Accommodation of aminoacyl-tRNA into the ribosome involves reversible excursions along multiple pathways

Supercomputing in the biological sciences: Toward Zettascale and Yottascale simulations

Molecular simulations of biological systems tend to be significantly more compute-intensive than those in materials science and astrophysics, due to important contributions of long-range electrostatic forces and large numbers of time steps (>1E9) required. Simulations of biomolecular complexes of microseconds to milliseconds are considered state-of-the-art today. However, these time scales are miniscule in comparison to physiological time scales relevant to molecular machine activity, drug action, and elongation cycles for protein synthesis, RNA synthesis, and DNA synthesis (seconds to days). While an exascale supercomputer has simulated an entire virus for nanoseconds, this supercomputer would need to be 10 billion times faster to simulate that virus for 3 hours of physiological time, demonstrating the insatiable need for computing power. With growing interest in computational drug design from the pharmaceutical sector, the biological sciences are positioned to be an industry driver in computing.

Temperature-Dependent tRNA Modifications in Bacillales

Transfer RNA (tRNA) modifications are essential for the temperature adaptation of thermophilic and psychrophilic organisms as they control the rigidity and flexibility of transcripts. To further understand how specific tRNA modifications are adjusted to maintain functionality in response to temperature fluctuations, we investigated whether tRNA modifications represent an adaptation of bacteria to different growth temperatures (minimal, optimal, and maximal), focusing on closely related psychrophilic (P. halocryophilus and E. sibiricum), mesophilic (B. subtilis), and thermophilic (G. stearothermophilus) Bacillales. Utilizing an RNA sequencing approach combined with chemical pre-treatment of tRNA samples, we systematically profiled dihydrouridine (D), 4-thiouridine (s4U), 7-methyl-guanosine (m7G), and pseudouridine (Ψ) modifications at single-nucleotide resolution. Despite their close relationship, each bacterium exhibited a unique tRNA modification profile. Our findings revealed increased tRNA modifications in the thermophilic bacterium at its optimal growth temperature, particularly showing elevated levels of s4U8 and Ψ55 modifications compared to non-thermophilic bacteria, indicating a temperature-dependent regulation that may contribute to thermotolerance. Furthermore, we observed higher levels of D modifications in psychrophilic and mesophilic bacteria, indicating an adaptive strategy for cold environments by enhancing local flexibility in tRNAs. Our method demonstrated high effectiveness in identifying tRNA modifications compared to an established tool, highlighting its potential for precise tRNA profiling studies.

Alternate conformational trajectories in ribosome translocation

Translocation in protein synthesis entails the efficient and accurate movement of the mRNA-[tRNA]2 substrate through the ribosome after peptide bond formation. An essential conformational change during this process is the swiveling of the small subunit head domain about two rRNA 'hinge' elements. Using iterative selection and molecular dynamics simulations, we derive alternate hinge elements capable of translocation in vitro and in vivo and describe their effects on the conformational trajectory of the EF-G-bound, translocating ribosome. In these alternate conformational pathways, we observe a diversity of swivel kinetics, hinge motions, three-dimensional head domain trajectories and tRNA dynamics. By finding alternate conformational pathways of translocation, we identify motions and intermediates that are essential or malleable in this process. These findings highlight the plasticity of protein synthesis and provide a more thorough understanding of the available sequence and conformational landscape of a central biological process.

Single-Molecule Imaging of Integral Membrane Protein Dynamics and Function

Integral membrane proteins (IMPs) play central roles in cellular physiology and represent the majority of known drug targets. Single-molecule fluorescence and fluorescence resonance energy transfer (FRET) methods have recently emerged as valuable tools for investigating structure-function relationships in IMPs. This review focuses on the practical foundations required for examining polytopic IMP function using single-molecule FRET (smFRET) and provides an overview of the technical and conceptual frameworks emerging from this area of investigation. In this context, we highlight the utility of smFRET methods to reveal transient conformational states critical to IMP function and the use of smFRET data to guide structural and drug mechanism-of-action investigations. We also identify frontiers where progress is likely to be paramount to advancing the field.

Engineering Ribosomal Machinery for Noncanonical Amino Acid Incorporation

The introduction of noncanonical amino acids into proteins has enabled researchers to modify fundamental physicochemical and functional properties of proteins. While the alteration of the genetic code, via the introduction of orthogonal aminoacyl-tRNA synthetase:tRNA pairs, has driven many of these efforts, the various components involved in the process of translation are important for the development of new genetic codes. In this review, we will focus on recent advances in engineering ribosomal machinery for noncanonical amino acid incorporation and genetic code modification. The engineering of the ribosome itself will be considered, as well as the many factors that interact closely with the ribosome, including both tRNAs and accessory factors, such as the all-important EF-Tu. Given the success of genome re-engineering efforts, future paths for radical alterations of the genetic code will require more expansive alterations in the translation machinery.

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.

Beyond the Anticodon: tRNA Core Modifications and Their Impact on Structure, Translation and Stress Adaptation

