The effect of chaperonin buffering on protein evolution

Protein Quality Control is a Master Modulator of Molecular Evolution in Bacteria

The bacterial protein quality control (PQC) network comprises a set of genes that promote proteostasis (proteome homeostasis) through proper protein folding and function via chaperones, proteases, and protein translational machinery. It participates in vital cellular processes and influences organismal development and evolution. In this review, we examine the mechanistic bases for how the bacterial PQC network influences molecular evolution. We discuss the relevance of PQC components to contemporary issues in evolutionary biology including epistasis, evolvability, and the navigability of protein space. We examine other areas where proteostasis affects aspects of evolution and physiology, including host-parasite interactions. More generally, we demonstrate that the study of bacterial systems can aid in broader efforts to understand the relationship between genotype and phenotype across the biosphere.

Expression divergence of BAG gene family in maize under heat stress

Heat stress poses a significant challenge for maize production, especially during the spring when high temperatures disrupt cellular processes, impeding plant growth and development. The B-cell lymphoma-2 (Bcl-2) associated athanogene (BAG) gene family is known to be relatively conserved across various species. It plays a crucial role as molecular chaperone cofactors that are responsible for programmed cell death and tumorigenesis. Once the plant is under heat stress, the BAG genes act as co-chaperones and modulate the molecular functions of HSP70/HSC70 saving the plant from the damage of high temperature stress. The study was planned to identify and characterize the BAG genes for heat stress responsiveness in maize. Twenty-one (21) BAG genes were identified in the latest maize genome. The evolutionary relationship of Zea mays BAGs (ZmBAGs) with Arabidopsis thaliana, Solanum lycopersicum, Theobroma cacao, Sorghum bicolor, Ananas comosus, Physcomitrium patens, Oryza sativa and Populus trichocarpa were represented by the phylogenetic analysis. Differential expressions of BAG gene family in leaf, endosperm, anther, silk, seed and developing embryo depict their contribution to the growth and development. The in-silico gene expression analysis indicated ZmBAG-8 (Zm00001eb170080), and ZmBAG-11 (Zm00001eb237960) showed higher expression under abiotic stresses (cold, heat and salinity). The RT-qPCR further confirmed the expression of ZmBAG-8 and ZmBAG-11 in plant leaf tissue across the contrasting inbred lines and their F1 hybrid (DR-139, UML-1 and DR-139 × UML-1) when exposed to heat stress. Furthermore, the protein-protein interaction networks of ZmBAG-8 and ZmBAG-11 further elucidated their role in stress tolerance related pathways. This research offers a roadmap to plan functional research and utilize ZmBAG genes to enhance heat tolerance in grasses.

The Hsp90 Molecular Chaperone as a Global Modifier of the Genotype-Phenotype-Fitness Map: An Evolutionary Perspective

Global modifier genes influence the mapping of genotypes onto phenotypes and fitness through their epistatic interactions with genetic variants on a massive scale. The first such factor to be identified, Hsp90, is a highly conserved molecular chaperone that plays a central role in protein homeostasis. Hsp90 is a "hub of hubs" that chaperones proteins engaged in many key cellular and developmental regulatory networks. These clients, which are enriched in kinases, transcription factors, and E3 ubiquitin ligases, drive diverse cellular functions and are themselves highly connected. By contrast to many other hub proteins, the abundance and activity of Hsp90 changes substantially in response to shifting environmental conditions. As a result, Hsp90 modifies the functional impact of many genetic variants simultaneously in a manner that depends on environmental stress. Studies in diverse organisms suggest that this coupling between Hsp90 function and challenging environments exerts a substantial impact on what parts of the genome are visible to natural selection, expanding adaptive opportunities when most needed. In this Perspective, we explore the multifaceted role of Hsp90 as global modifier of the genotype-phenotype-fitness map as well as its implications for evolution in nature and the clinic.

Dominant-negative chaperonin mutation ptCPN60α1S57F uncovers redundancy in chloroplast rRNA processing

The biogenesis of functional forms of chloroplast ribosomal RNAs (rRNAs) is crucial for the translation of chloroplast mRNAs into polypeptides. However, the molecular mechanisms underlying the proper processing and maturation of chloroplast rRNA species are poorly understood. Through a genetic approach, we isolated and characterized an Arabidopsis mutant, α1-4, harboring a missense mutation in the plastid chaperonin-60α1 gene. Using allelism tests and transgenic manipulation, we determined functional redundancy among ptCPN60 subunits. The ptCPN60α1S57F mutation caused specific defects in the formation of chloroplast rRNA species, including 23S, 5S, and 4.5S rRNAs, but not 16S rRNAs. Allelism tests suggested that the dysfunctional ptCPN60α1S57F competes with other members of the ptCPN60 family. Indeed, overexpression of the ptCPN60α1S57F protein in wild-type plants mimicked the phenotypes observed in the α1-4 mutant, while increasing the endogenous transcriptional levels of ptCPN60α2, β1, β2, and β3 in the α1-4 mutant partially mitigated the abnormal fragmentation processing of chloroplast 23S, 5S, and 4.5S rRNAs. Furthermore, we demonstrated functional redundancy between ptCPN60β1 and ptCPN60β2 in chloroplast rRNA processing through double-mutant analysis. Collectively, our data reveal a novel physiological role of ptCPN60 subunits in generating the functional rRNA species of the large 50S ribosomal subunit in Arabidopsis chloroplasts.

The ribosome lowers the entropic penalty of protein folding

Most proteins fold during biosynthesis on the ribosome1, and co-translational folding energetics, pathways and outcomes of many proteins have been found to differ considerably from those in refolding studies2-10. The origin of this folding modulation by the ribosome has remained unknown. Here we have determined atomistic structures of the unfolded state of a model protein on and off the ribosome, which reveal that the ribosome structurally expands the unfolded nascent chain and increases its solvation, resulting in its entropic destabilization relative to the peptide chain in isolation. Quantitative 19F NMR experiments confirm that this destabilization reduces the entropic penalty of folding by up to 30 kcal mol-1 and promotes formation of partially folded intermediates on the ribosome, an observation that extends to other protein domains and is obligate for some proteins to acquire their active conformation. The thermodynamic effects also contribute to the ribosome protecting the nascent chain from mutation-induced unfolding, which suggests a crucial role of the ribosome in supporting protein evolution. By correlating nascent chain structure and dynamics to their folding energetics and post-translational outcomes, our findings establish the physical basis of the distinct thermodynamics of co-translational protein folding.

Mutational robustness and the role of buffer genes in evolvability

Organisms rely on mutations to fuel adaptive evolution. However, many mutations impose a negative effect on fitness. Cells may have therefore evolved mechanisms that affect the phenotypic effects of mutations, thus conferring mutational robustness. Specifically, so-called buffer genes are hypothesized to interact directly or indirectly with genetic variation and reduce its effect on fitness. Environmental or genetic perturbations can change the interaction between buffer genes and genetic variation, thereby unmasking the genetic variation's phenotypic effects and thus providing a source of variation for natural selection to act on. This review provides an overview of our understanding of mutational robustness and buffer genes, with the chaperone gene HSP90 as a key example. It discusses whether buffer genes merely affect standing variation or also interact with de novo mutations, how mutational robustness could influence evolution, and whether mutational robustness might be an evolved trait or rather a mere side-effect of complex genetic interactions.

