Atomic model of the Thermus thermophilus 70S ribosome developed in silico

Molecular basis of the pleiotropic effects by the antibiotic amikacin on the ribosome

Aminoglycosides are a class of antibiotics that bind to ribosomal RNA and exert pleiotropic effects on ribosome function. Amikacin, the semisynthetic derivative of kanamycin, is commonly used for treating severe infections with multidrug-resistant, aerobic Gram-negative bacteria. Amikacin carries the 4-amino-2-hydroxy butyrate (AHB) moiety at the N1 amino group of the central 2-deoxystreptamine (2-DOS) ring, which may confer amikacin a unique ribosome inhibition profile. Here we use in vitro fast kinetics combined with X-ray crystallography and cryo-EM to dissect the mechanisms of ribosome inhibition by amikacin and the parent compound, kanamycin. Amikacin interferes with tRNA translocation, release factor-mediated peptidyl-tRNA hydrolysis, and ribosome recycling, traits attributed to the additional interactions amikacin makes with the decoding center. The binding site in the large ribosomal subunit proximal to the 3'-end of tRNA in the peptidyl (P) site lays the groundwork for rational design of amikacin derivatives with improved antibacterial properties.

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.

Getting to the bottom of lncRNA mechanism: structure-function relationships

While long non-coding RNAs are known to play key roles in disease and development, relatively few structural studies have been performed for this important class of RNAs. Here, we review functional studies of long non-coding RNAs and expose the need for high-resolution 3-D structural studies, discussing the roles of long non-coding RNAs in the cell and how structure–function relationships might be used to elucidate further understanding. We then describe structural studies of other classes of RNAs using chemical probing, nuclear magnetic resonance, small-angle X-ray scattering, X-ray crystallography, and cryogenic electron microscopy (cryo-EM). Next, we review early structural studies of long non-coding RNAs to date and describe the way forward for the structural biology of long non-coding RNAs in terms of cryo-EM.

Comparative Analysis of Aggregation of Thermus thermophilus Ribosomal Protein bS1 and Its Stable Fragment

Functionally important multidomain bacterial protein bS1 is the largest ribosomal protein of subunit 30S. It interacts with both mRNA and proteins and is prone to aggregation, although this process has not been studied in detail. Here, we obtained bacterial strains overproducing ribosomal bS1 protein from Thermus thermophilus and its stable fragment bS1(49) and purified these proteins. Using fluorescence spectroscopy, dynamic light scattering, and high-performance liquid chromatography combined with mass spectrometric analysis of products of protein limited proteolysis, we demonstrated that disordered regions at the N- and C-termini of bS1 can play a key role in the aggregation of this protein. The truncated fragment bS1(49) was less prone to aggregation compared to the full-size bS1. The revealed properties of the studied proteins can be used to obtain protein crystals for elucidating the structure of the bS1 stable fragment.

Modeling the Influenza A NP-vRNA-Polymerase Complex in Atomic Detail

Seasonal flu is an acute respiratory disease that exacts a massive toll on human populations, healthcare systems and economies. The disease is caused by an enveloped Influenza virus containing eight ribonucleoprotein (RNP) complexes. Each RNP incorporates multiple copies of nucleoprotein (NP), a fragment of the viral genome (vRNA), and a viral RNA-dependent RNA polymerase (POL), and is responsible for packaging the viral genome and performing critical functions including replication and transcription. A complete model of an Influenza RNP in atomic detail can elucidate the structural basis for viral genome functions, and identify potential targets for viral therapeutics. In this work we construct a model of a complete Influenza A RNP complex in atomic detail using multiple sources of structural and sequence information and a series of homology-modeling techniques, including a motif-matching fragment assembly method. Our final model provides a rationale for experimentally-observed changes to viral polymerase activity in numerous mutational assays. Further, our model reveals specific interactions between the three primary structural components of the RNP, including potential targets for blocking POL-binding to the NP-vRNA complex. The methods developed in this work open the possibility of elucidating other functionally-relevant atomic-scale interactions in additional RNP structures and other biomolecular complexes.

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.

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.

Analysis of ribosomal inter-subunit sites as targets for complementary oligonucleotides