Transfer RNAs (tRNAs) are heavily decorated with post-transcriptional chemical modifications. Approximately 100 different modifications have been identified in tRNAs, and each tRNA typically contains 5-15 modifications that are incorporated at specific sites along the tRNA sequence. These modifications may be classified into two groups according to their position in the three-dimensional tRNA structure, i.e., modifications in the tRNA core and modifications in the anticodon-loop (ACL) region. Since many modified nucleotides in the tRNA core are involved in the formation of tertiary interactions implicated in tRNA folding, these modifications are key to tRNA stability and resistance to RNA decay pathways. In comparison to the extensively studied ACL modifications, tRNA core modifications have generally received less attention, although they have been shown to play important roles beyond tRNA stability. Here, we review and place in perspective selected data on tRNA core modifications. We present their impact on tRNA structure and stability and report how these changes manifest themselves at the functional level in translation, fitness and stress adaptation.

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.

RlmQ: a newly discovered rRNA modification enzyme bridging RNA modification and virulence traits in Staphylococcus aureus

rRNA modifications play crucial roles in fine-tuning the delicate balance between translation speed and accuracy, yet the underlying mechanisms remain elusive. Comparative analyses of the rRNA modifications in taxonomically distant bacteria could help define their general, as well as species-specific, roles. In this study, we identified a new methyltransferase, RlmQ, in Staphylococcus aureus responsible for the Gram-positive specific m7G2601, which is not modified in Escherichia coli (G2574). We also demonstrate the absence of methylation on C1989, equivalent to E. coli C1962, which is methylated at position 5 by the Gram-negative specific RlmI methyltransferase, a paralog of RlmQ. Both modifications (S. aureus m7G2601 and E. coli m5C1962) are situated within the same tRNA accommodation corridor, hinting at a potential shared function in translation. Inactivation of S. aureus rlmQ causes the loss of methylation at G2601 and significantly impacts growth, cytotoxicity, and biofilm formation. These findings unravel the intricate connections between rRNA modifications, translation, and virulence in pathogenic Gram-positive bacteria.

Dynamics and Function of sRNA/mRNAs Under the Scrutiny of Computational Simulation Methods

Molecular dynamics simulations have proved extremely useful in investigating the functioning of proteins with atomic-scale resolution. Many applications to the study of RNA also exist, and their number increases by the day. However, implementing MD simulations for RNA molecules in solution faces challenges that the MD practitioner must be aware of for the appropriate use of this tool. In this chapter, we present the fundamentals of MD simulations, in general, and the peculiarities of RNA simulations, in particular. We discuss the strengths and limitations of the technique and provide examples of its application to elucidate small RNA's performance.

Ribosome Abundance Control in Prokaryotes

Cell growth is an essential phenotype of any unicellular organism and it crucially depends on precise control of protein synthesis. We construct a model of the feedback mechanisms that regulate abundance of ribosomes in E. coli, a prototypical prokaryotic organism. Since ribosomes are needed to produce more ribosomes, the model includes a positive feedback loop central to the control of cell growth. Our analysis of the model shows that there can be only two coexisting equilibrium states across all 23 parameters. This precludes the existence of hysteresis, suggesting that the ribosome abundance changes continuously with parameters. These states are related by a transcritical bifurcation, and we provide an analytic formula for parameters that admit either state.

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.

Characterizing the Folding Transition-State Ensembles in the Energy Landscape of an RNA Tetraloop

Molecular dynamics (MD) simulations have become increasingly powerful and can now describe the folding/unfolding of small biomolecules in atomic detail. However, a major challenge in MD simulations is to represent the complex energy landscape of biomolecules using a small number of reaction coordinates. In this study, we investigate the folding pathways of an RNA tetraloop, gcGCAAgc, using five classical MD simulations with a combined simulation time of approximately 120 μs. Our approach involves analyzing the tetraloop dynamics, including the folding transition state ensembles, using the energy landscape visualization method (ELViM). The ELViM is an approach that uses internal distances to compare any two conformations, allowing for a detailed description of the folding process without requiring root mean square alignment of structures. This method has previously been applied to describe the energy landscape of disordered β-amyloid peptides and other proteins. The ELViM results in a non-linear projection of the multidimensional space, providing a comprehensive representation of the tetraloop's energy landscape. Our results reveal four distinct transition-state regions and establish the paths that lead to the folded tetraloop structure. This detailed analysis of the tetraloop's folding process has important implications for understanding RNA folding, and the ELViM approach can be used to study other biomolecules.

Decoding and Recoding of mRNA Sequences by the Ribosome

Faithful translation of messenger RNA (mRNA) into protein is essential to maintain protein homeostasis in the cell. Spontaneous translation errors are very rare due to stringent selection of cognate aminoacyl transfer RNAs (tRNAs) and the tight control of the mRNA reading frame by the ribosome. Recoding events, such as stop codon readthrough, frameshifting, and translational bypassing, reprogram the ribosome to make intentional mistakes and produce alternative proteins from the same mRNA. The hallmark of recoding is the change of ribosome dynamics. The signals for recoding are built into the mRNA, but their reading depends on the genetic makeup of the cell, resulting in cell-specific changes in expression programs. In this review, I discuss the mechanisms of canonical decoding and tRNA-mRNA translocation; describe alternative pathways leading to recoding; and identify the links among mRNA signals, ribosome dynamics, and recoding.

