In vivo translation rates can substantially delay the cotranslational folding of the Escherichia coli cytosolic proteome

Emerging frontiers in protein structure prediction following the AlphaFold revolution

Models of protein structures enable molecular understanding of biological processes. Current protein structure prediction tools lie at the interface of biology, chemistry and computer science. Millions of protein structure models have been generated in a very short space of time through a revolution in protein structure prediction driven by deep learning, led by AlphaFold. This has provided a wealth of new structural information. Interpreting these predictions is critical to determining where and when this information is useful. But proteins are not static nor do they act alone, and structures of proteins interacting with other proteins and other biomolecules are critical to a complete understanding of their biological function at the molecular level. This review focuses on the application of state-of-the-art protein structure prediction to these advanced applications. We also suggest a set of guidelines for reporting AlphaFold predictions.

Native Fold Delay and its implications for co-translational chaperone binding and protein aggregation

Because of vectorial protein translation, residues that interact in the native protein structure but are distantly separated in the primary sequence are unavailable simultaneously. Instead, there is a temporal delay during which the N-terminal interaction partner is unsatisfied and potentially vulnerable to non-native interactions. We introduce "Native Fold Delay" (NFD), a metric that integrates protein topology with translation kinetics to quantify such delays. We found that many proteins exhibit residues with NFDs in the range of tens of seconds. These residues, predominantly in well-structured, buried regions, often coincide with aggregation-prone regions. NFD correlates with co-translational engagement by the yeast Hsp70 chaperone Ssb, suggesting that native fold-delayed regions have a propensity to misfold. Supporting this, we show that proteins with long NFDs are more frequently co-translationally ubiquitinated and prone to aggregate upon Ssb deletion.

Three Stages of Nascent Protein Translocation Through the Ribosome Exit Tunnel

All proteins in living organisms are produced in ribosomes that facilitate the translation of genetic information into a sequence of amino acid residues. During translation, the ribosome undergoes initiation, elongation, termination, and recycling. In fact, peptide bonds are formed only during the elongation phase, which comprises periodic association of transfer RNAs and multiple auxiliary proteins with the ribosome and the addition of an amino acid to the nascent polypeptide one at a time. The protein spends a considerable amount of time attached to the ribosome. Here, we conceptually divide this portion of the protein lifetime into three stages. We define each stage on the basis of the position of the N-terminus of the nascent polypeptide within the ribosome exit tunnel and the context of the catalytic center. We argue that nascent polypeptides experience a variety of forces that determine how they translocate through the tunnel and interact with the tunnel walls. We review current knowledge about nascent polypeptide translocation and identify several white spots in our understanding of the birth of proteins.

Disruption of the nascent polypeptide-associated complex leads to reduced polyglutamine aggregation and toxicity

The nascent polypeptide-associate complex (NAC) is a heterodimeric chaperone complex that binds near the ribosome exit tunnel and is the first point of chaperone contact for newly synthesized proteins. Deletion of the NAC induces embryonic lethality in many multi-cellular organisms. Previous work has shown that the deletion of the NAC rescues cells from prion-induced cytotoxicity. This counterintuitive result led us to hypothesize that NAC disruption would improve viability in cells expressing human misfolding proteins. Here, we show that NAC disruption improves viability in cells expressing expanded polyglutamine and also leads to delayed and reduced aggregation of expanded polyglutamine and changes in polyglutamine aggregate morphology. Moreover, we show that NAC disruption leads to changes in de novo yeast prion induction. These results indicate that the NAC plays a critical role in aggregate organization as a potential therapeutic target in neurodegenerative disorders.

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.

Mapping Protein-Protein Interactions at Birth: Single-Particle Cryo-EM Analysis of a Ribosome-Nascent Globin Complex

Interactions between ribosome-bound nascent chains (RNCs) and ribosomal components are critical to elucidate the mechanism of cotranslational protein folding. Nascent protein-ribosome contacts within the ribosomal exit tunnel were previously assessed mostly in the presence of C-terminal stalling sequences, yet little is known about contacts taking place in the absence of these strongly interacting motifs. Further, there is nearly no information about ribosomal proteins (r-proteins) interacting with nascent chains within the outer surface of the ribosome. Here, we combine chemical cross-linking, single-particle cryo-EM, and fluorescence anisotropy decays to determine the structural features of ribosome-bound apomyoglobin (apoMb). Within the ribosomal exit tunnel core, interactions are similar to those identified in previous reports. However, once the RNC enters the tunnel vestibule, it becomes more dynamic and interacts with ribosomal RNA (rRNA) and the L23 r-protein. Remarkably, on the outer surface of the ribosome, RNCs interact mainly with a highly conserved nonpolar patch of the L23 r-protein. RNCs also comprise a compact and dynamic N-terminal region lacking contact with the ribosome. In all, apoMb traverses the ribosome and interacts with it via its C-terminal region, while N-terminal residues sample conformational space and form a compact subdomain before the entire nascent protein sequence departs from the ribosome.

Specific codons control cellular resources and fitness

As cellular engineering progresses from simply overexpressing proteins to imparting complex phenotypes through multigene expression, judicious appropriation of cellular resources is essential. Since codon use is degenerate and biased, codons may control cellular resources at a translational level. We investigate how partitioning transfer RNA (tRNA) resources by incorporating dissimilar codon usage can drastically alter interdependence of expression level and burden on the host. By isolating the effect of individual codons' use during translation elongation while eliminating confounding factors, we show that codon choice can trans-regulate fitness of the host and expression of other heterologous or native genes. We correlate specific codon usage patterns with host fitness and derive a coding scheme for multigene expression called the Codon Health Index (CHI, χ). This empirically derived coding scheme (χ) enables the design of multigene expression systems that avoid catastrophic cellular burden and is robust across several proteins and conditions.

