New scenarios of protein folding can occur on the ribosome

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.

Modulating co-translational protein folding by rational design and ribosome engineering

Co-translational folding is a fundamental process for the efficient biosynthesis of nascent polypeptides that emerge through the ribosome exit tunnel. To understand how this process is modulated by the shape and surface of the narrow tunnel, we have rationally engineered three exit tunnel protein loops (uL22, uL23 and uL24) of the 70S ribosome by CRISPR/Cas9 gene editing, and studied the co-translational folding of an immunoglobulin-like filamin domain (FLN5). Our thermodynamics measurements employing ¹⁹F/¹⁵N/methyl-TROSY NMR spectroscopy together with cryo-EM and molecular dynamics simulations reveal how the variations in the lengths of the loops present across species exert their distinct effects on the free energy of FLN5 folding. A concerted interplay of the uL23 and uL24 loops is sufficient to alter co-translational folding energetics, which we highlight by the opposite folding outcomes resulting from their extensions. These subtle modulations occur through a combination of the steric effects relating to the shape of the tunnel, the dynamic interactions between the ribosome surface and the unfolded nascent chain, and its altered exit pathway within the vestibule. These results illustrate the role of the exit tunnel structure in co-translational folding, and provide principles for how to remodel it to elicit a desired folding outcome.

The narrow exit tunnel of the ribosome is important for cotranslational protein folding. Here, authors show that their rationally designed and engineered exit tunnel protein loops modulate the free energy of nascent chain dynamics and folding.

The driving force for co-translational protein folding is weaker in the ribosome vestibule due to greater water ordering

Interactions between the ribosome and nascent chain can destabilize folded domains in the ribosome exit tunnel's vestibule, the last 3 nm of the exit tunnel where tertiary folding can occur. Here, we test if a contribution to this destabilization is a weakening of hydrophobic association, the driving force for protein folding. Using all-atom molecular dynamics simulations, we calculate the potential-of-mean force between two methane molecules along the center line of the ribosome exit tunnel and in bulk solution. Associated methanes, we find, are half as stable in the ribosome's vestibule as compared to bulk solution, demonstrating that the hydrophobic effect is weakened by the presence of the ribosome. This decreased stability arises from a decrease in the amount of water entropy gained upon the association of the methanes. And this decreased entropy gain originates from water molecules being more ordered in the vestibule as compared to bulk solution. Therefore, the hydrophobic effect is weaker in the vestibule because waters released from the first solvation shell of methanes upon association do not gain as much entropy in the vestibule as they do upon release in bulk solution. These findings mean that nascent proteins pass through a ribosome vestibule environment that can destabilize folded structures, which has the potential to influence co-translational protein folding pathways, energetics, and kinetics.

The ins and outs of Bacillus proteases: activities, functions and commercial significance

Because the majority of bacterial species divide by binary fission, and do not have distinguishable somatic and germline cells, they could be considered to be immortal. However, bacteria ‘age’ due to damage to vital cell components such as DNA and proteins. DNA damage can often be repaired using efficient DNA repair mechanisms. However, many proteins have a functional ‘shelf life’; some are short lived, while others are relatively stable. Specific degradation processes are built into the life span of proteins whose activities are required to fulfil a specific function during a prescribed period of time (e.g. cell cycle, differentiation process, stress response). In addition, proteins that are irreparably damaged or that have come to the end of their functional life span need to be removed by quality control proteases. Other proteases are involved in performing a variety of specific functions that can be broadly divided into three categories: processing, regulation and feeding. This review presents a systematic account of the proteases of Bacillus subtilis and their activities. It reviews the proteases found in, or associated with, the cytoplasm, the cell membrane, the cell wall and the external milieu. Where known, the impacts of the deletion of particular proteases are discussed, particularly in relation to industrial applications.

Nascent Folding of Proteins Across the Three Domains of Life

We study the nascent behavior of three model coarse-grained proteins in six rigid all-atom structures representing ribosomes that come from three domains of life. The synthesis of the proteins is implemented as a growth process. The geometry of the exit tunnel is quantified and shown to differ between the domains of life: both in volume and the size of constriction sites. This results in different characteristic times of capture within the tunnel and various probabilities of the escape. One of the proteins studied is the bacterial YibK which is knotted in its native state. A fraction of the trajectories results in knotting and the probability of doing so is largest for the bacterial ribosomes. Relaxing the condition of the rigidness of the ribosomes should result in a better avoidance of trapping and better proper folding.

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 ribosome modulates folding inside the ribosomal exit tunnel

Proteins commonly fold co-translationally at the ribosome, while the nascent chain emerges from the ribosomal exit tunnel. Protein domains that are sufficiently small can even fold while still located inside the tunnel. However, the effect of the tunnel on the folding dynamics of these domains is not well understood. Here, we combine optical tweezers with single-molecule FRET and molecular dynamics simulations to investigate folding of the small zinc-finger domain ADR1a inside and at the vestibule of the ribosomal tunnel. The tunnel is found to accelerate folding and stabilize the folded state, reminiscent of the effects of chaperonins. However, a simple mechanism involving stabilization by confinement does not explain the results. Instead, it appears that electrostatic interactions between the protein and ribosome contribute to the observed folding acceleration and stabilization of ADR1a.

Protein escape at the ribosomal exit tunnel: Effect of the tunnel shape