mRNA decoding in human is kinetically and structurally distinct from bacteria

In all species, ribosomes synthesize proteins by faithfully decoding messenger RNA (mRNA) nucleotide sequences using aminoacyl-tRNA substrates. Current knowledge of the decoding mechanism derives principally from studies on bacterial systems1. Although key features are conserved across evolution2, eukaryotes achieve higher-fidelity mRNA decoding than bacteria3. In human, changes in decoding fidelity are linked to ageing and disease and represent a potential point of therapeutic intervention in both viral and cancer treatment4-6. Here we combine single-molecule imaging and cryogenic electron microscopy methods to examine the molecular basis of human ribosome fidelity to reveal that the decoding mechanism is both kinetically and structurally distinct from that of bacteria. Although decoding is globally analogous in both species, the reaction coordinate of aminoacyl-tRNA movement is altered on the human ribosome and the process is an order of magnitude slower. These distinctions arise from eukaryote-specific structural elements in the human ribosome and in the elongation factor eukaryotic elongation factor 1A (eEF1A) that together coordinate faithful tRNA incorporation at each mRNA codon. The distinct nature and timing of conformational changes within the ribosome and eEF1A rationalize how increased decoding fidelity is achieved and potentially regulated in eukaryotic species.

Computationally-guided design and selection of high performing ribosomal active site mutants

Understanding how modifications to the ribosome affect function has implications for studying ribosome biogenesis, building minimal cells, and repurposing ribosomes for synthetic biology. However, efforts to design sequence-modified ribosomes have been limited because point mutations in the ribosomal RNA (rRNA), especially in the catalytic active site (peptidyl transferase center; PTC), are often functionally detrimental. Moreover, methods for directed evolution of rRNA are constrained by practical considerations (e.g. library size). Here, to address these limitations, we developed a computational rRNA design approach for screening guided libraries of mutant ribosomes. Our method includes in silico library design and selection using a Rosetta stepwise Monte Carlo method (SWM), library construction and in vitro testing of combined ribosomal assembly and translation activity, and functional characterization in vivo. As a model, we apply our method to making modified ribosomes with mutant PTCs. We engineer ribosomes with as many as 30 mutations in their PTCs, highlighting previously unidentified epistatic interactions, and show that SWM helps identify sequences with beneficial phenotypes as compared to random library sequences. We further demonstrate that some variants improve cell growth in vivo, relative to wild type ribosomes. We anticipate that SWM design and selection may serve as a powerful tool for rRNA engineering.

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.

Didemnin B and ternatin-4 differentially inhibit conformational changes in eEF1A required for aminoacyl-tRNA accommodation into mammalian ribosomes

Rapid and accurate mRNA translation requires efficient codon-dependent delivery of the correct aminoacyl-tRNA (aa-tRNA) to the ribosomal A site. In mammals, this fidelity-determining reaction is facilitated by the GTPase elongation factor-1 alpha (eEF1A), which escorts aa-tRNA as an eEF1A(GTP)-aa-tRNA ternary complex into the ribosome. The structurally unrelated cyclic peptides didemnin B and ternatin-4 bind to the eEF1A(GTP)-aa-tRNA ternary complex and inhibit translation but have different effects on protein synthesis in vitro and in vivo. Here, we employ single-molecule fluorescence imaging and cryogenic electron microscopy to determine how these natural products inhibit translational elongation on mammalian ribosomes. By binding to a common site on eEF1A, didemnin B and ternatin-4 trap eEF1A in an intermediate state of aa-tRNA selection, preventing eEF1A release and aa-tRNA accommodation on the ribosome. We also show that didemnin B and ternatin-4 exhibit distinct effects on the dynamics of aa-tRNA selection that inform on observed disparities in their inhibition efficacies and physiological impacts. These integrated findings underscore the value of dynamics measurements in assessing the mechanism of small-molecule inhibition and highlight potential of single-molecule methods to reveal how distinct natural products differentially impact the human translation mechanism.

Exploring the Energy Landscape of Riboswitches Using Collective Variables Based on Tertiary Contacts

Messenger RNA regulatory elements, such as riboswitches, can display a high degree of flexibility. By characterizing their energy landscapes, and corresponding distributions of 3D configurations, structure-function relationships can be elucidated. Molecular dynamics simulation with enhanced sampling is an important strategy used to computationally access free energy landscapes characterizing the accessible 3D conformations of RNAs. While tertiary contacts are thought to play important roles in RNA dynamics, it is difficult, in explicit solvent, to sample the formation and breakage of tertiary contacts, such as helix-helix interactions, pseudoknot interactions, and junction interactions, while maintaining intact secondary structure elements. To this end, we extend previously developed collective variables and metadynamics efforts, to establish a simple metadynamics protocol, which utilizes only one collective variable, based on multiple tertiary contacts, to characterize the underlying free energy landscape of any RNA molecule. We develop a modified collective variable, the tertiary contacts distance (Q_TC), which can probe the formation and breakage of all or selectively chosen tertiary contacts of the RNA. The SAM-I riboswitch in the presence of three ionic and substrate conditions was investigated and validated against the structure ensemble previously generated using SAXS experiments. This efficient and easy to implement all-atom MD simulation based approach incorporating metadynamics to study RNA conformational dynamics can also be transferred to any other type of biomolecule.

