The complex impact of cancer-related missense mutations on the stability and on the biophysical and biochemical properties of MAPK1 and MAPK3 somatic variants.

Mitogen-activated protein kinases 1 and 3 (MAPK1 and MAPK3), also called extracellular regulated kinases (ERK2 and ERK1), are serine/threonine kinase activated downstream by the Ras/Raf/MEK/ERK signal transduction cascade that regulates a variety of cellular processes. A dysregulation of MAPK cascade is frequently associated to missense mutations on its protein components and may be related to many pathologies, including cancer. In this study we selected from COSMIC database a set of MAPK1 and MAPK3 somatic variants found in cancer tissues carrying missense mutations distributed all over the MAPK1 and MAPK3 sequences. The proteins were expressed as pure recombinant proteins, and their biochemical and biophysical properties have been studied in comparison with the wild type. The missense mutations lead to changes in the tertiary arrangements of all the variants. The thermodynamic stability of the wild type and variants has been investigated in the non-phosphorylated and in the phosphorylated form. Significant differences in the thermal stabilities of most of the variants have been observed, as well as changes in the catalytic efficiencies. The energetics of the catalytic reaction is affected for all the variants for both the MAPK proteins. The stability changes and the variation in the enzyme catalysis observed for most of MAPK1/3 variants suggest that a local change in a residue, distant from the catalytic site, may have long-distance effects that reflect globally on enzyme stability and functions.

Mutation in the Common Docking Domain Affects MAP Kinase ERK2 Catalysis and Stability.

The extracellular-signal-regulated kinase 2 (ERK2), a mitogen-activated protein kinase (MAPK) located downstream of the Ras-Raf-MEK-ERK signal transduction cascade, is involved in the regulation of a large variety of cellular processes. The ERK2, activated by phosphorylation, is the principal effector of a central signaling cascade that converts extracellular stimuli into cells. Deregulation of the ERK2 signaling pathway is related to many human diseases, including cancer. This study reports a comprehensive biophysical analysis of structural, function, and stability data of pure, recombinant human non-phosphorylated (NP-) and phosphorylated (P-) ERK2 wild-type and missense variants in the common docking site (CD-site) found in cancer tissues. Because the CD-site is involved in interaction with protein substrates and regulators, a biophysical characterization of missense variants adds information about the impact of point mutations on the ERK2 structure-function relationship. Most of the P-ERK2 variants in the CD-site display a reduced catalytic efficiency, and for the P-ERK2 D321E, D321N, D321V and E322K, changes in thermodynamic stability are observed. The thermal stability of NP-ERK2 and P-ERK2 D321E, D321G, and E322K is decreased with respect to the wild-type. In general, a single residue mutation in the CD-site may lead to structural local changes that reflects in alterations in the global ERK2 stability and catalysis.

Zn-Induced Interactions Between SARS-CoV-2 orf7a and BST2/Tetherin.

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.

Analysis and Interpretation of the Impact of Missense Variants in Cancer.

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.

A Glimpse into the Structural Properties of the Intermediate and Transition State in the Folding of Bromodomain 2 Domain 2 by Φ Value Analysis.

Bromodomains (BRDs) are small protein interaction modules of about 110 amino acids that selectively recognize acetylated lysine in histones and other proteins. These domains have been identified in a variety of multi-domain proteins involved in transcriptional regulation or chromatin remodeling in eukaryotic cells. BRD inhibition is considered an attractive therapeutic approach in epigenetic disorders, particularly in oncology. Here, we present a Φ value analysis to investigate the folding pathway of the second domain of BRD2 (BRD2(2)). Using an extensive mutational analysis based on 25 site-directed mutants, we provide structural information on both the intermediate and late transition state of BRD2(2). The data reveal that the C-terminal region represents part of the initial folding nucleus, while the N-terminal region of the domain consolidates its structure only later in the folding process. Furthermore, only a small number of native-like interactions have been identified, suggesting the presence of a non-compact, partially folded state with scarce native-like characteristics. Taken together, these results indicate that, in BRD2(2), a hierarchical mechanism of protein folding can be described with non-native interactions that play a significant role in folding.

Dealing with Cu reduction in X-ray absorption spectroscopy experiments.