We study the post-translational escape of nascent proteins at the ribosomal exit tunnel with the consideration of a real shape atomistic tunnel based on the Protein Data Bank structure of the large ribosome subunit of archeon Haloarcula marismortui. Molecular dynamics simulations employing the Go-like model for the proteins show that at intermediate and high temperatures, including a presumable physiological temperature, the protein escape process at the atomistic tunnel is quantitatively similar to that at a cylinder tunnel of length L = 72 Å and diameter d = 16 Å. At low temperatures, the atomistic tunnel, however, yields an increased probability of protein trapping inside the tunnel, while the cylinder tunnel does not cause the trapping. All-β proteins tend to escape faster than all-α proteins, but this difference is blurred on increasing the protein's chain length. A 29-residue zinc-finger domain is shown to be severely trapped inside the tunnel. Most of the single-domain proteins considered, however, can escape efficiently at the physiological temperature with the escape time distribution following the diffusion model proposed in our previous works. An extrapolation of the simulation data to a realistic value of the friction coefficient for amino acids indicates that the escape times of globular proteins are at the sub-millisecond scale. It is argued that this time scale is short enough for the smooth functioning of the ribosome by not allowing nascent proteins to jam the ribosome tunnel.

Water-Mediated Interactions Determine Helix Formation of Peptides in Open Nanotubes

Water-mediated interactions (WMIs) play diverse roles in molecular biology. They are particularly relevant in geometrically confined spaces such as the interior of the chaperonin, at the interface between ligands and their binding partners, and in the ribosome tunnel. Inspired in part by the geometry of the ribosome tunnel, we consider confinement effects on the stability of peptides. We describe results from replica exchange molecular dynamics simulations of a system containing a 23-alanine or 23-serine polypeptide confined to nonpolar and polar nanotubes in the gas phase and when open to a water reservoir. We quantify the effect of water in determining the preferred conformational states of these polypeptides by calculating the difference in the solvation free energy for the helix and coil states in the open nanotube in the two phases. Our simulations reveal several possibilities. We find that nanoscopic confinement preferentially stabilizes the helical state of polypeptides with hydrophobic side chains, which is explained by the entropic stabilization mechanism proposed on the basis of polymer physics. Polypeptide chains with hydrophilic side chains can adopt helical structures within nanotubes, but helix formation is sensitive to the nature of the nanotube due to WMIs. We elaborate on the potential implications of our findings to the stability of peptides in the ribosome tunnel.

Interactions between nascent proteins and the ribosome surface inhibit co-translational folding

Most proteins begin to fold during biosynthesis on the ribosome. It has been suggested that interactions between the emerging polypeptide and the ribosome surface might allow the ribosome itself to modulate co-translational folding. Here we combine protein engineering and NMR spectroscopy to characterize a series of interactions between the ribosome surface and unfolded nascent chains of the immunoglobulin-like FLN5 filamin domain. The strongest interactions are found for a C-terminal segment that is essential for folding, and we demonstrate quantitative agreement between the strength of this interaction and the energetics of the co-translational folding process itself. Mutations in this region that reduce the extent of binding result in a shift in the co-translational folding equilibrium towards the native state. Our results therefore demonstrate that a competition between folding and binding provides a simple, dynamic mechanism for the modulation of co-translational folding by the ribosome.

Hidden kinetic traps in multidomain folding highlight the presence of a misfolded but functionally competent intermediate

Much of our current knowledge on protein folding is based on work focused on isolated domains. In this study, using a combination of NMR and kinetic experiments, we depict the folding pathway of a multidomain construct comprising two PDZ domains in tandem, belonging to the protein Whirlin. We demonstrate the presence of a misfolded intermediate that competes with productive folding. Interestingly, we show that, unexpectedly, this misfolded state retains the native-like functional ability to bind its physiological ligand, representing a clear example of a functionally competent misfolded state. On the basis of these results and a comparative analysis of the amino acidic sequences of Whirlin from different species, we propose a possible physiological role of the misfolded intermediate.

Although more than 75% of the proteome is composed of multidomain proteins, current knowledge of protein folding is based primarily on studies of isolated domains. In this work, we describe the folding mechanism of a multidomain tandem construct comprising two distinct covalently bound PDZ domains belonging to a protein called Whirlin, a scaffolding protein of the hearing apparatus. In particular, via a synergy between NMR and kinetic experiments, we demonstrate the presence of a misfolded intermediate that competes with productive folding. In agreement with the view that tandem domain swapping is a potential source of transient misfolding, we demonstrate that such a kinetic trap retains native-like functional activity, as shown by the preserved ability to bind its physiological ligand. Thus, despite the general knowledge that protein misfolding is intimately associated with dysfunction and diseases, we provide a direct example of a functionally competent misfolded state. Remarkably, a bioinformatics analysis of the amino acidic sequence of Whirlin from different species suggests that the tendency to perform tandem domain swapping between PDZ1 and PDZ2 is highly conserved, as demonstrated by their unexpectedly high sequence identity. On the basis of these observations, we discuss on a possible physiological role of such misfolded intermediate.

Cotranslational folding cooperativity of contiguous domains of α-spectrin