Viral Evolution Shaped by Host Proteostasis Networks

Understanding the factors that shape viral evolution is critical for developing effective antiviral strategies, accurately predicting viral evolution, and preventing pandemics. One fundamental determinant of viral evolution is the interplay between viral protein biophysics and the host machineries that regulate protein folding and quality control. Most adaptive mutations in viruses are biophysically deleterious, resulting in a viral protein product with folding defects. In cells, protein folding is assisted by a dynamic system of chaperones and quality control processes known as the proteostasis network. Host proteostasis networks can determine the fates of viral proteins with biophysical defects, either by assisting with folding or by targeting them for degradation. In this review, we discuss and analyze new discoveries revealing that host proteostasis factors can profoundly shape the sequence space accessible to evolving viral proteins. We also discuss the many opportunities for research progress proffered by the proteostasis perspective on viral evolution and adaptation.

Hsp90 is a potential risk factor for ovarian cancer prognosis: an evidence of a Chinese clinical center

BACKGROUND

The potential treatment effects of heat shock protein 90 (Hsp90) inhibitors in ovarian cancer (OC) are controversial. This research aims to investigate the relationship between the level of Hsp90 in peripheral blood and the prognosis of OC patients, as well as the clinicopathological indicators.

MATERIALS AND METHODS

We retrospectively collected the clinicopathological indicators of OC patients who were admitted to the Department of Obstetrics and Gynecology of the First Affiliated Hospital of Anhui Medical University from 2017 to 2022. Hsp90 level in patient blood was detected by enzyme-linked immunosorbent assay, and the correlation between Hsp90 level and OC prognosis was systematically investigated. Kaplan-Meier method was used to draw the survival curve, and the average survival time and survival rate were calculated. The log-rank test and Cox model were used for univariate survival analysis, and the Cox proportional hazards model was applied for multivariate survival analysis. Based on the TCGA dataset of OC obtained by cBioPortal, Pearson's correlation coefficients between Hsp90 level values and other mRNA expression values were calculated to further conduct bioinformatics analysis. GSEA and GSVA analysis were also conducted for gene functional enrichment. The expression of Hsp90 in OC tissues were evaluated and compared by Immunohistochemical staining.

RESULTS

According to the established screening criteria, 106 patients were selected. The enzyme-linked immunosorbent assay results showed that 50.94% OC patients with abnormal Hsp90 level. According to the outcome of Kaplan-Meier curves, the results revealed that the abnormal level of Hsp90 was suggested to poor prognosis (P = 0.001) of OC patients. Furthermore, the result of multivariate Cox proportional hazards regression model analysis also predicted that abnormal Hsp90 level (HR = 2.838, 95%CI = 1.139-7.069, P = 0.025) was linked to poor prognosis, which could be an independent prognostic factor for the prognosis of OC patients. Moreover, top 100 genes screened by Pearson's value associated with Hsp90, indicating that Hsp90 participated in the regulation of ATF5 target genes, PRAGC1A target genes and BANP target genes and also enriched in the metabolic processes of cell response to DNA damage stimulus, response to heat and protein folding.

CONCLUSION

Hsp90 level is positively associated with OC mortality and is a potential prognostic indicator of OC.

DnaK response to expression of protein mutants is dependent on translation rate and stability

Chaperones play a central part in the quality control system in cells by clearing misfolded and aggregated proteins. The chaperone DnaK acts as a sensor for molecular stress by recognising short hydrophobic stretches of misfolded proteins. As the level of unfolded protein is a function of protein stability, we hypothesised that the level of DnaK response upon overexpression of recombinant proteins would be correlated to stability. Using a set of mutants of the λ-repressor with varying thermal stabilities and a fluorescent reporter system, the effect of stability on DnaK response and protein abundance was investigated. Our results demonstrate that the initial DnaK response is largely dependent on protein synthesis rate but as the recombinantly expressed protein accumulates and homeostasis is approached the response correlates strongly with stability. Furthermore, we observe a large degree of cell-cell variation in protein abundance and DnaK response in more stable proteins.

Non-adaptive complexity and biochemical function

Intricate biochemical structures are usually thought to be useful, because natural selection preserves them from degradation by a constant hail of destructive mutations. Biochemists therefore often deliberately disrupt them to understand how complexity improves protein function or fitness. However, evolutionary theory suggests that even useless complexity that never improved fitness can become completely essential if a simple set of evolutionary conditions is fulfilled. We review evidence that stable protein complexes, protein-chaperone interactions, and complexes consisting of several paralogs all fulfill these conditions. This makes reverse genetics or destructive mutagenesis unsuitable for assigning functions to these kinds of complexity. Instead, we advocate that incorporating evolutionary approaches into biochemistry overcomes this difficulty and allows us to distinguish useless from useful biochemical complexity.

The endoplasmic reticulum proteostasis network profoundly shapes the protein sequence space accessible to HIV envelope

The sequence space accessible to evolving proteins can be enhanced by cellular chaperones that assist biophysically defective clients in navigating complex folding landscapes. It is also possible, at least in theory, for proteostasis mechanisms that promote strict quality control to greatly constrain accessible protein sequence space. Unfortunately, most efforts to understand how proteostasis mechanisms influence evolution rely on artificial inhibition or genetic knockdown of specific chaperones. The few experiments that perturb quality control pathways also generally modulate the levels of only individual quality control factors. Here, we use chemical genetic strategies to tune proteostasis networks via natural stress response pathways that regulate the levels of entire suites of chaperones and quality control mechanisms. Specifically, we upregulate the unfolded protein response (UPR) to test the hypothesis that the host endoplasmic reticulum (ER) proteostasis network shapes the sequence space accessible to human immunodeficiency virus-1 (HIV-1) envelope (Env) protein. Elucidating factors that enhance or constrain Env sequence space is critical because Env evolves extremely rapidly, yielding HIV strains with antibody- and drug-escape mutations. We find that UPR-mediated upregulation of ER proteostasis factors, particularly those controlled by the IRE1-XBP1s UPR arm, globally reduces Env mutational tolerance. Conserved, functionally important Env regions exhibit the largest decreases in mutational tolerance upon XBP1s induction. Our data indicate that this phenomenon likely reflects strict quality control endowed by XBP1s-mediated remodeling of the ER proteostasis environment. Intriguingly, and in contrast, specific regions of Env, including regions targeted by broadly neutralizing antibodies, display enhanced mutational tolerance when XBP1s is induced, hinting at a role for host proteostasis network hijacking in potentiating antibody escape. These observations reveal a key function for proteostasis networks in decreasing instead of expanding the sequence space accessible to client proteins, while also demonstrating that the host ER proteostasis network profoundly shapes the mutational tolerance of Env in ways that could have important consequences for HIV adaptation.

The host cell’s endoplasmic reticulum proteostasis network has a profound, constraining impact on the protein sequence space accessible to HIV’s envelope protein, which is a major target of the host’s adaptive immune system; in particular, upregulation of stringent quality control pathways appears to restrict the viability of destabilizing envelope variants.

