Identification of important amino acid residues that modulate binding of Escherichia coli GroEL to its various cochaperones

DeMAG predicts the effects of variants in clinically actionable genes by integrating structural and evolutionary epistatic features

Despite the increasing use of genomic sequencing in clinical practice, the interpretation of rare genetic variants remains challenging even in well-studied disease genes, resulting in many patients with Variants of Uncertain Significance (VUSs). Computational Variant Effect Predictors (VEPs) provide valuable evidence in variant assessment, but they are prone to misclassifying benign variants, contributing to false positives. Here, we develop Deciphering Mutations in Actionable Genes (DeMAG), a supervised classifier for missense variants trained using extensive diagnostic data available in 59 actionable disease genes (American College of Medical Genetics and Genomics Secondary Findings v2.0, ACMG SF v2.0). DeMAG improves performance over existing VEPs by reaching balanced specificity (82%) and sensitivity (94%) on clinical data, and includes a novel epistatic feature, the 'partners score', which leverages evolutionary and structural partnerships of residues. The 'partners score' provides a general framework for modeling epistatic interactions, integrating both clinical and functional information. We provide our tool and predictions for all missense variants in 316 clinically actionable disease genes (demag.org) to facilitate the interpretation of variants and improve clinical decision-making.

Ring Separation Highlights the Protein-Folding Mechanism Used by the Phage EL-Encoded Chaperonin

Chaperonins are ubiquitous, ATP-dependent protein-folding molecular machines that are essential for all forms of life. Bacteriophage φEL encodes its own chaperonin to presumably fold exceedingly large viral proteins via profoundly different nucleotide-binding conformations. Our structural investigations indicate that ATP likely binds to both rings simultaneously and that a misfolded substrate acts as the trigger for ATP hydrolysis. More importantly, the φEL complex dissociates into two single rings resulting from an evolutionarily altered residue in the highly conserved ATP-binding pocket. Conformational changes also more than double the volume of the single-ring internal chamber such that larger viral proteins are accommodated. This is illustrated by the fact that φEL is capable of folding β-galactosidase, a 116-kDa protein. Collectively, the architecture and protein-folding mechanism of the φEL chaperonin are significantly different from those observed in group I and II chaperonins.

Allosteric Mechanisms in Chaperonin Machines

Chaperonins are nanomachines that facilitate protein folding by undergoing energy (ATP)-dependent movements that are coordinated in time and space owing to complex allosteric regulation. They consist of two back-to-back stacked oligomeric rings with a cavity at each end where protein substrate folding can take place. Here, we focus on the GroEL/GroES chaperonin system from Escherichia coli and, to a lesser extent, on the more poorly characterized eukaryotic chaperonin CCT/TRiC. We describe their various functional (allosteric) states and how they are affected by substrates and allosteric effectors that include ATP, ADP, nonfolded protein substrates, potassium ions, and GroES (in the case of GroEL). We also discuss the pathways of intra- and inter-ring allosteric communication by which they interconvert and the coupling between allosteric transitions and protein folding reactions.

GroEL2 of Mycobacterium tuberculosis Reveals the Importance of Structural Pliability in Chaperonin Function

Intracellular protein folding is mediated by molecular chaperones, the best studied among which are the chaperonins GroEL and GroES. Conformational changes and allosteric transitions between different metastable states are hallmarks of the chaperonin mechanism. These conformational transitions between three structural domains of GroEL are anchored at two hinges. Although hinges are known to be critical for mediating the communication between different domains of GroEL, the relative importance of hinges on GroEL oligomeric assembly, ATPase activity, conformational changes, and functional activity is not fully characterized. We have exploited the inability of Mycobacterium tuberculosis GroEL2 to functionally complement an Escherichia coli groEL mutant to address the importance of hinge residues in the GroEL mechanism. Various chimeras of M. tuberculosis GroEL2 and E. coli GroEL allowed us to understand the role of hinges and dissect the consequences of oligomerization and substrate binding capability on conformational transitions. The present study explains the concomitant conformational changes observed with GroEL hinge variants and is best supported by the normal mode analysis.

IMPORTANCE

