The role of L1 stalk-tRNA interaction in the ribosome elongation cycle

On the Power and Challenges of Atomistic Molecular Dynamics to Investigate RNA Molecules

RNA molecules play a vital role in biological processes within the cell, with significant implications for science and medicine. Notably, the biological functions exerted by specific RNA molecules are often linked to the RNA conformational ensemble. However, the experimental characterization of such three-dimensional RNA structures is challenged by the structural heterogeneity of RNA and by its multiple dynamic interactions with binding partners such as small molecules, proteins, and metal ions. Consequently, our current understanding of the structure-function relationship of RNA molecules is still limited. In this context, we highlight molecular dynamics (MD) simulations as a powerful tool to complement experimental efforts on RNAs. Despite the recognized limitations of current force fields for RNA MD simulations, examining the dynamics of selected RNAs has provided valuable functional insights into their structures.

The Beak of Eukaryotic Ribosomes: Life, Work and Miracles

Ribosomes are not totally globular machines. Instead, they comprise prominent structural protrusions and a myriad of tentacle-like projections, which are frequently made up of ribosomal RNA expansion segments and N- or C-terminal extensions of ribosomal proteins. This is more evident in higher eukaryotic ribosomes. One of the most characteristic protrusions, present in small ribosomal subunits in all three domains of life, is the so-called beak, which is relevant for the function and regulation of the ribosome's activities. During evolution, the beak has transitioned from an all ribosomal RNA structure (helix h33 in 16S rRNA) in bacteria, to an arrangement formed by three ribosomal proteins, eS10, eS12 and eS31, and a smaller h33 ribosomal RNA in eukaryotes. In this review, we describe the different structural and functional properties of the eukaryotic beak. We discuss the state-of-the-art concerning its composition and functional significance, including other processes apparently not related to translation, and the dynamics of its assembly in yeast and human cells. Moreover, we outline the current view about the relevance of the beak's components in human diseases, especially in ribosomopathies and cancer.

Extraribosomal Functions of Bacterial Ribosomal Proteins-An Update, 2023

Ribosomal proteins (r-proteins) are abundant, highly conserved, and multifaceted cellular proteins in all domains of life. Most r-proteins have RNA-binding properties and can form protein-protein contacts. Bacterial r-proteins govern the co-transcriptional rRNA folding during ribosome assembly and participate in the formation of the ribosome functional sites, such as the mRNA-binding site, tRNA-binding sites, the peptidyl transferase center, and the protein exit tunnel. In addition to their primary role in a cell as integral components of the protein synthesis machinery, many r-proteins can function beyond the ribosome (the phenomenon known as moonlighting), acting either as individual regulatory proteins or in complexes with various cellular components. The extraribosomal activities of r-proteins have been studied over the decades. In the past decade, our understanding of r-protein functions has advanced significantly due to intensive studies on ribosomes and gene expression mechanisms not only in model bacteria like Escherichia coli or Bacillus subtilis but also in little-explored bacterial species from various phyla. The aim of this review is to update information on the multiple functions of r-proteins in bacteria.

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.

Minimization of the E. coli ribosome, aided and optimized by community science

The ribosome is a ribonucleoprotein complex found in all domains of life. Its role is to catalyze protein synthesis, the messenger RNA (mRNA)-templated formation of amide bonds between α-amino acid monomers. Amide bond formation occurs within a highly conserved region of the large ribosomal subunit known as the peptidyl transferase center (PTC). Here we describe the step-wise design and characterization of mini-PTC 1.1, a 284-nucleotide RNA that recapitulates many essential features of the Escherichia coli PTC. Mini-PTC 1.1 folds into a PTC-like structure under physiological conditions, even in the absence of r-proteins, and engages small molecule analogs of A- and P-site tRNAs. The sequence of mini-PTC 1.1 differs from the wild type E. coli ribosome at 12 nucleotides that were installed by a cohort of citizen scientists using the on-line video game Eterna. These base changes improve both the secondary structure and tertiary folding of mini-PTC 1.1 as well as its ability to bind small molecule substrate analogs. Here, the combined input from Eterna citizen-scientists and RNA structural analysis provides a robust workflow for the design of a minimal PTC that recapitulates many features of an intact ribosome.

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.

Features of smaller ribosomes in candidate phyla radiation (CPR) bacteria revealed with a molecular evolutionary analysis

The candidate phyla radiation (CPR) is a large bacterial group consisting mainly of uncultured lineages. They have small cells and small genomes, and they often lack ribosomal proteins uL1, bL9, and/or uL30, which are basically ubiquitous in non-CPR bacteria. Here, we comprehensively analyzed the genomic information on CPR bacteria and identified their unique properties. The distribution of protein lengths in CPR bacteria peaks at around 100-150 amino acids, whereas the position of the peak varies in the range of 100-300 amino acids in free-living non-CPR bacteria, and at around 100-200 amino acids in most symbiotic non-CPR bacteria. These results show that the proteins of CPR bacteria are smaller, on average, than those of free-living non-CPR bacteria, like those of symbiotic non-CPR bacteria. We found that ribosomal proteins bL28, uL29, bL32, and bL33 have been lost in CPR bacteria in a taxonomic lineage-specific manner. Moreover, the sequences of approximately half of all ribosomal proteins of CPR differ, in part, from those of non-CPR bacteria, with missing regions or specifically added regions. We also found that several regions in the 16S, 23S, and 5S rRNAs of CPR bacteria are lacking, which presumably caused the total predicted lengths of the three rRNAs of CPR bacteria to be smaller than those of non-CPR bacteria. The regions missing in the CPR ribosomal proteins and rRNAs are located near the surface of the ribosome, and some are close to one another. These observations suggest that ribosomes are smaller in CPR bacteria than those in free-living non-CPR bacteria, with simplified surface structures.

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.

