Structure-based energetics of mRNA decoding on the ribosome

Why Does the E1219V Mutation Expand T-Rich PAM Recognition in Cas9 from Streptococcus pyogenes?

Popular RNA-guided DNA endonuclease Cas9 from Streptococcus pyogenes (SpCas9) recognizes the canonical 5'-NGG-3' protospacer adjacent motif (PAM) and triggers double-stranded DNA cleavage activity. Mutations in SpCas9 were demonstrated to expand the PAM readability and hold promise for therapeutic and genome editing applications. However, the energetics of the PAM recognition and its relation to the atomic structure remain unknown. Using the X-ray structure (precatalytic SpCas9:sgRNA:dsDNA) as a template, we calculated the change in the PAM binding affinity in response to SpCas9 mutations using computer simulations. The E1219V mutation in SpCas9 fine-tunes the water accessibility in the PAM binding pocket and promotes new interactions in the SpCas9:noncanonical T-rich PAM, thus weakening the PAM stringency. The nucleotide-specific interaction of two arginine residues (i.e., R1333 and R1335 of SpCas9) ensured stringent 5'-NGG-3' PAM recognition. R1335A substitution (SpCas9R1335A) completely disrupts the direct interaction between SpCas9 and PAM sequences (canonical or noncanonical), accounting for the loss of editing activity. Interestingly, the double mutant (SpCas9R1335A,E1219V) boosts DNA binding affinity by favoring protein:PAM electrostatic contact in a desolvated pocket. The underlying thermodynamics explain the varied DNA cleavage activity of SpCas9 variants. A direct link between the energetics, structures, and activity is highlighted, which can aid in the rational design of improved SpCas9-based genome editing tools.

The energy landscape of the ribosome

The ribosome is a prototypical assembly that can be used to establish general principles and techniques for the study of biological molecular machines. Motivated by the fact that the dynamics of every biomolecule is governed by an underlying energy landscape, there has been great interest to understand and quantify ribosome energetics. In the present review, we will focus on theoretical and computational strategies for probing the interactions that shape the energy landscape of the ribosome, with an emphasis on more recent studies of the elongation cycle. These efforts include the application of quantum mechanical methods for describing chemical kinetics, as well as classical descriptions to characterize slower (microsecond to millisecond) large-scale (10-100 Å) rearrangements, where motion is described in terms of diffusion across an energy landscape. Together, these studies provide broad insights into the factors that control a diverse range of dynamical processes in this assembly.

Importance of residue 248 in Escherichia coli RNase P RNA mediated cleavage

tRNA genes are transcribed as precursors and RNase P generates the matured 5' end of tRNAs. It has been suggested that residue - 1 (the residue immediately 5' of the scissile bond) in the pre-tRNA interacts with the well-conserved bacterial RNase P RNA (RPR) residue A248 (Escherichia coli numbering). The way A248 interacts with residue - 1 is not clear. To gain insight into the role of A248, we analyzed cleavage as a function of A248 substitutions and N-1 nucleobase identity by using pre-tRNA and three model substrates. Our findings are consistent with a model where the structural topology of the active site varies and depends on the identity of the nucleobases at, and in proximity to, the cleavage site and their potential to interact. This leads to positioning of Mg2+ that activates the water that acts as the nucleophile resulting in efficient and correct cleavage. We propose that in addition to be involved in anchoring the substrate the role of A248 is to exclude bulk water from access to the amino acid acceptor stem, thereby preventing non-specific hydrolysis of the pre-tRNA. Finally, base stacking is discussed as a way to protect functionally important base-pairing interactions from non-specific hydrolysis, thereby ensuring high fidelity during RNA processing and the decoding of mRNA.

Machines on Genes through the Computational Microscope

Macromolecular machines acting on genes are at the core of life's fundamental processes, including DNA replication and repair, gene transcription and regulation, chromatin packaging, RNA splicing, and genome editing. Here, we report the increasing role of computational biophysics in characterizing the mechanisms of "machines on genes", focusing on innovative applications of computational methods and their integration with structural and biophysical experiments. We showcase how state-of-the-art computational methods, including classical and ab initio molecular dynamics to enhanced sampling techniques, and coarse-grained approaches are used for understanding and exploring gene machines for real-world applications. As this review unfolds, advanced computational methods describe the biophysical function that is unseen through experimental techniques, accomplishing the power of the "computational microscope", an expression coined by Klaus Schulten to highlight the extraordinary capability of computer simulations. Pushing the frontiers of computational biophysics toward a pragmatic representation of large multimegadalton biomolecular complexes is instrumental in bridging the gap between experimentally obtained macroscopic observables and the molecular principles playing at the microscopic level. This understanding will help harness molecular machines for medical, pharmaceutical, and biotechnological purposes.