Conformational changes and allosteric transitions are hallmarks of the chaperonin mechanism. We have exploited the inability of M. tuberculosis GroEL2 to functionally complement a strain of E. coli in which groEL expression is repressed to address the importance of hinges. The significance of conservation at the hinge regions stands out as a prominent feature of the GroEL mechanism in binding to GroES and substrate polypeptides. The hinge residues play a significant role in the chaperonin activity in vivo and in vitro.

G673 could be a novel mutational hot spot for intragenic suppressors of pheS5 lesion in Escherichia coli

The pheS5 Ts mutant of Escherichia coli defined by a G293 → A293 transition, which is responsible for thermosensitive Phenylalanyl-tRNA synthetase has been well studied at both biochemical and molecular level but genetic analyses pertaining to suppressors of pheS5 were hard to come by. Here we have systematically analyzed a spectrum of Temperature-insensitive derivatives isolated from pheS5 Ts mutant and identified two intragenic suppressors affecting the same base pair coordinate G673 (pheS19 defines G673 → T673 ; Gly225 → Cys225 and pheS28 defines G673 → C673 ; Gly225 → Arg225). In fact in the third derivative, the intragenic suppressor originally named pheS43 (G673 → C673 transversion) is virtually same as pheS28. In the fourth case, the very pheS5 lesion itself has got changed from A293 → T293 (named pheS40). Cloning of pheS(+), pheS5, pheS5-pheS19, pheS5-pheS28 alleles into pBR322 and introduction of these clones into pheS5 mutant revealed that excess of double mutant protein is not at all good for the survival of cells at 42°C. These results clearly indicate a pivotal role for Gly225 in the structural/functional integrity of alpha subunit of E. coli PheRS enzyme and it is proposed that G673 might define a hot spot for intragenic suppressors of pheS5.

Identification of elements that dictate the specificity of mitochondrial Hsp60 for its co-chaperonin

Type I chaperonins (cpn60/Hsp60) are essential proteins that mediate the folding of proteins in bacteria, chloroplast and mitochondria. Despite the high sequence homology among chaperonins, the mitochondrial chaperonin system has developed unique properties that distinguish it from the widely-studied bacterial system (GroEL and GroES). The most relevant difference to this study is that mitochondrial chaperonins are able to refold denatured proteins only with the assistance of the mitochondrial co-chaperonin. This is in contrast to the bacterial chaperonin, which is able to function with the help of co-chaperonin from any source. The goal of our work was to determine structural elements that govern the specificity between chaperonin and co-chaperonin pairs using mitochondrial Hsp60 as model system. We used a mutagenesis approach to obtain human mitochondrial Hsp60 mutants that are able to function with the bacterial co-chaperonin, GroES. We isolated two mutants, a single mutant (E321K) and a double mutant (R264K/E358K) that, together with GroES, were able to rescue an E. coli strain, in which the endogenous chaperonin system was silenced. Although the mutations are located in the apical domain of the chaperonin, where the interaction with co-chaperonin takes place, none of the residues are located in positions that are directly responsible for co-chaperonin binding. Moreover, while both mutants were able to function with GroES, they showed distinct functional and structural properties. Our results indicate that the phenotype of the E321K mutant is caused mainly by a profound increase in the binding affinity to all co-chaperonins, while the phenotype of R264K/E358K is caused by a slight increase in affinity toward co-chaperonins that is accompanied by an alteration in the allosteric signal transmitted upon nucleotide binding. The latter changes lead to a great increase in affinity for GroES, with only a minor increase in affinity toward the mammalian mitochondrial co-chaperonin.

An ORFan no more: the bacteriophage T4 39.2 gene product, NwgI, modulates GroEL chaperone function