Structural Studies Reveal the Role of Helix 68 in the Elongation Step of Protein Biosynthesis

The ribosome, a multicomponent assembly consisting of RNA and proteins, is a pivotal macromolecular machine that translates the genetic code into proteins. The large ribosomal subunit rRNA helix 68 (H68) is a key element in the protein synthesis process, as it coordinates the coupled movements of the actors involved in translocation, including the tRNAs and L1 stalk. Examination of cryo-electron microscopy (cryo-EM) structures of ribosomes incubated for various time durations at physiological temperatures led to the identification of functionally relevant H68 movements. These movements assist the transition of the L1 stalk between its open and closed states. H68 spatial flexibility and its significance to the protein synthesis process were confirmed through its effective targeting with antisense PNA oligomers. Our results suggest that H68 is actively involved in ribosome movements that are central to the elongation process. IMPORTANCE The mechanism that regulates the translocation step in ribosomes during protein synthesis is not fully understood. In this work, cryo-EM techniques used to image ribosomes from Staphylococcus aureus after incubation at physiological temperature allowed the identification of a conformation of the helix 68 that has never been observed so far. We then propose a mechanism in which such helix, switching between two different conformations, actively coordinates the translocation step, shedding light on the dynamics of ribosomal components. In addition, the relevance of helix 68 to ribosome function and its potential as an antibiotic target was proved by inhibiting Staphylococcus aureus ribosomes activity in vitro using oligomers with sequence complementarity.

Types and Functions of Mitoribosome-Specific Ribosomal Proteins across Eukaryotes

Mitochondria are key organelles that combine features inherited from their bacterial endosymbiotic ancestor with traits that arose during eukaryote evolution. These energy producing organelles have retained a genome and fully functional gene expression machineries including specific ribosomes. Recent advances in cryo-electron microscopy have enabled the characterization of a fast-growing number of the low abundant membrane-bound mitochondrial ribosomes. Surprisingly, mitoribosomes were found to be extremely diverse both in terms of structure and composition. Still, all of them drastically increased their number of ribosomal proteins. Interestingly, among the more than 130 novel ribosomal proteins identified to date in mitochondria, most of them are composed of a-helices. Many of them belong to the nuclear encoded super family of helical repeat proteins. Here we review the diversity of functions and the mode of action held by the novel mitoribosome proteins and discuss why these proteins that share similar helical folds were independently recruited by mitoribosomes during evolution in independent eukaryote clades.

Exploring Allosteric Signaling in the Exit Tunnel of the Bacterial Ribosome by Molecular Dynamics Simulations and Residue Network Model

The bacterial ribosomal tunnel is equipped with numerous sites highly sensitive to the course of the translation process. This study investigates allosteric pathways linking distant functional sites that collaboratively play a role either in translation regulation or recruitment of chaperones. We apply perturbation response scanning (PRS) analysis to 700 ns long and 500 ns long coarse-grained molecular dynamics simulations of E. coli and T. thermophilus large subunits, respectively, to reveal nucleotides/residues with the ability to transmit perturbations by dynamic rationale. We also use the residue network model with the k-shortest pathways method to calculate suboptimal pathways based on the contact topology of the ribosomal tunnel of E. coli crystal structure and 101 ClustENM generated conformers of T. thermophilus large subunit. In the upper part of the tunnel, results suggest that A2062 and A2451 can communicate in both directions for translation stalling, mostly through dynamically coupled C2063, C2064, and A2450. For a similar purpose, U2585 and U2586 are coupled with A2062, while they are also sensitive to uL4 and uL22 at the constriction region through two different pathways at the opposite sides of the tunnel wall. In addition, the constriction region communicates with the chaperone binding site on uL23 at the solvent side but through few nucleotides. Potential allosteric communication pathways between the lower part of the tunnel and chaperone binding site mostly use the flexible loop of uL23, while A1336–G1339 provide a suboptimal pathway. Both species seem to employ similar mechanisms in the long tunnel, where a non-conserved cavity at the bacterial uL23 and 23S rRNA interface is proposed as a novel drug target.

Structural insights into mammalian mitochondrial translation elongation catalyzed by mtEFG1