Folding stabilities of ribosome-bound nascent polypeptides probed by mass spectrometry

The folding of most proteins occurs during the course of their translation while their tRNA-bound C termini are embedded in the ribosome. How the close proximity of nascent proteins to the ribosome influences their folding thermodynamics remains poorly understood. Here, we have developed a mass spectrometry-based approach for determining the stabilities of nascent polypeptide chains using methionine oxidation as a folding probe. This approach enables quantitative measurement subglobal folding stabilities of ribosome nascent chains within complex protein mixtures and extracts. To validate the methodology, we analyzed the folding thermodynamics of three model proteins (dihydrofolate reductase, chemotaxis protein Y, and DNA polymerase IV) in soluble and ribosome-bound states. The data indicate that the ribosome can significantly alter the stability of nascent polypeptides. Ribosome-induced stability modulations were highly variable among different folding domains and were dependent on localized charge distributions within nascent polypeptides. The results implicated electrostatic interactions between the ribosome surface and nascent polypeptides as the cause of ribosome-induced stability modulations. The study establishes a robust proteomic methodology for analyzing localized stabilities within ribosome-bound nascent polypeptides and sheds light on how the ribosome influences the thermodynamics of protein folding.

Is Posttranslational Folding More Efficient Than Refolding from a Denatured State: A Computational Study

The folding of proteins into their native conformation is a complex process that has been extensively studied over the past half-century. The ribosome, the molecular machine responsible for protein synthesis, is known to interact with nascent proteins, adding further complexity to the protein folding landscape. Consequently, it is unclear whether the folding pathways of proteins are conserved on and off the ribosome. The main question remains: to what extent does the ribosome help proteins fold? To address this question, we used coarse-grained molecular dynamics simulations to compare the mechanisms by which the proteins dihydrofolate reductase, type III chloramphenicol acetyltransferase, and d-alanine-d-alanine ligase B fold during and after vectorial synthesis on the ribosome to folding from the full-length unfolded state in bulk solution. Our results reveal that the influence of the ribosome on protein folding mechanisms varies depending on the size and complexity of the protein. Specifically, for a small protein with a simple fold, the ribosome facilitates efficient folding by helping the nascent protein avoid misfolded conformations. However, for larger and more complex proteins, the ribosome does not promote folding and may contribute to the formation of intermediate misfolded states cotranslationally. These misfolded states persist posttranslationally and do not convert to the native state during the 6 μs runtime of our coarse-grain simulations. Overall, our study highlights the complex interplay between the ribosome and protein folding and provides insight into the mechanisms of protein folding on and off the ribosome.

Hold the fold: how delayed folding aids protein secretion

In bacteria, N-terminal signal peptides mark proteins for transport across the plasma membrane. A recent study by Smets et al (2022) followed the folding of a pair of structural twins to shed light on how evolution has optimised the secretory process.

Ribosome profiling enhances understanding of mycobacterial translation

A recent addition to the -omics toolkit, ribosome profiling, enables researchers to gain insight into the process and regulation of translation by mapping fragments of mRNA protected from nuclease digestion by ribosome binding. In this review, we discuss how ribosome profiling applied to mycobacteria has led to discoveries about translational regulation. Using case studies, we show that the traditional view of “canonical” translation mechanisms needs expanding to encompass features of mycobacterial translation that are more widespread than previously recognized. We also discuss the limitations of the method and potential future developments that could yield further insight into the fundamental biology of this important human pathogen.

Large protein complex interfaces have evolved to promote cotranslational assembly

Assembly pathways of protein complexes should be precise and efficient to minimise misfolding and unwanted interactions with other proteins in the cell. One way to achieve this efficiency is by seeding assembly pathways during translation via the cotranslational assembly of subunits. While recent evidence suggests that such cotranslational assembly is widespread, little is known about the properties of protein complexes associated with the phenomenon. Here, using a combination of proteome-specific protein complex structures and publicly available ribosome profiling data, we show that cotranslational assembly is particularly common between subunits that form large intermolecular interfaces. To test whether large interfaces have evolved to promote cotranslational assembly, as opposed to cotranslational assembly being a non-adaptive consequence of large interfaces, we compared the sizes of first and last translated interfaces of heteromeric subunits in bacterial, yeast, and human complexes. When considering all together, we observe the N-terminal interface to be larger than the C-terminal interface 54% of the time, increasing to 64% when we exclude subunits with only small interfaces, which are unlikely to cotranslationally assemble. This strongly suggests that large interfaces have evolved as a means to maximise the chance of successful cotranslational subunit binding.

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.

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.

Engineering of the Recombinant Expression and PEGylation Efficiency of the Therapeutic Enzyme Human Thymidine Phosphorylase