Hyper-swivel head domain motions are required for complete mRNA-tRNA translocation and ribosome resetting

Translocation of messenger RNA (mRNA) and transfer RNA (tRNA) substrates through the ribosome during protein synthesis, an exemplar of directional molecular movement in biology, entails a complex interplay of conformational, compositional, and chemical changes. The molecular determinants of early translocation steps have been investigated rigorously. However, the elements enabling the ribosome to complete translocation and reset for subsequent protein synthesis reactions remain poorly understood. Here, we have combined molecular simulations with single-molecule fluorescence resonance energy transfer imaging to gain insights into the rate-limiting events of the translocation mechanism. We find that diffusive motions of the ribosomal small subunit head domain to hyper-swivelled positions, governed by universally conserved rRNA, can maneuver the mRNA and tRNAs to their fully translocated positions. Subsequent engagement of peptidyl-tRNA and disengagement of deacyl-tRNA from mRNA, within their respective small subunit binding sites, facilitate the ribosome resetting mechanism after translocation has occurred to enable protein synthesis to resume.

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.

Exploring cryo-electron microscopy with molecular dynamics

Single particle analysis cryo-electron microscopy (EM) and molecular dynamics (MD) have been complimentary methods since cryo-EM was first applied to the field of structural biology. The relationship started by biasing structural models to fit low-resolution cryo-EM maps of large macromolecular complexes not amenable to crystallization. The connection between cryo-EM and MD evolved as cryo-EM maps improved in resolution, allowing advanced sampling algorithms to simultaneously refine backbone and sidechains. Moving beyond a single static snapshot, modern inferencing approaches integrate cryo-EM and MD to generate structural ensembles from cryo-EM map data or directly from the particle images themselves. We summarize the recent history of MD innovations in the area of cryo-EM modeling. The merits for the myriad of MD based cryo-EM modeling methods are discussed, as well as, the discoveries that were made possible by the integration of molecular modeling with cryo-EM. Lastly, current challenges and potential opportunities are reviewed.

The Structural Dynamics of Translation

Accurate protein synthesis (translation) relies on translation factors that rectify ribosome fluctuations into a unidirectional process. Understanding this process requires structural characterization of the ribosome and translation-factor dynamics. In the 2000s, crystallographic studies determined high-resolution structures of ribosomes stalled with translation factors, providing a starting point for visualizing translation. Recent progress in single-particle cryogenic electron microscopy (cryo-EM) has enabled near-atomic resolution of numerous structures sampled in heterogeneous complexes (ensembles). Ensemble and time-resolved cryo-EM have now revealed unprecedented views of ribosome transitions in the three principal stages of translation: initiation, elongation, and termination. This review focuses on how translation factors help achieve high accuracy and efficiency of translation by monitoring distinct ribosome conformations and by differentially shifting the equilibria of ribosome rearrangements for cognate and near-cognate substrates. Expected final online publication date for the Annual Review of Biochemistry, Volume 91 is June 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Interplay between an ATP-binding cassette F protein and the ribosome from Mycobacterium tuberculosis

EttA, energy-dependent translational throttle A, is a ribosomal factor that gates ribosome entry into the translation elongation cycle. A detailed understanding of its mechanism of action is limited due to the lack of high-resolution structures along its ATPase cycle. Here we present the cryo-electron microscopy (cryo-EM) structures of EttA from Mycobacterium tuberculosis (Mtb), referred to as MtbEttA, in complex with the Mtb 70S ribosome initiation complex (70SIC) at the pre-hydrolysis (ADPNP) and transition (ADP-VO4) states, and the crystal structure of MtbEttA alone in the post-hydrolysis (ADP) state. We observe that MtbEttA binds the E-site of the Mtb 70SIC, remodeling the P-site tRNA and the ribosomal intersubunit bridge B7a during the ribosomal ratcheting. In return, the rotation of the 30S causes conformational changes in MtbEttA, forcing the two nucleotide-binding sites (NBSs) to alternate to engage each ADPNP in the pre-hydrolysis states, followed by complete engagements of both ADP-VO4 molecules in the ATP-hydrolysis transition states. In the post-hydrolysis state, the conserved ATP-hydrolysis motifs of MtbEttA dissociate from both ADP molecules, leaving two nucleotide-binding domains (NBDs) in an open conformation. These structures reveal a dynamic interplay between MtbEttA and the Mtb ribosome, providing insights into the mechanism of translational regulation by EttA-like proteins.

The energetics of subunit rotation in the ribosome

