Catalysis of protein folding by chaperones accelerates evolutionary dynamics in adapting cell populations

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.

A Theoretical Framework on the Biology of Prion Diseases

Background

Prion diseases or TSE diseases are a group of neurodegenerative disorders that manifest in several forms in humans, such as Kuru disease, Creutzfeldt‒Jakob disease (CJD), Gerstmann-Sträussler-Scheinker syndrome (GSS) and fatal familial insomnia.

Objective

In this study, we propose a multimodular hypothesis of prion diseases. According to this hypothesis, a prion disease manifests because of the interaction of two genetic modules, such as the PRNP gene module and that of the gene or genes responsible for one or more chaperones, with one or some chemical module on whose structure the products of the genes or genetic modules interact.

Methods

This study presents the perspective that modular thinking can allow us to overcome conceptual obstacles in the understanding and interpretation of prion diseases.

Results and Discussion

The structure of the chemical module or modules is directly responsible for the folding or misfolding of the PrPC protein. The etiology of acquired prion diseases is explained based on this hypothesis. Hence, it has been proposed that (g) CJD involves the PRNP gene mutant and one or more mutant genes for one or more chaperone genes. In contrast, sCJD has one or more mutant chaperone genes. When does iCJD occur? Healthy individuals manifest acquired prion disease through contamination when infected with one or more mutant chaperones. The mutant chaperones interact with the prion protein, and PrPC is converted to its isoform PrPSc. In a recent study, there was a case of an individual with CJD after COVID-19 infection.

Conclusion

This case emphasizes the link between neuroinflammation and protein misfolding and provides proof that chemical module formation is a necessary condition for the manifestation of prion diseases.

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.

Could Small Heat Shock Protein HSP27 Be a First-Line Target for Preventing Protein Aggregation in Parkinson's Disease?

Small heat shock proteins (HSPs), such as HSP27, are ubiquitously expressed molecular chaperones and are essential for cellular homeostasis. The major functions of HSP27 include chaperoning misfolded or unfolded polypeptides and protecting cells from toxic stress. Dysregulation of stress proteins is associated with many human diseases including neurodegenerative diseases, such as Parkinson's disease (PD). PD is characterized by the presence of aggregates of α-synuclein in the central and peripheral nervous system, which induces the degeneration of dopaminergic neurons in the substantia nigra pars compacta (SNpc) and in the autonomic nervous system. Autonomic dysfunction is an important non-motor phenotype of PD, which includes cardiovascular dysregulation, among others. Nowadays, the therapies for PD focus on dopamine (DA) replacement. However, certain non-motor symptoms with a great impact on quality of life do not respond to dopaminergic drugs; therefore, the development and testing of new treatments for non-motor symptoms of PD remain a priority. Since small HSP27 was shown to prevent α-synuclein aggregation and cytotoxicity, this protein might constitute a suitable target to prevent or delay the motor and non-motor symptoms of PD. In the first part of our review, we focus on the cardiovascular dysregulation observed in PD patients. In the second part, we present data on the possible role of HSP27 in preventing the accumulation of amyloid fibrils and aggregated forms of α-synuclein. We also include our own studies, highlighting the possible protective cardiac effects induced by L-DOPA treatment through the enhancement of HSP27 levels and activity.

Weak Chemical Interactions That Drive Protein Evolution: Crowding, Sticking, and Quinary Structure in Folding and Function

In recent years, better instrumentation and greater computing power have enabled the imaging of elusive biomolecule dynamics in cells, driving many advances in understanding the chemical organization of biological systems. The focus of this Review is on interactions in the cell that affect both biomolecular stability and function and modulate them. The same protein or nucleic acid can behave differently depending on the time in the cell cycle, the location in a specific compartment, or the stresses acting on the cell. We describe in detail the crowding, sticking, and quinary structure in the cell and the current methods to quantify them both in vitro and in vivo. Finally, we discuss protein evolution in the cell in light of current biophysical evidence. We describe the factors that drive protein evolution and shape protein interaction networks. These interactions can significantly affect the free energy, ΔG, of marginally stable and low-population proteins and, due to epistasis, direct the evolutionary pathways in an organism. We finally conclude by providing an outlook on experiments to come and the possibility of collaborative evolutionary biology and biophysical efforts.