Proteins synthesized in the cell can begin to fold during translation before the entire polypeptide has been produced, which may be particularly relevant to the folding of multidomain proteins. Here, we study the cotranslational folding of adjacent domains from the cytoskeletal protein α-spectrin using force profile analysis (FPA). Specifically, we investigate how the cotranslational folding behavior of the R15 and R16 domains are affected by their neighboring R14 and R16, and R15 and R17 domains, respectively. Our results show that the domains impact each other's folding in distinct ways that may be important for the efficient assembly of α-spectrin, and may reduce its dependence on chaperones. Furthermore, we directly relate the experimentally observed yield of full-length protein in the FPA assay to the force exerted by the folding protein in piconewtons. By combining pulse-chase experiments to measure the rate at which the arrested protein is converted into full-length protein with a Bell model of force-induced rupture, we estimate that the R16 domain exerts a maximal force on the nascent chain of ∼15 pN during cotranslational folding.

Kinetic trapping in protein folding

The founding principles of protein folding introduced by Christian Anfinsen, together with the numerous mechanistic investigations that followed, assume that protein folding is a thermodynamically controlled process. On the other hand, this review underscores the fact that thermodynamic control is far from being the norm in protein folding, as long as one considers an extended chemical-potential landscape encompassing aggregates, in addition to native, unfolded and intermediate states. Here, we highlight the key role of kinetic trapping of the protein native state relative to unfolded, intermediate and, most importantly, aggregated states. We propose that kinetic trapping serves an important role in biology by protecting the bioactive states of a large number of proteins from deleterious aggregation. In the event that undesired aggregates were somehow formed, specialized intracellular disaggregation machines have evolved to convert any aberrant populations back to the native state, thus restoring a fully bioactive and aggregation-protected protein cohort.

Templated folding of intrinsically disordered proteins

Much of our current knowledge of biological chemistry is founded in the structure-function relationship, whereby sequence determines structure that determines function. Thus, the discovery that a large fraction of the proteome is intrinsically disordered, while being functional, has revolutionized our understanding of proteins and raised new and interesting questions. Many intrinsically disordered proteins (IDPs) have been determined to undergo a disorder-to-order transition when recognizing their physiological partners, suggesting that their mechanisms of folding are intrinsically different from those observed in globular proteins. However, IDPs also follow some of the classic paradigms established for globular proteins, pointing to important similarities in their behavior. In this review, we compare and contrast the folding mechanisms of globular proteins with the emerging features of binding-induced folding of intrinsically disordered proteins. Specifically, whereas disorder-to-order transitions of intrinsically disordered proteins appear to follow rules of globular protein folding, such as the cooperative nature of the reaction, their folding pathways are remarkably more malleable, due to the heterogeneous nature of their folding nuclei, as probed by analysis of linear free-energy relationship plots. These insights have led to a new model for the disorder-to-order transition in IDPs termed "templated folding," whereby the binding partner dictates distinct structural transitions en route to product, while ensuring a cooperative folding.

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.

Differences in the path to exit the ribosome across the three domains of life

The ribosome exit tunnel is an important structure involved in the regulation of translation and other essential functions such as protein folding. By comparing 20 recently obtained cryo-EM and X-ray crystallography structures of the ribosome from all three domains of life, we here characterize the key similarities and differences of the tunnel across species. We first show that a hierarchical clustering of tunnel shapes closely reflects the species phylogeny. Then, by analyzing the ribosomal RNAs and proteins, we explain the observed geometric variations and show direct association between the conservations of the geometry, structure and sequence. We find that the tunnel is more conserved in the upper part close to the polypeptide transferase center, while in the lower part, it is substantially narrower in eukaryotes than in bacteria. Furthermore, we provide evidence for the existence of a second constriction site in eukaryotic exit tunnels. Overall, these results have several evolutionary and functional implications, which explain certain differences between eukaryotes and prokaryotes in their translation mechanisms. In particular, they suggest that major co-translational functions of bacterial tunnels were externalized in eukaryotes, while reducing the tunnel size provided some other advantages, such as facilitating the nascent chain elongation and enabling antibiotic resistance.

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.

The shape of the bacterial ribosome exit tunnel affects cotranslational protein folding

The E. coli ribosome exit tunnel can accommodate small folded proteins, while larger ones fold outside. It remains unclear, however, to what extent the geometry of the tunnel influences protein folding. Here, using E. coli ribosomes with deletions in loops in proteins uL23 and uL24 that protrude into the tunnel, we investigate how tunnel geometry determines where proteins of different sizes fold. We find that a 29-residue zinc-finger domain normally folding close to the uL23 loop folds deeper in the tunnel in uL23 Δloop ribosomes, while two ~ 100 residue proteins normally folding close to the uL24 loop near the tunnel exit port fold at deeper locations in uL24 Δloop ribosomes, in good agreement with results obtained by coarse-grained molecular dynamics simulations. This supports the idea that cotranslational folding commences once a protein domain reaches a location in the exit tunnel where there is sufficient space to house the folded structure.

Folding pathway of an Ig domain is conserved on and off the ribosome

Most proteins need to fold into a specific 3D structure to function. The mechanism by which isolated proteins fold has been thoroughly studied by experiment and theory. However, in the cell proteins do not fold in isolation but are synthesized as linear chains by the ribosome during translation. It is therefore natural to ask at which point during synthesis proteins fold, and whether this differs from the folding of isolated protein molecules. By studying folding of a well-characterized protein domain, titin I27, stalled at different points during translation, we show that it already folds in the mouth of the ribosome exit tunnel and that the mechanism is almost identical to that of the isolated protein.

