The influence of base pair tautomerism on single point mutations in aqueous DNA

The relationship between base pair hydrogen bond proton transfer and the rate of spontaneous single point mutations at ambient temperatures and pressures in aqueous DNA is investigated. By using an ensemble-based multiscale computational modelling method, statistically robust rates of proton transfer for the A:T and G:C base pairs within a solvated DNA dodecamer are calculated. Several different proton transfer pathways are observed within the same base pair. It is shown that, in G:C, the double proton transfer tautomer is preferred, while the single proton transfer process is favoured in A:T. The reported range of rate coefficients for double proton transfer is consistent with recent experimental data. Notwithstanding the approximately 1000 times more common presence of single proton transfer products from A:T, observationally there is bias towards G:C to A:T mutations in a wide range of living organisms. We infer that the double proton transfer reactions between G:C base pairs have a negligible contribution towards this bias for the following reasons: (i) the maximum half-life of the G*:C* tautomer is in the range of picoseconds, which is significantly smaller than the milliseconds it takes for DNA to unwind during replication, (ii) statistically, the majority of G*:C* tautomers revert back to their canonical forms through a barrierless process, and (iii) the thermodynamic instability of the tautomers with respect to the canonical base pairs. Through similar reasoning, we also deduce that proton transfer in the A:T base pair does not contribute to single point mutations in DNA.

AT Base Pair Anions versus (9-Methyl-A)(1-Methyl-T) Base Pair Anions

The anionic base pairs of adenine and thymine, (AT)(-), and 9-methyladenine and 1-methylthymine, (MAMT)(-), have been investigated both theoretically and experimentally in a complementary, synergistic study. Calculations on (AT)(-) found that it had undergone a barrier-free proton transfer (BFPT) similar to that seen in other dimer anion systems and that its structural configuration was neither Watson-Crick (WC) nor Hoogsteen (HS). The vertical detachment energy (VDE) of (AT)(-) was determined by anion photoelectron spectroscopy and found to be in agreement with the VDE value predicted by theory for the BFPT mechanism. An AT pair in DNA is structurally immobilized into the WC configuration, in part, by being bonded to the sugars of the double helix. This circumstance was mimicked by methylating the sites on both A and T where these sugars would have been tied, viz., 9-methyladenine and 1-methylthymine. Calculations found no BFPT in (MAMT)(-) and a resulting (MAMT)(-) configuration that was either HS or WC, with the configurations differing in stability by ca. 2 kcal/mol. The photoelectron spectrum of (MAMT)(-) occurred at a completely different electron binding energy than had (AT)(-). Moreover, the VDE value of (MAMT)(-) was in agreement with that predicted by theory. The configuration of (MAMT)(-) and its lack of electron-induced proton transfer are inter-related. While there may be other pathways for electron-induced DNA alterations, BFPT in the WC/HS configurations of (AT)(-) is not feasible.