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.

K-Pro: Kinetics Data on Proteins and Mutants.

The study of protein folding plays a crucial role in improving our understanding of protein function and of the relationship between genetics and phenotypes. In particular, understanding the thermodynamics and kinetics of the folding process is important for uncovering the mechanisms behind human disorders caused by protein misfolding. To address this issue, it is essential to collect and curate experimental kinetic and thermodynamic data on protein folding. K-Pro is a new database designed for collecting and storing experimental kinetic data on monomeric proteins, with a two-state folding mechanism. With 1,529 records from 62 proteins corresponding to 65 structures, K-Pro contains various kinetic parameters such as the logarithm of the folding and unfolding rates, Tanford's ß and the ϕ values. When available, the database also includes thermodynamic parameters associated with the kinetic data. K-Pro features a user-friendly interface that allows browsing and downloading kinetic data of interest. The graphical interface provides a visual representation of the protein and mutants, and it is cross-linked to key databases such as PDB, UniProt, and PubMed. K-Pro is open and freely accessible through https://folding.biofold.org/k-pro and supports the latest versions of popular browsers.

Resources and tools for rare disease variant interpretation.

Collectively, rare genetic disorders affect a substantial portion of the world's population. In most cases, those affected face difficulties in receiving a clinical diagnosis and genetic characterization. The understanding of the molecular mechanisms of these diseases and the development of therapeutic treatments for patients are also challenging. However, the application of recent advancements in genome sequencing/analysis technologies and computer-aided tools for predicting phenotype-genotype associations can bring significant benefits to this field. In this review, we highlight the most relevant online resources and computational tools for genome interpretation that can enhance the diagnosis, clinical management, and development of treatments for rare disorders. Our focus is on resources for interpreting single nucleotide variants. Additionally, we present use cases for interpreting genetic variants in clinical settings and review the limitations of these results and prediction tools. Finally, we have compiled a curated set of core resources and tools for analyzing rare disease genomes. Such resources and tools can be utilized to develop standardized protocols that will enhance the accuracy and effectiveness of rare disease diagnosis.

Modulation of the H+/ATP coupling ratio by ADP and ATP as a possible regulatory feature in the F-type ATP synthases.

F-type ATP synthases are transmembrane enzymes, which play a central role in the metabolism of all aerobic and photosynthetic cells and organisms, being the major source of their ATP synthesis. Catalysis occurs via a rotary mechanism, in which the free energy of a transmembrane electrochemical ion gradient is converted into the free energy of ATP phosphorylation from ADP and Pi, and vice versa. An ADP, tightly bound to one of the three catalytic sites on the stator head, is associated with catalysis inhibition, which is relieved by the transmembrane proton gradient and by ATP. By preventing wasteful ATP hydrolysis in times of low osmotic energy and low ATP/ADP ratio, such inhibition constitutes a classical regulatory feedback effect, likely to be an integral component of in vivo regulation. The present miniview focuses on an additional putative regulatory phenomenon, which has drawn so far little attention, consisting in a substrate-induced tuning of the H+/ATP coupling ratio during catalysis, which might represent an additional key to energy homeostasis in the cell. Experimental pieces of evidence in support of such a phenomenon are reviewed.

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.

ThermoScan: Semi-automatic Identification of Protein Stability Data From PubMed.

During the last years, the increasing number of DNA sequencing and protein mutagenesis studies has generated a large amount of variation data published in the biomedical literature. The collection of such data has been essential for the development and assessment of tools predicting the impact of protein variants at functional and structural levels. Nevertheless, the collection of manually curated data from literature is a highly time consuming and costly process that requires domain experts. In particular, the development of methods for predicting the effect of amino acid variants on protein stability relies on the thermodynamic data extracted from literature. In the past, such data were deposited in the ProTherm database, which however is no longer maintained since 2013. For facilitating the collection of protein thermodynamic data from literature, we developed the semi-automatic tool ThermoScan. ThermoScan is a text mining approach for the identification of relevant thermodynamic data on protein stability from full-text articles. The method relies on a regular expression searching for groups of words, including the most common conceptual words appearing in experimental studies on protein stability, several thermodynamic variables, and their units of measure. ThermoScan analyzes full-text articles from the PubMed Central Open Access subset and calculates an empiric score that allows the identification of manuscripts reporting thermodynamic data on protein stability. The method was optimized on a set of publications included in the ProTherm database, and tested on a new curated set of articles, manually selected for presence of thermodynamic data. The results show that ThermoScan returns accurate predictions and outperforms recently developed text-mining algorithms based on the analysis of publication abstracts. Availability: The ThermoScan server is freely accessible online at https://folding.biofold.org/thermoscan. The ThermoScan python code and the Google Chrome extension for submitting visualized PMC web pages to the ThermoScan server are available at https://github.com/biofold/ThermoScan.

