Functional conformations of the L11-ribosomal RNA complex revealed by correlative analysis of cryo-EM and molecular dynamics simulations

Ribosomal Protein L11 Selectively Stabilizes a Tertiary Structure of the GTPase Center rRNA Domain

The GTPase Center (GAC) RNA domain in bacterial 23S rRNA is directly bound by ribosomal protein L11, and this complex is essential to ribosome function. Previous cocrystal structures of the 58-nucleotide GAC RNA bound to L11 revealed the intricate tertiary fold of the RNA domain, with one monovalent and several divalent ions located in specific sites within the structure. Here, we report a new crystal structure of the free GAC that is essentially identical to the L11-bound structure, which retains many common sites of divalent ion occupation. This new structure demonstrates that RNA alone folds into its tertiary structure with bound divalent ions. In solution, we find that this tertiary structure is not static, but rather is best described as an ensemble of states. While L11 protein cannot bind to the GAC until the RNA has adopted its tertiary structure, new experimental data show that L11 binds to Mg²⁺-dependent folded states, which we suggest lie along the folding pathway of the RNA. We propose that L11 stabilizes a specific GAC RNA tertiary state, corresponding to the crystal structure, and that this structure reflects the functionally critical conformation of the rRNA domain in the fully assembled ribosome.

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.

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.

YcaO-Dependent Posttranslational Amide Activation: Biosynthesis, Structure, and Function

With advances in sequencing technology, uncharacterized proteins and domains of unknown function (DUFs) are rapidly accumulating in sequence databases and offer an opportunity to discover new protein chemistry and reaction mechanisms. The focus of this review, the formerly enigmatic YcaO superfamily (DUF181), has been found to catalyze a unique phosphorylation of a ribosomal peptide backbone amide upon attack by different nucleophiles. Established nucleophiles are the side chains of Cys, Ser, and Thr which gives rise to azoline/azole biosynthesis in ribosomally synthesized and posttranslationally modified peptide (RiPP) natural products. However, much remains unknown about the potential for YcaO proteins to collaborate with other nucleophiles. Recent work suggests potential in forming thioamides, macroamidines, and possibly additional post-translational modifications. This review covers all knowledge through mid-2016 regarding the biosynthetic gene clusters (BGCs), natural products, functions, mechanisms, and applications of YcaO proteins and outlines likely future research directions for this protein superfamily.

Chemical footprinting reveals conformational changes of 18S and 28S rRNAs at different steps of translation termination on the human ribosome

Translation termination in eukaryotes is mediated by release factors: eRF1, which is responsible for stop codon recognition and peptidyl-tRNA hydrolysis, and GTPase eRF3, which stimulates peptide release. Here, we have utilized ribose-specific probes to investigate accessibility of rRNA backbone in complexes formed by association of mRNA- and tRNA-bound human ribosomes with eRF1•eRF3•GMPPNP, eRF1•eRF3•GTP, or eRF1 alone as compared with complexes where the A site is vacant or occupied by tRNA. Our data show which rRNA ribose moieties are protected from attack by the probes in the complexes with release factors and reveal the rRNA regions increasing their accessibility to the probes after the factors bind. These regions in 28S rRNA are helices 43 and 44 in the GTPase associated center, the apical loop of helix 71, and helices 89, 92, and 94 as well as 18S rRNA helices 18 and 34. Additionally, the obtained data suggest that eRF3 neither interacts with the rRNA ribose-phosphate backbone nor dissociates from the complex after GTP hydrolysis. Taken together, our findings provide new information on architecture of the eRF1 binding site on mammalian ribosome at various translation termination steps and on conformational rearrangements induced by binding of the release factors.

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.

Fluctuations between multiple EF-G-induced chimeric tRNA states during translocation on the ribosome

The coupled translocation of transfer RNA and messenger RNA through the ribosome entails large-scale structural rearrangements, including step-wise movements of the tRNAs. Recent structural work has visualized intermediates of translocation induced by elongation factor G (EF-G) with tRNAs trapped in chimeric states with respect to 30S and 50S ribosomal subunits. The functional role of the chimeric states is not known. Here we follow the formation of translocation intermediates by single-molecule fluorescence resonance energy transfer. Using EF-G mutants, a non-hydrolysable GTP analogue, and fusidic acid, we interfere with either translocation or EF-G release from the ribosome and identify several rapidly interconverting chimeric tRNA states on the reaction pathway. EF-G engagement prevents backward transitions early in translocation and increases the fraction of ribosomes that rapidly fluctuate between hybrid, chimeric and posttranslocation states. Thus, the engagement of EF-G alters the energetics of translocation towards a flat energy landscape, thereby promoting forward tRNA movement.