Mitochondria are eukaryotic organelles of bacterial origin where respiration takes place to produce cellular chemical energy. These reactions are catalyzed by the respiratory chain complexes located in the inner mitochondrial membrane. Notably, key components of the respiratory chain complexes are encoded on the mitochondrial chromosome and their expression relies on a dedicated mitochondrial translation machinery. Defects in the mitochondrial gene expression machinery lead to a variety of diseases in humans mostly affecting tissues with high energy demand such as the nervous system, the heart, or the muscles. The mitochondrial translation system has substantially diverged from its bacterial ancestor, including alterations in the mitoribosomal architecture, multiple changes to the set of translation factors and striking reductions in otherwise conserved tRNA elements. Although a number of structures of mitochondrial ribosomes from different species have been determined, our mechanistic understanding of the mitochondrial translation cycle remains largely unexplored. Here, we present two cryo-EM reconstructions of human mitochondrial elongation factor G1 bound to the mammalian mitochondrial ribosome at two different steps of the tRNA translocation reaction during translation elongation. Our structures explain the mechanism of tRNA and mRNA translocation on the mitoribosome, the regulation of mtEFG1 activity by the ribosomal GTPase-associated center, and the basis of decreased susceptibility of mtEFG1 to the commonly used antibiotic fusidic acid.

Structural analysis of 70S ribosomes by cross-linking/mass spectrometry reveals conformational plasticity

The ribosome is not only a highly complex molecular machine that translates the genetic information into proteins, but also an exceptional specimen for testing and optimizing cross-linking/mass spectrometry (XL-MS) workflows. Due to its high abundance, ribosomal proteins are frequently identified in proteome-wide XL-MS studies of cells or cell extracts. Here, we performed in-depth cross-linking of the E. coli ribosome using the amine-reactive cross-linker disuccinimidyl diacetic urea (DSAU). We analyzed 143 E. coli ribosomal structures, mapping a total of 10,771 intramolecular distances for 126 cross-link-pairs and 3,405 intermolecular distances for 97 protein pairs. Remarkably, 44% of intermolecular cross-links covered regions that have not been resolved in any high-resolution E. coli ribosome structure and point to a plasticity of cross-linked regions. We systematically characterized all cross-links and discovered flexible regions, conformational changes, and stoichiometric variations in bound ribosomal proteins, and ultimately remodeled 2,057 residues (15,794 atoms) in total. Our working model explains more than 95% of all cross-links, resulting in an optimized E. coli ribosome structure based on the cross-linking data obtained. Our study might serve as benchmark for conducting biochemical experiments on newly modeled protein regions, guided by XL-MS. Data are available via ProteomeXchange with identifier PXD018935.

Resolution-exchanged structural modeling and simulations jointly unravel that subunit rolling underlies the mechanism of programmed ribosomal frameshifting

MOTIVATION

Programmed ribosomal frameshifting (PRF) is widely used by viruses and bacteria to produce different proteins from a single mRNA template. How steric hindrance of a PRF-stimulatory mRNA structure transiently modifies the conformational dynamics of the ribosome, and thereby allows tRNA slippage, remains elusive.

RESULTS

Here, we leverage linear response theories and resolution-exchanged simulations to construct a structural/dynamics model that connects and rationalizes existing structural, single-molecule and mutagenesis data by resolution-exchanged structural modelling and simulations. Our combined theoretical techniques provide a temporal and spatial description of PRF with unprecedented mechanistic details. We discover that ribosomal unfolding of the PRF-stimulating pseudoknot exerts resistant forces on the mRNA entrance of the ribosome, and thereby drives 30S subunit rolling. Such motion distorts tRNAs, leads to tRNA slippage, and in turn serves as a delicate control of cis-element's unwinding forces over PRF.

AVAILABILITY AND IMPLEMENTATION

All the simulation scripts and computational implementations of our methods/analyses (including linear response theory) are included in the bioStructureM suite, provided through GitHub at https://github.com/Yuan-Yu/bioStructureM.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

An Efficient Timer and Sizer of Biomacromolecular Motions

Life ticks as fast as how proteins move. Computationally expensive molecular dynamics simulation has been the only theoretical tool to gauge the time and sizes of these motions, though barely to their slowest ends. Here, we convert a computationally cheap elastic network model (ENM) into a molecular timer and sizer to gauge the slowest functional motions of structured biomolecules. Quasi-harmonic analysis, fluctuation profile matching, and the Wiener-Khintchine theorem are used to define the "time periods," t, for anharmonic principal components (PCs), which are validated by nuclear magnetic resonance (NMR) order parameters. The PCs with their respective "time periods" are mapped to the eigenvalues (λ_ENM) of the corresponding ENM modes. Thus, the power laws t(ns) = 56.1λ_ENM^-1.6 and σ²(Ų) = 32.7λ_ENM^-3.0 can be established allowing the characterization of the timescales of NMR-resolved conformers, crystallographic anisotropic displacement parameters, and important ribosomal motions, as well as motional sizes of the latter.

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.

Bayesian-Estimated Hierarchical HMMs Enable Robust Analysis of Single-Molecule Kinetic Heterogeneity

