Stabilization Energies of the Hydrogen-Bonded and Stacked Structures of Nucleic Acid Base Pairs in the Crystal Geometries of CG, AT, and AC DNA Steps and in the NMR Geometry of the 5‘-d(GCGAAGC)-3‘ Hairpin: Complete Basis Set Calculations at the MP2 and CCSD(T) Levels
- 26 January 2005
- journal article
- Published by American Chemical Society (ACS) in The Journal of Physical Chemistry A
- Vol. 109 (6) , 1131-1136
- https://doi.org/10.1021/jp046738a
Abstract
Stabilization energies of the H-bonded and stacked structures of a DNA base pair were studied in the crystal structures of adenine-thymine, cytosine-guanine, and adenine-cytosine steps as well as in the 5'-d(GCGAAGC)-3' hairpin (utilizing the NMR geometry). Stabilization energies were determined as the sum of the complete basis set (CBS) limit of MP2 stabilization energies and the Delta E(CCSD(T)) - Delta E(MP2) correction term evaluated with the 6-31G*(0.25) basis set. The CBS limit was determined by a two-point extrapolation using the aug-cc-pVXZ basis sets for X = D and T. While the H-bonding energies are comparable to those of base pairs in a crystal and a vacuum, the stacking energies are considerably smaller in a crystal. Despite this, the stacking is still important and accounts for a significant part of the overall stabilization. It contributes equally to the stability of DNA as does H-bonding for AT-rich DNAs, while in the case of GC-rich DNAs it forms about one-third of the total stabilization. Interstrand stacking reaches surprisingly large values, well comparable to the intrastrand ones, and thus contributes significantly to the overall stabilization. The hairpin structure is characterized by significant stacking, and both guanine...cytosine pairs possess stacking energies larger than 11.5 kcal/mol. A high portion of stabilization in the studied hairpin comes from stacking (similar to that found for AT-rich DNAs) despite the fact that it contains two GC Watson-Crick pairs having very large H-bonding stabilization. The DFT/B3LYP/6-31G** method yields satisfactory values of interaction energies for H-bonded structures, while it fails completely for stacking.Keywords
This publication has 21 references indexed in Scilit:
- 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 ExperimentJournal of the American Chemical Society, 2003
- Solvent Organization in an Oligonucleotide CrystalJournal of Biological Chemistry, 2000
- Noncovalent Interactions: A Challenge for Experiment and TheoryChemical Reviews, 1999
- Preorganization of DNA: Design Principles for Improving Nucleic Acid Recognition by Synthetic OligonucleotidesChemical Reviews, 1997
- Experimental Measurement of Aromatic Stacking Affinities in the Context of Duplex DNAJournal of the American Chemical Society, 1996
- Nature of Nucleic Acid−Base Stacking: Nonempirical ab Initio and Empirical Potential Characterization of 10 Stacked Base Dimers. Comparison of Stacked and H-Bonded Base PairsThe Journal of Physical Chemistry, 1996
- Theoretical Investigation of the Hydrogen Bond Strengths in Guanine-Cytosine and Adenine-Thymine Base PairsJournal of the American Chemical Society, 1994
- Structure of the B-DNA decamer C-C-A-A-C-G-T-T-G-G and comparison with isomorphous decamers C-C-A-A-G-A-T-T-G-G and C-C-A-G-G-C-C-T-G-GJournal of Molecular Biology, 1991
- Stability of XGCGCp, GCGCYp, and XGCGCYp helixes: an empirical estimate of the energetics of hydrogen bonds in nucleic acidsBiochemistry, 1986
- Experimental studies of molecular interactions between nitrogen bases of nucleic acidsBiopolymers, 1979