The bacterial ribosome has many functional ribosomal RNA (rRNA) sites. We have computationally analyzed the rRNA regions involved in the interactions between the 30S and 50S subunits. Various properties of rRNA such as solvent accessibility, opening energy, hydrogen bonding pattern, van der Waals energy, thermodynamic stability were determined. Based on these properties we selected rRNA targets for hybridization with complementary 2'-O-methyl oligoribonucleotides (2'-OMe RNAs). Further, the inhibition efficiencies of the designed ribosome-interfering 2'-OMe RNAs were tested using a β-galactosidase assay in a translation system based on the E. coli extract. Several of the oligonucleotides displayed IC₅₀ values below 1 μM, which were in a similar range as those determined for known ribosome inhibitors, tetracycline and pactamycin. The calculated opening and van der Waals stacking energies of the rRNA targets correlated best with the inhibitory efficiencies of 2'-OMe RNAs. Moreover, the binding affinities of several oligonucleotides to both 70S ribosomes and isolated 30S and 50S subunits were measured using a double-filter retention assay. Further, we applied heat-shock chemical transformation to introduce 2'-OMe RNAs to E. coli cells and verify inhibition of bacterial growth. We observed high correlation between IC₅₀ in the cell-free extract and bacterial growth inhibition. Overall, the results suggest that the computational analysis of potential rRNA targets within the conformationally dynamic regions of inter-subunit bridges can help design efficient antisense oligomers to probe the ribosome function.

Ribosome Mechanics Informs about Mechanism

The essential aspects of the ribosome's mechanism can be extracted from coarse-grained simulations, including the ratchet motion, the movement together of critical bases at the decoding center, and movements of the peptide tunnel lining that assist in the expulsion of the synthesized peptide. Because of its large size, coarse graining helps to simplify and to aid in the understanding of its mechanism. Results presented here utilize coarse-grained elastic network modeling to extract the dynamics, and both RNAs and proteins are coarse grained. We review our previous results, showing the well-known ratchet motions and the motions in the peptide tunnel and in the mRNA tunnel. The motions of the lining of the peptide tunnel appear to assist in the expulsion of the growing peptide chain, and clamps at the ends of the mRNA tunnel with three proteins ensure that the mRNA is held tightly during decoding and essential for the helicase activity at the entrance. The entry clamp may also assist in base recognition to ensure proper selection of the incoming tRNA. The overall precision of the ribosome machine-like motions is remarkable.

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.

Integrating molecular dynamics simulations with chemical probing experiments using SHAPE-FIT

Integration and calibration of molecular dynamics simulations with experimental data remain a challenging endeavor. We have developed a novel method to integrate chemical probing experiments with molecular simulations of RNA molecules by using a native structure-based model. Selective 2'-hydroxyl acylation by primer extension (SHAPE) characterizes the mobility of each residue in the RNA. Our method, SHAPE-FIT, automatically optimizes the potential parameters of the force field according to measured reactivities from SHAPE. The optimized parameter set allows simulations of dynamics highly consistent with SHAPE probing experiments. Such atomistic simulations, thoroughly grounded in experiment, can open a new window on RNA structure-function relations.

Cryo-EM of ribosomal 80S complexes with termination factors reveals the translocated cricket paralysis virus IRES

The cricket paralysis virus (CrPV) uses an internal ribosomal entry site (IRES) to hijack the ribosome. In a remarkable RNA-based mechanism involving neither initiation factor nor initiator tRNA, the CrPV IRES jumpstarts translation in the elongation phase from the ribosomal A site. Here, we present cryoelectron microscopy (cryo-EM) maps of 80S⋅CrPV-STOP ⋅ eRF1 ⋅ eRF3 ⋅ GMPPNP and 80S⋅CrPV-STOP ⋅ eRF1 complexes, revealing a previously unseen binding state of the IRES and directly rationalizing that an eEF2-dependent translocation of the IRES is required to allow the first A-site occupation. During this unusual translocation event, the IRES undergoes a pronounced conformational change to a more stretched conformation. At the same time, our structural analysis provides information about the binding modes of eRF1 ⋅ eRF3 ⋅ GMPPNP and eRF1 in a minimal system. It shows that neither eRF3 nor ABCE1 are required for the active conformation of eRF1 at the intersection between eukaryotic termination and recycling.