Single-molecule kinetic experiments allow the reaction trajectories of individual biomolecules to be directly observed, eliminating the effects of population averaging and providing a powerful approach for elucidating the kinetic mechanisms of biomolecular processes. A major challenge to the analysis and interpretation of these experiments, however, is the kinetic heterogeneity that almost universally complicates the recorded single-molecule signal versus time trajectories (i.e., signal trajectories). Such heterogeneity manifests as changes and/or differences in the transition rates that are observed within individual signal trajectories or across a population of signal trajectories. Because characterizing kinetic heterogeneity can provide critical mechanistic information, we have developed a computational method that effectively and comprehensively enables such analysis. To this end, we have developed a computational algorithm and software program, hFRET, that uses the variational approximation for Bayesian inference to estimate the parameters of a hierarchical hidden Markov model, thereby enabling robust identification and characterization of kinetic heterogeneity. Using simulated signal trajectories, we demonstrate the ability of hFRET to accurately and precisely characterize kinetic heterogeneity. In addition, we use hFRET to analyze experimentally recorded signal trajectories reporting on the conformational dynamics of ribosomal pre-translocation (PRE) complexes. The results of our analyses demonstrate that PRE complexes exhibit kinetic heterogeneity, reveal the physical origins of this heterogeneity, and allow us to expand the current model of PRE complex dynamics. The methods described here can be applied to signal trajectories generated using any type of signal and can be easily extended to the analysis of signal trajectories exhibiting more complex kinetic behaviors. Moreover, variations of our approach can be easily developed to integrate kinetic data obtained from different experimental constructs and/or from molecular dynamics simulations of a biomolecule of interest.

The T-Box Riboswitch: tRNA as an Effector to Modulate Gene Regulation

The T-box riboswitch is a unique, RNA-based regulatory mechanism that modulates expression of a wide variety of amino acid-related genes, predominantly in Firmicutes. RNAs of this class selectively bind a specific cognate tRNA, utilizing recognition of the tRNA anticodon and other tRNA features. The riboswitch monitors the aminoacylation status of the tRNA to induce expression of the regulated downstream gene(s) at the level of transcription antitermination or derepression of translation initiation in response to reduced tRNA charging via stabilization of an antiterminator or antisequestrator. Recent biochemical and structural studies have revealed new features of tRNA recognition that extend beyond the initially identified Watson-Crick base-pairing of a codon-like sequence in the riboswitch with the tRNA anticodon, and residues in the antiterminator or antisequestrator with the tRNA acceptor end. These studies have revealed new tRNA contacts and new modes of riboswitch function and ligand recognition that expand our understanding of RNA-RNA recognition and the biological roles of tRNA.

Evolutionary shift toward protein-based architecture in trypanosomal mitochondrial ribosomes

Ribosomal RNA (rRNA) plays key functional and architectural roles in ribosomes. Using electron microscopy, we determined the atomic structure of a highly divergent ribosome found in mitochondria of Trypanosoma brucei, a unicellular parasite that causes sleeping sickness in humans. The trypanosomal mitoribosome features the smallest rRNAs and contains more proteins than all known ribosomes. The structure shows how the proteins have taken over the role of architectural scaffold from the rRNA: They form an autonomous outer shell that surrounds the entire particle and stabilizes and positions the functionally important regions of the rRNA. Our results also reveal the "minimal" set of conserved rRNA and protein components shared by all ribosomes that help us define the most essential functional elements.

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.

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

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

The ribosome moves: RNA mechanics and translocation

During protein synthesis, mRNA and tRNAs must be moved rapidly through the ribosome while maintaining the translational reading frame. This process is coupled to large- and small-scale conformational rearrangements in the ribosome, mainly in its rRNA. The free energy from peptide-bond formation and GTP hydrolysis is probably used to impose directionality on those movements. We propose that the free energy is coupled to two pawls, namely tRNA and EF-G, which enable two ratchet mechanisms to act separately and sequentially on the two ribosomal subunits.

Ribosome structural dynamics in translocation: yet another functional role for ribosomal RNA

Ribosomes are remarkable ribonucleoprotein complexes that are responsible for protein synthesis in all forms of life. They polymerize polypeptide chains programmed by nucleotide sequences in messenger RNA in a mechanism mediated by transfer RNA. One of the most challenging problems in the ribosome field is to understand the mechanism of coupled translocation of mRNA and tRNA during the elongation phase of protein synthesis. In recent years, the results of structural, biophysical and biochemical studies have provided extensive evidence that translocation is based on the structural dynamics of the ribosome itself. Detailed structural analysis has shown that ribosome dynamics, like aminoacyl-tRNA selection and catalysis of peptide bond formation, is made possible by the properties of ribosomal RNA.

Anisotropic Fluctuations in the Ribosome Determine tRNA Kinetics

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

Quantitative Characterization of Domain Motions in Molecular Machines

The work of molecular machines such as the ribosome is accompanied by conformational changes, often characterized by relative motions of their domains. The method we have developed seeks to quantify these motions in a general way, facilitating comparisons of results obtained by different researchers. Typically there are multiple snapshots of a structure in the form of pdb coordinates resulting from flexible fitting of low-resolution density maps, from X-ray crystallography, or from molecular dynamics simulation trajectories. Our objective is to characterize the motion of each domain as a coordinate transformation using moments of inertia tensor, a method we developed earlier. What has been missing until now are ancillary tools that make this task practical, general, and biologically informative. We have provided a comprehensive solution to this task with a set of tools implemented on the VMD platform. These tools address the need for reproducible segmentation of domains, and provide a generalized description of their motions using principal axes of inertia. Although this methodology has been specifically developed for studying ribosome motion, it is applicable to any molecular machine.

The Renormalization Group and Its Applications to Generating Coarse-Grained Models of Large Biological Molecular Systems