Proteins that fold cotranslationally may do so in a restricted configurational space, due to the volume occupied by the ribosome. How does this environment, coupled with the close proximity of the ribosome, affect the folding pathway of a protein? Previous studies have shown that the cotranslational folding process for many proteins, including small, single domains, is directly affected by the ribosome. Here, we investigate the cotranslational folding of an all-β Ig domain, titin I27. Using an arrest peptide-based assay and structural studies by cryo-EM, we show that I27 folds in the mouth of the ribosome exit tunnel. Simulations that use a kinetic model for the force dependence of escape from arrest accurately predict the fraction of folded protein as a function of length. We used these simulations to probe the folding pathway on and off the ribosome. Our simulations—which also reproduce experiments on mutant forms of I27—show that I27 folds, while still sequestered in the mouth of the ribosome exit tunnel, by essentially the same pathway as free I27, with only subtle shifts of critical contacts from the C to the N terminus.

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.

Effects of protein size, thermodynamic stability, and net charge on cotranslational folding on the ribosome

During the last five decades, studies of protein folding in dilute buffer solutions have produced a rich picture of this complex process. In the cell, however, proteins can start to fold while still attached to the ribosome (cotranslational folding) and it is not yet clear how the ribosome affects the folding of protein domains of different sizes, thermodynamic stabilities, and net charges. Here, by using arrest peptides as force sensors and on-ribosome pulse proteolysis, we provide a comprehensive picture of how the distance from the peptidyl transferase center in the ribosome at which proteins fold correlates with protein size. Moreover, an analysis of a large collection of mutants of the Escherichia coli ribosomal protein S6 shows that the force exerted on the nascent chain by protein folding varies linearly with the thermodynamic stability of the folded state, and that the ribosome environment disfavors folding of domains of high net-negative charge.

Charge Interactions Can Dominate Coupled Folding and Binding on the Ribosome

Interactions between emerging nascent polypeptide chains and the ribosome can modulate cotranslational protein folding. However, it has remained unclear how such interactions can affect the binding of nascent chains to their cellular targets. We thus investigated on the ribosome the interaction between two intrinsically disordered proteins of opposite charge, ACTR and NCBD, which form a high-affinity complex in a coupled folding-and-binding reaction. Using fluorescence correlation spectroscopy and arrest-peptide-mediated force measurements in vitro and in vivo, we find that the ACTR-NCBD complex can form cotranslationally but only with ACTR as the nascent chain and NCBD free in solution, not vice versa. We show that this surprising asymmetry in behavior is caused by pronounced charge interactions: attraction of the positively charged nascent chain of NCBD to the negatively charged ribosomal surface competes with complex formation and prevents ACTR binding. In contrast, the negatively charged nascent ACTR is repelled by the ribosomal surface and thus remains available for productively binding its partner. Electrostatic interactions may thus be more important for cotranslational folding and binding than previously thought.

Co-translational folding of α-helical proteins: structural studies of intermediate-length variants of the λ repressor

Nascent polypeptide chains fold cotranslationally, but the atomic‐level details of this process remain unknown. Here, we report crystallographic, de novo modeling, and spectroscopic studies of intermediate‐length variants of the λ repressor N‐terminal domain. Although the ranges of helical regions of the half‐length variant were almost identical to those of the full‐length protein, the relative orientations of these helices in the intermediate‐length variants differed. Our results suggest that cotranslational folding of the λ repressor initially forms a helical structure with a transient conformation, as in the case of a molten globule state. This conformation subsequently matures during the course of protein synthesis.

Database

Structural data are available in the PDB under the accession numbers http://www.rcsb.org/pdb/search/structidSearch.do?structureId=5ZCA and http://www.rcsb.org/pdb/search/structidSearch.do?structureId=3WOA.

Kinetic and structural comparison of a protein's cotranslational folding and refolding pathways

Precise protein folding is essential for the survival of all cells, and protein misfolding causes a number of diseases that lack effective therapies, yet the general principles governing protein folding in the cell remain poorly understood. In vivo, folding can begin cotranslationally and protein quality control at the ribosome is essential for cellular proteostasis. We directly characterize and compare the refolding and cotranslational folding trajectories of the protein HaloTag. We introduce new techniques for both measuring folding kinetics and detecting the conformations of partially folded intermediates during translation in real time. We find that, although translation does not affect the rate-limiting step of HaloTag folding, a key aggregation-prone intermediate observed during in vitro refolding experiments is no longer detectable. This rerouting of the folding pathway increases HaloTag's folding efficiency and may serve as a general chaperone-independent mechanism of quality control by the ribosome.

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.

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.

Rotational restriction of nascent peptides as an essential element of co-translational protein folding: possible molecular players and structural consequences

Background

A basic tenet of protein science is that all information about the spatial structure of proteins is present in their sequences. Nonetheless, many proteins fail to attain native structure upon experimental denaturation and refolding in vitro, raising the question of the specific role of cellular machinery in protein folding in vivo. Recently, we hypothesized that energy-dependent twisting of the protein backbone is an unappreciated essential factor guiding the protein folding process in vivo. Torque force may be applied by the ribosome co-translationally, and when accompanied by simultaneous restriction of the rotational mobility of the distal part of the growing chain, the resulting tension in the protein backbone would facilitate the formation of local secondary structure and direct the folding process.

Results