Human thymidine phosphorylase (HsTP) is an enzyme with important implications in the field of rare metabolic diseases. Defective mutations of HsTP lead to mitochondrial neurogastrointestinal encephalomyopathy (MNGIE), a disease with a high unmet medical need that is associated with severe neurological and gastrointestinal complications. Current efforts focus on the development of an enzyme replacement therapy (ERT) using the Escherichia coli ortholog (EcTP). However, bacterial enzymes are counter-indicated for human therapeutic applications because they are recognized as foreign by the human immune system, thereby eliciting adverse immune responses and raising significant safety and efficacy risks. Thus, it is critical to utilize the HsTP enzyme as starting scaffold for pre-clinical drug development, thus de-risking the safety concerns associated with the use of bacterial enzymes. However, HsTP expresses very poorly in E. coli, whereas its PEGylation, a crucial chemical modification for achieving long serum persistence of therapeutic enzymes, is highly inefficient and negatively affects its catalytic activity. Here we focused on the engineering of the recombinant expression profile of HsTP in E. coli cells, as well as on the optimization of its PEGylation efficiency aiming at the development of an alternative therapeutic approach for MNGIE. We show that phylogenetic and structural analysis of proteins can provide important insights for the rational design of N’-terminus-truncation constructs which exhibit significantly improved recombinant expression levels. In addition, we developed and implemented a criteria-driven rational surface engineering strategy for the substitution of arginine-to-lysine and lysine-to-arginine residues to achieve more efficient, homogeneous and reproducible PEGylation without negatively affecting the enzymatic catalytic activity upon PEGylation. Collectively, our proposed strategies provide an effective way to optimize enzyme PEGylation and E. coli recombinant expression and are likely applicable for other proteins and enzymes.

Translational autoregulation of RF2 protein in E. coli through programmed frameshifting

Various feedback mechanisms regulate the expression of different genes to ensure the required protein levels inside a cell. In this paper, we develop a kinetic model for one such mechanism that autoregulates RF2 protein synthesis in E. coli through programmed frameshifting. The model finds that the programmed frameshifting autoregulates RF2 protein synthesis by two independent mechanisms. First, it increases the rate of RF2 synthesis from each mRNA transcript at low RF2 concentration. Second, programmed frameshifting can dramatically increase the lifetime of RF2 transcripts when RF2 protein levels are lower than a threshold. This sharp increase in mRNA lifetime is caused by a first-order phase transition from a low to a high ribosome density on an RF2 transcript. The high ribosome density prevents the transcript's degradation by shielding it from nucleases, which increases its average lifetime and hence RF2 protein levels. Our study identifies this quality control mechanism that regulates the cellular protein levels by breaking the hierarchy of processes involved in gene expression.

How Does the Ribosome Fold the Proteome?

Folding of polypeptides begins during their synthesis on ribosomes. This process has evolved as a means for the cell to maintain proteostasis, by mitigating the risk of protein misfolding and aggregation. The capacity to now depict this cellular feat at increasingly higher resolution is providing insight into the mechanistic determinants that promote successful folding. Emerging from these studies is the intimate interplay between protein translation and folding, and within this the ribosome particle is the key player. Its unique structural properties provide a specialized scaffold against which nascent polypeptides can begin to form structure in a highly coordinated, co-translational manner. Here, we examine how, as a macromolecular machine, the ribosome modulates the intrinsic dynamic properties of emerging nascent polypeptide chains and guides them toward their biologically active structures.

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.

Riboexp: an interpretable reinforcement learning framework for ribosome density modeling

Translation elongation is a crucial phase during protein biosynthesis. In this study, we develop a novel deep reinforcement learning-based framework, named Riboexp, to model the determinants of the uneven distribution of ribosomes on mRNA transcripts during translation elongation. In particular, our model employs a policy network to perform a context-dependent feature selection in the setting of ribosome density prediction. Our extensive tests demonstrated that Riboexp can significantly outperform the state-of-the-art methods in predicting ribosome density by up to 5.9% in terms of per-gene Pearson correlation coefficient on the datasets from three species. In addition, Riboexp can indicate more informative sequence features for the prediction task than other commonly used attribution methods in deep learning. In-depth analyses also revealed the meaningful biological insights generated by the Riboexp framework. Moreover, the application of Riboexp in codon optimization resulted in an increase of protein production by around 31% over the previous state-of-the-art method that models ribosome density. These results have established Riboexp as a powerful and useful computational tool in the studies of translation dynamics and protein synthesis. Availability: The data and code of this study are available on GitHub: https://github.com/Liuxg16/Riboexp. Contact:  zengjy321@tsinghua.edu.cn; songsen@tsinghua.edu.cn.

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.

Interplay between Conformational Strain and Intramolecular Interaction in Protein Structures: Which of Them Is Evolutionarily Conserved?

By computing strain energies of peptide fragments within protein structures and their intramolecular interaction energies, we attempt to reveal general biophysical trends behind the secondary structure formation in the context of protein evolution. Our "protein basis set" consisted of 1143 representatives of different folds obtained from curated SCOPe database, and for each member of the set, the strain and intramolecular energy was calculated on the "rolling tripeptide" basis, employing the DFT-D3/COSMO-RS method for the former and the QM-calibrated force field method (MM) for the latter. The calculated data, strain and interactions, were correlated with the conservation of amino acid residues in secondary structure elements and also with the level of the residue burial within the protein three-dimensional structure. It allowed us to formulate several observations concerning fundamental differences between two main secondary structure motifs: α-helices and β-strands. We have shown that a strong interaction is one of the determining characteristics of the β-sheet formation, at least at the level of tripeptides (and likely penta- or heptapeptides, too), and that the β-strand is a prevailing secondary structure in the strongly-interacting regions of the protein folds conserved by evolution. On the other hand, low strain was neither proven to be an important physicochemical property conserved by evolution nor does it correlate with the propensity for the α-helix and β-strand. Finally, it has been demonstrated that the strong interaction has a certain level of connection with residue burial; however, we demonstrate that these two characteristics should be rather regarded as two complementary factors. These findings represent an important contribution to understanding protein folding from first principles, which is a complementary approach to ongoing efforts to solve the protein folding problem by knowledge-based approaches and machine-learning.

Dynamics of Co-translational Membrane Protein Integration and Translocation via the Sec Translocon