Thermophilic Adaptation in Prokaryotes Is Constrained by Metabolic Costs of Proteostasis

Prokaryotes evolved to thrive in an extremely diverse set of habitats, and their proteomes bear signatures of environmental conditions. Although correlations between amino acid usage and environmental temperature are well-documented, understanding of the mechanisms of thermal adaptation remains incomplete. Here, we couple the energetic costs of protein folding and protein homeostasis to build a microscopic model explaining both the overall amino acid composition and its temperature trends. Low biosynthesis costs lead to low diversity of physical interactions between amino acid residues, which in turn makes proteins less stable and drives up chaperone activity to maintain appropriate levels of folded, functional proteins. Assuming that the cost of chaperone activity is proportional to the fraction of unfolded client proteins, we simulated thermal adaptation of model proteins subject to minimization of the total cost of amino acid synthesis and chaperone activity. For the first time, we predicted both the proteome-average amino acid abundances and their temperature trends simultaneously, and found strong correlations between model predictions and 402 genomes of bacteria and archaea. The energetic constraint on protein evolution is more apparent in highly expressed proteins, selected by codon adaptation index. We found that in bacteria, highly expressed proteins are similar in composition to thermophilic ones, whereas in archaea no correlation between predicted expression level and thermostability was observed. At the same time, thermal adaptations of highly expressed proteins in bacteria and archaea are nearly identical, suggesting that universal energetic constraints prevail over the phylogenetic differences between these domains of life.

Host proteostasis modulates influenza evolution

Predicting and constraining RNA virus evolution require understanding the molecular factors that define the mutational landscape accessible to these pathogens. RNA viruses typically have high mutation rates, resulting in frequent production of protein variants with compromised biophysical properties. Their evolution is necessarily constrained by the consequent challenge to protein folding and function. We hypothesized that host proteostasis mechanisms may be significant determinants of the fitness of viral protein variants, serving as a critical force shaping viral evolution. Here, we test that hypothesis by propagating influenza in host cells displaying chemically-controlled, divergent proteostasis environments. We find that both the nature of selection on the influenza genome and the accessibility of specific mutational trajectories are significantly impacted by host proteostasis. These findings provide new insights into features of host–pathogen interactions that shape viral evolution, and into the potential design of host proteostasis-targeted antiviral therapeutics that are refractory to resistance.

Influenza viruses, commonly called flu, can evade our immune system and develop resistance to treatments by changing frequently. Specifically, mutations in their genome cause influenza proteins to change in ways that can help the virus evade our defences. However, these mutations come at a cost and can prevent the viral proteins from forming functional and stable three-dimensional shapes – a process known as protein folding – thereby hampering the virus’ ability to replicate.

In human cells, proteins called chaperones can help our other proteins fold properly. Influenza viruses do not have their own chaperones and, instead, hijack those of their host. Host chaperones are therefore crucial to the virus’ ability to replicate. However, until now, it was not known if host chaperones can influence how these viruses evolve.

Here, Phillips et al. used mammalian cells to study how host chaperones affect an evolving influenza population. First, cells were engineered to either have normal chaperone levels, elevated chaperone levels, or inactive chaperones. Next, the H3N2 influenza strain was grown in these different conditions for nearly 200 generations and sequenced to determine how the virus evolved in each distinctive host chaperone environment.

Phillips et al. discovered that host chaperones affect the rate at which mutations accumulate in the influenza population, and also the types of mutations in the influenza genome. For instance, when a chaperone called Hsp90 was inactivated, mutations became prevalent in the viral population more slowly than in cells with normal or elevated chaperone levels. Moreover, some specific mutations fared better in cells with high chaperone levels, whilst others worked better in cells with inactivated chaperones.

These results suggest that influenza evolution is affected by host chaperone levels in complex and important ways. Moreover, whether chaperones will promote or hinder the effects of any single mutation is difficult to predict ahead of time. This discovery is significant, as the chaperones available to influenza can vary in different tissues, organisms and infectious conditions, and may therefore influence the virus' ability to change and evolve in a context-specific manner.

