Regulation of Aerobic Energy Metabolism in Podospora anserina by Two Paralogous Genes Encoding Structurally Different c-Subunits of ATP Synthase.

Most of the ATP in living cells is produced by an F-type ATP synthase. This enzyme uses the energy of a transmembrane electrochemical proton gradient to synthesize ATP from ADP and inorganic phosphate. Proton movements across the membrane domain (FO) of the ATP synthase drive the rotation of a ring of 8-15 c-subunits, which induces conformational changes in the catalytic part (F1) of the enzyme that ultimately promote ATP synthesis. Two paralogous nuclear genes, called Atp9-5 and Atp9-7, encode structurally different c-subunits in the filamentous fungus Podospora anserina. We have in this study identified differences in the expression pattern for the two genes that correlate with the mitotic activity of cells in vegetative mycelia: Atp9-7 is transcriptionally active in non-proliferating (stationary) cells while Atp9-5 is expressed in the cells at the extremity (apex) of filaments that divide and are responsible for mycelium growth. When active, the Atp9-5 gene sustains a much higher rate of c-subunit synthesis than Atp9-7. We further show that the ATP9-7 and ATP9-5 proteins have antagonist effects on the longevity of P. anserina. Finally, we provide evidence that the ATP9-5 protein sustains a higher rate of mitochondrial ATP synthesis and yield in ATP molecules per electron transferred to oxygen than the c-subunit encoded by Atp9-7. These findings reveal that the c-subunit genes play a key role in the modulation of ATP synthase production and activity along the life cycle of P. anserina. Such a degree of sophistication for regulating aerobic energy metabolism has not been described before.

Multidimensional mutual information methods for the analysis of covariation in multiple sequence alignments.

BACKGROUND: Several methods are available for the detection of covarying positions from a multiple sequence alignment (MSA). If the MSA contains a large number of sequences, information about the proximities between residues derived from covariation maps can be sufficient to predict a protein fold. However, in many cases the structure is already known, and information on the covarying positions can be valuable to understand the protein mechanism and dynamic properties. RESULTS: In this study we have sought to determine whether a multivariate (multidimensional) extension of traditional mutual information (MI) can be an additional tool to study covariation. The performance of two multidimensional MI (mdMI) methods, designed to remove the effect of ternary/quaternary interdependencies, was tested with a set of 9 MSAs each containing <400 sequences, and was shown to be comparable to that of the newest methods based on maximum entropy/pseudolikelyhood statistical models of protein sequences. However, while all the methods tested detected a similar number of covarying pairs among the residues separated by < 8 Å in the reference X-ray structures, there was on average less than 65% overlap between the top scoring pairs detected by methods that are based on different principles. CONCLUSIONS: Given the large variety of structure and evolutionary history of different proteins it is possible that a single best method to detect covariation in all proteins does not exist, and that for each protein family the best information can be derived by merging/comparing results obtained with different methods. This approach may be particularly valuable in those cases in which the size of the MSA is small or the quality of the alignment is low, leading to significant differences in the pairs detected by different methods.

Assembly-dependent translation of subunits 6 (Atp6) and 9 (Atp9) of ATP synthase in yeast mitochondria.

The yeast mitochondrial ATP synthase is an assembly of 28 subunits of 17 types of which 3 (subunits 6, 8, and 9) are encoded by mitochondrial genes, while the 14 others have a nuclear genetic origin. Within the membrane domain (FO) of this enzyme, the subunit 6 and a ring of 10 identical subunits 9 transport protons across the mitochondrial inner membrane coupled to ATP synthesis in the extra-membrane structure (F1) of ATP synthase. As a result of their dual genetic origin, the ATP synthase subunits are synthesized in the cytosol and inside the mitochondrion. How they are produced in the proper stoichiometry from two different cellular compartments is still poorly understood. The experiments herein reported show that the rate of translation of the subunits 9 and 6 is enhanced in strains with mutations leading to specific defects in the assembly of these proteins. These translation modifications involve assembly intermediates interacting with subunits 6 and 9 within the final enzyme and cis-regulatory sequences that control gene expression in the organelle. In addition to enabling a balanced output of the ATP synthase subunits, these assembly-dependent feedback loops are presumably important to limit the accumulation of harmful assembly intermediates that have the potential to dissipate the mitochondrial membrane electrical potential and the main source of chemical energy of the cell.

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

