Defining the pathogenesis of the human Atp12p W94R mutation using a Saccharomyces cerevisiae yeast model

Predicting the effects of mutations on protein solubility using graph convolution network and protein language model representation

Solubility is one of the most important properties of protein. Protein solubility can be greatly changed by single amino acid mutations and the reduced protein solubility could lead to diseases. Since experimental methods to determine solubility are time-consuming and expensive, in-silico methods have been developed to predict the protein solubility changes caused by mutations mostly through protein evolution information. However, these methods are slow since it takes long time to obtain evolution information through multiple sequence alignment. In addition, these methods are of low performance because they do not fully utilize protein 3D structures due to a lack of experimental structures for most proteins. Here, we proposed a sequence-based method DeepMutSol to predict solubility change from residual mutations based on the Graph Convolutional Neural Network (GCN), where the protein graph was initiated according to predicted protein structure from Alphafold2, and the nodes (residues) were represented by protein language embeddings. To circumvent the small data of solubility changes, we further pretrained the model over absolute protein solubility. DeepMutSol was shown to outperform state-of-the-art methods in benchmark tests. In addition, we applied the method to clinically relevant genes from the ClinVar database and the predicted solubility changes were shown able to separate pathogenic mutations. All of the data sets and the source code are available at https://github.com/biomed-AI/DeepMutSol.

Gene co-expression network analysis reveals the positive impact of endocytosis and mitochondria-related genes over nitrogen metabolism in Saccharomyces cerevisiae

Nitrogen metabolism is essential for most cellular activities. Therefore, a deep understanding of its regulatory mechanisms is necessary for the efficient utilization of nitrogen sources for Saccharomyces cerevisiae. In this study, a gene co-expression network was constructed for S. cerevisiae S288C with different nitrogen sources. From this, a key gene co-expression module related to nitrogen source preference utilization was obtained, and 10 hub genes centrally located in the co-expression network were identified. Functional studies verified that the endocytosis-related genes CAP1 and END3 significantly increased the utilization of multiple non-preferred amino acids and reduced the accumulation of the harmful nitrogen metabolite precursor urea by regulating amino acid transporters and TOR pathway. The mitochondria-related gene ATP12, MRPL22, MRP1 and NAM9 significantly increased the utilization of multiple non-preferred amino acids and reduced accumulation of the urea by coordinately regulating nitrogen catabolism repression, Ssy1p-Ptr3p-Ssy5p signaling sensor system, amino acid transporters, TOR pathway and urea metabolism-related pathways. Furthermore, these data revealed the potential positive effects of endocytosis and mitochondrial ribosomes protein translation on nitrogen source preference. This study provides new analytical perspectives for complex regulatory networks involving nitrogen metabolism in S. cerevisiae.

PON-Sol2: Prediction of Effects of Variants on Protein Solubility

Genetic variations have a multitude of effects on proteins. A substantial number of variations affect protein-solvent interactions, either aggregation or solubility. Aggregation is often related to structural alterations, whereas solubilizable proteins in the solid phase can be made again soluble by dilution. Solubility is a central protein property and when reduced can lead to diseases. We developed a prediction method, PON-Sol2, to identify amino acid substitutions that increase, decrease, or have no effect on the protein solubility. The method is a machine learning tool utilizing gradient boosting algorithm and was trained on a large dataset of variants with different outcomes after the selection of features among a large number of tested properties. The method is fast and has high performance. The normalized correct prediction rate for three states is 0.656, and the normalized GC2 score is 0.312 in 10-fold cross-validation. The corresponding numbers in the blind test were 0.545 and 0.157. The performance was superior in comparison to previous methods. The PON-Sol2 predictor is freely available. It can be used to predict the solubility effects of variants for any organism, even in large-scale projects.

The ATP Synthase Deficiency in Human Diseases

