Molecular Dynamics Simulations Reveal an Interplay between SHAPE Reagent Binding and RNA Flexibility

The function of RNA molecules usually depends on their overall fold and on the presence of specific structural motifs. Chemical probing methods are routinely used in combination with nearest-neighbor models to determine RNA secondary structure. Among the available methods, SHAPE is relevant due to its capability to probe all RNA nucleotides and the possibility to be used in vivo. However, the structural determinants for SHAPE reactivity and its mechanism of reaction are still unclear. Here molecular dynamics simulations and enhanced sampling techniques are used to predict the accessibility of nucleotide analogs and larger RNA structural motifs to SHAPE reagents. We show that local RNA reconformations are crucial in allowing reagents to reach the 2′-OH group of a particular nucleotide and that sugar pucker is a major structural factor influencing SHAPE reactivity.

[1]  Kyle E. Watters,et al.  Using in-cell SHAPE-Seq and simulations to probe structure–function design principles of RNA transcriptional regulators , 2016, RNA.

[2]  D. Lilley.  How RNA acts as a nuclease: some mechanistic comparisons in the nucleolytic ribozymes. , 2017, Biochemical Society transactions.

[3]  Howard Y. Chang,et al.  RNA SHAPE analysis in living cells. , 2013, Nature chemical biology.

[4]  Giovanni Bussi,et al.  The role of nucleobase interactions in RNA structure and dynamics , 2014, Nucleic acids research.

[5]  Giovanni Bussi,et al.  Unraveling Mg2+–RNA binding with atomistic molecular dynamics , 2016, RNA.

[6]  Laura R. Ganser,et al.  Resolving sugar puckers in RNA excited states exposes slow modes of repuckering dynamics , 2017, Nucleic acids research.

[7]  J. Onuchic,et al.  Integrating molecular dynamics simulations with chemical probing experiments using SHAPE-FIT. , 2015, Methods in enzymology.

[8]  K. Weeks,et al.  RNA structure analysis at single nucleotide resolution by selective 2'-hydroxyl acylation and primer extension (SHAPE). , 2005, Journal of the American Chemical Society.

[9]  Kevin Y. Yip,et al.  Improved prediction of RNA secondary structure by integrating the free energy model with restraints derived from experimental probing data , 2015, Nucleic acids research.

[10]  Peter F. Stadler,et al.  SHAPE directed RNA folding , 2015, bioRxiv.

[11]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[12]  P. Walter,et al.  NMR studies of the most conserved RNA domain of the mammalian signal recognition particle (SRP). , 1996, RNA.

[13]  Peter M. Kasson,et al.  GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit , 2013, Bioinform..

[14]  Ryan A. Flynn,et al.  Comparative Analysis Reveals Furoyl in Vivo Selective Hydroxyl Acylation Analyzed by Primer Extension Reagents Form Stable Ribosyl Ester Adducts. , 2017, Biochemistry.

[15]  Gaurav Sharma,et al.  Modeling RNA Secondary Structure with Sequence Comparison and Experimental Mapping Data. , 2017, Biophysical journal.

[16]  J. Šponer,et al.  Noncanonical α/γ Backbone Conformations in RNA and the Accuracy of Their Description by the AMBER Force Field. , 2017, The journal of physical chemistry. B.

[17]  J. Doudna,et al.  Crystal structure of the ribonucleoprotein core of the signal recognition particle. , 2000, Science.

[18]  A. Laio,et al.  A bias-exchange approach to protein folding. , 2007, The journal of physical chemistry. B.

[19]  K. Weeks,et al.  Direct identification of base-paired RNA nucleotides by correlated chemical probing , 2017, RNA.

[20]  Kristen K. Dang,et al.  Architecture and Secondary Structure of an Entire HIV-1 RNA Genome , 2009, Nature.

[21]  K. Morris,et al.  The rise of regulatory RNA , 2014, Nature Reviews Genetics.

[22]  D. Mathews,et al.  Accurate SHAPE-directed RNA secondary structure modeling, including pseudoknots , 2013, Proceedings of the National Academy of Sciences.

[23]  Kyle E. Watters,et al.  Cotranscriptional Folding of a Riboswitch at Nucleotide Resolution , 2016, Nature Structural &Molecular Biology.

[24]  K. Weeks,et al.  The mechanisms of RNA SHAPE chemistry. , 2012, Journal of the American Chemical Society.

