Diffusion of tRNA inside the ribosome is position-dependent

Mitoribosome structure with cofactors and modifications reveals mechanism of ligand binding and interactions with L1 stalk

The mitoribosome translates mitochondrial mRNAs and regulates energy conversion that is a signature of aerobic life forms. We present a 2.2 Å resolution structure of human mitoribosome together with validated mitoribosomal RNA (rRNA) modifications, including aminoacylated CP-tRNAVal. The structure shows how mitoribosomal proteins stabilise binding of mRNA and tRNA helping to align it in the decoding center, whereas the GDP-bound mS29 stabilizes intersubunit communication. Comparison between different states, with respect to tRNA position, allowed us to characterize a non-canonical L1 stalk, and molecular dynamics simulations revealed how it facilitates tRNA transitions in a way that does not require interactions with rRNA. We also report functionally important polyamines that are depleted when cells are subjected to an antibiotic treatment. The structural, biochemical, and computational data illuminate the principal functional components of the translation mechanism in mitochondria and provide a description of the structure and function of the human mitoribosome.

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.

Structural insights into the mechanism of GTP initiation of microtubule assembly

In eukaryotes, the dynamic assembly of microtubules (MT) plays an important role in numerous cellular processes. The underlying mechanism of GTP triggering MT assembly is still unknown. Here, we present cryo-EM structures of tubulin heterodimer at their GTP- and GDP-bound states, intermediate assembly states of GTP-tubulin, and final assembly stages of MT. Both GTP- and GDP-tubulin heterodimers adopt similar curved conformations with subtle flexibility differences. In head-to-tail oligomers of tubulin heterodimers, the inter-dimer interface of GDP-tubulin exhibits greater flexibility, particularly in tangential bending. Cryo-EM of the intermediate assembly states reveals two types of tubulin lateral contacts, "Tube-bond" and "MT-bond". Further, molecular dynamics (MD) simulations show that GTP triggers lateral contact formation in MT assembly in multiple sequential steps, gradually straightening the curved tubulin heterodimers. Therefore, we propose a flexible model of GTP-initiated MT assembly, including the formation of longitudinal and lateral contacts, to explain the nucleation and assembly of MT.

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.

Structure of mitoribosome reveals mechanism of mRNA binding, tRNA interactions with L1 stalk, roles of cofactors and rRNA modifications

The mitoribosome translates mitochondrial mRNAs and regulates energy conversion that is a signature of aerobic life forms. We present a 2.2 Å resolution structure of human mitoribosome together with validated mitoribosomal RNA (rRNA) modifications, including aminoacylated CP-tRNA Val . The structure shows how mitoribosomal proteins stabilise binding of mRNA and tRNA helping to align it in the decoding center, whereas the GDP-bound mS29 stabilizes intersubunit communication. Comparison between different states, with respect to tRNA position, allowed to characterize a non-canonical L1 stalk, and molecular dynamics simulations revealed how it facilitates tRNA transition in a way that does not require interactions with rRNA. We also report functionally important polyamines that are depleted when cells are subjected to an antibiotic treatment. The structural, biochemical, and computational data illuminate the principal functional components of the translation mechanism in mitochondria and provide the most complete description so far of the structure and function of the human mitoribosome.

Coordinate-Dependent Drift-Diffusion Reveals the Kinetic Intermediate Traps of Top7-Based Proteins

The computer-designed Top7 served as a scaffold to produce immunoreactive proteins by grafting of the 2F5 HIV-1 antibody epitope (Top7-2F5) followed by biotinylation (Top7-2F5-biotin). The resulting nonimmunoglobulin affinity proteins were effective in inducing and detecting the HIV-1 antibody. However, the grafted Top7-2F5 design led to protein aggregation, as opposed to the soluble biotinylated Top7-2F5-biotin. The structure-based model predicted that the thermodynamic cooperativity of Top7 increases after grafting and biotin-labeling, reducing their intermediate state populations. In this work, the folding kinetic traps that might contribute to the aggregation propensity are investigated by the diffusion theory. Since the engineered proteins have similar sequence and structural homology, they served as protein models to study the kinetic intermediate traps that were uncovered by characterizing the position-dependent drift-velocity (v(Q)) and the diffusion (D(Q)) coefficients. These coordinate-dependent coefficients were taken into account to obtain the folding and transition path times over the free energy transition states containing the intermediate kinetic traps. This analysis may be useful to predict the aggregated kinetic traps of scaffold-epitope proteins that might compose novel diagnostic and therapeutic platforms.

