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

Tuning tRNAs for improved translation

Transfer RNAs have been extensively explored as the molecules that translate the genetic code into proteins. At this interface of genetics and biochemistry, tRNAs direct the efficiency of every major step of translation by interacting with a multitude of binding partners. However, due to the variability of tRNA sequences and the abundance of diverse post-transcriptional modifications, a guidebook linking tRNA sequences to specific translational outcomes has yet to be elucidated. Here, we review substantial efforts that have collectively uncovered tRNA engineering principles that can be used as a guide for the tuning of translation fidelity. These principles have allowed for the development of basic research, expansion of the genetic code with non-canonical amino acids, and tRNA therapeutics.

Conformational changes of ribosomes during translation elongation resolved by molecular dynamics simulations

Molecular dynamics simulations have emerged as a powerful set of tools to unravel the intricate dynamics of ribosomes during protein synthesis. Recent advancements in this field have enabled simulations to delve deep into the conformational rearrangements of ribosomes and associated factors, providing invaluable insights into the intricacies of translation. Emphasis on simulations has recently been on translation elongation, such as tRNA selection, translocation, and ribosomal head-swivel motions. These studies have offered crucial structural interpretations of how genetic information is faithfully translated into proteins. This review outlines recent discoveries concerning ribosome conformational changes occurring during translation elongation, as elucidated through molecular dynamics simulations.

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.

tRNA shape is an identity element for an archaeal pyrrolysyl-tRNA synthetase from the human gut

Protein translation is orchestrated through tRNA aminoacylation and ribosomal elongation. Among the highly conserved structure of tRNAs, they have distinguishing features which promote interaction with their cognate aminoacyl tRNA synthetase (aaRS). These key features are referred to as identity elements. In our study, we investigated the tRNA:aaRS pair that installs the 22nd amino acid, pyrrolysine (tRNAPyl:PylRS). Pyrrolysyl-tRNA synthetases (PylRSs) are naturally encoded in some archaeal and bacterial genomes to acylate tRNAPyl with pyrrolysine. Their large amino acid binding pocket and poor recognition of the tRNA anticodon have been instrumental in incorporating >200 noncanonical amino acids. PylRS enzymes can be divided into three classes based on their genomic structure. Two classes contain both an N-terminal and C-terminal domain, however the third class (ΔpylSn) lacks the N-terminal domain. In this study we explored the tRNA identity elements for a ΔpylSn tRNAPyl from Candidatus Methanomethylophilus alvus which drives the orthogonality seen with its cognate PylRS (MaPylRS). From aminoacylation and translation assays we identified five key elements in ΔpylSn tRNAPyl necessary for MaPylRS activity. The absence of a base (position 8) and a G-U wobble pair (G28:U42) were found to affect the high-resolution structure of the tRNA, while molecular dynamic simulations led us to acknowledge the rigidity imparted from the G-C base pairs (G3:C70 and G5:C68).

Geometric alignment of aminoacyl-tRNA relative to catalytic centers of the ribosome underpins accurate mRNA decoding

Accurate protein synthesis is determined by the two-subunit ribosome's capacity to selectively incorporate cognate aminoacyl-tRNA for each mRNA codon. The molecular basis of tRNA selection accuracy, and how fidelity can be affected by antibiotics, remains incompletely understood. Using molecular simulations, we find that cognate and near-cognate tRNAs delivered to the ribosome by Elongation Factor Tu (EF-Tu) can follow divergent pathways of motion into the ribosome during both initial selection and proofreading. Consequently, cognate aa-tRNAs follow pathways aligned with the catalytic GTPase and peptidyltransferase centers of the large subunit, while near-cognate aa-tRNAs follow pathways that are misaligned. These findings suggest that differences in mRNA codon-tRNA anticodon interactions within the small subunit decoding center, where codon-anticodon interactions occur, are geometrically amplified over distance, as a result of this site's physical separation from the large ribosomal subunit catalytic centers. These insights posit that the physical size of both tRNA and ribosome are key determinants of the tRNA selection fidelity mechanism.

