Theoretical investigation of hydrogen transfer mechanism in the guanine–cytosine base pair

A review on point mutations via proton transfer in DNA base pairs in the absence and presence of electric fields

This comprehensive review focuses on spontaneous mutations that may occur during DNA replication, the fundamental process responsible for transferring genetic information. In 1963, Löwdin postulated that these mutations are primarily a result of proton transfer reactions within the hydrogen-bonded DNA base pairs. The single and double proton transfer reactions within the base pairs in DNA result in zwitterions and rare tautomers, respectively. For persistent mutations, these products must be generated at high rates and should be thermodynamically stable. This review covers the proton transfer reactions studied experimentally and computationally. The review also examines the influence of externally applied electric fields on the thermodynamics and kinetics of proton transfer reactions within DNA base pairs, and their biological implications.

Proton transfer during DNA strand separation as a source of mutagenic guanine-cytosine tautomers

Proton transfer between the DNA bases can lead to mutagenic Guanine-Cytosine tautomers. Over the past several decades, a heated debate has emerged over the biological impact of tautomeric forms. Here, we determine that the energy required for generating tautomers radically changes during the separation of double-stranded DNA. Density Functional Theory calculations indicate that the double proton transfer in Guanine-Cytosine follows a sequential, step-like mechanism where the reaction barrier increases quasi-linearly with strand separation. These results point to increased stability of the tautomer when the DNA strands unzip as they enter the helicase, effectively trapping the tautomer population. In addition, molecular dynamics simulations indicate that the relevant strand separation time is two orders of magnitude quicker than previously thought. Our results demonstrate that the unwinding of DNA by the helicase could simultaneously slow the formation but significantly enhance the stability of tautomeric base pairs and provide a feasible pathway for spontaneous DNA mutations.

A theoretical study on the excited-state deactivation paths for the A-5FU dimer

The photostability of DNA plays a key role in the normal function of organisms. A-5FU is a base pair derivative of the A-T dimer where the methyl group is replaced by a F atom. Here, accurate static TDDFT calculations and non-adiabatic dynamic simulations are used to systematically investigate the excited-state decay paths of the A-5FU dimer related to the proton transfer and the out-of-plane twisting deformation motion of A and 5FU in the 1ππ* and 1nπ* states. CC2 is used to check the accuracy of the current TDDFT calculations. Our results show that the deformation of the C[double bond, length as m-dash]C or C[double bond, length as m-dash]N double bond in A and 5FU provides an efficient pathway for the depopulation of the lowest excited states, which can compete with the excited-state proton transfer paths in the dimer. This finding indicates that monomer-like decay paths could be important for the photostability of weakly hydrogen-bonded DNA base pairs and provide a new insight into the excited-state decay paths in base pairs and their analogues.

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.

Proton transfer in guanine-cytosine base pair analogues studied by NMR spectroscopy and PIMD simulations

It has been hypothesised that proton tunnelling between paired nucleobases significantly enhances the formation of rare tautomeric forms and hence leads to errors in DNA replication. Here, we study nuclear quantum effects (NQEs) using deuterium isotope-induced changes of nitrogen NMR chemical shifts in a model base pair consisting of two tautomers of isocytosine, which form hydrogen-bonded dimers in the same way as the guanine-cytosine base pair. Isotope effects in NMR are consequences of NQEs, because ro-vibrational averaging of different isotopologues gives rise to different magnetic shielding of the nuclei. The experimental deuterium-induced chemical shift changes are compared with those calculated by a combination of path integral molecular dynamics (PIMD) simulations with DFT calculations of nuclear shielding. These calculations can directly link the observable isotope-induced shifts with NQEs. A comparison of the deuterium-induced changes of 15N chemical shifts with those predicted by PIMD simulations shows that inter-base proton transfer reactions do not take place in this system. We demonstrate, however, that NMR isotope shifts provide a unique possibility to study NQEs and to evaluate the accuracy of the computational methods used for modelling quantum effects in molecules. Calculations based on the PBE functional from the general-gradient-approximation family provided significantly worse predictions of deuterium isotope shifts than those with the hybrid B3LYP functional.

Influence of Coupling and Embedding Schemes on QM Size Convergence in QM/MM Approaches for the Example of a Proton Transfer in DNA

The influence of embedding and coupling schemes on the convergence of the QM size in the QM/MM approach is investigated for the transfer of a proton in a DNA base pair. We find that the embedding scheme (mechanical or electrostatic) has a much greater impact on the convergence behavior than the coupling scheme (additive QM/MM or subtractive ONIOM). To achieve size convergence, QM regions with up to 6000 atoms are necessary for pure QM or mechanical embedding. In contrast, electrostatic embedding converges faster: for the example of the transfer of a proton between DNA base pairs, we recommend including at least five base pairs and 5 Å of solvent (including counterions) into the QM region, i.e., a total of 1150 atoms.

Exposing the G-quadruplex to electric fields: the role played by telomeres in the propagation of DNA errors

We use quantum calculations to assess the impact of external electric fields on the stability of G-quadruplex, a key structure in telomere functionality.

Assessing the importance of proton transfer reactions in DNA