Precise Steric Features Control Aminoacyl-tRNA Accommodation on the Ribosome

Protein synthesis involves a complex series of large-scale conformational changes in the ribosome. While long-lived intermediate states of these processes can be characterized by experiments, computational methods can be used to identify the interactions that contribute to the rate-limiting free-energy barriers. To this end, we use a simplified energetic model to perform molecular dynamics (MD) simulations of aminoacyl-tRNA (aa-tRNA) accommodation on the ribosome. While numerous studies have probed the energetics of the early stages of accommodation, we focus on the final stage of accommodation, where the 3'-CCA tail of aa-tRNA enters the peptidyl transferase center (PTC). These simulations show how a distinct intermediate is induced by steric confinement of the tail, immediately before it completes accommodation. Multiple pathways for 3'-CCA tail accommodation can be quantitatively distinguished, where the tail enters the PTC by moving past a pocket enclosed by Helix 89, 90, and 92, or through an alternate route formed by Helix 93 and the P-site tRNA. C2573, located within Helix 90, is shown to provide the largest contribution to this late-accommodation steric barrier, such that sub-Å perturbations to this residue can alter the time scale of tail accommodation by nearly an order of magnitude. In terms of biological function, these calculations suggest how this late-stage sterically induced barrier may contribute to tRNA proofreading by the ribosome.

Identifying Strategies to Experimentally Probe Multidimensional Dynamics in the Ribosome

The ribosome is a complex biomolecular machine that utilizes large-scale conformational rearrangements to synthesize proteins. For example, during the elongation cycle, the "head" domain of the ribosomal small subunit (SSU) is known to undergo transient rotation events that allow for movement of tRNA molecules (i.e., translocation). While the head may exhibit rigid-body-like properties, the precise relationship between experimentally accessible probes and multidimensional rotations has yet to be established. To address this gap, we perform molecular dynamics simulations of the translocation step of the elongation cycle in the ribosome, where the SSU head spontaneously undergoes rotation and tilt-like motions. With this data set (1250 simulated events), we used statistical and information-theory-based measures to identify possible single-molecule probes that can isolate SSU head rotation and head tilting. This analysis provides a molecular interpretation for previous single-molecule measurements, while establishing a framework for the design of next-generation experiments that may precisely probe the mechanistic and kinetic aspects of the ribosome.

Biotin-painted proteins have thermodynamic stability switched by kinetic folding routes

Biotin-labeled proteins are widely used as tools to study protein-protein interactions and proximity in living cells. Proteomic methods broadly employ proximity-labeling technologies based on protein biotinylation in order to investigate the transient encounters of biomolecules in subcellular compartments. Biotinylation is a post-translation modification in which the biotin molecule is attached to lysine or tyrosine residues. So far, biotin-based technologies proved to be effective instruments as affinity and proximity tags. However, the influence of biotinylation on aspects such as folding, binding, mobility, thermodynamic stability, and kinetics needs to be investigated. Here, we selected two proteins [biotin carboxyl carrier protein (BCCP) and FKBP3] to test the influence of biotinylation on thermodynamic and kinetic properties. Apo (without biotin) and holo (biotinylated) protein structures were used separately to generate all-atom structure-based model simulations in a wide range of temperatures. Holo BCCP contains one biotinylation site, and FKBP3 was modeled with up to 23 biotinylated lysines. The two proteins had their estimated thermodynamic stability changed by altering their energy landscape. In all cases, after comparison between the apo and holo simulations, differences were observed on the free-energy profiles and folding routes. Energetic barriers were altered with the density of states clearly showing changes in the transition state. This study suggests that analysis of large-scale datasets of biotinylation-based proximity experiments might consider possible alterations in thermostability and folding mechanisms imposed by the attached biotins.

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.