Uncovering translation roadblocks during the development of a synthetic tRNA

Ribosomes are remarkable in their malleability to accept diverse aminoacyl-tRNA substrates from both the same organism and other organisms or domains of life. This is a critical feature of the ribosome that allows the use of orthogonal translation systems for genetic code expansion. Optimization of these orthogonal translation systems generally involves focusing on the compatibility of the tRNA, aminoacyl-tRNA synthetase, and a non-canonical amino acid with each other. As we expand the diversity of tRNAs used to include non-canonical structures, the question arises as to the tRNA suitability on the ribosome. Specifically, we investigated the ribosomal translation of allo-tRNAUTu1, a uniquely shaped (9/3) tRNA exploited for site-specific selenocysteine insertion, using single-molecule fluorescence. With this technique we identified ribosomal disassembly occurring from translocation of allo-tRNAUTu1 from the A to the P site. Using cryo-EM to capture the tRNA on the ribosome, we pinpointed a distinct tertiary interaction preventing fluid translocation. Through a single nucleotide mutation, we disrupted this tertiary interaction and relieved the translation roadblock. With the continued diversification of genetic code expansion, our work highlights a targeted approach to optimize translation by distinct tRNAs as they move through the ribosome.

Diffuse Ions Coordinate Dynamics in a Ribonucleoprotein Assembly

Proper ionic concentrations are required for the functional dynamics of RNA and ribonucleoprotein (RNP) assemblies. While experimental and computational techniques have provided many insights into the properties of chelated ions, less is known about the energetic contributions of diffuse ions to large-scale conformational rearrangements. To address this, we present a model that is designed to quantify the influence of diffuse monovalent and divalent ions on the dynamics of biomolecular assemblies. This model employs all-atom (non-H) resolution and explicit ions, where effective potentials account for hydration effects. We first show that the model accurately predicts the number of excess Mg²⁺ ions for prototypical RNA systems, at a level comparable to modern coarse-grained models. We then apply the model to a complete ribosome and show how the balance between diffuse Mg²⁺ and K⁺ ions can control the dynamics of tRNA molecules during translation. The model predicts differential effects of diffuse ions on the free-energy barrier associated with tRNA entry and the energy of tRNA binding to the ribosome. Together, this analysis reveals the direct impact of diffuse ions on the dynamics of an RNP assembly.

The dynamics of subunit rotation in a eukaryotic ribosome

Protein synthesis by the ribosome is coordinated by an intricate series of large-scale conformational rearrangements. Structural studies can provide information about long-lived states, however biological kinetics are controlled by the intervening free-energy barriers. While there has been progress describing the energy landscapes of bacterial ribosomes, very little is known about the energetics of large-scale rearrangements in eukaryotic systems. To address this topic, we constructed an all-atom model with simplified energetics and performed simulations of subunit rotation in the yeast ribosome. In these simulations, the small subunit (SSU; ~1MDa) undergoes spontaneous and reversible rotations (~8°). By enabling the simulation of this rearrangement under equilibrium conditions, these calculations provide initial insights into the molecular factors that control dynamics in eukaryotic ribosomes. Through this, we are able to identify specific inter-subunit interactions that have a pronounced influence on the rate-limiting free-energy barrier. We also show that, as a result of changes in molecular flexibility, the thermodynamic balance between the rotated and unrotated states is temperature-dependent. This effect may be interpreted in terms of differential molecular flexibility within the rotated and unrotated states. Together, these calculations provide a foundation, upon which the field may begin to dissect the energetics of these complex molecular machines.

A steric gate controls P/E hybrid-state formation of tRNA on the ribosome