Regulation of the mammalian elongation cycle by subunit rolling: a eukaryotic-specific ribosome rearrangement

The extent to which bacterial ribosomes and the significantly larger eukaryotic ribosomes share the same mechanisms of ribosomal elongation is unknown. Here, we present subnanometer resolution cryoelectron microscopy maps of the mammalian 80S ribosome in the posttranslocational state and in complex with the eukaryotic eEF1A⋅Val-tRNA⋅GMPPNP ternary complex, revealing significant differences in the elongation mechanism between bacteria and mammals. Surprisingly, and in contrast to bacterial ribosomes, a rotation of the small subunit around its long axis and orthogonal to the well-known intersubunit rotation distinguishes the posttranslocational state from the classical pretranslocational state ribosome. We term this motion "subunit rolling." Correspondingly, a mammalian decoding complex visualized in substates before and after codon recognition reveals structural distinctions from the bacterial system. These findings suggest how codon recognition leads to GTPase activation in the mammalian system and demonstrate that in mammalia subunit rolling occurs during tRNA selection.

Fitness costs of rifampicin resistance in Mycobacterium tuberculosis are amplified under conditions of nutrient starvation and compensated by mutation in the β' subunit of RNA polymerase

Rifampicin resistance, a defining attribute of multidrug-resistant tuberculosis, is conferred by mutations in the β subunit of RNA polymerase. Sequencing of rifampicin-resistant (RIF-R) clinical isolates of Mycobacterium tuberculosis revealed, in addition to RIF-R mutations, enrichment of potential compensatory mutations around the double-psi β-barrel domain of the β' subunit comprising the catalytic site and the exit tunnel for newly synthesized RNA. Sequential introduction of the resistance allele followed by the compensatory allele in isogenic Mycobacterium smegmatis showed that these mutations respectively caused and compensated a starvation enhanced growth defect by altering RNA polymerase activity. While specific combinations of resistance and compensatory alleles converged in divergent lineages, other combinations recurred among related isolates suggesting transmission of compensated RIF-R strains. These findings suggest nutrient poor growth conditions impose larger selective pressure on RIF-R organisms that results in the selection of compensatory mutations in a domain involved in catalysis and starvation control of RNA polymerase transcription.

Tackling structures of long noncoding RNAs

RNAs are important catalytic machines and regulators at every level of gene expression. A new class of RNAs has emerged called long non-coding RNAs, providing new insights into evolution, development and disease. Long non-coding RNAs (lncRNAs) predominantly found in higher eukaryotes, have been implicated in the regulation of transcription factors, chromatin-remodeling, hormone receptors and many other processes. The structural versatility of RNA allows it to perform various functions, ranging from precise protein recognition to catalysis and metabolite sensing. While major housekeeping RNA molecules have long been the focus of structural studies, lncRNAs remain the least characterized class, both structurally and functionally. Here, we review common methodologies used to tackle RNA structure, emphasizing their potential application to lncRNAs. When considering the complexity of lncRNAs and lack of knowledge of their structure, chemical probing appears to be an indispensable tool, with few restrictions in terms of size, quantity and heterogeneity of the RNA molecule. Probing is not constrained to in vitro analysis and can be adapted to high-throughput sequencing platforms. Significant efforts have been applied to develop new in vivo chemical probing reagents, new library construction protocols for sequencing platforms and improved RNA prediction software based on the experimental evidence.

Computational studies of molecular machines: the ribosome

The past decade has produced an avalanche of experimental data on the structure and dynamics of the ribosome. Groundbreaking studies in structural biology and kinetics have placed important constraints on ribosome structural dynamics. However, a gulf remains between static structures and time dependent data. In particular, X-ray crystallography and cryo-EM studies produce static models of the ribosome in various states, but lack dynamic information. Single molecule studies produce information on the rates of transitions between these states but do not have high-resolution spatial information. Computational studies have aided in bridging this gap by providing atomic resolution simulations of structural fluctuations and transitions between configurations.

Key intermolecular interactions in the E. coli 70S ribosome revealed by coarse-grained analysis