Chaperonin Abundance Enhances Bacterial Fitness

The ability of chaperonins to buffer mutations that affect protein folding pathways suggests that their abundance should be evolutionarily advantageous. Here, we investigate the effect of chaperonin overproduction on cellular fitness in Escherichia coli. We demonstrate that chaperonin abundance confers 1) an ability to tolerate higher temperatures, 2) improved cellular fitness, and 3) enhanced folding of metabolic enzymes, which is expected to lead to enhanced energy harvesting potential.

Chaperone client proteins evolve slower than non-client proteins

Chaperones are important molecular machinery that assists proteins to attain their native three-dimensional structure crucial for function. Earlier studies using experimental evolution showed that chaperones impose a relaxation of sequence constraints on their "client" proteins, which may lead to the fixation of slightly deleterious mutations on the latter. However, we hypothesized that such a phenomenon might be harmful to the organism in a natural physiological condition. In this study, we investigated the evolutionary rates of chaperone client and non-client proteins in five model organisms from both prokaryotic and eukaryotic lineages. Our study reveals a slower evolutionary rate of chaperone client proteins in all five organisms. Additionally, the slower folding rate and lower aggregation propensity of chaperone client proteins reveal that the chaperone may play an essential role in rescuing the slightly disadvantageous effects due to random mutations and subsequent protein misfolding. However, the fixation of such mutations is less likely to be selected in the natural population.

Molecular Chaperones Accelerate the Evolution of Their Protein Clients in Yeast

Protein stability is a major constraint on protein evolution. Molecular chaperones, also known as heat-shock proteins, can relax this constraint and promote protein evolution by diminishing the deleterious effect of mutations on protein stability and folding. This effect, however, has only been stablished for a few chaperones. Here, we use a comprehensive chaperone–protein interaction network to study the effect of all yeast chaperones on the evolution of their protein substrates, that is, their clients. In particular, we analyze how yeast chaperones affect the evolutionary rates of their clients at two very different evolutionary time scales. We first study the effect of chaperone-mediated folding on protein evolution over the evolutionary divergence of Saccharomyces cerevisiae and S. paradoxus. We then test whether yeast chaperones have left a similar signature on the patterns of standing genetic variation found in modern wild and domesticated strains of S. cerevisiae. We find that genes encoding chaperone clients have diverged faster than genes encoding non-client proteins when controlling for their number of protein–protein interactions. We also find that genes encoding client proteins have accumulated more intraspecific genetic diversity than those encoding non-client proteins. In a number of multivariate analyses, controlling by other well-known factors that affect protein evolution, we find that chaperone dependence explains the largest fraction of the observed variance in the rate of evolution at both evolutionary time scales. Chaperones affecting rates of protein evolution mostly belong to two major chaperone families: Hsp70s and Hsp90s. Our analyses show that protein chaperones, by virtue of their ability to buffer destabilizing mutations and their role in modulating protein genotype–phenotype maps, have a considerable accelerating effect on protein evolution.

Chaperonin overproduction and metabolic erosion caused by mutation accumulation in Escherichia coli

Bacterial cells adapting to a constant environment tend to accumulate mutations in portions of their genome that are not maintained by selection. This process has been observed in bacteria evolving under strong genetic drift, and especially in bacterial endosymbionts of insects. Here, we study this process in hypermutable Escherichia coli populations evolved through 250 single-cell bottlenecks on solid rich medium in a mutation accumulation experiment that emulates the evolution of bacterial endosymbionts. Using phenotype microarrays monitoring metabolic activity in 95 environments distinguished by their carbon sources, we observe how mutation accumulation has decreased the ability of cells to metabolize most carbon sources. We study if the chaperonin GroEL, which is naturally overproduced in bacterial endosymbionts, can ameliorate the process of metabolic erosion, because of its known ability to buffer destabilizing mutations in metabolic enzymes. Our results indicate that GroEL can slow down the negative phenotypic consequences of genome decay in some environments.

Metabolic Determinants of Enzyme Evolution in a Genome-Scale Bacterial Metabolic Network

Different genes and proteins evolve at very different rates. To identify the factors that explain these differences is an important aspect of research in molecular evolution. One such factor is the role a protein plays in a large molecular network. Here, we analyze the evolutionary rates of enzyme-coding genes in the genome-scale metabolic network of Escherichia coli to find the evolutionary constraints imposed by the structure and function of this complex metabolic system. Central and highly connected enzymes appear to evolve more slowly than less connected enzymes, but we find that they do so as a by-product of their high abundance, and not because of their position in the metabolic network. In contrast, enzymes catalyzing reactions with high metabolic flux-high substrate to product conversion rates-evolve slowly even after we account for their abundance. Moreover, enzymes catalyzing reactions that are difficult to by-pass through alternative pathways, such that they are essential in many different genetic backgrounds, also evolve more slowly. Our analyses show that an enzyme's role in the function of a metabolic network affects its evolution more than its place in the network's structure. They highlight the value of a system-level perspective for studies of molecular evolution.

Inhibition of HSP90α protects cultured neurons from oxygen-glucose deprivation induced necroptosis by decreasing RIP3 expression

Heat shock protein 90α (HSP90α) maintains cell stabilization and regulates cell death, respectively. Recent studies have shown that HSP90α is involved in receptor interacting protein 3 (RIP3)-mediated necroptosis in HT29 cells. It is known that oxygen and glucose deprivation (OGD) can induce necroptosis, which is regulated by RIP3 in neurons. However, it is still unclear whether HSP90α participates in the process of OGD-induced necroptosis in cultured neurons via the regulation of RIP3. Our study found that necroptosis occurs in primary cultured cortical neurons and PC-12 cells following exposure to OGD insult. Additionally, the expression of RIP3/p-RIP3, MLKL/p-MLKL, and the RIP1/RIP3 complex (necrosome) significantly increased following OGD, as measured through immunofluorescence (IF) staining, Western blotting (WB), and immunoprecipitation (IP) assay. Additionally, data from computer simulations and IP assays showed that HSP90α interacts with RIP3. In addition, HSP90α was overexpressed following OGD in cultured neurons, as measured through WB and IF staining. Inhibition of HSP90α in cultured neurons, using the specific inhibitor, geldanamycin (GA), and siRNA/shRNA of HSP90α, protected cultured neurons from necrosis. Our study showed that the inhibitor of HSP90α, GA, rescued cultured neurons not only by decreasing the expression of total RIP3/MLKL, but also by decreasing the expression of p-RIP3/p-MLKL and the RIP1/RIP3 necrosome. In this study, we reveal that inhibition of HSP90α protects primary cultured cortical neurons and PC-12 cells from OGD-induced necroptosis through the modulation of RIP3 expression.

Hsp90 shapes protein and RNA evolution to balance trade-offs between protein stability and aggregation