In this paper we prove in the exemplary case of the amyloid-ß peptide in complex with Cu(ii) that at the current low temperatures employed in XAS experiments, the time needed for collecting a good quality XAS spectrum is significantly shorter than the time after which structural damage becomes appreciable. Our method takes advantage of the well-known circumstance that the transition of Cu from the oxidized to the reduced form under ionizing radiation can be quantified by monitoring a characteristic peak in the pre-edge region. We show that there exists a sufficiently large time window in which good XAS spectra can be acquired before the structure around the oxidized Cu(ii) ion reorganizes itself into the reduced Cu(i) "resting" structure. We suggest that similar considerations apply to other cases of biological interest, especially when dealing with macromolecules in complex with transition metal ions.

Evaluating the predictions of the protein stability change upon single amino acid substitutions for the FXN CAGI5 challenge.

Frataxin (FXN) is a highly conserved protein found in prokaryotes and eukaryotes that is required for efficient regulation of cellular iron homeostasis. Experimental evidence associates amino acid substitutions of the FXN to Friedreich Ataxia, a neurodegenerative disorder. Recently, new thermodynamic experiments have been performed to study the impact of somatic variations identified in cancer tissues on protein stability. The Critical Assessment of Genome Interpretation (CAGI) data provider at the University of Rome measured the unfolding free energy of a set of variants (FXN challenge data set) with far-UV circular dichroism and intrinsic fluorescence spectra. These values have been used to calculate the change in unfolding free energy between the variant and wild-type proteins at zero concentration of denaturant (ΔΔGH2O) . The FXN challenge data set, composed of eight amino acid substitutions, was used to evaluate the performance of the current computational methods for predicting the ΔΔGH2O value associated with the variants and to classify them as destabilizing and not destabilizing. For the fifth edition of CAGI, six independent research groups from Asia, Australia, Europe, and North America submitted 12 sets of predictions from different approaches. In this paper, we report the results of our assessment and discuss the limitations of the tested algorithms.

Characterization of human frataxin missense variants in cancer tissues.

Human frataxin is an iron-binding protein involved in the mitochondrial iron-sulfur (Fe-S) clusters assembly, a process fundamental for the functional activity of mitochondrial proteins. Decreased level of frataxin expression is associated with the neurodegenerative disease Friedreich ataxia. Defective function of frataxin may cause defects in mitochondria, leading to increased tumorigenesis. Tumor-initiating cells show higher iron uptake, a decrease in iron storage and a reduced Fe-S clusters synthesis and utilization. In this study, we selected, from COSMIC database, the somatic human frataxin missense variants found in cancer tissues p.D104G, p.A107V, p.F109L, p.Y123S, p.S161I, p.W173C, p.S181F, and p.S202F to analyze the effect of the single amino acid substitutions on frataxin structure, function, and stability. The spectral properties, the thermodynamic and the kinetic stability, as well as the molecular dynamics of the frataxin missense variants found in cancer tissues point to local changes confined to the environment of the mutated residues. The global fold of the variants is not altered by the amino acid substitutions; however, some of the variants show a decreased stability and a decreased functional activity in comparison with that of the wild-type protein.

The phosphoglycerate kinase 1 variants found in carcinoma cells display different catalytic activity and conformational stability compared to the native enzyme.

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.

Unveiling the folding mechanism of the Bromodomains.

Bromodomains (BRDs) are small protein domains often present in large multidomain proteins involved in transcriptional regulation in eukaryotic cells. They currently represent valuable targets for the development of inhibitors of aberrant transcriptional processes in a variety of human diseases. Here we report urea-induced equilibrium unfolding experiments monitored by circular dichroism (CD) and fluorescence on two structurally similar BRDs: BRD2(2) and BRD4(1), showing that BRD4(1) is more stable than BRD2(2). Moreover, we report a description of their kinetic folding mechanism, as obtained by careful analysis of stopped-flow and temperature-jump data. The presence of a high energy intermediate for both proteins, suggested by the non-linear dependence of the folding rate on denaturant concentration in the millisec time regime, has been experimentally observed by temperature-jump experiments. Quantitative global analysis of all the rate constants obtained over a wide range of urea concentrations, allowed us to propose a common, three-state, folding mechanism for these two BRDs. Interestingly, the intermediate of BRD4(1) appears to be more stable and structurally native-like than that populated by BRD2(2). Our results underscore the role played by structural topology and sequence in determining and tuning the folding mechanism.

Single-Nucleotide Polymorphism of PPARγ, a Protein at the Crossroads of Physiological and Pathological Processes.

