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

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.

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.

The protein escape process at the ribosomal exit tunnel has conserved mechanisms across the domains of life

The ribosomal exit tunnel is the primary structure affecting the release of nascent proteins at the ribosome. The ribosomal exit tunnels from different species have elements of conservation and differentiation in structural and physico-chemical properties. In this study, by simulating the elongation and escape processes of nascent proteins at the ribosomal exit tunnels of four different organisms, we show that the escape process has conserved mechanisms across the domains of life. Specifically, it is found that the escape process of proteins follows the diffusion mechanism given by a simple diffusion model, and the median escape time positively correlates with the number of hydrophobic residues and the net charge of a protein for all the exit tunnels considered. These properties hold for 12 distinct proteins considered in two slightly different and improved Gō-like models. It is also found that the differences in physico-chemical properties of the tunnels lead to quantitative differences in the protein escape times. In particular, the relatively strong hydrophobicity of E. coli's tunnel and the unusually high number of negatively charged amino acids on the tunnel's surface of H. marismortui lead to substantially slower escapes of proteins at these tunnels than at those of S. cerevisiae and H. sapiens.

Extracting time series matching a small-angle X-ray scattering profile from trajectories of molecular dynamics simulations

Solving structural ensembles of flexible biomolecules is a challenging research area. Here, we propose a method to obtain possible structural ensembles of a biomolecule based on small-angle X-ray scattering (SAXS) and molecular dynamics simulations. Our idea is to clip a time series that matches a SAXS profile from a simulation trajectory. To examine its practicability, we applied our idea to a multi-domain protein ER-60 and successfully extracted time series longer than 1 micro second from trajectories of coarse-grained molecular dynamics simulations. In the extracted time series, the domain conformation was distributed continuously and smoothly in a conformational space. Preferred domain conformations were also observed. Diversity among scattering curves calculated from each ER-60 structure was interpreted to reflect an open-close motion of the protein. Although our approach did not provide a unique solution for the structural ensemble of the biomolecule, each extracted time series can be an element of the real behavior of ER-60. Considering its low computational cost, our approach will play a key role to identify biomolecular dynamics by integrating SAXS, simulations, and other experiments.

Hydrophobic and electrostatic interactions modulate protein escape at the ribosomal exit tunnel

After translation, nascent proteins must escape the ribosomal exit tunnel to attain complete folding to their native states. This escape process also frees up the ribosome tunnel for a new translation job. In this study, we investigate the impacts of energetic interactions between the ribosomal exit tunnel and nascent proteins on the protein escape process by molecular dynamics simulations using partially coarse-grained models that incorporate hydrophobic and electrostatic interactions of the ribosome tunnel of Haloarcula marismortui with nascent proteins. We find that, in general, attractive interactions slow down the protein escape process, whereas repulsive interactions speed it up. For the small globular proteins considered, the median escape time correlates with both the number of hydrophobic residues, N_h, and the net charge, Q, of a nascent protein. A correlation coefficient exceeding 0.96 is found for the relation between the median escape time and a combined quantity of N_h + 5.9Q, suggesting that it is ∼6 times more efficient to modulate the escape time by changing the total charge than the number of hydrophobic residues. The estimated median escape times are found in the submillisecond-to-millisecond range, indicating that the escape does not delay the ribosome recycling. For various types of the tunnel model, with and without hydrophobic and electrostatic interactions, the escape time distribution always follows a simple diffusion model that describes the escape process as a downhill drift of a Brownian particle, suggesting that nascent proteins escape along barrier-less pathways at the ribosome tunnel.

Multiscale Coarse-Grained Model for the Stepping of Molecular Motors with Application to Kinesin

Conventional kinesin, a motor protein that transports cargo within cells, walks by taking multiple steps toward the plus end of the microtubule (MT). While significant progress has been made in understanding the details of the walking mechanism of kinesin, there are many unresolved issues. From a computational perspective, a central challenge is the large size of the system, which limits the scope of time scales accessible in standard computer simulations. Here, we create a general multiscale coarse-grained model for motors that enables us to simulate the stepping process of motors on polar tracks (actin and MT) with a focus on kinesin. Our approach greatly shortens the computation times without a significant loss in detail, thus allowing us to better describe the molecular basis of the stepping kinetics. The small number of parameters, which are determined by fitting to experimental data, allows us to develop an accurate method that may be adopted to simulate stepping in other molecular motors. The model enables us to simulate a large number of steps, which was not possible previously. We show in agreement with experiments that due to the docking of the neck linker (NL) of kinesin, sometimes deemed as the power stroke, the space explored diffusively by the tethered head is severely restricted, allowing the step to be completed in tens of microseconds. We predict that increasing the interaction strength between the NL and the motor head, achievable by mutations in the NL, decreases the stepping time but reaches a saturation value. Furthermore, the full three-dimensional dynamics of the cargo are fully resolved in our model, contributing to the predictive power and allowing us to study the important aspects of cargo-motor interactions.