Bacteriophages are the most abundant biological entities in our biosphere, characterized by their hyperplasticity, mosaic composition, and the many unknown functions (ORFans) encoded by their immense genetic repertoire. These genes are potentially maintained by the bacteriophage to allow efficient propagation on hosts encountered in nature. To test this hypothesis, we devised a selection to identify bacteriophage-encoded gene(s) that modulate the host Escherichia coli GroEL/GroES chaperone machine, which is essential for the folding of certain host and bacteriophage proteins. As a result, we identified the bacteriophage RB69 gene 39.2, of previously unknown function and showed that homologs of 39.2 in bacteriophages T4, RB43, and RB49 similarly modulate GroEL/GroES. Production of wild-type bacteriophage T4 Gp39.2, a 58-amino-acid protein, (a) enables diverse bacteriophages to plaque on the otherwise nonpermissive groES or groEL mutant hosts in an allele-specific manner, (b) suppresses the temperature-sensitive phenotype of both groES and groEL mutants, (c) suppresses the defective UV-induced PolV function (UmuCD) of the groEL44 mutant, and (d) is lethal to the host when overproduced. Finally, as proof of principle that Gp39.2 is essential for bacteriophage growth on certain bacterial hosts, we constructed a T4 39.2 deletion strain and showed that, unlike the isogenic wild-type parent, it is incapable of propagating on certain groEL mutant hosts. We propose a model of how Gp39.2 modulates GroES/GroEL function.

Manipulation of conformational change in proteins by single-residue perturbations

Using the perturbation-response scanning (PRS) technique, we study a set of 25 proteins that display a variety of conformational motions upon ligand binding (e.g., shear, hinge, allosteric). In most cases, PRS determines single residues that may be manipulated to achieve the resulting conformational change. PRS reveals that for some proteins, binding-induced conformational change may be achieved through the perturbation of residues scattered throughout the protein, whereas in others, perturbation of specific residues confined to a highly specific region is necessary. Overlaps between the experimental and PRS-calculated atomic displacement vectors are usually more descriptive of the conformational change than those obtained from a modal analysis of elastic network models. Furthermore, the largest overlaps obtained by the latter approach do not always appear in the most collective modes; there are cases where more than one mode yields comparable overlap sizes. We show that success of the modal analysis depends on an absence of redundant paths in the protein. PRS thus demonstrates that several relevant modes can be induced simultaneously by perturbing a single select residue on the protein. We also illustrate the biological relevance of applying PRS to the GroEL, adenylate kinase, myosin, and kinesin structures in detail by showing that the residues whose perturbation leads to precise conformational changes usually correspond to those experimentally determined to be functionally important.

Epistatic interactions modulate the evolution of mammalian mitochondrial respiratory complex components

BACKGROUND

The deleterious effect of a mutation can be reverted by a second-site interacting residue. This is an epistatic compensatory process explaining why mutations that are deleterious in some species are tolerated in phylogenetically related lineages, rendering evident that those mutations are, by all means, only deleterious in the species-specific context. Although an extensive and refined theoretical framework on compensatory evolution does exist, the supporting evidence remains limited, especially for protein models. In this current study, we focused on the molecular mechanism underlying the epistatic compensatory process in mammalian mitochondrial OXPHOS proteins using a combination of in-depth structural and sequence analyses.

RESULTS

Modeled human structures were used in this study to predict the structural impairment and recovery of deleterious mutations alone and combined with an interacting compensatory partner, respectively. In two cases, COI and COIII, intramolecular interactions between spatially linked residues restore the folding pattern impaired by the deleterious mutation. In a third case, intermolecular contact between mitochondrial CYB and nuclear CYT1 encoded components of the cytochrome bc1 complex are likely to restore protein binding. Moreover, we observed different modes of compensatory evolution that have resulted in either a quasi-simultaneous occurrence of a mutation and corresponding compensatory partner, or in independent occurrences of mutations in distinct lineages that were always preceded by the compensatory site.

CONCLUSION

Epistatic interactions between individual replacements involving deleterious mutations seems to follow a parsimonious model of evolution in which genomes hold pre-compensating states that subsequently tolerate deleterious mutations. This phenomenon is likely to have been constraining the variability at coevolving sites and shaping the interaction between the mitochondrial and the nuclear genome.

Conserved inserts in the Hsp60 (GroEL) and Hsp70 (DnaK) proteins are essential for cellular growth