SMOG 2 and OpenSMOG: Extending the limits of structure-based models

Applying simulations with structure-based G o ¯ - like models has proven to be an effective strategy for investigating the factors that control biomolecular dynamics. The common element of these models is that some (or all) of the intra/inter-molecular interactions are explicitly defined to stabilize an experimentally determined structure. To facilitate the development and application of this broad class of models, we previously released the SMOG 2 software package. This suite allows one to easily customize and distribute structure-based (i.e., SMOG) models for any type of polymer-ligand system. The force fields generated by SMOG 2 may then be used to perform simulations in highly optimized MD packages, such as Gromacs, NAMD, LAMMPS, and OpenMM. Here, we describe extensions to the software and demonstrate the capabilities of the most recent version (SMOG v2.4.2). Changes include new tools that aid user-defined customization of force fields, as well as an interface with the OpenMM simulation libraries (OpenSMOG v1.1.0). The OpenSMOG module allows for arbitrary user-defined contact potentials and non-bonded potentials to be employed in SMOG models, without source-code modifications. To illustrate the utility of these advances, we present applications to systems with millions of atoms, long polymers and explicit ions, as well as models that include non-structure-based (e.g., AMBER-based) energetic terms. Examples include large-scale rearrangements of the SARS-CoV-2 Spike protein, the HIV-1 capsid with explicit ions, and crystallographic lattices of ribosomes and proteins. In summary, SMOG 2 and OpenSMOG provide robust support for researchers who seek to develop and apply structure-based models to large and/or intricate biomolecular systems.

Sterically confined rearrangements of SARS-CoV-2 Spike protein control cell invasion

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is highly contagious, and transmission involves a series of processes that may be targeted by vaccines and therapeutics. During transmission, host cell invasion is controlled by a large-scale (200–300 Å) conformational change of the Spike protein. This conformational rearrangement leads to membrane fusion, which creates transmembrane pores through which the viral genome is passed to the host. During Spike-protein-mediated fusion, the fusion peptides must be released from the core of the protein and associate with the host membrane. While infection relies on this transition between the prefusion and postfusion conformations, there has yet to be a biophysical characterization reported for this rearrangement. That is, structures are available for the endpoints, though the intermediate conformational processes have not been described. Interestingly, the Spike protein possesses many post-translational modifications, in the form of branched glycans that flank the surface of the assembly. With the current lack of data on the pre-to-post transition, the precise role of glycans during cell invasion has also remained unclear. To provide an initial mechanistic description of the pre-to-post rearrangement, an all-atom model with simplified energetics was used to perform thousands of simulations in which the protein transitions between the prefusion and postfusion conformations. These simulations indicate that the steric composition of the glycans can induce a pause during the Spike protein conformational change. We additionally show that this glycan-induced delay provides a critical opportunity for the fusion peptides to capture the host cell. In contrast, in the absence of glycans, the viral particle would likely fail to enter the host. This analysis reveals how the glycosylation state can regulate infectivity, while providing a much-needed structural framework for studying the dynamics of this pervasive pathogen.

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.

RNA Pore Translocation with Static and Periodic Forces: Effect of Secondary and Tertiary Elements on Process Activation and Duration

We use MD simulations to study the pore translocation properties of a pseudoknotted viral RNA. We consider the 71-nucleotide-long xrRNA from the Zika virus and establish how it responds when driven through a narrow pore by static or periodic forces applied to either of the two termini. Unlike the case of fluctuating homopolymers, the onset of translocation is significantly delayed with respect to the application of static driving forces. Because of the peculiar xrRNA architecture, activation times can differ by orders of magnitude at the two ends. Instead, translocation duration is much smaller than activation times and occurs on time scales comparable at the two ends. Periodic forces amplify significantly the differences at the two ends, for both activation times and translocation duration. Finally, we use a waiting-times analysis to examine the systematic slowing downs in xrRNA translocations and associate them to the hindrance of specific secondary and tertiary elements of xrRNA. The findings provide a useful reference to interpret and design future theoretical and experimental studies of RNA translocation.

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.

