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.

• MeSH
• DNA chemie   MeSH
• katalýza MeSH
• konformace nukleové kyseliny * MeSH
• počítačová simulace MeSH
• RNA chemie   MeSH
• simulace molekulární dynamiky * MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• přehledy MeSH
• Research Support, N.I.H., Extramural MeSH
• Názvy látek
• DNA MeSH
• RNA MeSH

We present extensive explicit solvent molecular dynamics analysis of three RNA three-way junctions (3WJs) from the large ribosomal subunit: the 3WJ formed by Helices 90-92 (H90-H92) of 23S rRNA; the 3WJ formed by H42-H44 organizing the GTPase associated center (GAC) of 23S rRNA; and the 3WJ of 5S rRNA. H92 near the peptidyl transferase center binds the 3'-CCA end of amino-acylated tRNA. The GAC binds protein factors and stimulates GTP hydrolysis driving protein synthesis. The 5S rRNA binds the central protuberance and A-site finger (ASF) involved in bridges with the 30S subunit. The simulations reveal that all three 3WJs possess significant anisotropic hinge-like flexibility between their stacked stems and dynamics within the compact regions of their adjacent stems. The A-site 3WJ dynamics may facilitate accommodation of tRNA, while the 5S 3WJ flexibility appears to be essential for coordinated movements of ASF and 5S rRNA. The GAC 3WJ may support large-scale dynamics of the L7/L12-stalk region. The simulations reveal that H42-H44 rRNA segments are not fully relaxed and in the X-ray structures they are bent towards the large subunit. The bending may be related to L10 binding and is distributed between the 3WJ and the H42-H97 contact.

• MeSH
• archeální RNA chemie   MeSH
• bakteriální RNA chemie   MeSH
• Escherichia coli genetika   MeSH
• fosfáty chemie   MeSH
• Haloarcula marismortui genetika   MeSH
• konformace nukleové kyseliny MeSH
• RNA ribozomální 23S chemie   MeSH
• RNA ribozomální 5S chemie   MeSH
• simulace molekulární dynamiky MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, N.I.H., Extramural MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Názvy látek
• archeální RNA MeSH
• bakteriální RNA MeSH
• fosfáty MeSH
• RNA ribozomální 23S MeSH
• RNA ribozomální 5S MeSH

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.

• MeSH
• chlorid draselný chemie   MeSH
• Deinococcus genetika   MeSH
• elektronová kryomikroskopie MeSH
• Escherichia coli genetika   MeSH
• Haloarcula marismortui genetika   MeSH
• konformace nukleové kyseliny MeSH
• RNA ribozomální 23S chemie   MeSH
• simulace molekulární dynamiky MeSH
• sodík chemie   MeSH
• Thermus thermophilus genetika   MeSH
• velké podjednotky ribozomu archebakteriální chemie   MeSH
• velké podjednotky ribozomu bakteriální chemie   MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, N.I.H., Extramural MeSH
• Názvy látek
• chlorid draselný MeSH
• RNA ribozomální 23S MeSH
• sodík MeSH

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.

• MeSH
• archeální RNA chemie   MeSH
• Haloarcula marismortui genetika   MeSH
• ionty chemie   MeSH
• konformace nukleové kyseliny MeSH
• konzervovaná sekvence MeSH
• molekulární modely * MeSH
• molekulární sekvence - údaje MeSH
• párování bází MeSH
• počítačová simulace MeSH
• pohyb těles MeSH
• RNA ribozomální 23S chemie   MeSH
• sekvence nukleotidů MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• archeální RNA MeSH
• ionty MeSH
• RNA ribozomální 23S MeSH

The hepatitis delta virus (HDV) ribozyme is an RNA enzyme from the human pathogenic HDV. Cations play a crucial role in self-cleavage of the HDV ribozyme, by promoting both folding and chemistry. Experimental studies have revealed limited but intriguing details on the location and structural and catalytic functions of metal ions. Here, we analyze a total of approximately 200 ns of explicit-solvent molecular dynamics simulations to provide a complementary atomistic view of the binding of monovalent and divalent cations as well as water molecules to reaction precursor and product forms of the HDV ribozyme. Our simulations find that an Mg2+ cation binds stably, by both inner- and outer-sphere contacts, to the electronegative catalytic pocket of the reaction precursor, in a position to potentially support chemistry. In contrast, protonation of the catalytically involved C75 in the precursor or artificial placement of this Mg2+ into the product structure result in its swift expulsion from the active site. These findings are consistent with a concerted reaction mechanism in which C75 and hydrated Mg2+ act as general base and acid, respectively. Monovalent cations bind to the active site and elsewhere assisted by structurally bridging long-residency water molecules, but are generally delocalized.

