Cotranslational folding allows misfolding-prone proteins to circumvent deep kinetic traps

Translation Rates and Protein Folding

The mRNA coding sequence defines not only the amino acid sequence of the protein, but also the speed at which the ribosomes move along the mRNA while making the protein. The non-uniform local kinetics - denoted as translational rhythm - is similar among mRNAs coding for related protein folds. Deviations from this conserved rhythm can result in protein misfolding. In this review we summarize the experimental evidence demonstrating how local translation rates affect cotranslational protein folding, with the focus on the synonymous codons and patches of charged residues in the nascent peptide as best-studied examples. Alterations in nascent protein conformations due to disturbed translational rhythm can persist off the ribosome, as demonstrated by the effects of synonymous codon variants of several disease-related proteins. Charged amino acid patches in nascent chains also modulate translation and cotranslational protein folding, and can abrogate translation when placed at the N-terminus of the nascent peptide. During cotranslational folding, incomplete nascent chains navigate through a unique conformational landscape in which earlier intermediate states become inaccessible as the nascent peptide grows. Precisely tuned local translation rates, as well as interactions with the ribosome, guide the folding pathway towards the native structure, whereas deviations from the natural translation rhythm may favor pathways leading to trapped misfolded states. Deciphering the 'folding code' of the mRNA will contribute to understanding the diseases caused by protein misfolding and to rational protein design.

Impact of Nonnative Interactions on the Binding Kinetics of Intrinsically Disordered p53 with MDM2: Insights from All-Atom Simulation and Markov State Model Analysis

Intrinsically disordered proteins (IDPs) lack a well-defined tertiary structure but are essential players in various biological processes. Their ability to undergo a disorder-to-order transition upon binding to their partners, known as the folding-upon-binding process, is crucial for their function. One classical example is the intrinsically disordered transactivation domain (TAD) of the tumor suppressor protein p53, which quickly forms a structured α-helix after binding to its partner MDM2, with clinical significance for cancer treatment. However, the contribution of nonnative interactions between the IDP and its partner to the rapid binding kinetics, as well as their interplay with native interactions, is not well understood at the atomic level. Here, we used molecular dynamics simulation and Markov state model (MSM) analysis to study the folding-upon-binding mechanism between p53-TAD and MDM2. Our results suggest that the system progresses from the nascent encounter complex to the well-structured encounter complex and finally reaches the native complex, following an induced-fit mechanism. We found that nonnative hydrophobic and hydrogen bond interactions, combined with native interactions, effectively stabilize the nascent and well-structured encounter complexes. Among the nonnative interactions, Leu25p53-Leu54MDM2 and Leu25p53-Phe55MDM2 are particularly noteworthy, as their interaction strength is close to the optimum. Evidently, strengthening or weakening these interactions could both adversely affect the binding kinetics. Overall, our findings suggest that nonnative interactions are evolutionarily optimized to accelerate the binding kinetics of IDPs in conjunction with native interactions.

How the ribosome shapes cotranslational protein folding

During protein synthesis, the growing nascent peptide chain moves inside the polypeptide exit tunnel of the ribosome from the peptidyl transferase center towards the exit port where it emerges into the cytoplasm. The ribosome defines the unique energy landscape of the pioneering round of protein folding. The spatial confinement and the interactions of the nascent peptide with the tunnel walls facilitate formation of secondary structures, such as α-helices. The vectorial nature of protein folding inside the tunnel favors local intra- and inter-molecular interactions, thereby inducing cotranslational folding intermediates that do not form upon protein refolding in solution. Tertiary structures start to fold in the lower part of the tunnel, where interactions with the ribosome destabilize native protein folds. The present review summarizes the recent progress in understanding the driving forces of nascent protein folding inside the tunnel and at the surface of the ribosome.

Nuclear Hsp104 safeguards the dormant translation machinery during quiescence

The resilience of cellular proteostasis declines with age, which drives protein aggregation and compromises viability. The nucleus has emerged as a key quality control compartment that handles misfolded proteins produced by the cytosolic protein biosynthesis system. Here, we find that age-associated metabolic cues target the yeast protein disaggregase Hsp104 to the nucleus to maintain a functional nuclear proteome during quiescence. The switch to respiratory metabolism and the accompanying decrease in translation rates direct cytosolic Hsp104 to the nucleus to interact with latent translation initiation factor eIF2 and to suppress protein aggregation. Hindering Hsp104 from entering the nucleus in quiescent cells results in delayed re-entry into the cell cycle due to compromised resumption of protein synthesis. In sum, we report that cytosolic-nuclear partitioning of the Hsp104 disaggregase is a critical mechanism to protect the latent protein synthesis machinery during quiescence in yeast, ensuring the rapid restart of translation once nutrients are replenished.