Studies in yeast have shown that a deficiency in Atp12p prevents assembly of the extrinsic domain (F(1)) of complex V and renders cells unable to make ATP through oxidative phosphorylation. De Meirleir et al. (De Meirleir, L., Seneca, S., Lissens, W., De Clercq, I., Eyskens, F., Gerlo, E., Smet, J., and Van Coster, R. (2004) J. Med. Genet. 41, 120-124) have reported that a homozygous missense mutation in the gene for human Atp12p (HuAtp12p), which replaces Trp-94 with Arg, was linked to the death of a 14-month-old patient. We have investigated the impact of the pathogenic W94R mutation on Atp12p structure/function. Plasmid-borne wild type human Atp12p rescues the respiratory defect of a yeast ATP12 deletion mutant (Deltaatp12). The W94R mutation alters the protein at the most highly conserved position in the Pfam sequence and renders HuAtp12p insoluble in the background of Deltaatp12. In contrast, the yeast protein harboring the corresponding mutation, ScAtp12p(W103R), is soluble in the background of Deltaatp12 but not in the background of Deltaatp12Deltafmc1, a strain that also lacks Fmc1p. Fmc1p is a yeast mitochondrial protein not found in higher eukaryotes. Tryptophan 94 (human) or 103 (yeast) is located in a positively charged region of Atp12p, and hence its mutation to arginine does not alter significantly the electrostatic properties of the protein. Instead, we provide evidence that the primary effect of the substitution is on the dynamic properties of Atp12p.

Consequences of the pathogenic T9176C mutation of human mitochondrial DNA on yeast mitochondrial ATP synthase.

Several human neurological disorders have been associated with various mutations affecting mitochondrial enzymes involved in cellular ATP production. One of these mutations, T9176C in the mitochondrial DNA (mtDNA), changes a highly conserved leucine residue into proline at position 217 of the mitochondrially encoded Atp6p (or a) subunit of the F1FO-ATP synthase. The consequences of this mutation on the mitochondrial ATP synthase are still poorly defined. To gain insight into the primary pathogenic mechanisms induced by T9176C, we have investigated the consequences of this mutation on the ATP synthase of yeast where Atp6p is also encoded by the mtDNA. In vitro, yeast atp6-T9176C mitochondria showed a 30% decrease in the rate of ATP synthesis. When forcing the F1FO complex to work in the reverse mode, i.e. F1-catalyzed hydrolysis of ATP coupled to proton transport out of the mitochondrial matrix, the mutant showed a normal proton-pumping activity and this activity was fully sensitive to oligomycin, an inhibitor of the ATP synthase proton channel. However, under conditions of maximal ATP hydrolytic activity, using non-osmotically protected mitochondria, the mutant ATPase activity was less efficiently inhibited by oligomycin (60% inhibition versus 85% for the wild type control). Blue Native Polyacrylamide Gel Electrophoresis analyses revealed that atp6-T9176C yeast accumulated rather good levels of fully assembled ATP synthase complexes. However, a number of sub-complexes (F1, Atp9p-ring, unassembled alpha-F1 subunits) could be detected as well, presumably because of a decreased stability of Atp6p within the ATP synthase. Although the oxidative phosphorylation capacity was reduced in atp6-T9176C yeast, the number of ATP molecules synthesized per electron transferred to oxygen was similar compared with wild type yeast. It can therefore be inferred that the coupling efficiency within the ATP synthase was mostly unaffected and that the T9176C mutation did not increase the proton permeability of the mitochondrial inner membrane.

Chaperones of F1-ATPase.

Mitochondrial F(1)-ATPase contains a hexamer of alternating alpha and beta subunits. The assembly of this structure requires two specialized chaperones, Atp11p and Atp12p, that bind transiently to beta and alpha. In the absence of Atp11p and Atp12p, the hexamer is not formed, and alpha and beta precipitate as large insoluble aggregates. An early model for the mechanism of chaperone-mediated F(1) assembly (Wang, Z. G., Sheluho, D., Gatti, D. L., and Ackerman, S. H. (2000) EMBO J. 19, 1486-1493) hypothesized that the chaperones themselves look very much like the alpha and beta subunits, and proposed an exchange of Atp11p for alpha and of Atp12p for beta; the driving force for the exchange was expected to be a higher affinity of alpha and beta for each other than for the respective chaperone partners. One important feature of this model was the prediction that as long as Atp11p is bound to beta and Atp12p is bound to alpha, the two F(1) subunits cannot interact at either the catalytic site or the noncatalytic site interface. Here we present the structures of Atp11p from Candida glabrata and Atp12p from Paracoccus denitrificans, and we show that some features of the Wang model are correct, namely that binding of the chaperones to alpha and beta prevents further interactions between these F(1) subunits. However, Atp11p and Atp12p do not resemble alpha or beta, and it is instead the F(1) gamma subunit that initiates the release of the chaperones from alpha and beta and their further assembly into the mature complex.