An important aspect of cellular function is the correct targeting and delivery of newly synthesized proteins. Central to this task is the machinery of the Sec translocon, a transmembrane channel that is involved in both the translocation of nascent proteins across cell membranes and the integration of proteins into the membrane. Considerable experimental and computational effort has focused on the Sec translocon and its role in nascent protein biosynthesis, including the correct folding and expression of integral membrane proteins. However, the use of molecular simulation methods to explore Sec-facilitated protein biosynthesis is hindered by the large system sizes and long (i.e., minute) time scales involved. In this work, we describe the development and application of a coarse-grained simulation approach that addresses these challenges and allows for direct comparison with both in vivo and in vitro experiments. The method reproduces a wide range of experimental observations, providing new insights into the underlying molecular mechanisms, predictions for new experiments, and a strategy for the rational enhancement of membrane protein expression levels.

Electrostatic Interactions Govern Extreme Nascent Protein Ejection Times from Ribosomes and Can Delay Ribosome Recycling

The ejection of nascent proteins out of the ribosome exit tunnel, after their covalent bond to transfer-RNA has been broken, has not been experimentally studied due to challenges in sample preparation. Here, we investigate this process using a combination of multiscale modeling, ribosome profiling, and gene ontology analyses. Simulating the ejection of a representative set of 122 E. coli proteins we find a greater than 1000-fold variation in ejection times. Nascent proteins enriched in negatively charged residues near their C-terminus eject the fastest, while nascent chains enriched in positively charged residues tend to eject much more slowly. More work is required to pull slowly ejecting proteins out of the exit tunnel than quickly ejecting proteins, according to all-atom simulations. An energetic decomposition reveals, for slowly ejecting proteins, that this is due to the strong attractive electrostatic interactions between the nascent chain and the negatively charged ribosomal-RNA lining the exit tunnel, and for quickly ejecting proteins, it is due to their repulsive electrostatic interactions with the exit tunnel. Ribosome profiling data from E. coli reveals that the presence of slowly ejecting sequences correlates with ribosomes spending more time at stop codons, indicating that the ejection process might delay ribosome recycling. Proteins that have the highest positive charge density at their C-terminus are overwhelmingly ribosomal proteins, suggesting the possibility that this sequence feature may aid in the cotranslational assembly of ribosomes by delaying the release of nascent ribosomal proteins into the cytosol. Thus, nascent chain ejection times from the ribosome can vary greatly between proteins due to differential electrostatic interactions, can influence ribosome recycling, and could be particularly relevant to the synthesis and cotranslational behavior of some proteins.

Codon Selection Affects Recruitment of Ribosome-Associating Factors during Translation

An intriguing aspect of protein synthesis is how cotranslational events are managed inside the cell. In this study, we developed an in vivo bimolecular fluorescence complementation assay coupled to SecM stalling (BiFC-SecM) to study how codon usage influences the interactions of ribosome-associating factors that occur cotranslationally. We profiled ribosomal associations of a number of proteins, and observed differential association of chaperone proteins TF, DnaK, GroEL, and translocation factor Ffh as a result of introducing synonymous codon substitutions that change the affinity of the translating sequence to the ribosomal anti-Shine-Dalgarno (aSD) sequence. The use of pausing sequences within proteins regulates their transit within the translating ribosome. Our results indicate that the dynamics between cellular factors and the new polypeptide chain are affected by how codon composition is designed. Furthermore, associating factors may play a role in processes including protein quality control (folding and degradation) and cellular respiration.

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

Many large proteins suffer from slow or inefficient folding in vitro. It has long been known that this problem can be alleviated in vivo if proteins start folding cotranslationally. However, the molecular mechanisms underlying this improvement have not been well established. To address this question, we use an all-atom simulation-based algorithm to compute the folding properties of various large protein domains as a function of nascent chain length. We find that for certain proteins, there exists a narrow window of lengths that confers both thermodynamic stability and fast folding kinetics. Beyond these lengths, folding is drastically slowed by nonnative interactions involving C-terminal residues. Thus, cotranslational folding is predicted to be beneficial because it allows proteins to take advantage of this optimal window of lengths and thus avoid kinetic traps. Interestingly, many of these proteins' sequences contain conserved rare codons that may slow down synthesis at this optimal window, suggesting that synthesis rates may be evolutionarily tuned to optimize folding. Using kinetic modeling, we show that under certain conditions, such a slowdown indeed improves cotranslational folding efficiency by giving these nascent chains more time to fold. In contrast, other proteins are predicted not to benefit from cotranslational folding due to a lack of significant nonnative interactions, and indeed these proteins' sequences lack conserved C-terminal rare codons. Together, these results shed light on the factors that promote proper protein folding in the cell and how biomolecular self-assembly may be optimized evolutionarily.

Cotranslational Folding of Proteins on the Ribosome

Many proteins in the cell fold cotranslationally within the restricted space of the polypeptide exit tunnel or at the surface of the ribosome. A growing body of evidence suggests that the ribosome can alter the folding trajectory in many different ways. In this review, we summarize the recent examples of how translation affects folding of single-domain, multiple-domain and oligomeric proteins. The vectorial nature of translation, the spatial constraints of the exit tunnel, and the electrostatic properties of the ribosome-nascent peptide complex define the onset of early folding events. The ribosome can facilitate protein compaction, induce the formation of intermediates that are not observed in solution, or delay the onset of folding. Examples of single-domain proteins suggest that early compaction events can define the folding pathway for some types of domain structures. Folding of multi-domain proteins proceeds in a domain-wise fashion, with each domain having its role in stabilizing or destabilizing neighboring domains. Finally, the assembly of protein complexes can also begin cotranslationally. In all these cases, the ribosome helps the nascent protein to attain a native fold and avoid the kinetic traps of misfolding.