Limitations and challenges in protein stability prediction upon genome variations: towards future applications in precision medicine.

Protein stability predictions are becoming essential in medicine to develop novel immunotherapeutic agents and for drug discovery. Despite the large number of computational approaches for predicting the protein stability upon mutation, there are still critical unsolved problems: 1) the limited number of thermodynamic measurements for proteins provided by current databases; 2) the large intrinsic variability of ΔΔG values due to different experimental conditions; 3) biases in the development of predictive methods caused by ignoring the anti-symmetry of ΔΔG values between mutant and native protein forms; 4) over-optimistic prediction performance, due to sequence similarity between proteins used in training and test datasets. Here, we review these issues, highlighting new challenges required to improve current tools and to achieve more reliable predictions. In addition, we provide a perspective of how these methods will be beneficial for designing novel precision medicine approaches for several genetic disorders caused by mutations, such as cancer and neurodegenerative diseases.

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.

Thermodynamics of proton transport coupled ATP synthesis.

The thermodynamic H(+)/ATP ratio of the H(+)-ATP synthase from chloroplasts was measured in proteoliposomes after energization of the membrane by an acid base transition (Turina et al. 2003 [13], 418-422). The method is discussed, and all published data obtained with this system are combined and analyzed as a single dataset. This meta-analysis led to the following results. 1) At equilibrium, the transmembrane ΔpH is energetically equivalent to the transmembrane electric potential difference. 2) The standard free energy for ATP synthesis (reference reaction) is ΔG°(ref)=33.8±1.3kJ/mol. 3) The thermodynamic H(+)/ATP ratio, as obtained from the shift of the ATP synthesis equilibrium induced by changing the transmembrane ΔpH (varying either pH(in) or pH(out)) is 4.0±0.1. The structural H(+)/ATP ratio, calculated from the ratio of proton binding sites on the c-subunit-ring in F(0) to the catalytic nucleotide binding sites on the ß-subunits in F(1), is c/ß=14/3=4.7. We infer that the energy of 0.7 protons per ATP that flow through the enzyme, but do not contribute to shifting the ATP/(ADP·Pi) ratio, is used for additional processes within the enzyme, such as activation, and/or energy dissipation, due e.g. to internal uncoupling. The ratio between the thermodynamic and the structural H(+)/ATP values is 0.85, and we conclude that this value represents the efficiency of the chemiosmotic energy conversion within the chloroplast H(+)-ATP synthase.

Dehydration affects the electronic structure of the primary electron donor in bacterial photosynthetic reaction centers: evidence from visible-NIR and light-induced difference FTIR spectroscopy.

The photosynthetic reaction center (RC) is a membrane pigment-protein complex that catalyzes the initial charge separation reactions of photosynthesis. Following photoexcitation, the RC undergoes conformational relaxations which stabilize the charge-separated state. Dehydration of the complex inhibits its conformational dynamics, providing a useful tool to gain insights into the relaxational processes. We analyzed the effects of dehydration on the electronic structure of the primary electron donor P, as probed by visible-NIR and light-induced FTIR difference spectroscopy, in RC films equilibrated at different relative humidities r. Previous FTIR and ENDOR spectroscopic studies revealed that P, an excitonically coupled dimer of bacteriochlorophylls, can be switched between two conformations, P866 and P850, which differ in the extent of delocalization of the unpaired electron between the two bacteriochlorophyll moieties (PL and PM) of the photo-oxidized radical P(+). We found that dehydration (at r = 11%) shifts the optical Qy band of P from 866 to 850-845 nm, a large part of the effect occurring already at r = 76%. Such a dehydration weakens light-induced difference FTIR marker bands, which probe the delocalization of charge distribution within the P(+) dimer (the electronic band of P(+) at 2700 cm(-1), and the associated phase-phonon vibrational modes at around 1300, 1480, and 1550 cm(-1)). From the analysis of the P(+) keto C[double bond, length as m-dash]O bands at 1703 and 1713-15 cm(-1), we inferred that dehydration induces a stronger localization of the unpaired electron on PL(+). The observed charge redistribution is discussed in relation to the dielectric relaxation of the photoexcited RC on a long (10(2) s) time scale.