Energy landscape remodeling mechanism of Hsp70-chaperone-accelerated protein folding

Hsp70 chaperone is one of the key protein machines responsible for the quality control of protein production in cells. Facilitating in vivo protein folding by counteracting misfolding and aggregation is the essence of its biological function. Although the allosteric cycle during its functional actions has been well characterized both experimentally and computationally, the mechanism by which Hsp70 assists protein folding is still not fully understood. In this work, we studied the Hsp70-mediated folding of model proteins with rugged energy landscape by using molecular simulations. Different from the canonical scenario of Hsp70 functioning, which assumes that folding of substrate proteins occurs spontaneously after releasing from chaperones, our results showed that the substrate protein remains in contacts with the chaperone during its folding process. The direct chaperone-substrate interactions in the open conformation of Hsp70 tend to shield the substrate sites prone to form non-native contacts, which therefore avoids the frustrated folding pathway, leading to a higher folding rate and less probability of misfolding. Our results suggest that in addition to the unfoldase and holdase functions widely addressed in previous studies, Hsp70 can facilitate the folding of its substrate proteins by remodeling the folding energy landscape and directing the folding processes, demonstrating the foldase scenario. These findings add new, to our knowledge, insights into the general molecular mechanisms of chaperone-mediated protein folding.

Solution structure of multi-domain protein ER-60 studied by aggregation-free SAXS and coarse-grained-MD simulation

Multi-domain proteins (MDPs) show a variety of domain conformations under physiological conditions, regulating their functions through such conformational changes. One of the typical MDPs, ER-60 which is a protein folding enzyme, has a U-shape with four domains and is thought to have different domain conformations in solution depending on the redox state at the active centres of the edge domains. In this work, an aggregation-free small-angle X-ray scattering revealed that the structures of oxidized and reduced ER-60 in solution are different from each other and are also different from those in the crystal. Furthermore, structural modelling with coarse-grained molecular dynamics simulation indicated that the distance between the two edge domains of oxidized ER-60 is longer than that of reduced ER-60. In addition, one of the edge domains has a more flexible conformation than the other.

Investigating the trade-off between folding and function in a multidomain Y-family DNA polymerase

The way in which multidomain proteins fold has been a puzzling question for decades. Until now, the mechanisms and functions of domain interactions involved in multidomain protein folding have been obscure. Here, we develop structure-based models to investigate the folding and DNA-binding processes of the multidomain Y-family DNA polymerase IV (DPO4). We uncover shifts in the folding mechanism among ordered domain-wise folding, backtracking folding, and cooperative folding, modulated by interdomain interactions. These lead to ‘U-shaped’ DPO4 folding kinetics. We characterize the effects of interdomain flexibility on the promotion of DPO4–DNA (un)binding, which probably contributes to the ability of DPO4 to bypass DNA lesions, which is a known biological role of Y-family polymerases. We suggest that the native topology of DPO4 leads to a trade-off between fast, stable folding and tight functional DNA binding. Our approach provides an effective way to quantitatively correlate the roles of protein interactions in conformational dynamics at the multidomain level.

Observing the nonvectorial yet cotranslational folding of a multidomain protein, LDL receptor, in the ER of mammalian cells

Proteins have evolved by incorporating several structural units within a single polypeptide. As a result, multidomain proteins constitute a large fraction of all proteomes. Their domains often fold to their native structures individually and vectorially as each domain emerges from the ribosome or the protein translocation channel, leading to the decreased risk of interdomain misfolding. However, some multidomain proteins fold in the endoplasmic reticulum (ER) nonvectorially via intermediates with nonnative disulfide bonds, which were believed to be shuffled to native ones slowly after synthesis. Yet, the mechanism by which they fold nonvectorially remains unclear. Using two-dimensional (2D) gel electrophoresis and a conformation-specific antibody that recognizes a correctly folded domain, we show here that shuffling of nonnative disulfide bonds to native ones in the most N-terminal region of LDL receptor (LDLR) started at a specific timing during synthesis. Deletion analysis identified a region on LDLR that assisted with disulfide shuffling in the upstream domain, thereby promoting its cotranslational folding. Thus, a plasma membrane-bound multidomain protein has evolved a sequence that promotes the nonvectorial folding of its upstream domains. These findings demonstrate that nonvectorial folding of a multidomain protein in the ER of mammalian cells is more coordinated and elaborated than previously thought. Thus, our findings alter our current view of how a multidomain protein folds nonvectorially in the ER of living cells.

Successes and challenges in simulating the folding of large proteins