Human diseases range from gene-associated to gene-non-associated disorders, including age-related diseases, neurodegenerative, neuromuscular, cardiovascular, diabetic diseases, neurocognitive disorders and cancer. Mitochondria participate to the cascades of pathogenic events leading to the onset and progression of these diseases independently of their association to mutations of genes encoding mitochondrial protein. Under physiological conditions, the mitochondrial ATP synthase provides the most energy of the cell via the oxidative phosphorylation. Alterations of oxidative phosphorylation mainly affect the tissues characterized by a high-energy metabolism, such as nervous, cardiac and skeletal muscle tissues. In this review, we focus on human diseases caused by altered expressions of ATP synthase genes of both mitochondrial and nuclear origin. Moreover, we describe the contribution of ATP synthase to the pathophysiological mechanisms of other human diseases such as cardiovascular, neurodegenerative diseases or neurocognitive disorders.

The Power of Yeast in Modelling Human Nuclear Mutations Associated with Mitochondrial Diseases

The increasing application of next generation sequencing approaches to the analysis of human exome and whole genome data has enabled the identification of novel variants and new genes involved in mitochondrial diseases. The ability of surviving in the absence of oxidative phosphorylation (OXPHOS) and mitochondrial genome makes the yeast Saccharomyces cerevisiae an excellent model system for investigating the role of these new variants in mitochondrial-related conditions and dissecting the molecular mechanisms associated with these diseases. The aim of this review was to highlight the main advantages offered by this model for the study of mitochondrial diseases, from the validation and characterisation of novel mutations to the dissection of the role played by genes in mitochondrial functionality and the discovery of potential therapeutic molecules. The review also provides a summary of the main contributions to the understanding of mitochondrial diseases emerged from the study of this simple eukaryotic organism.

Human Mitochondrial Pathologies of the Respiratory Chain and ATP Synthase: Contributions from Studies of Saccharomyces cerevisiae

The ease with which the unicellular yeast Saccharomyces cerevisiae can be manipulated genetically and biochemically has established this organism as a good model for the study of human mitochondrial diseases. The combined use of biochemical and molecular genetic tools has been instrumental in elucidating the functions of numerous yeast nuclear gene products with human homologs that affect a large number of metabolic and biological processes, including those housed in mitochondria. These include structural and catalytic subunits of enzymes and protein factors that impinge on the biogenesis of the respiratory chain. This article will review what is currently known about the genetics and clinical phenotypes of mitochondrial diseases of the respiratory chain and ATP synthase, with special emphasis on the contribution of information gained from pet mutants with mutations in nuclear genes that impair mitochondrial respiration. Our intent is to provide the yeast mitochondrial specialist with basic knowledge of human mitochondrial pathologies and the human specialist with information on how genes that directly and indirectly affect respiration were identified and characterized in yeast.

Homoplasmic deleterious MT-ATP6/8 mutations in adult patients

To address the frequency of complex V defects, we systematically sequenced MT-ATP6/8 genes in 512 consecutive patients. We performed functional analysis in muscle or fibroblasts for 12 out of 27 putative homoplasmic mutations and in cybrids for four. Fibroblasts, muscle and cybrids with known deleterious mutations underwent parallel analysis. It included oxidative phosphorylation spectrophotometric assays, western blots, structural analysis, ATP production, glycolysis and cell proliferation evaluation. We demonstrated the deleterious nature of three original mutations. Striking gradation in severity of the mutations consequences and differences between muscle, fibroblasts and cybrids implied a likely under-diagnosis of human complex V defects.

Modular biogenesis of mitochondrial respiratory complexes

Mitochondrial function relies on the activity of oxidative phosphorylation to synthesise ATP and generate an electrochemical gradient across the inner mitochondrial membrane. These coupled processes are mediated by five multi-subunit complexes that reside in this inner membrane. These complexes are the product of both nuclear and mitochondrial gene products. Defects in the function or assembly of these complexes can lead to mitochondrial diseases due to deficits in energy production and mitochondrial functions. Appropriate biogenesis and function are mediated by a complex number of assembly factors that promote maturation of specific complex subunits to form the active oxidative phosphorylation complex. The understanding of the biogenesis of each complex has been informed by studies in both simple eukaryotes such as Saccharomyces cerevisiae and human patients with mitochondrial diseases. These studies reveal each complex assembles through a pathway using specific subunits and assembly factors to form kinetically distinct but related assembly modules. The current understanding of these complexes has embraced the revolutions in genomics and proteomics to further our knowledge on the impact of mitochondrial biology in genetics, medicine, and evolution.