The findings are likely to extend to other viruses such as HIV and Ebola, which also hijack host chaperones for the same purpose. More work is now needed to systematically quantify these effects so that we can better predict how specific chaperones will affect the ability of viruses to adapt, especially in pathologically relevant conditions like fever or viral host-switching. In the future, such insights could help shape the design of treatments to which viruses do not evolve resistance.

Bridging the physical scales in evolutionary biology: from protein sequence space to fitness of organisms and populations

Bridging the gap between the molecular properties of proteins and organismal/population fitness is essential for understanding evolutionary processes. This task requires the integration of the several physical scales of biological organization, each defined by a distinct set of mechanisms and constraints, into a single unifying model. The molecular scale is dominated by the constraints imposed by the physico-chemical properties of proteins and their substrates, which give rise to trade-offs and epistatic (non-additive) effects of mutations. At the systems scale, biological networks modulate protein expression and can either buffer or enhance the fitness effects of mutations. The population scale is influenced by the mutational input, selection regimes, and stochastic changes affecting the size and structure of populations, which eventually determine the evolutionary fate of mutations. Here, we summarize the recent advances in theory, computer simulations, and experiments that advance our understanding of the links between various physical scales in biology.

Massively parallel sampling of lattice proteins reveals foundations of thermal adaptation

Evolution of proteins in bacteria and archaea living in different conditions leads to significant correlations between amino acid usage and environmental temperature. The origins of these correlations are poorly understood, and an important question of protein theory, physics-based prediction of types of amino acids overrepresented in highly thermostable proteins, remains largely unsolved. Here, we extend the random energy model of protein folding by weighting the interaction energies of amino acids by their frequencies in protein sequences and predict the energy gap of proteins designed to fold well at elevated temperatures. To test the model, we present a novel scalable algorithm for simultaneous energy calculation for many sequences in many structures, targeting massively parallel computing architectures such as graphics processing unit. The energy calculation is performed by multiplying two matrices, one representing the complete set of sequences, and the other describing the contact maps of all structural templates. An implementation of the algorithm for the CUDA platform is available at http://www.github.com/kzeldovich/galeprot and calculates protein folding energies over 250 times faster than a single central processing unit. Analysis of amino acid usage in 64-mer cubic lattice proteins designed to fold well at different temperatures demonstrates an excellent agreement between theoretical and simulated values of energy gap. The theoretical predictions of temperature trends of amino acid frequencies are significantly correlated with bioinformatics data on 191 bacteria and archaea, and highlight protein folding constraints as a fundamental selection pressure during thermal adaptation in biological evolution.

Regulation of human 3-beta-hydroxysteroid dehydrogenase type-2 (3βHSD2) by molecular chaperones and the mitochondrial environment affects steroidogenesis

Human 3-β-hydroxysteroid dehydrogenase/isomerase types 1 and 2 (3βHSD1 and 3βHSD2, respectively) are expressed in a tissue-specific pattern by different genes. Site-directed mutagenesis studies have confirmed the function of the catalytic amino acids (Tyr154, Lys 158, Ser124 in both isoenzymes), substrate/inhibitor isoform-specific residues (His156 and Arg195 in 3βHSD1) and cofactor binding residues (Asp36 provides NAD(+) specificity in both isoenzymes). However, detailed analysis of isoform-specific organelle localization and characterization is difficult due to the 93% amino acid identity between the two isoforms. With recent advances in the knowledge of mitochondrial architecture and localization of the various translocases, our laboratory has studied the mechanisms regulating mitochondrial 3βHSD2 localization. The mitochondrial N-terminal leader sequence of 3βHSD2 directs its entry into the mitochondria where it is localized to the intermembrane space (IMS). Unlike other mitochondrial proteins, the N-terminal signal sequence of 3βHSD2 is not cleaved upon mitochondrial import. 3βHSD2 interacts with the mitochondrial translocase, Tim50, to regulate progesterone and androstenedione formation. Our studies suggest that its activity at the IMS is facilitated in a partially unfolded "molten globule" conformation by the proton pump between the matrix and IMS. The unfolded protein is refolded by the mitochondrial chaperones. The protons at the IMS are absorbed by the lipid vesicles, to maintain the proton pump and recycle 3βHSD2. As a result, one molecule of 3βHSD2 may participate in multiple catalytic reactions. In summary, the steroidogenic cell recycles 3βHSD2 to catalyze the reactions needed to produce androstenedione, progesterone and 17α-hydroxyprogesterone on demand in coordination with the mitochondrial translocase, Tim50. This article is part of a Special Issue entitled 'Steroid/Sterol signaling'.