The ribosome is a very large complex that consists of many RNA and protein molecules and plays a central role in protein biosynthesis in all organisms. Extensive interactions between different molecules are critical to ribosomal functional dynamics. In this work, intermolecular interactions in the Escherichia coli 70S ribosome are investigated by coarse-grained (CG) analysis. CG models are defined to preserve dynamic domains in RNAs and proteins and to capture functional motions in the ribosome, and then the CG sites are connected by harmonic springs, and spring constants are obtained by matching the computed fluctuations to those of an all-atom molecular dynamics (MD) simulation. Those spring constants indicate how strong the interactions are between the ribosomal components, and they are in good agreement with various experimental data. Nearly all the bridges between the small and large ribosomal subunits are indicated by CG interactions with large spring constants. The head of the small subunit is very mobile because it has minimal CG interactions with the rest of the subunit; however, a large number of small subunit proteins bind to maintain the internal structure of the head. The results show a clear connection between the intermolecular interactions and the structural and functional properties of the ribosome because of the reduced complexity in domain-based CG models. The present approach also provides a useful strategy to map interactions between molecules within large biomolecular complexes since it is not straightforward to investigate these by either atomistic MD simulations or residue-based elastic network models.

ModeRNA: a tool for comparative modeling of RNA 3D structure

RNA is a large group of functionally important biomacromolecules. In striking analogy to proteins, the function of RNA depends on its structure and dynamics, which in turn is encoded in the linear sequence. However, while there are numerous methods for computational prediction of protein three-dimensional (3D) structure from sequence, with comparative modeling being the most reliable approach, there are very few such methods for RNA. Here, we present ModeRNA, a software tool for comparative modeling of RNA 3D structures. As an input, ModeRNA requires a 3D structure of a template RNA molecule, and a sequence alignment between the target to be modeled and the template. It must be emphasized that a good alignment is required for successful modeling, and for large and complex RNA molecules the development of a good alignment usually requires manual adjustments of the input data based on previous expertise of the respective RNA family. ModeRNA can model post-transcriptional modifications, a functionally important feature analogous to post-translational modifications in proteins. ModeRNA can also model DNA structures or use them as templates. It is equipped with many functions for merging fragments of different nucleic acid structures into a single model and analyzing their geometry. Windows and UNIX implementations of ModeRNA with comprehensive documentation and a tutorial are freely available.

Simulating activity of the bacterial ribosome

Computational modeling studies that investigate activity of the bacterial ribosome were reviewed. Computational approaches became possible with the availability of three-dimensional atomic resolution structures of the ribosomal subunits. However, due to the enormous size of the system, theoretical efforts to study the ribosome are few and challenging. For example, to extend the simulation timescales to biologically relevant ones, often, reduced models that require tedious parameterizations need to be applied. To that end, modeling of the ribosome focused on its internal dynamics, electrostatic properties, inhibition by antibiotics, polypeptide folding in the ribosome tunnel and assembly mechanisms driving the formation of the small ribosomal subunit.

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

The ribosome is a massive ribonucleoprotein complex ( approximately 2.4 MDa) that utilizes large-scale structural fluctuations to produce unidirectional protein synthesis. Accommodation is a key conformational change during transfer RNA (tRNA) selection that allows movement of tRNA into the ribosome. Here, we address the structure-function relationship that governs accommodation using all-atom molecular simulations and single-molecule fluorescence resonance energy transfer (smFRET). Simulations that employ an all-atom, structure-based (Gō-like) model illuminate the interplay between configurational entropy and effective enthalpy during the accommodation process. This delicate balance leads to spontaneous reversible accommodation attempts, which are corroborated by smFRET measurements. The dynamics about the endpoints of accommodation (the A/T and A/A conformations) obtained from structure-based simulations are validated by multiple 100-200 ns explicit-solvent simulations (3.2 million atoms for a cumulative 1.4 micros), and previous crystallographic analysis. We find that the configurational entropy of the 3'-CCA end of aminoacyl-tRNA resists accommodation, leading to a multistep accommodation process that encompasses a distribution of parallel pathways. The calculated mechanism is robust across simulation methods and protocols, suggesting that the structure of the accommodation corridor imposes stringent limitations on the accessible pathways. The identified mechanism and observed parallel pathways establish an atomistic framework for interpreting a large body of biochemical data and demonstrate that conformational changes during translation occur through a stochastic trial-and-error process, rather than in concerted lock-step motions.