Mechanisms and pathology of protein misfolding and aggregation

Despite advances in machine learning-based protein structure prediction, we are still far from fully understanding how proteins fold into their native conformation. The conventional notion that polypeptides fold spontaneously to their biologically active states has gradually been replaced by our understanding that cellular protein folding often requires context-dependent guidance from molecular chaperones in order to avoid misfolding. Misfolded proteins can aggregate into larger structures, such as amyloid fibrils, which perpetuate the misfolding process, creating a self-reinforcing cascade. A surge in amyloid fibril structures has deepened our comprehension of how a single polypeptide sequence can exhibit multiple amyloid conformations, known as polymorphism. The assembly of these polymorphs is not a random process but is influenced by the specific conditions and tissues in which they originate. This observation suggests that, similar to the folding of native proteins, the kinetics of pathological amyloid assembly are modulated by interactions specific to cells and tissues. Here, we review the current understanding of how intrinsic protein conformational propensities are modulated by physiological and pathological interactions in the cell to shape protein misfolding and aggregation pathology.

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.

Co-translational formation of disulfides guides folding of the SARS-CoV-2 receptor binding domain

Many secreted proteins, including viral proteins, contain multiple disulfide bonds. How disulfide formation is coupled to protein folding in the cell remains poorly understood at the molecular level. Here, we combine experiment and simulation to address this question as it pertains to the SARS-CoV-2 receptor binding domain (RBD). We show that the RBD can only refold reversibly if its native disulfides are present before folding. But in their absence, the RBD spontaneously misfolds into a nonnative, molten-globule-like state that is structurally incompatible with complete disulfide formation and that is highly prone to aggregation. Thus, the RBD native structure represents a metastable state on the protein's energy landscape with reduced disulfides, indicating that nonequilibrium mechanisms are needed to ensure native disulfides form before folding. Our atomistic simulations suggest that this may be achieved via co-translational folding during RBD secretion into the endoplasmic reticulum. Namely, at intermediate translation lengths, native disulfide pairs are predicted to come together with high probability, and thus, under suitable kinetic conditions, this process may lock the protein into its native state and circumvent highly aggregation-prone nonnative intermediates. This detailed molecular picture of the RBD folding landscape may shed light on SARS-CoV-2 pathology and molecular constraints governing SARS-CoV-2 evolution.

Chromosome and protein folding: In search for unified principles

Structural biology has traditionally focused on the structures of proteins, short nucleic acids, small molecules, and their complexes. However, it is now widely recognized that the 3D organization of chromosomes should also be included in this list, despite significant differences in scale and complexity of organization. Here we highlight some notable similarities between the folding processes that shape proteins and chromosomes. Both biomolecules are folded by two types of processes: the affinity-mediated interactions, and by active (ATP-dependent) processes. Both chromosome and proteins in vivo can have partially unstructured and non-equilibrium ensembles with yet to be understood functional roles. By analyzing these biological systems in parallel, we can uncover universal principles of biomolecular organization that transcend specific biopolymers.

Learning to learn by using nonequilibrium training protocols for adaptable materials

Evolution in time-varying environments naturally leads to adaptable biological systems that can easily switch functionalities. Advances in the synthesis of environmentally responsive materials therefore open up the possibility of creating a wide range of synthetic materials which can also be trained for adaptability. We consider high-dimensional inverse problems for materials where any particular functionality can be realized by numerous equivalent choices of design parameters. By periodically switching targets in a given design algorithm, we can teach a material to perform incompatible functionalities with minimal changes in design parameters. We exhibit this learning strategy for adaptability in two simulated settings: elastic networks that are designed to switch deformation modes with minimal bond changes and heteropolymers whose folding pathway selections are controlled by a minimal set of monomer affinities. The resulting designs can reveal physical principles, such as nucleation-controlled folding, that enable such adaptability.

How synonymous mutations alter enzyme structure and function over long timescales

