RNA molecules play a key role in countless biochemical processes. RNA interactions, which are of highly diverse nature, are determined by the fact that RNA is a highly negatively charged polyelectrolyte, which leads to intimate interactions with an ion atmosphere. Although RNA molecules are formally single-stranded, canonical (Watson-Crick) duplexes are key components of folded RNAs. A double-stranded (ds) RNA is also important for the design of RNA-based nanostructures and assemblies. Despite the fact that the description of canonical dsRNA is considered the least problematic part of RNA modeling, the imperfect shape and flexibility of dsRNA can lead to imbalances in the simulations of larger RNAs and RNA-containing assemblies. We present a comprehensive set of molecular dynamics (MD) simulations of four canonical A-RNA duplexes. Our focus was directed toward the characterization of the influence of varying ion concentrations and of the size of the solvation box. We compared several water models and four RNA force fields. The simulations showed that the A-RNA shape was most sensitive to the RNA force field, with some force fields leading to a reduced inclination of the A-RNA duplexes. The ions and water models played a minor role. The effect of the box size was negligible, and even boxes with a small fraction of the bulk solvent outside the RNA hydration sphere were sufficient for the simulation of the dsRNA.

• MeSH
• ionty chemie   MeSH
• konformace nukleové kyseliny MeSH
• RNA * chemie   MeSH
• simulace molekulární dynamiky * MeSH
• voda chemie   MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• ionty MeSH
• RNA * MeSH
• voda MeSH

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

The 22-mer c-kit promoter sequence folds into a parallel-stranded quadruplex with a unique structure, which has been elucidated by crystallographic and NMR methods and shows a high degree of structural conservation. We have carried out a series of extended (up to 10 μs long, ∼50 μs in total) molecular dynamics simulations to explore conformational stability and loop dynamics of this quadruplex. Unfolding no-salt simulations are consistent with a multi-pathway model of quadruplex folding and identify the single-nucleotide propeller loops as the most fragile part of the quadruplex. Thus, formation of propeller loops represents a peculiar atomistic aspect of quadruplex folding. Unbiased simulations reveal μs-scale transitions in the loops, which emphasizes the need for extended simulations in studies of quadruplex loops. We identify ion binding in the loops which may contribute to quadruplex stability. The long lateral-propeller loop is internally very stable but extensively fluctuates as a rigid entity. It creates a size-adaptable cleft between the loop and the stem, which can facilitate ligand binding. The stability gain by forming the internal network of GA base pairs and stacks of this loop may be dictating which of the many possible quadruplex topologies is observed in the ground state by this promoter quadruplex.

• MeSH
• denaturace nukleových kyselin MeSH
• draslík chemie   MeSH
• G-kvadruplexy * MeSH
• kationty MeSH
• párování bází MeSH
• promotorové oblasti (genetika) * MeSH
• protoonkogenní proteiny c-kit genetika   MeSH
• simulace molekulární dynamiky MeSH
• sodík chemie   MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• draslík MeSH
• kationty MeSH
• protoonkogenní proteiny c-kit MeSH
• sodík 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

Hybrid QM/MM methods combine the rigor of quantum mechanical (QM) calculations with the low computational cost of empirical molecular mechanical (MM) treatment allowing to capture dynamic properties to probe critical atomistic details of enzyme reactions. Catalysis by RNA enzymes (ribozymes) has only recently begun to be addressed with QM/MM approaches and is thus still a field under development. This review surveys methodology as well as recent advances in QM/MM applications to RNA mechanisms, including those of the HDV, hairpin, and hammerhead ribozymes, as well as the ribosome. We compare and correlate QM/MM results with those from QM and/or molecular dynamics (MD) simulations, and discuss scope and limitations with a critical eye on current shortcomings in available methodologies and computer resources. We thus hope to foster mutual appreciation and facilitate collaboration between experimentalists and theorists to jointly advance our understanding of RNA catalysis at an atomistic level.

• MeSH
• biofyzika metody   MeSH
• fosfáty chemie   MeSH
• fosforylace MeSH
• hořčík chemie   MeSH
• katalýza MeSH
• konformace nukleové kyseliny MeSH
• kvantová teorie MeSH
• lidé MeSH
• molekulární modely MeSH
• počítačová simulace MeSH
• ribozomy chemie   MeSH
• RNA katalytická chemie   MeSH
• RNA virová chemie   MeSH
• RNA chemie   MeSH
• software MeSH
• Check Tag
• lidé 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
• fosfáty MeSH
• hammerhead ribozyme MeSH Prohlížeč
• hořčík MeSH
• RNA katalytická MeSH
• RNA virová MeSH
• RNA MeSH

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.

• MeSH
• chemické modely * MeSH
• dimerizace MeSH
• HIV-1 chemie   genetika   MeSH
• konformace nukleové kyseliny MeSH
• molekulární modely * MeSH
• párování bází genetika   MeSH
• počátek transkripce * MeSH
• počítačová simulace MeSH
• RNA virová chemie   MeSH
• vazebná místa MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• RNA virová 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č