[25]  Howard Y. Chang,et al.  Structural imprints in vivo decode RNA regulatory mechanisms , 2015, Nature.

[26]  N. Walter.  The blessing and curse of RNA dynamics: past, present, and future. , 2009, Methods.

[27]  K. Weeks,et al.  Slow conformational dynamics at C2'-endo nucleotides in RNA. , 2008, Journal of the American Chemical Society.

[28]  A. Laederach,et al.  Transcending the prediction paradigm: novel applications of SHAPE to RNA function and evolution , 2016, Wiley interdisciplinary reviews. RNA.

[29]  Niel M. Henriksen,et al.  Highly sampled tetranucleotide and tetraloop motifs enable evaluation of common RNA force fields , 2015, RNA.

[30]  Kyle E. Watters,et al.  Characterizing RNA structures in vitro and in vivo with selective 2'-hydroxyl acylation analyzed by primer extension sequencing (SHAPE-Seq). , 2016, Methods.

[31]  Giovanni Bussi,et al.  Understanding in-line probing experiments by modeling cleavage of nonreactive RNA nucleotides. , 2017, RNA.

[32]  Giovanni Bussi,et al.  Elastic network models for RNA: a comparative assessment with molecular dynamics and SHAPE experiments , 2015, Nucleic acids research.

[33]  B. Shapiro,et al.  Correlating SHAPE signatures with three-dimensional RNA structures. , 2011, RNA.

[34]  K. Weeks,et al.  Ribosome RNA assembly intermediates visualized in living cells. , 2014, Biochemistry.

[35]  S. Acharya,et al.  The pKa's of 2‘-Hydroxyl Group in Nucleosides and Nucleotides , 2001 .

[36]  K. Weeks,et al.  Fingerprinting noncanonical and tertiary RNA structures by differential SHAPE reactivity. , 2012, Journal of the American Chemical Society.

[37]  Christine E. Heitsch,et al.  Evaluating the accuracy of SHAPE-directed RNA secondary structure predictions , 2013, Nucleic acids research.

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

[39]  Julius B. Lucks,et al.  An RNA Mapping DataBase for curating RNA structure mapping experiments , 2012, Bioinform..

[40]  Gabriele Varani,et al.  Strong correlation between SHAPE chemistry and the generalized NMR order parameter (S2) in RNA. , 2008, Journal of the American Chemical Society.

[41]  S. Harvey,et al.  Molecular paleontology: a biochemical model of the ancestral ribosome , 2013, Nucleic acids research.

[42]  K. Weeks,et al.  Exploring RNA structural codes with SHAPE chemistry. , 2011, Accounts of chemical research.

[43]  David H Mathews,et al.  Modeling RNA secondary structure folding ensembles using SHAPE mapping data , 2017, Nucleic acids research.

[44]  Christopher A. Lavender,et al.  In-cell SHAPE reveals that free 30S ribosome subunits are in the inactive state , 2015, Proceedings of the National Academy of Sciences.

[45]  Kaoru Inoue,et al.  SHAPE reveals transcript-wide interactions, complex structural domains, and protein interactions across the Xist lncRNA in living cells , 2016, Proceedings of the National Academy of Sciences.

[46]  Siqi Tian,et al.  High-throughput mutate-map-rescue evaluates SHAPE-directed RNA structure and uncovers excited states. , 2014, RNA.

[47]  Kevin M Weeks,et al.  Selective 2′-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP) for direct, versatile and accurate RNA structure analysis , 2015, Nature Protocols.

[48]  Giovanni Bussi,et al.  Free Energy Landscape of GAGA and UUCG RNA Tetraloops. , 2016, The journal of physical chemistry letters.

[49]  Morgan C. Giddings,et al.  Influence of nucleotide identity on ribose 2'-hydroxyl reactivity in RNA. , 2009, RNA.

[50]  K. Weeks,et al.  Selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE): quantitative RNA structure analysis at single nucleotide resolution , 2006, Nature Protocols.

[51]  K. Weeks.  Advances in RNA structure analysis by chemical probing. , 2010, Current opinion in structural biology.

[52]  K. Weeks,et al.  Detection of RNA-Protein Interactions in Living Cells with SHAPE. , 2015, Biochemistry.

[53]  Giovanni Bussi,et al.  Computer Folding of RNA Tetraloops: Identification of Key Force Field Deficiencies. , 2016, Journal of chemical theory and computation.