The specific activity of enzymes can be altered over long timescales in cells by synonymous mutations that alter a messenger RNA molecule's sequence but not the encoded protein's primary structure. How this happens at the molecular level is unknown. Here, we use multiscale modelling of three Escherichia coli enzymes (type III chloramphenicol acetyltransferase, D-alanine-D-alanine ligase B and dihydrofolate reductase) to understand experimentally measured changes in specific activity due to synonymous mutations. The modelling involves coarse-grained simulations of protein synthesis and post-translational behaviour, all-atom simulations to test robustness and quantum mechanics/molecular mechanics calculations to characterize enzymatic function. We show that changes in codon translation rates induced by synonymous mutations cause shifts in co-translational and post-translational folding pathways that kinetically partition molecules into subpopulations that very slowly interconvert to the native, functional state. Structurally, these states resemble the native state, with localized misfolding near the active sites of the enzymes. These long-lived states exhibit reduced catalytic activity, as shown by their increased activation energies for the reactions they catalyse.

Cotranslational formation of disulfides guides folding of the SARS COV-2 receptor binding domain

Many secreted proteins contain multiple disulfide bonds. How disulfide formation is coupled to protein folding in the cell remains poorly understood at the molecular level. Here, we combine experiment and simulation to address this question as it pertains to the SARS-CoV-2 receptor binding domain (RBD). We show that, whereas RBD can refold reversibly when its disulfides are intact, their disruption causes misfolding into a nonnative molten-globule state that is highly prone to aggregation and disulfide scrambling. Thus, non-equilibrium mechanisms are needed to ensure disulfides form prior to folding in vivo. Our simulations suggest that co-translational folding may accomplish this, as native disulfide pairs are predicted to form with high probability at intermediate lengths, ultimately committing the RBD to its metastable native state and circumventing nonnative intermediates. This detailed molecular picture of the RBD folding landscape may shed light on SARS-CoV-2 pathology and molecular constraints governing SARS-CoV-2 evolution.

Universal protein misfolding intermediates can bypass the proteostasis network and remain soluble and less functional

Some misfolded protein conformations can bypass proteostasis machinery and remain soluble in vivo. This is an unexpected observation, as cellular quality control mechanisms should remove misfolded proteins. Three questions, then, are: how do long-lived, soluble, misfolded proteins bypass proteostasis? How widespread are such misfolded states? And how long do they persist? We address these questions using coarse-grain molecular dynamics simulations of the synthesis, termination, and post-translational dynamics of a representative set of cytosolic E. coli proteins. We predict that half of proteins exhibit misfolded subpopulations that bypass molecular chaperones, avoid aggregation, and will not be rapidly degraded, with some misfolded states persisting for months or longer. The surface properties of these misfolded states are native-like, suggesting they will remain soluble, while self-entanglements make them long-lived kinetic traps. In terms of function, we predict that one-third of proteins can misfold into soluble less-functional states. For the heavily entangled protein glycerol-3-phosphate dehydrogenase, limited-proteolysis mass spectrometry experiments interrogating misfolded conformations of the protein are consistent with the structural changes predicted by our simulations. These results therefore provide an explanation for how proteins can misfold into soluble conformations with reduced functionality that can bypass proteostasis, and indicate, unexpectedly, this may be a wide-spread phenomenon.

Co-Translational Folding of Multi-Domain Proteins

The majority of proteins in nature are composed of multiple domains connected in a single polypeptide. How these long sequences fold into functional structures without forming toxic misfolds or aggregates is poorly understood. Their folding is inextricably linked to protein synthesis and interactions with cellular machinery, making mechanistic studies challenging. Recent progress has revealed critical features of multi-domain protein folding in isolation and in the context of translation by the ribosome. In this review, we discuss challenges and progress in understanding multi-domain protein folding, and highlight how molecular interactions shape folding and misfolding pathways. With the development of new approaches and model systems, the stage is now set for mechanistically exploring the folding of large multi-domain proteins.

Thermodynamics of co-translational folding and ribosome-nascent chain interactions

Proteins can begin the conformational search for their native structure in parallel with biosynthesis on the ribosome, in a process termed co-translational folding. In contrast to the reversible folding of isolated domains, as a nascent chain emerges from the ribosome exit tunnel during translation the free energy landscape it explores also evolves as a function of chain length. While this presents a substantially more complex measurement problem, this review will outline the progress that has been made recently in understanding, quantitatively, the process by which a nascent chain attains its full native stability, as well as the mechanisms through which interactions with the nearby ribosome surface can perturb or modulate this process.

