Constructing RNA dynamical ensembles by combining MD and motionally decoupled NMR RDCs: new insights into RNA dynamics and adaptive ligand recognition

We describe a strategy for constructing atomic resolution dynamical ensembles of RNA molecules, spanning up to millisecond timescales, that combines molecular dynamics (MD) simulations with NMR residual dipolar couplings (RDC) measured in elongated RNA. The ensembles are generated via a Monte Carlo procedure by selecting snap-shot from an MD trajectory that reproduce experimentally measured RDCs. Using this approach, we construct ensembles for two variants of the transactivation response element (TAR) containing three (HIV-1) and two (HIV-2) nucleotide bulges. The HIV-1 TAR ensemble reveals significant mobility in bulge residues C24 and U25 and to a lesser extent U23 and neighboring helical residue A22 that give rise to large amplitude spatially correlated twisting and bending helical motions. Omission of bulge residue C24 in HIV-2 TAR leads to a significant reduction in both the local mobility in and around the bulge and amplitude of inter-helical bending motions. In contrast, twisting motions of the helices remain comparable in amplitude to HIV-1 TAR and spatial correlations between them increase significantly. Comparison of the HIV-1 TAR dynamical ensemble and ligand bound TAR conformations reveals that several features of the binding pocket and global conformation are dynamically preformed, providing support for adaptive recognition via a ‘conformational selection’ type mechanism.

[1]  I. Andricioaei,et al.  Impact of static and dynamic A-form heterogeneity on the determination of RNA global structural dynamics using NMR residual dipolar couplings , 2006, Journal of biomolecular NMR.

[2]  Modesto Orozco,et al.  Recent advances in the study of nucleic acid flexibility by molecular dynamics. , 2008, Current opinion in structural biology.

[3]  J. Williamson,et al.  Base flexibility in HIV-2 TAR RNA mapped by solution (15)N, (13)C NMR relaxation. , 2002, Journal of molecular biology.

[4]  D. Patel,et al.  Concerted motions in HIV-1 TAR RNA may allow access to bound state conformations: RNA dynamics from NMR residual dipolar couplings. , 2002, Journal of molecular biology.

[5]  H. Sticht,et al.  Structural Rearrangements of HIV-1 Tat-responsive RNA upon Binding of Neomycin B* , 2000, The Journal of Biological Chemistry.

[6]  A. Szabó,et al.  Model-free approach to the interpretation of nuclear magnetic resonance relaxation in macromolecules. 1. Theory and range of validity , 1982 .

[7]  Gabriele Varani,et al.  Rational design of inhibitors of HIV-1 TAR RNA through the stabilisation of electrostatic "hot spots". , 2004, Journal of molecular biology.

[8]  Zhihua Du,et al.  Structure of TAR RNA complexed with a Tat-TAR interaction nanomolar inhibitor that was identified by computational screening. , 2002, Chemistry & biology.

[9]  J H Prestegard,et al.  Order matrix analysis of residual dipolar couplings using singular value decomposition. , 1999, Journal of magnetic resonance.

[10]  T. Steitz,et al.  A 1.3-A resolution crystal structure of the HIV-1 trans-activation response region RNA stem reveals a metal ion-dependent bulge conformation. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[11]  M. Karplus,et al.  Locally accessible conformations of proteins: Multiple molecular dynamics simulations of crambin , 1998, Protein science : a publication of the Protein Society.

[12]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[13]  J. Williamson Induced fit in RNA–protein recognition , 2000, Nature Structural Biology.

[14]  R Lavery,et al.  The definition of generalized helicoidal parameters and of axis curvature for irregular nucleic acids. , 1988, Journal of biomolecular structure & dynamics.

[15]  Nathan A. Baker,et al.  Iron Responsive Element RNA Flexibility Described by NMR and Isotropic Reorientational Eigenmode Dynamics , 2005, Journal of biomolecular NMR.

