Performance of Molecular Mechanics Force Fields for RNA Simulations: Stability of UUCG and GNRA Hairpins.

The RNA hairpin loops represent important RNA topologies with indispensable biological functions in RNA folding and tertiary interactions. 5′-UNCG-3′ and 5′-GNRA-3′ RNA tetraloops are the most important classes of RNA hairpin loops. Both tetraloops are highly structured with characteristic signature three-dimensional features and are recurrently seen in functional RNAs and ribonucleoprotein particles. Explicit solvent molecular dynamics (MD) simulation is a computational technique which can efficiently complement the experimental data and provide unique structural dynamics information on the atomic scale. Nevertheless, the outcome of simulations is often compromised by imperfections in the parametrization of simplified pairwise additive empirical potentials referred to also as force fields. We have pointed out in several recent studies that a force field description of single-stranded hairpin segments of nucleic acids may be particularly challenging for the force fields. In this paper, we report a critica...

[1]  P. Moore,et al.  The sarcin/ricin loop, a modular RNA. , 1995, Journal of molecular biology.

[2]  Joachim Frank,et al.  Dynamics of the base of ribosomal A-site finger revealed by molecular dynamics simulations and Cryo-EM , 2009, Nucleic acids research.

[3]  H. Heus,et al.  A network of heterogeneous hydrogen bonds in GNRA tetraloops. , 1996, Journal of molecular biology.

[4]  E. Westhof,et al.  Analysis of RNA motifs. , 2003, Current opinion in structural biology.

[5]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[6]  J. Šponer,et al.  Single Stranded Loops of Quadruplex DNA As Key Benchmark for Testing Nucleic Acids Force Fields. , 2009, Journal of chemical theory and computation.

[7]  Hirotaka Ode,et al.  Force field parameters for rotation around χ torsion axis in nucleic acids , 2008, J. Comput. Chem..

[8]  Seema Chauhan,et al.  Tertiary interactions determine the accuracy of RNA folding. , 2008, Journal of the American Chemical Society.

[9]  M. Bansal,et al.  A self-consistent formulation for analysis and generation of non-uniform DNA structures. , 1989, Journal of biomolecular structure & dynamics.

[10]  J. Feigon,et al.  Site-directed spin labeling studies reveal solution conformational changes in a GAAA tetraloop receptor upon Mg(2+)-dependent docking of a GAAA tetraloop. , 2005, Journal of molecular biology.

[11]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

[12]  Harry A. Stern,et al.  Reparameterization of RNA χ Torsion Parameters for the AMBER Force Field and Comparison to NMR Spectra for Cytidine and Uridine , 2010, Journal of chemical theory and computation.

[13]  T. Cheatham,et al.  Determination of Alkali and Halide Monovalent Ion Parameters for Use in Explicitly Solvated Biomolecular Simulations , 2008, The journal of physical chemistry. B.

[14]  Jaroslav Koča,et al.  Conformations of Flanking Bases in HIV-1 RNA DIS Kissing Complexes Studied by Molecular Dynamics , 2007, Biophysical journal.

[15]  D. Crothers,et al.  Bent helix formation between RNA hairpins with complementary loops. , 1995, Science.

[16]  Jiří Šponer,et al.  Molecular dynamics simulations of sarcin–ricin rRNA motif , 2006, Nucleic acids research.

[17]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[18]  Pavel Hobza,et al.  Stabilization and structure calculations for noncovalent interactions in extended molecular systems based on wave function and density functional theories. , 2010, Chemical reviews.

[19]  T. Xia,et al.  Direct revelation of multiple conformations in RNA by femtosecond dynamics. , 2007, Journal of the American Chemical Society.

[20]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[21]  T. Steitz,et al.  Crystal structure of the ribosomal RNA domain essential for binding elongation factors. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[22]  M. Klein,et al.  Constant pressure molecular dynamics algorithms , 1994 .

[23]  T. Xia.  Taking femtosecond snapshots of RNA conformational dynamics and complexity. , 2008, Current opinion in chemical biology.

[24]  Eric J. Sorin,et al.  Equilibrium conformational dynamics in an RNA tetraloop from massively parallel molecular dynamics , 2010, Nucleic acids research.

[25]  K. Hall,et al.  Thermodynamics of 2'-ribose substitutions in UUCG tetraloops. , 2001, RNA.

[26]  J. Šponer,et al.  Different intrastrand and interstrand contributions to stacking account for roll variations at the alternating purine-pyrimidine sequences in A-DNA and A-RNA. , 1991, Journal of molecular biology.

[27]  Philip C. Bevilacqua,et al.  Contribution of the closing base pair to exceptional stability in RNA tetraloops: roles for molecular mimicry and electrostatic factors. , 2009, Journal of the American Chemical Society.

[28]  K. Hall,et al.  Unrestrained stochastic dynamics simulations of the UUCG tetraloop using an implicit solvation model. , 1999, Biophysical journal.