Mechanochemistry in Translation

As the influence of translation rates on protein folding and function has come to light, the mechanisms by which translation speed is modulated have become an important issue. One mechanism entails the generation of force by the nascent protein. Cotranslational processes, such as nascent protein folding, the emergence of unfolded nascent chain segments from the ribosome's exit tunnel, and insertion of the nascent chain into or translocation of the nascent chain through membranes, can generate forces that are transmitted back to the peptidyl transferase center and affect translation rates. In this Perspective, we examine the processes that generate these forces, the mechanisms of transmission along the ribosomal exit tunnel to the peptidyl transferase center, and the effects of force on the ribosome's catalytic cycle. We also discuss the physical models that have been developed to predict and explain force generation for individual processes and speculate about other processes that may generate forces that have yet to be tested.

Timing and specificity of cotranslational nascent protein modification in bacteria

The nascent polypeptide exit site of the ribosome is a crowded environment where multiple ribosome-associated protein biogenesis factors (RPBs) compete for the nascent polypeptide to influence their localization, folding, or quality control. Here we address how N-terminal methionine excision (NME), a ubiquitous process crucial for the maturation of over 50% of the bacterial proteome, occurs in a timely and selective manner in this crowded environment. In bacteria, NME is mediated by 2 essential enzymes, peptide deformylase (PDF) and methionine aminopeptidase (MAP). We show that the reaction of MAP on ribosome-bound nascent chains approaches diffusion-limited rates, allowing immediate methionine excision of optimal substrates after deformylation. Specificity is achieved by kinetic competition of NME with translation elongation and by regulation from other RPBs, which selectively narrow the processing time window for suboptimal substrates. A mathematical model derived from the data accurately predicts cotranslational NME efficiency in the cytosol. Our results demonstrate how a fundamental enzymatic activity is reshaped by its associated macromolecular environment to optimize both efficiency and selectivity, and provides a platform to study other cotranslational protein biogenesis pathways.

Ligand Binding Site Structure Shapes Folding, Assembly and Degradation of Homomeric Protein Complexes

Ligand binding site structure has profound consequences for the evolution of function of protein complexes, particularly in homomers—complexes comprising multiple copies of the same protein. Previously, we have shown that homomers with multichain binding sites (MBSs) are characterized by more conserved binding sites and quaternary structure, and qualitatively different allosteric pathways than homomers with single-chain binding sites (SBSs) or monomers. Here, using computational methods, we show that the folds of single-domain MBS and SBS homomers are different, and SBS homomers are likely to be folded cotranslationally, while MBS homomers are more likely to form post-translationally and rely on more advanced folding-assistance and quality control mechanisms, which include chaperonins. In addition, our findings demonstrate that MBS homomers are qualitatively different from monomers, while SBS homomers are much less distinct, supporting the hypothesis that the evolution of quaternary structure in SBS homomers is significantly influenced by stochastic processes.

Nature and Regulation of Protein Folding on the Ribosome

Co-translational protein folding is an essential process by which cells ensure the safe and efficient production and assembly of new proteins in their functional native states following biosynthesis on the ribosome. In this review, we describe recent progress in probing the changes during protein synthesis of the free energy landscapes that underlie co-translational folding and discuss the critical coupling between these landscapes and the rate of translation that ultimately determines the success or otherwise of the folding process. Recent developments have revealed a variety of mechanisms by which both folding and translation can be modulated or regulated, and we discuss how these effects are utilised by the cell to optimise the outcome of protein biosynthesis.

Non-equilibrium dynamics of a nascent polypeptide during translation suppress its misfolding

Protein folding can begin co-translationally. Due to the difference in timescale between folding and synthesis, co-translational folding is thought to occur at equilibrium for fast-folding domains. In this scenario, the folding kinetics of stalled ribosome-bound nascent chains should match the folding of nascent chains in real time. To test if this assumption is true, we compare the folding of a ribosome-bound, multi-domain calcium-binding protein stalled at different points in translation with the nascent chain as is it being synthesized in real-time, via optical tweezers. On stalled ribosomes, a misfolded state forms rapidly (1.5 s). However, during translation, this state is only attained after a long delay (63 s), indicating that, unexpectedly, the growing polypeptide is not equilibrated with its ensemble of accessible conformations. Slow equilibration on the ribosome can delay premature folding until adequate sequence is available and/or allow time for chaperone binding, thus promoting productive folding.

Co-translational protein folding is thought to occur at equilibrium for fast-folding domains. Here authors use optical tweezers to show that the folding kinetics of stalled ribosome-bound nascent chains do not match the folding of nascent chains in real time.

Mechanisms of Cotranslational Maturation of Newly Synthesized Proteins

The timely production of functional proteins is of critical importance for the biological activity of cells. To reach the functional state, newly synthesized polypeptides have to become enzymatically processed, folded, and assembled into oligomeric complexes and, for noncytosolic proteins, translocated across membranes. Key activities of these processes occur cotranslationally, assisted by a network of machineries that transiently engage nascent polypeptides at distinct phases of translation. The sequence of events is tuned by intrinsic features of the nascent polypeptides and timely association of factors with the translating ribosome. Considering the dynamics of translation, the heterogeneity of cellular proteins, and the diversity of interaction partners, it is a major cellular achievement that these processes are temporally and spatially so precisely coordinated, minimizing the generation of damaged proteins. This review summarizes the current progress we have made toward a comprehensive understanding of the cotranslational interactions of nascent chains, which pave the way to their functional state.