Genome polymorphisms are responsible for phenotypic differences between humans and for individual susceptibility to genetic diseases and therapeutic responses. Non-synonymous single-nucleotide polymorphisms (nsSNPs) lead to protein variants with a change in the amino acid sequence that may affect the structure and/or function of the protein and may be utilized as efficient structural and functional markers of association to complex diseases. This study is focused on nsSNP variants of the ligand binding domain of PPARγ a nuclear receptor in the superfamily of ligand inducible transcription factors that play an important role in regulating lipid metabolism and in several processes ranging from cellular differentiation and development to carcinogenesis. Here we selected nine nsSNPs variants of the PPARγ ligand binding domain, V290M, R357A, R397C, F360L, P467L, Q286P, R288H, E324K, and E460K, expressed in cancer tissues and/or associated with partial lipodystrophy and insulin resistance. The effects of a single amino acid change on the thermodynamic stability of PPARγ, its spectral properties, and molecular dynamics have been investigated. The nsSNPs PPARγ variants show alteration of dynamics and tertiary contacts that impair the correct reciprocal positioning of helices 3 and 12, crucially important for PPARγ functioning.

PEGylated ß-Sheet Breaker Peptides as Inhibitors of ß-Amyloid Fibrillization.

Three PEGylated ß-sheet breaker peptides are designed as new inhibitors of ß-amyloid fibrillization. The peptide Ac-Leu-Pro-Phe-Phe-Asp-NH2 , considered the lead compound, and hexamers in which taurine and ß-alanine substitute the acetyl group, are conjugated to poly(ethylene glycol); this conjugates self-assemble into nanoparticles. The activity of the PEGylated peptides as inhibitors of amyloid fibrillization are tested in vitro using circular dichroism spectroscopy and scanning electron microscopy. The experimental results indicate that PEGylation does not impair the ability of the ß-sheet breaker peptides to inhibit fibrillogenesis in vitro. Moreover, microscopy images of ß-amyloid incubated for 6 days with the taurine-containing peptide, suggest that this conjugate has major anti-fibrillogenesis activity and demonstrate the important role of the sulfonamide function against the amyloid aggregation.

Production of functional, stable, unmutated recombinant human papillomavirus E6 oncoprotein: implications for HPV-tumor diagnosis and therapy.

BACKGROUND: High-risk human papillomaviruses (HR-HPVs) types 16 and 18 are the main etiological agents of cervical cancer, with more than 550,000 new cases each year worldwide. HPVs are also associated with other ano-genital and head-and-neck tumors. The HR-HPV E6 and E7 oncoproteins are responsible for onset and maintenance of the cell transformation state, and they represent appropriate targets for development of diagnostic and therapeutic tools. METHODS: The unmutated E6 gene from HPV16 and HPV18 and from low-risk HPV11 was cloned in a prokaryotic expression vector for expression of the Histidine-tagged E6 protein (His6-E6), according to a novel procedure. The structural properties were determined using circular dichroism and fluorescence spectroscopy. His6-E6 oncoprotein immunogenicity was assessed in a mouse model, and its functionality was determined using in vitro GST pull-down and protein degradation assays. RESULTS: The His6-tagged E6 proteins from HPV16, HPV18, and HPV11 E6 genes, without any further modification in the amino-acid sequence, were produced in bacteria as soluble and stable molecules. Structural analyses of HPV16 His6-E6 suggests that it maintains correct folding and conformational properties. C57BL/6 mice immunized with HPV16 His6-E6 developed significant humoral immune responses. The E6 proteins from HPV16, HPV18, and HPV11 were purified according to a new procedure, and investigated for protein-protein interactions. HR-HPV His6-E6 bound p53, the PDZ1 motif from MAGI-1 proteins, the human discs large tumor suppressor, and the human ubiquitin ligase E6-associated protein, thus suggesting that it is biologically active. The purified HR-HPV E6 proteins also targeted the MAGI-3 and p53 proteins for degradation. CONCLUSIONS: This new procedure generates a stable, unmutated HPV16 E6 protein, which maintains the E6 properties in in vitro binding assays. This will be useful for basic studies, and for development of diagnostic kits and immunotherapies in preclinical mouse models of HPV-related tumorigenesis.

Effect of BET Missense Mutations on Bromodomain Function, Inhibitor Binding and Stability.

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.

Structural basis of the transactivation deficiency of the human PPARγ F360L mutant associated with familial partial lipodystrophy.