The Hsp60 and Hsp70 chaperones contain a number of conserved inserts that are restricted to particular phyla of bacteria. A one aa insert in the E. coli GroEL and a 21-23 insert in the DnaK proteins are specific for most Gram-negative bacteria. Two other inserts in DnaK are limited to certain groups of proteobacteria. The requirement of these inserts for cellular growth was examined by carrying out complementation studies with temperature-sensitive (T(s)) mutants of E. coli groEL or dnaK. Our results demonstrate that deletion or most changes in these inserts completely abolished the complementation ability of the mutant proteins. Studies with GroEL and DnaK from some other species that either lacked or contained these inserts also indicated that these inserts are essential for growth of E. coli. The DnaK from some bacteria contains a two aa insert that is not found in E. coli. Introduction of this insert into the E. coli DnaK also led to its inactivation, indicating that these inserts are specific for different groups. We postulate that these conserved inserts that are localized in loop regions on protein surfaces, are involved in some ancillary functions that are essential for the groups of bacteria where they are found.

Allostery wiring diagrams in the transitions that drive the GroEL reaction cycle

Determining the network of residues that transmit allosteric signals is crucial to understanding the function of biological nanomachines. During the course of a reaction cycle, biological machines in general, and Escherichia coli chaperonin GroEL in particular, undergo large-scale conformational changes in response to ligand binding. Normal mode analyses, based on structure-based coarse-grained models where each residue is represented by an alpha carbon atom, have been widely used to describe the motions encoded in the structures of proteins. Here, we propose a new Calpha-side chain elastic network model of proteins that includes information about the physical identity of each residue and accurately accounts for the side-chain topology and packing within the structure. Using the Calpha-side chain elastic network model and the structural perturbation method, which probes the response of a local perturbation at a given site at all other sites in the structure, we determine the network of key residues (allostery wiring diagram) responsible for the T-->R and R''-->T transitions in GroEL. A number of residues, both within a subunit and at the interface of two adjacent subunits, are found to be at the origin of the positive cooperativity in the ATP-driven T-->R transition. Of particular note are residues G244, R58, D83, E209, and K327. Of these, R38, D83, and K327 are highly conserved. G244 is located in the apical domain at the interface between two subunits; E209 and K327 are located in the apical domain, toward the center of a subunit; R58 and D83 are equatorial domain residues. The allostery wiring diagram shows that the network of residues are interspersed throughout the structure. Residues D83, V174, E191, and D359 play a critical role in the R''-->T transition, which implies that mutations of these residues would compromise the ATPase activity. D83 and E191 are also highly conserved; D359 is moderately conserved. The negative cooperativity between the rings in the R''-->T transition is orchestrated through several interface residues within a single ring, including N10, E434, D435, and E451. Signal from the trans ring that is transmitted across the interface between the equatorial domains is responsible for the R''-->T transition. The cochaperonin GroES plays a passive role in the R''-->T transition. Remarkably, the binding affinity of GroES for GroEL is allosterically linked to GroEL residues 350-365 that span helices K and L. The movements of helices K and L alter the polarity of the cavity throughout the GroEL functional cycle and undergo large-scale motions that are anticorrelated with the other apical domain residues. The allostery wiring diagrams for the T-->R and R''-->T transitions of GroEL provide a microscopic foundation for the cooperativity (anticooperativity) within (between) the ring (rings). Using statistical coupling analysis, we extract evolutionarily linked clusters of residues in GroEL and GroES. We find that several substrate protein binding residues as well as sites related to ATPase activity belong to a single functional network in GroEL. For GroES, the mobile loop residues and GroES/GroES interface residues are linked.

Allosteric transitions in the chaperonin GroEL are captured by a dominant normal mode that is most robust to sequence variations

The Escherichia coli chaperonin GroEL, which helps proteins to fold, consists of two heptameric rings stacked back-to-back. During the reaction cycle GroEL undergoes a series of allosteric transitions triggered by ligand (substrate protein, ATP, and the cochaperonin GroES) binding. Based on an elastic network model of the bullet-shaped double-ring chaperonin GroEL-(ADP)(7)-GroES structure (R''T state), we perform a normal mode analysis to explore the energetically favorable collective motions encoded in the R''T structure. By comparing each normal mode with the observed conformational changes in the R''T --> TR'' transition, a single dominant normal mode provides a simple description of this highly intricate allosteric transition. A detailed analysis of this relatively high-frequency mode describes the structural and dynamic changes that underlie the positive intra-ring and negative inter-ring cooperativity. The dynamics embedded in the dominant mode entails highly concerted structural motions with approximate preservation of sevenfold symmetry within each ring and negatively correlated ones between the two rings. The dominant normal mode (in comparison with the other modes) is robust to parametric perturbations caused by sequence variations, which validates its functional importance. Response of the dominant mode to local changes that mimic mutations using the structural perturbation method technique leads to a wiring diagram that identifies a network of key residues that regulate the allosteric transitions. Many of these residues are located in intersubunit interfaces, and may therefore play a critical role in transmitting allosteric signals between subunits.