Methods to determine the status of mitochondrial ATP synthase assembly.

The adenosine triphosphate (ATP) synthase (F1-F0 complex) of the mitochondrial inner membrane is responsible for making nearly all of the ATP utilized by eukaryotic organisms. The enzyme is an oligomer of more than 20 different subunits, 14 of which are essential for its catalytic activity. The other subunits function in the regulation and structure of the complex. Subunits essential for catalytic activity make up the proton pore, the bulk of the F1 headpiece, and the two stalks that physically and functionally couple the catalytic and proton-translocating activities of the ATP synthase. Saccharomyces cerevisiae provides an excellent model system for studying mutations that affect assembly of the complex because of the ability of this organism to survive on the ATP produced from fermentation in the absence of mitochondrial respiration or oxidative phosphorylation. Studies of such mutants have been instrumental in identifying novel molecular chaperones that act at discrete steps of F1-F0 assembly. Here, we describe some experimental approaches useful in assessing the status of F1-F0 assembly.

Atp11p and Atp12p are chaperones for F(1)-ATPase biogenesis in mitochondria.

The bioenergetic needs of aerobic cells are met principally through the action of the F(1)F(0) ATP synthase, which catalyzes ATP synthesis during oxidative phosphorylation. The catalytic unit of the enzyme (F(1)) is a multimeric protein of the subunit composition alpha(3)beta(3)(gamma)(delta) epsilon. Our work, which employs the yeast Saccharomyces cerevisiae as a model system for studies of mitochondrial function, has provided evidence that assembly of the mitochondrial alpha and beta subunits into the F(1) oligomer requires two molecular chaperone proteins called Atp11p and Atp12p. Comprehensive knowledge of Atp11p and Atp12p activities in mitochondria bears relevance to human physiology and disease as these chaperone actions are now known to exist in mitochondria of human cells.

Defining the impact on yeast ATP synthase of two pathogenic human mitochondrial DNA mutations, T9185C and T9191C.

Mutations in the human mitochondrial ATP6 gene encoding ATP synthase subunit a/6 (referred to as Atp6p in yeast) are at the base of neurodegenerative disorders like Neurogenic Ataxia and Retinitis Pigmentosa (NARP), Leigh syndrome (LS), Charcot-Marie-Tooth (CMT), and ataxia telangiectasia. In previous studies, using the yeast Saccharomyces cerevisiae as a model we were able to better define how several of these mutations impact the ATP synthase. Here we report the construction of yeast models of two other ATP6 pathogenic mutations, T9185C and T9191C. The first one was reported as conferring a mild, sometimes reversible, CMT clinical phenotype; the second one has been described in a patient presenting with severe LS. We found that an equivalent of the T9185C mutation partially impaired the functioning of yeast ATP synthase, with only a 30% deficit in mitochondrial ATP production. An equivalent of the mutation T9191C had much more severe effects, with a nearly complete block in yeast Atp6p assembly and an >95% drop in the rate of ATP synthesis. These findings provide a molecular basis for the relative severities of the diseases induced by T9185C and T9191C.

The depletion of F1 subunit ε in yeast leads to an uncoupled respiratory phenotype that is rescued by mutations in the proton-translocating subunits of F0.

The central stalk of the ATP synthase is an elongated hetero-oligomeric structure providing a physical connection between the catalytic sites in F1 and the proton translocation channel in F0 for energy transduction between the two subdomains. The shape of the central stalk and relevance to energy coupling are essentially the same in ATP synthases from all forms of life, yet the protein composition of this domain changed during evolution of the mitochondrial enzyme from a two- to a three-subunit structure (γ, δ, ε). Whereas the mitochondrial γ- and δ-subunits are homologues of the bacterial central stalk proteins, the deliberate addition of subunit ε is poorly understood. Here we report that down-regulation of the gene (ATP15) encoding the ε-subunit rapidly leads to lethal F0-mediated proton leaks through the membrane because of the loss of stability of the ATP synthase. The ε-subunit is thus essential for oxidative phosphorylation. Moreover, mutations in F0 subunits a and c, which slow the proton translocation rate, are identified that prevent ε-deficient ATP synthases from dissipating the electrochemical potential. Cumulatively our data lead us to propose that the ε-subunit evolved to permit operation of the central stalk under the torque imposed at the normal speed of proton movement through mitochondrial F0.

