Theoretical investigation of hydrogen transfer mechanism in the adenine–thymine 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.

Multiscale simulations reveal the role of PcrA helicase in protecting against spontaneous point mutations in DNA

Proton transfer across hydrogen bonds in DNA can produce non-canonical nucleobase dimers and is a possible source of single-point mutations when these forms mismatch under replication. Previous computational studies have revealed this process to be energetically feasible for the guanine-cytosine (GC) base pair, but the tautomeric product (G[Formula: see text]C[Formula: see text]) is short-lived. In this work we reveal, for the first time, the direct effect of the replisome enzymes on proton transfer, rectifying the shortcomings of existing models. Multi-scale quantum mechanical/molecular dynamics (QM/MM) simulations reveal the effect of the bacterial PcrA Helicase on the double proton transfer in the GC base pair. It is shown that the local protein environment drastically increases the activation and reaction energies for the double proton transfer, modifying the tautomeric equilibrium. We propose a regime in which the proton transfer is dominated by tunnelling, taking place instantaneously and without atomic rearrangement of the local environment. In this paradigm, we can reconcile the metastable nature of the tautomer and show that ensemble averaging methods obscure detail in the reaction profile. Our results highlight the importance of explicit environmental models and suggest that asparagine N624 serves a secondary function of reducing spontaneous mutations in PcrA Helicase.

Tautomerisation Mechanisms in the Adenine-Thymine Nucleobase Pair during DNA Strand Separation

The adenine-thymine tautomer (A*-T*) has previously been discounted as a spontaneous mutagenesis mechanism due to the energetic instability of the tautomeric configuration. We study the stability of A*-T* while the nucleobases undergo DNA strand separation. Our calculations indicate an increase in the stability of A*-T* as the DNA strands unzip and the hydrogen bonds between the bases stretch. Molecular Dynamics simulations reveal the time scales and dynamics of DNA strand separation and the statistical ensemble of opening angles present in a biological environment. Our results demonstrate that the unwinding of DNA, an inherently out-of-equilibrium process facilitated by helicase, will change the energy landscape of the adenine-thymine tautomerization reaction. We propose that DNA strand separation allows the stable tautomerization of adenine-thymine, providing a feasible pathway for genetic point mutations via proton transfer between the A-T bases.

Time delay during intra-base proton tunneling in the guanine base of the single stranded DNA

Spontaneous point mutations are one of the main actors in evolution, and the tautomerization of organic bases in the DNA is hypothesized to be the underlying mechanism of this crucial process. Tautomerization can be induced by proton tunneling, and if it occurs in single-stranded DNA (ssDNA) during the replication process, tautomerized bases might give rise to a mismatch, which will eventually defect the genetic code. In the present work, we report on the tautomerization time in the guanine base of the ssDNA. The model we use includes an intra-base tunneling mechanism such that time tunneling (delay time) is estimated to be around a few pico-seconds. The time delay is found to be biologically relevant which indicates that it is long enough to induce point mutations. Our results close a gap in the literature and sheds light on the importance of point mutations originating from quantum effects in the ssDNA.

Non-poissonian Distribution of Point Mutations in DNA

In general, for chemical reactions occurring in systems, where fluctuations are not negligibly small, it is necessary to introduce a master equation for distribution of probability of fluctuations. It has been established that the monomolecular reactions of a type as A ↔ X are described by the master equation, which leads to a Poisson distribution with the variance equal to the average value N 0. However, the consideration of the Löwdin mechanism as autocatalytic non-linear chemical reactions such as A + X ↔ 2X and the corresponding master equation lead to a non-Poissonian probability distribution of fluctuations. In the presented work, first-order autocatalysis has been applied to the Löwdin's mechanism of spontaneous mutations in DNA. Describing double proton transfers between complimentary nucleotide bases along the chain by first-order autocatalytic reactions, the corresponding master equation for protons in tautomeric states becomes non-linear, and at non-equilibrium conditions this leads to the non-Poissonian distribution of spontaneous mutations in DNA. It is also suggested that the accumulation of large fluctuations of successive cooperative concerted protons along the chain may produce higher non-linearities which could have a significant impact on some biochemical processes, occurring in DNA.

Modeling the action of environment on proton tunneling in the adenine-thymine base pair