Our model of the early stages of protein folding in vivo postulates that the free motion of both terminal regions of the protein during its synthesis and maturation is restricted. The long-known but unexplained phenomenon of statistical overrepresentation of protein termini on the surfaces of the protein structures may be an indication of the backbone twist-based folding mechanism; sustained maintenance of a twist requires that both ends of the protein chain are anchored in space, and if the ends are released only after the majority of folding is complete, they are much more likely to remain on the surface of the molecule. We identified the molecular components that are likely to play a role in the twisting of the nascent protein chain and in the anchoring of its N-terminus. The twist may be induced at the C-terminus of the nascent polypeptide by the peptidyltransferase center of the ribosome. Several ribosome-associated proteins, including the trigger factor in bacteria and the nascent polypeptide-associated complex in archaea and eukaryotes, may restrict the rotational mobility of the N-proximal regions of the peptides.

Conclusions

Many experimental observations are consistent with the hypothesis of co-translational twisting of the protein backbone. Several molecular players in this hypothetical mechanism of protein folding can be suggested. In addition, the new view of protein folding in vivo opens the possibility of novel potential drug targets to combat human protein folding diseases.

Reviewers

This article was reviewed by Lakshminarayan Iyer and István Simon.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-017-0186-1) contains supplementary material, which is available to authorized users.

Translation and folding of single proteins in real time

Protein biosynthesis is inherently coupled to cotranslational protein folding. Folding of the nascent chain already occurs during synthesis and is mediated by spatial constraints imposed by the ribosomal exit tunnel as well as self-interactions. The polypeptide's vectorial emergence from the ribosomal tunnel establishes the possible folding pathways leading to its native tertiary structure. How cotranslational protein folding and the rate of synthesis are linked to a protein's amino acid sequence is still not well defined. Here, we follow synthesis by individual ribosomes using dual-trap optical tweezers and observe simultaneous folding of the nascent polypeptide chain in real time. We show that observed stalling during translation correlates with slowed peptide bond formation at successive proline sequence positions and electrostatic interactions between positively charged amino acids and the ribosomal tunnel. We also determine possible cotranslational folding sites initiated by hydrophobic collapse for an unstructured and two globular proteins while directly measuring initial cotranslational folding forces. Our study elucidates the intricate relationship among a protein's amino acid sequence, its cotranslational nascent-chain elongation rate, and folding.

Effect of Nascent Peptide Steric Bulk on Elongation Kinetics in the Ribosome Exit Tunnel

All proteins are synthesized by the ribosome, a macromolecular complex that accomplishes the life-sustaining tasks of faithfully decoding mRNA and catalyzing peptide bond formation at the peptidyl transferase center (PTC). The ribosome has evolved an exit tunnel to host the elongating new peptide, protect it from proteolytic digestion, and guide its emergence. It is here that the nascent chain begins to fold. This folding process depends on the rate of translation at the PTC. We report here that besides PTC events, translation kinetics depend on steric constraints on nascent peptide side chains and that confined movements of cramped side chains within and through the tunnel fine-tune elongation rates.

The role of the backbone torsion in protein folding

BACKGROUND

The set of forces and sequence of events that govern the transition from an unfolded polypeptide chain to a functional protein with correct spatial structure remain incompletely known, despite the importance of the problem and decades of theory development, computer simulations, and laboratory experiments. Information about the correctly folded state of most proteins is likely to be present in their sequences, and yet many proteins fail to attain native structure after overexpression in a non-native environment or upon experimental denaturation and refolding.

PRESENTATION OF THE HYPOTHESIS

We hypothesize that correct protein folding in vivo is an active, energy-dependent process that most likely applies torque force co-translationally to all proteins and possibly also post-translationally to many proteins in every cell. When a site on an unfolded polypeptide is rotationally constrained, torsion applied at another site would induce twisting of the main chain, which would initiate the formation of a local secondary structure, such as an alpha-helical turn or a beta-turn/beta-hairpin. The nucleation of structural elements is a rate-limiting, energetically unfavorable step in the process of protein folding, and energy-dependent chain torsion is likely to help overcome this barrier in vivo. Several molecular machines in a cell, primarily ribosomes, but also possibly signal recognition particles and chaperone systems, may play a role in applying torque to an unfolded protein chain, using the energy of GTP or ATP hydrolysis. Lack of such force in the in vitro systems may be the main reason of the failure of many longer proteins to attain the correct functional conformation.

TESTING THE HYPOTHESIS

The hypothesis can be tested using single-molecule approaches, by measuring directly the forces applied to polypeptide chains under controlled conditions in vitro, and in bulk, by assessing folding rates and extent of misfolding in proteins that are engineered to experience transient spatial constraint during their synthesis.

IMPLICATIONS OF THE HYPOTHESIS

Learning about the role of main chain torsion in protein folding will improve our understanding of folding mechanisms and may lead to bioengineering solutions that would enhance the yields of correctly folded proteins in heterologous expression systems.

REVIEWERS

This article was reviewed by Frank Eisenhaber, Igor Berezovsky and Michael Gromiha.

Quantitative determination of ribosome nascent chain stability

