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.
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.
Dabkowska
Iwona
I
Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Center for Complex Molecular Systems and Biomolecules, Flemingovo nam. 2, 166 10 Prague 6, Czech Republic.
Gonzalez
Haydee Valdes
HV
Jurecka
Petr
P
Hobza
Pavel
P
D000225:Adenine / Q000737:chemistry* / Q000378:metabolism; D000465:Algorithms; D020029:Base Pairing; D001483:Base Sequence; D003596:Cytosine / Q000737:chemistry*; D006147:Guanine / Q000737:chemistry*; D006860:Hydrogen Bonding; D009682:Magnetic Resonance Spectroscopy; D008958:Models, Molecular; D008969:Molecular Sequence Data; D009690:Nucleic Acid Conformation; D013816:Thermodynamics; D013941:Thymine / Q000737:chemistry* / Q000378:metabolism
D006147:Guanine; D003596:Cytosine; D000225:Adenine; D013941:Thymine