Collective dynamics of the ribosomal tunnel revealed by elastic network modeling

The collective dynamics of the nascent polypeptide exit tunnel are investigated with the computationally efficient elastic network model using normal mode analysis. The calculated normal modes are considered individually and in linear combinations with different coefficients mimicking the phase angles between modes, in order to follow the mechanistic motions of tunnel wall residues. The low frequency fluctuations indicate three distinct regions along the tunnel-the entrance, the neck, and the exit-each having distinctly different domain motions. Generally, the lining of the entrance region moves in the exit direction, with the exit region having significantly larger motions, but in a perpendicular direction, whereas the confined neck region has rotational motions. Especially the universally conserved extensions of ribosomal proteins L4 and L22 located at the narrowest and mechanistically strategic region of tunnel undergo generally anti- or non-correlated motions, which may have an important role in nascent polypeptide gating mechanism. These motions appear to be sufficiently robust so as to be unaffected by the presence of a peptide chain in the tunnel.

Path of nascent polypeptide in exit tunnel revealed by molecular dynamics simulation of ribosome

Molecular dynamics simulations were carried out on Thermus thermophilus 70S ribosome with and without a nascent polypeptide inside the exit tunnel. Modeling of the polypeptide in the tunnel revealed two possible paths: one over Arg92 of L22 and one under (from the viewpoint of 50S on top of 30S). A strong interaction between L4 and Arg92 was observed without the polypeptide and when it passed over Arg92. However, when the polypeptide passed under, Arg92 repositioned to interact with Ade2059 of 23S rRNA. Using steered molecular dynamics the polypeptide could be pulled through the L4-L22 constriction when situated under Arg92, but did not move when over. These results suggest that the tunnel is closed by the Arg92-L4 interaction before elongation of the polypeptide and the tunnel leads the entering polypeptide from the peptidyl transferase center to the passage under Arg92, causing Arg92 to switch to an open position. It is possible, therefore, that Arg92 plays the role of a gate, opening and closing the tunnel at L4-L22. There is some disagreement over whether the tunnel is dynamic or rigid. At least within the timescale of our simulations conformational analysis showed that global motions mainly involve relative movement of the 50S and 30S subunits and seem not to affect the conformation of the tunnel.

Cross-crystal averaging reveals that the structure of the peptidyl-transferase center is the same in the 70S ribosome and the 50S subunit

Recently, two crystal structures of the Thermus thermophilus 70S ribosome in the same functional state were determined at 2.8 and 3.7 A resolution but were different throughout. The most functionally significant structural differences are in the conformation of the peptidyl-transferase center (PTC) and the interface between the PTC and the CCA end of the P-site tRNA. Likewise, the 3.7 A PTC differed from the functionally equivalent structure of the Haloarcula marismortui 50S subunit. To ascertain whether the 3.7 A model does indeed differ from the other two, we performed cross-crystal averaging of the two 70S data sets. The unbiased maps suggest that the conformation of the PTC-CCA in the two 70S crystal forms is identical to that of the 2.8 A 70S model as well as that of the H. marismortui 50S subunit. We conclude that the structure of the PTC is the same in the functionally equivalent 70S ribosome and the 50S subunit.

High performance computing in biology: multimillion atom simulations of nanoscale systems

Computational methods have been used in biology for sequence analysis (bioinformatics), all-atom simulation (molecular dynamics and quantum calculations), and more recently for modeling biological networks (systems biology). Of these three techniques, all-atom simulation is currently the most computationally demanding, in terms of compute load, communication speed, and memory load. Breakthroughs in electrostatic force calculation and dynamic load balancing have enabled molecular dynamics simulations of large biomolecular complexes. Here, we report simulation results for the ribosome, using approximately 2.64 million atoms, the largest all-atom biomolecular simulation published to date. Several other nano-scale systems with different numbers of atoms were studied to measure the performance of the NAMD molecular dynamics simulation program on the Los Alamos National Laboratory Q Machine. We demonstrate that multimillion atom systems represent a 'sweet spot' for the NAMD code on large supercomputers. NAMD displays an unprecedented 85% parallel scaling efficiency for the ribosome system on 1024 CPUs. We also review recent targeted molecular dynamics simulations of the ribosome that prove useful for studying conformational changes of this large biomolecular complex in atomic detail.

