Highly sampled tetranucleotide and tetraloop motifs enable evaluation of common RNA force fields

RNA

Volumn 21, Issue 9, 2015, Pages 1578-1590

  Bergonzo, Christina ^a   Henriksen, Niel M ^a,b   Roe, Daniel R ^a   Cheatham, Thomas E ^a  

^a UNIVERSITY OF UTAH   (United States)

^b UNIVERSITY OF CALIFORNIA SAN DIEGO   (United States)

Author keywords

AMBER; CHARMM; Enhanced sampling; Force fields; Molecular dynamics; Replica exchange; RNA

Indexed keywords

RNA;

ARTICLE; DNA HELIX; HYDROGEN BOND; MOLECULAR DYNAMICS; MULTIDIMENSIONAL REPLICA EXCHANGE MOLECULAR DYNAMICS; NUCLEOTIDE MOTIF; PRIORITY JOURNAL; RNA CONFORMATION; RNA STRUCTURE; TETRANUCLEOTIDE MOTIF; THERMODYNAMICS; UUCG TETRALOOP MOTIF; BIOLOGY; CHEMISTRY; CONFORMATION; NUCLEAR MAGNETIC RESONANCE SPECTROSCOPY; PROCEDURES;

COMPUTATIONAL BIOLOGY; MAGNETIC RESONANCE SPECTROSCOPY; MOLECULAR DYNAMICS SIMULATION; NUCLEIC ACID CONFORMATION; NUCLEOTIDE MOTIFS; RNA; THERMODYNAMICS;

EID: 84939789620     PISSN: 13558382     EISSN: 14699001     Source Type: Journal    
DOI: 10.1261/rna.051102.115     Document Type: Article

References (56)

