HIV-1 TAR RNA Spontaneously Undergoes Relevant Apo-to-Holo Conformational Transitions in Molecular Dynamics and Constrained Geometrical Simulations

We report all-atom molecular dynamics and replica exchange molecular dynamics simulations on the unbound human immunodeficiency virus type-1 (HIV-1) transactivation responsive region (TAR) RNA structure and three TAR RNA structures in bound conformations of, in total, approximately 250 ns length. We compare the extent of observed conformational sampling with that of the conceptually simpler and computationally much cheaper constrained geometrical simulation approach framework rigidity optimized dynamic algorithm (FRODA). Atomic fluctuations obtained by replica-exchange molecular dynamics (REMD) simulations agree quantitatively with those obtained by molecular dynamics (MD) and FRODA simulations for the unbound TAR structure. Regarding the stereochemical quality of the generated conformations, backbone torsion angles and puckering modes of the sugar-phosphate backbone were reproduced equally well by MD and REMD simulations, but further improvement is needed in the case of FRODA simulations. Essential dynamics analysis reveals that all three simulation approaches show a tendency to sample bound conformations when starting from the unbound TAR structure, with MD and REMD simulations being superior with respect to FRODA. These results are consistent with the experimental view that bound TAR RNA conformations are transiently sampled in the free ensemble, following a conformation selection model. The simulation-generated TAR RNA conformations have been successfully used as receptor structures for docking. This finding has important implications for RNA-ligand docking in that docking into an ensemble of simulation-generated RNA structures is shown to be a valuable means to cope with large apo-to-holo conformational transitions of the receptor structure.

[1]  G. Bodenhausen,et al.  Cross correlations between 13C-1H dipolar interactions and 15N chemical shift anisotropy in nucleic acids , 2003, Journal of biomolecular NMR.

[2]  H. Gohlke,et al.  Determinants of the Unexpected Stability of RNA Fluorobenzene Self Pairs , 2008, Chembiochem : a European journal of chemical biology.

[3]  M. Karplus,et al.  Collective motions in proteins: A covariance analysis of atomic fluctuations in molecular dynamics and normal mode simulations , 1991, Proteins.

[4]  Valentina Tozzini,et al.  Coarse-grained models for proteins. , 2005, Current opinion in structural biology.

[5]  R. Nussinov,et al.  Folding and binding cascades: shifts in energy landscapes. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[6]  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.

[7]  D. Jacobs,et al.  Protein flexibility predictions using graph theory , 2001, Proteins.

[8]  R. Abagyan,et al.  Flexible ligand docking to multiple receptor conformations: a practical alternative. , 2008, Current opinion in structural biology.

[9]  Jacobs,et al.  Generic rigidity percolation: The pebble game. , 1995, Physical review letters.

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

[11]  Holger Gohlke,et al.  Target flexibility: an emerging consideration in drug discovery and design. , 2008, Journal of medicinal chemistry.

[12]  Ly Le,et al.  Explicitly solvated ligand contribution to continuum solvation models for binding free energies: selectivity of theophylline binding to an RNA aptamer. , 2010, The journal of physical chemistry. B.

[13]  X. Zou,et al.  Ensemble docking of multiple protein structures: Considering protein structural variations in molecular docking , 2006, Proteins.

[14]  Claudio N. Cavasotto,et al.  Protein flexibility in ligand docking and virtual screening to protein kinases. , 2004, Journal of molecular biology.

[15]  Viktor Hornak,et al.  HIV-1 protease flaps spontaneously open and reclose in molecular dynamics simulations. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[16]  G. Varani,et al.  Targeting RNA with small-molecule drugs: therapeutic promise and chemical challenges. , 2001, Accounts of chemical research.

[17]  I. Bahar,et al.  Coarse-grained normal mode analysis in structural biology. , 2005, Current opinion in structural biology.

[18]  P. Kollman,et al.  A modified version of the Cornell et al. force field with improved sugar pucker phases and helical repeat. , 1999, Journal of biomolecular structure & dynamics.

[19]  Markus Wagener,et al.  A flexible approach to induced fit docking. , 2007, Journal of medicinal chemistry.

[20]  V. Helms,et al.  Transient pockets on protein surfaces involved in protein-protein interaction. , 2007, Journal of medicinal chemistry.

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

[22]  Petra Fromme,et al.  Fitting low-resolution cryo-EM maps of proteins using constrained geometric simulations. , 2008, Biophysical journal.

[23]  Claudio N. Cavasotto,et al.  Representing receptor flexibility in ligand docking through relevant normal modes. , 2005, Journal of the American Chemical Society.

[24]  Bryan R. Cullen,et al.  Trans-activation of human immunodeficiency virus occurs via a bimodal mechanism , 1986, Cell.

[25]  Holger Gohlke,et al.  Elastic Potential Grids: Accurate and Efficient Representation of Intermolecular Interactions for Fully Flexible Docking , 2009, ChemMedChem.

[26]  M Feig,et al.  Conformations of an adenine bulge in a DNA octamer and its influence on DNA structure from molecular dynamics simulations. , 2001, Biophysical journal.

[27]  Holger Gohlke,et al.  Analyzing the flexibility of RNA structures by constraint counting. , 2008, Biophysical journal.

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

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

[30]  Holger Gohlke,et al.  Change in protein flexibility upon complex formation: Analysis of Ras‐Raf using molecular dynamics and a molecular framework approach , 2004, Proteins.

[31]  Nils G Walter,et al.  Molecular dynamics simulations of RNA: an in silico single molecule approach. , 2007, Biopolymers.

[32]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[33]  Aqeel Ahmed,et al.  Protein Flexibility and Mobility in Structure-Based Drug Design , 2007 .

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