[16]  W. Olson,et al.  3DNA: a software package for the analysis, rebuilding and visualization of three-dimensional nucleic acid structures. , 2003, Nucleic acids research.

[17]  A. Bax,et al.  Direct measurement of distances and angles in biomolecules by NMR in a dilute liquid crystalline medium. , 1997, Science.

[18]  Qi Zhang,et al.  Probing motions between equivalent RNA domains using magnetic field induced residual dipolar couplings: accounting for correlations between motions and alignment. , 2003, Journal of the American Chemical Society.

[19]  M. DePristo,et al.  Simultaneous determination of protein structure and dynamics , 2005, Nature.

[20]  Alexander D. MacKerell,et al.  Development and current status of the CHARMM force field for nucleic acids , 2000, Biopolymers.

[21]  R. Lavery,et al.  Defining the structure of irregular nucleic acids: conventions and principles. , 1989, Journal of biomolecular structure & dynamics.

[22]  S. Nosé A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[23]  R. Breaker,et al.  Riboswitches as versatile gene control elements. , 2005, Current opinion in structural biology.

[24]  D. Patel,et al.  Characterizing the relative orientation and dynamics of RNA A-form helices using NMR residual dipolar couplings , 2007, Nature Protocols.

[25]  E. Westhof,et al.  H-bond stability in the tRNA(Asp) anticodon hairpin: 3 ns of multiple molecular dynamics simulations. , 1996, Biophysical journal.

[26]  H. Schwalbe,et al.  NMR and MD studies of the temperature-dependent dynamics of RNA YNMG-tetraloops , 2008, Nucleic acids research.

[27]  K. Hall,et al.  Isotropic reorientational eigenmode dynamics complements NMR relaxation measurements for RNA. , 2005, Methods in enzymology.

[28]  D. Patel,et al.  Evidence that electrostatic interactions dictate the ligand-induced arrest of RNA global flexibility. , 2005, Angewandte Chemie.

[29]  J. Karn,et al.  The structure of the human immunodeficiency virus type-1 TAR RNA reveals principles of RNA recognition by Tat protein. , 1995, Journal of molecular biology.

[30]  Alexander D. MacKerell,et al.  Molecular dynamics simulations of nucleic acid-protein complexes. , 2008, Current opinion in structural biology.

[31]  H. Al‐Hashimi,et al.  Dynamics of large elongated RNA by NMR carbon relaxation. , 2007, Journal of the American Chemical Society.

[32]  J. Karn,et al.  Structure of HIV-1 TAR RNA in the absence of ligands reveals a novel conformation of the trinucleotide bulge. , 1996, Nucleic acids research.

[33]  Qi Zhang,et al.  Resolving the Motional Modes That Code for RNA Adaptation , 2006, Science.

[34]  R. Brüschweiler,et al.  Quantitative molecular ensemble interpretation of NMR dipolar couplings without restraints. , 2007, Journal of the American Chemical Society.

[35]  L. Nilsson,et al.  Cytosine ribose flexibility in DNA: a combined NMR 13C spin relaxation and molecular dynamics simulation study , 2008, Nucleic acids research.

[36]  J. Karn,et al.  Structure-based drug design targeting an inactive RNA conformation: exploiting the flexibility of HIV-1 TAR RNA. , 2004, Journal of molecular biology.

[37]  A. Saupe Recent Results in the Field of Liquid Crystals , 1968 .

[38]  D. Kern,et al.  Dynamic personalities of proteins , 2007, Nature.

[39]  Nikolay V Dokholyan,et al.  Deciphering protein dynamics from NMR data using explicit structure sampling and selection. , 2007, Biophysical journal.

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

[41]  H. Al‐Hashimi,et al.  Probing Na(+)-induced changes in the HIV-1 TAR conformational dynamics using NMR residual dipolar couplings: new insights into the role of counterions and electrostatic interactions in adaptive recognition. , 2007, Biochemistry.

[42]  Charles K. Fisher,et al.  Visualizing spatially correlated dynamics that directs RNA conformational transitions , 2007, Nature.