Understanding the dynamics of biomolecules is the key to understanding their biological activities. Computational methods ranging from all-atom molecular dynamics simulations to coarse-grained normal-mode analyses based on simplified elastic networks provide a general framework to studying these dynamics. Despite recent successes in studying very large systems with up to a 100,000,000 atoms, those methods are currently limited to studying small- to medium-sized molecular systems due to computational limitations. One solution to circumvent these limitations is to reduce the size of the system under study. In this paper, we argue that coarse-graining, the standard approach to such size reduction, must define a hierarchy of models of decreasing sizes that are consistent with each other, i.e., that each model contains the information of the dynamics of its predecessor. We propose a new method, Decimate, for generating such a hierarchy within the context of elastic networks for normal-mode analysis. This method is based on the concept of the renormalization group developed in statistical physics. We highlight the details of its implementation, with a special focus on its scalability to large systems of up to millions of atoms. We illustrate its application on two large systems, the capsid of a virus and the ribosome translation complex. We show that highly decimated representations of those systems, containing down to 1% of their original number of atoms, still capture qualitatively and quantitatively their dynamics. Decimate is available as an OpenSource resource.

Recurring RNA structural motifs underlie the mechanics of L1 stalk movement

The L1 stalk of the large ribosomal subunit undergoes large-scale movements coupled to the translocation of deacylated tRNA during protein synthesis. We use quantitative comparative structural analysis to localize the origins of L1 stalk movement and to understand its dynamic interactions with tRNA and other structural elements of the ribosome. Besides its stacking interactions with the tRNA elbow, stalk movement is directly linked to intersubunit rotation, rotation of the 30S head domain and contact of the acceptor arm of deacylated tRNA with helix 68 of 23S rRNA. Movement originates from pivoting at stacked non-canonical base pairs in a Family A three-way junction and bending in an internal G-U-rich zone. Use of these same motifs as hinge points to enable such dynamic events as rotation of the 30S subunit head domain and in flexing of the anticodon arm of tRNA suggests that they represent general strategies for movement of functional RNAs.

Translocation of the tRNA on the ribosome is associated with large-scale molecular movements of the ribosomal L1 stalk. Here the authors identify the key determinants that allow these dramatic movements, and suggest they represent general strategies used to enable large-scale motions in functional RNAs.

Virus-induced gene silencing of the RPC5-like subunit of RNA polymerase III caused pleiotropic effects in Nicotiana benthamiana

In eukaryotic cells, RNA polymerase III is highly conserved and transcribes housekeeping genes such as ribosomal 5S rRNA, tRNA and other small RNAs. The RPC5-like subunit is one of the 17 subunits forming RNAPIII and its exact functional roles in the transcription are poorly understood. In this work, we report that virus-induced gene silencing of transcripts encoding a putative RPC5-like subunit of the RNA Polymerase III in a model species Nicotiana benthamiana had pleiotropic effects, including but not limited to severe dwarfing appearance, chlorosis, nearly complete reduction of internodes and abnormal leaf shape. Using transcriptomic analysis, we identified genes and pathways affected by RPC5 silencing and thus presumably related to the cellular roles of the subunit as well as to the downstream cascade of reactions in response to partial loss of RNA Polymerase III function. Our results suggest that silencing of the RPC5L in N. benthamiana disrupted not only functions commonly associated with the core RNA Polymerase III transcripts, but also more diverse cellular processes, including responses to stress. We believe this is the first demonstration that activity of the RPC5 subunit is critical for proper functionality of RNA Polymerase III and normal plant development.

Computational Methodologies for Real-Space Structural Refinement of Large Macromolecular Complexes

The rise of the computer as a powerful tool for model building and refinement has revolutionized the field of structure determination for large biomolecular systems. Despite the wide availability of robust experimental methods capable of resolving structural details across a range of spatiotemporal resolutions, computational hybrid methods have the unique ability to integrate the diverse data from multimodal techniques such as X-ray crystallography and electron microscopy into consistent, fully atomistic structures. Here, commonly employed strategies for computational real-space structural refinement are reviewed, and their specific applications are illustrated for several large macromolecular complexes: ribosome, virus capsids, chemosensory array, and photosynthetic chromatophore. The increasingly important role of computational methods in large-scale structural refinement, along with current and future challenges, is discussed.

Advances in the molecular dynamics flexible fitting method for cryo-EM modeling

Molecular Dynamics Flexible Fitting (MDFF) is an established technique for fitting all-atom structures of molecules into corresponding cryo-electron microscopy (cryo-EM) densities. The practical application of MDFF is simple but requires a user to be aware of and take measures against a variety of possible challenges presented by each individual case. Some of these challenges arise from the complexity of a molecular structure or the limited quality of available structural models and densities to be interpreted, while others stem from the intricacies of MDFF itself. The current article serves as an overview of the strategies that have been developed since MDFF's inception to overcome common challenges and successfully perform MDFF simulations.

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

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

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

Nascent peptide assists the ribosome in recognizing chemically distinct small molecules

