A nucleotide-level coarse-grained model of RNA.

We present a new, nucleotide-level model for RNA, oxRNA, based on the coarse-graining methodology recently developed for the oxDNA model of DNA. The model is designed to reproduce structural, mechanical, and thermodynamic properties of RNA, and the coarse-graining level aims to retain the relevant physics for RNA hybridization and the structure of single- and double-stranded RNA. In order to explore its strengths and weaknesses, we test the model in a range of nanotechnological and biological settings. Applications explored include the folding thermodynamics of a pseudoknot, the formation of a kissing loop complex, the structure of a hexagonal RNA nanoring, and the unzipping of a hairpin motif. We argue that the model can be used for efficient simulations of the structure of systems with thousands of base pairs, and for the assembly of systems of up to hundreds of base pairs. The source code implementing the model is released for public use.

[1]  A. Turberfield,et al.  Mechanism for a directional, processive, and reversible DNA motor. , 2009, Small.

[2]  Wade W Grabow,et al.  Self-assembling RNA nanorings based on RNAI/II inverse kissing complexes. , 2011, Nano letters.

[3]  P. Rothemund.  Folding DNA to create nanoscale shapes and patterns , 2006, Nature.

[4]  Song Cao,et al.  Predicting RNA folding thermodynamics with a reduced chain representation model. , 2005, RNA.

[5]  S. Neidle.  Oxford handbook of nucleic acid structure , 1998 .

[6]  J. Doye,et al.  DNA hybridization kinetics: zippering, internal displacement and sequence dependence , 2013, Nucleic acids research.

[7]  Adelene Y. L. Sim,et al.  Modeling nucleic acids. , 2012, Current opinion in structural biology.

[8]  A. E. Walter,et al.  Coaxial stacking of helixes enhances binding of oligoribonucleotides and improves predictions of RNA folding. , 1994, Proceedings of the National Academy of Sciences of the United States of America.

[9]  P. Derreumaux,et al.  HiRE-RNA: a high resolution coarse-grained energy model for RNA. , 2010, The journal of physical chemistry. B.

[10]  F. Ritort,et al.  Force-dependent fragility in RNA hairpins. , 2006, Physical review letters.

[11]  J. Šponer,et al.  Comparison of intrinsic stacking energies of ten unique dinucleotide steps in A-RNA and B-DNA duplexes. Can we determine correct order of stability by quantum-chemical calculations? , 2010, The journal of physical chemistry. B.

[12]  David H. Mathews,et al.  RNAstructure: software for RNA secondary structure prediction and analysis , 2010, BMC Bioinformatics.

[13]  Michael Zuker,et al.  UNAFold: software for nucleic acid folding and hybridization. , 2008, Methods in molecular biology.

[14]  Peng Yin,et al.  Conditional Dicer Substrate Formation via Shape and Sequence Transduction with Small Conditional RNAs , 2013, Journal of the American Chemical Society.

[15]  P. Hagerman,et al.  Flexibility of RNA. , 1997, Annual review of biophysics and biomolecular structure.

[16]  Feng Ding,et al.  RNA-Puzzles: a CASP-like evaluation of RNA three-dimensional structure prediction. , 2012, RNA.

[17]  Yanga Byun,et al.  PseudoViewer: web application and web service for visualizing RNA pseudoknots and secondary structures , 2006, Nucleic Acids Res..

[18]  E Rivas,et al.  A dynamic programming algorithm for RNA structure prediction including pseudoknots. , 1998, Journal of molecular biology.

[19]  Syma Khalid,et al.  Single-stranded DNA within nanopores: conformational dynamics and implications for sequencing; a molecular dynamics simulation study. , 2012, Biophysical journal.

[20]  D. Jost,et al.  Prediction of RNA multiloop and pseudoknot conformations from a lattice-based, coarse-grain tertiary structure model. , 2010, The Journal of chemical physics.

[21]  D. Baker,et al.  Atomic accuracy in predicting and designing non-canonical RNA structure , 2010, Nature Methods.

[22]  G. Weber.  Mesoscopic model parametrization of hydrogen bonds and stacking interactions of RNA from melting temperatures , 2012, Nucleic acids research.

[23]  Peter Minary,et al.  Dynamical Spatial Warping: A Novel Method for the Conformational Sampling of Biophysical Structure , 2008, SIAM J. Sci. Comput..

[24]  B. Shapiro,et al.  Coarse-graining RNA nanostructures for molecular dynamics simulations , 2010, Physical biology.

[25]  D. Turner,et al.  Thermodynamic parameters for an expanded nearest-neighbor model for formation of RNA duplexes with Watson-Crick base pairs. , 1998, Biochemistry.

[26]  Walter Fontana,et al.  Fast folding and comparison of RNA secondary structures , 1994 .

[27]  Flavio Romano,et al.  Coarse-grained simulations of DNA overstretching. , 2012, The Journal of chemical physics.

[28]  F. Ding,et al.  Ab initio RNA folding by discrete molecular dynamics: from structure prediction to folding mechanisms. , 2008, RNA.