The ribosome is a biomolecular machine that undergoes multiple large-scale structural rearrangements during protein elongation. Here, we focus on a conformational rearrangement during translocation, known as P/E hybrid-state formation. Using a model that explicitly represents all non-hydrogen atoms, we simulated more than 120 spontaneous transitions, where the tRNA molecule is displaced between the P and E sites of the large subunit. In addition to predicting a free-energy landscape that is consistent with previous experimental observations, the simulations reveal how a six-residue gate-like region can limit P/E formation, where sub-angstrom structural perturbations lead to an order-of-magnitude change in kinetics. Thus, this precisely defined set of residues represents a novel target that may be used to control functional dynamics in bacterial ribosomes. This theoretical analysis establishes a direct relationship between ribosome structure and large-scale dynamics, and it suggests how next-generation experiments may precisely dissect the energetics of hybrid formation on the ribosome.

The ribosome undergoes multiple large-scale structural rearrangements during protein elongation. Here the authors present an all-atom model of the ribosome to study the energetics of P/E hybrid-state formation, an early conformational rearrangement occurring during translocation.

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.

How Nanopore Translocation Experiments Can Measure RNA Unfolding

Electrokinetic translocation of biomolecules through solid-state nanopores represents a label-free single-molecule technique that may be used to measure biomolecular structure and dynamics. Recent investigations have attempted to distinguish individual transfer RNA (tRNA) species based on the associated pore translocation times, ion-current noise, and blockage currents. By manufacturing sufficiently smaller pores, each tRNA is required to undergo a deformation to translocate. Accordingly, differences in nanopore translocation times and distributions may be used to infer the mechanical properties of individual tRNA molecules. To bridge our understanding of tRNA structural dynamics and nanopore measurements, we apply molecular dynamics simulations using a simplified "structure-based" energetic model. Calculating the free-energy landscape for distinct tRNA species implicates transient unfolding of the terminal RNA helix during nanopore translocation. This provides a structural and energetic framework for interpreting current experiments, which can aid the design of methods for identifying macromolecules using nanopores.

The Israeli acute paralysis virus IRES captures host ribosomes by mimicking a ribosomal state with hybrid tRNAs

Colony collapse disorder (CCD) is a multi-faceted syndrome decimating bee populations worldwide, and a group of viruses of the widely distributed Dicistroviridae family have been identified as a causing agent of CCD. This family of viruses employs non-coding RNA sequences, called internal ribosomal entry sites (IRESs), to precisely exploit the host machinery for viral protein production. Using single-particle cryo-electron microscopy (cryo-EM), we have characterized how the IRES of Israeli acute paralysis virus (IAPV) intergenic region captures and redirects translating ribosomes toward viral RNA messages. We reconstituted two in vitro reactions targeting a pre-translocation and a post-translocation state of the IAPV-IRES in the ribosome, allowing us to identify six structures using image processing classification methods. From these, we reconstructed the trajectory of IAPV-IRES from the early small subunit recruitment to the final post-translocated state in the ribosome. An early commitment of IRES/ribosome complexes for global pre-translocation mimicry explains the high efficiency observed for this IRES. Efforts directed toward fighting CCD by targeting the IAPV-IRES using RNA-interference technology are underway, and the structural framework presented here may assist in further refining these approaches.

Diffusion of tRNA inside the ribosome is position-dependent

In recent years, there has been a growing interest to quantify the energy landscape that governs ribosome dynamics. However, in order to quantitatively integrate theoretical predictions and experimental measurements, it is essential that one has a detailed understanding of the associated diffusive properties. Here, all-atom explicit-solvent simulations (50 μs of aggregate sampling) predict that the diffusion coefficient of a tRNA molecule will depend on its position within the ribosome. Specifically, during aa-tRNA accommodation (i.e., the process by which tRNA enters the ribosome), the apparent diffusion coefficient decreases by approximately an order of magnitude. By comparing these to values obtained with an energetically "smooth" model, we show that the observed nonuniform behavior likely arises from electrostatic and solvation interactions between the tRNA and ribosome. These calculations also reveal the hierarchical character of ribosomal energetics, where steric interactions induce a large-scale free-energy barrier, and short-scale roughness determines the rate of diffusive movement across the landscape.