Within-Gene Shine-Dalgarno Sequences Are Not Selected for Function

The Shine-Dalgarno (SD) sequence motif facilitates translation initiation and is frequently found upstream of bacterial start codons. However, thousands of instances of this motif occur throughout the middle of protein coding genes in a typical bacterial genome. Here, we use comparative evolutionary analysis to test whether SD sequences located within genes are functionally constrained. We measure the conservation of SD sequences across Enterobacteriales, and find that they are significantly less conserved than expected. Further, the strongest SD sequences are the least conserved whereas we find evidence of conservation for the weakest possible SD sequences given amino acid constraints. Our findings indicate that most SD sequences within genes are likely to be deleterious and removed via selection. To illustrate the origin of these deleterious costs, we show that ATG start codons are significantly depleted downstream of SD sequences within genes, highlighting the constraint that these sequences impose on the surrounding nucleotides to minimize the potential for erroneous translation initiation.

Novel insights into gene expression regulation during meiosis revealed by translation elongation dynamics

The ability to dynamically control mRNA translation has a great impact on many intracellular processes. Whereas it is believed that translational control in eukaryotes occurs mainly at initiation, the condition-specific changes at the elongation level and their potential regulatory role remain unclear. Using computational approaches applied to ribosome profiling data, we show that elongation rate is dynamic and can change considerably during the yeast meiosis to facilitate the selective translation of stage-specific transcripts. We observed unique elongation changes during meiosis II, including a global inhibition of translation elongation at the onset of anaphase II accompanied by a sharp shift toward increased elongation for genes required at this meiotic stage. We also show that ribosomal proteins counteract the global decreased elongation by maintaining high initiation rates. Our findings provide new insights into gene expression regulation during meiosis and demonstrate that codon usage evolved, among others, to optimize timely translation.

Domain topology, stability, and translation speed determine mechanical force generation on the ribosome

The concomitant folding of a nascent protein domain with its synthesis can generate mechanical forces that act on the ribosome and alter translation speed. Such changes in speed can affect the structure and function of the newly synthesized protein as well as cellular phenotype. The domain properties that govern force generation have yet to be identified and understood, and the influence of translation speed is unknown because all reported measurements have been carried out on arrested ribosomes. Here, using coarse-grained molecular simulations and statistical mechanical modeling of protein synthesis, we demonstrate that force generation is determined by a domain's stability and topology, as well as translation speed. The statistical mechanical models we create predict how force profiles depend on these properties. These results indicate that force measurements on arrested ribosomes will not always accurately reflect what happens in a cell, especially for slow-folding domains, and suggest the possibility that certain domain properties may be enriched or depleted across the structural proteome of organisms through evolutionary selection pressures to modulate protein synthesis speed and posttranslational protein behavior.

Structural Origins of FRET-Observed Nascent Chain Compaction on the Ribosome

A fluorescence signal arising from a Förster resonance energy transfer process was used to monitor conformational changes of a domain within the E. coli protein HemK during its synthesis by the ribosome. An increase in fluorescence was observed to begin 10 s after translation was initiated, indicating the domain became more compact in size. Since fluorescence only reports a single value at each time point it contains very little information about the structural ensemble that gives rise to it. Here, we supplement this experimental information with coarse-grained simulations that describe protein conformations and transitions at a spatial resolution of 3.8 Å. We use these simulations to test three hypotheses that might explain the cause of domain compaction: (1) that poor solvent quality conditions drive the unfolded state to compact, (2) that a change in the dimension of the space the domain occupies upon moving outside the exit tunnel causes compaction, or (3) that domain folding causes compaction. We find that domain folding and dimensional collapse are both consistent with the experimental data, while poor-solvent collapse is inconsistent. We identify alternative dye labeling positions on HemK that upon fluorescence can differentiate between the domain folding and dimensional collapse mechanisms. Partial folding of domains has been observed in C-terminally truncated forms of proteins. Therefore, it is likely that the experimentally observed compact state is a partially folded intermediate consisting, according to our simulations, of the first three helices of the HemK N-terminal domain adopting a native, tertiary configuration. With these simulations we also identify the possible cotranslational folding pathways of HemK.

Investigating the Effect of Chain Connectivity on the Folding of a Beta-Sheet Protein On and Off the Ribosome

Determining the relationship between protein folding pathways on and off the ribosome remains an important area of investigation in biology. Studies on isolated domains have shown that alteration of the separation of residues in a polypeptide chain, while maintaining their spatial contacts, may affect protein stability and folding pathway. Due to the vectorial emergence of the polypeptide chain from the ribosome, chain connectivity may have an important influence upon cotranslational folding. Using MATH, an all β-sandwich domain, we investigate whether the connectivity of residues and secondary structure elements is a key determinant of when cotranslational folding can occur on the ribosome. From Φ-value analysis, we show that the most structured region of the transition state for folding in MATH includes the N and C terminal strands, which are located adjacent to each other in the structure. However, arrest peptide force-profile assays show that wild-type MATH is able to fold cotranslationally, while some C-terminal residues remain sequestered in the ribosome, even when destabilized by 2–3 kcal mol⁻¹. We show that, while this pattern of Φ-values is retained in two circular permutants in our studies of the isolated domains, one of these permutants can fold only when fully emerged from the ribosome. We propose that in the case of MATH, onset of cotranslational folding is determined by the ability to form a sufficiently stable folding nucleus involving both β-sheets, rather than by the location of the terminal strands in the ribosome tunnel.

•

Adjacent N and C terminal strands are most structured region in the transition state.

•

Two circular permutants retain the same folding pathway as wild-type MATH.

•