Acquisition of mutations is central to evolution; however, the detrimental effects of most mutations on protein folding and stability limit protein evolvability. Molecular chaperones, which suppress aggregation and facilitate polypeptide folding, may alleviate the effects of destabilizing mutations thus promoting sequence diversification. To illuminate how chaperones can influence protein evolution, we examined the effect of reduced activity of the chaperone Hsp90 on poliovirus evolution. We find that Hsp90 offsets evolutionary trade-offs between protein stability and aggregation. Lower chaperone levels favor variants of reduced hydrophobicity and protein aggregation propensity but at a cost to protein stability. Notably, reducing Hsp90 activity also promotes clusters of codon-deoptimized synonymous mutations at inter-domain boundaries, likely to facilitate cotranslational domain folding. Our results reveal how a chaperone can shape the sequence landscape at both the protein and RNA levels to harmonize competing constraints posed by protein stability, aggregation propensity, and translation rate on successful protein biogenesis.

Comparative genomic analysis of mollicutes with and without a chaperonin system

The GroE chaperonin system, which comprises GroEL and GroES, assists protein folding in vivo and in vitro. It is conserved in all prokaryotes except in most, but not all, members of the class of mollicutes. In Escherichia coli, about 60 proteins were found to be obligatory clients of the GroE system. Here, we describe the properties of the homologs of these GroE clients in mollicutes and the evolution of chaperonins in this class of bacteria. Comparing the properties of these homologs in mollicutes with and without chaperonins enabled us to search for features correlated with the presence of GroE. Interestingly, no sequence-based features of proteins such as average length, amino acid composition and predicted folding/disorder propensity were found to be affected by the absence of GroE. Other properties such as genome size and number of proteins were also found to not differ between mollicute species with and without GroE. Our data suggest that two clades of mollicutes re-acquired the GroE system, thereby supporting the view that gaining the system occurred polyphyletically and not monophyletically, as previously debated. Our data also suggest that there might have been three isolated cases of lateral gene transfer from specific bacterial sources. Taken together, our data indicate that loss of GroE does not involve crossing a high evolutionary barrier and can be compensated for by a small number of changes within the few dozen client proteins.

Disulfide Bonds Enable Accelerated Protein Evolution

The different proteins of any proteome evolve at enormously different rates. What factors contribute to this variability, and to what extent, is still a largely open question. We hypothesized that disulfide bonds, by increasing protein stability, should make proteins' structures relatively independent of their amino acid sequences, thus acting as buffers of deleterious mutations and enabling accelerated sequence evolution. In agreement with this hypothesis, we observed that membrane proteins with disulfide bonds evolved 88% faster than those without disulfide bonds, and that extracellular proteins with disulfide bonds evolved 49% faster than those without disulfide bonds. In addition, genes encoding proteins with disulfide bonds exhibit an increased likelihood of showing signatures of positive selection. Multivariate analyses indicate that the trend is independent of a number of potentially confounding factors. The effect, however, is not observed among the longest proteins, which can become stabilized by mechanisms other than disulfide bonds.

The versatile mutational "repertoire" of Escherichia coli GroEL, a multidomain chaperonin nanomachine

The bacterial chaperonins are highly sophisticated molecular nanomachines, controlled by the hydrolysis of ATP to dynamically trap and remove from the environment unstable protein molecules that are susceptible to denaturation and aggregation. Chaperonins also act to assist in the refolding of these unstable proteins, providing a means by which these proteins may return in active form to the complex environment of the cell. The Escherichia coli GroE chaperonin system is one of the largest protein supramolecular complexes known, whose quaternary structure is required for segregating aggregation-prone proteins. Over the course of more than two decades of research on GroE, it has become accepted that GroE, more specifically the GroEL subunit, is a "high-tolerance" molecular system, capable of accommodating numerous mutations, while retaining its molecular integrity. In some cases, a given site of mutation was revealed to be absolutely required for GroEL function, providing hints regarding the network of signals and triggers that propel this unique system. In other instances, however, a mutation has produced a more delicate response, altering only part of, or in some cases, only a single facet of, the molecular mechanism, and these mutants have often provided invaluable hints on the extent of the complexity underlying chaperonin-assisted protein folding. In this review, we highlight some examples of the latter type of GroEL mutants which compose the unique "mutational repertoire" of GroEL and touch upon the important clues that each mutant provided to the overall effort to elucidate the details of GroE action.

Interactome evolution: insights from genome-wide analyses of protein-protein interactions

We highlight new evolutionary insights enabled by recent genome-wide studies on protein-protein interaction (PPI) networks ('interactomes'). While most PPIs are mediated by a single sequence region promoting or inhibiting interactions, many PPIs are mediated by multiple sequence regions acting cooperatively. Most PPIs perform important functions maintained by negative selection: we estimate that less than ∼10% of the human interactome is effectively neutral upon perturbation (i.e. 'junk' PPIs), and the rest are deleterious upon perturbation; interfacial sites evolve more slowly than other sites; many conserved PPIs show signatures of co-evolution at the interface; PPIs evolve more slowly than protein sequence. At the same time, many PPIs undergo rewiring during evolution for lineage-specific adaptation. Finally, chaperone-protein and host-pathogen interactomes are governed by distinct evolutionary principles.

Evolution of Chaperonin Gene Duplication in Stigonematalean Cyanobacteria (Subsection V)

Chaperonins promote protein folding and are known to play a role in the maintenance of cellular stability under stress conditions. The group I bacterial chaperonin complex comprises GroEL, that forms a barrel-like oligomer, and GroES that forms the lid. In most eubacteria the GroES/GroEL chaperonin is encoded by a single-copy bicistronic operon, whereas in cyanobacteria up to three groES/groEL paralogs have been documented. Here we study the evolution and functional diversification of chaperonin paralogs in the heterocystous, multi-seriate filament forming cyanobacterium Chlorogloeopsis fritschii PCC 6912. The genome of C. fritschii encodes two groES/groEL operons (groESL1, groESL1.2) and a monocistronic groEL gene (groEL2). A phylogenetic reconstruction reveals that the groEL2 duplication is as ancient as cyanobacteria, whereas the groESL1.2 duplication occurred at the ancestor of heterocystous cyanobacteria. A comparison of the groEL paralogs transcription levels under different growth conditions shows that they have adapted distinct transcriptional regulation. Our results reveal that groEL1 and groEL1.2 are upregulated during diazotrophic conditions and the localization of their promoter activity points towards a role in heterocyst differentiation. Furthermore, protein–protein interaction assays suggest that paralogs encoded in the two operons assemble into hybrid complexes. The monocistronic encoded GroEL2 is not forming oligomers nor does it interact with the co-chaperonins. Interaction between GroES1.2 and GroEL1.2 could not be documented, suggesting that the groESL1.2 operon does not encode a functional chaperonin complex. Functional complementation experiments in Escherichia coli show that only GroES1/GroEL1 and GroES1/GroEL1.2 can substitute the native operon. In summary, the evolutionary consequences of chaperonin duplication in cyanobacteria include the retention of groESL1 as a housekeeping gene, subfunctionalization of groESL1.2 and neofunctionalization of the monocistronic groEL2 paralog.

Mechanisms of Hsp90 regulation