Comparison of the H+/ATP ratios of the H+-ATP synthases from yeast and from chloroplast.

F(0)F(1)-ATP synthases use the free energy derived from a transmembrane proton transport to synthesize ATP from ADP and inorganic phosphate. The number of protons translocated per ATP (H(+)/ATP ratio) is an important parameter for the mechanism of the enzyme and for energy transduction in cells. Current models of rotational catalysis predict that the H(+)/ATP ratio is identical to the stoichiometric ratio of c-subunits to ß-subunits. We measured in parallel the H(+)/ATP ratios at equilibrium of purified F(0)F(1)s from yeast mitochondria (c/ß = 3.3) and from spinach chloroplasts (c/ß = 4.7). The isolated enzymes were reconstituted into liposomes and, after energization of the proteoliposomes with acid-base transitions, the initial rates of ATP synthesis and hydrolysis were measured as a function of ΔpH. The equilibrium ΔpH was obtained by interpolation, and from its dependency on the stoichiometric ratio, [ATP]/([ADP]·[P(i)]), finally the thermodynamic H(+)/ATP ratios were obtained: 2.9 ± 0.2 for the mitochondrial enzyme and 3.9 ± 0.3 for the chloroplast enzyme. The data show that the thermodynamic H(+)/ATP ratio depends on the stoichiometry of the c-subunit, although it is not identical to the c/ß ratio.

Quantitative evaluation of the intrinsic uncoupling modulated by ADP and P(i) in the reconstituted ATP synthase of Escherichia coli.

The ATP synthase from Escherichia coli was isolated and reconstituted into liposomes. The ATP hydrolysis by these proteoliposomes was coupled to proton pumping, and the ensuing inner volume acidification was measured by the fluorescent probe 9-amino-6-chloro-2-methoxyacridine (ACMA). The ACMA response was calibrated by acid-base transitions, and converted into internal pH values. The rates of internal acidification and of ATP hydrolysis were measured in parallel, as a function of P(i) or ADP concentration. Increasing P(i) monotonically inhibited the hydrolysis rate with a half-maximal effect at 510µM, whereas it stimulated the acidification rate up to 100-200µM, inhibiting it only at higher concentrations. The ADP concentration in the assay, due both to contaminant ADP in ATP and to the hydrolysis reaction, was progressively decreased by means of increasing pyruvate kinase activities. Decreasing ADP stimulated the hydrolysis rate, whereas it inhibited the internal acidification rate. The quantitative analysis showed that the relative number of translocated protons per hydrolyzed ATP, i.e. the relative coupling ratio, depended on the concentrations of P(i) and ADP with apparent K(d) values of 220µM and 27nM respectively. At the smallest ADP concentrations reached, and in the absence of P(i), the coupling ratio dropped down to 15% relative to the value observed at the highest ADP and P(i) concentrations tested. In addition, the data indicate the presence of two ADP and P(i) binding sites, of which only the highest affinity one is related to changes in the coupling ratio.

Proton transport coupled ATP synthesis by the purified yeast H+ -ATP synthase in proteoliposomes.

The H(+)/ATP synthase from yeast mitochondria, MF0F1, was purified and reconstituted into liposomes prepared from phosphatidylcholine and phosphatidic acid. Analysis by mass spectrometry revealed the presence of all subunits of the yeast enzyme with the exception of the K-subunit. The MF0F1 liposomes were energized by acid-base transitions (DeltapH) and a K(+)/valinomycin diffusion potential (Deltaphi). ATP synthesis was completely abolished by the addition of uncouplers as well as by the inhibitor oligomycin. The rate of ATP synthesis was optimized as a function of various parameters and reached a maximum value (turnover number) of 120s⁻¹ at a transmembrane pH difference of 3.2 units (at pH(in)=4.8 and pH(out)=8.0) and a Deltaphi of 133mV (Nernst potential). Functional studies showed that the monomeric MF0F1, was fully active in ATP synthesis. The turnover increased in a sigmoidal way with increasing internal and decreasing external proton concentration. The dependence of the turnover on the phosphate concentration and the dependence of K(M) on pH(out) indicated that the substrate for ATP synthesis is the monoanionic phosphate species H2PO⁻4.

