Biopolym. Cell. 2005; 21(1):70-79.
Molecular Biophysics
The molecular mechanism of the spontaneous substitution mutations caused by tautomerism of bases: Post Hartree-Fock study of the DNA rare base pairs
1Danilov V. I., 1Hovorun D. M., 2Kurita N.
  1. Institute of Molecular Biology and Genetics, NAS of Ukraine
    150, Akademika Zabolotnoho Str., Kyiv, Ukraine, 03680
  2. Toyohashi University of Technology
    1-1 Hibarigaoka Tenpaku-cho, Toyohashi, Aichi 441-8580


Gas-phase gradient optimization of the DNA rare base pairs containing lactam-lactim and amino-imino tautomers was carried oat using the Hartree-Fock (HF), Density Functional Theory (DFT) and the second-order Moller-Plesset perturbation (MP2) methods at the 6-31G(d, p) basis set. It is shown that full geometry optimization at the MP2 level leads to an intrinsically nonplanar propeller-twisted and buckled geometry of G*-T and G-T* base pairs. The nonplanarity of the pairs is caused by pyramidalization of the amino nitrogen atoms, which is underestimated by the HF and DFT methods. This justifies the importance of geometry optimization at the MP2 level for obtaining reliable prediction of the charge distribution, molecular dipole moments and geometrical structure of the base pairs. The comparison of the formation energies for the rare base pairs shows the energetical preference of the G*-T and A-C* base pairs as compared with the G-T* and A*-C ones, respectively. It is detected that the stabilization energies of the G-T* and A*-C base pairs describing the interaction between monomers are essentially larger than those of the G*-T and A-C* base pairs, respectively. An analysis of the decomposition members for molecular HF interaction energies by Morokuma-Kitaura (MK) and the Reduced Variational Space (RVS) methods showed that the nature of a larger stability of the G-T* and A*-C base pairs as compared to the G*-T and A-C* ones is due to the electrostatic interactions by 60–65 % and the polarization and charge transfer interactions by 35–40 %.
Keywords: spontaneous mutations, tautomerism, rare base pairs, Hartree-Fock, DFT, Moller-Plesset, geometry optimization, nonplanarity of amino group, propeller-twisted and buckled geometry, Morokuma-Kitaura analysis, interaction energy, stability energy


[1] Watson JD, Crick FH. The structure of DNA. Cold Spring Harb Symp Quant Biol. 1953;18:123-31.
[2] Lowdin PO. Quantum genetics and the aperiodic solid. Some aspects of the biological problems of heredity, mutations, aging and tumors in view of the quantum theory of the DNA molecule. Adv. Quant. Chem. 1965; 2:213-360.
[3] Goodman MF. DNA models. Mutations caught in the act. Nature. 1995;378(6554):237-8.
[4] Jiang SP, Raghunathan G, Ting KL, Xuan JC, Jernigan RL. Geometries, charges, dipole moments and interaction energies of normal, tautomeric and novel bases. J Biomol Struct Dyn. 1994;12(2):367-82.
[5] Krechkivska OM, Kosach DA, Sudakov OO, Hovorun DM. Is the complementarity of canonical nucleotide bases in DNA their inherent physical-chemical characteristics? Results of the non-empirical quantum-chemical investigation. Biopolym Cell. 2003; 19(4):382-5.
[6] Hehre WJ, Radom L, Schleyer PVR, Pople JA. Ab initio molecular orbital theory. New York: J. Wiley, 1986.
[7] Danilov VI, Anisimov VM. Post Hartree-Fock studies of the canonical Watson-Crick DNA base pairs: molecular structure and the nature of stability. J Biomol Struct Dyn. 2005;22(4):471-82.
[8] Chalasinski G, Szczesniak MM. Origins of Structure and Energetics of van der Waals Clusters from ab Initio Calculations. Chem Rev.1994;94(7):1723–65.
[9] M?ller C, Plesset MS. Note on an approximation treatment for many-electron systems. Phys Rev. 1934;46(7):618–22.
[10] Morokuma K, Kitaura K. Energy decomposition analysis. Chemical applications of atomic and molecular electrostatic potentials. Eds P. Politzer, D. G. Truhlar. New York: Plenum press, 1981: 215-42.
[11] Stevens WJ, Fink WH. Frozen fragment reduced variational space analysis of hydrogen bonding interactions. Application to the water dimer. Chem Phys Lett. 1987;139(1):15–22.
[12] Gordon MS, Jensen JH. Wavefunctions and chemical bonding: interpretation. The Encyclopedia of Computational Chemistry. New York: J. Wiley, 1998: 3198-3214.
[13] Schmidt MW, Baldridge KK, Boatz JA, Elbert ST, Gordon MS, Jensen JH, et al. General atomic and molecular electronic structure system. J Comp Chem. 1993;14(11):1347–63.
[14] Hobza P, Sponer J. Toward true DNA base-stacking energies: MP2, CCSD(T), and complete basis set calculations. J Am Chem Soc. 2002;124(39):11802-8.
[15] Sponer J, Leszczynski J, Hobza P. Electronic properties, hydrogen bonding, stacking, and cation binding of DNA and RNA bases. Biopolymers. 2001-2002;61(1):3-31.
[16] Diekmann S. Definitions and nomenclature of nucleic acid structure parameters. J Mol Biol. 1989;205(4):787-91.
[17] Komarov VM, Polozov RV, Konoplev GG. Non-planar structure of nitrous bases and non-coplanarity of Watson-Crick pairs. J Theor Biol. 1992;155(3):281-94.
[18] Sloane DL, Goodman MF, Echols H. The fidelity of base selection by the polymerase subunit of DNA polymerase III holoenzyme. Nucleic Acids Res. 1988;16(14A):6465-75.
[19] Echols H, Goodman MF. Fidelity mechanisms in DNA replication. Annu Rev Biochem. 1991;60:477-511.
[20] Goodman MF. Hydrogen bonding revisited: geometric selection as a principal determinant of DNA replication fidelity. Proc Natl Acad Sci U S A. 1997;94(20):10493-5.
[21] Topal MD, Fresco JR. Complementary base pairing and the origin of substitution mutations. Nature. 1976;263(5575):285-9.