Heat shock protein 90 (Hsp90) is a molecular chaperone that is involved in the activation of disparate client proteins. This implicates Hsp90 in diverse biological processes that require a variety of co-ordinated regulatory mechanisms to control its activity. Perhaps the most important regulator is heat shock factor 1 (HSF1), which is primarily responsible for upregulating Hsp90 by binding heat shock elements (HSEs) within Hsp90 promoters. HSF1 is itself subject to a variety of regulatory processes and can directly respond to stress. HSF1 also interacts with a variety of transcriptional factors that help integrate biological signals, which in turn regulate Hsp90 appropriately. Because of the diverse clientele of Hsp90 a whole variety of co-chaperones also regulate its activity and some are directly responsible for delivery of client protein. Consequently, co-chaperones themselves, like Hsp90, are also subject to regulatory mechanisms such as post translational modification. This review, looks at the many different levels by which Hsp90 activity is ultimately regulated.

Biogenesis and Metabolic Maintenance of Rubisco

Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) mediates the fixation of atmospheric CO₂ in photosynthesis by catalyzing the carboxylation of the 5-carbon sugar ribulose-1,5-bisphosphate (RuBP). Rubisco is a remarkably inefficient enzyme, fixing only 2-10 CO₂ molecules per second. Efforts to increase crop yields by bioengineering Rubisco remain unsuccessful, owing in part to the complex cellular machinery required for Rubisco biogenesis and metabolic maintenance. The large subunit of Rubisco requires the chaperonin system for folding, and recent studies have shown that assembly of hexadecameric Rubisco is mediated by specific assembly chaperones. Moreover, Rubisco function can be inhibited by a range of sugar-phosphate ligands, including RuBP. Metabolic repair depends on remodeling of Rubisco by the ATP-dependent Rubisco activase and hydrolysis of inhibitory sugar phosphates by specific phosphatases. Here, we review our present understanding of the structure and function of these auxiliary factors and their utilization in efforts to engineer more catalytically efficient Rubisco enzymes.

Myxococcus xanthus DK1622 Coordinates Expressions of the Duplicate groEL and Single groES Genes for Synergistic Functions of GroELs and GroES

Chaperonin GroEL (Cpn60) requires cofactor GroES (Cpn10) for protein refolding in bacteria that possess single groEL and groES genes in a bicistronic groESL operon. Among 4,861 completely-sequenced prokaryotic genomes, 884 possess duplicate groEL genes and 770 possess groEL genes with no neighboring groES. It is unclear whether stand-alone groEL requires groES in order to function and, if required, how duplicate groEL genes and unequal groES genes balance their expressions. In Myxococcus xanthus DK1622, we determined that, while duplicate groELs were alternatively deletable, the single groES that clusters with groEL1 was essential for cell survival. Either GroEL1 or GroEL2 required interactions with GroES for in vitro and in vivo functions. Deletion of groEL1 or groEL2 resulted in decreased expressions of both groEL and groES; and ectopic complementation of groEL recovered not only the groEL but also groES expressions. The addition of an extra groES gene upstream groEL2 to form a bicistronic operon had almost no influence on groES expression and the cell survival rate, whereas over-expression of groES using a self-replicating plasmid simultaneously increased the groEL expressions. The results indicated that M. xanthus DK1622 cells coordinate expressions of the duplicate groEL and single groES genes for synergistic functions of GroELs and GroES. We proposed a potential regulation mechanism for the expression coordination.

Chance and necessity in the genome evolution of endosymbiotic bacteria of insects

An open question in evolutionary biology is how does the selection–drift balance determine the fates of biological interactions. We searched for signatures of selection and drift in genomes of five endosymbiotic bacterial groups known to evolve under strong genetic drift. Although most genes in endosymbiotic bacteria showed evidence of relaxed purifying selection, many genes in these bacteria exhibited stronger selective constraints than their orthologs in free-living bacterial relatives. Remarkably, most of these highly constrained genes had no role in the host–symbiont interactions but were involved in either buffering the deleterious consequences of drift or other host-unrelated functions, suggesting that they have either acquired new roles or their role became more central in endosymbiotic bacteria. Experimental evolution of Escherichia coli under strong genetic drift revealed remarkable similarities in the mutational spectrum, genome reduction patterns and gene losses to endosymbiotic bacteria of insects. Interestingly, the transcriptome of the experimentally evolved lines showed a generalized deregulation of the genome that affected genes encoding proteins involved in mutational buffering, regulation and amino acid biosynthesis, patterns identical to those found in endosymbiotic bacteria. Our results indicate that drift has shaped endosymbiotic associations through a change in the functional landscape of bacterial genes and that the host had only a small role in such a shift.

The Molecular Chaperone DnaK Is a Source of Mutational Robustness

Molecular chaperones, also known as heat-shock proteins, refold misfolded proteins and help other proteins reach their native conformation. Thanks to these abilities, some chaperones, such as the Hsp90 protein or the chaperonin GroEL, can buffer the deleterious phenotypic effects of mutations that alter protein structure and function. Hsp70 chaperones use a chaperoning mechanism different from that of Hsp90 and GroEL, and it is not known whether they can also buffer mutations. Here, we show that they can. To this end, we performed a mutation accumulation experiment in Escherichia coli, followed by whole-genome resequencing. Overexpression of the Hsp70 chaperone DnaK helps cells cope with mutational load and completely avoid the extinctions we observe in lineages evolving without chaperone overproduction. Additionally, our sequence data show that DnaK overexpression increases mutational robustness, the tolerance of its clients to nonsynonymous nucleotide substitutions. We also show that this elevated mutational buffering translates into differences in evolutionary rates on intermediate and long evolutionary time scales. Specifically, we studied the evolutionary rates of DnaK clients using the genomes of E. coli, Salmonella enterica, and 83 other gamma-proteobacteria. We find that clients that interact strongly with DnaK evolve faster than weakly interacting clients. Our results imply that all three major chaperone classes can buffer mutations and affect protein evolution. They illustrate how an individual protein like a chaperone can have a disproportionate effect on the evolution of a proteome.

Molecular Determinants of Mutant Phenotypes, Inferred from Saturation Mutagenesis Data

Understanding how mutations affect protein activity and organismal fitness is a major challenge. We used saturation mutagenesis combined with deep sequencing to determine mutational sensitivity scores for 1,664 single-site mutants of the 101 residue Escherichia coli cytotoxin, CcdB at seven different expression levels. Active-site residues could be distinguished from buried ones, based on their differential tolerance to aliphatic and charged amino acid substitutions. At nonactive-site positions, the average mutational tolerance correlated better with depth from the protein surface than with accessibility. Remarkably, similar results were observed for two other small proteins, PDZ domain (PSD95^pdz3) and IgG-binding domain of protein G (GB1). Mutational sensitivity data obtained with CcdB were used to derive a procedure for predicting functional effects of mutations. Results compared favorably with those of two widely used computational predictors. In vitro characterization of 80 single, nonactive-site mutants of CcdB showed that activity in vivo correlates moderately with thermal stability and solubility. The inability to refold reversibly, as well as a decreased folding rate in vitro, is associated with decreased activity in vivo. Upon probing the effect of modulating expression of various proteases and chaperones on mutant phenotypes, most deleterious mutants showed an increased in vivo activity and solubility only upon over-expression of either Trigger factor or SecB ATP-independent chaperones. Collectively, these data suggest that folding kinetics rather than protein stability is the primary determinant of activity in vivo. This study enhances our understanding of how mutations affect phenotype, as well as the ability to predict fitness effects of point mutations.