Accurate protein folding is essential for proper cellular and organismal function. In the cell, protein folding is carefully regulated; changes in folding homeostasis (proteostasis) can disrupt many cellular processes and have been implicated in various neurodegenerative diseases and other pathologies. For many proteins, the initial folding process begins during translation while the protein is still tethered to the ribosome; however, most biophysical studies of a protein's energy landscape are carried out in isolation under idealized, dilute conditions and may not accurately report on the energy landscape in vivo. Thus, the energy landscape of ribosome nascent chains and the effect of the tethered ribosome on nascent chain folding remain unclear. Here we have developed a general assay for quantitatively measuring the folding stability of ribosome nascent chains, and find that the ribosome exerts a destabilizing effect on the polypeptide chain. This destabilization decreases as a function of the distance away from the peptidyl transferase center. Thus, the ribosome may add an additional layer of robustness to the protein-folding process by avoiding the formation of stable partially folded states before the protein has completely emerged from the ribosome.

Basic units of protein structure, folding, and function

Study of the hierarchy of domain structure with alternative sets of domains and analysis of discontinuous domains, consisting of remote segments of the polypeptide chain, raised a question about the minimal structural unit of the protein domain. The hypothesis on the decisive role of the polypeptide backbone in determining the elementary units of globular proteins have led to the discovery of closed loops. It is reviewed here how closed loops form the loop-n-lock structure of proteins, providing the foundation for stability and designability of protein folds/domain and underlying their co-translational folding. Simplified protein sequences are considered here with the aim to explore the basic principles that presumably dominated the folding and stability of proteins in the early stages of structural evolution. Elementary functional loops (EFLs), closed loops with one or few catalytic residues, are, in turn, units of the protein function. They are apparent descendants of the prebiotic ring-like peptides, which gave rise to the first functional folds/domains being fused in the beginning of the evolution of protein structure. It is also shown how evolutionary relations between protein functional superfamilies and folds delineated with the help of EFLs can contribute to establishing the rules for design of desired enzymatic functions. Generalized descriptors of the elementary functions are proposed to be used as basic units in the future computational design.

In vivo aspects of protein folding and quality control

Most proteins must fold into unique three-dimensional structures to perform their biological functions. In the crowded cellular environment, newly synthesized proteins are at risk of misfolding and forming toxic aggregate species. To ensure efficient folding, different classes of molecular chaperones receive the nascent protein chain emerging from the ribosome and guide it along a productive folding pathway. Because proteins are structurally dynamic, constant surveillance of the proteome by an integrated network of chaperones and protein degradation machineries is required to maintain protein homeostasis (proteostasis). The capacity of this proteostasis network declines during aging, facilitating neurodegeneration and other chronic diseases associated with protein aggregation. Understanding the proteostasis network holds the promise of identifying targets for pharmacological intervention in these pathologies.

Insights into Cotranslational Nascent Protein Behavior from Computer Simulations

Regulation of protein stability and function in vivo begins during protein synthesis, when the ribosome translates a messenger RNA into a nascent polypeptide. Cotranslational processes involving a nascent protein include folding, binding to other macromolecules, enzymatic modification, and secretion through membranes. Experiments have shown that the rate at which the ribosome adds amino acids to the elongating nascent chain influences the efficiency of these processes, with alterations to these rates possibly contributing to diseases, including some types of cancer. In this review, we discuss recent insights into cotranslational processes gained from molecular simulations, how different computational approaches have been combined to understand cotranslational processes at multiple scales, and the new scenarios illuminated by these simulations. We conclude by suggesting interesting questions that computational approaches in this research area can address over the next few years.

Ribosomal Stalk Protein Silencing Partially Corrects the ΔF508-CFTR Functional Expression Defect

The most common cystic fibrosis (CF) causing mutation, deletion of phenylalanine 508 (ΔF508 or Phe508del), results in functional expression defect of the CF transmembrane conductance regulator (CFTR) at the apical plasma membrane (PM) of secretory epithelia, which is attributed to the degradation of the misfolded channel at the endoplasmic reticulum (ER). Deletion of phenylalanine 670 (ΔF670) in the yeast oligomycin resistance 1 gene (YOR1, an ABC transporter) of Saccharomyces cerevisiae phenocopies the ΔF508-CFTR folding and trafficking defects. Genome-wide phenotypic (phenomic) analysis of the Yor1-ΔF670 biogenesis identified several modifier genes of mRNA processing and translation, which conferred oligomycin resistance to yeast. Silencing of orthologues of these candidate genes enhanced the ΔF508-CFTR functional expression at the apical PM in human CF bronchial epithelia. Although knockdown of RPL12, a component of the ribosomal stalk, attenuated the translational elongation rate, it increased the folding efficiency as well as the conformational stability of the ΔF508-CFTR, manifesting in 3-fold augmented PM density and function of the mutant. Combination of RPL12 knockdown with the corrector drug, VX-809 (lumacaftor) restored the mutant function to ~50% of the wild-type channel in primary CFTR^{ΔF508/ΔF508} human bronchial epithelia. These results and the observation that silencing of other ribosomal stalk proteins partially rescue the loss-of-function phenotype of ΔF508-CFTR suggest that the ribosomal stalk modulates the folding efficiency of the mutant and is a potential therapeutic target for correction of the ΔF508-CFTR folding defect.

Reducing the rate of translational elongation by silencing ribosomal stalk proteins ameliorates the folding and stability defect of the cystic fibrosis mutant protein ΔF508-CFTR, partially restoring the plasma membrane chloride conductance.