Protein synthesis in the cell is controlled by an elaborate sequence of conformational rearrangements in the ribosome. The composition of a ribosome varies by species, though they typically contain ∼ 50-100 RNA and protein molecules. While advances in structural techniques have revolutionized our understanding of long-lived conformational states, a vast range of transiently visited configurations can not be directly observed. In these cases, computational/simulation methods can be used to understand the mechanical properties of the ribosome. Insights from these approaches can then help guide next-generation experimental measurements. In this short review, we discuss theoretical strategies that have been deployed to quantitatively describe the energetics of collective rearrangements in the ribosome. We focus on efforts to probe large-scale subunit rotation events, which involve the coordinated displacement of large numbers of atoms (tens of thousands). These investigations are revealing how the molecular structure of the ribosome encodes the mechanical properties that control large-scale dynamics.

Exploring the Minimum-Energy Pathways and Free-Energy Profiles of Enzymatic Reactions with QM/MM Calculations

Understanding molecular mechanisms of enzymatic reactions is of vital importance in biochemistry and biophysics. Here, we introduce new functions of hybrid quantum mechanical/molecular mechanical (QM/MM) calculations in the GENESIS program to compute the minimum-energy pathways (MEPs) and free-energy profiles of enzymatic reactions. For this purpose, an interface in GENESIS is developed to utilize a highly parallel electronic structure program, QSimulate-QM (https://qsimulate.com), calling it as a shared library from GENESIS. Second, algorithms to search the MEP are implemented, combining the string method (E et al. J. Chem. Phys. 2007, 126, 164103) with the energy minimization of the buffer MM region. The method implemented in GENESIS is applied to an enzyme, triosephosphate isomerase, which converts dihyroxyacetone phosphate to glyceraldehyde 3-phosphate in four proton-transfer processes. QM/MM-molecular dynamics simulations show performances of greater than 1 ns/day with the density functional tight binding (DFTB), and 10-30 ps/day with the hybrid density functional theory, B3LYP-D3. These performances allow us to compute not only MEP but also the potential of mean force (PMF) of the enzymatic reactions using the QM/MM calculations. The barrier height obtained as 13 kcal mol-1 with B3LYP-D3 in the QM/MM calculation is in agreement with the experimental results. The impact of conformational sampling in PMF calculations and the level of electronic structure calculations (DFTB vs B3LYP-D3) suggests reliable computational protocols for enzymatic reactions without high computational costs.

Dynamics of uS19 C-Terminal Tail during the Translation Elongation Cycle in Human Ribosomes

Ribosomes undergo multiple conformational transitions during translation elongation. Here, we report the high-resolution cryoelectron microscopy (cryo-EM) structure of the human 80S ribosome in the post-decoding pre-translocation state (classical-PRE) at 3.3-Å resolution along with the rotated (hybrid-PRE) and the post-translocation states (POST). The classical-PRE state ribosome structure reveals a previously unobserved interaction between the C-terminal region of the conserved ribosomal protein uS19 and the A- and P-site tRNAs and the mRNA in the decoding site. In addition to changes in the inter-subunit bridges, analysis of different ribosomal conformations reveals the dynamic nature of this domain and suggests a role in tRNA accommodation and translocation during elongation. Furthermore, we show that disease-associated mutations in uS19 result in increased frameshifting. Together, this structure-function analysis provides mechanistic insights into the role of the uS19 C-terminal tail in the context of mammalian ribosomes.

FRET-based dynamic structural biology: Challenges, perspectives and an appeal for open-science practices

Single-molecule FRET (smFRET) has become a mainstream technique for studying biomolecular structural dynamics. The rapid and wide adoption of smFRET experiments by an ever-increasing number of groups has generated significant progress in sample preparation, measurement procedures, data analysis, algorithms and documentation. Several labs that employ smFRET approaches have joined forces to inform the smFRET community about streamlining how to perform experiments and analyze results for obtaining quantitative information on biomolecular structure and dynamics. The recent efforts include blind tests to assess the accuracy and the precision of smFRET experiments among different labs using various procedures. These multi-lab studies have led to the development of smFRET procedures and documentation, which are important when submitting entries into the archiving system for integrative structure models, PDB-Dev. This position paper describes the current 'state of the art' from different perspectives, points to unresolved methodological issues for quantitative structural studies, provides a set of 'soft recommendations' about which an emerging consensus exists, and lists openly available resources for newcomers and seasoned practitioners. To make further progress, we strongly encourage 'open science' practices.

Multi-Scale Flexible Fitting of Proteins to Cryo-EM Density Maps at Medium Resolution

Structure determination using cryo-electron microscopy (cryo-EM) medium-resolution density maps is often facilitated by flexible fitting. Avoiding overfitting, adjusting force constants driving the structure to the density map, and emulating complex conformational transitions are major concerns in the fitting. To address them, we develop a new method based on a three-step multi-scale protocol. First, flexible fitting molecular dynamics (MD) simulations with coarse-grained structure-based force field and replica-exchange scheme between different force constants replicas are performed. Second, fitted Cα atom positions guide the all-atom structure in targeted MD. Finally, the all-atom flexible fitting refinement in implicit solvent adjusts the positions of the side chains in the density map. Final models obtained via the multi-scale protocol are significantly better resolved and more reliable in comparison with long all-atom flexible fitting simulations. The protocol is useful for multi-domain systems with intricate structural transitions as it preserves the secondary structure of single domains.

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.

Intrinsically Disordered Regions of the DNA-Binding Domain of Human FoxP1 Facilitate Domain Swapping

Forkhead box P (FoxP) proteins are unique transcription factors that spatiotemporally regulate gene expression by tethering two chromosome loci together via functional domain-swapped dimers formed through their DNA-binding domains. Further, the differential kinetics on this dimerization mechanism underlie an intricate gene regulation network at physiological conditions. Nonetheless, poor understanding of the structural dynamics and steps of the association process impedes to link the functional domain swapping to human-associated diseases. Here, we have characterized the DNA-binding domain of human FoxP1 by integrating single-molecule Förster resonance energy transfer and hydrogen-deuterium exchange mass spectrometry data with molecular dynamics simulations. Our results confirm the formation of a previously postulated domain-swapped (DS) FoxP1 dimer in solution and reveal the presence of highly populated, heterogeneous, and locally disordered dimeric intermediates along the dimer dissociation pathway. The unique features of FoxP1 provide a glimpse of how intrinsically disordered regions can facilitate domain swapping oligomerization and other tightly regulated association mechanisms relevant in biological processes.

Structure and function of yeast Lso2 and human CCDC124 bound to hibernating ribosomes

Cells adjust to nutrient deprivation by reversible translational shutdown. This is accompanied by maintaining inactive ribosomes in a hibernation state, in which they are bound by proteins with inhibitory and protective functions. In eukaryotes, such a function was attributed to suppressor of target of Myb protein 1 (Stm1; SERPINE1 mRNA-binding protein 1 [SERBP1] in mammals), and recently, late-annotated short open reading frame 2 (Lso2; coiled-coil domain containing short open reading frame 124 [CCDC124] in mammals) was found to be involved in translational recovery after starvation from stationary phase. Here, we present cryo-electron microscopy (cryo-EM) structures of translationally inactive yeast and human ribosomes. We found Lso2/CCDC124 accumulating on idle ribosomes in the nonrotated state, in contrast to Stm1/SERBP1-bound ribosomes, which display a rotated state. Lso2/CCDC124 bridges the decoding sites of the small with the GTPase activating center (GAC) of the large subunit. This position allows accommodation of the duplication of multilocus region 34 protein (Dom34)-dependent ribosome recycling system, which splits Lso2-containing, but not Stm1-containing, ribosomes. We propose a model in which Lso2 facilitates rapid translation reactivation by stabilizing the recycling-competent state of inactive ribosomes.

Cryo-EM of elongating ribosome with EF-Tu•GTP elucidates tRNA proofreading

Ribosomes accurately decode mRNA by proofreading each aminoacyl-tRNA delivered by elongation factor EF-Tu¹. Understanding the molecular mechanism of proofreading requires visualizing GTP-catalyzed elongation, which has remained a challenge^2–4. Here, time-resolved cryo-EM revealed 33 states following aminoacyl-tRNA delivery by EF-Tu•GTP. Instead of locking cognate tRNA upon initial recognition, the ribosomal decoding center (DC) dynamically monitors codon-anticodon interactions before and after GTP hydrolysis. GTP hydrolysis allows EF-Tu’s GTPase domain to extend away, releasing EF-Tu from tRNA. Then, the 30S subunit locks cognate tRNA in the DC, and rotates, enabling the tRNA to bypass 50S protrusions during accommodation into the peptidyl transferase center. By contrast, the DC fails to lock near-cognate tRNA, allowing dissociation of near-cognate tRNA during both initial selection (before GTP hydrolysis) and proofreading (after GTP hydrolysis). These findings reveal structural similarity between initial selection^5,6 and the previously unseen proofreading, which together govern efficient rejection of incorrect tRNA.

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.

Elongation factor-Tu can repetitively engage aminoacyl-tRNA within the ribosome during the proofreading stage of tRNA selection

Elongation factor Tu (EF-Tu) facilitates rapid and accurate selection of aminoacyl-tRNA (aa-tRNA) by the bacterial ribosome during protein synthesis. We show that EF-Tu dissociates from the ribosome as aa-tRNA navigates the accommodation corridor en route to peptide bond formation. We find that EF-Tu’s release from the ribosome during aa-tRNA selection can be reversible. We also demonstrate that new ternary complex formation, accompanied by futile cycles of GTP hydrolysis, can occur on aa-tRNA bound within the ribosome. These findings inform on the decoding mechanism, the contributions of EF-Tu to the fidelity of translation, and the potential consequences of reduced rates of peptide bond formation on cellular physiology.

The substrate for ribosomes actively engaged in protein synthesis is a ternary complex of elongation factor Tu (EF-Tu), aminoacyl-tRNA (aa-tRNA), and GTP. EF-Tu plays a critical role in mRNA decoding by increasing the rate and fidelity of aa-tRNA selection at each mRNA codon. Here, using three-color single-molecule fluorescence resonance energy transfer imaging and molecular dynamics simulations, we examine the timing and role of conformational events that mediate the release of aa-tRNA from EF-Tu and EF-Tu from the ribosome after GTP hydrolysis. Our investigations reveal that conformational changes in EF-Tu coordinate the rate-limiting passage of aa-tRNA through the accommodation corridor en route to the peptidyl transferase center of the large ribosomal subunit. Experiments using distinct inhibitors of the accommodation process further show that aa-tRNA must at least partially transit the accommodation corridor for EF-Tu⋅GDP to release. aa-tRNAs failing to undergo peptide bond formation at the end of accommodation corridor passage after EF-Tu release can be reengaged by EF-Tu⋅GTP from solution, coupled to GTP hydrolysis. These observations suggest that additional rounds of ternary complex formation can occur on the ribosome during proofreading, particularly when peptide bond formation is slow, which may serve to increase both the rate and fidelity of protein synthesis at the expense of GTP hydrolysis.

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.

tRNA Dissociation from EF-Tu after GTP Hydrolysis: Primary Steps and Antibiotic Inhibition

In each round of ribosomal translation, the translational GTPase elongation factor Tu (EF-Tu) delivers a transfer RNA (tRNA) to the ribosome. After successful decoding, EF-Tu hydrolyzes GTP, which triggers a conformational change that ultimately results in the release of the tRNA from EF-Tu. To identify the primary steps of these conformational changes and how they are prevented by the antibiotic kirromycin, we employed all-atom explicit-solvent molecular dynamics simulations of the full ribosome-EF-Tu complex. Our results suggest that after GTP hydrolysis and Pi release, the loss of interactions between the nucleotide and the switch 1 loop of EF-Tu allows domain D1 of EF-Tu to rotate relative to domains D2 and D3 and leads to an increased flexibility of the switch 1 loop. This rotation induces a closing of the D1-D3 interface and an opening of the D1-D2 interface. We propose that the opening of the D1-D2 interface, which binds the CCA tail of the tRNA, weakens the crucial EF-Tu-tRNA interactions, which lowers tRNA binding affinity, representing the first step of tRNA release. Kirromycin binds within the D1-D3 interface, sterically blocking its closure, but does not prevent hydrolysis. The resulting increased flexibility of switch 1 explains why it is not resolved in kirromycin-bound structures.

Gō model revisited

This review discusses Gō models broadly used in biomolecular simulations. I start with a brief description of the original lattice model study by Nobuhiro Gō. Then, the theory of protein folding behind Gō model, free energy approaches, and off-lattice Gō models are reviewed. I also mention a stringent test for the assumption in Gō models given from all-atom molecular dynamics simulations. Subsequently, I move to application of Gō models to protein dynamical functions. Various extension of Gō models is also reviewed. Finally, some publicly available tools to use Gō models are listed.

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.

Large-scale simulations of nucleoprotein complexes: ribosomes, nucleosomes, chromatin, chromosomes and CRISPR

Recent advances in biotechnology such as Hi-C, CRISPR/Cas9 and ribosome display have placed nucleoprotein complexes at center stage. Understanding the structural dynamics of these complexes aids in optimizing protocols and interpreting data for these new technologies. The integration of simulation and experiment has helped advance mechanistic understanding of these systems. Coarse-grained simulations, reduced-description models, and explicit solvent molecular dynamics simulations yield useful complementary perspectives on nucleoprotein complex structural dynamics. When combined with Hi-C, cryo-EM, and single molecule measurements, these simulations integrate disparate forms of experimental data into a coherent mechanism.

Rapid interpretation of small-angle X-ray scattering data

The fundamental aim of structural analyses in biophysics is to reveal a mutual relation between a molecule’s dynamic structure and its physiological function. Small-angle X-ray scattering (SAXS) is an experimental technique for structural characterization of macromolecules in solution and enables time-resolved analysis of conformational changes under physiological conditions. As such experiments measure spatially averaged low-resolution scattering intensities only, the sparse information obtained is not sufficient to uniquely reconstruct a three-dimensional atomistic model. Here, we integrate the information from SAXS into molecular dynamics simulations using computationally efficient native structure-based models. Dynamically fitting an initial structure towards a scattering intensity, such simulations produce atomistic models in agreement with the target data. In this way, SAXS data can be rapidly interpreted while retaining physico-chemical knowledge and sampling power of the underlying force field. We demonstrate our method’s performance using the example of three protein systems. Simulations are faster than full molecular dynamics approaches by more than two orders of magnitude and consistently achieve comparable accuracy. Computational demands are reduced sufficiently to run the simulations on commodity desktop computers instead of high-performance computing systems. These results underline that scattering-guided structure-based simulations provide a suitable framework for rapid early-stage refinement of structures towards SAXS data with particular focus on minimal computational resources and time.

Proteins are the molecular nanomachines in biological cells and thus vital to any known form of life. From the evolutionary perspective, viable protein structure emerges on the basis of a ‘form-follows-function’ principle. A protein’s designated function is inextricably linked to dynamic conformational changes, which can be observed by small-angle X-ray scattering. Intensities from SAXS contain low-resolution information on the protein’s shape at different steps of its functional cycle. We are interested in directly getting an atomistic model of this encoded structure. One powerful approach is to include the experimental data into computational simulations of the protein’s function-related physical motions. We combine scattering intensities with coarse-grained native structure-based models. These models are computationally highly efficient yet describe the system’s dynamics realistically. Here, we present our method for rapid interpretation of scattering intensities from SAXS to derive structural models, using minimal computational resources and time.

Structure-Based Model of RNA Pseudoknot Captures Magnesium-Dependent Folding Thermodynamics

We develop a simple, coarse-grained approach for simulating the folding of the Beet Western Yellow Virus (BWYV) pseudoknot toward the goal of creating a transferable model that can be used to study other small RNA molecules. This approach combines a structure-based model (SBM) of RNA with an electrostatic scheme that has previously been shown to correctly reproduce ionic condensation in the native basin. Mg²⁺ ions are represented explicitly, directly incorporating ion-ion correlations into the system, and K⁺ is represented implicitly, through the mean-field generalized Manning counterion condensation theory. Combining the electrostatic scheme with a SBM enables the electrostatic scheme to be tested beyond the native basin. We calibrate the SBM to reproduce experimental BWYV unfolding data by eliminating overstabilizing backbone interactions from the molecular contact map and by strengthening base pairing and stacking contacts relative to other native contacts, consistent with the experimental observation that relative helical stabilities are central determinants of the RNA unfolding sequence. We find that this approach quantitatively captures the Mg²⁺ dependence of the folding temperature and generates intermediate states that better approximate those revealed by experiment. Finally, we examine how our model captures Mg²⁺ condensation about the BWYV pseudoknot and a U-tail variant, for which the nine 3' end nucleotides are replaced with uracils, and find our results to be consistent with experimental condensation measurements. This approach can be easily transferred to other RNA molecules by eliminating and strengthening the same classes of contacts in the SBM and including generalized Manning counterion condensation.

Using SMOG 2 to Simulate Complex Biomolecular Assemblies

Over the last 20 years, the application of structure-based (Gō-like) models has ranged from protein folding with coarse-grained models to all-atom representations of large-scale molecular assemblies. While there are many variants that may be employed, the common feature of these models is that some (or all) of the stabilizing energetic interactions are defined based on the knowledge of a particular experimentally obtained conformation. With the generality of this approach, there was a need for a versatile computational platform for designing and implementing this class of models. To this end, the SMOG 2 software package provides an easy-to-use interface, where the user has full control of the model parameters. This software allows the user to edit XML-formatted files in order to provide definitions of new structure-based models. SMOG 2 reads these "template" files and maps the interactions onto specific structures, which are provided in PDB format. The force field files produced by SMOG 2 may then be used to perform simulations with a variety of popular molecular dynamics suites. In this chapter, we describe some of the key features of the SMOG 2 package, while providing examples and strategies for applying these techniques to complex (often large-scale) molecular assemblies, such as 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.

Ligand Modulates Cross-Coupling between Riboswitch Folding and Transcriptional Pausing

Numerous classes of riboswitches have been found to regulate bacterial gene expression in response to physiological cues, offering new paths to antibacterial drugs. As common studies of isolated riboswitches lack the functional context of the transcription machinery, we here combine single-molecule, biochemical, and simulation approaches to investigate the coupling between co-transcriptional folding of the pseudoknot-structured preQ1 riboswitch and RNA polymerase (RNAP) pausing. We show that pausing at a site immediately downstream of the riboswitch requires a ligand-free pseudoknot in the nascent RNA, a precisely spaced sequence resembling the pause consensus, and electrostatic and steric interactions with the RNAP exit channel. While interactions with RNAP stabilize the native fold of the riboswitch, binding of the ligand signals RNAP release from the pause. Our results demonstrate that the nascent riboswitch and its ligand actively modulate the function of RNAP and vice versa, a paradigm likely to apply to other cellular RNA transcripts.

RNA as a Complex Polymer with Coupled Dynamics of Ions and Water in the Outer Solvation Sphere

We unravel the internal and collective modes of a widely studied 58-nucleotide rRNA fragment in solvent using atomically detailed molecular dynamics simulations. The variation of lifetimes for water hydrogen bonds with nucleotide groups indicates heterogeneity of water dynamics on the RNA surface. The time scales of interactions of the discrete water molecules with RNA nucleotides extend from several hundred picoseconds to a few nanoseconds. We determine all of the association sites and the spatial distribution of residence times for Mg²⁺, K⁺, and water molecules in those sites. We provide insights into the population of Mg²⁺ and K⁺ ions and water molecules in the outer sphere and how their fluctuations are intricately linked with the kinetics of the 58-mer. We find that many of the long-lived Mg²⁺ sites identified from the simulations agree with the locations of ions in the X-ray structure. We determine the excess ion atmosphere around the rRNA and compare it with experimental data. We investigate the collective behavior of RNA, ions, and water, by performing a joint principle component analysis for the Cartesian coordinates of the RNA phosphorus atoms and for the occupation counts of the association sites. Our results indicate that the 58-mer system is a complex polymer, composed of RNA that is encased by a fluctuating network of associated counterions, co-ions, and water.

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.