Computational simulations of protein folding can be used to interpret experimental folding results, to design new folding experiments, and to test the effects of mutations and small molecules on folding. However, whereas major experimental and computational progress has been made in understanding how small proteins fold, research on larger, multidomain proteins, which comprise the majority of proteins, is less advanced. Specifically, large proteins often fold via long-lived partially folded intermediates, whose structures, potentially toxic oligomerization, and interactions with cellular chaperones remain poorly understood. Molecular dynamics based folding simulations that rely on knowledge of the native structure can provide critical, detailed information on folding free energy landscapes, intermediates, and pathways. Further, increases in computational power and methodological advances have made folding simulations of large proteins practical and valuable. Here, using serpins that inhibit proteases as an example, we review native-centric methods for simulating the folding of large proteins. These synergistic approaches range from Gō and related structure-based models that can predict the effects of the native structure on folding to all-atom-based methods that include side-chain chemistry and can predict how disease-associated mutations may impact folding. The application of these computational approaches to serpins and other large proteins highlights the successes and limitations of current computational methods and underscores how computational results can be used to inform experiments. These powerful simulation approaches in combination with experiments can provide unique insights into how large proteins fold and misfold, expanding our ability to predict and manipulate protein folding.

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.

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.

Molecular simulations of cellular processes

It is, nowadays, possible to simulate biological processes in conditions that mimic the different cellular compartments. Several groups have performed these calculations using molecular models that vary in performance and accuracy. In many cases, the atomistic degrees of freedom have been eliminated, sacrificing both structural complexity and chemical specificity to be able to explore slow processes. In this review, we will discuss the insights gained from computer simulations on macromolecule diffusion, nuclear body formation, and processes involving the genetic material inside cell-mimicking spaces. We will also discuss the challenges to generate new models suitable for the simulations of biological processes on a cell scale and for cell-cycle-long times, including non-equilibrium events such as the co-translational folding, misfolding, and aggregation of proteins. A prominent role will be played by the wise choice of the structural simplifications and, simultaneously, of a relatively complex energetic description. These challenging tasks will rely on the integration of experimental and computational methods, achieved through the application of efficient algorithms. Graphical abstract.

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.

Transiently disordered tails accelerate folding of globular proteins

Numerous biological proteins exhibit intrinsic disorder at their termini, which are associated with multifarious functional roles. Here, we show the surprising result that an increased percentage of terminal short transiently disordered regions with enhanced flexibility (TstDREF) is associated with accelerated folding rates of globular proteins. Evolutionary conservation of predicted disorder at TstDREFs and drastic alteration of folding rates upon point-mutations suggest critical regulatory role(s) of TstDREFs in shaping the folding kinetics. TstDREFs are associated with long-range intramolecular interactions and the percentage of native secondary structural elements physically contacted by TstDREFs exhibit another surprising positive correlation with folding kinetics. These results allow us to infer probable molecular mechanisms behind the TstDREF-mediated regulation of folding kinetics that challenge protein biochemists to assess by direct experimental testing.

Differential bicodon usage in lowly and highly abundant proteins

Degeneracy in the genetic code implies that different codons can encode the same amino acid. Usage preference of synonymous codons has been observed in all domains of life. There is much evidence suggesting that this bias has a major role on protein elongation rate, contributing to differential expression and to co-translational folding. In addition to codon usage bias, other preference variations have been observed such as codon pairs. In this paper, I report that codon pairs have significant different frequency usage for coding either lowly or highly abundant proteins. These usage preferences cannot be explained by the frequency usage of the single codons. The statistical analysis of coding sequences of nine organisms reveals that in many cases bicodon preferences are shared between related organisms. Furthermore, it is observed that misfolding in the drug-transport protein, encoded by MDR1 gene, is better explained by a big change in the pause propensity due to the synonymous bicodon variant, rather than by a relatively small change in codon usage. These findings suggest that codon pair usage can be a more powerful framework to understand translation elongation rate, protein folding efficiency, and to improve protocols to optimize heterologous gene expression.

Protein folding in the cell envelope of Escherichia coli

While the entire proteome is synthesized on cytoplasmic ribosomes, almost half associates with, localizes in or crosses the bacterial cell envelope. In Escherichia coli a variety of mechanisms are important for taking these polypeptides into or across the plasma membrane, maintaining them in soluble form, trafficking them to their correct cell envelope locations and then folding them into the right structures. The fidelity of these processes must be maintained under various environmental conditions including during stress; if this fails, proteases are called in to degrade mislocalized or aggregated proteins. Various soluble, diffusible chaperones (acting as holdases, foldases or pilotins) and folding catalysts are also utilized to restore proteostasis. These responses can be general, dealing with multiple polypeptides, with functional overlaps and operating within redundant networks. Other chaperones are specialized factors, dealing only with a few exported proteins. Several complex machineries have evolved to deal with binding to, integration in and crossing of the outer membrane. This complex protein network is responsible for fundamental cellular processes such as cell wall biogenesis; cell division; the export, uptake and degradation of molecules; and resistance against exogenous toxic factors. The underlying processes, contributing to our fundamental understanding of proteostasis, are a treasure trove for the development of novel antibiotics, biopharmaceuticals and vaccines.