Adaptive evolution in GroEL from distantly related endosymbiotic bacteria of insects

Many symbioses between bacteria and insects resulted from ancient infections followed by strict vertical transmission within host lineages. The strong bottlenecks under which this transmission occurs promote the neutral fixation of slightly deleterious mutations by genetic drift. As predicted by Muller's ratchet, this fixation will drive endosymbiotic bacteria through an irreversible dynamics of fitness loss. The chaperonin GroEL has been proposed as a compensatory mechanism whereby endosymbiotic bacteria of aphids persist. Here, we show that endosymbiotic bacteria of insects from two phylogenetically very distant bacterial phyla have fixed amino acid substitutions by positive selection in functionally important GroEL regions involved in either GroES/peptide binding or in the en bloc movement of the GroEL apical domain. These results, together with the high levels of constitutive expression of GroEL in these endosymbionts, provide valuable insights into the evolution of a molecular mechanism responsible for the maintenance of the symbiotic lifestyle.

The N Terminus of the Head Protein of T4 Bacteriophage Directs Proteins to the GroEL Chaperonin

The head protein of T4 bacteriophage requires the GroEL chaperonin for its insertion into a growing T4 head. Hundreds of thousands of copies of this protein must pass through the chaperonin in a limited time later in infection, indicating that the protein must use GroEL very efficiently and may contain sequences that bind tightly to GroEL. We show that green fluorescent protein (GFP) fused to the N terminus of the head protein can fold at temperatures higher than those at which the GFP protein can fold well by itself. We present evidence that this folding is promoted by the strong binding of N-terminal head protein sequences to GroEL. This binding is so strong that some fusion proteins can apparently deplete the cell of the GroEL needed for other cellular functions, altering the cellular membranes and slowing growth.

A mobile loop order-disorder transition modulates the speed of chaperonin cycling

Molecular machines order and disorder polypeptides as they form and dissolve large intermolecular interfaces, but the biological significance of coupled ordering and binding has been established in few, if any, macromolecular systems. The ordering and binding of GroES co-chaperonin mobile loops accompany an ATP-dependent conformational change in the GroEL chaperonin that promotes client protein folding. Following ATP hydrolysis, disordering of the mobile loops accompanies co-chaperonin dissociation, reversal of the GroEL conformational change, and release of the client protein. "High-affinity" GroEL mutants were identified by their compatibility with "low-affinity" co-chaperonin mutants and incompatibility with high-affinity co-chaperonin mutants. Analysis of binding kinetics using the intrinsic fluorescence of tryptophan-containing co-chaperonin variants revealed that excessive affinity causes the chaperonin to stall in a conformation that forms in the presence of ATP. Destabilizing the beta-hairpins formed by the mobile loops restores the normal rate of dissociation. Thus, the free energy of mobile-loop ordering and disordering acts like the inertia of an engine's flywheel by modulating the speed of chaperonin conformational changes.

Regulation of Translation of the Head Protein of T4 Bacteriophage by Specific Binding of EF-Tu to a Leader Sequence

Recent evidence indicates that translation elongation factor Tu (EF-Tu) has a role in the cell in addition to its well established role in translation. The translation factor binds to a specific region called the Gol region close to the N terminus of the T4 bacteriophage major head protein as the head protein emerges from the ribosome. This binding was discovered because EF-Tu bound to Gol peptide is the specific substrate of the Lit protease that cleaves the EF-Tu between amino acid residues Gly59 and lle60, blocking phage development. These experiments raised the question of why the Gol region of the incipient head protein binds to EF-Tu, as binding to incipient proteins is not expected from the canonical role of EF-Tu. Here, we use gol-lacZ translational fusions to show that cleavage of EF-Tu in the complex with Gol peptide can block translation of a lacZ reporter gene fused translationally downstream of the Gol peptide that activated the cleavage. We propose a model to explain how binding of EF-Tu to the emerging Gol peptide could cause translation to pause temporarily and allow time for the leader polypeptide to bind to the GroEL chaperonin before translation continues, allowing cotranslation of the head protein with its insertion into the GroEL chaperonin chamber, and preventing premature synthesis and precipitation of the head protein. Cleavage of EF-Tu in the complex would block translation of the head protein and therefore development of the infecting phage. Experiments are presented that confirm two predictions of this model. Considering the evolutionary conservation of the components of this system, this novel regulatory mechanism could be used in other situations, both in bacteria and eukaryotes, where proteins are cotranslated with their insertion into cellular structures.