Protein folding in vitro and in the cell: From a solitary journey to a team effort

Correct protein folding is essential for the health and function of living organisms. Yet, it is not well understood how unfolded proteins reach their native state and avoid aggregation, especially within the cellular milieu. Some proteins, especially small, single-domain and apparent two-state folders, successfully attain their native state upon dilution from denaturant. Yet, many more proteins undergo misfolding and aggregation during this process, in a concentration-dependent fashion. Once formed, native and aggregated states are often kinetically trapped relative to each other. Hence, the early stages of protein life are absolutely critical for proper kinetic channeling to the folded state and for long-term solubility and function. This review summarizes current knowledge on protein folding/aggregation mechanisms in buffered solution and within the bacterial cell, highlighting early stages. Remarkably, teamwork between nascent chain, ribosome, trigger factor and Hsp70 molecular chaperones enables all proteins to overcome aggregation propensities and reach a long-lived bioactive state.

A small molecule high throughput screening platform to profile conformational properties of nascent, ribosome-bound proteins

Genetic mutations cause a wide spectrum of human disease by disrupting protein folding, both during and after synthesis. Transient de-novo folding intermediates therefore represent potential drug targets for pharmacological correction of protein folding disorders. Here we develop a FRET-based high-throughput screening (HTS) assay in 1,536-well format capable of identifying small molecules that interact with nascent polypeptides and correct genetic, cotranslational folding defects. Ribosome nascent chain complexes (RNCs) containing donor and acceptor fluorophores were isolated from cell free translation reactions, immobilized on Nickel-NTA/IDA beads, and imaged by high-content microscopy. Quantitative FRET measurements obtained from as little as 0.4 attomole of protein/bead enabled rapid assessment of conformational changes with a high degree of reproducibility. Using this assay, we performed a pilot screen of ~ 50,000 small molecules to identify compounds that interact with RNCs containing the first nucleotide-binding domain (NBD1) of the cystic fibrosis transmembrane conductance regulator (CFTR) harboring a disease-causing mutation (A455E). Screen results yielded 133 primary hits and 1 validated hit that normalized FRET values of the mutant nascent peptide. This system provides a scalable, tractable, structure-based discovery platform for screening small molecules that bind to or impact the folding of protein substrates that are not amenable to traditional biochemical analyses.

Towards a generic prototyping approach for therapeutically-relevant peptides and proteins in a cell-free translation system

Advances in peptide and protein therapeutics increased the need for rapid and cost-effective polypeptide prototyping. While in vitro translation systems are well suited for fast and multiplexed polypeptide prototyping, they suffer from misfolding, aggregation and disulfide-bond scrambling of the translated products. Here we propose that efficient folding of in vitro produced disulfide-rich peptides and proteins can be achieved if performed in an aggregation-free and thermodynamically controlled folding environment. To this end, we modify an E. coli-based in vitro translation system to allow co-translational capture of translated products by affinity matrix. This process reduces protein aggregation and enables productive oxidative folding and recycling of misfolded states under thermodynamic control. In this study we show that the developed approach is likely to be generally applicable for prototyping of a wide variety of disulfide-constrained peptides, macrocyclic peptides with non-native bonds and antibody fragments in amounts sufficient for interaction analysis and biological activity assessment.

Generic approach for rapid prototyping is essential for the progress of synthetic biology. Here the authors modify the cell-free translation system to control protein aggregation and folding and validate the approach by using single conditions for prototyping of various disulfide-constrained polypeptides.

Slowest-first protein translation scheme: Structural asymmetry and co-translational folding

Proteins are translated from the N to the C terminus, raising the basic question of how this innate directionality affects their evolution. To explore this question, we analyze 16,200 structures from the Protein Data Bank (PDB). We find remarkable enrichment of α helices at the C terminus and β strands at the N terminus. Furthermore, this α-β asymmetry correlates with sequence length and contact order, both determinants of folding rate, hinting at possible links to co-translational folding (CTF). Hence, we propose the "slowest-first" scheme, whereby protein sequences evolved structural asymmetry to accelerate CTF: the slowest of the cooperatively folding segments are positioned near the N terminus so they have more time to fold during translation. A phenomenological model predicts that CTF can be accelerated by asymmetry in folding rate, up to double the rate, when folding time is commensurate with translation time; analysis of the PDB predicts that structural asymmetry is indeed maximal in this regime. This correspondence is greater in prokaryotes, which generally require faster protein production. Altogether, this indicates that accelerating CTF is a substantial evolutionary force whose interplay with stability and functionality is encoded in secondary structure asymmetry.