SODA: prediction of protein solubility from disorder and aggregation propensity

Solubility is an important, albeit not well understood, feature determining protein behavior. It is of paramount importance in protein engineering, where similar folded proteins may behave in very different ways in solution. Here we present SODA, a novel method to predict the changes of protein solubility based on several physico-chemical properties of the protein. SODA uses the propensity of the protein sequence to aggregate as well as intrinsic disorder, plus hydrophobicity and secondary structure preferences to estimate changes in solubility. It has been trained and benchmarked on two different datasets. The comparison to other recently published methods shows that SODA has state-of-the-art performance and is particularly well suited to predict mutations decreasing solubility. The method is fast, returning results for single mutations in seconds. A usage example estimating the full repertoire of mutations for a human germline antibody highlights several solubility hotspots on the surface. The web server, complete with RESTful interface and extensive help, can be accessed from URL: http://protein.bio.unipd.it/soda.

Biallelic Mutations in ATP5F1D, which Encodes a Subunit of ATP Synthase, Cause a Metabolic Disorder

ATP synthase, H+ transporting, mitochondrial F1 complex, δ subunit (ATP5F1D; formerly ATP5D) is a subunit of mitochondrial ATP synthase and plays an important role in coupling proton translocation and ATP production. Here, we describe two individuals, each with homozygous missense variants in ATP5F1D, who presented with episodic lethargy, metabolic acidosis, 3-methylglutaconic aciduria, and hyperammonemia. Subject 1, homozygous for c.245C>T (p.Pro82Leu), presented with recurrent metabolic decompensation starting in the neonatal period, and subject 2, homozygous for c.317T>G (p.Val106Gly), presented with acute encephalopathy in childhood. Cultured skin fibroblasts from these individuals exhibited impaired assembly of F1FO ATP synthase and subsequent reduced complex V activity. Cells from subject 1 also exhibited a significant decrease in mitochondrial cristae. Knockdown of Drosophila ATPsynδ, the ATP5F1D homolog, in developing eyes and brains caused a near complete loss of the fly head, a phenotype that was fully rescued by wild-type human ATP5F1D. In contrast, expression of the ATP5F1D c.245C>T and c.317T>G variants rescued the head-size phenotype but recapitulated the eye and antennae defects seen in other genetic models of mitochondrial oxidative phosphorylation deficiency. Our data establish c.245C>T (p.Pro82Leu) and c.317T>G (p.Val106Gly) in ATP5F1D as pathogenic variants leading to a Mendelian mitochondrial disease featuring episodic metabolic decompensation.

Insight into the Role of the Hv1 C-Terminal Domain in Dimer Stabilization

The voltage-gated proton-selective channel (Hv1) conducts protons in response to changes in membrane potential. The Hv1 protein forms dimers in the membrane. Crystal structures of Hv1 channels have revealed that the primary contacts between the two monomers are in the C-terminal domain (CTD), which forms a coiled-coil structure. The role of Hv1-CTD in channel assembly and activity is not fully understood. Here, molecular dynamics (MD) simulations of full-length and truncated CTD models of human and mouse Hv1 channels reveal a strong contribution of the CTD to the packing of the transmembrane domains. Simulations of the CTD models highlight four fundamental interactions of the key residues contributing to dimer stability. These include salt bridges, hydrophobic interactions, hydrogen bonds, and a disulfide bond across the dimer interface. At neutral pH, salt-bridge interactions increase dimer stability and the dimer becomes less stable at acidic pH. Hydrophobic core packing of the heptad pattern is important for stability, as shown by favorable nonpolar binding free energies rather than by electrostatic components. Moreover, free-energy calculations indicate that a more uniform hydrophobic core in the coiled-coil structure of the Hv1-NIN, a channel carrying the triple mutation M234N-N235I-V236N, leads to an increase in dimer stability with respect to the wild-type. A Cys disulfide bond has a strong impact on dimer stability by holding the dimer together and facilitating the interactions described above. These results are consistent with dissociative temperatures and energy barriers of dimer dissociation obtained from the temperature-accelerated MD.