Functional characterization of an archaeal GroEL/GroES chaperonin system: significance of substrate encapsulation

In all three kingdoms of life chaperonins assist the folding of a range of newly synthesized proteins. As shown recently, Archaea of the genus Methanosarcina contain both group I (GroEL/GroES) and group II (thermosome) chaperonins in the cytosol. Here we report on a detailed functional analysis of the archaeal GroEL/GroES system of Methanosarcina mazei (Mm) in comparison to its bacterial counterpart from Escherichia coli (Ec). We find that the groESgroEL operon of M. mazei is unable to functionally replace groESgroEL in E. coli. However, the MmGroES protein can largely complement a mutant EcGroES protein in vivo. The ATPase rate of MmGroEL is very low and the dissociation of MmGroES from MmGroEL is 15 times slower than for the EcGroEL/GroES system. This slow ATPase cycle results in a prolonged enclosure time for model substrate proteins, such as rhodanese, in the MmGroEL:GroES folding cage before their release into the medium. Interestingly, optimal functionality of MmGroEL/GroES and its ability to encapsulate larger proteins, such as malate dehydrogenase, requires the presence of ammonium sulfate in vitro. In the absence of ammonium sulfate, malate dehydrogenase fails to be encapsulated by GroES and rather cycles on and off the GroEL trans ring in a non-productive reaction. These results indicate that the archaeal GroEL/GroES system has preserved the basic encapsulation mechanism of bacterial GroEL and suggest that it has adjusted the length of its reaction cycle to the slower growth rates of Archaea. Additionally, the release of only the folded protein from the GroEL/GroES cage may prevent adverse interactions of the GroEL substrates with the thermosome, which is not normally located within the same compartment.

Dobzhansky-Muller incompatibilities in protein evolution

We study fitness landscape in the space of protein sequences by relating sets of human pathogenic missense mutations in 32 proteins to amino acid substitutions that occurred in the course of evolution of these proteins. On average, approximately 10% of deviations of a nonhuman protein from its human ortholog are compensated pathogenic deviations (CPDs), i.e., are caused by an amino acid substitution that, at this site, would be pathogenic to humans. Normal functioning of a CPD-containing protein must be caused by other, compensatory deviations of the nonhuman species from humans. Together, a CPD and the corresponding compensatory deviation form a Dobzhansky-Muller incompatibility that can be visualized as the corner on a fitness ridge. Thus, proteins evolve along fitness ridges which contain only approximately 10 steps between successive corners. The fraction of CPDs among all deviations of a protein from its human ortholog does not increase with the evolutionary distance between the proteins, indicating that substitutions that carry evolving proteins around these corners occur in rapid succession, driven by positive selection. Data on fitness of interspecies hybrids suggest that the compensatory change that makes a CPD fit usually occurs within the same protein. Data on protein structures and on cooccurrence of amino acids at different sites of multiple orthologous proteins often make it possible to provisionally identify the substitution that compensates a particular CPD.

Mapping pathways of allosteric communication in GroEL by analysis of correlated mutations

An interesting example of an allosteric protein is the chaperonin GroEL. It undergoes adenosine 5'-triphosphate-induced conformational changes that are reflected in binding of adenosine 5'-triphosphate with positive cooperativity within rings and negative cooperativity between rings. Herein, correlated mutations in chaperonins are analyzed to unravel routes of allosteric communication in GroEL and in its complex with its co-chaperonin GroES. It is shown that analysis of correlated mutations in the chaperonin family can provide information about pathways of allosteric communication within GroEL and between GroEL and GroES. The results are discussed in the context of available structural, genetic, and biochemical data concerning short- and long-range interactions in the GroE system.