Folding Intermediates, Heterogeneous Native Ensembles and Protein Function

Single domain proteins fold via diverse mechanisms emphasizing the intricate relationship between energetics and structure, which is a direct consequence of functional constraints and demands imposed at the level of sequence. On the other hand, elucidating the interplay between folding mechanisms and function is challenging in large proteins, given the inherent shortcomings in identifying metastable states experimentally and the sampling limitations associated with computational methods. Here, we show that free energy profiles and surfaces of large systems (>150 residues), as predicted by a statistical mechanical model, display a wide array of folding mechanisms with ubiquitous folding intermediates and heterogeneous native ensembles. Importantly, residues around the ligand binding or enzyme active site display a larger tendency to partially unfold and this manifests as intermediates or excited states along the folding coordinate in ligand binding domains, transcription repressors, and representative enzymes from all the six classes, including the SARS-CoV-2 receptor binding domain (RBD) of the spike protein and the protease M^pro. It thus appears that it is relatively easier to distill the imprints of function on the folding landscape of larger proteins as opposed to smaller systems. We discuss how an understanding of energetic-entropic features in ordered proteins can pinpoint specific avenues through which folding mechanisms, populations of partially structured states and function can be engineered.

Ribosome Elongation Kinetics of Consecutively Charged Residues Are Coupled to Electrostatic Force

The speed of protein synthesis can dramatically change when consecutively charged residues are incorporated into an elongating nascent protein by the ribosome. The molecular origins of this class of allosteric coupling remain unknown. We demonstrate, using multiscale simulations, that positively charged residues generate large forces that move the P-site amino acid away from the A-site amino acid. Negatively charged residues generate forces of similar magnitude but move the A- and P-sites closer together. These conformational changes, respectively, increase and decrease the transition state barrier height to peptide bond formation, explaining how charged residues mechanochemically alter translation speed. This mechanochemical mechanism is consistent with in vivo ribosome profiling data exhibiting proportionality between translation speed and the number of charged residues, experimental data characterizing nascent chain conformations, and a previously published cryo-EM structure of a ribosome-nascent chain complex containing consecutive lysines. These results expand the role of mechanochemistry in translation and provide a framework for interpreting experimental results on translation speed.

The Protein Folding Problem: The Role of Theory

The protein folding problem was first articulated as question of how order arose from disorder in proteins: How did the various native structures of proteins arise from interatomic driving forces encoded within their amino acid sequences, and how did they fold so fast? These matters have now been largely resolved by theory and statistical mechanics combined with experiments. There are general principles. Chain randomness is overcome by solvation-based codes. And in the needle-in-a-haystack metaphor, native states are found efficiently because protein haystacks (conformational ensembles) are funnel-shaped. Order-disorder theory has now grown to encompass a large swath of protein physical science across biology.

Effects of electrostatic interactions on global folding and local conformational dynamics of a multidomain Y-family DNA polymerase

The Y-family DNA polymerases specialize in translesion DNA synthesis, which is essential for replicating damaged DNA. The Y-family polymerases, which are made up of four stable domains, exhibit extensive distributions of charged residues, and are responsible for the tight formation of the protein-DNA complex. However, it is still unclear how the electrostatic interactions influence the conformational dynamics of the polymerases. Here, we focus on the case of a prototype Y-family DNA polymerase, Dpo4. Using coarse-grained models including a salt-dependent electrostatic potential, we investigate the effects of the electrostatic interactions on the folding process of Dpo4. Our simulations show that strong electrostatic interactions result in a three-state folding of Dpo4, consistent with the experimental observations. This folding process exhibits low cooperativity led by low salt concentration, where the individual domains fold one by one through one single pathway. Since the refined folding order of domains in multidomain proteins can shrink the configurational space, we suggest that the electrostatic interactions facilitate the Dpo4 folding. In addition, we study the local conformational dynamics of Dpo4 in terms of fluctuation and frustration analyses. We show that the electrostatic interactions can exaggerate the local conformational properties, which are in favor of the large-scale conformational transition of Dpo4 during the functional DNA binding. Our results underline the importance of electrostatic interactions in the conformational dynamics of Dpo4 at both the global and local scale, providing useful guidance in protein engineering at the multidomain level.

Protein structure and aggregation: a marriage of necessity ruled by aggregation gatekeepers