Genetic code degeneracy is established by the decoding center of the ribosome

The degeneracy of the genetic code confers a wide array of properties to coding sequences. Yet, its origin is still unclear. A structural analysis has shown that the stability of the Watson–Crick base pair at the second position of the anticodon–codon interaction is a critical parameter controlling the extent of non-specific pairings accepted at the third position by the ribosome, a flexibility at the root of degeneracy. Based on recent cryo-EM analyses, the present work shows that residue A1493 of the decoding center provides a significant contribution to the stability of this base pair, revealing that the ribosome is directly involved in the establishment of degeneracy. Building on existing evolutionary models, we show the evidence that the early appearance of A1493 and A1492 established the basis of degeneracy when an elementary kinetic scheme of translation was prevailing. Logical considerations on the expansion of this kinetic scheme indicate that the acquisition of the peptidyl transferase center was the next major evolutionary step, while the induced-fit mechanism, that enables a sharp selection of the tRNAs, necessarily arose later when G530 was acquired by the decoding center.

Constraints on error rate revealed by computational study of G•U tautomerization in translation

In translation, G•U mismatch in codon-anticodon decoding is an error hotspot likely due to transition of G•U from wobble (wb) to Watson-Crick (WC) geometry, which is governed by keto/enol tautomerization (wb-WC reaction). Yet, effects of the ribosome on the wb-WC reaction and its implications for decoding mechanism remain unclear. Employing quantum-mechanical/molecular-mechanical umbrella sampling simulations using models of the ribosomal decoding site (A site) we determined that the wb-WC reaction is endoergic in the open, but weakly exoergic in the closed A-site state. We extended the classical ‘induced-fit’ model of initial selection by incorporating wb-WC reaction parameters in open and closed states. For predicted parameters, the non-equilibrium exoergic wb-WC reaction is kinetically limited by the decoding rates. The model explains early observations of the WC geometry of G•U from equilibrium structural studies and reveals discrimination capacity for the working ribosome operating at non-equilibrium conditions. The equilibration of the exoergic wb-WC reaction counteracts the equilibration of the open-closed transition of the A site, constraining the decoding accuracy and potentially explaining the persistence of the G•U as an error hotspot. Our results unify structural and mechanistic views of codon-anticodon decoding and generalize the ‘induced-fit’ model for flexible substrates.

Pyrophosphate and Irreversibility in Evolution, or why PPi Is Not an Energy Currency and why Nature Chose Triphosphates

The possible evolutionary significance of pyrophosphate (PPi) has been discussed since the early 1960s. Lipmann suggested that PPi could have been an ancient currency or a possible environmental source of metabolic energy at origins, while Kornberg proposed that PPi vectorializes metabolism because ubiquitous pyrophosphatases render PPi forming reactions kinetically irreversible. To test those ideas, we investigated the reactions that consume phosphoanhydride bonds among the 402 reactions of the universal biosynthetic core that generates amino acids, nucleotides, and cofactors from H2, CO2, and NH3. We find that 36% of the core's phosphoanhydride hydrolyzing reactions generate PPi, while no reactions use PPi as an energy currency. The polymerization reactions that generate ~80% of cell mass - protein, RNA, and DNA synthesis - all generate PPi, while none use PPi as an energy source. In typical prokaryotic cells, aminoacyl tRNA synthetases (AARS) underlie ~80% of PPi production. We show that the irreversibility of the AARS reaction is a kinetic, not a thermodynamic effect. The data indicate that PPi is not an ancient energy currency and probably never was. Instead, PPi hydrolysis is an ancient mechanism that imparts irreversibility, as Kornberg suggested, functioning like a ratchet's pawl to vectorialize the life process toward growth. The two anhydride bonds in nucleoside triphosphates offer ATP-cleaving enzymes an option to impart either thermodynamic control (Pi formation) or kinetic control (PPi formation) upon reactions. This dual capacity explains why nature chose the triphosphate moiety of ATP as biochemistry's universal energy currency.

Thermodynamic control of -1 programmed ribosomal frameshifting

mRNA contexts containing a 'slippery' sequence and a downstream secondary structure element stall the progression of the ribosome along the mRNA and induce its movement into the -1 reading frame. In this study we build a thermodynamic model based on Bayesian statistics to explain how -1 programmed ribosome frameshifting can work. As training sets for the model, we measured frameshifting efficiencies on 64 dnaX mRNA sequence variants in vitro and also used 21 published in vivo efficiencies. With the obtained free-energy difference between mRNA-tRNA base pairs in the 0 and -1 frames, the frameshifting efficiency of a given sequence can be reproduced and predicted from the tRNA-mRNA base pairing in the two frames. Our results further explain how modifications in the tRNA anticodon modulate frameshifting and show how the ribosome tunes the strength of the base-pair interactions.