[29]  C R Woese,et al.  Architecture of ribosomal RNA: constraints on the sequence of "tetra-loops". , 1990, Proceedings of the National Academy of Sciences of the United States of America.

[30]  A. Garcia,et al.  Simulation of the pressure and temperature folding/unfolding equilibrium of a small RNA hairpin. , 2008, Journal of the American Chemical Society.

[31]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[32]  T. Steitz,et al.  The kink‐turn: a new RNA secondary structure motif , 2001, The EMBO journal.

[33]  A. Klug,et al.  The crystal structure of an AII-RNAhammerhead ribozyme: A proposed mechanism for RNA catalytic cleavage , 1995, Cell.

[34]  Takeshi Sakata,et al.  Studies on the structure and stabilizing factor of the CUUCGG hairpin RNA using chemically synthesized oligonucleotides , 1990, Nucleic Acids Res..

[35]  E Westhof,et al.  Involvement of a GNRA tetraloop in long-range RNA tertiary interactions. , 1994, Journal of molecular biology.

[36]  Daniel Herschlag,et al.  RNA simulations: probing hairpin unfolding and the dynamics of a GNRA tetraloop. , 2002, Journal of molecular biology.

[37]  F. Javier Luque,et al.  Towards a molecular dynamics consensus view of B-DNA flexibility , 2008, Nucleic acids research.

[38]  G. Varani.  Exceptionally stable nucleic acid hairpins. , 1995, Annual review of biophysics and biomolecular structure.

[39]  P A Kollman,et al.  Theoretical studies of an exceptionally stable RNA tetraloop: observation of convergence from an incorrect NMR structure to the correct one using unrestrained molecular dynamics. , 1997, Journal of molecular biology.

[40]  Y. Chan,et al.  The common and the distinctive features of the bulged-G motif based on a 1.04 A resolution RNA structure. , 2003, Nucleic acids research.

[41]  Leonidas J. Guibas,et al.  Structural Insight into RNA Hairpin Folding Intermediates , 2008, Journal of the American Chemical Society.

[42]  Alexander D. MacKerell,et al.  All‐atom empirical force field for nucleic acids: II. Application to molecular dynamics simulations of DNA and RNA in solution , 2000 .

[43]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[44]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[45]  C. Ehresmann,et al.  Detailed analysis of RNA-protein interactions within the ribosomal protein S8-rRNA complex from the archaeon Methanococcus jannaschii. , 2001, Journal of molecular biology.

[46]  G. Stock,et al.  Molecular dynamics simulation of the structure, dynamics, and thermostability of the RNA hairpins uCACGg and cUUCGg. , 2008, The journal of physical chemistry. B.

[47]  P. Bevilacqua,et al.  Structures, kinetics, thermodynamics, and biological functions of RNA hairpins. , 2008, Annual review of physical chemistry.

[48]  C. G. Hoogstraten,et al.  Extensive backbone dynamics in the GCAA RNA tetraloop analyzed using 13C NMR spin relaxation and specific isotope labeling. , 2008, Journal of the American Chemical Society.

[49]  K. Réblová,et al.  Conformational transitions of flanking purines in HIV‐1 RNA dimerization initiation site kissing complexes studied by CHARMM explicit solvent molecular dynamics , 2008, Biopolymers.

[50]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[51]  G. Varani,et al.  Structure of the P1 helix from group I self-splicing introns. , 1995, Journal of molecular biology.

[52]  A. Serganov,et al.  The crystal structure of UUCG tetraloop. , 2000, Journal of molecular biology.

[53]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[54]  J. Šponer,et al.  Theoretical studies of RNA catalysis: hybrid QM/MM methods and their comparison with MD and QM. , 2009, Methods.

[55]  Piotr Cieplak,et al.  Free energy profile of RNA hairpins: a molecular dynamics simulation study. , 2010, Biophysical journal.

[56]  K. Réblová,et al.  Extensive molecular dynamics simulations showing that canonical G8 and protonated A38H+ forms are most consistent with crystal structures of hairpin ribozyme. , 2010, The journal of physical chemistry. B.

[57]  H. Heus,et al.  Structural features that give rise to the unusual stability of RNA hairpins containing GNRA loops. , 1991, Science.

[58]  F. Eckstein,et al.  Dynamics of the RNA hairpin GNRA tetraloop. , 2000, Biochemistry.

[59]  Wei Wang,et al.  Folding of a small RNA hairpin based on simulation with replica exchange molecular dynamics. , 2010, The journal of physical chemistry. B.

[60]  David H Mathews,et al.  Prediction of RNA secondary structure by free energy minimization. , 2006, Current opinion in structural biology.

[61]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[62]  A. Tannenbaum,et al.  Single nucleotide RNA choreography , 2006, Nucleic acids research.

[63]  Gerhard Stock,et al.  Free-energy landscape of RNA hairpins constructed via dihedral angle principal component analysis. , 2009, The journal of physical chemistry. B.

[64]  V. Ramakrishnan,et al.  Functional insights from the structure of the 30S ribosomal subunit and its interactions with antibiotics , 2000, Nature.