On the ribosome, early emergence of terminal strands does not promote earlier folding.

•

Formation of both β-sheets is energetically critical for folding on the ribosome.

•

Folding pathway minimizes formation of partly structured states prone to mis-folding.

Protein Sequences Recapitulate Genetic Code Evolution

Several hypotheses predict ranks of amino acid assignments to genetic code's codons. Analyses here show that average positions of amino acid species in proteins correspond to assignment ranks, in particular as predicted by Juke's neutral mutation hypothesis for codon assignments. In all tested protein groups, including co- and post-translationally folding proteins, 'recent' amino acids are on average closer to gene 5' extremities than 'ancient' ones. Analyses of pairwise residue contact energies matrices suggest that early amino acids stereochemically selected late ones that stablilize residue interactions within protein cores, presumably producing 5'-late-to-3'-early amino acid protein sequence gradients. The gradient might reduce protein misfolding, also after mutations, extending principles of neutral mutations to protein folding. Presumably, in self-perpetuating and self-correcting systems like the genetic code, initial conditions produce similarities between evolution of the process (the genetic code) and 'ontogeny' of resulting structures (here proteins), producing apparent teleonomy between process and product.

Does mRNA structure contain genetic information for regulating co-translational protein folding?

Currently many facets of genetic information are illdefined. In particular, how protein folding is genetically regulated has been a long-standing issue for genetics and protein biology. And a generic mechanistic model with supports of genomic data is still lacking. Recent technological advances have enabled much needed genome-wide experiments. While putting the effect of codon optimality on debate, these studies have supplied mounting evidence suggesting a role of mRNA structure in the regulation of protein folding by modulating translational elongation rate. In conjunctions with previous theories, this mechanistic model of protein folding guided by mRNA structure shall expand our understandings of genetic information and offer new insights into various biomedical puzzles.

Non-equilibrium coupling of protein structure and function to translation-elongation kinetics

Protein folding research has been dominated by the assumption that thermodynamics determines protein structure and function. And that when the folding process is compromised in vivo the proteostasis machinery-chaperones, deaggregases, the proteasome-work to restore proteins to their soluble, functional form or degrade them to maintain the cellular pool of proteins in a quasi-equilibrium state. During the past decade, however, more and more proteins have been identified for which altering only their speed of synthesis alters their structure and function, the efficiency of the down-stream processes they take part in, and cellular phenotype. Indeed, evidence has emerged that evolutionary selection pressures have encoded translation-rate information into mRNA molecules to coordinate diverse co-translational processes. Thus, non-equilibrium physics can play a fundamental role in influencing nascent protein behavior, mRNA sequence evolution, and disease. Here, we discuss how our understanding of this phenomenon is being advanced by the application of theoretical tools from the physical sciences.

Evidence of evolutionary selection for cotranslational folding

Recent experiments and simulations have demonstrated that proteins can fold on the ribosome. However, the extent and generality of fitness effects resulting from cotranslational folding remain open questions. Here we report a genome-wide analysis that uncovers evidence of evolutionary selection for cotranslational folding. We describe a robust statistical approach to identify loci within genes that are both significantly enriched in slowly translated codons and evolutionarily conserved. Surprisingly, we find that domain boundaries can explain only a small fraction of these conserved loci. Instead, we propose that regions enriched in slowly translated codons are associated with cotranslational folding intermediates, which may be smaller than a single domain. We show that the intermediates predicted by a native-centric model of cotranslational folding account for the majority of these loci across more than 500 Escherichia coli proteins. By making a direct connection to protein folding, this analysis provides strong evidence that many synonymous substitutions have been selected to optimize translation rates at specific locations within genes. More generally, our results indicate that kinetics, and not just thermodynamics, can significantly alter the efficiency of self-assembly in a biological context.

Protein folding and tRNA biology

Polypeptides can fold into tertiary structures while they are synthesized by the ribosome. In addition to the amino acid sequence, protein folding is determined by several factors within the cell. Among others, the folding pathway of a nascent polypeptide can be affected by transient interactions with other proteins, ligands, or the ribosome, as well as by the translocation through membrane pores. Particularly, the translation machinery and the population of tRNA under different physiological or adaptive responses can dramatically affect protein folding. This review summarizes the scientific evidence describing the role of translation kinetics and tRNA populations on protein folding and addresses current efforts to better understand tRNA biology. It is organized into three main parts, which are focused on: (i) protein folding in the cellular context; (ii) tRNA biology and the complexity of the tRNA population; and (iii) available methods and technical challenges in the characterization of tRNA pools. In this manner, this work illustrates the ways by which functional properties of proteins may be modulated by cellular tRNA populations.

Miscoding-induced stalling of substrate translocation on the bacterial ribosome

Directional transit of the ribosome along the messenger RNA (mRNA) template is a key determinant of the rate and processivity of protein synthesis. Imaging of the multistep translocation mechanism using single-molecule FRET has led to the hypothesis that substrate movements relative to the ribosome resolve through relatively long-lived late intermediates wherein peptidyl-tRNA enters the P site of the small ribosomal subunit via reversible, swivel-like motions of the small subunit head domain within the elongation factor G (GDP)-bound ribosome complex. Consistent with translocation being rate-limited by recognition and productive engagement of peptidyl-tRNA within the P site, we now show that base-pairing mismatches between the peptidyl-tRNA anticodon and the mRNA codon dramatically delay this rate-limiting, intramolecular process. This unexpected relationship between aminoacyl-tRNA decoding and translocation suggests that miscoding antibiotics may impact protein synthesis by impairing the recognition of peptidyl-tRNA in the small subunit P site during EF-G-catalyzed translocation. Strikingly, we show that elongation factor P (EF-P), traditionally known to alleviate ribosome stalling at polyproline motifs, can efficiently rescue translocation defects arising from miscoding. These findings help reveal the nature and origin of the rate-limiting steps in substrate translocation on the bacterial ribosome and indicate that EF-P can aid in resuming translation elongation stalled by miscoding errors.

