A molecular dynamics simulation of the flavin mononucleotide-RNA aptamer complex.

We report on an unrestrained molecular dynamics simulation of the flavin mononucleotide (FMN)-RNA aptamer. The simulated average structure maintains both cross-strand and intermolecular FMN-RNA nuclear Overhauser effects from the nmr experiments and has all qualitative features of the nmr structure including the G10-U12-A25 base triple and the A13-G24, A8-G28, and G9-G27 mismatches. However, the relative orientation of the hairpin loop to the remaining part of the molecule differs from the nmr structure. The simulation predicts that the flexible phosphoglycerol part of FMN moves toward G27 and forms hydrogen bonds. There are structurally long-lived water molecules in the FMN binding pocket forming hydrogen bonds within FMN and between FMN and RNA. In addition, long-lived water is found bridging primarily RNA backbone atoms. A general feature of the environment of long-lived "structural" water is at least two and in most cases three or four potential acceptor atoms. The 2'-OH group of RNA usually acts as an acceptor in interactions with the solvent. There are almost no intrastrand O2'H(n) vertical ...O4'(n + 1) hydrogen bonds within the RNA backbone. In the standard case the preferred orientation of the 2'-OH hydrogen atoms is approximately toward O3' of the same nucleotide. However, a relatively large number of conformations with the backbone torsional angle gamma in the trans orientation is found. A survey of all experimental RNA x-ray structures shows that this backbone conformation occurs but is less frequent than found in the simulation. Experimental nmr RNA aptamer structures have a higher fraction of this conformation as compared to the x-ray structures. The backbone conformation of nucleotide n + 1 with the torsional angle gamma in the trans orientation leads to a relatively short distance between 2'-OH(n) and O5'(n + 1), enabling hydrogen-bond formation. In this case the preferred orientation of the 2'-OH hydrogen atom is approximately toward O5'(n + 1). We find two relatively short and dynamically stable types of backbone-backbone next-neighbor contacts, namely C2'(H)(n) vertical ...O4'(n + 1) and C5'(H)(n + 1) vertical ...O2'(n). These interactions may affect both backbone rigidity and thermodynamic stability of RNA helical structures.

[1]  Peter A. Kollman,et al.  Molecular dynamics simulations highlight the structural differences among DNA: DNA, RNA:RNA, and DNA:RNA hybrid duplexes , 1997 .

[2]  G. Varani.  Exceptionally stable nucleic acid hairpins. , 1995, Annual review of biophysics and biomolecular structure.

[3]  E. Westhof,et al.  Molecular dynamics investigations of hammerhead ribozyme RNA , 1998, European Biophysics Journal.

[4]  I. Tinoco,et al.  The structure of the L3 loop from the hepatitis delta virus ribozyme: a syn cytidine. , 1998, Nucleic acids research.

[5]  W. Olson.  Syn–Anti effects on the spatial configuration of polynucleotide chains , 1973, Biopolymers.

[6]  J. Feigon,et al.  Solution structure of an ATP-binding RNA aptamer reveals a novel fold. , 1997, RNA.

[7]  M. Billeter,et al.  MOLMOL: a program for display and analysis of macromolecular structures. , 1996, Journal of molecular graphics.

[8]  J. Sühnel,et al.  C–H...O and C–H...N interactions in RNA structures , 1999 .

[9]  G. Varani,et al.  Structure of an unusually stable RNA hairpin. , 1991, Biochemistry.

[10]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[11]  T. Steitz,et al.  Metals, Motifs, and Recognition in the Crystal Structure of a 5S rRNA Domain , 1997, Cell.

[12]  J A McCammon,et al.  Phenylalanine transfer RNA: molecular dynamics simulation. , 1984, Science.

[13]  L. Kay,et al.  α Helix-RNA Major Groove Recognition in an HIV-1 Rev Peptide-RRE RNA Complex , 1996, Science.

[14]  L. Gold,et al.  Systematic evolution of ligands by exponential enrichment: RNA ligands to bacteriophage T4 DNA polymerase. , 1990, Science.

[15]  A. R. Srinivasan,et al.  The nucleic acid database. A comprehensive relational database of three-dimensional structures of nucleic acids. , 1992, Biophysical journal.

[16]  Tom Brown,et al.  The orientation and dynamics of the C2'-OH and hydration of RNA and DNA.RNA hybrids , 1998, Nucleic Acids Res..

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

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

[19]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

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

[21]  Eric Westhof,et al.  MULTIPLE MOLECULAR DYNAMICS SIMULATIONS OF THE ANTICODON LOOP OF TRNAASP IN AQUEOUS SOLUTION WITH COUNTERIONS , 1995 .

[22]  C. Sich,et al.  Structure of an RNA hairpin loop with a 5'-CGUUUCG-3' loop motif by heteronuclear NMR spectroscopy and distance geometry. , 1997, Biochemistry.

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