[29]  N H Dekker,et al.  Single-molecule measurements of the persistence length of double-stranded RNA. , 2005, Biophysical journal.

[30]  D. Giedroc,et al.  Contribution of the intercalated adenosine at the helical junction to the stability of the gag-pro frameshifting pseudoknot from mouse mammary tumor virus. , 2000, RNA.

[31]  Magdalena A. Jonikas,et al.  Coarse-grained modeling of large RNA molecules with knowledge-based potentials and structural filters. , 2009, RNA.

[32]  J. Sabina,et al.  Expanded sequence dependence of thermodynamic parameters improves prediction of RNA secondary structure. , 1999, Journal of molecular biology.

[33]  Christian Laing,et al.  Computational approaches to 3D modeling of RNA , 2010, Journal of physics. Condensed matter : an Institute of Physics journal.

[34]  J. Doye,et al.  DNA nanotweezers studied with a coarse-grained model of DNA. , 2009, Physical review letters.

[35]  T. Schlick,et al.  Computational approaches to RNA structure prediction, analysis, and design. , 2011, Current opinion in structural biology.

[36]  Helmut Grubmüller,et al.  Determining equilibrium constants for dimerization reactions from molecular dynamics simulations , 2011, J. Comput. Chem..

[37]  김삼묘,et al.  “Bioinformatics” 특집을 내면서 , 2000 .

[38]  J. Doye,et al.  The effect of topology on the structure and free energy landscape of DNA kissing complexes. , 2012, The Journal of chemical physics.

[39]  Conrad Steenberg,et al.  NUPACK: Analysis and design of nucleic acid systems , 2011, J. Comput. Chem..

[40]  L. Jaeger,et al.  In vitro Assembly of Cubic RNA-Based Scaffolds Designed in silico , 2010, Nature nanotechnology.

[41]  I. Tinoco,et al.  The structure of an RNA pseudoknot that causes efficient frameshifting in mouse mammary tumor virus. , 1995, Journal of molecular biology.

[42]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[43]  Christian Matek,et al.  DNA cruciform arms nucleate through a correlated but asynchronous cooperative mechanism. , 2012, The journal of physical chemistry. B.

[44]  Jonathan Bath,et al.  Optimizing DNA nanotechnology through coarse-grained modeling: a two-footed DNA walker. , 2013, ACS nano.

[45]  J. SantaLucia,et al.  The thermodynamics of DNA structural motifs. , 2004, Annual review of biophysics and biomolecular structure.

[46]  David H. Mathews,et al.  RNAstructure: web servers for RNA secondary structure prediction and analysis , 2013, Nucleic Acids Res..

[47]  L M Adleman,et al.  Molecular computation of solutions to combinatorial problems. , 1994, Science.

[48]  宁北芳,et al.  疟原虫var基因转换速率变化导致抗原变异[英]／Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A , 2005 .

[49]  D. Turner,et al.  Incorporating chemical modification constraints into a dynamic programming algorithm for prediction of RNA secondary structure. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[50]  D Thirumalai,et al.  Coarse-grained model for predicting RNA folding thermodynamics. , 2013, The journal of physical chemistry. B.

[51]  S. Wereley,et al.  Soft Matter , 2014 .

[52]  E. Herrero-Galán,et al.  Mechanical identities of RNA and DNA double helices unveiled at the single-molecule level. , 2013, Journal of the American Chemical Society.

[53]  D. Turner,et al.  A set of nearest neighbor parameters for predicting the enthalpy change of RNA secondary structure formation , 2006, Nucleic acids research.

[54]  Yaroslava G. Yingling,et al.  Computational design of an RNA hexagonal nanoring and an RNA nanotube. , 2007, Nano letters.

[55]  Adelene Y. L. Sim,et al.  Modeling and design by hierarchical natural moves , 2012, Proceedings of the National Academy of Sciences.

[56]  H. Hansma,et al.  Building Programmable Jigsaw Puzzles with RNA , 2004, Science.

[57]  Dmitrii V Pyshnyi,et al.  The Influence of Nearest Neighbours on the Efficiency of Coaxial Stacking at Contiguous Stacking Hybridization of Oligodeoxyribonucleotides , 2004, Nucleosides, nucleotides & nucleic acids.

[58]  Aya Kojima,et al.  fRNAdb: a platform for mining/annotating functional RNA candidates from non-coding RNA sequences , 2006, Nucleic Acids Res..

[59]  P. Schuster,et al.  RNA folding at elementary step resolution. , 1999, RNA.

[60]  Changbong Hyeon,et al.  Pathways and kinetic barriers in mechanical unfolding and refolding of RNA and proteins. , 2006, Structure.

[61]  Michal Otyepka,et al.  Computer Folding of RNA Tetraloops? Are We There Yet? , 2013, Journal of chemical theory and computation.

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

[63]  J. Doye,et al.  Structural, mechanical, and thermodynamic properties of a coarse-grained DNA model. , 2010, The Journal of chemical physics.

[64]  BMC Bioinformatics , 2005 .

[65]  P. Derreumaux,et al.  Coarse-grained simulations of RNA and DNA duplexes. , 2013, The journal of physical chemistry. B.

