Are the hydrogen bonds of RNA (AU) stronger than those of DNA (AT)? A quantum mechanics study.

The intrinsic stability of Watson-Crick d(AT) and r(AU) hydrogen bonds was analyzed by employing a variety of quantum-mechanical techniques, such as energy calculations, determination of reactivity indexes, and analysis of electron density topology. The analyses were performed not only for equilibrium gas-phase geometries, but also on hundreds of conformations derived from molecular dynamics (MD) and database analysis. None of our results support the idea that r(AU) hydrogen bonds are intrinsically more stable than those of d(AT). Instead, our data are in accordance with the traditional view that the greater stability of RNA relative to DNA is attributable to a variety of effects (e.g., stacking, sugar puckering, solvation) rather than to a significant difference in the hydrogen bonding of DNA and RNA base pairs.

[1]  José Elguero,et al.  Non-conventional hydrogen bonds , 1998 .

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

[3]  F. Javier Luque,et al.  Polarization effects in generalized molecular interaction potential: New Hamiltonian for reactivity studies and mixed QM/MM calculations , 1998 .

[4]  R. Bader Atoms in molecules : a quantum theory , 1990 .

[5]  F. J. Luque,et al.  Is polarization important in cation-pi interactions? , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[6]  Trygve Helgaker,et al.  A priori calculation of molecular properties to chemical accuracy , 2004 .

[7]  Michael B Hursthouse,et al.  Halogen bonding: a new interaction for liquid crystal formation. , 2004, Journal of the American Chemical Society.

[8]  M. Riley,et al.  Physical and chemical characterization of two- and three-stranded adenine-thymine and adenine-uracil homopolymer complexes. , 1966, Journal of molecular biology.

[9]  Pavel Hobza,et al.  Interactions of Hydrated Mg2+ Cation with Bases, Base Pairs, and Nucleotides. Electron Topology, Natural Bond Orbital, Electrostatic, and Vibrational Study , 2001 .

[10]  I. Tinoco,et al.  Calorimetric and spectroscopic investigation of the helix-to-coil transition of a ribo-oligonucleotide: rA7U7. , 1975, Journal of molecular biology.

[11]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[12]  A. LiWang,et al.  N1...N3 hydrogen bonds of A:U base pairs of RNA are stronger than those of A:T base pairs of DNA. , 2004, Journal of the American Chemical Society.

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

[14]  Richard F. W. Bader A quantum theory of molecular structure and its applications , 1991 .

[15]  Modesto Orozco,et al.  Relative flexibility of DNA and RNA: a molecular dynamics study. , 2004, Journal of molecular biology.

[16]  Pavel Hobza,et al.  On the convergence of the (ΔECCSD(T)−ΔEMP2) term for complexes with multiple H-bonds , 2002 .

[17]  Pavel Hobza,et al.  True stabilization energies for the optimal planar hydrogen-bonded and stacked structures of guanine...cytosine, adenine...thymine, and their 9- and 1-methyl derivatives: complete basis set calculations at the MP2 and CCSD(T) levels and comparison with experiment. , 2003, Journal of the American Chemical Society.

[18]  F. J. Luque,et al.  The relative flexibility of B-DNA and A-RNA duplexes: database analysis. , 2004, Nucleic acids research.

[19]  S. Acharya,et al.  Measurement of nucleobase pKa values in model mononucleotides shows RNA-RNA duplexes to be more stable than DNA-DNA duplexes. , 2004, Journal of the American Chemical Society.

[20]  Hans W. Horn,et al.  ELECTRONIC STRUCTURE CALCULATIONS ON WORKSTATION COMPUTERS: THE PROGRAM SYSTEM TURBOMOLE , 1989 .

[21]  N. Sugimoto,et al.  Influences of ribonucleotide on a duplex conformation and its thermal stability: study with the chimeric RNA-DNA strands. , 2004, Journal of the American Chemical Society.

[22]  P Hobza,et al.  Structure, energetics, and dynamics of the nucleic Acid base pairs: nonempirical ab initio calculations. , 1999, Chemical reviews.

[23]  J. Elguero,et al.  Bond Length–Electron Density Relationships: From Covalent Bonds to Hydrogen Bond Interactions , 1998 .

[24]  F. Weigend,et al.  RI-MP2: first derivatives and global consistency , 1997 .

[25]  Trygve Helgaker,et al.  Basis-set convergence of the energy in molecular Hartree–Fock calculations , 1999 .

[26]  M. Swart,et al.  Hydrogen bonds of RNA are stronger than those of DNA, but NMR monitors only presence of methyl substituent in uracil/thymine. , 2004, Journal of the American Chemical Society.

[27]  Trygve Helgaker,et al.  Basis-set convergence in correlated calculations on Ne, N2, and H2O , 1998 .

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

[29]  J. G. Snijders,et al.  Hydrogen Bonding in DNA Base Pairs: Reconciliation of Theory and Experiment , 2000 .

[30]  Friedrich Biegler-König,et al.  Calculation of the average properties of atoms in molecules. II , 1982 .