Protein aggregation propensity is a pervasive and seemingly inescapable property of proteomes. Strikingly, a significant fraction of the proteome is supersaturated, meaning that, for these proteins, their native conformation is less stable than the aggregated state. Maintaining the integrity of a proteome under such conditions is precarious and requires energy-consuming proteostatic regulation. Why then is aggregation propensity maintained at such high levels over long evolutionary timescales? Here, we argue that the conformational stability of the native and aggregated states are correlated thermodynamically and that codon usage strengthens this correlation. As a result, the folding of stable proteins requires kinetic control to avoid aggregation, provided by aggregation gatekeepers. These unique residues are evolutionarily selected to kinetically favor native folding, either on their own or by coopting chaperones.

Mechanisms of Cotranslational Protein Maturation in Bacteria

Growing cells invest a significant part of their biosynthetic capacity into the production of proteins. To become functional, newly-synthesized proteins must be N-terminally processed, folded and often translocated to other cellular compartments. A general strategy is to integrate these protein maturation processes with translation, by cotranslationally engaging processing enzymes, chaperones and targeting factors with the nascent polypeptide. Precise coordination of all factors involved is critical for the efficiency and accuracy of protein synthesis and cellular homeostasis. This review provides an overview of the current knowledge on cotranslational protein maturation, with a focus on the production of cytosolic proteins in bacteria. We describe the role of the ribosome and the chaperone network in protein folding and how the dynamic interplay of all cotranslationally acting factors guides the sequence of cotranslational events. Finally, we discuss recent data demonstrating the coupling of protein synthesis with the assembly of protein complexes and end with a brief discussion of outstanding questions and emerging concepts in the field of cotranslational protein maturation.

The effects of codon bias and optimality on mRNA and protein regulation

The central dogma of molecular biology entails that genetic information is transferred from nucleic acid to proteins. Notwithstanding retro-transcribing genetic elements, DNA is transcribed to RNA which in turn is translated into proteins. Recent advancements have shown that each stage is regulated to control protein abundances for a variety of essential physiological processes. In this regard, mRNA regulation is essential in fine-tuning or calibrating protein abundances. In this review, we would like to discuss one of several mRNA-intrinsic features of mRNA regulation that has been gaining traction of recent-codon bias and optimality. Specifically, we address the effects of codon bias with regard to codon optimality in several biological processes centred on translation, such as mRNA stability and protein folding among others. Finally, we examine how different organisms or cell types, through this system, are able to coordinate physiological pathways to respond to a variety of stress or growth conditions.

Energy-dependent protein folding: modeling how a protein folding machine may work

Background: Proteins fold robustly and reproducibly in vivo, but many cannot fold in vitro in isolation from cellular components. Despite the remarkable progress that has been achieved by the artificial intelligence approaches in predicting the protein native conformations, the pathways that lead to such conformations, either in vitro or in vivo, remain largely unknown. The slow progress in recapitulating protein folding pathways in silico may be an indication of the fundamental deficiencies in our understanding of folding as it occurs in nature. Here we consider the possibility that protein folding in living cells may not be driven solely by the decrease in Gibbs free energy and propose that protein folding in vivo should be modeled as an active energy-dependent process. The mechanism of action of such a protein folding machine might include direct manipulation of the peptide backbone. Methods: To show the feasibility of a protein folding machine, we conducted molecular dynamics simulations that were augmented by the application of mechanical force to rotate the C-terminal amino acid while simultaneously limiting the N-terminal amino acid movements. Results: Remarkably, the addition of this simple manipulation of peptide backbones to the standard molecular dynamics simulation indeed facilitated the formation of native structures in five diverse alpha-helical peptides. Steric clashes that arise in the peptides due to the forced directional rotation resulted in the behavior of the peptide backbone no longer resembling a freely jointed chain. Conclusions: These simulations show the feasibility of a protein folding machine operating under the conditions when the movements of the polypeptide backbone are restricted by applying external forces and constraints. Further investigation is needed to see whether such an effect may play a role during co-translational protein folding in vivo and how it can be utilized to facilitate folding of proteins in artificial environments.

Protein folding stability and binding interactions through the lens of evolution: a dynamical perspective