EF-G enhances the rate of tRNA–mRNA translocation on the ribosome. Here the authors use single-molecule FRET to follow tRNA translocation in real time, identifying new chimeric intermediates and suggesting how EF-G binding and GTP hydrolysis change the energetic landscape of translocation to accelerate forward tRNA movement.

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.

Influence of thiostrepton binding on the ribosomal GTPase associated region characterized by molecular dynamics simulation

The thiostrepton antibiotic inhibits bacterial protein synthesis by binding to a cleft formed by the ribosomal protein L11 and 23S's rRNA helices 43-44 on the 70S ribosome. It was proposed from crystal structures that the ligand restricts L11's N-terminal movement and thus prevents proper translation factor binding. An exact understanding of thiostrepton's impact on the binding site's dynamics at atomistic resolution is still missing. Here we report an all-atom molecular dynamics simulations of the binary L11·rRNA and the ternary L11·rRNA·thiostrepton complex (rRNA = helices 43-44). We demonstrate that thiostrepton directly impacts the binding site's atomic and biomacromolecular dynamics.

Common chaperone activity in the G-domain of trGTPase protects L11-L12 interaction on the ribosome

Translational GTPases (trGTPases) regulate all phases of protein synthesis. An early event in the interaction of a trGTPase with the ribosome is the contact of the G-domain with the C-terminal domain (CTD) of ribosomal protein L12 (L12-CTD) and subsequently interacts with the N-terminal domain of L11 (L11-NTD). However, the structural and functional relationships between L12-CTD and L11-NTD remain unclear. Here, we performed mutagenesis, biochemical and structural studies to identify the interactions between L11-NTD and L12-CTD. Mutagenesis of conserved residues in the interaction site revealed their role in the docking of trGTPases. During docking, loop62 of L11-NTD protrudes into a cleft in L12-CTD, leading to an open conformation of this domain and exposure of hydrophobic core. This unfavorable situation for L12-CTD stability is resolved by a chaperone-like activity of the contacting G-domain. Our results suggest that all trGTPases-regardless of their different specific functions-use a common mechanism for stabilizing the L11-NTD•L12-CTD interactions.

Principal component and clustering analysis on molecular dynamics data of the ribosomal L11·23S subdomain

With improvements in computer speed and algorithm efficiency, MD simulations are sampling larger amounts of molecular and biomolecular conformations. Being able to qualitatively and quantitatively sift these conformations into meaningful groups is a difficult and important task, especially when considering the structure-activity paradigm. Here we present a study that combines two popular techniques, principal component (PC) analysis and clustering, for revealing major conformational changes that occur in molecular dynamics (MD) simulations. Specifically, we explored how clustering different PC subspaces effects the resulting clusters versus clustering the complete trajectory data. As a case example, we used the trajectory data from an explicitly solvated simulation of a bacteria’s L11·23S ribosomal subdomain, which is a target of thiopeptide antibiotics. Clustering was performed, using K-means and average-linkage algorithms, on data involving the first two to the first five PC subspace dimensions. For the average-linkage algorithm we found that data-point membership, cluster shape, and cluster size depended on the selected PC subspace data. In contrast, K-means provided very consistent results regardless of the selected subspace. Since we present results on a single model system, generalization concerning the clustering of different PC subspaces of other molecular systems is currently premature. However, our hope is that this study illustrates a) the complexities in selecting the appropriate clustering algorithm, b) the complexities in interpreting and validating their results, and c) by combining PC analysis with subsequent clustering valuable dynamic and conformational information can be obtained.

A conserved proline switch on the ribosome facilitates the recruitment and binding of trGTPases

When elongation factor G (EF-G) binds to the ribosome, it first makes contact with the C-terminal domain (CTD) of L12 before interacting with the N-terminal domain (NTD) of L11. Here we have identified a universally conserved residue, Pro22 of L11, that functions as a proline switch (PS22), as well as the corresponding center of peptidyl-prolyl cis-trans isomerase (PPIase) activity on EF-G that drives the cis-trans isomerization of PS22. Only the cis configuration of PS22 allows direct contact between the L11 NTD and the L12 CTD. Mutational analyses of both PS22 and the residues of the EF-G PPIase center reveal their function in translational GTPase (trGTPase) activity, protein synthesis and cell survival in Escherichia coli. Finally, we demonstrate that all known universal trGTPases contain an active PPIase center. Our observations suggest that the cis-trans isomerization of the L11 PS22 is a universal event required for an efficient turnover of trGTPases throughout the translation process.

Structural insights into cognate versus near-cognate discrimination during decoding