Intrinsic uncoupling in the ATP synthase of Escherichia coli.

The ATP hydrolysis activity and proton pumping of the ATP synthase of Escherichia coli in isolated native membranes have been measured and compared as a function of ADP and Pi concentration. The ATP hydrolysis activity was inhibited by Pi with an half-maximal effect at 140 microM, which increased progressively up in the millimolar range when the ADP concentration was progressively decreased by increasing amounts of an ADP trap. In addition, the relative extent of this inhibition decreased with decreasing ADP. The half-maximal inhibition by ADP was found in the submicromolar range, and the extent of inhibition was enhanced by the presence of Pi. The parallel measurement of ATP hydrolysis activity and proton pumping indicated that, while the rate of ATP hydrolysis was decreased as a function of either ligand, the rate of proton pumping increased. The latter showed a biphasic response to the concentration of Pi, in which an inhibition followed the initial stimulation. Similarly as previously found for the ATP synthase from Rhodobacter caspulatus [P. Turina, D. Giovannini, F. Gubellini, B.A. Melandri, Physiological ligands ADP and Pi modulate the degree of intrinsic coupling in the ATP synthase of the photosynthetic bacterium Rhodobacter capsulatus, Biochemistry 43 (2004) 11126-11134], these data indicate that the E. coli ATP synthase can operate at different degrees of energetic coupling between hydrolysis and proton transport, which are modulated by ADP and Pi.

The thermodynamic H+/ATP ratios of the H+-ATPsynthases from chloroplasts and Escherichia coli.

The H(+)/ATP ratio is an important parameter for the energy balance of all cells and for the coupling mechanism between proton transport and ATP synthesis. A straightforward interpretation of rotational catalysis predicts that the H(+)/ATP coincides with the ratio of the c-subunits to beta-subunits, implying that, for the chloroplast and Escherichia coli ATPsynthases, numbers of 4.7 and 3.3 are expected. Here, the energetics described by the chemiosmotic theory was used to determine the H(+)/ATP ratio for the two enzymes. The isolated complexes were reconstituted into liposomes, and parallel measurements were performed under identical conditions. The internal phase of the liposomes was equilibrated with the acidic medium during reconstitution, allowing to measure the internal pH with a glass electrode. An acid-base transition was carried out and the initial rates of ATP synthesis or ATP hydrolysis were measured with luciferin/luciferase as a function of DeltapH at constant Q = [ATP]/([ADP][P(i)]). From the shift of the equilibrium DeltapH as a function of Q the standard Gibbs free energy for phosphorylation, DeltaG(p)(0)'; and the H(+)/ATP ratio were determined. It resulted DeltaG(p)(0)' = 38 +/- 3 kJ.mol(-1) and H(+)/ATP = 4.0 +/- 0.2 for the chloroplast and H(+)/ATP = 4.0 +/- 0.3 for the E. coli enzyme, indicating that the thermodynamic H(+)/ATP ratio is the same for both enzymes and that it is different from the subunit stoichiometric ratio.

Heterogeneity of photosynthetic membranes from Rhodobacter capsulatus: size dispersion and ATP synthase distribution.

The density distribution of photosynthetic membrane vesicles (chromatophores) from Rhodobacter capsulatus has been studied by isopicnic centrifugation. The average vesicle diameters, examined by electron microscopy, varied between 61 and 72 nm in different density fractions (70 nm in unfractionated chromatophores). The ATP synthase catalytic activities showed maxima displaced toward the higher density fractions relative to bacteriochlorophyll, resulting in higher specific activities in those fractions (about threefold). The amount of ATP synthase, measured by quantitative Western blotting, paralleled the catalytic activities. The average number of ATP synthases per chromatophore, evaluated on the basis of the Western blotting data and of vesicle density analysis, ranged between 8 and 13 (10 in unfractionated chromatophores). Poisson distribution analysis indicated that the probability of chromatophores devoid of ATP synthase was negligible. The effects of ATP synthase inhibition by efrapeptin on the time course of the transmembrane electric potential (evaluated as carotenoid electrochromic response) and on ATP synthesis were studied comparatively. The ATP produced after a flash and the total charge associated with the proton flow coupled to ATP synthesis were more resistant to efrapeptin than the initial value of the phosphorylating currents, indicating that several ATP synthases are fed by protons from the same vesicle.