Crystal Structure of a 70S Ribosome-tRNA Complex Reveals Functional Interactions and Rearrangements

Our understanding of the mechanism of protein synthesis has undergone rapid progress in recent years as a result of low-resolution X-ray and cryo-EM structures of ribosome functional complexes and high-resolution structures of ribosomal subunits and vacant ribosomes. Here, we present the crystal structure of the Thermus thermophilus 70S ribosome containing a model mRNA and two tRNAs at 3.7 A resolution. Many structural details of the interactions between the ribosome, tRNA, and mRNA in the P and E sites and the ways in which tRNA structure is distorted by its interactions with the ribosome are seen. Differences between the conformations of vacant and tRNA-bound 70S ribosomes suggest an induced fit of the ribosome structure in response to tRNA binding, including significant changes in the peptidyl-transferase catalytic site.

Structure, dynamics, and elasticity of free 16s rRNA helix 44 studied by molecular dynamics simulations

Molecular dynamics (MD) simulations were employed to investigate the structure, dynamics, and local base-pair step deformability of the free 16S ribosomal helix 44 from Thermus thermophilus and of a canonical A-RNA double helix. While helix 44 is bent in the crystal structure of the small ribosomal subunit, the simulated helix 44 is intrinsically straight. It shows, however, substantial instantaneous bends that are isotropic. The spontaneous motions seen in simulations achieve large degrees of bending seen in the X-ray structure and would be entirely sufficient to allow the dynamics of the upper part of helix 44 evidenced by cryo-electron microscopic studies. Analysis of local base-pair step deformability reveals a patch of flexible steps in the upper part of helix 44 and in the area proximal to the bulge bases, suggesting that the upper part of helix 44 has enhanced flexibility. The simulations identify two conformational substates of the second bulge area (bottom part of the helix) with distinct base pairing. In agreement with nuclear magnetic resonance (NMR) and X-ray studies, a flipped out conformational substate of conserved 1492A is seen in the first bulge area. Molecular dynamics (MD) simulations reveal a number of reversible alpha-gamma backbone flips that correspond to transitions between two known A-RNA backbone families. The flipped substates do not cumulate along the trajectory and lead to a modest transient reduction of helical twist with no significant influence on the overall geometry of the duplexes. Despite their considerable flexibility, the simulated structures are very stable with no indication of substantial force field inaccuracies.

Ribose 2'-hydroxyl groups in the 5' strand of the acceptor arm of P-site tRNA are not essential for EF-G catalyzed translocation

The coupled movement of tRNA-mRNA complex through the ribosome is a fundamental step during the protein elongation process. We demonstrate that the ribosome will translocate a P-site-bound tRNA(Met) with a break in the phosphodiester backbone between positions 17 and 18 in the D-loop. Crystallographic data showed that the acceptor arms of P- and E-site tRNA interact extensively with the ribosomal large subunit. Therefore, we used this fragmented P-site-bound tRNA(Met) to investigate the contributions of single 2'-hydroxyl groups in the 5' strand of the acceptor arm for translocation into the ribosomal E-site. EF-G-dependent translocation of the tRNAs was monitored using a toeprinting assay and a fluorescence-based rapid kinetic method. Surprisingly, our results show that none of the 2'-hydroxyl groups in the 5' strand of the acceptor arm of P-site-bound tRNA(Met) between positions 1-17 play a critical role during translocation. This suggests that either these 2'-hydroxyl groups are not important for translocation or they are redundant and the three-dimensional shape of the P-site tRNA is more important for translocation.

Molecular dynamics simulations of sarcin-ricin rRNA motif