Biapenem inactivation by B2 metallo ß-lactamases: energy landscape of the hydrolysis reaction.

BACKGROUND: A general mechanism has been proposed for metallo ß-lactamases (MßLs), in which deprotonation of a water molecule near the Zn ion(s) results in the formation of a hydroxide ion that attacks the carbonyl oxygen of the ß-lactam ring. However, because of the absence of X-ray structures that show the exact position of the antibiotic in the reactant state (RS) it has been difficult to obtain a definitive validation of this mechanism. METHODOLOGY/PRINCIPAL FINDINGS: We have employed a strategy to identify the RS, which does not rely on substrate docking and/or molecular dynamics. Starting from the X-ray structure of the enzyme:product complex (the product state, PS), a QM/MM scan was used to drive the reaction uphill from product back to reactant. Since in this process also the enzyme changes from PS to RS, we actually generate the enzyme:substrate complex from product and avoid the uncertainties associated with models of the reactant state. We used this strategy to study the reaction of biapenem hydrolysis by B2 MßL CphA. QM/MM simulations were carried out under 14 different ionization states of the active site, in order to generate potential energy surfaces (PESs) corresponding to a variety of possible reaction paths. CONCLUSIONS/SIGNIFICANCE: The calculations support a model for biapenem hydrolysis by CphA, in which the nucleophile that attacks the ß-lactam ring is not the water molecule located in proximity of the active site Zn, but a second water molecule, hydrogen bonded to the first one, which is used up in the reaction, and thus is not visible in the X-ray structure of the enzyme:product complex.

Accurate simulation and detection of coevolution signals in multiple sequence alignments.

BACKGROUND: While the conserved positions of a multiple sequence alignment (MSA) are clearly of interest, non-conserved positions can also be important because, for example, destabilizing effects at one position can be compensated by stabilizing effects at another position. Different methods have been developed to recognize the evolutionary relationship between amino acid sites, and to disentangle functional/structural dependencies from historical/phylogenetic ones. METHODOLOGY/PRINCIPAL FINDINGS: We have used two complementary approaches to test the efficacy of these methods. In the first approach, we have used a new program, MSAvolve, for the in silico evolution of MSAs, which records a detailed history of all covarying positions, and builds a global coevolution matrix as the accumulated sum of individual matrices for the positions forced to co-vary, the recombinant coevolution, and the stochastic coevolution. We have simulated over 1600 MSAs for 8 protein families, which reflect sequences of different sizes and proteins with widely different functions. The calculated coevolution matrices were compared with the coevolution matrices obtained for the same evolved MSAs with different coevolution detection methods. In a second approach we have evaluated the capacity of the different methods to predict close contacts in the representative X-ray structures of an additional 150 protein families using only experimental MSAs. CONCLUSIONS/SIGNIFICANCE: Methods based on the identification of global correlations between pairs were found to be generally superior to methods based only on local correlations in their capacity to identify coevolving residues using either simulated or experimental MSAs. However, the significant variability in the performance of different methods with different proteins suggests that the simulation of MSAs that replicate the statistical properties of the experimental MSA can be a valuable tool to identify the coevolution detection method that is most effective in each case.

The contribution of coevolving residues to the stability of KDO8P synthase.

BACKGROUND: The evolutionary tree of 3-deoxy-D-manno-octulosonate 8-phosphate (KDO8P) synthase (KDO8PS), a bacterial enzyme that catalyzes a key step in the biosynthesis of bacterial endotoxin, is evenly divided between metal and non-metal forms, both having similar structures, but diverging in various degrees in amino acid sequence. Mutagenesis, crystallographic and computational studies have established that only a few residues determine whether or not KDO8PS requires a metal for function. The remaining divergence in the amino acid sequence of KDO8PSs is apparently unrelated to the underlying catalytic mechanism. METHODOLOGY/PRINCIPAL FINDINGS: The multiple alignment of all known KDO8PS sequences reveals that several residue pairs coevolved, an indication of their possible linkage to a structural constraint. In this study we investigated by computational means the contribution of coevolving residues to the stability of KDO8PS. We found that about 1/4 of all strongly coevolving pairs probably originated from cycles of mutation (decreasing stability) and suppression (restoring it), while the remaining pairs are best explained by a succession of neutral or nearly neutral covarions. CONCLUSIONS/SIGNIFICANCE: Both sequence conservation and coevolution are involved in the preservation of the core structure of KDO8PS, but the contribution of coevolving residues is, in proportion, smaller. This is because small stability gains or losses associated with selection of certain residues in some regions of the stability landscape of KDO8PS are easily offset by a large number of possible changes in other regions. While this effect increases the tolerance of KDO8PS to deleterious mutations, it also decreases the probability that specific pairs of residues could have a strong contribution to the thermodynamic stability of the protein.