DnaK-Dependent Accelerated Evolutionary Rate in Prokaryotes

Many proteins depend on an interaction with molecular chaperones in order to fold into a functional tertiary structure. Previous studies showed that protein interaction with the GroEL/GroES chaperonine and Hsp90 chaperone can buffer the impact of slightly deleterious mutations in the protein sequence. This capacity of GroEL/GroES to prevent protein misfolding has been shown to accelerate the evolution of its client proteins. Whether other bacterial chaperones have a similar effect on their client proteins is currently unknown. Here, we study the impact of DnaK (Hsp70) chaperone on the evolution of its client proteins. Evolutionary parameters were derived from comparison of the Escherichia coli proteome to 1,808,565 orthologous proteins in 1,149 proteobacterial genomes. Our analysis reveals a significant positive correlation between protein binding frequency with DnaK and evolutionary rate. Proteins with high binding affinity to DnaK evolve on average 4.3-fold faster than proteins in the lowest binding affinity class at the genus resolution. Differences in evolutionary rates of DnaK interactor classes are still significant after adjusting for possible effects caused by protein expression level. Furthermore, we observe an additive effect of DnaK and GroEL chaperones on the evolutionary rates of their common interactors. Finally, we found pronounced similarities in the physicochemical profiles that characterize proteins belonging to DnaK and GroEL interactomes. Our results thus implicate DnaK-mediated folding as a major component in shaping protein evolutionary dynamics in bacteria and supply further evidence for the long-term manifestation of chaperone-mediated folding on genome evolution.

Transcriptomic profiling of host-parasite interactions in the microsporidian Trachipleistophora hominis

Background

Trachipleistophora hominis was isolated from an HIV/AIDS patient and is a member of a highly successful group of obligate intracellular parasites.

Methods

Here we have investigated the evolution of the parasite and the interplay between host and parasite gene expression using transcriptomics of T. hominis-infected rabbit kidney cells.

Results

T. hominis has about 30 % more genes than small-genome microsporidians. Highly expressed genes include those involved in growth, replication, defence against oxidative stress, and a large fraction of uncharacterised genes. Chaperones are also highly expressed and may buffer the deleterious effects of the large number of non-synonymous mutations observed in essential T. hominis genes. Host expression suggests a general cellular shutdown upon infection, but ATP, amino sugar and nucleotide sugar production appear enhanced, potentially providing the parasite with substrates it cannot make itself. Expression divergence of duplicated genes, including transporters used to acquire host metabolites, demonstrates ongoing functional diversification during microsporidian evolution. We identified overlapping transcription at more than 100 loci in the sparse T. hominis genome, demonstrating that this feature is not caused by genome compaction. The detection of additional transposons of insect origin strongly suggests that the natural host for T. hominis is an insect.

Conclusions

Our results reveal that the evolution of contemporary microsporidian genomes is highly dynamic and innovative. Moreover, highly expressed T. hominis genes of unknown function include a cohort that are shared among all microsporidians, indicating that some strongly conserved features of the biology of these enormously successful parasites remain uncharacterised.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1989-z) contains supplementary material, which is available to authorized users.

The GroEL-GroES Chaperonin Machine: A Nano-Cage for Protein Folding

The bacterial chaperonin GroEL and its cofactor GroES constitute the paradigmatic molecular machine of protein folding. GroEL is a large double-ring cylinder with ATPase activity that binds non-native substrate protein (SP) via hydrophobic residues exposed towards the ring center. Binding of the lid-shaped GroES to GroEL displaces the bound protein into an enlarged chamber, allowing folding to occur unimpaired by aggregation. GroES and SP undergo cycles of binding and release, regulated allosterically by the GroEL ATPase. Recent structural and functional studies are providing insights into how the physical environment of the chaperonin cage actively promotes protein folding, in addition to preventing aggregation. Here, we review different models of chaperonin action and discuss issues of current debate.

Fitness Trade-Offs Determine the Role of the Molecular Chaperonin GroEL in Buffering Mutations

Molecular chaperones fold many proteins and their mutated versions in a cell and can sometimes buffer the phenotypic effect of mutations that affect protein folding. Unanswered questions about this buffering include the nature of its mechanism, its influence on the genetic variation of a population, the fitness trade-offs constraining this mechanism, and its role in expediting evolution. Answering these questions is fundamental to understand the contribution of buffering to increase genetic variation and ecological diversification. Here, we performed experimental evolution, genome resequencing, and computational analyses to determine the trade-offs and evolutionary trajectories of Escherichia coli expressing high levels of the essential chaperonin GroEL. GroEL is abundantly present in bacteria, particularly in bacteria with large loads of deleterious mutations, suggesting its role in mutational buffering. We show that groEL overexpression is costly to large populations evolving in the laboratory, leading to groE expression decline within 66 generations. In contrast, populations evolving under the strong genetic drift characteristic of endosymbiotic bacteria avoid extinction or can be rescued in the presence of abundant GroEL. Genomes resequenced from cells evolved under strong genetic drift exhibited significantly higher tolerance to deleterious mutations at high GroEL levels than at native levels, revealing that GroEL is buffering mutations in these cells. GroEL buffered mutations in a highly diverse set of proteins that interact with the environment, including substrate and ion membrane transporters, hinting at its role in ecological diversification. Our results reveal the fitness trade-offs of mutational buffering and how genetic variation is maintained in populations.

Role of auxiliary proteins in Rubisco biogenesis and function

Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) catalyses the conversion of atmospheric CO2 into organic compounds during photosynthesis. Despite its pivotal role in plant metabolism, Rubisco is an inefficient enzyme and has therefore been a key target in bioengineering efforts to improve crop yields. Much has been learnt about the complex cellular machinery involved in Rubisco assembly and metabolic repair over recent years. The simple form of Rubisco found in certain bacteria and dinoflagellates comprises two large subunits, and generally requires the chaperonin system for folding. However, the evolution of hexadecameric Rubisco, which comprises eight large and eight small subunits, from its dimeric precursor has rendered Rubisco in most plants, algae, cyanobacteria and proteobacteria dependent on an array of additional factors. These auxiliary factors include several chaperones for assembly as well as ATPases of the AAA+ family for functional maintenance. An integrated view of the pathways underlying Rubisco biogenesis and repair will pave the way for efforts to improve the enzyme with the goal of increasing crop yields.

The origins of mutational robustness

Biological systems are resistant to genetic changes; a property known as mutational robustness, the origin of which remains an open question. In recent years, researchers have explored emergent properties of biological systems and mechanisms of genetic redundancy to reveal how mutational robustness emerges and persists. Several mechanisms have been proposed to explain the origin of mutational robustness, including molecular chaperones and gene duplication. The latter has received much attention, but its role in robustness remains controversial. Here, I examine recent findings linking genetic redundancy through gene duplication and mutational robustness. Experimental evolution and genome resequencing have made it possible to test the role of gene duplication in tolerating mutations at both the coding and regulatory levels. This evidence as well as previous findings on regulatory reprogramming of duplicates support the role of gene duplication in the origin of robustness.