Protein Folding Mediated by Trigger Factor and Hsp70: New Insights from Single-Molecule Approaches

Chaperones assist in protein folding, but what this common phrase means in concrete terms has remained surprisingly poorly understood. We can readily measure chaperone binding to unfolded proteins, but how they bind and affect proteins along folding trajectories has remained obscure. Here we review recent efforts by our labs and others that are beginning to pry into this issue, with a focus on the chaperones trigger factor and Hsp70. Single-molecule methods are central, as they allow the stepwise process of folding to be followed directly. First results have already revealed contrasts with long-standing paradigms: rather than acting only "early" by stabilizing unfolded chain segments, these chaperones can bind and stabilize partially folded structures as they grow to their native state. The findings suggest a fundamental redefinition of the protein folding problem and a more extensive functional repertoire of chaperones than previously assumed.

Genetic Code Optimization for Cotranslational Protein Folding: Codon Directional Asymmetry Correlates with Antiparallel Betasheets, tRNA Synthetase Classes

A new codon property, codon directional asymmetry in nucleotide content (CDA), reveals a biologically meaningful genetic code dimension: palindromic codons (first and last nucleotides identical, codon structure XZX) are symmetric (CDA = 0), codons with structures ZXX/XXZ are 5'/3' asymmetric (CDA = - 1/1; CDA = - 0.5/0.5 if Z and X are both purines or both pyrimidines, assigning negative/positive (-/+) signs is an arbitrary convention). Negative/positive CDAs associate with (a) Fujimoto's tetrahedral codon stereo-table; (b) tRNA synthetase class I/II (aminoacylate the 2'/3' hydroxyl group of the tRNA's last ribose, respectively); and (c) high/low antiparallel (not parallel) betasheet conformation parameters. Preliminary results suggest CDA-whole organism associations (body temperature, developmental stability, lifespan). Presumably, CDA impacts spatial kinetics of codon-anticodon interactions, affecting cotranslational protein folding. Some synonymous codons have opposite CDA sign (alanine, leucine, serine, and valine), putatively explaining how synonymous mutations sometimes affect protein function. Correlations between CDA and tRNA synthetase classes are weaker than between CDA and antiparallel betasheet conformation parameters. This effect is stronger for mitochondrial genetic codes, and potentially drives mitochondrial codon-amino acid reassignments. CDA reveals information ruling nucleotide-protein relations embedded in reversed (not reverse-complement) sequences (5'-ZXX-3'/5'-XXZ-3').

Fast Protein Translation Can Promote Co- and Posttranslational Folding of Misfolding-Prone Proteins

Chemical kinetic modeling has previously been used to predict that fast-translating codons can enhance cotranslational protein folding by helping to avoid misfolded intermediates. Consistent with this prediction, protein aggregation in yeast and worms was observed to increase when translation was globally slowed down, possibly due to increased cotranslational misfolding. Observation of similar behavior in molecular simulations would confirm predictions from the simpler chemical kinetic model and provide a molecular perspective on cotranslational folding, misfolding, and the impact of translation speed on these processes. All-atom simulations cannot reach the timescales relevant to protein synthesis, and most conventional structure-based coarse-grained models do not allow for nonnative structure formation. Here, we introduce a protocol to incorporate misfolding using the functional forms of publicly available force fields. With this model we create two artificial proteins that are capable of undergoing structural transitions between a native and a misfolded conformation and simulate their synthesis by the ribosome. Consistent with the chemical kinetic predictions, we find that rapid synthesis of misfolding-prone nascent-chain segments increases the fraction of folded proteins by kinetically partitioning more molecules through on-pathway intermediates, decreasing the likelihood of sampling misfolded conformations. Novel to this study, to our knowledge, we observe that differences in protein dynamics, arising from different translation-elongation schedules, can persist long after the nascent protein has been released from the ribosome, and that a sufficient level of energetic frustration is needed for fast-translating codons to be beneficial for folding. These results provide further evidence that fast-translating codons can be as biologically important as pause sites in coordinating cotranslational folding.

Increasing Protein Production Rates Can Decrease the Rate at Which Functional Protein Is Produced and Their Steady-State Levels

The rate at which soluble, functional protein is produced by the ribosome has recently been found to vary in complex and unexplained ways as various translation-associated rates are altered through synonymous codon substitutions. To understand this phenomenon, here, we combine a well-established ribosome-traffic model with a master-equation model of cotranslational domain folding to explore the scenarios that are possible for the protein production rate, J, and the functional-nascent protein production rate, F, as the rates of various translation processes are altered for five different E. coli proteins. We find that while J monotonically increases as the rates of translation-initiation, -elongation, and -termination increase, F can either increase or decrease. We show that F's nonmonotonic behavior arises within the model from two opposing trends: the tendency for increased translation rates to produce more total protein but less cotranslationally folded protein. We further demonstrate that under certain conditions these nonmonotonic changes in F can result in nonmonotonic variations in post-translational, steady-state levels of functional protein. These results provide a potential explanation for recent experimental observations in which the specific activity of enzymatic proteins decreased with increased synthesis rates. Additionally our model has the potential to be used to rationally design transcripts to maximize the production of functional nascent protein by simultaneously optimizing translation initiation, elongation, and termination rates.