[24]  I. Tinoco,et al.  Crystal structure of an RNA double helix incorporating a track of non-Watson–Crick base pairs , 1991, Nature.

[25]  J. Szostak,et al.  In vitro selection of RNA molecules that bind specific ligands , 1990, Nature.

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

[27]  D. Patel,et al.  Deep penetration of an α-helix into a widened RNA major groove in the HIV-1 rev peptide–RNA aptamer complex , 1996, Nature Structural Biology.

[28]  G. Stormo,et al.  CUUCGG hairpins: extraordinarily stable RNA secondary structures associated with various biochemical processes. , 1988, Proceedings of the National Academy of Sciences of the United States of America.

[29]  D. Patel,et al.  RNA folding topology and intermolecular contacts in the AMP-RNA aptamer complex. , 1996, Biochemistry.

[30]  E Westhof,et al.  Isoalloxazine derivatives promote photocleavage of natural RNAs at G.U base pairs embedded within helices. , 1997, Nucleic acids research.

[31]  E. Kool,et al.  Origins of the large differences in stability of DNA and RNA helices: C-5 methyl and 2'-hydroxyl effects. , 1995, Biochemistry.

[32]  E. Westhof,et al.  Exploration of metal ion binding sites in RNA folds by Brownian-dynamics simulations. , 1998, Structure.

[33]  C R Woese,et al.  Architecture of ribosomal RNA: constraints on the sequence of "tetra-loops". , 1990, Proceedings of the National Academy of Sciences of the United States of America.

[34]  E. Westhof,et al.  Hydration of C-H groups in tRNA. , 1996, Faraday discussions.

[35]  H. Berman,et al.  Geometric Parameters in Nucleic Acids: Sugar and Phosphate Constituents , 1996 .

[36]  P. Burgstaller,et al.  Structural probing and damage selection of citrulline- and arginine-specific RNA aptamers identify base positions required for binding. , 1995, Nucleic acids research.

[37]  D. Patel,et al.  Saccharide-RNA recognition in an aminoglycoside antibiotic-RNA aptamer complex. , 1997, Chemistry & biology.

[38]  Grant R. Zimmermann,et al.  Interlocking structural motifs mediate molecular discrimination by a theophylline-binding RNA , 1997, Nature Structural Biology.

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

[40]  M. Sundaralingam,et al.  C-H...O hydrogen bonding in biology. , 1997, Trends in biochemical sciences.

[41]  P. Kraulis.  A program to produce both detailed and schematic plots of protein structures , 1991 .

[42]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[43]  Alexander D. MacKerell.  Influence of Magnesium Ions on Duplex DNA Structural, Dynamic, and Solvation Properties , 1997 .

[44]  E Westhof,et al.  Water: an integral part of nucleic acid structure. , 1988, Annual review of biophysics and biophysical chemistry.

[45]  P A Kollman,et al.  Theoretical studies of an exceptionally stable RNA tetraloop: observation of convergence from an incorrect NMR structure to the correct one using unrestrained molecular dynamics. , 1997, Journal of molecular biology.

[46]  D. Patel,et al.  Molecular recognition in the FMN-RNA aptamer complex. , 1996, Journal of molecular biology.

[47]  E. Westhof,et al.  RNA hydration: three nanoseconds of multiple molecular dynamics simulations of the solvated tRNA(Asp) anticodon hairpin. , 1997, Journal of molecular biology.

[48]  E Westhof,et al.  Structural Basis of Ligand Discrimination by Two Related RNA Aptamers Resolved by NMR Spectroscopy , 1996, Science.

[49]  T Hermann,et al.  Determination of nucleotide distances in RNA by means of copper phenanthroline-generated hydroxyl radical cleavage pattern. , 1995, RNA.

[50]  P. Burgstaller,et al.  Isolation of RNA Aptamers for Biological Cofactors by In Vitro Selection , 1994 .

[51]  E. Westhof,et al.  Molecular Dynamics Simulations of the Anticodon Hairpin of tRNAAsp: Structuring Effects of C−H···O Hydrogen Bonds and of Long-Range Hydration Forces , 1996 .

[52]  N. Usman,et al.  RNA hydration: a detailed look. , 1996, Biochemistry.

[53]  D. Zichi.  Molecular Dynamics of RNA with the OPLS Force Field. Aqueous Simulation of a Hairpin Containing a Tetranucleotide Loop , 1995 .

[54]  T. Darden,et al.  The effect of long‐range electrostatic interactions in simulations of macromolecular crystals: A comparison of the Ewald and truncated list methods , 1993 .

[55]  J. Åqvist,et al.  Ion-water interaction potentials derived from free energy perturbation simulations , 1990 .