Principles of tRNAAla Selection by Alanyl-tRNA Synthetase Based on the Critical G3·U70 Base Pair

Throughout evolution, the presence of a single G3·U70 mismatch in the acceptor stem of tRNA^Ala is the major determinant for aminoacylation with alanine by alanyl-tRNA synthetase (AlaRS). Recently reported crystal structures of the complexes AlaRS-tRNA^Ala/G3·U70 and AlaRS-tRNA^Ala/A3·U70 suggest two very different conformations, representing a reactive and a nonreactive state, respectively. On the basis of these structures, it has been proposed that the G3·U70 base pair guides the -CCA end of the tRNA acceptor stem into the active site of AlaRS, thereby enabling aminoacylation. The crystal structures open up the possibility of directly computing the energetics of tRNA specificity by AlaRS. We have carried out molecular dynamics free-energy simulations to quantitatively estimate tRNA discrimination by AlaRS, focusing on the mutations of the single critical base pair G3·U70 to uncover the energetics underlying the accuracy of tRNA selection. The calculations show that the reactive complex is highly selective in favor of the cognate tRNA^Ala/G3·U70 over its noncognate analogues (A3·U70/G3·C70/A3·C70). In contrast, the nonreactive complex is predicted to be unselective between tRNA^Ala/G3·U70 and tRNA^Ala/A3·U70. Utilizing our calculated relative binding free energies, we show how a simple three-step kinetic scheme for aminoacylation, involving both an initial nonspecific binding step and a subsequent transition to a selective reactive complex, accounts for the observed kinetics of the process.

Cryo-EM shows stages of initial codon selection on the ribosome by aa-tRNA in ternary complex with GTP and the GTPase-deficient EF-TuH84A

The GTPase EF-Tu in ternary complex with GTP and aminoacyl-tRNA (aa-tRNA) promotes rapid and accurate delivery of cognate aa-tRNAs to the ribosomal A site. Here we used cryo-EM to study the molecular origins of the accuracy of ribosome-aided recognition of a cognate ternary complex and the accuracy-amplifying role of the monitoring bases A1492, A1493 and G530 of the 16S rRNA. We used the GTPase-deficient EF-Tu variant H84A with native GTP, rather than non-cleavable GTP analogues, to trap a near-cognate ternary complex in high-resolution ribosomal complexes of varying codon-recognition accuracy. We found that ribosome complexes trapped by GTPase-deficicent ternary complex due to the presence of EF-Tu^H84A or non-cleavable GTP analogues have very similar structures. We further discuss speed and accuracy of initial aa-tRNA selection in terms of conformational changes of aa-tRNA and stepwise activation of the monitoring bases at the decoding center of the ribosome.

Accuracy of genetic code translation and its orthogonal corruption by aminoglycosides and Mg2+ ions

We studied the effects of aminoglycosides and changing Mg2+ ion concentration on the accuracy of initial codon selection by aminoacyl-tRNA in ternary complex with elongation factor Tu and GTP (T3) on mRNA programmed ribosomes. Aminoglycosides decrease the accuracy by changing the equilibrium constants of 'monitoring bases' A1492, A1493 and G530 in 16S rRNA in favor of their 'activated' state by large, aminoglycoside-specific factors, which are the same for cognate and near-cognate codons. Increasing Mg2+ concentration decreases the accuracy by slowing dissociation of T3 from its initial codon- and aminoglycoside-independent binding state on the ribosome. The distinct accuracy-corrupting mechanisms for aminoglycosides and Mg2+ ions prompted us to re-interpret previous biochemical experiments and functional implications of existing high resolution ribosome structures. We estimate the upper thermodynamic limit to the accuracy, the 'intrinsic selectivity' of the ribosome. We conclude that aminoglycosides do not alter the intrinsic selectivity but reduce the fraction of it that is expressed as the accuracy of initial selection. We suggest that induced fit increases the accuracy and speed of codon reading at unaltered intrinsic selectivity of the ribosome.

How Messenger RNA and Nascent Chain Sequences Regulate Translation Elongation