Met23Lys mutation in subunit gamma of F(O)F(1)-ATP synthase from Rhodobacter capsulatus impairs the activation of ATP hydrolysis by protonmotive force.

H(+)-F(O)F(1)-ATP synthase couples proton flow through its membrane portion, F(O), to the synthesis of ATP in its headpiece, F(1). Upon reversal of the reaction the enzyme functions as a proton pumping ATPase. Even in the simplest bacterial enzyme the ATPase activity is regulated by several mechanisms, involving inhibition by MgADP, conformational transitions of the epsilon subunit, and activation by protonmotive force. Here we report that the Met23Lys mutation in the gamma subunit of the Rhodobacter capsulatus ATP synthase significantly impaired the activation of ATP hydrolysis by protonmotive force. The impairment in the mutant was due to faster enzyme deactivation that was particularly evident at low ATP/ADP ratio. We suggest that the electrostatic interaction of the introduced gammaLys23 with the DELSEED region of subunit beta stabilized the ADP-inhibited state of the enzyme by hindering the rotation of subunit gamma rotation which is necessary for the activation.

Modulation of proton pumping efficiency in bacterial ATP synthases.

The ATP synthase in chromatophores of Rhodobacter caspulatus can effectively generate a transmembrane pH difference coupled to the hydrolysis of ATP. The rate of hydrolysis was rather insensitive to the depletion of ADP in the assay medium by an ATP regenerating system (phospho-enol-pyruvate (PEP) and pyruvate kinase (PK)). The steady state values of DeltapH were however drastically reduced as a consequence of ADP depletion. The clamped concentrations of ADP obtained using different PK activities in the assay medium could be calculated and an apparent Kd approximately 0.5 microM was estimated. The extent of proton uptake was also strongly dependent on the addition of phosphate to the assay medium. The Kd for this effect was about 70 microM. Analogous experiments were performed in membrane fragment from Escherichia coli. In this case, however, the hydrolysis rate was strongly inhibited by Pi, added up to 3 mM. Inhibition by Pi was nearly completely suppressed following depletion of ADP. The Kd's for the ADP and Pi were in the micromolar range and submillimolar range, respectively, and were mutually dependent from the concentration of the other ligand. Contrary to hydrolysis, the pumping of protons was rather insensitive to changes in the concentrations of the two ligands. At intermediate concentrations, proton pumping was actually stimulated, while the hydrolysis was inhibited. It is concluded that, in these two bacterial organisms, ADP and phosphate induce a functional state of the ATP synthase competent for a tightly coupled proton pumping, while the depletion of either one of these two ligands favors an inefficient (slipping) functional state. The switch between these states can probably be related to a structural change in the C-terminal alpha-helical hairpin of the epsilon-subunit, from an extended conformation, in which ATP hydrolysis is tightly coupled to proton pumping, to a retracted one, in which ATP hydrolysis and proton pumping are loosely coupled.

Physiological ligands ADP and Pi modulate the degree of intrinsic coupling in the ATP synthase of the photosynthetic bacterium Rhodobacter capsulatus.

The proton-pumping and the ATP hydrolysis activities of the ATP synthase of Rhodobacter capsulatus have been compared as a function of the ADP and P(i) concentrations. The proton pumping was measured either with the transmembrane pH difference probe, 9-amino-6-chloro-2-methoxyacridine, or with the transmembrane electric potential difference probe, bis(3-propyl-5-oxoisoxazol-4-yl)pentamethine oxonol, obtaining consistent results. The comparison indicates that an intrinsic uncoupling of ATP synthase is induced when the concentration of either ligand is decreased. The half-maximal effect was found in the submicromolar range for ADP and at about 70 microM for P(i). It is proposed that a switch from a partially uncoupled state of ATP synthase to the coupled state is induced by the simultaneous binding of ADP and P(i).