Cystic fibrosis (CF) is one of the most common autosomal recessive diseases in Caucasians. It is caused by mutations in the CF transmembrane conductance regulator (CFTR), which functions as an anion channel at the apical plasma membrane of secretory epithelia. The most common CF mutation, a deletion of the phenylalanine residue at position 508 (ΔF508), results in the channel misfolding and subsequent intracellular degradation. Our previous genome-wide phenotypic screens, using a yeast variant, have predicted modifier genes for ΔF508-CFTR biogenesis. Here, we show that silencing of one of these candidate genes, RPL12, a component of the ribosomal stalk, increased the folding and stabilization of ΔF508-CFTR, resulting in its increased plasma membrane expression and function. Our data suggest that reducing the translational elongation rate via RPL12 silencing can partially reverse the ΔF508-CFTR folding defect. Importantly, RPL12 silencing in combination with the corrector drug VX-809 (lumacaftor), increased the mutant function to 50% of the wild-type CFTR channel, suggesting that the ribosomal stalk perturbation may represent a therapeutic target for rescuing the ΔF508-CFTR biogenesis defect.

Modeling the effect of codon translation rates on co-translational protein folding mechanisms of arbitrary complexity

In a cell, the folding of a protein molecule into tertiary structure can begin while it is synthesized by the ribosome. The rate at which individual amino acids are incorporated into the elongating nascent chain has been shown to affect the likelihood that proteins will populate their folded state, indicating that co-translational protein folding is a far from equilibrium process. Developing a theoretical framework to accurately describe this process is, therefore, crucial for advancing our understanding of how proteins acquire their functional conformation in living cells. Current state-of-the-art computational approaches, such as molecular dynamics simulations, are very demanding in terms of the required computer resources, making the simulation of co-translational protein folding difficult. Here, we overcome this limitation by introducing an efficient approach that predicts the effects that variable codon translation rates have on co-translational folding pathways. Our approach is based on Markov chains. By using as an input a relatively small number of molecular dynamics simulations, it allows for the computation of the probability that a nascent protein is in any state as a function of the translation rate of individual codons along a mRNA's open reading frame. Due to its computational efficiency and favorable scalability with the complexity of the folding mechanism, this approach could enable proteome-wide computational studies of the influence of translation dynamics on co-translational folding.

A structural ensemble of a ribosome-nascent chain complex during cotranslational protein folding

Although detailed pictures of ribosome structures are emerging, little is known about the structural and cotranslational folding properties of nascent polypeptide chains at the atomic level. Here we used solution-state NMR spectroscopy to define a structural ensemble of a ribosome-nascent chain complex (RNC) formed during protein biosynthesis in Escherichia coli, in which a pair of immunoglobulin-like domains adopts a folded N-terminal domain (FLN5) and a disordered but compact C-terminal domain (FLN6). To study how FLN5 acquires its native structure cotranslationally, we progressively shortened the RNC constructs. We found that the ribosome modulates the folding process, because the complete sequence of FLN5 emerged well beyond the tunnel before acquiring native structure, whereas FLN5 in isolation folded spontaneously, even when truncated. This finding suggests that regulating structure acquisition during biosynthesis can reduce the probability of misfolding, particularly of homologous domains.

Accurate prediction of cellular co-translational folding indicates proteins can switch from post- to co-translational folding

The rates at which domains fold and codons are translated are important factors in determining whether a nascent protein will co-translationally fold and function or misfold and malfunction. Here we develop a chemical kinetic model that calculates a protein domain's co-translational folding curve during synthesis using only the domain's bulk folding and unfolding rates and codon translation rates. We show that this model accurately predicts the course of co-translational folding measured in vivo for four different protein molecules. We then make predictions for a number of different proteins in yeast and find that synonymous codon substitutions, which change translation-elongation rates, can switch some protein domains from folding post-translationally to folding co-translationally--a result consistent with previous experimental studies. Our approach explains essential features of co-translational folding curves and predicts how varying the translation rate at different codon positions along a transcript's coding sequence affects this self-assembly process.

Identification of fibrillogenic regions in human triosephosphate isomerase

Background. Amyloid secondary structure relies on the intermolecular assembly of polypeptide chains through main-chain interaction. According to this, all proteins have the potential to form amyloid structure, nevertheless, in nature only few proteins aggregate into toxic or functional amyloids. Structural characteristics differ greatly among amyloid proteins reported, so it has been difficult to link the fibrillogenic propensity with structural topology. However, there are ubiquitous topologies not represented in the amyloidome that could be considered as amyloid-resistant attributable to structural features, such is the case of TIM barrel topology.

Methods. This work was aimed to study the fibrillogenic propensity of human triosephosphate isomerase (HsTPI) as a model of TIM barrels. In order to do so, aggregation of HsTPI was evaluated under native-like and destabilizing conditions. Fibrillogenic regions were identified by bioinformatics approaches, protein fragmentation and peptide aggregation.

Results. We identified four fibrillogenic regions in the HsTPI corresponding to the β3, β6, β7 y α8 of the TIM barrel. From these, the β3-strand region (residues 59–66) was highly fibrillogenic. In aggregation assays, HsTPI under native-like conditions led to amorphous assemblies while under partially denaturing conditions (urea 3.2 M) formed more structured aggregates. This slightly structured aggregates exhibited residual cross-β structure, as demonstrated by the recognition of the WO1 antibody and ATR-FTIR analysis.