Translation elongation is a highly coordinated, multistep, multifactor process that ensures accurate and efficient addition of amino acids to a growing nascent-peptide chain encoded in the sequence of translated messenger RNA (mRNA). Although translation elongation is heavily regulated by external factors, there is clear evidence that mRNA and nascent-peptide sequences control elongation dynamics, determining both the sequence and structure of synthesized proteins. Advances in methods have driven experiments that revealed the basic mechanisms of elongation as well as the mechanisms of regulation by mRNA and nascent-peptide sequences. In this review, we highlight how mRNA and nascent-peptide elements manipulate the translation machinery to alter the dynamics and pathway of elongation.

Substrate-Induced Formation of Ribosomal Decoding Center for Accurate and Rapid Genetic Code Translation

Accurate translation of genetic information is crucial for synthesis of functional proteins in all organisms. We use recent experimental data to discuss how induced fit affects accuracy of initial codon selection on the ribosome by aminoacyl transfer RNA in ternary complex ( T₃) with elongation factor Tu (EF-Tu) and guanosine-5'-triphosphate (GTP). We define actual accuracy ([Formula: see text]) of a particular protein synthesis system as its current accuracy and the effective selectivity ([Formula: see text]) as [Formula: see text] in the limit of zero ribosomal binding affinity for T₃. Intrinsic selectivity ([Formula: see text]), defined as the upper thermodynamic limit of [Formula: see text], is determined by the free energy difference between near-cognate and cognate T₃ in the pre-GTP hydrolysis state on the ribosome. [Formula: see text] is much larger than [Formula: see text], suggesting the possibility of a considerable increase in [Formula: see text] and [Formula: see text] at negligible kinetic cost. Induced fit increases [Formula: see text] and [Formula: see text] without affecting [Formula: see text], and aminoglycoside antibiotics reduce [Formula: see text] and [Formula: see text] at unaltered [Formula: see text].

Roles of specific aminoglycoside-ribosome interactions in the inhibition of translation

Aminoglycosides containing a 2-deoxystreptamine core (AGs) represent a large family of antibiotics that target the ribosome. These compounds promote miscoding, inhibit translocation, and inhibit ribosome recycling. AG binding to helix h44 of the small subunit induces rearrangement of A-site nucleotides A1492 and A1493, which promotes a key open-to-closed conformational change of the subunit and thereby increases miscoding. Mechanisms by which AGs inhibit translocation and recycling remain less clear. Structural studies have revealed a secondary AG binding site in H69 of the large subunit, and it has been proposed that interaction at this site is crucial for inhibition of translocation and recycling. Here, we analyze ribosomes with mutations targeting either or both AG binding sites. Assaying translocation, we find that ablation of the h44 site increases the IC₅₀ values for AGs dramatically, while removal of the H69 site increases these values modestly. This suggests that the AG-h44 interaction is primarily responsible for inhibition, with H69 playing a minor role. Assaying recycling, we find that mutation of h44 has no effect on AG inhibition, consistent with a primary role for AG-H69 interaction. Collectively, these findings help clarify the roles of the two AG binding sites in mechanisms of inhibition by these compounds.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

2'-O-methylation in mRNA disrupts tRNA decoding during translation elongation

Chemical modifications of messenger RNA (mRNA) may regulate many aspects of mRNA processing and protein synthesis. Recently, 2′-O-methylation of nucleotides was identified as a frequent modification in translated regions of human mRNA, showing enrichment in codons for certain amino acid. Here, using single-molecule, bulk kinetics and structural methods, we show that 2′-O-methylation within coding regions of mRNA disrupts key steps in codon reading during cognate transfer RNA (tRNA) selection. Our results suggest that 2′-O-methylation sterically perturbs interactions of ribosomal monitoring bases (G530, A1492 and A1493) with cognate codon-anticodon helices, thereby inhibiting downstream GTP-hydrolysis by elongation factor Tu (EF-Tu) and A-site tRNA accommodation, leading to excessive rejection of cognate aminoacylated-tRNAs in initial selection and proofreading. Our current and prior findings highlight how chemical modifications of mRNA tune the dynamics of protein synthesis at different steps of translation elongation.

Roles of elusive translational GTPases come to light and inform on the process of ribosome biogenesis in bacteria

Protein synthesis relies on several translational GTPases (trGTPases), related proteins that couple the hydrolysis of GTP to specific molecular events on the ribosome. Most bacterial trGTPases, including IF2, EF-Tu, EF-G and RF3, play well-known roles in translation. The cellular functions of LepA (also termed EF4) and BipA (also termed TypA), conversely, have remained enigmatic. Recent studies provide compelling in vivo evidence that LepA and BipA function in biogenesis of the 30S and 50S subunit respectively. These findings have important implications for ribosome biogenesis in bacteria. Because the GTPase activity of each of these proteins depends on interactions with both ribosomal subunits, some portion of 30S and 50S assembly must occur in the context of the 70S ribosome. In this review, we introduce the trGTPases of bacteria, describe the new functional data on LepA and BipA, and discuss the how these findings shape our current view of ribosome biogenesis in bacteria.