Simulations of Phage T7 Capsid Expansion Reveal the Role of Molecular Sterics on Dynamics

Molecular dynamics techniques provide numerous strategies for investigating biomolecular energetics, though quantitative analysis is often only accessible for relatively small (frequently monomeric) systems. To address this limit, we use simulations in combination with a simplified energetic model to study complex rearrangements in a large assembly. We use cryo-EM reconstructions to simulate the DNA packaging-associated 3 nm expansion of the protein shell of an initially assembled phage T7 capsid (called procapsid or capsid I). This is accompanied by a disorder-order transition and expansion-associated externalization displacement of the 420 N-terminal tails of the shell proteins. For the simulations, we use an all-atom structure-based model (1.07 million atoms), which is specifically designed to probe the influence of molecular sterics on dynamics. We find that the rate at which the N-terminal tails undergo translocation depends heavily on their position within hexons and pentons. Specifically, trans-shell displacements of the hexon E subunits are the most frequent and hexon A subunits are the least frequent. The simulations also implicate numerous tail translocation intermediates during tail translocation that involve topological traps, as well as sterically induced barriers. The presented study establishes a foundation for understanding the precise relationship between molecular structure and phage maturation.

Quantitative comparison between sub-millisecond time resolution single-molecule FRET measurements and 10-second molecular simulations of a biosensor protein

Molecular Dynamics (MD) simulations seek to provide atomic-level insights into conformationally dynamic biological systems at experimentally relevant time resolutions, such as those afforded by single-molecule fluorescence measurements. However, limitations in the time scales of MD simulations and the time resolution of single-molecule measurements have challenged efforts to obtain overlapping temporal regimes required for close quantitative comparisons. Achieving such overlap has the potential to provide novel theories, hypotheses, and interpretations that can inform idealized experimental designs that maximize the detection of the desired reaction coordinate. Here, we report MD simulations at time scales overlapping with in vitro single-molecule Förster (fluorescence) resonance energy transfer (smFRET) measurements of the amino acid binding protein LIV-BP^SS at sub-millisecond resolution. Computationally efficient all-atom structure-based simulations, calibrated against explicit solvent simulations, were employed for sampling multiple cycles of LIV-BP^SS clamshell-like conformational changes on the time scale of seconds, examining the relationship between these events and those observed by smFRET. The MD simulations agree with the smFRET measurements and provide valuable information on local dynamics of fluorophores at their sites of attachment on LIV-BP^SS and the correlations between fluorophore motions and large-scale conformational changes between LIV-BP^SS domains. We further utilize the MD simulations to inform the interpretation of smFRET data, including Förster radius (R₀) and fluorophore orientation factor (κ²) determinations. The approach we describe can be readily extended to distinct biochemical systems, allowing for the interpretation of any FRET system conjugated to protein or ribonucleoprotein complexes, including those with more conformational processes, as well as those implementing multi-color smFRET.

Förster (fluorescence) resonance energy transfer (FRET) has been used extensively by biophysicists as a molecular-scale ruler that yields fundamental structural and kinetic insights into transient processes including complex formation and conformational rearrangements required for biological function. FRET techniques require the identification of informative fluorophore labeling sites, spaced at defined distances to inform on a reaction coordinate of interest and consideration of noise sources that have the potential to obscure quantitative interpretations. Here, we describe an approach to leverage advancements in computationally efficient all-atom structure-based molecular dynamics simulations in which structural dynamics observed via FRET can be interpreted in full atomistic detail on commensurate time scales. We demonstrate the potential of this approach using a model FRET system, the amino acid binding protein LIV-BP^SS conjugated to self-healing organic fluorophores. LIV-BP^SS exhibits large scale, sub-millisecond clamshell-like conformational changes between open and closed conformations associated with ligand unbinding and binding, respectively. Our findings inform on the molecular basis of the dynamics observed by smFRET and on strategies to optimize fluorophore labeling sites, the manner of fluorophore attachment, and fluorophore composition.