Regulation of gene expression in response to the changing environment is critical for cell survival. For instance, binding of macrolide antibiotics to the ribosome promote the translation arrest at the leader ORFs ermCL and ermBL necessary for inducing antibiotic resistance genes ermC and ermB. Cladinose-containing macrolides, like erythromycin (ERY), but not ketolides e.g., telithromycin (TEL), arrest translation of ermCL, while either ERY or TEL stall ermBL translation. How the ribosome distinguishes between chemically similar small molecules is unknown.

We show that single amino acid changes in the leader peptide switch the specificity of recognition of distinct molecules, triggering gene activation in response to only ERY, only TEL, to both antibiotics, or preventing stalling altogether. Thus, the ribosomal response to chemical signals can be modulated by minute changes in the nascent peptide, suggesting that protein sequences could have been optimized for rendering translation sensitive to environmental cues.

Functional Importance of Mobile Ribosomal Proteins

Although the dynamic motions and peptidyl transferase activity seem to be embedded in the rRNAs, the ribosome contains more than 50 ribosomal proteins (r-proteins), whose functions remain largely elusive. Also, the precise forms of some of these r-proteins, as being part of the ribosome, are not structurally solved due to their high flexibility, which hinders the efforts in their functional elucidation. Owing to recent advances in cryo-electron microscopy, single-molecule techniques, and theoretical modeling, much has been learned about the dynamics of these r-proteins. Surprisingly, allosteric regulations have been found in between spatially separated components as distant as those in the opposite sides of the ribosome. Here, we focus on the functional roles and intricate regulations of the mobile L1 and L12 stalks and L9 and S1 proteins. Conformational flexibility also enables versatile functions for r-proteins beyond translation. The arrangement of r-proteins may be under evolutionary pressure that fine-tunes mass distributions for optimal structural dynamics and catalytic activity of the ribosome.

Universal and domain-specific sequences in 23S-28S ribosomal RNA identified by computational phylogenetics

Comparative analysis of ribosomal RNA (rRNA) sequences has elucidated phylogenetic relationships. However, this powerful approach has not been fully exploited to address ribosome function. Here we identify stretches of evolutionarily conserved sequences, which correspond with regions of high functional importance. For this, we developed a structurally aligned database, FLORA (full-length organismal rRNA alignment) to identify highly conserved nucleotide elements (CNEs) in 23S-28S rRNA from each phylogenetic domain (Eukarya, Bacteria, and Archaea). Universal CNEs (uCNEs) are conserved in sequence and structural position in all three domains. Those in regions known to be essential for translation validate our approach. Importantly, some uCNEs reside in areas of unknown function, thus identifying novel sequences of likely great importance. In contrast to uCNEs, domain-specific CNEs (dsCNEs) are conserved in just one phylogenetic domain. This is the first report of conserved sequence elements in rRNA that are domain-specific; they are largely a eukaryotic phenomenon. The locations of the eukaryotic dsCNEs within the structure of the ribosome suggest they may function in nascent polypeptide transit through the ribosome tunnel and in tRNA exit from the ribosome. Our findings provide insights and a resource for ribosome function studies.

The ribosome's energy landscape: Recent insights from computation

The ever-increasing capacity of computing resources has extended ribosome calculations from the study of small-scale fluctuations to large-scale barrier-crossing processes. As the field of computational/theoretical biophysics shifts focus to large-scale conformational transitions, there is a growing need for a systematic framework to interpret and analyze ribosome dynamics. To this end, energy landscape principles, largely developed for the study of biomolecular folding, have proven to be invaluable. These tools not only provide a foundation for describing simulations but can be used to reconcile experimental results, as well. In this review, I will discuss recent efforts to employ computational methods to reveal the characteristics of the ribosome's landscape and how these studies can help guide a new generation of experiments that more closely probe the underlying energetics. As a result of these investigations, general principles about ribosome function are beginning to emerge, including that: (1) small-scale fluctuations are the result of structure, rather than detailed energetics, (2) molecular flexibility leads to entropically favored rearrangements, and (3) tRNA dynamics may be accurately described as diffusive movement across an energy landscape.

Dynamic contact network between ribosomal subunits enables rapid large-scale rotation during spontaneous translocation

During ribosomal translation, the two ribosomal subunits remain associated through intersubunit bridges, despite rapid large-scale intersubunit rotation. The absence of large barriers hindering rotation is a prerequisite for rapid rotation. Here, we investigate how such a flat free-energy landscape is achieved, in particular considering the large shifts the bridges undergo at the periphery. The dynamics and energetics of the intersubunit contact network are studied using molecular dynamics simulations of the prokaryotic ribosome in intermediate states of spontaneous translocation. Based on observed occupancies of intersubunit contacts, residues were grouped into clusters. In addition to the central contact clusters, peripheral clusters were found to maintain strong steady interactions by changing contacts in the course of rotation. The peripheral B1 bridges are stabilized by a changing contact pattern of charged residues that adapts to the rotational state. In contrast, steady strong interactions of the B4 bridge are ensured by the flexible helix H34 following the movement of protein S15. The tRNAs which span the subunits contribute to the intersubunit binding enthalpy to an almost constant degree, despite their different positions in the ribosome. These mechanisms keep the intersubunit interaction strong and steady during rotation, thereby preventing dissociation and enabling rapid rotation.

Designed Regular Tetragon-Shaped RNA-Protein Complexes with Ribosomal Protein L1 for Bionanotechnology and Synthetic Biology