Failure to assemble the alpha 3 beta 3 subcomplex of the ATP synthase leads to accumulation of the alpha and beta subunits within inclusion bodies and the loss of mitochondrial cristae in Saccharomyces cerevisiae.

The F(1) component of mitochondrial ATP synthase is an oligomeric assembly of five different subunits, alpha, beta, gamma, delta, and epsilon. In terms of mass, the bulk of the structure ( approximately 90%) is provided by the alpha and beta subunits, which form an (alphabeta)(3) hexamer with adenine nucleotide binding sites at the alpha/beta interfaces. We report here ultrastructural and immunocytochemical analyses of yeast mutants that are unable to form the alpha(3)beta(3) oligomer, either because the alpha or the beta subunit is missing or because the cells are deficient for proteins that mediate F assembly (e.g. Atp11p, Atp12p, or Fmc1p). The F(1) alpha(1) and beta subunits of such mutant strains are detected within large electron-dense particles in the mitochondrial matrix. The composition of the aggregated species is principally full-length F(1) alpha and/or beta subunit protein that has been processed to remove the amino-terminal targeting peptide. To our knowledge this is the first demonstration of mitochondrial inclusion bodies that are formed largely of one particular protein species. We also show that yeast mutants lacking the alpha(3)beta(3) oligomer are devoid of mitochondrial cristae and are severely deficient for respiratory complexes III and IV. These observations are in accord with other studies in the literature that have pointed to a central role for the ATP synthase in biogenesis of the mitochondrial inner membrane.

The molecular chaperone, Atp12p, from Homo sapiens. In vitro studies with purified wild type and mutant (E240K) proteins.

Work in Saccharomyces cerevisiae has shown that Atp12p binds to unassembled alpha subunits of F(1) and in so doing prevents the alpha subunit from associating with itself in non-productive complexes during assembly of the F(1) moiety of the mitochondrial ATP synthase. We have developed a method to prepare recombinant Atp12p after expression of its human cDNA in bacterial cells. The molecular chaperone activity of HuAtp12p was studied using citrate synthase as a model substrate. Wild type HuAtp12p suppresses the aggregation of thermally inactivated citrate synthase. In contrast, the mutant protein HuAtp12p(E240K), which harbors a lysine at the position of the highly conserved Glu-240, fails to prevent citrate synthase aggregation at 43 degrees C. No significant differences were observed between the wild type and the mutant proteins as judged by sedimentation analysis, cysteine titration, tryptophan emission spectra, or limited proteolysis, which suggests that the E240K mutation alters the activity of HuAtp12p with minimal effects on the physical integrity of the protein. An additional important finding of this work is that the equilibrium chemical denaturation curve of HuAtp12p shows two components, the first of which is associated with protein aggregation. This result is consistent with a model for Atp12p structure in which there is a hydrophobic chaperone domain that is buried within the protein interior.

A purified subfragment of yeast Atp11p retains full molecular chaperone activity.

Atp11p is a molecular chaperone of the mitochondrial matrix that participates in the biogenesis pathway to form F1, the catalytic unit of the ATP synthase. Affinity tag pull-down assays and yeast two-hybrid screens have shown that Atp11p binds to free beta subunits of F1 (Wang, Z. G., and Ackerman, S. H. (2000) J. Biol. Chem. 275, 5767-5772). This binding action prevents the beta subunit from associating with itself in non-productive complexes and fosters the formation of a (alpha beta)3 hexamer. Following the premise that Atp11p action is mediated primarily through a surface (as opposed to specific amino acids, as in an enzyme active site), solving its three-dimensional structure so that we may learn how the shape of the protein influences its function is a high priority. Recombinant yeast Atp11p has proven refractory for such analysis because of the presence of a disordered region in the protein. In this article, we show that removal of 67 residues from the amino terminus of recombinant Atp11p yields a subfragment of the protein (called Atp11pTRNC) that retains molecular chaperone function as determined in vitro with both a surrogate substrate (reduced insulin) and the natural substrate (F1 beta). Moreover, preliminary 15N-1H heteronuclear single quantum coherence spectra obtained with Atp11pTRNC indicate that the truncated protein is well ordered and amenable to structure determination by nuclear magnetic resonance.