The peroxisome proliferator-activated receptors (PPARs) are transcription factors that regulate glucose and lipid metabolism. The role of PPARs in several chronic diseases such as type 2 diabetes, obesity and atherosclerosis is well known and, for this reason, they are the targets of antidiabetic and hypolipidaemic drugs. In the last decade, some rare mutations in human PPARγ that might be associated with partial lipodystrophy, dyslipidaemia, insulin resistance and colon cancer have emerged. In particular, the F360L mutant of PPARγ (PPARγ2 residue 388), which is associated with familial partial lipodystrophy, significantly decreases basal transcriptional activity and impairs stimulation by synthetic ligands. To date, the structural reason for this defective behaviour is unclear. Therefore, the crystal structure of PPARγ F360L together with the partial agonist LT175 has been solved and the mutant has been characterized by circular-dichroism spectroscopy (CD) in order to compare its thermal stability with that of the wild-type receptor. The X-ray analysis showed that the mutation induces dramatic conformational changes in the C-terminal part of the receptor ligand-binding domain (LBD) owing to the loss of van der Waals interactions made by the Phe360 residue in the wild type and an important salt bridge made by Arg357, with consequent rearrangement of loop 11/12 and the activation function helix 12 (H12). The increased mobility of H12 makes the binding of co-activators in the hydrophobic cleft less efficient, thereby markedly lowering the transactivation activity. The spectroscopic analysis in solution and molecular-dynamics (MD) simulations provided results which were in agreement and consistent with the mutant conformational changes observed by X-ray analysis. Moreover, to evaluate the importance of the salt bridge made by Arg357, the crystal structure of the PPARγ R357A mutant in complex with the agonist rosiglitazone has been solved.

Computational and experimental studies on ß-sheet breakers targeting Aß1-40 fibrils.

In this work we present and compare the results of extensive molecular dynamics simulations of model systems comprising an Aß1-40 peptide in water in interaction with short peptides (ß-sheet breakers) mimicking the 17-21 region of the Aß1-40 sequence. Various systems differing in the customized ß-sheet breaker structure have been studied. Specifically we have considered three kinds of ß-sheet breakers, namely Ac-LPFFD-NH2 and two variants thereof, one obtained by substituting the acetyl group with the sulfonic amino acid taurine (Tau-LPFFD-NH2) and a second novel one in which the aspartic acid is substituted by an asparagine (Ac-LPFFN-NH2). Thioflavin T fluorescence, circular dichroism, and mass spectrometry experiments have been performed indicating that ß-sheet breakers are able to inhibit in vitro fibril formation and prevent the ß sheet folding of portions of the Aß1-40 peptide. We show that molecular dynamics simulations and far UV circular dichroism provide consistent evidence that the new Ac-LPFFN-NH2 ß-sheet breaker is more effective than the other two in stabilizing the native α-helix structure of Aß1-40. In agreement with these results thioflavin T fluorescence experiments confirm the higher efficiency in inhibiting Aß1-40 aggregation. Furthermore, mass spectrometry data and molecular dynamics simulations consistently identified the 17-21 Aß1-40 portion as the location of the interaction region between peptide and the Ac-LPFFN-NH2 ß-sheet breaker.

Effect of single amino acid substitution observed in cancer on Pim-1 kinase thermodynamic stability and structure.

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.

Structural stability of human protein tyrosine phosphatase ρ catalytic domain: effect of point mutations.

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.

Inhibition of amyloid peptide fragment Aß25-35 fibrillogenesis and toxicity by N-terminal ß-amino acid-containing esapeptides: is taurine moiety essential for in vivo effects?

We report the synthesis and fibrillogenesis inhibiting activity of the new peptide derivatives 1-4, related to the pentapeptide Ac-LPFFD-NH(2) (iAß5p), proposed by Soto and co-workers and widely recognized as one of the most active ß-sheet breaker agents. The Aß(25-35) fragment of the parent full-length Aß(1-42) was used as fibrillogenesis model. The activity of peptide derivatives 1-4 was tested in vitro by thioflavin T binding assay, far UV CD spectroscopy, and scanning electron microscopy. Their ability to hinder the toxic effect of Aß(25-35) in vivo was studied by monitoring the viability of human SH-SY5Y neuroblastoma cells and the prevention of superoxide anion radical release from BV2 microglial cells. The results point to a favourable role in the fibrillogenesis inhibitory activity of the sulphonamide junction for compounds 1 and 2, containing an N,N-dimethyltaurine and a taurine amino-terminal moiety, respectively. Furthermore, compounds 1 and 2 show a significant protective effect on cell viability, rescuing the cells from the toxicity exerted by Aß(25-35) treatment.