Although engineered by millions of years of evolution, the cellular machinery is not flawless, and errors regularly appear during DNA replication. The subsequent alteration of the stored genetic message results in a mutation and might be the starting point of important health disorders. The question therefore is what causes DNA mutations? All living organisms are constantly exposed to a number of external agents such as free radicals and to radiation, which may lead to induced mutations. There are also mutations happening without invoking the action of any exogenous element, the so-called spontaneous mutations. The former can be partially controlled by avoiding exposure to high-risk environments, while the latter are more intriguing because their origin is unclear and difficult to determine. As noted by Watson and Crick when they first discovered the DNA structure, the correct replication of DNA rests on the assumption that the base pairs remain in their most stable, canonical form. However, protons along the interbase hydrogen-bond network are not static entities. They can in fact interchange their positions in DNA bases through proton transfer (PT) reactions before strands unwind, giving rise to noncanonical structures defined as rare tautomers. The importance of these rare tautomers was also cleverly anticipated by Watson and Crick and some years later claimed by Löwdin to be a source of spontaneous mutations. In Watson and Crick's words: "It would be of interest to know the precise difference in free energy between the various tautomeric forms under physiological conditions." Unfortunately, rare tautomeric forms are very difficult to detect, so no direct and accurate free energy measure has been discerned. In contrast, theoretical chemistry is making good progress toward the quantification of PT reactions in DNA and their biological consequences. This Account touches upon the theoretical studies devoted to appraising the importance of rare tautomers as promoters of spontaneous mutations. We focus in particular on the crucial role played by the biological environment on DNA stability. It has now been demonstrated that valuable macroscopic predictions require not only highly accurate theories but also refined chemical models. Hybrid quantum mechanics/molecular mechanics (QM/MM) simulations performed on short but complete DNA sequence fragments emerge in this context as the most adequate tools. In addition, these methods can be used to quantify the effect of different external agents on the PT tautomeric equilibria and, eventually, to conveniently handle them. This is the case for the possible alteration of the naturally observed mutation rate by exposure to intense electric fields. Theoretical predictions envision in this respect promising applications of ultrashort electric pulses in medicine to selectively modify the mutated/canonical ratio in DNA.

Coupling between hydrogen atoms transfer and stacking interaction in adenine-thymine/guanine-cytosine complexes: a theoretical study

Four different complexes of two base pairs, an adenine-thymine and a guanine-cytosine one, have been studied in order to understand the modifications induced by the staking interaction between the two base pairs on the hydrogen atoms transfers between the bases in either base pair. The inclusion of these two kinds of interactions allows us to clarify if some properties, as the mechanism of hydrogen transfer, is exclusively a local effect of a base pair or can be modified by a more long-range interaction between the base pairs. The results on these four complexes are compared with those of the monomeric systems, the A-T and G-C base pair, and with those of the A-T and G-C dimers. The specificity of each complex and of each hydrogen bond has been analyzed.

Hydrogen-bonded proton transfer in the hydrated adenine–thymine anion

The proton transfer processes of microhydrated adenine–thymine anions are studied using density functional theory with the B3LYP method and DZP++ basis set. The microhydration effects on the geometrical structures, adsorption site, and the proton transfer reaction of the adenine–thymine anion are investigated. The site N10 atom of the adenine moiety has a larger proton affinity than the site O24 atom of thymine, which facilitates the proton H26 transfers from the N25 site of thymine to the N10 site of adenine. Therefore, the first single-proton transfer pathway (SPT1) is found for the all of the monohydrated adenine–thymine anions (AN4T)−·H2O, (AN13T)−·H2O, (ATO24)−·H2O, and (ATO28)−·H2O and tetrahydrated adenine–thymine anions (AT)−·4H2O. The proton H9 at the N7 site of adenine is also found to transfer to the O24 site of thymine for (AN4T)−·H2O and (AN13T)−·H2O in the gas phase. The double-proton transferred pathway is found when one water molecule interacts with the O28 atom of thymine. The reactant structures before the proton transfer are more stable than the product structures, and the structural changes mainly occur in thymine. The reaction energies of the microhydrated adenine–thymine anion in the gas phase and in the aqueous environment predict that the proton transfer process of the microhydrated adenine–thymine anion are more favorable in the gas phase than in aqueous solution.

Theoretical investigation of the coupling between hydrogen-atom transfer and stacking interaction in adenine-thymine dimers

Three different dimers of the adenine-thymine (A-T) base pair are studied to point out the changes of important properties (structure, atomic charge, energy and so on) induced by coupling between the movement of the atoms in the hydrogen bonds and the stacking interaction. The comparison of these results with those for the A-T monomer system explains the role of the stacking interaction in the hydrogen-atom transfer in this biologically important base pair. The results support the idea that this coupling depends on the exact dimer considered and is different for the N-N and N-O hydrogen bonds. In particular, the correlation between the hydrogen transfer and the stacking interaction is more relevant for the N-N bridge than for the N-O one. Also, the two different mechanisms of two-hydrogen transfer (step by step and concerted) can be modified by the stacking interaction between the base pairs.

Enol tautomers of Watson-Crick base pair models are metastable because of nuclear quantum effects

Intermolecular enol tautomers of Watson-Crick base pairs could emerge spontaneously via interbase double proton transfer. It has been hypothesized that their formation could be facilitated by thermal fluctuations and proton tunneling, and possibly be relevant to DNA damage. Theoretical and computational studies, assuming classical nuclei, have confirmed the dynamic stability of these rare tautomers. However, by accounting for nuclear quantum effects explicitly through Car-Parrinello path integral molecular dynamics calculations, we find the tautomeric enol form to be dynamically metastable, with lifetimes too insignificant to be implicated in DNA damage.