[35]  Holger Gohlke,et al.  DrugScoreRNAKnowledge-Based Scoring Function To Predict RNA-Ligand Interactions , 2007, J. Chem. Inf. Model..

[36]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[37]  Stephen Hanessian,et al.  Docking of aminoglycosides to hydrated and flexible RNA. , 2006, Journal of medicinal chemistry.

[38]  Richard J. Hall,et al.  Protein-Ligand Docking against Non-Native Protein Conformers , 2008, J. Chem. Inf. Model..

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

[40]  J. Berg,et al.  Molecular dynamics simulations of biomolecules , 2002, Nature Structural Biology.

[41]  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.

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

[43]  Ming Lei,et al.  Sampling protein conformations and pathways , 2004, J. Comput. Chem..

[44]  M. Zacharias,et al.  Conformational transitions in RNA single uridine and adenosine bulge structures: a molecular dynamics free energy simulation study. , 2006, Biophysical journal.

[45]  G. Varani,et al.  Intermediate rate atomic trajectories of RNA by solid-state NMR spectroscopy. , 2010, Journal of the American Chemical Society.

[46]  Wolfram Saenger,et al.  Principles of Nucleic Acid Structure , 1983 .

[47]  Gabriele Cruciani,et al.  Targeting the Conformational Transitions of MDM2 and MDMX: Insights into Dissimilarities and Similarities of p53 Recognition , 2008, J. Chem. Inf. Model..

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

[49]  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.

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

[51]  H. Gohlke,et al.  Multi-scale modelling of macromolecular conformational changes , 2009 .

[52]  Xian Zhao,et al.  Conformational Transition Map of an RNA GCAA Tetraloop Explored by Replica-Exchange Molecular Dynamics Simulation. , 2009, Journal of chemical theory and computation.

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

[54]  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.

[55]  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.

[56]  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.

[57]  R. Pellicciari,et al.  Molecular docking and spatial coarse graining simulations as tools to investigate substrate recognition, enhancer binding and conformational transitions in indoleamine-2,3-dioxygenase (IDO). , 2007, Biochimica et biophysica acta.

[58]  Thanyada Rungrotmongkol,et al.  Dynamic Behavior of Avian Influenza A Virus Neuraminidase Subtype H5N1 in Complex with Oseltamivir, Zanamivir, Peramivir, and Their Phosphonate Analogues , 2009, J. Chem. Inf. Model..

[59]  G. Klebe,et al.  Expect the unexpected or caveat for drug designers: multiple structure determinations using aldose reductase crystals treated under varying soaking and co-crystallisation conditions. , 2006, Journal of molecular biology.

[60]  D. J. Price,et al.  Assessing scoring functions for protein-ligand interactions. , 2004, Journal of medicinal chemistry.

[61]  Michal Otyepka,et al.  Dependence of A-RNA simulations on the choice of the force field and salt strength. , 2009, Physical chemistry chemical physics : PCCP.

[62]  P A Kollman,et al.  Molecular dynamics studies of the HIV-1 TAR and its complex with argininamide. , 2000, Nucleic acids research.

[63]  D. Turner,et al.  An RNA Molecular Switch: Intrinsic Flexibility of 23S rRNA Helices 40 and 68 5'-UAA/5'-GAN Internal Loops Studied by Molecular Dynamics Methods. , 2010, Journal of chemical theory and computation.

[64]  Holger Gohlke,et al.  A natural coarse graining for simulating large biomolecular motion. , 2006, Biophysical journal.

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

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

[67]  Richard A. Lewis,et al.  Lessons in molecular recognition: the effects of ligand and protein flexibility on molecular docking accuracy. , 2004, Journal of medicinal chemistry.

[68]  J. Šponer,et al.  Refinement of the AMBER Force Field for Nucleic Acids: Improving the Description of α/γ Conformers , 2007 .

[69]  J. Šponer,et al.  Molecular Dynamics and Quantum Mechanics of RNA: Conformational and Chemical Change We Can Believe In , 2009, Accounts of chemical research.

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

[71]  P. Kollman,et al.  Automatic atom type and bond type perception in molecular mechanical calculations. , 2006, Journal of molecular graphics & modelling.

[72]  M. Thorpe,et al.  Constrained geometric simulation of diffusive motion in proteins , 2005, Physical biology.

[73]  David S. Goodsell,et al.  Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function , 1998, J. Comput. Chem..

[74]  Thomas L. James,et al.  Docking to RNA via Root-Mean-Square-Deviation-Driven Energy Minimization with Flexible Ligands and Flexible Targets , 2008, J. Chem. Inf. Model..

[75]  Holger Gohlke,et al.  Constraint counting on RNA structures: linking flexibility and function. , 2009, Methods.

[76]  E Westhof,et al.  Rules governing the orientation of the 2'-hydroxyl group in RNA. , 1997, Journal of molecular biology.

[77]  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.

[78]  M. Zentgraf Characterization of Binding Pocket Flexibility of Aldose Reductase , 2006 .

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

[80]  A. W. Schüttelkopf,et al.  PRODRG: a tool for high-throughput crystallography of protein-ligand complexes. , 2004, Acta crystallographica. Section D, Biological crystallography.

[81]  Urban Bren,et al.  DNA duplex stability: the role of preorganized electrostatics. , 2010, The journal of physical chemistry. B.

[82]  H. Gohlke,et al.  Multiscale modeling of macromolecular conformational changes combining concepts from rigidity and elastic network theory , 2006, Proteins.

[83]  R. Nussinov,et al.  Folding funnels, binding funnels, and protein function , 1999, Protein science : a publication of the Protein Society.

[84]  Holger Gohlke,et al.  Molecular recognition of RNA: challenges for modelling interactions and plasticity , 2009, Journal of molecular recognition : JMR.