While the function of a protein depends heavily on its ability to fold into a correct 3D structure, billions of years of evolution have tailored proteins from highly stable objects to flexible molecules as they adapted to environmental changes. Nature maintains the fine balance of protein folding and stability while still evolving towards new function through generations of fine-tuning necessary interactions with other proteins and small molecules. Here we focus on recent computational and experimental studies that shed light onto how evolution molds protein folding and the functional landscape from a conformational dynamics' perspective. Particularly, we explore the importance of dynamic allostery throughout protein evolution and discuss how the protein anisotropic network can give rise to allosteric and epistatic interactions.

Piecewise All-Atom SMD Simulations Reveal Key Secondary Structures in Luciferase Unfolding Pathway

Although the folding of single-domain proteins is well characterized theoretically and experimentally, the folding of large multidomain proteins is less well known. Firefly luciferase, a 550 residue three-domain protein, has been commonly used as a substrate to study chaperone reactions and as a model system for the study of folding of long polypeptide chains, including related phenomena such as cotranslational folding. Despite being characterized by various experimental techniques, the atomic-level contributions of various secondary structures of luciferase to its fold's mechanical stability remain unknown. Here, we developed a piecewise approach for all-atom steered molecular dynamics simulations to examine specific secondary structures that resist mechanical unfolding while minimizing the amount of computational resources required by the large water box of standard all-atom steered molecular dynamics simulations. We validated the robustness of this approach with a small NI3C protein and used our approach to elucidate the specific secondary structures that provide the largest contributions to luciferase mechanostability. In doing so, we show that piecewise all-atom steered molecular dynamics simulations can provide novel atomic resolution details regarding mechanostability and can serve as a platform for novel mutagenesis studies as well as a point for comparison with high-resolution force spectroscopy experiments.

The Uniqueness of Tryptophan in Biology: Properties, Metabolism, Interactions and Localization in Proteins

Tryptophan (Trp) holds a unique place in biology for a multitude of reasons. It is the largest of all twenty amino acids in the translational toolbox. Its side chain is indole, which is aromatic with a binuclear ring structure, whereas those of Phe, Tyr, and His are single-ring aromatics. In part due to these elaborate structural features, the biosynthetic pathway of Trp is the most complex and the most energy-consuming among all amino acids. Essential in the animal diet, Trp is also the least abundant amino acid in the cell, and one of the rarest in the proteome. In most eukaryotes, Trp is the only amino acid besides Met, which is coded for by a single codon, namely UGG. Due to the large and hydrophobic π-electron surface area, its aromatic side chain interacts with multiple other side chains in the protein, befitting its strategic locations in the protein structure. Finally, several Trp derivatives, namely tryptophylquinone, oxitriptan, serotonin, melatonin, and tryptophol, have specialized functions. Overall, Trp is a scarce and precious amino acid in the cell, such that nature uses it parsimoniously, for multiple but selective functions. Here, the various aspects of the uniqueness of Trp are presented in molecular terms.

Validation of DBFOLD: An efficient algorithm for computing folding pathways of complex proteins

Atomistic simulations can provide valuable, experimentally-verifiable insights into protein folding mechanisms, but existing ab initio simulation methods are restricted to only the smallest proteins due to severe computational speed limits. The folding of larger proteins has been studied using native-centric potential functions, but such models omit the potentially crucial role of non-native interactions. Here, we present an algorithm, entitled DBFOLD, which can predict folding pathways for a wide range of proteins while accounting for the effects of non-native contacts. In addition, DBFOLD can predict the relative rates of different transitions within a protein’s folding pathway. To accomplish this, rather than directly simulating folding, our method combines equilibrium Monte-Carlo simulations, which deploy enhanced sampling, with unfolding simulations at high temperatures. We show that under certain conditions, trajectories from these two types of simulations can be jointly analyzed to compute unknown folding rates from detailed balance. This requires inferring free energies from the equilibrium simulations, and extrapolating transition rates from the unfolding simulations to lower, physiologically-reasonable temperatures at which the native state is marginally stable. As a proof of principle, we show that our method can accurately predict folding pathways and Monte-Carlo rates for the well-characterized Streptococcal protein G. We then show that our method significantly reduces the amount of computation time required to compute the folding pathways of large, misfolding-prone proteins that lie beyond the reach of existing direct simulation. Our algorithm, which is available online, can generate detailed atomistic models of protein folding mechanisms while shedding light on the role of non-native intermediates which may crucially affect organismal fitness and are frequently implicated in disease.