RNA chaperones buffer deleterious mutations in E. coli

Both proteins and RNAs can misfold into non-functional conformations. Protein chaperones promote native folding of nascent polypeptides and refolding of misfolded species, thereby buffering mutations that compromise protein structure and function. Here, we show that RNA chaperones can also act as mutation buffers that enhance organismal fitness. Using competition assays, we demonstrate that overexpression of select RNA chaperones, including three DEAD box RNA helicases (DBRHs) (CsdA, SrmB, RhlB) and the cold shock protein CspA, improves fitness of two independently evolved Escherichia coli mutator strains that have accumulated deleterious mutations during short- and long-term laboratory evolution. We identify strain-specific mutations that are deleterious and subject to buffering when introduced individually into the ancestral genotype. For DBRHs, we show that buffering requires helicase activity, implicating RNA structural remodelling in the buffering process. Our results suggest that RNA chaperones might play a fundamental role in RNA evolution and evolvability.

Is catalytic activity of chaperones a selectable trait for the emergence of heat shock response?

Although heat shock response is ubiquitous in bacterial cells, the underlying physical chemistry behind heat shock response remains poorly understood. To study the response of cell populations to heat shock we employ a physics-based ab initio model of living cells where protein biophysics (i.e., folding and protein-protein interactions in crowded cellular environments) and important aspects of proteins homeostasis are coupled with realistic population dynamics simulations. By postulating a genotype-phenotype relationship we define a cell division rate in terms of functional concentrations of proteins and protein complexes, whose Boltzmann stabilities of folding and strengths of their functional interactions are exactly evaluated from their sequence information. We compare and contrast evolutionary dynamics for two models of chaperon action. In the active model, foldase chaperones function as nonequilibrium machines to accelerate the rate of protein folding. In the passive model, holdase chaperones form reversible complexes with proteins in their misfolded conformations to maintain their solubility. We find that only cells expressing foldase chaperones are capable of genuine heat shock response to the increase in the amount of unfolded proteins at elevated temperatures. In response to heat shock, cells' limited resources are redistributed differently for active and passive models. For the active model, foldase chaperones are overexpressed at the expense of downregulation of high abundance proteins, whereas for the passive model; cells react to heat shock by downregulating their high abundance proteins, as their low abundance proteins are upregulated.

Survival and innovation: The role of mutational robustness in evolution

Biological systems are resistant to perturbations caused by the environment and by the intrinsic noise of the system. Robustness to mutations is a particular aspect of robustness in which the phenotype is resistant to genotypic variation. Mutational robustness has been linked to the ability of the system to generate heritable genetic variation (a property known as evolvability). It is known that greater robustness leads to increased evolvability. Therefore, mechanisms that increase mutational robustness fuel evolvability. Two such mechanisms, molecular chaperones and gene duplication, have been credited with enormous importance in generating functional diversity through the increase of system's robustness to mutational insults. However, the way in which such mechanisms regulate robustness remains largely uncharacterized. In this review, I provide evidence in support of the role of molecular chaperones and gene duplication in innovation. Specifically, I present evidence that these mechanisms regulate robustness allowing unstable systems to survive long periods of time, and thus they provide opportunity for other mutations to compensate the destabilizing effects of functionally innovative mutations. The findings reported in this study set new questions with regards to the synergy between robustness mechanisms and how this synergy can alter the adaptive landscape of proteins. The ideas proposed in this article set the ground for future research in the understanding of the role of robustness in evolution.

Hsp90 promotes kinase evolution

Heat-shock protein 90 (Hsp90) promotes the maturation and stability of its client proteins, including many kinases. In doing so, Hsp90 may allow its clients to accumulate mutations as previously proposed by the capacitor hypothesis. If true, Hsp90 clients should show increased evolutionary rate compared with nonclients; however, other factors, such as gene expression and protein connectivity, may confound or obscure the chaperone’s putative contribution. Here, we compared the evolutionary rates of many Hsp90 clients and nonclients in the human protein kinase superfamily. We show that Hsp90 client status promotes evolutionary rate independently of, but in a small magnitude similar to that of gene expression and protein connectivity. Hsp90’s effect on kinase evolutionary rate was detected across mammals, specifically relaxing purifying selection. Hsp90 clients also showed increased nucleotide diversity and harbored more damaging variation than nonclient kinases across humans. These results are consistent with the central argument of the capacitor hypothesis that interaction with the chaperone allows its clients to harbor genetic variation. Hsp90 client status is thought to be highly dynamic with as few as one amino acid change rendering a protein dependent on the chaperone. Contrary to this expectation, we found that across protein kinase phylogeny Hsp90 client status tends to be gained, maintained, and shared among closely related kinases. We also infer that the ancestral protein kinase was not an Hsp90 client. Taken together, our results suggest that Hsp90 played an important role in shaping the kinase superfamily.

Reduction and expansion in microsporidian genome evolution: new insights from comparative genomics

Microsporidia are an abundant group of obligate intracellular parasites of other eukaryotes, including immunocompromised humans, but the molecular basis of their intracellular lifestyle and pathobiology are poorly understood. New genomes from a taxonomically broad range of microsporidians, complemented by published expression data, provide an opportunity for comparative analyses to identify conserved and lineage-specific patterns of microsporidian genome evolution that have underpinned this success. In this study, we infer that a dramatic bottleneck in the last common microsporidian ancestor (LCMA) left a small conserved core of genes that was subsequently embellished by gene family expansion driven by gene acquisition in different lineages. Novel expressed protein families represent a substantial fraction of sequenced microsporidian genomes and are significantly enriched for signals consistent with secretion or membrane location. Further evidence of selection is inferred from the gain and reciprocal loss of functional domains between paralogous genes, for example, affecting transport proteins. Gene expansions among transporter families preferentially affect those that are located on the plasma membrane of model organisms, consistent with recruitment to plug conserved gaps in microsporidian biosynthesis and metabolism. Core microsporidian genes shared with other eukaryotes are enriched in orthologs that, in yeast, are highly expressed, highly connected, and often essential, consistent with strong negative selection against further reduction of the conserved gene set since the LCMA. Our study reveals that microsporidian genome evolution is a highly dynamic process that has balanced constraint, reductive evolution, and genome expansion during adaptation to an extraordinarily successful obligate intracellular lifestyle.

A comprehensive analysis of the Omp85/TpsB protein superfamily structural diversity, taxonomic occurrence, and evolution