Explicit solvent molecular dynamics (MD) simulations were carried out for sarcin-ricin domain (SRD) motifs from 23S (Escherichia coli) and 28S (rat) rRNAs. The SRD motif consists of GAGA tetraloop, G-bulged cross-strand A-stack, flexible region and duplex part. Detailed analysis of the overall dynamics, base pairing, hydration, cation binding and other SRD features is presented. The SRD is surprisingly static in multiple 25 ns long simulations and lacks any non-local motions, with root mean square deviation (r.m.s.d.) values between averaged MD and high-resolution X-ray structures of 1-1.4 A. Modest dynamics is observed in the tetraloop, namely, rotation of adenine in its apex and subtle reversible shift of the tetraloop with respect to the adjacent base pair. The deformed flexible region in low-resolution rat X-ray structure is repaired by simulations. The simulations reveal few backbone flips, which do not affect positions of bases and do not indicate a force field imbalance. Non-Watson-Crick base pairs are rigid and mediated by long-residency water molecules while there are several modest cation-binding sites around SRD. In summary, SRD is an unusually stiff rRNA building block. Its intrinsic structural and dynamical signatures seen in simulations are strikingly distinct from other rRNA motifs such as Loop E and Kink-turns.

Modeling the highly specific ribotoxin recognition of ribosomes

The three-dimensional structures of the alpha-sarcin ribotoxin and its delta(7-22) deletion mutant, both complexed with a 20-mer oligonucleotide mimicking the sarcin/ricin loop (SRL) of the ribosome, have been docked into the structure of the Halobacterium marismortui ribosome by fitting the nucleotide atomic coordinates into those of the ribosomal SRL. This study has revealed that two regions of the ribotoxin, residues 11-16 and 84-85, contact the ribosomal proteins L14 (residues 99-105) and L6 (residues 88-92), respectively. The first of these two ribotoxin regions appears to be crucial for its specific ribosome recognition.

Simulating movement of tRNA into the ribosome during decoding

Decoding is the key step during protein synthesis that enables information transfer from RNA to protein, a process critical for the survival of all organisms. We have used large-scale (2.64 x 10(6) atoms) all-atom simulations of the entire ribosome to understand a critical step of decoding. Although the decoding problem has been studied for more than four decades, the rate-limiting step of cognate tRNA selection has only recently been identified. This step, known as accommodation, involves the movement inside the ribosome of the aminoacyl-tRNA from the partially bound "A/T" state to the fully bound "A/A" state. Here, we show that a corridor of 20 universally conserved ribosomal RNA bases interacts with the tRNA during the accommodation movement. Surprisingly, the tRNA is impeded by the A-loop (23S helix 92), instead of enjoying a smooth transition to the A/A state. In particular, universally conserved 23S ribosomal RNA bases U2492, C2556, and C2573 act as a 3D gate, causing the acceptor stem to pause before allowing entrance into the peptidyl transferase center. Our simulations demonstrate that the flexibility of the acceptor stem of the tRNA, in addition to flexibility of the anticodon arm, is essential for tRNA selection. This study serves as a template for simulating conformational changes in large (>10(6) atoms) biological and artificial molecular machines.

Homology modeling based solution structure of Hoxc8-DNA complex: role of context bases outside TAAT stretch

The 3D structure of neither Hoxc8 nor Hoxc8-DNA complex is known. The repressor protein Hoxc8 binds to the TAAT stretch of the promoter of the osteopontin gene and modulates its expression. Over expression of the osteopontin gene is related to diseases like osteoporosis, multiple sclerosis, cancer et cetera. In this paper we have proposed a 3D structure of Hoxc8-DNA complex obtained by Homology modeling and molecular dynamics (MD) simulation in explicit water. The crystal structure (9ant.pdb) of Antennapedia homeodomain in complex with its DNA sequence was chosen as the template based on (i) high sequence identity (85% for the protein and 60% for the DNA) and (ii) the presence of the TAAT stretch in interaction with the protein. The resulting model was refined by MD simulation for 2.0ns in explicit water. This refined model was then characterized in terms of the structural and the interactional features to improve our understanding of the mechanism of Hoxc8-DNA recognition. The interaction pattern shows that the residues Ile(195), Gln(198), and Asn(199), and the bases S2-(4)TAATG(8) are most important for recognition suggesting the stretch TAATG as the "true recognition element" in the present case. A strong and long-lived water bridge connecting Gln(198) and the base of S1-C(7) complementary to S2-G(8) was observed. Our predicted model of Hoxc8-DNA complex provides us with features that are consistent with the available experimental data on Hoxc8 and the general features of other homeodomain-DNA complexes. The predictions based on the model are also amenable to experimental verification.