The structural basis of the tRNA selection process is investigated by cryo-electron microscopy of ribosomes programmed with UGA codons and incubated with ternary complex (TC) containing the near-cognate Trp-tRNA(Trp) in the presence of kirromycin. Going through more than 350 000 images and employing image classification procedures, we find ∼8% in which the TC is bound to the ribosome. The reconstructed 3D map provides a means to characterize the arrangement of the near-cognate aa-tRNA with respect to elongation factor Tu (EF-Tu) and the ribosome, as well as the domain movements of the ribosome. One of the interesting findings is that near-cognate tRNA's acceptor stem region is flexible and CCA end becomes disordered. The data bring direct structural insights into the induced-fit mechanism of decoding by the ribosome, as the analysis of the interactions between small and large ribosomal subunit, aa-tRNA and EF-Tu and comparison with the cognate case (UGG codon) offers clues on how the conformational signals conveyed to the GTPase differ in the two cases.

Molecular dynamics of EF-G during translocation

Elongation factor G (EF-G) plays a crucial role in two stages of mRNA-(tRNA)(2) translocation. First, EF-G•GTP enters the pre-translocational ribosome in its intersubunit-rotated state, with tRNAs in their hybrid (P/E and A/P) positions. Second, a conformational change in EF-G's Domain IV induced by GTP hydrolysis disengages the mRNA-anticodon stem-loops of the tRNAs from the decoding center to advance relative to the small subunit when the ribosome undergoes a backward inter-subunit rotation. These events take place as EF-G undergoes a series of large conformational changes as visualized by cryo-EM and X-ray studies. The number and variety of these structures leave open many questions on how these different configurations form during the dynamic translocation process. To understand the molecular mechanism of translocation, we examined the molecular motions of EF-G in solution by means of molecular dynamics simulations. Our results show: (1) rotations of the super-domain formed by Domains III-V with respect to the super-domain formed by I-II, and rotations of Domain IV with respect to Domain III; (2) flexible conformations of both 503- and 575-loops; (3) large conformational variability in the bound form caused by the interaction between Domain V and the GTPase-associated center; (4) after GTP hydrolysis, the Switch I region seems to be instrumental for effecting the conformational change at the end of Domain IV implicated in the disengagement of the codon-anticodon helix from the decoding center.

Structural dynamics of the box C/D RNA kink-turn and its complex with proteins: the role of the A-minor 0 interaction, long-residency water bridges, and structural ion-binding sites revealed by molecular simulations