Molecular simulations of the ribosome and associated translation factors

The ribosome is a macromolecular complex which is responsible for protein synthesis in all living cells according to their transcribed genetic information. Using X-ray crystallography and, more recently, cryo-electron microscopy (cryo-EM), the structure of the ribosome was resolved at atomic resolution in many functional and conformational states. Molecular dynamics simulations have added information on dynamics and energetics to the available structural information, thereby have bridged the gap to the kinetics obtained from single-molecule and bulk experiments. Here, we review recent computational studies that brought notable insights into ribosomal structure and function.

A recent intermezzo at the Ribosome Club

Two sets of ribosome structures have recently led to two different interpretations of what limits the accuracy of codon translation by transfer RNAs. In this review, inspired by this intermezzo at the Ribosome Club, we briefly discuss accuracy amplification by energy driven proofreading and its implementation in genetic code translation. We further discuss general ways by which the monitoring bases of 16S rRNA may enhance the ultimate accuracy (d-values) and how the codon translation accuracy is reduced by the actions of Mg²⁺ ions and the presence of error inducing aminoglycoside antibiotics. We demonstrate that complete freezing-in of cognate-like tautomeric states of ribosome-bound nucleotide bases in transfer RNA or messenger RNA is not compatible with recent experiments on initial codon selection by transfer RNA in ternary complex with elongation factor Tu and GTP. From these considerations, we suggest that the sets of 30S subunit structures from the Ramakrishnan group and 70S structures from the Yusupov/Yusupova group may, after all, reflect two sides of the same coin and how the structurally based intermezzo at the Ribosome Club may be resolved simply by taking the dynamic aspects of ribosome function into account.

This article is part of the themed issue ‘Perspectives on the ribosome’.

Anisotropic Fluctuations in the Ribosome Determine tRNA Kinetics

The ribosome is a large ribonucleoprotein complex that is responsible for the production of proteins in all organisms. Accommodation is the process by which an incoming aminoacyl-tRNA (aa-tRNA) molecule binds the ribosomal A site, and its kinetics has been implicated in the accuracy of tRNA selection. In addition to rearrangements in the aa-tRNA molecule, the L11 stalk can undergo large-scale anisotropic motions during translation. To explore the potential impact of this protruding region on the rate of aa-tRNA accommodation, we used molecular dynamics simulations with a simplified model to evaluate the free energy as a function of aa-tRNA position. Specifically, these calculations describe the transition between A/T and elbow-accommodated (EA) configurations (∼20 Å displacement). We find that the free-energy barrier associated with elbow accommodation is proportional to the degree of mobility exhibited by the L11 stalk. That is, when L11 is more rigid, the free-energy barrier height is decreased. This effect arises from the ability of L11 to confine, and thereby destabilize, the A/T ensemble. In addition, when elongation factor Tu (EF-Tu) is present, the A/T ensemble is further destabilized in an L11-dependent manner. These results provide a framework that suggests how next-generation experiments may precisely control the dynamics of the ribosome.

Structure-Based Energetics of Stop Codon Recognition by Eukaryotic Release Factor

In translation termination, the eukaryotic release factor (eRF1) recognizes mRNA stop codons (UAA, UAG, or UGA) in a ribosomal A site and triggers release of the nascent polypeptide chain from P-site tRNA. eRF1 is highly selective for U in the first position and a combination of purines (except two consecutive guanines, i.e., GG) in the second and third positions. Eukaryotes decode all three stop codons with a single release factor eRF1, instead of two (RF1 and RF2), in bacteria. Furthermore, unlike bacterial RF1/RF2, eRF1 stabilizes the compact U-turn mRNA configuration in the ribosomal A site by accommodating four nucleotides instead of three. Despite the available cryo-EM structures (resolution ∼3.5-3.8 Å), the energetic principle for eRF1 selectivity toward a stop codon remains a fundamentally unsolved problem. Using cryo-EM structures of eukaryotic translation termination complexes as templates, we carried out molecular dynamics free energy simulations of cognate and near-cognate complexes to quantitatively address the energetics of stop codon recognition by eRF1. Our results suggest that eRF1 has a higher discriminatory power against sense codons, compared to that reported earlier for RF1/RF2. The compact mRNA formed specific intra-mRNA interactions, which itself contributed to stop codon specificity. Furthermore, the specificity is enhanced by the loss of protein-mRNA interactions and, most importantly, by desolvation of the incorrect codons in the near-cognate complexes. Our work provides a clue to how eRF1 discriminates between cognate and near-cognate codons during protein synthesis.