[43]  M. Karplus,et al.  Deformable stochastic boundaries in molecular dynamics , 1983 .

[44]  Ioan Andricioaei,et al.  Characterizing complex dynamics in the transactivation response element apical loop and motional correlations with the bulge by NMR, molecular dynamics, and mutagenesis. , 2008, Biophysical journal.

[45]  D. Patel,et al.  Towards structural genomics of RNA: rapid NMR resonance assignment and simultaneous RNA tertiary structure determination using residual dipolar couplings. , 2002, Journal of molecular biology.

[46]  Martin Blackledge,et al.  Conformational distributions of unfolded polypeptides from novel NMR techniques. , 2008, The Journal of chemical physics.

[47]  Lewis E. Kay,et al.  New Tools Provide New Insights in NMR Studies of Protein Dynamics , 2006, Science.

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

[49]  G. Varani,et al.  Recent advances in RNA-protein recognition. , 2001, Current opinion in structural biology.

[50]  H. Al‐Hashimi,et al.  RNA dynamics: it is about time. , 2008, Current opinion in structural biology.

[51]  Ioan Andricioaei,et al.  iRED analysis of TAR RNA reveals motional coupling, long-range correlations, and a dynamical hinge. , 2007, Biophysical journal.

[52]  K. Hall,et al.  RNA in motion. , 2008, Current opinion in chemical biology.

[53]  D. Crothers,et al.  Characterization of the solution conformations of unbound and Tat peptide-bound forms of HIV-1 TAR RNA. , 1999, Biochemistry.

[54]  G. Varani,et al.  Current topics in RNA-protein recognition: control of specificity and biological function through induced fit and conformational capture. , 2001, Biochemistry.

[55]  Harald Schwalbe,et al.  Time-resolved NMR studies of RNA folding. , 2007, Biopolymers.

[56]  A. Serganov,et al.  Argininamide binding arrests global motions in HIV-1 TAR RNA: comparison with Mg2+-induced conformational stabilization. , 2004, Journal of molecular biology.

[57]  J. Williamson,et al.  Solution structure of the HIV-2 TAR-argininamide complex. , 1997, Journal of molecular biology.

[58]  H. Al‐Hashimi Dynamics‐Based Amplification of RNA Function and Its Characterization by Using NMR Spectroscopy , 2005, Chembiochem : a European journal of chemical biology.

[59]  J. Puglisi,et al.  Conformation of the TAR RNA-arginine complex by NMR spectroscopy. , 1992, Science.

[60]  Evgeny Nudler,et al.  Flipping Riboswitches , 2006, Cell.

[61]  P. Kadeřávek,et al.  NMR 13C-relaxation Study of Base and Sugar Dynamics in GCAA RNA Hairpin Tetraloop , 2007, Journal of biomolecular structure & dynamics.

[62]  G. Bouvignies,et al.  Exploring multiple timescale motions in protein GB3 using accelerated molecular dynamics and NMR spectroscopy. , 2007, Journal of the American Chemical Society.

[63]  J H Prestegard,et al.  Nuclear magnetic dipole interactions in field-oriented proteins: information for structure determination in solution. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[64]  H. Al‐Hashimi,et al.  Extending the NMR spatial resolution limit for RNA by motional couplings , 2008, Nature Methods.

[65]  Charles K. Fisher,et al.  Ultrahigh resolution characterization of domain motions and correlations by multialignment and multireference residual dipolar coupling NMR. , 2008, The journal of physical chemistry. B.

[66]  Giovanni Lipari,et al.  MODEL-FREE APPROACH TO THE INTERPRETATION OF NUCLEAR MAGNETIC RESONANCE RELAXATION IN MACROMOLECULES. 1. THEORY AND RANGE OF VALIDITY , 1982 .

[67]  D. Patel,et al.  Mg2+-induced variations in the conformation and dynamics of HIV-1 TAR RNA probed using NMR residual dipolar couplings. , 2003, Journal of molecular biology.