Kink-turns (K-turns) are recurrent elbow-like RNA motifs that participate in protein-assisted RNA folding and contribute to RNA dynamics. We carried out a set of molecular dynamics (MD) simulations using parm99 and parmbsc0 force fields to investigate structural dynamics of the box C/D RNA and its complexes with two proteins: native archaeal L7ae protein and human 15.5 kDa protein, originally bound to very similar structure of U4 snRNA. The box C/D RNA forms K-turn with A-minor 0 tertiary interaction between its canonical (C) and noncanonical (NC) stems. The local K-turn architecture is thus different from the previously studied ribosomal K-turns 38 and 42 having A-minor I interaction. The simulations reveal visible structural dynamics of this tertiary interaction involving altogether six substates which substantially contribute to the elbow-like flexibility of the K-turn. The interaction can even temporarily shift to the A-minor I type pattern; however, this is associated with distortion of the G/A base pair in the NC-stem of the K-turn. The simulations show reduction of the K-turn flexibility upon protein binding. The protein interacts with the apex of the K-turn and with the NC-stem. The protein-RNA interface includes long-residency hydration sites. We have also found long-residency hydration sites and major ion-binding sites associated with the K-turn itself. The overall topology of the K-turn remains stable in all simulations. However, in simulations of free K-turn, we observed instability of the key C16(O2')-A7(N1) H-bond, which is a signature interaction of K-turns and which was visibly more stable in simulations of K-turns possessing A-minor I interaction. It may reflect either some imbalance of the force field or it may be a correct indication of early stages of unfolding since this K-turn requires protein binding for its stabilization. Interestingly, the 16(O2')-7(N1) H- bond is usually not fully lost since it is replaced by a water bridge with a tightly bound water, which is adenine-specific similarly as the original interaction. The 16(O2')-7(N1) H-bond is stabilized by protein binding. The stabilizing effect is more visible with the human 15.5 kDa protein, which is attributed to valine to arginine substitution in the binding site. The behavior of the A-minor interaction is force-field-dependent because the parmbsc0 force field attenuates the A-minor fluctuations compared to parm99 simulations. Behavior of other regions of the box C/D RNA is not sensitive to the force field choice. Simulation with net-neutralizing Na(+) and 0.2 M excess salt conditions appear in all aspects equivalent. The simulations show loss of a hairpin tetraloop, which is not part of the K-turn. This was attributed to force field limitations.

An RNA molecular switch: Intrinsic flexibility of 23S rRNA Helices 40 and 68 5'-UAA/5'-GAN internal loops studied by molecular dynamics methods

Functional RNA molecules such as ribosomal RNAs frequently contain highly conserved internal loops with a 5'-UAA/5'-GAN (UAA/GAN) consensus sequence. The UAA/GAN internal loops adopt distinctive structure inconsistent with secondary structure predictions. The structure has a narrow major groove and forms a trans Hoogsteen/Sugar edge (tHS) A/G base pair followed by an unpaired stacked adenine, a trans Watson-Crick/Hoogsteen (tWH) U/A base pair and finally by a bulged nucleotide (N). The structure is further stabilized by a three-adenine stack and base-phosphate interaction. In the ribosome, the UAA/GAN internal loops are involved in extensive tertiary contacts, mainly as donors of A-minor interactions. Further, this sequence can adopt an alternative 2D/3D pattern stabilized by a four-adenine stack involved in a smaller number of tertiary interactions. The solution structure of an isolated UAA/GAA internal loop shows substantially rearranged base pairing with three consecutive non-Watson-Crick base pairs. Its A/U base pair adopts an incomplete cis Watson-Crick/Sugar edge (cWS) A/U conformation instead of the expected Watson-Crick arrangement. We performed 3.1 µs of explicit solvent molecular dynamics (MD) simulations of the X-ray and NMR UAA/GAN structures, supplemented by MM-PBSA free energy calculations, locally enhanced sampling (LES) runs, targeted MD (TMD) and nudged elastic band (NEB) analysis. We compared parm99 and parmbsc0 force fields and net-neutralizing Na(+) vs. excess salt KCl ion environments. Both force fields provide a similar description of the simulated structures, with the parmbsc0 leading to modest narrowing of the major groove. The excess salt simulations also cause a similar effect. While the NMR structure is entirely stable in simulations, the simulated X-ray structure shows considerable widening of the major groove, loss of base-phosphate interaction and other instabilities. The alternative X-ray geometry even undergoes conformational transition towards the solution 2D structure. Free energy calculations confirm that the X-ray arrangement is less stable than the solution structure. LES, TMD and NEB provide a rather consistent pathway for interconversion between the X-ray and NMR structures. In simulations, the incomplete cWS A/U base pair of the NMR structure is water mediated and alternates with the canonical A-U base pair, which is not indicated by the NMR data. Completion of full cWS A/U base pair is prevented by the overall internal loop arrangement. In summary, the simulations confirm that the UAA/GAN internal loop is a molecular switch RNA module that adopts its functional geometry upon specific tertiary contexts.

Dynamics of the base of ribosomal A-site finger revealed by molecular dynamics simulations and Cryo-EM

Helix 38 (H38) of the large ribosomal subunit, with a length of 110 A, reaches the small subunit through intersubunit bridge B1a. Previous cryo-EM studies revealed that the tip of H38 moves by more than 10 A from the non-ratcheted to the ratcheted state of the ribosome while mutational studies implicated a key role of flexible H38 in attenuation of translocation and in dynamical signaling between ribosomal functional centers. We investigate a region including the elbow-shaped kink-turn (Kt-38) in the Haloarcula marismortui archaeal ribosome, and equivalently positioned elbows in three eubacterial species, located at the H38 base. We performed explicit solvent molecular dynamics simulations on the H38 elbows in all four species. They are formed by at first sight unrelated sequences resulting in diverse base interactions but built with the same overall topology, as shown by X-ray crystallography. The elbows display similar fluctuations and intrinsic flexibilities in simulations indicating that the eubacterial H38 elbows are structural and dynamical analogs of archaeal Kt-38. We suggest that this structural element plays a pivotal role in the large motions of H38 and may act as fulcrum for the abovementioned tip motion. The directional flexibility inferred from simulations correlates well with the cryo-EM results.

Interaction of IF2 with the ribosomal GTPase-associated center during 70S initiation complex formation

Addition of an Escherichia coli 50S subunit (50S(Cy5)) containing a Cy5-labeled L11 N-terminal domain (L11-NTD) within the GTPase-associated center (GAC) to an E. coli 30S initiation complex (30SIC(Cy3)) containing Cy3-labeled initiation factor 2 complexed with GTP leads to rapid development of a FRET signal during formation of the 70S initiation complex (70SIC). Initiation factor 2 (IF2) and elongation factor G (EF-G) induce similar changes in ribosome structure. Here we show that such similarities are maintained on a dynamic level as well. Thus, movement of IF2 toward L11-NTD after initial 70S ribosome formation follows GTP hydrolysis and precedes P(i) release, paralleling movement of EF-G following its binding to the ribosome [Seo, H., et al. (2006) Biochemistry 45, 2504-2514], and in both cases, the rate of such movement is slowed if GTP hydrolysis is prevented. The 30SIC(Cy3):50S(Cy5) FRET signal also provides a sensitive probe of the ability of initiation factor 3 to discriminate between a canonical and a noncanonical initiation codon during 70SIC formation. We employ Bacillus stearothermophilus IF2 as a substitute for E. coli IF2 to take advantage of the higher stability of the complexes it forms with E. coli ribosomes. While Bst-IF2 is fully functional in formation of E. coli 70SIC, relative reactivities toward dipeptide formation of 70SICs formed with the two IF2s suggest that the Bst-IF2.GDP complex is more difficult to displace from the GAC than the E. coli IF2.GDP complex.

Molecular dynamics suggest multifunctionality of an adenine imino group in acid-base catalysis of the hairpin ribozyme

Despite numerous structural and biochemical investigations, the catalytic mechanism of hairpin ribozyme self-cleavage remains elusive. To gain insight into the coupling of active site dynamics with activity of this small catalytic RNA, we analyzed a total of approximately 300 ns of molecular dynamics (MD) simulations. Our simulations predict improved global stability for an in vitro selected "gain of function" mutation, which is validated by native gel electrophoretic mobility shift assay. We observe that active site nucleobases and water molecules stabilize a geometry favorable to catalysis through a dynamic hydrogen bonding network. Simulations in which A38 is unprotonated show its N1 move into close proximity of the active site 2'-OH, indicating that A38 may act as a general base during cleavage, a role that has generally been discounted due to the longer distances observed in crystal structures involving inactivating substrate analogs. By contrast, simulations in which N1 of A38 is protonated place N1 in close proximity to the 5'-oxygen leaving group, which supports the proposal that A38 serves as a general acid. In analogy to protein enzymes, we discuss a plausible mechanism in which A38 acts bifunctionally and shuttles a proton directly from the 2'-OH to the 5'-oxygen. Furthermore, our simulations suggest an important role for protonation of N1 of A38 in promoting a favorable geometry similar to that observed in transition-state analog crystal structures, and support previously proposed roles of A38, G8, and long residency water molecules in transition-state stabilization.

Ribosome-induced changes in elongation factor Tu conformation control GTP hydrolysis

In translation, elongation factor Tu (EF-Tu) molecules deliver aminoacyl-tRNAs to the mRNA-programmed ribosome. The GTPase activity of EF-Tu is triggered by ribosome-induced conformational changes of the factor that play a pivotal role in the selection of the cognate aminoacyl-tRNAs. We present a 6.7-A cryo-electron microscopy map of the aminoacyl-tRNA x EF-Tu x GDP x kirromycin-bound Escherichia coli ribosome, together with an atomic model of the complex obtained through molecular dynamics flexible fitting. The model reveals the conformational changes in the conserved GTPase switch regions of EF-Tu that trigger hydrolysis of GTP, along with key interactions, including those between the sarcin-ricin loop and the P loop of EF-Tu, and between the effector loop of EF-Tu and a conserved region of the 16S rRNA. Our data suggest that GTP hydrolysis on EF-Tu is controlled through a hydrophobic gate mechanism.

Recognition of aminoacyl-tRNA: a common molecular mechanism revealed by cryo-EM

The accuracy of ribosomal translation is achieved by an initial selection and a proofreading step, mediated by EF-Tu, which forms a ternary complex with aminoacyl(aa)-tRNA. To study the binding modes of different aa-tRNAs, we compared cryo-EM maps of the kirromycin-stalled ribosome bound with ternary complexes containing Phe-tRNA^Phe, Trp-tRNA^Trp, or Leu-tRNA^LeuI. The three maps suggest a common binding manner of cognate aa-tRNAs in their specific binding with both the ribosome and EF-Tu. All three aa-tRNAs have the same ‘loaded spring' conformation with a kink and twist between the D-stem and anticodon stem. The three complexes are similarly integrated in an interaction network, extending from the anticodon loop through h44 and protein S12 to the EF-Tu-binding CCA end of aa-tRNA, proposed to signal cognate codon–anticodon interaction to the GTPase centre and tune the accuracy of aa-tRNA selection.

RNA in motion

Although RNA duplex regions are highly structured and inflexible, other elements of an RNA molecule are capable of dynamic motions. These flexible regions are the sites of interactions with small molecules, proteins, and other RNAs, yet there are few descriptions of these regions that include the timescale and amplitude of their motions. No one technique is sufficient to accurately describe these motions, but the combination of in vitro methods, particularly NMR relaxation methods, and more robust in silico methods, is beginning to yield the type of data that can be used to understand RNA function. Very few RNAs have been described by both techniques, and here one such RNA and one RNA:protein complex are reviewed.

Multiple-site trimethylation of ribosomal protein L11 by the PrmA methyltransferase

Ribosomal protein L11 is a universally conserved component of the large subunit, and plays a significant role during initiation, elongation, and termination of protein synthesis. In Escherichia coli, the lysine methyltransferase PrmA trimethylates the N-terminal alpha-amino group and the epsilon-amino groups of Lys3 and Lys39. Here, we report four PrmA-L11 complex structures in different orientations with respect to the PrmA active site. Two structures capture the L11 N-terminal alpha-amino group in the active site in a trimethylated post-catalytic state and in a dimethylated state with bound S-adenosyl-L-homocysteine. Two other structures show L11 in a catalytic orientation to modify Lys39 and in a noncatalytic orientation. The comparison of complex structures in different orientations with a minimal substrate recognition complex shows that the binding mode remains conserved in all L11 orientations, and that substrate orientation is brought about by the unusual interdomain flexibility of PrmA.

Perturbation of the tRNA tertiary core differentially affects specific steps of the elongation cycle

The tRNA tertiary core region is important for both tRNA stability and activity in the translation elongation cycle. Here we report the effects of mutating each of two highly conserved base pairs in the tertiary core of Phe-tRNA(Phe), 18-55 and 19-56, on rate and equilibrium constants for specific steps of this cycle, beginning with formation of aminoacyl-tRNA.EF-Tu.GTP ternary complexs and culminating with translocation of A-site-bound peptidyl-tRNA into the P-site. We find that codon-dependent binding of aminoacyl-tRNA to the A/T-site and proofreading of near-cognate tRNA are sensitive to perturbation of either base pair; formation of the ternary complex and accommodation from the A/T to the A-site are sensitive to 18-55 perturbation only, and translocation of peptidyl-tRNA from the A- to P-site is insensitive to perturbation of either. These results underline the importance of the core region in promoting the efficiency and accuracy of translation, and they likely reflect different requirements for structural integrity of the core during specific steps of the elongation cycle.

Molecular dynamics simulations of nucleic acid-protein complexes

Molecular dynamics simulation studies of protein-nucleic acid complexes are more complicated than studies of either component alone-the force field has to be properly balanced, the systems tend to become very large, and a careful treatment of solvent and of electrostatic interactions is necessary. Recent investigations into several protein-DNA and protein-RNA systems have shown the feasibility of the simulation approach, yielding results of biological interest not readily accessible to experimental methods.

Transfer RNA in the hybrid P/E state: correlating molecular dynamics simulations with cryo-EM data

Transfer RNA (tRNA) transiently occupies the hybrid P/E state (P/E-tRNA) when mRNA-tRNA are translocated in the ribosome. In this study, we characterize the structure of P/E-tRNA and its interactions with the ribosome by correlating the results from molecular dynamics simulations on free tRNA with the cryo-EM map of P/E-tRNA. In our approach, we show that the cryo-EM map may be interpreted as a conformational average. Along the molecular dynamics trajectories (44 ns, 18 ns, and 18 ns), some of the snapshots prove to be quite close to the observed density. In a representative structure, the CCA (3') arm is uniquely twisted, and the anticodon stem loop is kinked at the junctions to both the anticodon loop and the D stem. In addition, the map shows that the P/E-tRNA is no longer bound to helix H69 of 23S rRNA and is flexible, and the conformations of helices H68 and h44 of 16S rRNA differ from those in the x-ray structure. Thus, our study presents structural and dynamic information on the P/E-tRNA and characterizes its interactions with the translocating ribosome.

In vivo assembling of bacterial ribosomal protein L11 into yeast ribosomes makes the particles sensitive to the prokaryotic specific antibiotic thiostrepton

Eukaryotic ribosomal stalk protein L12 and its bacterial orthologue L11 play a central role on ribosomal conformational changes during translocation. Deletion of the two genes encoding L12 in Saccharomyces cerevisiae resulted in a very slow-growth phenotype. Gene RPL12B, but not the RPL12A, cloned in centromeric plasmids fully restored control protein level and the growth rate when expressed in a L12-deprived strain. The same strain has been transformed to express Escherichia coli protein EcL11 under the control of yeast RPL12B promoter. The bacterial protein has been found in similar amounts in washed ribosomes from the transformed yeast strain and from control E. coli cells, however, EcL11 was unable to restore the defective acidic protein stalk composition caused by the absence of ScL12 in the yeast ribosome. Protein EcL11 induced a 10% increase in L12-defective cell growth rate, although the in vitro polymerizing capacity of the EcL11-containing ribosomes is restored in a higher proportion, and, moreover, the particles became partially sensitive to the prokaryotic specific antibiotic thiostrepton. Molecular dynamic simulations using modelled complexes support the correct assembly of bacterial L11 into the yeast ribosome and confirm its direct implication of its CTD in the binding of thiostrepton to ribosomes.

Conformations of flanking bases in HIV-1 RNA DIS kissing complexes studied by molecular dynamics

Explicit solvent molecular dynamics simulations (in total almost 800 ns including locally enhanced sampling runs) were applied with different ion conditions and with two force fields (AMBER and CHARMM) to characterize typical geometries adopted by the flanking bases in the RNA kissing-loop complexes. We focus on flanking base positions in multiple x-ray and NMR structures of HIV-1 DIS kissing complexes and kissing complex from the large ribosomal subunit of Haloarcula marismortui. An initial x-ray open conformation of bulged-out bases in HIV-1 DIS complexes, affected by crystal packing, tends to convert to a closed conformation formed by consecutive stretch of four stacked purine bases. This is in agreement with those recent crystals where the packing is essentially avoided. We also observed variants of the closed conformation with three stacked bases, while nonnegligible populations of stacked geometries with bulged-in bases were detected, too. The simulation results reconcile differences in positions of the flanking bases observed in x-ray and NMR studies. Our results suggest that bulged-out geometries are somewhat more preferred, which is in accord with recent experiments showing that they may mediate tertiary contacts in biomolecular assemblies or allow binding of aminoglycoside antibiotics.

Methylation of proteins involved in translation

Methylation is one of the most common protein modifications. Many different prokaryotic and eukaryotic proteins are methylated, including proteins involved in translation, including ribosomal proteins (RPs) and translation factors (TFs). Positions of the methylated residues in six Escherichia coli RPs and two Saccharomyces cerevisiae RPs have been determined. At least two RPs, L3 and L12, are methylated in both organisms. Both prokaryotic and eukaryotic elongation TFs (EF1A) are methylated at lysine residues, while both release factors are methylated at glutamine residues. The enzymes catalysing methylation reactions, protein methyltransferases (MTases), generally use S-adenosylmethionine as the methyl donor to add one to three methyl groups that, in case of arginine, can be asymetrically positioned. The biological significance of RP and TF methylation is poorly understood, and deletions of the MTase genes usually do not cause major phenotypes. Apparently methylation modulates intra- or intermolecular interactions of the target proteins or affects their affinity for RNA, and, thus, influences various cell processes, including transcriptional regulation, RNA processing, ribosome assembly, translation accuracy, protein nuclear trafficking and metabolism, and cellular signalling. Differential methylation of specific RPs and TFs in a number of organisms at different physiological states indicates that this modification may play a regulatory role.

Elastic properties of ribosomal RNA building blocks: molecular dynamics of the GTPase-associated center rRNA

Explicit solvent molecular dynamics (MD) was used to describe the intrinsic flexibility of the helix 42–44 portion of the 23S rRNA (abbreviated as Kt-42+rGAC; kink-turn 42 and GTPase-associated center rRNA). The bottom part of this molecule consists of alternating rigid and flexible segments. The first flexible segment (Hinge1) is the highly anharmonic kink of Kt-42. The second one (Hinge2) is localized at the junction between helix 42 and helices 43/44. The rigid segments are the two arms of helix 42 flanking the kink. The whole molecule ends up with compact helices 43/44 (Head) which appear to be modestly compressed towards the subunit in the Haloarcula marismortui X-ray structure. Overall, the helix 42–44 rRNA is constructed as a sophisticated intrinsically flexible anisotropic molecular limb. The leading flexibility modes include bending at the hinges and twisting. The Head shows visible internal conformational plasticity, stemming from an intricate set of base pairing patterns including dynamical triads and tetrads. In summary, we demonstrate how rRNA building blocks with contrasting intrinsic flexibilities can form larger architectures with highly specific patterns of preferred low-energy motions and geometries.

Structure of the base of the L7/L12 stalk of the Haloarcula marismortui large ribosomal subunit: analysis of L11 movements

Initiation factors, elongation factors, and release factors all interact with the L7/L12 stalk of the large ribosomal subunit during their respective GTP-dependent cycles on the ribosome. Electron density corresponding to the stalk is not present in previous crystal structures of either 50 S subunits or 70 S ribosomes. We have now discovered conditions that result in a more ordered factor-binding center in the Haloarcula marismortui (H.ma) large ribosomal subunit crystals and consequently allows the visualization of the full-length L11, the N-terminal domain (NTD) of L10 and helices 43 and 44 of 23 S rRNA. The resulting model is currently the most complete reported structure of a L7/L12 stalk in the context of a ribosome. This region contains a series of intermolecular interfaces that are smaller than those typically seen in other ribonucleoprotein interactions within the 50 S subunit. Comparisons of the L11 NTD position between the current structure, which is has an NTD splayed out with respect to previous structures, and other structures of ribosomes in different functional states demonstrates a dynamic range of L11 NTD movements. We propose that the L11 NTD moves through three different relative positions during the translational cycle: apo-ribosome, factor-bound pre-GTP hydrolysis and post-GTP hydrolysis. These positions outline a pathway for L11 NTD movements that are dependent on the specific nucleotide state of the bound ligand. These three states are represented by the orientations of the L11 NTD relative to the ribosome and suggest that L11 may play a more specialized role in the factor binding cycle than previously appreciated.

Molecular dynamics simulations of RNA: an in silico single molecule approach

RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application. Several types of simulations are discussed in detail, including those of structured RNA molecules and their interactions with the surrounding solvent and ions, catalytic RNAs, and RNA-small molecule and RNA-protein complexes. Increased cooperation between theorists and experimentalists will allow expanded judicious use of MD simulations to complement conceptually related single molecule experiments. Such cooperation will open the door to a fundamental understanding of the structure-function relationships in diverse and complex RNA molecules. .

The structure of free L11 and functional dynamics of L11 in free, L11-rRNA(58 nt) binary and L11-rRNA(58 nt)-thiostrepton ternary complexes

The L11 binding site is one of the most important functional sites in the ribosome. The N-terminal domain of L11 has been implicated as a "reversible switch" in facilitating the coordinated movements associated with EF-G-driven GTP hydrolysis. The reversible switch mechanism has been hypothesized to require conformational flexibility involving re-orientation and re-positioning of the two L11 domains, and warrants a close examination of the structure and dynamics of L11. Here we report the solution structure of free L11, and relaxation studies of free L11, L11 complexed to its 58 nt RNA recognition site, and L11 in a ternary complex with the RNA and thiostrepton antibiotic. The binding site of thiostrepton on L11 was also defined by analysis of structural and dynamics data and chemical shift mapping. The conclusions of this work are as follows: first, the binding of L11 to RNA leads to sizable conformation changes in the regions flanking the linker and in the hinge area that links a beta-sheet and a 3(10)-helix-turn-helix element in the N terminus. Concurrently, the change in the relative orientation may lead to re-positioning of the N terminus, as implied by a decrease of radius of gyration from 18.5 A to 16.2 A. Second, the regions, which undergo large conformation changes, exhibit motions on milliseconds-microseconds or nanoseconds-picoseconds time scales. Third, binding of thiostrepton results in more rigid conformations near the linker (Thr71) and near its putative binding site (Leu12). Lastly, conformational changes in the putative thiostrepton binding site are implicated by the re-emergence of cross-correlation peaks in the spectrum of the ternary complex, which were missing in that of the binary complex. Our combined analysis of both the chemical shift perturbation and dynamics data clearly indicates that thiostrepton binds to a pocket involving residues in the 3(10)-helix in L11.

The flexible N-terminal domain of ribosomal protein L11 from Escherichia coli is necessary for the activation of stringent factor

The stringent response is activated by the binding of stringent factor to stalled ribosomes that have an unacylated tRNA in the ribosomal aminoacyl-site. Ribosomes lacking ribosomal protein L11 are deficient in stimulating stringent factor. L11 consists of a dynamic N-terminal domain (amino acid residues 1-72) connected to an RNA-binding C-terminal domain (amino acid residues 76-142) by a flexible linker (amino acid residues 73-75). In vivo data show that mutation of proline 22 in the N-terminal domain is important for initiation of the stringent response. Here, six different L11 point and deletion-mutants have been constructed to determine which regions of L11 are necessary for the activation of stringent factor. The different mutants were reconstituted with programmed 70 S(DeltaL11) ribosomes and tested for their ability to stimulate stringent factor in a sensitive in vitro pppGpp synthesis assay. It was found that a single-site mutation at proline 74 in the linker region between the two domains did not affect the stimulatory activity of the reconstituted ribosomes, whereas the single-site mutation at proline 22 reduced the activity of SF to 33% compared to ribosomes reconstituted with wild-type L11. Removal of the entire linker between the N and C-terminal domains or removal of the entire proline-rich helix beginning at proline 22 in L11 resulted in an L11 protein, which was unable to stimulate stringent factor in the ribosome-dependent assay. Surprisingly, the N-terminal domain of L11 on its own activated stringent factor in a ribosome-dependent manner without restoring the L11 footprint in 23 S rRNA in the 50 S subunit. This suggests that the N-terminal domain can activate stringent factor in trans. It is also shown that this activation is dependent on unacylated tRNA.