RNA nanotechnology has been established by employing the molecular architecture of RNA structural motifs. Here, we report two designed RNA-protein complexes (RNPs) composed of ribosomal protein L1 (RPL1) and its RNA-binding motif that are square-shaped nano-objects. The formation and the shape of the objects were confirmed by gel electrophoresis analysis and atomic force microscopy, respectively. Any protein can be attached to the RNA via a fusion protein with RPL1, indicating that it can be used as a scaffold for loading a variety of functional proteins or for building higher-order structures. In summary, the RNP object will serve as a useful tool in the fields of bionanotechnology and synthetic biology. Moreover, the RNP interaction enhances the RNA stability against nucleases, rendering these complexes stable in cells.

Using Molecular Simulation to Model High-Resolution Cryo-EM Reconstructions

An explosion of new data from high-resolution cryo-electron microscopy (cryo-EM) studies has produced a large number of data sets for many species of ribosomes in various functional states over the past few years. While many methods exist to produce structural models for lower resolution cryo-EM reconstructions, high-resolution reconstructions are often modeled using crystallographic techniques and extensive manual intervention. Here, we present an automated fitting technique for high-resolution cryo-EM data sets that produces all-atom models highly consistent with the EM density. Using a molecular dynamics approach, atomic positions are optimized with a potential that includes the cross-correlation coefficient between the structural model and the cryo-EM electron density, as well as a biasing potential preserving the stereochemistry and secondary structure of the biomolecule. Specifically, we use a hybrid structure-based/ab initio molecular dynamics potential to extend molecular dynamics fitting. In addition, we find that simulated annealing integration, as opposed to straightforward molecular dynamics integration, significantly improves performance. We obtain atomistic models of the human ribosome consistent with high-resolution cryo-EM reconstructions of the human ribosome. Automated methods such as these have the potential to produce atomistic models for a large number of ribosome complexes simultaneously that can be subsequently refined manually.

Quantitative Connection between Ensemble Thermodynamics and Single-Molecule Kinetics: A Case Study Using Cryogenic Electron Microscopy and Single-Molecule Fluorescence Resonance Energy Transfer Investigations of the Ribosome

At equilibrium, thermodynamic and kinetic information can be extracted from biomolecular energy landscapes by many techniques. However, while static, ensemble techniques yield thermodynamic data, often only dynamic, single-molecule techniques can yield the kinetic data that describe transition-state energy barriers. Here we present a generalized framework based upon dwell-time distributions that can be used to connect such static, ensemble techniques with dynamic, single-molecule techniques, and thus characterize energy landscapes to greater resolutions. We demonstrate the utility of this framework by applying it to cryogenic electron microscopy (cryo-EM) and single-molecule fluorescence resonance energy transfer (smFRET) studies of the bacterial ribosomal pre-translocation complex. Among other benefits, application of this framework to these data explains why two transient, intermediate conformations of the pre-translocation complex, which are observed in a cryo-EM study, may not be observed in several smFRET studies.

Studying the properties of domain I of the ribosomal protein l1: incorporation into ribosome and regulation of the l1 operon expression

L1 is a conserved protein of the large ribosomal subunit. This protein binds strongly to the specific region of the high molecular weight rRNA of the large ribosomal subunit, thus forming a conserved flexible structural element--the L1 stalk. L1 protein also regulates translation of the operon that comprises its own gene. Crystallographic data suggest that L1 interacts with RNA mainly by means of its domain I. We show here for the first time that the isolated domain I of the bacterial protein L1 of Thermus thermophilus and Escherichia coli is able to incorporate in vivo into the E. coli ribosome. Furthermore, domain I of T. thermophilus L1 can regulate expression of the L1 gene operon of Archaea in the coupled transcription-translation system in vitro, as well as the intact protein. We have identified the structural elements of domain I of the L1 protein that may be responsible for its regulatory properties.

Can We Execute Stable Microsecond-Scale Atomistic Simulations of Protein-RNA Complexes?

We report over 30 μs of unrestrained molecular dynamics simulations of six protein-RNA complexes in explicit solvent. We utilize the AMBER ff99bsc0χ(OL3) RNA force field combined with the ff99SB protein force field and its more recent ff12SB version with reparametrized side-chain dihedrals. The simulations show variable behavior, ranging from systems that are essentially stable to systems with progressive deviations from the experimental structure, which we could not stabilize anywhere close to the starting experimental structure. For some systems, microsecond-scale simulations are necessary to achieve stabilization after initial sizable structural perturbations. The results show that simulations of protein-RNA complexes are challenging and every system should be treated individually. The simulations are affected by numerous factors, including properties of the starting structures (the initially high force field potential energy, resolution limits, conformational averaging, crystal packing, etc.), force field imbalances, and real flexibility of the studied systems. These factors, and thus the simulation behavior, differ from system to system. The structural stability of simulated systems does not correlate with the size of buried interaction surface or experimentally determined binding affinities but reflects the type of protein-RNA recognition. Protein-RNA interfaces involving shape-specific recognition of RNA are more stable than those relying on sequence-specific RNA recognition. The differences between the protein force fields are considerably smaller than the uncertainties caused by sampling and starting structures. The ff12SB improves description of the tyrosine side-chain group, which eliminates some problems associated with tyrosine dynamics.