A close-up view of codon selection in eukaryotic initiation

When given an option to choose among a set of alternatives and only one selection is right, one might stop and reflect over which one is best. However, the ribosome has no time to stop and make such reflections, proteins need to be produced and very fast. Eukaryotic translation initiation is an example of such a conundrum. Here, scanning for the correct codon match must be fast, efficient and accurate. We highlight our recent computational findings, which show how the initiation machinery manages to recognize one specific codon among many possible challengers, by fine-tuning the energetic landscape of base-pairing with the aid of the initiation factors eIF1 and eIF1A. Using a recent 3-dimensional structure of the eukaryotic initiation complex we have performed simulations of codon recognition in atomic detail. These calculations provide an in-depth energetic and structural view of how discrimination against near-cognate codons is achieved by the initiation complex.

How the Ribosomal A-Site Finger Can Lead to tRNA Species-Dependent Dynamics

Proteins are synthesized by the joint action of the ribosome and tRNA molecules, where the rate of synthesis can be affected by numerous factors, such as the concentration of tRNA, the binding affinity of tRNA for the ribosome, or post-transcriptional modifications. Here, we expand this range of contributors by demonstrating how differences in tRNA structure can give rise to tRNA species-specific dynamics in the ribosome. To show this, we perform simulations of A/P hybrid-state formation for two tRNA species (tRNA^Phe and tRNA^Leu), which differ in the size of their variable loops (VLs). These calculations reveal that steric interactions between the VL and the ribosomal A-site finger (ASF, i.e., H38 of 23S rRNA) can directly modulate the free-energy landscape for each tRNA species. We also find that tRNA and ASF motions are highly correlated, where fluctuations of the ASF are predictive of tRNA transition events. Finally, by introducing perturbations to the model, we demonstrate that ASF flexibility is a determinant of the rate of A/P hybrid-state formation.

Atomic mutagenesis at the ribosomal decoding site

Ribosomal decoding is an essential process in every living cell. During protein synthesis the 30S ribosomal subunit needs to accomplish binding and accurate decoding of mRNAs. From mutational studies and high-resolution crystal structures nucleotides G530, A1492 and A1493 of the 16S rRNA came into focus as important elements for the decoding process. Recent crystallographic data challenged the so far accepted model for the decoding mechanism. To biochemically investigate decoding in greater detail we applied an in vitro reconstitution approach to modulate single chemical groups at A1492 and A1493. The modified ribosomes were subsequently tested for their ability to efficiently decode the mRNA. Unexpectedly, the ribosome was rather tolerant toward modifications of single groups either at the base or at the sugar moiety in terms of translation activity. Concerning translation fidelity, the elimination of single chemical groups involved in a hydrogen bonding network between the tRNA, mRNA and rRNA did not change the accuracy of the ribosome. These results indicate that the contribution of those chemical groups and the formed hydrogen bonds are not crucial for ribosomal decoding.

How EF-Tu can contribute to efficient proofreading of aa-tRNA by the ribosome

It has long been recognized that the thermodynamics of mRNA–tRNA base pairing is insufficient to explain the high fidelity and efficiency of aminoacyl-tRNA (aa-tRNA) selection by the ribosome. To rationalize this apparent inconsistency, Hopfield proposed that the ribosome may improve accuracy by utilizing a multi-step kinetic proofreading mechanism. While biochemical, structural and single-molecule studies have provided a detailed characterization of aa-tRNA selection, there is a limited understanding of how the physical–chemical properties of the ribosome enable proofreading. To this end, we probe the role of EF-Tu during aa-tRNA accommodation (the proofreading step) through the use of energy landscape principles, molecular dynamics simulations and kinetic models. We find that the steric composition of EF-Tu can reduce the free-energy barrier associated with the first step of accommodation: elbow accommodation. We interpret this effect within an extended kinetic model of accommodation and show how EF-Tu can contribute to efficient and accurate proofreading.

The translation of mRNA by the ribosome is governed by a series of large-scale conformational transitions. Here the authors use MD simulations to demonstrate how the rate of dissociation of elongation factor Tu affects the dynamics of tRNA accommodation and proofreading.

Principles of start codon recognition in eukaryotic translation initiation