Effect of environment coupling on the quantum-biological phenomenon of proton tunneling in the hydrogen bonds of the adenine-thymine base pair in DNA was modeled within the framework of quantum statistics and perturbation theory. A number of important thermodynamic indicators including partition function, free energy, and entropy were then calculated and examined. The proton was then assumed to be subject to an attraction represented by a double-well potential energy surface with a small asymmetry, which was considered as the perturbation introduced to the system. The action of environment manifested itself in the form of a global minimum in the free energy curve, as an implicit implication of the tendency of the system toward randomness and disorder, at which no spontaneous change such as quantum tunneling will accordingly occur. Furthermore, assuming the free energy to be in a close neighborhood of its minimum truly explained the smallness of the contribution of environment coupling to the tunneling probability reported in the literature based on the fact that the closer the free energy to its minimum, the less the transition probability to this point.

Solvent effect on the intermolecular proton transfer of the Watson and Crick guanine-cytosine and adenine-thymine base pairs: a polarizable continuum model study

Herein we present our results on the study of the double proton transfer (DPT) mechanism in the adenine-thymine (AT) and guanine-cytosine (GC) base pairs, both in gas phase and in solution. The latter was modeled using the polarizable continuum method (PCM) in different solvents. According to our DFT calculations, the DPT may occur for both complexes in a stepwise mechanism in condensate phase. In gas phase only the GC base pair exhibits a concerted DPT mechanism. Using the Wigner's tunneling corrections to the transition state theory we demonstrate that such corrections are important for the prediction of the rate constants of both systems in gas and in condensate phase. We also show that (i) as the polarity of the medium decreases the equilibrium constant of the DPT reaction increases in both complexes, and (ii) that the equilibrium constant in the GC complex is four orders of magnitude larger than in AT. This observation suggests that the spontaneous mutations in DNA base pairs are more probable in GC than in AT.

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.

Competition between stacked and hydrogen bonded structures of cytosine aggregates

The four bases of DNA constitute what is known as the “alphabet of life”.

Modelling proton tunnelling in the adenine–thymine base pair

Environment-assisted quantum tunnelling is unlikely to play a significant role in adenine–thymine tautomerisation.

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.

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.

Investigation of C–H…OC and N–H…OC hydrogen-bonding interactions in crystalline thymine by DFT calculations of O-17, N-14 and H-2 NQR parameters

A computational study at the level of density functional theory (DFT) was carried out to investigate C-H...O=C and N-H...O=C hydrogen-bonding interactions (HBs) in the real crystalline cluster of thymine by O-17, N-14 and H-2 calculated nuclear quadrupole resonance (NQR) parameters. To perform the calculations, a hydrogen-bonded pentameric cluster of thymine was created using X-ray coordinates where the hydrogen atoms positions are optimized and the electric field gradient (EFG) tensors were calculated for the target molecule. Additional EFG calculations were also performed for crystalline monomer and an optimized isolated gas-phase thymine. The calculated EFG tensors at the level of B3LYP and B3PW91 DFT methods and 6-311++G**and CC-pVTZ basis sets were converted to those experimentally measurable NQR parameters, quadrupole coupling constants and asymmetry parameters. The results reveal that because of strong contribution to N-H...O=C HBs, NQR parameters of O2, N1 and N3 undergo significant changes from monomer to the target molecule in cluster. Furthermore, the NQR parameters of O2 also undergo some changes because of non-classical C-H...O=C HBs.

Role of Electron-Driven Proton-Transfer Processes in the Excited-State Deactivation of the Adenine−Thymine Base Pair

Exploratory electronic structure calculations have been performed with the CC2 (simplified singles and doubles coupled-cluster) method for two conformers of the adenine (A)-thymine (T) base pair, with emphasis on excited-state proton-transfer reactions. The Watson-Crick conformer and the most stable (in the gas-phase) conformer of the A-T base pair have been considered. The equilibrium geometries of the ground state and of the lowest excited electronic states have been determined with the MP2 (second-order Møller-Plesset) and CC2 methods, respectively. Vertical and adiabatic excitation energies, oscillator strengths, and dipole moments of the excited states are reported. Of particular relevance for the photochemistry of the A-T base pair are optically dark (1)pipi* states of charge-transfer character. Although rather high in energy at the ground-state equilibrium geometry, these states are substantially lowered in energy by the transfer of a proton, which thus neutralizes the charge separation. A remarkable difference of the energetics of the proton-transfer reaction is predicted for the two tautomers of A-T: in the Watson-Crick conformer, but not in the most stable conformer, a sequence of conical intersections connects the UV-absorbing (1)pipi* state in a barrierless manner with the electronic ground state. These conical intersections allow a very fast deactivation of the potentially reactive excited states in the Watson-Crick conformer. The results provide evidence that the specific hydrogen-bonding pattern of the Watson-Crick conformer endows this structure with a greatly enhanced photostability. This property of the Watson-Crick conformer of A-T may have been essential for the selection of this species as carrier of genetic information in early stages of the biological evolution.