Particle Filter Method to Integrate High-Speed Atomic Force Microscopy Measurements with Biomolecular Simulations

High-speed atomic force microscopy (HS-AFM) can be used to observe the structural dynamics of biomolecules at the single-molecule level in real time under near-physiological conditions; however, its spatiotemporal resolution is limited. Complementarily, molecular dynamics (MD) simulations have higher spatiotemporal resolutions, albeit with some artifacts. Here, to integrate HS-AFM data and coarse-grained molecular dynamics (CG-MD) simulations, we develop a particle filter method that implements a sequential Bayesian data assimilation approach. We test the method in a twin experiment. First, we generate a reference HS-AFM movie from the CG-MD trajectory of a test molecule, a nucleosome; this serves as the "experimental measurement". Then, we perform a particle filter simulation with 512 particles, which captures the large-scale nucleosome structural dynamics compatible with the AFM movie. Comparing particle filter simulations with 8-8192 particles, we find that using greater numbers of particles consistently increases the likelihood of the whole AFM movie. By comparing the likelihoods for different ionic concentrations and time scale mappings, we find that the "true" concentration and time scale mapping can be inferred as the largest likelihood of the whole AFM movie but not that of each AFM image. The particle filter method provides a general approach for integrating HS-AFM data with MD simulations.

Frustration and inhomogeneous environments in relaxation of open chains with Ising-type interactions

Frustration can contribute to very slow relaxation times in large open chains, as in spin glasses and in biopolymers. However, frustration may not be sufficient to produce broken ergodicity in finite systems. Here we employ a system-plus-reservoir approach to investigate how strongly inhomogeneous environments and frustration compete in the relaxation of finite open chains. We find a sufficient condition for our inhomogeneous environments to break ergodicity. We use the microscopic model to derive a Markovian quantum master equation for a generic chain with ultrastrong intrachain couplings. We show that this microscopic model avoids a spurious broken ergodicity we find in the phenomenological model. We work out an explicit example of broken ergodicity due to the inhomogeneous environment of an unfrustrated spin chain as far as simulating a recent experiment on protein denaturation (where environment inhomogeneity is especially relevant). We finally show that an inhomogeneous environment can mitigate the effects of frustration-induced degeneracies.

Directional translocation resistance of Zika xrRNA

xrRNAs from flaviviruses survive in host cells because of their exceptional dichotomic response to the unfolding action of different enzymes. They can be unwound, and hence copied, by replicases, and yet can resist degradation by exonucleases. How the same stretch of xrRNA can encode such diverse responses is an open question. Here, by using atomistic models and translocation simulations, we uncover an elaborate and directional mechanism for how stress propagates when the two xrRNA ends, [Formula: see text] and [Formula: see text], are driven through a pore. Pulling the [Formula: see text] end, as done by replicases, elicits a progressive unfolding; pulling the [Formula: see text] end, as done by exonucleases, triggers a counterintuitive molecular tightening. Thus, in what appears to be a remarkable instance of intra-molecular tensegrity, the very pulling of the [Formula: see text] end is what boosts resistance to translocation and consequently to degradation. The uncovered mechanistic principle might be co-opted to design molecular meta-materials.

Molecular dynamics simulations of DNA-DNA and DNA-protein interactions

The all-atom molecular dynamics method can characterize the molecular-level interactions in DNA and DNA-protein systems with unprecedented resolution. Recent advances in computational technologies have allowed the method to reveal the unbiased behavior of such systems at the microseconds time scale, whereas enhanced sampling approaches have matured enough to characterize the interaction free energy with quantitative precision. Here, we describe recent progress toward increasing the realism of such simulations by refining the accuracy of the molecular dynamics force field, and we highlight recent application of the method to systems of outstanding biological interest.

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.