Many proteins must adopt a specific structure in order to function. Computational simulations have been used to shed light on the mechanisms of protein folding, but unfortunately, realistic simulations can typically only be run for small proteins, due to severe limits in computational speed. Here, we present a method to solve this problem, whereby instead of directly simulating folding from an unfolded state, we run simulations that allow for computation of equilibrium folding free energies, alongside high temperature simulations to compute unfolding rates. From these quantities, folding rates can be computed using detailed balance. Importantly, our method can account for the effects of nonnative contacts which transiently form during folding and must be broken prior to adoption of the native state. Such contacts, which are often excluded from simple models of folding, may crucially affect real protein folding pathways and are often observed in folding intermediates implicated in disease.

Metabolic Specialization and Codon Preference of Lignocellulolytic Genes in the White Rot Basidiomycete Ceriporiopsis subvermispora

Ceriporiopsis subvermispora is a white-rot fungus with a high specificity towards lignin mineralization when colonizing dead wood or lignocellulosic compounds. Its lignocellulose degrading system is formed by cellulose hydrolytic enzymes, manganese peroxidases, and laccases that catalyze the efficient depolymerization and mineralization of lignocellulose. To determine if this metabolic specialization has modified codon usage of the lignocellulolytic system, improving its adaptation to the fungal translational machine, we analyzed the adaptation to host codon usage (CAI), tRNA pool (tAI, and AAtAI), codon pair bias (CPB), and the number of effective codons (Nc). These indexes were correlated with gene expression of C. subvermispora, in the presence of glucose and Aspen wood. General gene expression was not correlated with the index values. However, in media containing Aspen wood, the induction of expression of lignocellulose-degrading genes, showed significantly (p < 0.001) higher values of CAI, AAtAI, CPB, tAI, and lower values of Nc than non-induced genes. Cellulose-binding proteins and manganese peroxidases presented the highest adaptation values. We also identified an expansion of genes encoding glycine and glutamic acid tRNAs. Our results suggest that the metabolic specialization to use wood as the sole carbon source has introduced a bias in the codon usage of genes involved in lignocellulose degradation. This bias reduces codon diversity and increases codon usage adaptation to the tRNA pool available in C. subvermispora. To our knowledge, this is the first study showing that codon usage is modified to improve the translation efficiency of a group of genes involved in a particular metabolic process.

CFTR trafficking mutations disrupt cotranslational protein folding by targeting biosynthetic intermediates

Protein misfolding causes a wide spectrum of human disease, and therapies that target misfolding are transforming the clinical care of cystic fibrosis. Despite this success, however, very little is known about how disease-causing mutations affect the de novo folding landscape. Here we show that inherited, disease-causing mutations located within the first nucleotide-binding domain (NBD1) of the cystic fibrosis transmembrane conductance regulator (CFTR) have distinct effects on nascent polypeptides. Two of these mutations (A455E and L558S) delay compaction of the nascent NBD1 during a critical window of synthesis. The observed folding defect is highly dependent on nascent chain length as well as its attachment to the ribosome. Moreover, restoration of the NBD1 cotranslational folding defect by second site suppressor mutations also partially restores folding of full-length CFTR. These findings demonstrate that nascent folding intermediates can play an important role in disease pathogenesis and thus provide potential targets for pharmacological correction.

Effect of Protein Structure on Evolution of Cotranslational Folding

Cotranslational folding depends on the folding speed and stability of the nascent protein. It remains difficult, however, to predict which proteins cotranslationally fold. Here, we simulate evolution of model proteins to investigate how native structure influences evolution of cotranslational folding. We developed a model that connects protein folding during and after translation to cellular fitness. Model proteins evolved improved folding speed and stability, with proteins adopting one of two strategies for folding quickly. Low contact order proteins evolve to fold cotranslationally. Such proteins adopt native conformations early on during the translation process, with each subsequently translated residue establishing additional native contacts. On the other hand, high contact order proteins tend not to be stable in their native conformations until the full chain is nearly extruded. We also simulated evolution of slowly translating codons, finding that slower translation speeds at certain positions enhances cotranslational folding. Finally, we investigated real protein structures using a previously published data set that identified evolutionarily conserved rare codons in Escherichia coli genes and associated such codons with cotranslational folding intermediates. We found that protein substructures preceding conserved rare codons tend to have lower contact orders, in line with our finding that lower contact order proteins are more likely to fold cotranslationally. Our work shows how evolutionary selection pressure can cause proteins with local contact topologies to evolve cotranslational folding.