[66]  M. Zuker,et al.  Combining temperature and force to study folding of an RNA hairpin. , 2014, Physical chemistry chemical physics : PCCP.

[67]  D. Turner,et al.  Predicting thermodynamic properties of RNA. , 1995, Methods in enzymology.

[68]  David Rueda,et al.  Thermodynamic and kinetic analysis of an RNA kissing interaction and its resolution into an extended duplex. , 2012, Biophysical journal.

[69]  R. Rosenfeld.  Nature , 2009, Otolaryngology--head and neck surgery : official journal of American Academy of Otolaryngology-Head and Neck Surgery.

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

[71]  G. Fredrickson.  The theory of polymer dynamics , 1996 .

[72]  I. Z. Reguly,et al.  A comparison between parallelization approaches in molecular dynamics simulations on GPUs , 2014, J. Comput. Chem..

[73]  Lorenzo Rovigatti,et al.  Coarse-graining DNA for simulations of DNA nanotechnology. , 2013, Physical chemistry chemical physics : PCCP.

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

[75]  Alain Xayaphoummine,et al.  Kinefold web server for RNA/DNA folding path and structure prediction including pseudoknots and knots , 2005, Nucleic Acids Res..

[76]  D. Thirumalai,et al.  Mechanical unfolding of RNA hairpins. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[77]  T. Ouldridge.  Inferring bulk self-assembly properties from simulations of small systems with multiple constituent species and small systems in the grand canonical ensemble. , 2012, The Journal of chemical physics.

[78]  Y. Benenson.  RNA-based computation in live cells. , 2009, Current opinion in biotechnology.

[79]  J. Šponer,et al.  Reference simulations of noncanonical nucleic acids with different χ variants of the AMBER force field: quadruplex DNA, quadruplex RNA and Z-DNA. , 2012, Journal of chemical theory and computation.

[80]  Henri Orland,et al.  TT2NE: a novel algorithm to predict RNA secondary structures with pseudoknots , 2010, Nucleic acids research.

[81]  A. Louis.  Beware of density dependent pair potentials , 2002, cond-mat/0205110.

[82]  F. Ritort,et al.  Non-specific binding of Na+ and Mg2+ to RNA determined by force spectroscopy methods , 2012, Nucleic acids research.

[83]  Stephen Neidle,et al.  Principles of nucleic acid structure , 2007 .

[84]  D. Thirumalai,et al.  Mechanical unfolding of RNA: from hairpins to structures with internal multiloops. , 2006, Biophysical journal.

[85]  John Russo,et al.  Reversible gels of patchy particles: role of the valence. , 2009, The Journal of chemical physics.

[86]  Pengyu Ren,et al.  RNA 3D structure prediction by using a coarse-grained model and experimental data. , 2013, The journal of physical chemistry. B.

[87]  F. Young.  Biochemistry , 1955, The Indian Medical Gazette.

[88]  Jejoong Yoo,et al.  In situ structure and dynamics of DNA origami determined through molecular dynamics simulations , 2013, Proceedings of the National Academy of Sciences.

[89]  S. Whitelam,et al.  The role of collective motion in examples of coarsening and self-assembly. , 2008, Soft Matter.

[90]  J. Doye,et al.  Sequence-dependent thermodynamics of a coarse-grained DNA model. , 2012, The Journal of chemical physics.

[91]  D. Draper,et al.  Persistence length of RNA. , 1995, Biochemistry.

[92]  T. Nyholm,et al.  Nucleosides Nucleotides Nucleic Acids , 2003 .

[93]  Ivo L. Hofacker,et al.  Vienna RNA secondary structure server , 2003, Nucleic Acids Res..

[94]  Tamar Schlick,et al.  Molecular Modeling and Simulation: An Interdisciplinary Guide , 2010 .

[95]  J. Doye,et al.  Extracting bulk properties of self-assembling systems from small simulations , 2009, Journal of physics. Condensed matter : an Institute of Physics journal.

[96]  Peixuan Guo.  The emerging field of RNA nanotechnology. , 2010, Nature nanotechnology.

[97]  Joseph M. Schaeffer,et al.  On the biophysics and kinetics of toehold-mediated DNA strand displacement , 2013, Nucleic acids research.

[98]  J. H. Maddocks,et al.  Conformational analysis of nucleic acids revisited: Curves+ , 2009, Nucleic acids research.

[99]  Pavel Hobza,et al.  Accurate interaction energies of hydrogen-bonded nucleic acid base pairs. , 2004, Journal of the American Chemical Society.

[100]  T. Odijk.  Stiff chains and filaments under tension , 1995 .

[101]  A. E. Walter,et al.  Sequence dependence of stability for coaxial stacking of RNA helixes with Watson-Crick base paired interfaces. , 1994, Biochemistry.

[102]  N. Seeman.  Nanomaterials based on DNA. , 2010, Annual review of biochemistry.

[103]  F. Major,et al.  The MC-Fold and MC-Sym pipeline infers RNA structure from sequence data , 2008, Nature.