Selection of the correct start codon during initiation of translation on the ribosome is a key event in protein synthesis. In eukaryotic initiation, several factors have to function in concert to ensure that the initiator tRNA finds the cognate AUG start codon during mRNA scanning. The two initiation factors eIF1 and eIF1A are known to provide important functions for the initiation process and codon selection. Here, we have used molecular dynamics free energy calculations to evaluate the energetics of initiator tRNA binding to different near-cognate codons on the yeast 40S ribosomal subunit, in the presence and absence of these two initiation factors. The results show that eIF1 and eIF1A together cause a relatively uniform and high discrimination against near-cognate codons. This works such that eIF1 boosts the discrimination against a first position near-cognate G-U mismatch, and also against a second position A-A base pair, while eIF1A mainly acts on third codon position. The computer simulations further reveal the structural basis of the increased discriminatory effect caused by binding of eIF1 and eIF1A to the 40S ribosomal subunit.

Epistasis analysis of 16S rRNA ram mutations helps define the conformational dynamics of the ribosome that influence decoding

The ribosome actively participates in decoding, with a tRNA-dependent rearrangement of the 30S A site playing a key role. Ribosomal ambiguity (ram) mutations have mapped not only to the A site but also to the h12/S4/S5 region and intersubunit bridge B8, implicating other conformational changes such as 30S shoulder rotation and B8 disruption in the mechanism of decoding. Recent crystallographic data have revealed that mutation G299A in helix h12 allosterically promotes B8 disruption, raising the question of whether G299A and/or other ram mutations act mainly via B8. Here, we compared the effects of each of several ram mutations in the absence and presence of mutation h8Δ2, which effectively takes out bridge B8. The data obtained suggest that a subset of mutations including G299A act in part via B8 but predominantly through another mechanism. We also found that G299A in h12 and G347U in h14 each stabilize tRNA in the A site. Collectively, these data support a model in which rearrangement of the 30S A site, inward shoulder rotation, and bridge B8 disruption are loosely coupled events, all of which promote progression along the productive pathway toward peptide bond formation.

Intersubunit Bridges of the Bacterial Ribosome

The ribosome is a large two-subunit ribonucleoprotein machine that translates the genetic code in all cells, synthesizing proteins according to the sequence of the mRNA template. During translation, the primary substrates, transfer RNAs, pass through binding sites formed between the two subunits. Multiple interactions between the ribosomal subunits, termed intersubunit bridges, keep the ribosome intact and at the same time govern dynamics that facilitate the various steps of translation such as transfer RNA-mRNA movement. Here, we review the molecular nature of these intersubunit bridges, how they change conformation during translation, and their functional roles in the process.

Steric interactions lead to collective tilting motion in the ribosome during mRNA-tRNA translocation

Translocation of mRNA and tRNA through the ribosome is associated with large-scale rearrangements of the head domain in the 30S ribosomal subunit. To elucidate the relationship between 30S head dynamics and mRNA–tRNA displacement, we apply molecular dynamics simulations using an all-atom structure-based model. Here we provide a statistical analysis of 250 spontaneous transitions between the A/P–P/E and P/P–E/E ensembles. Consistent with structural studies, the ribosome samples a chimeric ap/P–pe/E intermediate, where the 30S head is rotated ∼18°. It then transiently populates a previously unreported intermediate ensemble, which is characterized by a ∼10° tilt of the head. To identify the origins of head tilting, we analyse 781 additional simulations in which specific steric features are perturbed. These calculations show that head tilting may be attributed to specific steric interactions between tRNA and the 30S subunit (PE loop and protein S13). Taken together, this study demonstrates how molecular structure can give rise to large-scale collective rearrangements.

During protein elongation, the translocation of mRNA and tRNA molecules across the 30S ribosomal subunit is associated with large-scale motions of the 30S head domain. Here the authors carry out MD simulations to probe the associated steric interactions and identify novel tilting motions during the late stages of translocation.

Theoretical Insights on the Mechanism of the GTP Hydrolysis Catalyzed by the Elongation Factor Tu (EF-Tu)

The purpose of this work is to have a better understanding of the mechanism of GTP hydrolysis catalyzed by the elongation factor Tu. Two main aspects are being discussed in the literature: the associative or dissociative character of the process and the nature of nucleophile activation. The calculations of the QM subsystem have been done by means of the M06-2X density functional and the split valence triple-ζ 6-311+G(d,p) basis set. The environmental effect has been introduced through the continuum SMD method. We have studied three models of increasing complexity in order to analyze the different factors that intervene in the catalytic action. The results obtained in this paper confirm that the protonated His84 plays a fundamental role in the catalytic mechanism, but we have also found that the crystallographic sodium ion has a notable effect in the catalysis. So, our work has permitted a new insight, complementary to those obtained with QM/MM calculations, into this very complex process.

Ribosome utilizes the minimum free energy changes to achieve the highest decoding rate and fidelity