PON-Sol: prediction of effects of amino acid substitutions on protein solubility

MOTIVATION

Solubility is one of the fundamental protein properties. It is of great interest because of its relevance to protein expression. Reduced solubility and protein aggregation are also associated with many diseases.

RESULTS

We collected from literature the largest experimentally verified solubility affecting amino acid substitution (AAS) dataset and used it to train a predictor called PON-Sol. The predictor can distinguish both solubility decreasing and increasing variants from those not affecting solubility. PON-Sol has normalized correct prediction ratio of 0.491 on cross-validation and 0.432 for independent test set. The performance of the method was compared both to solubility and aggregation predictors and found to be superior. PON-Sol can be used for the prediction of effects of disease-related substitutions, effects on heterologous recombinant protein expression and enhanced crystallizability. One application is to investigate effects of all possible AASs in a protein to aid protein engineering.

AVAILABILITY AND IMPLEMENTATION

PON-Sol is freely available at http://structure.bmc.lu.se/PON-Sol The training and test data are available at http://structure.bmc.lu.se/VariBench/ponsol.php

CONTACT

mauno.vihinen@med.lu.se

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Power(2): the power of yeast genetics applied to the powerhouse of the cell

The budding yeast Saccharomyces cerevisiae has served as a remarkable model organism for numerous seminal discoveries in biology. This paradigm extends to the mitochondria, a central hub for cellular metabolism, where studies in yeast have helped to reinvigorate the field and launch an exciting new era in mitochondrial biology. Here we discuss a few recent examples in which yeast research has laid a foundation for our understanding of evolutionarily conserved mitochondrial processes and functions, from key factors and pathways involved in the assembly of oxidative phosphorylation (OXPHOS) complexes to metabolite transport, lipid metabolism, and interorganelle communication. We also highlight new areas of yeast mitochondrial biology that are likely to aid in our understanding of the mitochondrial etiology of disease in the future.

Nuclear Genetic Defects of Mitochondrial ATP Synthase

Disorders of ATP synthase, the key enzyme of mitochondrial energy provision belong to the most severe metabolic diseases presenting as early-onset mitochondrial encephalo-cardiomyopathies. Up to now, mutations in four nuclear genes were associated with isolated deficiency of ATP synthase. Two of them, ATP5A1 and ATP5E encode enzyme’s structural subunits α and ε, respectively, while the other two ATPAF2 and TMEM70 encode specific ancillary factors that facilitate the biogenesis of ATP synthase. All these defects share a similar biochemical phenotype with pronounced decrease in the content of fully assembled and functional ATP synthase complex. However, substantial differences can be found in their frequency, molecular mechanism of pathogenesis, clinical manifestation as well as the course of the disease progression. While for TMEM70 the number of reported patients as well as spectrum of the mutations is steadily increasing, mutations in ATP5A1, ATP5E and ATPAF2 genes are very rare. Apparently, TMEM70 gene is highly prone to mutagenesis and this type of a rare mitochondrial disease has a rather frequent incidence. Here we present overview of individual reported cases of nuclear mutations in ATP synthase and discuss, how their analysis can improve our understanding of the enzyme biogenesis.

Mitochondrial ATP synthase: architecture, function and pathology

Human mitochondrial (mt) ATP synthase, or complex V consists of two functional domains: F(1), situated in the mitochondrial matrix, and F(o), located in the inner mitochondrial membrane. Complex V uses the energy created by the proton electrochemical gradient to phosphorylate ADP to ATP. This review covers the architecture, function and assembly of complex V. The role of complex V di-and oligomerization and its relation with mitochondrial morphology is discussed. Finally, pathology related to complex V deficiency and current therapeutic strategies are highlighted. Despite the huge progress in this research field over the past decades, questions remain to be answered regarding the structure of subunits, the function of the rotary nanomotor at a molecular level, and the human complex V assembly process. The elucidation of more nuclear genetic defects will guide physio(patho)logical studies, paving the way for future therapeutic interventions.