Discussion. Despite the fibrillogenic regions present in HsTPI, the enzyme maintained under native-favoring conditions displayed low fibrillogenic propensity. This amyloid-resistance can be attributed to the three-dimensional arrangement of the protein, where β-strands, susceptible to aggregation, are protected in the core of the molecule. Destabilization of the protein structure may expose inner regions promoting β-aggregation, as well as the formation of hydrophobic disordered aggregates. Being this last pathway kinetically favored over the thermodynamically more stable fibril aggregation pathway.

Physical Origins of Codon Positions That Strongly Influence Cotranslational Folding: A Framework for Controlling Nascent-Protein Folding

An emerging paradigm in the field of in vivo protein biophysics is that nascent-protein behavior is a type of nonequilibrium phenomenon, where translation-elongation kinetics can be more important in determining nascent-protein behavior than the thermodynamic properties of the protein. Synonymous codon substitutions, which change the translation rate at select codon positions along a transcript, have been shown to alter cotranslational protein folding, suggesting that evolution may have shaped synonymous codon usage in the genomes of organisms in part to increase the amount of folded and functional nascent protein. Here, we develop a Monte Carlo-master-equation method that allows for the control of nascent-chain folding during translation through the rational design of mRNA sequences to guide the cotranslational folding process. We test this framework using coarse-grained molecular dynamics simulations and find it provides optimal mRNA sequences to control the simulated, cotranslational folding of a protein in a user-prescribed manner. With this approach we discover that some codon positions in a transcript can have a much greater impact on nascent-protein folding than others because they tend to be positions where the nascent chain populates states that are far from equilibrium, as well as being dependent on a complex ratio of time scales. As a consequence, different cotranslational profiles of the same protein can have different critical codon positions and different numbers of synonymous mRNA sequences that encode for them. These findings explain that there is a fundamental connection between the nonequilibrium nature of cotranslational processes, nascent-protein behavior, and synonymous codon usage.

Cotranslational protein folding on the ribosome monitored in real time

Protein domains can fold into stable tertiary structures while they are synthesized on the ribosome. We used a high-performance, reconstituted in vitro translation system to investigate the folding of a small five-helix protein domain-the N-terminal domain of Escherichia coli N5-glutamine methyltransferase HemK-in real time. Our observations show that cotranslational folding of the protein, which folds autonomously and rapidly in solution, proceeds through a compact, non-native conformation that forms within the peptide tunnel of the ribosome. The compact state rearranges into a native-like structure immediately after the full domain sequence has emerged from the ribosome. Both folding transitions are rate-limited by translation, allowing for quasi-equilibrium sampling of the conformational space restricted by the ribosome. Cotranslational folding may be typical of small, intrinsically rapidly folding protein domains.

Cotranslational Protein Folding inside the Ribosome Exit Tunnel

At what point during translation do proteins fold? It is well established that proteins can fold cotranslationally outside the ribosome exit tunnel, whereas studies of folding inside the exit tunnel have so far detected only the formation of helical secondary structure and collapsed or partially structured folding intermediates. Here, using a combination of cotranslational nascent chain force measurements, inter-subunit fluorescence resonance energy transfer studies on single translating ribosomes, molecular dynamics simulations, and cryoelectron microscopy, we show that a small zinc-finger domain protein can fold deep inside the vestibule of the ribosome exit tunnel. Thus, for small protein domains, the ribosome itself can provide the kind of sheltered folding environment that chaperones provide for larger proteins.

•

Cotranslational folding is studied by arrest-peptide-mediated force measurements

•

Single-molecule measurements show that a pulling force prevents ribosome stalling

•

A ribosome-tethered zinc-finger domain is visualized by cryo-EM (electron microscopy)

•

The zinc-finger domain is shown to fold deep inside the ribosome exit tunnel

How Co-translational Folding of Multi-domain Protein Is Affected by Elongation Schedule: Molecular Simulations

Co-translational folding (CTF) facilitates correct folding in vivo, but its precise mechanism remains elusive. For the CTF of a three-domain protein SufI, it was reported that the translational attenuation is obligatory to acquire the functional state. Here, to gain structural insights on the underlying mechanisms, we performed comparative molecular simulations of SufI that mimic CTF as well as refolding schemes. A CTF scheme that relied on a codon-based prediction of translational rates exhibited folding probability markedly higher than that by the refolding scheme. When the CTF schedule is speeded up, the success rate dropped. These agree with experiments. Structural investigation clarified that misfolding of the middle domain was much more frequent in the refolding scheme than that in the codon-based CTF scheme. The middle domain is less stable and can fold via interactions with the folded N-terminal domain. Folding pathway networks showed the codon-based CTF gives narrower pathways to the native state than the refolding scheme.

Proteins are synthesized in vivo by ribosome from their N-termini. When N-terminal fragments of nascent proteins get out of the ribosome exit, they start folding, which is called co-translational folding. It has been suggested that well-scheduled co-translational folding schemes would facilitate correct acquisition of their native structures for some multi-domain proteins. In particular, an un-ambiguous experiment was recently reported for a model protein, SufI where pauses at certain positions in the translational elongation are obligatory for efficient folding. Here, for the first time to our knowledge, we performed molecular dynamics simulations of SufI with co-translational folding as well as re-folding schemes. We found a co-translational folding shceme with rare codon-based pauses indeed increased the success ratio of folding, which is consistent with recent experiments. On top, molecular simulations provided much of structural insights on the folding routes and misfolding in the case of re-folding scheme. This explains why pauses in the translational elongation rescue SufI from misfolding.