• MeSH
• hořčík chemie   MeSH
• kationty dvojmocné chemie   MeSH
• kationty jednomocné chemie   MeSH
• konformace nukleové kyseliny MeSH
• molekulární modely MeSH
• molekulární sekvence - údaje MeSH
• RNA katalytická chemie   MeSH
• sekvence nukleotidů MeSH
• sodík chemie   MeSH
• vazebná místa MeSH
• virus hepatitidy delta enzymologie   MeSH
• voda chemie   MeSH
• vodíková vazba MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, N.I.H., Extramural MeSH
• Názvy látek
• hořčík MeSH
• kationty dvojmocné MeSH
• kationty jednomocné MeSH
• RNA katalytická MeSH
• sodík MeSH
• voda MeSH

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.

• MeSH
• endoribonukleasy metabolismus   MeSH
• Escherichia coli genetika   MeSH
• fungální proteiny metabolismus   MeSH
• kationty chemie   MeSH
• konformace nukleové kyseliny MeSH
• krysa rodu Rattus MeSH
• krystalografie rentgenová MeSH
• molekulární modely * MeSH
• párování bází MeSH
• počítačová simulace MeSH
• ricin metabolismus   MeSH
• RNA ribozomální 23S chemie   metabolismus   MeSH
• RNA ribozomální 28S chemie   metabolismus   MeSH
• sacharidy chemie   MeSH
• vazebná místa MeSH
• voda chemie   MeSH
• vodíková vazba MeSH
• zvířata MeSH
• Check Tag
• krysa rodu Rattus MeSH
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Názvy látek
• alpha-sarcin MeSH Prohlížeč
• endoribonukleasy MeSH
• fungální proteiny MeSH
• kationty MeSH
• ricin MeSH
• RNA ribozomální 23S MeSH
• RNA ribozomální 28S MeSH
• sacharidy MeSH
• voda MeSH

Kink-turn (K-turn) motifs are asymmetric internal loops found at conserved positions in diverse RNAs, with sharp bends in phosphodiester backbones producing V-shaped structures. Explicit-solvent molecular dynamics simulations were carried out for three K-turns from 23S rRNA, i.e., Kt-38 located at the base of the A-site finger, Kt-42 located at the base of the L7/L12 stalk, and Kt-58 located in domain III, and for the K-turn of human U4 snRNA. The simulations reveal hinge-like K-turn motions on the nanosecond timescale. The first conserved A-minor interaction between the K-turn stems is entirely stable in all simulations. The angle between the helical arms of Kt-38 and Kt-42 is regulated by local variations of the second A-minor (type I) interaction between the stems. Its variability ranges from closed geometries to open ones stabilized by insertion of long-residency waters between adenine and cytosine. The simulated A-minor geometries fully agree with x-ray data. Kt-58 and Kt-U4 exhibit similar elbow-like motions caused by conformational change of the adenosine from the nominally unpaired region. Despite the observed substantial dynamics of K-turns, key tertiary interactions are stable and no sign of unfolding is seen. We suggest that some K-turns are flexible elements mediating large-scale ribosomal motions during the protein synthesis cycle.

• MeSH
• adenin chemie   MeSH
• aminokyselinové motivy MeSH
• biofyzika metody   MeSH
• časové faktory MeSH
• cytosin chemie   MeSH
• elongační faktor G chemie   MeSH
• katalýza MeSH
• konformace nukleové kyseliny MeSH
• konformace proteinů MeSH
• krystalografie rentgenová MeSH
• makromolekulární látky MeSH
• molekulární konformace MeSH
• molekulární modely MeSH
• molekulární sekvence - údaje MeSH
• oscilometrie MeSH
• párování bází * MeSH
• počítačová simulace MeSH
• rentgenové záření MeSH
• ribozomy chemie   MeSH
• RNA malá jaderná chemie   MeSH
• RNA ribozomální 23S chemie   MeSH
• RNA transferová chemie   MeSH
• RNA chemie   MeSH
• sekundární struktura proteinů MeSH
• sekvence nukleotidů MeSH
• software MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, N.I.H., Extramural MeSH
• Research Support, U.S. Gov't, P.H.S. MeSH
• Názvy látek
• adenin MeSH
• cytosin MeSH
• elongační faktor G MeSH
• makromolekulární látky MeSH
• RNA malá jaderná MeSH
• RNA ribozomální 23S MeSH
• RNA transferová MeSH
• RNA MeSH
• U4 small nuclear RNA MeSH Prohlížeč