Protein-RNA affinity of ribosomal protein L1 mutants does not correlate with the number of intermolecular interactions

Ribosomal protein L1, as part of the L1 stalk of the 50S ribosomal subunit, is implicated in directing tRNA movement through the ribosome during translocation. High-resolution crystal structures of four mutants (T217V, T217A, M218L and G219V) of the ribosomal protein L1 from Thermus thermophilus (TthL1) in complex with a specific 80 nt fragment of 23S rRNA and the structures of two of these mutants (T217V and G219V) in the RNA-unbound form are reported in this work. All mutations are located in the highly conserved triad Thr-Met-Gly, which is responsible for about 17% of all protein-RNA hydrogen bonds and 50% of solvent-inaccessible intermolecular hydrogen bonds. In the mutated proteins without bound RNA the RNA-binding regions show substantial conformational changes. On the other hand, in the complexes with RNA the structures of the RNA-binding surfaces in all studied mutants are very similar to the structure of the wild-type protein in complex with RNA. This shows that formation of the RNA complexes restores the distorted surfaces of the mutant proteins to a conformation characteristic of the wild-type protein complex. Domain I of the mutated TthL1 and helix 77 of 23S rRNA form a rigid body identical to that found in the complex of wild-type TthL1 with RNA, suggesting that the observed relative orientation is conserved and is probably important for ribosome function. Analysis of the complex structures and the kinetic data show that the number of intermolecular contacts and hydrogen bonds in the RNA-protein contact area does not correlate with the affinity of the protein for RNA and cannot be used as a measure of affinity.

Flexibility in the insulin receptor ectodomain enables docking of insulin in crystallographic conformation observed in a hormone-bound microreceptor

Insulin binding to the insulin receptor (IR) is the first key step in initiating downstream signaling cascades for glucose homeostasis in higher organisms. The molecular details of insulin recognition by IR are not yet completely understood, but a picture of hormone/receptor interactions at one of the epitopes (Site 1) is beginning to emerge from recent structural evidence. However, insulin-bound structures of truncated IR suggest that crystallographic conformation of insulin cannot be accommodated in the full IR ectodomain due to steric overlap of insulin with the first two type III fibronectin domains (F1 and F2), which are contributed to the insulin binding-pocket by the second subunit in the IR homodimer. A conformational change in the F1-F2 pair has thus been suggested. In this work, we present an all-atom structural model of complex of insulin and the IR ectodomain, where no structural overlap of insulin with the receptor domains (F1 and F2) is observed. This structural model was arrived at by flexibly fitting parts of our earlier insulin/IR all-atom model into the simulated density maps of crystallized constructs combined with conformational sampling from apo-IR solution conformations. Importantly, our experimentally-consistent model helps rationalize yet unresolved Site.

xMDFF: molecular dynamics flexible fitting of low-resolution X-ray structures

X-ray crystallography remains the most dominant method for solving atomic structures. However, for relatively large systems, the availability of only medium-to-low-resolution diffraction data often limits the determination of all-atom details. A new molecular dynamics flexible fitting (MDFF)-based approach, xMDFF, for determining structures from such low-resolution crystallographic data is reported. xMDFF employs a real-space refinement scheme that flexibly fits atomic models into an iteratively updating electron-density map. It addresses significant large-scale deformations of the initial model to fit the low-resolution density, as tested with synthetic low-resolution maps of D-ribose-binding protein. xMDFF has been successfully applied to re-refine six low-resolution protein structures of varying sizes that had already been submitted to the Protein Data Bank. Finally, via systematic refinement of a series of data from 3.6 to 7 Å resolution, xMDFF refinements together with electrophysiology experiments were used to validate the first all-atom structure of the voltage-sensing protein Ci-VSP.

Synchronous tRNA movements during translocation on the ribosome are orchestrated by elongation factor G and GTP hydrolysis

The translocation of tRNAs through the ribosome proceeds through numerous small steps in which tRNAs gradually shift their positions on the small and large ribosomal subunits. The most urgent questions are: (i) whether these intermediates are important; (ii) how the ribosomal translocase, the GTPase elongation factor G (EF-G), promotes directed movement; and (iii) how the energy of GTP hydrolysis is coupled to movement. In the light of recent advances in biophysical and structural studies, we argue that intermediate states of translocation are snapshots of dynamic fluctuations that guide the movement. In contrast to current models of stepwise translocation, kinetic evidence shows that the tRNAs move synchronously on the two ribosomal subunits in a rapid reaction orchestrated by EF-G and GTP hydrolysis. EF-G combines the energy regimes of a GTPase and a motor protein and facilitates tRNA movement by a combination of directed Brownian ratchet and power stroke mechanisms.

New trends in aminoglycosides use

Despite their inherent toxicity and the acquired bacterial resistance that continuously threaten their long-term clinical use, aminoglycosides (AGs) still remain valuable components of the antibiotic armamentarium. Recent literature shows that the AGs' role has been further expanded as multi-tasking players in different areas of study. This review aims at presenting some of the new trends observed in the use of AGs in the past decade, along with the current understanding of their mechanisms of action in various bacterial and eukaryotic cellular processes.