Dissecting the Energetics of Subunit Rotation in the Ribosome

The accurate expression of proteins requires the ribosome to efficiently undergo elaborate conformational rearrangements. The most dramatic of these motions is subunit rotation, which is necessary for tRNA molecules to transition between ribosomal binding sites. While rigid-body descriptions provide a qualitative picture of the process, obtaining quantitative mechanistic insights requires one to account for the relationship between molecular flexibility and collective dynamics. Using simulated rotation events, we assess the quality of experimentally accessible measures for describing the collective displacement of the ∼4000-residue small subunit. For this, we ask whether each coordinate is able to identify the underlying free-energy barrier and transition state ensemble (TSE). We find that intuitive structurally motivated coordinates (e.g., rotation angle, interprotein distances) can distinguish between the endpoints, though they are poor indicators of barrier-crossing events, and they underestimate the free-energy barrier. In contrast, coordinates based on intersubunit bridges can identify the TSE. We additionally verify that the committor probability for the putative TSE configurations is 0.5, a hallmark feature of any transition state. In terms of structural properties, these calculations implicate a transition state in which flexibility allows for asynchronous rearrangements of the bridges, as the ribosome adopts a partially rotated orientation. This provides a theoretical foundation, upon which experimental techniques may precisely quantify the energy landscape of the ribosome.

Ribosome Structure, Function, and Early Evolution

Ribosomes are among the largest and most dynamic molecular motors. The structure and dynamics of translation initiation and elongation are reviewed. Three ribosome motions have been identified for initiation and translocation. A swivel motion between the head/beak and the body of the 30S subunit was observed. A tilting dynamic of the head/beak versus the body of the 30S subunit was detected using simulations. A reversible ratcheting motion was seen between the 30S and the 50S subunits that slide relative to one another. The 30S⁻50S intersubunit contacts regulate translocation. IF2, EF-Tu, and EF-G are homologous G-protein GTPases that cycle on and off the same site on the ribosome. The ribosome, aminoacyl-tRNA synthetase (aaRS) enzymes, transfer ribonucleic acid (tRNA), and messenger ribonucleic acid (mRNA) form the core of information processing in cells and are coevolved. Surprisingly, class I and class II aaRS enzymes, with distinct and incompatible folds, are homologs. Divergence of class I and class II aaRS enzymes and coevolution of the genetic code are described by analysis of ancient archaeal species.

Quantifying the Relationship between Single-Molecule Probes and Subunit Rotation in the Ribosome

A major challenge in the study of biomolecular assemblies is to identify reaction coordinates that precisely monitor conformational rearrangements. This is central to the interpretation of single-molecule fluorescence resonance energy transfer measurements, where the observed dynamics depends on the labeling strategy. As an example, different probes of subunit rotation in the ribosome have provided qualitatively distinct descriptions. In one study, changes in fluorescence suggested that the 30S body undergoes a single rotation/back-rotation cycle during the process of mRNA-tRNA translocation. In contrast, an alternate assay implicated the presence of reversible rotation events before completing translocation. For future single-molecule experiments to unambiguously measure the relationship between subunit rotation and translocation, it is necessary to rationalize these conflicting descriptions. To this end, we have simulated hundreds of spontaneous subunit rotation events (≈8°) using a residue-level coarse-grained model of the ribosome. We analyzed nine different reaction coordinates and found that the apparently inconsistent measurements are likely a consequence of ribosomal flexibility. Further, we propose a metric for quantifying the degree of energetic coupling between experimentally measured degrees of freedom and subunit rotation. This analysis provides a physically grounded framework that can guide the development of more precise single-molecule techniques.

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.