Evolution of specificity in protein-protein interactions

Hub proteins are proteins that maintain promiscuous molecular recognition. Because they are reported to play essential roles in cellular control, there has been a special interest in the study of their structural and functional properties, yet the mechanisms by which they evolve to maintain functional interactions are poorly understood. By combining biophysical simulations of coarse-grained proteins and analysis of proteins-complex crystallographic structures, we seek to elucidate those mechanisms. We focus on two types of hub proteins: Multi hubs, which interact with their partners through different interfaces, and Singlish hubs, which do so through a single interface. We show that loss of structural stability is required for the evolution of protein-protein-interaction (PPI) networks, and it is more profound in Singlish hub systems. In addition, different ratios of hydrophobic to electrostatic interfacial amino acids are shown to support distinct network topologies (i.e., Singlish and Multi systems), and therefore underlie a fundamental design principle of PPI in a crowded environment. We argue that the physical nature of hydrophobic and electrostatic interactions, in particular, their favoring of either same-type interactions (hydrophobic-hydrophobic), or opposite-type interactions (negatively-positively charged) plays a key role in maintaining the network topology while allowing the protein amino acid sequence to evolve.

Biophysics of protein evolution and evolutionary protein biophysics

The study of molecular evolution at the level of protein-coding genes often entails comparing large datasets of sequences to infer their evolutionary relationships. Despite the importance of a protein's structure and conformational dynamics to its function and thus its fitness, common phylogenetic methods embody minimal biophysical knowledge of proteins. To underscore the biophysical constraints on natural selection, we survey effects of protein mutations, highlighting the physical basis for marginal stability of natural globular proteins and how requirement for kinetic stability and avoidance of misfolding and misinteractions might have affected protein evolution. The biophysical underpinnings of these effects have been addressed by models with an explicit coarse-grained spatial representation of the polypeptide chain. Sequence-structure mappings based on such models are powerful conceptual tools that rationalize mutational robustness, evolvability, epistasis, promiscuous function performed by 'hidden' conformational states, resolution of adaptive conflicts and conformational switches in the evolution from one protein fold to another. Recently, protein biophysics has been applied to derive more accurate evolutionary accounts of sequence data. Methods have also been developed to exploit sequence-based evolutionary information to predict biophysical behaviours of proteins. The success of these approaches demonstrates a deep synergy between the fields of protein biophysics and protein evolution.

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.

HSP90 empowers evolution of resistance to hormonal therapy in human breast cancer models

The efficacy of hormonal therapies for advanced estrogen receptor-positive breast cancers is limited by the nearly inevitable development of acquired resistance. Efforts to block the emergence of resistance have met with limited success, largely because the mechanisms underlying it are so varied and complex. Here, we investigate a new strategy aimed at the very processes by which cancers evolve resistance. From yeast to vertebrates, heat shock protein 90 (HSP90) plays a unique role among molecular chaperones by promoting the evolution of heritable new traits. It does so by regulating the folding of a diverse portfolio of metastable client proteins, many of which mediate adaptive responses that allow organisms to adapt and thrive in the face of diverse challenges, including those posed by drugs. Guided by our previous work in pathogenic fungi, in which very modest HSP90 inhibition impairs resistance to mechanistically diverse antifungals, we examined the effect of similarly modest HSP90 inhibition on the emergence of resistance to antiestrogens in breast cancer models. Even though this degree of inhibition fell below the threshold for proteotoxic activation of the heat-shock response and had no overt anticancer activity on its own, it dramatically impaired the emergence of resistance to hormone antagonists both in cell culture and in mice. Our findings strongly support the clinical testing of combined hormone antagonist-low-level HSP90 inhibitor regimens in the treatment of metastatic estrogen receptor-positive breast cancer. At a broader level, they also provide promising proof of principle for a generalizable strategy to combat the pervasive problem of rapidly emerging resistance to molecularly targeted therapeutics.

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.