Large scale genome sequencing allowed the identification of a massive number of genetic variations, whose impact on human health is still unknown. In this review we analyze, by an in silico-based strategy, the impact of missense variants on cancer-related genes, whose effect on protein stability and function was experimentally determined. We collected a set of 164 variants from 11 proteins to analyze the impact of missense mutations at structural and functional levels, and to assess the performance of state-of-the-art methods (FoldX and Meta-SNP) for predicting protein stability change and pathogenicity. The result of our analysis shows that a combination of experimental data on protein stability and in silico pathogenicity predictions allowed the identification of a subset of variants with a high probability of having a deleterious phenotypic effect, as confirmed by the significant enrichment of the subset in variants annotated in the COSMIC database as putative cancer-driving variants. Our analysis suggests that the integration of experimental and computational approaches may contribute to evaluate the risk for complex disorders and develop more effective treatment strategies.
Pim-1 kinase, a serine/threonine protein kinase encoded by the pim proto-oncogene, is involved in several signalling pathways such as the regulation of cell cycle progression and apoptosis. Many cancer types show high expression levels of Pim kinases and particularly Pim-1 has been linked to the initiation and progression of the malignant phenotype. In several cancer tissues somatic Pim-1 mutants have been identified. These natural variants are nonsynonymous single nucleotide polymorphisms, variations of a single nucleotide occurring in the coding region and leading to amino acid substitutions. In this study we investigated the effect of amino acid substitution on the structural stability and on the activity of Pim-1 kinase. We expressed and purified some of the mutants of Pim-1 kinase that are expressed in cancer tissues and reported in the single nucleotide polymorphisms database. The point mutations in the variants significantly affect the conformation of the native state of Pim-1. All the mutants, expressed as soluble recombinant proteins, show a decreased thermal and thermodynamic stability and a lower activation energy values for kinase activity. The decreased stability accompanied by an increased flexibility suggests that Pim-1 variants may be involved in a wider network of protein interactions. All mutants bound ATP and ATP mimetic inhibitors with comparable IC50 values suggesting that the studied Pim-1 kinase mutants can be efficiently targeted with inhibitors developed for the wild type protein.
Cancer cells are able to survive in difficult conditions, reprogramming their metabolism according to their requirements. Under hypoxic conditions they shift from oxidative phosphorylation to aerobic glycolysis, a behavior known as Warburg effect. In the last years, glycolytic enzymes have been identified as potential targets for alternative anticancer therapies. Recently, phosphoglycerate kinase 1 (PGK1), an ubiquitous enzyme expressed in all somatic cells that catalyzes the seventh step of glycolysis which consists of the reversible phosphotransfer reaction from 1,3-bisphosphoglycerate to ADP, has been discovered to be overexpressed in many cancer types. Moreover, several somatic variants of PGK1 have been identified in tumors. In this study we analyzed the effect of the single nucleotide variants found in cancer tissues on the PGK1 structure and function. Our results clearly show that the variants display a decreased catalytic efficiency and/or thermodynamic stability and an altered local tertiary structure, as shown by the solved X-ray structures. The changes in the catalytic properties and in the stability of the PGK1 variants, mainly due to the local changes evidenced by the X-ray structures, suggest also changes in the functional role of PGK to support the biosynthetic need of the growing and proliferating tumour cells.
Protein tyrosine phosphatase ρ (PTPρ) belongs to the classical receptor type IIB family of protein tyrosine phosphatase, the most frequently mutated tyrosine phosphatase in human cancer. There are evidences to suggest that PTPρ may act as a tumor suppressor gene and dysregulation of Tyr phosphorylation can be observed in diverse diseases, such as diabetes, immune deficiencies and cancer. PTPρ variants in the catalytic domain have been identified in cancer tissues. These natural variants are nonsynonymous single nucleotide polymorphisms, variations of a single nucleotide occurring in the coding region and leading to amino acid substitutions. In this study we investigated the effect of amino acid substitution on the structural stability and on the activity of the membrane-proximal catalytic domain of PTPρ. We expressed and purified as soluble recombinant proteins some of the mutants of the membrane-proximal catalytic domain of PTPρ identified in colorectal cancer and in the single nucleotide polymorphisms database. The mutants show a decreased thermal and thermodynamic stability and decreased activation energy relative to phosphatase activity, when compared to wild- type. All the variants show three-state equilibrium unfolding transitions similar to that of the wild- type, with the accumulation of a folding intermediate populated at ∼4.0 M urea.
Lysine acetylation is an important epigenetic mark regulating gene transcription and chromatin structure. Acetylated lysine residues are specifically recognized by bromodomains, small protein interaction modules that read these modification in a sequence and acetylation dependent way regulating the recruitment of transcriptional regulators and chromatin remodelling enzymes to acetylated sites in chromatin. Recent studies revealed that bromodomains are highly druggable protein interaction domains resulting in the development of a large number of bromodomain inhibitors. BET bromodomain inhibitors received a lot of attention in the oncology field resulting in the rapid translation of early BET bromodomain inhibitors into clinical studies. Here we investigated the effects of mutations present as polymorphism or found in cancer on BET bromodomain function and stability and the influence of these mutants on inhibitor binding. We found that most BET missense mutations localize to peripheral residues in the two terminal helices. Crystal structures showed that the three dimensional structure is not compromised by these mutations but mutations located in close proximity to the acetyl-lysine binding site modulate acetyl-lysine and inhibitor binding. Most mutations affect significantly protein stability and tertiary structure in solution, suggesting new interactions and an alternative network of protein-protein interconnection as a consequence of single amino acid substitution. To our knowledge this is the first report studying the effect of mutations on bromodomain function and inhibitor binding.
We present in this work a first X‐ray Absorption Spectroscopy study of the interactions of Zn with human BST2/tetherin and SARS‐CoV‐2 orf7a proteins as well as with some of their complexes. The analysis of the XANES region of the measured spectra shows that Zn binds to BST2, as well as to orf7a, thus resulting in the formation of BST2‐orf7a complexes. This structural information confirms the the conjecture, recently put forward by some of the present Authors, according to which the accessory orf7a (and possibly also orf8) viral protein are capable of interfering with the BST2 antiviral activity. Our explanation for this behavior is that, when BST2 gets in contact with Zn bound to the orf7a Cys15 ligand, it has the ability of displacing the metal owing to the creation of a new disulfide bridge across the two proteins. The formation of this BST2‐orf7a complex destabilizes BST2 dimerization, thus impairing the antiviral activity of the latter.