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

Synonymous Mutations Can Alter Protein Dimerization Through Localized Interface Misfolding Involving Self-entanglements

Synonymous mutations in messenger RNAs (mRNAs) can reduce protein-protein binding substantially without changing the protein's amino acid sequence. Here, we use coarse-grain simulations of protein synthesis, post-translational dynamics, and dimerization to understand how synonymous mutations can influence the dimerization of two E. coli homodimers, oligoribonuclease and ribonuclease T. We synthesize each protein from its wildtype, fastest- and slowest-translating synonymous mRNAs in silico and calculate the ensemble-averaged interaction energy between the resulting dimers. We find synonymous mutations alter oligoribonuclease's dimer properties. Relative to wildtype, the dimer interaction energy becomes 4% and 10% stronger, respectively, when translated from its fastest- and slowest-translating mRNAs. Ribonuclease T dimerization, however, is insensitive to synonymous mutations. The structural and kinetic origin of these changes are misfolded states containing non-covalent lasso-entanglements, many of which structurally perturb the dimer interface, and whose probability of occurrence depends on translation speed. These entangled states are kinetic traps that persist for long time scales. Entanglements cause altered dimerization energies for oligoribonuclease, as there is a large association (odds ratio: 52) between the co-occurrence of non-native self-entanglements and weak-binding dimer conformations. Simulated at all-atom resolution, these entangled structures persist for long timescales, indicating the conclusions are independent of model resolution. Finally, we show that regions of the protein we predict to have changes in entanglement are also structurally perturbed during refolding, as detected by limited-proteolysis mass spectrometry. Thus, non-native changes in entanglement at dimer interfaces is a mechanism through which oligomer structure and stability can be altered.

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

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

Thermodynamic Cost, Speed, Fluctuations, and Error Reduction of Biological Copy Machines

Due to large fluctuations in cellular environments, transfer of information in biological processes without regulation is error-prone. The mechanistic details of error-reducing mechanisms in biological copying processes have been a subject of active research; however, how error reduction of a process is balanced with its thermodynamic cost and dynamical properties remain largely unexplored. Here, we study the error reducing strategies in light of the recently discovered thermodynamic uncertainty relation (TUR) that sets a physical bound to the cost-precision trade-off for dissipative processes. We found that the two representative copying processes, DNA replication by the exonuclease-deficient T7 DNA polymerase and mRNA translation by the E. coli ribosome, reduce the error rates to biologically acceptable levels while also optimizing the processes close to the physical limit dictated by TUR.

Elongation Factor Tu Switch I Element is a Gate for Aminoacyl-tRNA Selection

Selection of correct aminoacyl (aa)-tRNA at the ribosomal A site is fundamental to maintaining translational fidelity. Aa-tRNA selection is a multistep process facilitated by the guanosine triphosphatase elongation factor (EF)-Tu. EF-Tu delivers aa-tRNA to the ribosomal A site and participates in tRNA selection. The structural mechanism of how EF-Tu is involved in proofreading remains to be fully resolved. Here, we provide evidence that switch I of EF-Tu facilitates EF-Tu's involvement during aa-tRNA selection. Using structure-based and explicit solvent molecular dynamics simulations based on recent cryo-electron microscopy reconstructions, we studied the conformational change of EF-Tu from the guanosine triphosphate to guanine diphosphate conformation during aa-tRNA accommodation. Switch I of EF-Tu rapidly converts from an α-helix into a β-hairpin and moves to interact with the acceptor stem of the aa-tRNA. In doing so, switch I gates the movement of the aa-tRNA during accommodation through steric interactions with the acceptor stem. Pharmacological inhibition of the aa-tRNA accommodation pathway prevents the proper positioning of switch I with the aa-tRNA acceptor stem, suggesting that the observed interactions are specific for cognate aa-tRNA substrates, and thus capable of contributing to the fidelity mechanism.

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

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

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

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

Translation in Prokaryotes

This review summarizes our current understanding of translation in prokaryotes, focusing on the mechanistic and structural aspects of each phase of translation: initiation, elongation, termination, and ribosome recycling. The assembly of the initiation complex provides multiple checkpoints for messenger RNA (mRNA) and start-site selection. Correct codon-anticodon interaction during the decoding phase of elongation results in major conformational changes of the small ribosomal subunit and shapes the reaction pathway of guanosine triphosphate (GTP) hydrolysis. The ribosome orchestrates proton transfer during peptide bond formation, but requires the help of elongation factor P (EF-P) when two or more consecutive Pro residues are to be incorporated. Understanding the choreography of transfer RNA (tRNA) and mRNA movements during translocation helps to place the available structures of translocation intermediates onto the time axis of the reaction pathway. The nascent protein begins to fold cotranslationally, in the constrained space of the polypeptide exit tunnel of the ribosome. When a stop codon is reached at the end of the coding sequence, the ribosome, assisted by termination factors, hydrolyzes the ester bond of the peptidyl-tRNA, thereby releasing the nascent protein. Following termination, the ribosome is dissociated into subunits and recycled into another round of initiation. At each step of translation, the ribosome undergoes dynamic fluctuations between different conformation states. The aim of this article is to show the link between ribosome structure, dynamics, and function.

Timing during translation matters: synonymous mutations in human pathologies influence protein folding and function

Ribosomes translate mRNAs with non-uniform speed. Translation velocity patterns are a conserved feature of mRNA and have evolved to fine-tune protein folding, expression and function. Synonymous single-nucleotide polymorphisms (sSNPs) that alter programmed translational speed affect expression and function of the encoded protein. Synergistic advances in next-generation sequencing have led to the identification of sSNPs associated with disease penetrance. Here, we draw on studies with disease-related proteins to enhance our understanding of mechanistic contributions of sSNPs to functional alterations of the encoded protein. We emphasize the importance of identification of sSNPs along with disease-causing mutations to understand genotype-phenotype relationships.

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.

Protein Sequences Recapitulate Genetic Code Evolution

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

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.

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.

Evidence of evolutionary selection for cotranslational folding

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

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

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

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

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