1. Banáš P, Hollas D, Zgarbová M, Jurecka P, Orozco M, Cheatham TE III, Sponer J, Otyepka M. 2010. Performance of molecular mechanics force fields for RNA simulations: stability of UUCG and GNRA hairpins. J Chem Theory Comput 6: 3836-3849.
2. Berendsen HJC, Postma JPM, van Gunsteren WF, DiNola A, Haak JR. 1984. Molecular dynamics with coupling to an external bath. J Chem Phys 81: 3684-3690.
3. Bergonzo C, Henriksen NM, Roe DR, Swails JM, Roitberg AE, Cheatham TE III. 2014. Multidimensional replica exchange molecular dynamics yields a converged ensemble of an RNA tetranucleo-tide. J Chem Theory Comput 10: 492-499.
4. Bevilacqua PC, Blose JM. 2008. Structures, kinetics, thermodynamics, and biological functions of RNA hairpins. Annu Rev Phys Chem 59: 79-103.
5. Case DA, Darden TA, Cheatham TE III, Simmerling CL, Wang J, Duke RE, Luo R, Walker RC, Zhang W, Merz KM, et al. 2012. AMBER 12. University of California, San Francisco, CA.
6. Case DA, Darden TA, Cheatham TE III, Simmerling CL, Roitberg A, Wang J, Duke RE, Luo R, Roe DR, Walker RC, et al. 2014. AMBER 14. University of California, San Francisco, CA.
7. Cheatham TE III, Case DA. 2013. Twenty-five years of nucleic acid simulations. Biopolymers 99: 969-977.
8. Cheatham TE III, Cieplak P, Kollman PA. 1999. A modified version of the Cornell et al. force field with improved sugar pucker phases and helical repeat. J Biomol Struct Dyn 16: 845-862.
9. Chen AA, García AE. 2013. High-resolution reversible folding of hyper-stable RNA tetraloops using molecular dynamics simulations. Proc Natl Acad Sci 110: 16820-16825.
10. Crowley MF, Williamson MJ, Walker RC. 2009. CHAMBER: comprehensive support for CHARMM force fields within the AMBER software. Int J Quant Chem 109: 3767-3772.
11. ′-hydroxyl sampling on the conformational properties of RNA: update of the CHARMM all-atom additive force field for RNA. J Comput Chem 32: 1929-1943.
12. Ennifar E, Nikulin A, Tishchenko S, Serganov A, Nevskaya N, Garber M, Ehresmann B, Ehresmann C, Nikonov S, Dumas P. 2000. The crystal structure of UUCG tetraloop. J Mol Biol 304: 35-42.
13. Feenstra KA, Hess B, Berendsen HJC. 1999. Improving efficiency of large time-scale molecular dynamics simulations of hydrogen-rich systems. J Comput Chem 20: 786-798.
14. Foloppe N, MacKerell AD Jr. 2000. All-atom empirical force field for nucleic acids: I. Parameter optimization based on small molecule and condensed phase macromolecular target data. J Comput Chem 21: 86-104.
15. Galindo-Murillo R, Roe DR, Cheatham TE III. 2014a. Convergence and reproducibility in molecular dynamics simulations of the DNA duplex d(GCACGAACGAACGAACGC). Biochim Biophys Acta 1850: 1041-1058.
16. Galindo-Murillo R, Roe DR, Cheatham TE III. 2014b. On the absence of intrahelical DNA dynamics on the μs to ms timescale. Nat Commun 5: 5152.
17. Giambaşu GM, York DM, Case DA. 2015. Structural fidelity and NMR relaxation analysis in a prototype RNA hairpin. RNA 21: 963-974.
18. Götz AW, Williamson MJ, Xu D, Poole D, Le Grand S, Walker RC. 2012. Routine microsecond molecular dynamics simulations with AMBER on GPUs. 1. Generalized born. J Chem Theory Comput 8: 1542-1555.
19. Hamelberg D, Mongan J, McCammon JA. 2004. Accelerated molecular dynamics: a promising and efficient simulation method for biomol-ecules. J Chem Phys 120: 11919-11929.
20. Hansmann UHE. 1997. Parallel tempering algorithm for conformation-al studies of biological molecules. Chem Phys Lett 281: 140-150.
21. Henriksen NM, Roe DR, Cheatham TE III. 2013. Reliable oligonucleo-tide conformational ensemble generation in explicit solvent for force field assessment using reservoir replica exchange molecular dynamics simulations. J Phys Chem B 117: 4014-4027.
22. Hopkins CW, Le Grand S, Walker RC, Roitberg AE. 2015. Long time step molecular dynamics through hydrogen mass repartitioning. J Chem Theory Comput 11: 1864-1874.
23. Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML. 1983. Comparison of simple potential functions for simulating liquid water. J Chem Phys 79: 926-935.
24. Joung IS, Cheatham TE III. 2008. Determination of alkali and halide monovalent ion parameters for use in explicitly solvated biomolec-ular simulations. J Phys Chem B 112: 9020-9041.
25. Kührová P, Banáš P, Best RB, Sponer J, Otyepka M. 2013. Computer folding of RNA tetraloops? Are we there yet? J Chem Theory Comput 9: 2115-2125.
26. Laing C, Schlick T. 2011. Computational approaches to RNA structure prediction, analysis, and design. Curr Opin Struct Biol 21: 306-318.
27. ′-methylamide. Biopolymers 32: 523-535.
28. Meister G. 2011. RNA biology: an introduction. Wiley-VCH, Weinheim, Germany.
29. Mlýnský V, Banás P, Hollas D, Réblová K, Walter NG, Sponer J, Otyepka M. 2010. Extensive molecular dynamics simulations showing that canonical G8 and protonated A38H+ forms are most consistent with crystal structures of hairpin ribozyme. J Phys Chem B 114: 6642-6652.
30. Nguyen H, Maier J, Huang H, Perrone V, Simmerling C. 2014. Folding simulations for proteins with diverse topologies are accessible in days with a physics-based force field and implicit solvent. J Am Chem Soc 136: 13959-13962.
31. Nozinovic S, Fürtig B, Jonker HR, Richter C, Schwalbe H. 2010. High-resolution NMR structure of an RNA model system: the 14-mer cUUCGg tetraloop hairpin RNA. Nucleic Acids Res 38: 683-694.
32. Patriksson A, van der Spoel D. 2008. A temperature predictor for parallel tempering simulations. Phys Chem Chem Phys 10: 2073-2077.
33. Pérez A, Luque FJ, Orozco M. 2007a. Dynamics of B-DNA on the microsecond time scale. J Am Chem Soc 101: 14771-14775.
34. Pérez A, Marchán I, Svozil D, Sponer J, Cheatham TE III, Laughton CA, Orozco M. 2007b. Refinement of the AMBER force field for nucleic acids: improving the description of α/γ conformers. Biophys J 92: 3817-3829.
35. Pérez A, Luque FJ, Orozco M. 2012. Frontiers in molecular dynamics simulations of DNA. Acc Chem Res 45: 196-205.
36. Richardson JS, Schneider B, Murray LW, Kapral GJ, Immormino RM, Headd JJ, Richardson DC, Ham D, Hershkovits ELI, Williams LD, et al. 2008. RNA backbone: consensus all-angle conformers and modular string nomenclature (an RNA Ontology Consortium contribution). RNA 14: 465-481.
37. Roe DR, Cheatham TE III. 2013. PTRAJ and CPPTRAJ: software for processing and analysis of molecular dynamics trajectory data. J Chem Theory Comput 9: 3084-3095.
38. Roe DR, Bergonzo C, Cheatham TE III. 2014. Evaluation of enhanced sampling provided by accelerated molecular dynamics with Hamil-tonian replica exchange methods. J Phys Chem B 118: 3543-3552.
39. Ryckaert J-P, Ciccotti G, Berendsen HJC. 1977. Numerical integration of the Cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes. J Comput Phys 23: 327-341.
40. Salomon-Ferrer R, Götz AW, Poole D, Le Grand S, Walker RC. 2013. Routine microsecond molecular dynamics simulations with AMBER on GPUs. 2. Explicit solvent particle mesh Ewald. J Chem Theory Comput 9: 3878-3888.
41. Shao J, Tanner SW, Thompson N, Cheatham TE III. 2007. Clustering molecular dynamics trajectories: 1. Characterizing the performance of different clustering algorithms. J Chem Theory Comput 3: 2312-2334.
42. Shaw DE, Deneroff MM, Dror RO, Kuskin JS, Larson RH, Salmon JK, Young C, Batson B, Bowers KJ, Chao JC, et al. 2008. Anton, a special-purpose machine for molecular dynamics simulations. Commun ACM 51: 91-97.
43. Simmerling C, Strockbine B, Roitberg AE. 2002. All-atom structure prediction and folding simulations of a stable protein. J Am Chem Soc 124: 11258-11259.
44. Sindhikara DJ, Kim S, Voter AF, Roitberg AE. 2009. Bad seeds sprout perilous dynamics: stochastic thermostat induced trajectory synchronization in biomolecules. J Chem Theory Comput 5: 1624-1631.
45. Sorin EJ, Engelhardt MA, Herschlag D, Pande VS. 2002. RNA simulations: probing hairpin unfolding and the dynamics of a GNRA tet-raloop. J Mol Biol 317: 493-506.
46. Steinbrecher T, Latzer J, Case DA. 2012. Revised AMBER parameters for bioorganic phosphates. J Chem Theory Comput 8: 4405-4412.
47. Sugita Y, Okamoto Y. 1999. Replica-exchange molecular dynamics method for protein folding. Chem Phys Lett 314: 141-151.
48. Towns J, Cockerill T, Dahan M, Foster I, Gaither K, Grimshaw A, Hazlewood V, Lathrop S, Lifka D, Peterson GD, et al. 2014. XSEDE: accelerating scientific discovery. Comput Sci Eng 16: 62-74.
49. Tubbs JD, Condon DE, Kennedy SD, Hauser M, Bevilacqua PC, Turner DH. 2013. The nuclear magnetic resonance of CCCC RNA reveals a right-handed helix, and revised parameters for AMBER force field torsions improve structural predictions from molecular dynamics. Biochemistry 52: 996-1010.
50. Uhlenbeck OC 1990. Tetraloops and RNA folding. Nature 346: 613-614.
51. Vaiana AC, Sanbonmatsu KY. 2009. Stochastic gating and drug-ribo-some interactions. J Mol Biol 386: 648-661.
52. Wang J, Cieplak P, Kollman PA. 2000. How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? J Comput Chem 21: 1049-1074.
53. Woese CR, Winkers S, Gutell RR. 1990. Architecture of ribosomal RNA: constraints on the sequence of "tetra-loops". Proc Natl Acad Sci 87: 8467-8471.
54. Yildirim I, Stern HA, Tubbs JD, Kennedy SD, Turner DH. 2011. Benchmarking AMBER force fields for RNA: comparisons to NMR spectra for single-stranded r(GACC) are improved by revised X torsions. J Phys Chem B 115: 9261-9270.
55. Yildirim I, Kennedy SD, Stern HA, Hart JM, Kierzek R, Turner DH. 2012. Revision of AMBER torsional parameters for RNA improves free energy predictions for tetramer duplexes with GC and iGiC base pairs. J Chem Theory Comput 8: 172-181.
56. Zgarbová M, Otyepka M, Sponer J, Mládek A, Banás P, Cheatham TE III, Jurecka P. 2011. Refinement of the Cornell et al. nucleic acids force field based on reference quantum chemical calculations of glycosidic torsion profiles. J Chem Theory Comput 7: 2886-2902.