Members of the Omp85/TpsB protein superfamily are ubiquitously distributed in Gram-negative bacteria, and function in protein translocation (e.g., FhaC) or the assembly of outer membrane proteins (e.g., BamA). Several recent findings are suggestive of a further level of variation in the superfamily, including the identification of the novel membrane protein assembly factor TamA and protein translocase PlpD. To investigate the diversity and the causal evolutionary events, we undertook a comprehensive comparative sequence analysis of the Omp85/TpsB proteins. A total of 10 protein subfamilies were apparent, distinguished in their domain structure and sequence signatures. In addition to the proteins FhaC, BamA, and TamA, for which structural and functional information is available, are families of proteins with so far undescribed domain architectures linked to the Omp85 β-barrel domain. This study brings a classification structure to a dynamic protein superfamily of high interest given its essential function for Gram-negative bacteria as well as its diverse domain architecture, and we discuss several scenarios of putative functions of these so far undescribed proteins.

How do eubacterial organisms manage aggregation-prone proteome?

Eubacterial genomes vary considerably in their nucleotide composition. The percentage of genetic material constituted by guanosine and cytosine (GC) nucleotides ranges from 20% to 70%.  It has been posited that GC-poor organisms are more dependent on protein folding machinery. Previous studies have ascribed this to the accumulation of mildly deleterious mutations in these organisms due to population bottlenecks. This phenomenon has been supported by protein folding simulations, which showed that proteins encoded by GC-poor organisms are more prone to aggregation than proteins encoded by GC-rich organisms. To test this proposition using a genome-wide approach, we classified different eubacterial proteomes in terms of their aggregation propensity and chaperone-dependence using multiple machine learning models. In contrast to the expected decrease in protein aggregation with an increase in GC richness, we found that the aggregation propensity of proteomes increases with GC content. A similar and even more significant correlation was obtained with the GroEL-dependence of proteomes: GC-poor proteomes have evolved to be less dependent on GroEL than GC-rich proteomes. We thus propose that a decrease in eubacterial GC content may have been selected in organisms facing proteostasis problems.

Interplay between chaperones and protein disorder promotes the evolution of protein networks

Evolution is driven by mutations, which lead to new protein functions but come at a cost to protein stability. Non-conservative substitutions are of interest in this regard because they may most profoundly affect both function and stability. Accordingly, organisms must balance the benefit of accepting advantageous substitutions with the possible cost of deleterious effects on protein folding and stability. We here examine factors that systematically promote non-conservative mutations at the proteome level. Intrinsically disordered regions in proteins play pivotal roles in protein interactions, but many questions regarding their evolution remain unanswered. Similarly, whether and how molecular chaperones, which have been shown to buffer destabilizing mutations in individual proteins, generally provide robustness during proteome evolution remains unclear. To this end, we introduce an evolutionary parameter λ that directly estimates the rate of non-conservative substitutions. Our analysis of λ in Escherichia coli, Saccharomyces cerevisiae, and Homo sapiens sequences reveals how co- and post-translationally acting chaperones differentially promote non-conservative substitutions in their substrates, likely through buffering of their destabilizing effects. We further find that λ serves well to quantify the evolution of intrinsically disordered proteins even though the unstructured, thus generally variable regions in proteins are often flanked by very conserved sequences. Crucially, we show that both intrinsically disordered proteins and highly re-wired proteins in protein interaction networks, which have evolved new interactions and functions, exhibit a higher λ at the expense of enhanced chaperone assistance. Our findings thus highlight an intricate interplay of molecular chaperones and protein disorder in the evolvability of protein networks. Our results illuminate the role of chaperones in enabling protein evolution, and underline the importance of the cellular context and integrated approaches for understanding proteome evolution. We feel that the development of λ may be a valuable addition to the toolbox applied to understand the molecular basis of evolution.

The C-terminal tails of the bacterial chaperonin GroEL stimulate protein folding by directly altering the conformation of a substrate protein

Many essential cellular proteins fold only with the assistance of chaperonin machines like the GroEL-GroES system of Escherichia coli. However, the mechanistic details of assisted protein folding by GroEL-GroES remain the subject of ongoing debate. We previously demonstrated that GroEL-GroES enhances the productive folding of a kinetically trapped substrate protein through unfolding, where both binding energy and the energy of ATP hydrolysis are used to disrupt the inhibitory misfolded states. Here, we show that the intrinsically disordered yet highly conserved C-terminal sequence of the GroEL subunits directly contributes to substrate protein unfolding. Interactions between the C terminus and the non-native substrate protein alter the binding position of the substrate protein on the GroEL apical surface. The C-terminal tails also impact the conformational state of the substrate protein during capture and encapsulation on the GroEL ring. Importantly, removal of the C termini results in slower overall folding, reducing the fraction of the substrate protein that commits quickly to a productive folding pathway and slowing several kinetically distinct folding transitions that occur inside the GroEL-GroES cavity. The conserved C-terminal tails of GroEL are thus important for protein folding from the beginning to the end of the chaperonin reaction cycle.

Integration of two ancestral chaperone systems into one: the evolution of eukaryotic molecular chaperones in light of eukaryogenesis

Eukaryotic genomes are mosaics of genes acquired from their prokaryotic ancestors, the eubacterial endosymbiont that gave rise to the mitochondrion and its archaebacterial host. Genomic footprints of the prokaryotic merger at the origin of eukaryotes are still discernable in eukaryotic genomes, where gene expression and function correlate with their prokaryotic ancestry. Molecular chaperones are essential in all domains of life as they assist the functional folding of their substrate proteins and protect the cell against the cytotoxic effects of protein misfolding. Eubacteria and archaebacteria code for slightly different chaperones, comprising distinct protein folding pathways. Here we study the evolution of the eukaryotic protein folding pathways following the endosymbiosis event. A phylogenetic analysis of all 64 chaperones encoded in the Saccharomyces cerevisiae genome revealed 25 chaperones of eubacterial ancestry, 11 of archaebacterial ancestry, 10 of ambiguous prokaryotic ancestry, and 18 that may represent eukaryotic innovations. Several chaperone families (e.g., Hsp90 and Prefoldin) trace their ancestry to only one prokaryote group, while others, such as Hsp40 and Hsp70, are of mixed ancestry, with members contributed from both prokaryotic ancestors. Analysis of the yeast chaperone–substrate interaction network revealed no preference for interaction between chaperones and substrates of the same origin. Our results suggest that the archaebacterial and eubacterial protein folding pathways have been reorganized and integrated into the present eukaryotic pathway. The highly integrated chaperone system of yeast is a manifestation of the central role of chaperone-mediated folding in maintaining cellular fitness. Most likely, both archaebacterial and eubacterial chaperone systems were essential at the very early stages of eukaryogenesis, and the retention of both may have offered new opportunities for expanding the scope of chaperone-mediated folding.

Diversity in the origins of proteostasis networks--a driver for protein function in evolution

Although the sequence of a protein largely determines its function, proteins can adopt different folding states in response to changes in the environment, some of which may be deleterious to the organism. All organisms--Bacteria, Archaea and Eukarya--have evolved a protein homeostasis, or proteostasis, network comprising chaperones and folding factors, degradation components, signalling pathways and specialized compartmentalized modules that manage protein folding in response to environmental stimuli and variation. Surveying the origins of proteostasis networks reveals that they have co-evolved with the proteome to regulate the physiological state of the cell, reflecting the unique stresses that different cells or organisms experience, and that they have a key role in driving evolution by closely managing the link between the phenotype and the genotype.