New thermodynamic and functional studies have been recently conducted to evaluate the impact of amino acid substitutions on the Mitogen Activated Protein Kinases 1 and 3 (MAPK1/3). The Critical Assessment of Genome Interpretation (CAGI) data provider, at Sapienza University of Rome, measured the unfolding free energy and the enzymatic activity of a set of variants (MAPK challenge dataset). Thermodynamic measurements for the denaturant-induced equilibrium unfolding of the phosphorylated and unphosphorylated forms of the MAPKs were obtained by monitoring the far-UV circular dichroism and intrinsic fluorescence changes as a function of denaturant concentration. 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 ( ΔΔ G^H_2O ). The enzymatic activity of the phosphorylated MAPKs variants was also measured using Chelation-Enhanced Fluorescence to monitor the phosphorylation of a peptide substrate. The MAPK challenge dataset, composed of a total of 23 single amino acid substitutions (11 and 12 for MAPK1 and MAPK3, respectively), was used to assess the effectiveness of the computational methods in predicting the ΔΔ G^H_2O values, associated with the variants, and categorize them as destabilizing and not destabilizing. The data on the enzymatic activity of the MAPKs mutants were used to assess the performance of the methods for predicting the functional impact of the variants. For the sixth edition of CAGI, thirteen independent research groups from four continents (Asia, Australia, Europe and North America) submitted > 80 sets of predictions, obtained from different approaches. In this manuscript, we summarized the results of our assessment to highlight the possible limitations of the available algorithms.

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.

In this paper, we provide evidence that Zn2+ ions play a role in the SARS-CoV-2 virus strategy to escape the immune response mediated by the BST2-tetherin host protein. This conclusion is based on sequence analysis and molecular dynamics simulations as well as X-ray absorption experiments [1].

Abstract Dynamic changes in the genome characterise the onset, progression, diffusion and metastasis of cancer. Genomic changes may confer a growth advantage: gain of function mutations in oncogenes and loss of function in tumour suppressor may drive the tumour invasion and the loss of control mechanism in the hosts. Genome instability and inflammation are at the basis of the genetic changes reported in the Cancer Genome Atlas during these last 20 years. Cancer genome landscape has been extensively analysed with the aim to distinguish, in the tumour samples, the cancer‐related mutations from those neutral or benign. Studies in cancer tissues indicate that the somatic mutations observed effect from 33 to 66 genes and that 95% of these mutations are single‐nucleotide variants, mainly missense, or nonsynonymous variants, that is a single nucleotide substitution leads to a change in the protein sequence. The vast amount of data collected is not accompanied yet by a meaningful interpretation, and a linear connection between clinics and biological data is not even on the way. We report on case studies of missense mutations in frequently mutated cancer genes, in relation to the effect of missense mutations on the protein structure, functions and stability. Key Concepts Cancer is a complex disease characterised by dynamic changes in the genome. Germline and somatic protein variants are present in cancer. Missense mutations in oncogenes or in tumour suppressor genes may promote uncontrolled cellular growth. Missense mutations may affect protein structure, stability and function. Structural and functional analysis of protein missense mutations in cancer may help in the development of personalised drugs.

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.

This work studies the stability of wild-type frataxin and some of its variants found in cancer tissues upon Co2+ binding. Although the physiologically involved metal ion in the frataxin enzymatic activity is Fe2+, as it is customarily done, Co2+ is most often used in experiments because Fe2+ is extremely unstable owing to the fast oxidation reaction Fe2+ → Fe3+. Protein stability is monitored following the conformational changes induced by Co2+ binding as measured by circular dichroism, fluorescence spectroscopy, and melting temperature measurements. The stability ranking among the wild-type frataxin and its variants obtained in this way is confirmed by a detailed comparative analysis of the XAS spectra of the metal-protein complex at the Co K-edge. In particular, a fit to the EXAFS region of the spectrum allows positively identifying the frataxin acidic ridge as the most likely location of the metal-binding sites. Furthermore, we can explain the surprising feature emerging from a detailed analysis of the XANES region of the spectrum, showing that the longer 81-210 frataxin fragment has a smaller propensity for Co2+ binding than the shorter 90-210 one. This fact is explained by the peculiar role of the N-terminal disordered tail in modulating the protein ability to interact with the metal.

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.

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.

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.

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.

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.

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.

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.

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.

Far-UV Circular Dichroism experiments and Atomic Force Microscopy tomography are employed to assess the impact of β-sheet breakers on the Aβ1–40 peptide aggregation process in the presence of Cu2+ or Zn2+ transition metals. In this work we focus on two specific 5-amino acids long β-sheet breakers, namely the LPFFD Soto peptide, already known in the literature, and the LPFFN peptide recently designed and studied by our team. We provide evidence that both β-sheet breakers are effective in reducing the Aβ1–40 aggregation propensity, even in the presence of metal ions.

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.

Three PEGylated beta-sheet breaker peptides are designed as new inhibitors of beta-amyloid fibrillization. The peptide Ac-Leu-Pro-Phe-Phe-Asp-NH2, considered the lead compound, and hexamers in which taurine and beta-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 beta-sheet breaker peptides to inhibit fibrillogenesis in vitro. Moreover, microscopy images of beta-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.

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.