The performance of ribosome translation can be characterized by two factors, the translation rate and fidelity. Here, we provide analytical studies of the effect of the near-cognate tRNAs on the two factors. It is shown that the increase of the concentration of the near-cognate tRNAs relative to that of the cognate tRNA has negative effects on the ribosome translation by reducing both the translation rate and the translation fidelity. The effect of the near-cognate ternary complexes on the translation rate results mainly from the initial selection phase, whereas the proofreading phase has a minor effect. By contrast, the effect of the near-cognate ternary complexes on the fidelity results almost equally from the two phases. By using two successive phases, the initial selection and the proofreading, the ribosome can achieve higher translation fidelity than the product of the fidelity when only the initial selection is included and when only the proofreading is included, especially at the large ratio of the concentration of the near-cognate tRNAs compared to that of the cognate tRNA. Moreover, we study the changes of the free energy landscape in the tRNA decoding step. It is found that the rate constants of the tRNA decoding step measured experimentally give the minimum energy changes for the ribosomal complex to attain the optimal performance with both the highest decoding rate and fidelity and/or with the maximum value of the decoding fitness function. This suggests that the ribosome has evolved to utilize the minimum free energy changes gained from the conformational changes of the ribosome, EF-Tu, and tRNA to achieve the optimal performance in the tRNA decoding.

Structural Insights into tRNA Dynamics on the Ribosome

High-resolution structures at different stages, as well as biochemical, single molecule and computational approaches have highlighted the elasticity of tRNA molecules when bound to the ribosome. It is well acknowledged that the inherent structural flexibility of the tRNA lies at the heart of the protein synthesis process. Here, we review the recent advances and describe considerations that the conformational changes of the tRNA molecules offer about the mechanisms grounded in translation.

DNA template dependent accuracy variation of nucleotide selection in transcription

It has been commonly assumed that the effect of erroneous transcription of DNA genes into messenger RNAs on peptide sequence errors are masked by much more frequent errors of mRNA translation to protein. We present a theoretical model of transcriptional accuracy. It uses experimentally estimated standard free energies of double-stranded DNA and RNA/DNA hybrids and predicts a DNA template dependent transcriptional accuracy variation spanning several orders of magnitude. The model also identifies high-error as well a high-accuracy transcription motifs. The source of the large accuracy span is the context dependent variation of the stacking free energy of pairs of correct and incorrect base pairs in the ever moving transcription bubble. Our model predictions have direct experimental support from recent single molecule based identifications of transcriptional errors in the C. elegans transcriptome. Our conclusions challenge the general view that amino acid substitution errors in proteins are mainly caused by translational errors. It suggests instead that transcriptional error hotspots are the dominating source of peptide sequence errors in some DNA template contexts, while mRNA translation is the major cause of protein errors in other contexts.

Modeling the mechanisms of biological GTP hydrolysis

Enzymes that hydrolyze GTP are currently in the spotlight, due to their molecular switch mechanism that controls many cellular processes. One of the best-known classes of these enzymes are small GTPases such as members of the Ras superfamily, which catalyze the hydrolysis of the γ-phosphate bond in GTP. In addition, the availability of an increasing number of crystal structures of translational GTPases such as EF-Tu and EF-G have made it possible to probe the molecular details of GTP hydrolysis on the ribosome. However, despite a wealth of biochemical, structural and computational data, the way in which GTP hydrolysis is activated and regulated is still a controversial topic and well-designed simulations can play an important role in resolving and rationalizing the experimental data. In this review, we discuss the contributions of computational biology to our understanding of GTP hydrolysis on the ribosome and in small GTPases.

Why base tautomerization does not cause errors in mRNA decoding on the ribosome

The structure of the genetic code implies strict Watson–Crick base pairing in the first two codon positions, while the third position is known to be degenerate, thus allowing wobble base pairing. Recent crystal structures of near-cognate tRNAs accommodated into the ribosomal A-site, however, show canonical geometry even with first and second position mismatches. This immediately raises the question of whether these structures correspond to tautomerization of the base pairs. Further, if unusual tautomers are indeed trapped why do they not cause errors in decoding? Here, we use molecular dynamics free energy calculations of ribosomal complexes with cognate and near-cognate tRNAs to analyze the structures and energetics of G-U mismatches in the first two codon positions. We find that the enol tautomer of G is almost isoenergetic with the corresponding ketone in the first position, while it is actually more stable in the second position. Tautomerization of U, on the other hand is highly penalized. The presence of the unusual enol form of G thus explains the crystallographic observations. However, the calculations also show that this tautomer does not cause high codon reading error frequencies, as the resulting tRNA binding free energies are significantly higher than for the cognate complex.