Research progress on the role and mechanism of DNA damage repair in germ cell development

In the complex and dynamic processes of replication, transcription, and translation of DNA molecules, a large number of replication errors or damage can occur which lead to obstacles in the development process of germ cells and result in a decreased reproductive rate. DNA damage repair has attracted widespread attention due to its important role in the maintenance and regulation of germ cells. This study reports on a systematic review of the role and mechanism of DNA damage repair in germline development. First, the causes, detection methods, and repair methods of DNA damage, and the mechanism of DNA damage repair are summarized. Second, a summary of the causes of abnormal DNA damage repair in germ cells is introduced along with common examples, and the relevant effects of germ cell damage. Third, we introduce the application of drugs related to DNA damage repair in the treatment of reproductive diseases and related surgical treatment of abnormal DNA damage, and summarize various applications of DNA damage repair in germ cells. Finally, a summary and discussion is given of the current deficiencies in DNA damage repair during germ cell development and future research development. The purpose of this paper is to provide researchers engaged in relevant fields with a further systematic understanding of the relevant applications of DNA damage repair in germ cells and to gain inspiration from it to provide new research ideas for related fields.

How proton transfer affects the helical parameters in DNA:DNA microhelices

Proton transfer reactions are a widespread phenomenon in many areas of the life sciences and it is one of the origins of the spontaneous point mutations during DNA replication. Because of its importance, many studies have been reported on these reactions. However, the present work is the first one focused on the structural geometrical changes by double proton transfer (DPT). Thus, different Watson-Crick (WC) pairs were optimized first in a simple model with one nucleoside base pair, and in a microhelix form with three nucleoside base pairs. The canonical and few tautomeric forms were considered in DNA:DNA microhelices with A-type and B-type helical forms. The stability of these structures and how the DPT process affects the main geometrical parameters was analyzed, in particular the deformation of the helical parameters. The M06-2X DFT method was used for this purpose. The purine/pyrimidine ring in the keto form appears easier to be deformed than when it is in the enol form. The weaker WC base pair formed with mixed microhelices than with nucleobases alone and the significant deformation of the helical and backbone parameters with the DPT appears to complicate this process in microhelices.Communicated by Ramaswamy H. Sarma.

Reduction of Np(VI) with hydrazinopropionitrile via water-mediated proton transfer

Effectively adjusting and controlling the valence state of neptunium (Np) is essential in its separation during spent fuel reprocessing. Hydrazine and its derivatives as free-salts can selectively reduce Np(VI) to Np(V). Reduction mechanisms of Np(VI) with hydrazine and four derivatives have been explored using multiple theoretical methods in our previous works. Herein, we examine the reduction mechanism of Np(VI) with hydrazinopropionitrile (NCCH₂N₂H₃) which exhibits faster kinetics than most other hydrazine derivatives probably due to its σ-π hyperconjugation effect. Free radical ion pathways I, II and III involving the three types of hydrazine H atoms were found that correspond to the experimentally established mechanism of reduction of two Np(VI) via initial oxidation to [NCCH₂N₂H₃]⁺˙, followed by conversion to NCCH₂N₂H (+2H₃O⁺) and ultimately to CH₃CN + N₂. Potential energy profiles suggest that the second redox stage is rate-determining for all three pathways. Pathway I with water-mediated proton transfer is energetically preferred for hydrazinopropionitrile. Analyses using the approaches of localized molecular orbitals (LMOs), quantum theory of atoms in molecules (QTAIM), and intrinsic reaction coordinate (IRC) elucidate the bonding evolution for the structures on the reaction pathways. The results of the spin density reveal that the reduction of the first Np(VI) ion is the outer-sphere electron transfer, while that of the second Np(VI) ion is the hydrogen transfer. This work offers new insights into the nature of reduction of Np(VI) by hydrazinopropionitrile via water-mediated proton transfer, and provides a basis for designing free-salt reductants for Np separations.

Long living excited state of protonated adenosine unveiled by ultrafast fluorescence spectroscopy and density functional theoretical study

Adenosine (Ado) possesses ultrafast nonradiative dynamics accounting for its remarkably high photostability. The deactivation dynamics of Ado after protonation in an aqueous solution remains an elusive issue. Herein we report an investigation of the excited state dynamics of protonated Ado (AdoH⁺) performed using ultrafast time-resolved fluorescence spectroscopy combined with density functional theoretical calculation. The result obtained from comparison of conformers with protonation at different sites revealed that the syn-conformer with protonation occurring at the N3 position (syn-N3) is the predominant form of AdoH⁺ in the ground state, similar to that of Ado. In contrast, the fluorescence of AdoH⁺ with maximum intensity at 385 nm, significantly red-shifted from that of Ado, displaying decay dynamics composed of an ultrafast component with the lifetime of ∼0.5 ps and a slower one of ∼2.9 ns. The former is because of the decay of the syn-N3 conformer, similar to that reported for AdoH⁺ under the gas phase condition. The latter is due to the syn-N1 conformer formed via ultrafast proton transfer of the syn-N3. The excited state of syn-N1 has a peculiar nonplanar conformation over the purine molecule, which is responsible for the substantial Stokes shift showed in the fluorescence spectrum and correlates with a large energy barrier for nonradiative decay likely involving a reversed proton transfer. This study demonstrates the importance of protonation and solvent environment in altering dramatically the excited states of Ado, providing insight for better understanding nonradiative dynamics of both the monomeric bases and the oligomeric or polymeric DNAs.

Insights into plasmon induced keto-enol isomerization

Chemical reactions that are driven by plasmon-induced hot carriers are a timely topic of interest to chemists and material scientists as they provide catalytic alternatives that may reduce cost and/or waste. Herein, we monitored the localized surface plasmon resonance-induced keto-enol isomerization process of 2-mercapto-4(3H)-quinazolinone (MQ) by time-dependent surface enhanced Raman scattering (SERS), where the MQ molecules are adsorbed on gold nanoparticles (GNP) surface by Au-S bonds. The mechanism of keto-enol isomerization has been successfully investigated, and it is found that the isomerization is induced by hot hole transfer from GNPs to the adsorbed molecules. The present investigation could provide significant insights into hot hole catalyzed chemical reactions via SERS spectra and theoretical calculations.

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.

Comparing the performances of various density functionals for modelling the mechanisms and kinetics of bimolecular free radical reactions in aqueous solution

We systematically tested the performances of 18 density functionals for the mechanisms and kinetics of reactions of the α-hydroxyisopropyl radical with 9 organic substrates.

Ultrafast Excited-State Dynamics of Hydrogen-Bonded Cytosine Microsolvated Clusters with Protic and Aprotic Polar Solvents

Microsolvation effects on the ultrafast excited-state deactivation dynamics of cytosine (Cy) were studied in hydrogen-bonded Cy clusters with protic and aprotic solvents using mass-resolved femtosecond pump-probe ionization spectroscopy. Two protic solvents, water (H₂O) and methanol (MeOH), and one aprotic solvent, tetrahydrofuran (THF), were investigated, and transients of Cy·(H₂O)_1-6, Cy·(MeOH)_1-3, and Cy·THF microsolvated clusters produced in supersonic expansions were measured. With the aid of electronic structure calculations, we assigned the observed dynamics to the low-energy isomers of various Cy clusters and discussed the microsolvation effect on the excited-state deactivation dynamics. With the protic solvents only the microsolvated clusters of Cy keto tautomer were observed. The observed decay time constants of Cy·(H₂O) _n are 0.5 ps for n = 1 and ∼0.2-0.25 ps for n = 2-6. For Cy·(MeOH) _n clusters, the decay time constant for n = 1 cluster is similar to that of the Cy monohydrate, but for n = 2 and 3 the decays are about a factor of 2 slower than the corresponding microhydrates. With the aprotic solvent, THF, hydrogen-bonded complexes of both keto and enol tautomers are present in the beam. The keto-Cy·THF shows a decay similar to that of the keto-Cy monomer, whereas the enol-Cy·THF exhibits a 2-fold slower decay than the enol-Cy monomer, suggesting an increase in the barrier to excited-state deactivation upon binding of one THF molecule to the enol form of Cy.

Catalytic Effect of Aqueous Solution in Water-Assisted Proton-Transfer Mechanism of 8-Hydroxy Guanine Radical

Water-assisted proton-transfer process is a key step in guanine damage reaction by hydroxyl radical in aqueous solution. In this article, we quantitatively determine the solvent effect in water-assisted proton-transfer mechanism of 8-hydroxy guanine radical using combined quantum mechanics and molecular mechanism with an explicit solvation model. Atomic-level reaction pathway was mapped, which shows a synchronized two-proton-transfer mechanism between the assistant water molecule and 8-hydroxy guanine radical. The transition-state dipole moment is the largest along the reaction pathway, which electrostatically stabilizes the proton-transfer transition-state complex. The free-energy reaction barrier for this water-assisted proton-transfer reaction was calculated at 19.2 kcal/mol with the density functional theory/M08-SO/cc-pVTZ+/molecular mechanics level of theory. The solvent effect not only has a big impact on geometries, but also dramatically changes the energetics along the reaction pathway. Among the solvent effect contributions to the transition state, the solvent energy contribution is -28.5 kcal/mol and the polarization effect contribution is 19.9 kcal/mol. In total, the solvent effect contributes -8.6 kcal/mol to the free-energy barrier height, which means that the presence of aqueous solution has a catalytic effect on the reaction mechanism and enhances the proton-transfer reactivity in aqueous solution.

Stacking of the mutagenic base analogue 5-bromouracil: energy landscapes of pyrimidine dimers in gas phase and water

The potential energy surfaces of stacked base pairs consisting of cytosine (C), thymine (T), uracil (U) and the mutagenic thymine analogue 5-bromouracil (BrU) have been searched to obtain all possible minima. Minima and transition states were optimised at the counterpoise-corrected M06-2X/6-31+G(d) level, both in the gas phase and in water, modelled by the polarizable continuum model. The stacked dimers studied are BrU/BrU, C/BrU, C/C, C/T, C/U, T/BrU and T/U. Both face-to-back and face-to-face structures were considered. Free energies were calculated at 298.15 K. Together with U/U, T/T and BrU/U results from previous work, these results complete the family consisting of every stacked dimer combination consisting of C, T, U and BrU. The results were used to assess the hypothesis suggested in the literature that BrU stacks stronger than T, which could stabilise the mispair formed by BrU and guanine. In the gas phase, structures of C/BrU, T/BrU and U/BrU with greater zero-point-corrected binding energies than C/T, T/T and U/T, respectively, were found, with differences in favour of BrU of 3.1 kcal mol(-1), 1.7 kcal mol(-1) and 0.5 kcal mol(-1), respectively. However, the structure of these dimers differed considerably from anything encountered in DNA. When only the dimers with the most "DNA-like" twist (±36°) were considered, C/BrU and T/BrU were still more strongly bound than C/T and T/T, by 0.5 kcal mol(-1) and 1.7 kcal mol(-1), respectively. However, when enthalpic and/or solvent contributions were taken into account, the stacking advantage of BrU was reversed in the gas phase and mostly nullified in water. Enhanced stacking therefore does not seem a plausible mechanism for the considerably greater ability of BrU-G mispairs over T-G mispairs to escape enzymatic repair.

A novel conception for spontaneous transversions caused by homo-pyrimidine DNA mismatches: a QM/QTAIM highlight

We have firstly shown that the T·T(w) and C·C(w) DNA mismatches with wobble (w) geometry stay in slow tautomeric equilibrium with short T·T*(WC) and C·C*(WC) Watson-Crick (WC) mispairs. These non-dissociative tautomeric rearrangements are controlled by the plane-symmetric, highly stable, highly polar and zwitterionic transition states. The obtained results allow us to understand in what way the T·T(w) and C·C(w) mismatches acquire enzymatically competent T·T*(WC) and C·C*(WC) conformations directly in the hydrophobic recognition pocket of a high-fidelity DNA-polymerase, thereby producing thermodynamically non-equilibrium spontaneous transversions. The simplest numerical estimation of the frequency ratio of the TT to CC spontaneous transversions satisfactorily agrees with experimental data.

Tautomeric transition between wobble A·C DNA base mispair and Watson-Crick-like A·C* mismatch: microstructural mechanism and biological significance

Here, we use MP2/DFT quantum-chemical methods combined with Quantum Theory of Atoms in Molecules to study the tautomeric transition between wobble A·C(w) mismatch and Watson-Crick-like A·C*(WC) base mispair, proceeding non-dissociatively via sequential proton transfer between bases through the planar, highly stable and zwitterionic TS(A∙C-)(A∙C(W)<-->A∙C&(WC)) transition state joined by the participation of (A)N6(+)H∙∙∙N4(-)(C), (A)N1(+)H∙∙∙N4(-)(C) and (A)C2(+)H∙∙∙N3(-)(C) H-bonds. Notably, the A·C(w) ↔ A·C*(WC) tautomerization reaction is accompanied by 10 unique patterns of the specific intermolecular interactions that consistently replace each other. Our data suggest that biologically significant A·C(w) → A·C*(WC) tautomerization is a kinetically controlled pathway for formation of the enzymatically competent Watson-Crick-like A·C*(WC) DNA base mispair in the essentially hydrophobic recognition pocket of the high-fidelity DNA-polymerase, responsible for the occurrence of spontaneous point AC/CA incorporation errors during DNA biosynthesis.

Theoretical study of enzymatically catalyzed tautomerization of carbon acids in aqueous solution: quantum calculations and steered molecular dynamics simulations

The thermodynamics and kinetics of enzymatically assisted reactions of carbon acids were studied theoretically in this work. Quantum electronic (QE) structure calculations and steered molecular dynamics (SMD) simulations were carried out. Three 3-butenal tautomerization reactions that proceed from the β,γ-unsaturated reactant (R) to the α,β-unsaturated carbon acid product (P) and occur in two elementary steps through an intermediate (I) were studied, ignoring or including the surrounding aqueous medium in the calculations. The Gibbs free energies of activation of the R ⇆ I enolization and I ⇆ P ketonization steps were found to decrease considerably when residues simulating enzymes were introduced into these processes. Although the processes became slightly more favorable thermodynamically when the solution was included in the simulations, they became less favorable kinetically. The results from SMD simulations of these reactions were qualitatively consistent with the values we obtained using QE as well as those found by other authors in similar studies. Our simulations also allowed us to perform a detailed study of these reactions in solution.

In silico studies with substituted adenines to achieve a remarkable stability of mispairs with thymine nucleobase

The modified adenine and thymine mispair achieves a remarkable stability, which can presumably help the DNA lesions to be less cytotoxic.

Whether the amino–imino tautomerism of 2-aminopurine is involved into its mutagenicity? Results of a thorough QM investigation

2AP* mutagenic tautomer is able to induce only one incorporation error – transversion – by pairing through the H-bonds into the G·2AP* mispair.

Wobble↔Watson-Crick tautomeric transitions in the homo-purine DNA mismatches: a key to the intimate mechanisms of the spontaneous transversions

The intrinsic capability of the homo-purine DNA base mispairs to perform wobble↔Watson-Crick/Topal-Fresco tautomeric transitions via the sequential intrapair double proton transfer was discovered for the first time using QM (MP2/DFT) and QTAIM methodologies that are crucial for understanding the microstructural mechanisms of the spontaneous transversions.

How many tautomerization pathways connect Watson-Crick-like G*·T DNA base mispair and wobble mismatches?

In this study, we have theoretically demonstrated the intrinsic ability of the wobble G·T(w)/G*·T*(w)/G·T(w1)/G·T(w2) and Watson-Crick-like G*·T(WC) DNA base mispairs to interconvert into each other via the DPT tautomerization. We have established that among all these transitions, only one single G·T(w) ↔ G*·T(WC) pathway is eligible from a biological perspective. It involves short-lived intermediate - the G·T*(WC) base mispair - and is governed by the planar, highly stable, and zwitterionic [Formula: see text] transition state stabilized by the participation of the unique pattern of the five intermolecular O6(+)H⋯O4(-), O6(+)H⋯N3(-), N1(+)H⋯N3(-), N1(+)H⋯O2(-), and N2(+)H⋯O2(-) H-bonds. This non-dissociative G·T(w) ↔ G*·T(WC) tautomerization occurs without opening of the pair: Bases within mispair remain connected by 14 different patterns of the specific intermolecular interactions that successively change each other along the IRC. Novel kinetically controlled mechanism of the thermodynamically non-equilibrium spontaneous point GT/TG incorporation errors has been suggested. The mutagenic effect of the analogues of the nucleotide bases, in particular 5-bromouracil, can be attributed to the decreasing of the barrier of the acquisition by the wobble pair containing these compounds of the enzymatically competent Watson-Crick's geometry via the intrapair mutagenic tautomerization directly in the essentially hydrophobic recognition pocket of the replication DNA-polymerase machinery. Proposed approaches are able to explain experimental data, namely growth of the rate of the spontaneous point incorporation errors during DNA biosynthesis with increasing temperature.

How does the long G·G* Watson-Crick DNA base mispair comprising keto and enol tautomers of the guanine tautomerise? The results of a QM/QTAIM investigation

The double proton transfer (DPT) in the long G·G* Watson-Crick base mispair (|C6N1(G*)N1C6(G)| = 36.4°; C1 symmetry), involving keto and enol tautomers of the guanine (G) nucleobase, along two intermolecular neighboring O6H···O6 (8.39) and N1···HN1 (6.14 kcal mol(-1)) H-bonds that were established to be slightly anti-cooperative, leads to its transformation into the G*·G base mispair through a single transition state (|C6N1N1C6| = 37.1°; C1), namely to the interconversion into itself. It was shown that the G·G* ↔ G*·G tautomerisation via the DPT is assisted by the third specific contact, that sequentially switches along the intrinsic reaction coordinate (IRC) in an original way: (G)N2H···N2(G*) H-bond (-25.13 to -10.37) → N2···N2 van der Waals contact (-10.37 to -9.23) → (G)N2···HN2(G*) H-bond (-9.23 to 0.79) → (G*)N2···HN2(G) H-bond (0.79 to 7.35 Bohr). The DPT tautomerisation was found to proceed through the asynchronous concerted mechanism by employing the QM/QTAIM approach and the methodology of the scans of the geometric, electron-topological, energetic, polar and NBO properties along the IRC. Nine key points, that can be considered as part of the tautomerisation repertoire, have been established and analyzed in detail. Furthermore, it was shown that the G·G* or G*·G base mispair is a thermodynamically and dynamically stable structure with a lifetime of 8.22 × 10(-10) s and all 6 low-frequency intermolecular vibrations are able to develop during this time span. Lastly, our results highlight the importance of the G·G* ↔ G*·G DPT tautomerisation, which can have implications for biological and chemical sensing applications.

DPT tautomerisation of the G·A(syn) and A*·G*(syn) DNA mismatches: a QM/QTAIM combined atomistic investigation

By applying a combined QM and QTAIM atomistic computational approach we have established for the first time that the G·A(syn) and A*·G*(syn) DNA mismatches (rare tautomers are marked with an asterisk), causing spontaneous transversions with substantially various probabilities, radically differ from each other in their ability to tautomerise through the double proton transfer (DPT). The A*·G*(syn) mismatch tautomerises quite easily (ΔΔG(TS) ≈ 4·kT at room temperature) into the A·G*(syn) mismatch through the asynchronous concerted mechanism, whereas the G·A(syn) base mispair does not tautomerise via the DPT at all, since there is no local minimum corresponding to the tautomerised G*·A*(syn) mismatch on the potential energy surface. It was established that the A·G*(syn) base mispair is a dynamically unstable H-bonded complex with an extremely short lifetime of 2.17 × 10(-13) s. Consequently, the obtained results allow us to believe that spontaneous or forced dissociation of both the G·A(syn) and A*·G*(syn) DNA mismatches by the DNA-polymerase occurs with the preservation of the tautomeric status of the bases.

Direct and solvent-assisted keto-enol tautomerism and hydrogen-bonding interactions in 4-(m-chlorobenzylamino)-3-phenyl-4,5-dihydro-1H-1,2,4-triazol-5-one: a quantum-chemical study

The tautomeric equilibrium of the title triazole compound was computationally analyzed at the B3LYP/6-311++G(d,p) and MP2/6-311++G(d,p) levels of theory. The solvent effect was considered for three solvents (chloroform, methanol, and water). Two distinct mechanisms were applied: a direct intramolecular transfer using the polarizable continuum model (PCM) and a solvent-assisted mechanism. The calculations indicated that the keto form is more stable in all cases. It was found that the barrier heights for the tautomerization reaction are very high, indicating a relatively disfavored process. Although the barrier heights for solvent-assisted reactions are significantly lower than those for the unassisted tautomerization reaction, implying the importance of the superior catalytic effect of the solvents, monosolvation was not found to be sufficient for the reaction to occur. Finally, the two intermolecular hydrogen-bonding interactions in the crystal structure were investigated in the gas phase; according to the calculated energies and structural parameters, the order of stability is N3-H3···O1 > N1-H1···O1.

Novel physico-chemical mechanism of the mutagenic tautomerisation of the Watson–Crick-like A·G and C·T DNA base mispairs: a quantum-chemical picture

Novel routes for the mutagenic tautomerisation of the long A·G and short C·T Watson–Crick DNA base mispairs via sequential DPT are reported, pursuing the goal of an estimation of their contribution into spontaneous point replication errors in DNA.

The nature of the transition mismatches with Watson-Crick architecture: the G*·T or G·T* DNA base mispair or both? A QM/QTAIM perspective for the biological problem

This study provides the first accurate investigation of the tautomerization of the biologically important guanine*·thymine (G*·T) DNA base mispair with Watson-Crick geometry, involving the enol mutagenic tautomer of the G and the keto tautomer of the T, into the G·T* mispair (∆G = .99 kcal mol(-1), population = 15.8% obtained at the MP2 level of quantum-mechanical theory in the continuum with ε = 4), formed by the keto tautomer of the G and the enol mutagenic tautomer of the T base, using DFT and MP2 methods in vacuum and in the weakly polar medium (ε = 4), characteristic for the hydrophobic interfaces of specific protein-nucleic acid interactions. We were first able to show that the G*·T↔G·T* tautomerization occurs through the asynchronous concerted double proton transfer along two antiparallel O6H···O4 and N1···HN3 H-bonds and is assisted by the third N2H···O2 H-bond, that exists along the entire reaction pathway. The obtained results indicate that the G·T* base mispair is stable from the thermodynamic point of view complex, while it is dynamically unstable structure in vacuum and dynamically stable structure in the continuum with ε = 4 with lifetime of 6.4·10(-12) s, that, on the one side, makes it possible to develop all six low-frequency intermolecular vibrations, but, on the other side, it is by three orders less than the time (several ns) required for the replication machinery to forcibly dissociate a base pair into the monomers during DNA replication. One of the more significant findings to emerge from this study is that the short-lived G·T* base mispair, which electronic interaction energy between the bases (-23.76 kcal mol(-1)) exceeds the analogical value for the G·C Watson-Crick nucleobase pair (-20.38 kcal mol(-1)), "escapes from the hands" of the DNA replication machinery by fast transforming into the G*·T mismatch playing an indirect role of its supplier during the DNA replication. So, exactly the G*·T mismatch was established to play the crucial role in the spontaneous point mutagenesis.

Atomistic nature of the DPT tautomerisation of the biologically important C·C* DNA base mispair containing amino and imino tautomers of cytosine: a QM and QTAIM approach

A theoretical study of tautomerisation of the biologically important cytosine·cytosine* (C·C*) DNA mismatch with a propeller-like structure (|C4N3N3C4| = 32.4°; C1 symmetry) and cis-oriented N1H glycosidic bonds, formed by the amino and imino tautomers of the C nucleobase, via the asynchronous concerted double proton transfer (DPT) along two H-bonds through the transition state (TSC·C*↔C*·C) (|C4N3N3C4| = 48.5°; C1 symmetry) into the C*·C mispair was carried out for the first time. It was established that the C·C*/C*·C DNA base mispair is associated by the antiparallel N4H···N4 (6.66 kcal mol(-1)), N3H···N3 (6.47 kcal mol(-1)) H-bonds and the O2···O2 van der Waals (vdW) contact (0.33 kcal mol(-1)), while the zwitterionic TSC·C*↔C*·C is stabilized by the parallel N4(+)H···N4(-) (13.55 kcal mol(-1)), N3(+)H···N3(-) (13.20 kcal mol(-1)) H-bonds and the O2(+)···O2(-) vdW contact (0.60 kcal mol(-1)). It was shown that the C·C* ↔ C*·C tautomerisation via the DPT is assisted by the O2···O2 vdW contact, that in contrast to the two others N4H···N4 and N3H···N3 H-bonds exists along the entire intrinsic reaction coordinate (IRC) range. The positive values of the Grunenberg's compliance constants (30.919 and 21.384 Å mdyn(-1) for C·C*/C*·C and TSC·C*↔C*·C, respectively) indicate that the O2···O2 vdW contact is a stabilizing closed-shell interaction. It was found that the middle N3H···N3 H-bond is anti-cooperative with the upper N4H···N4 H-bond and cooperative with the lower O2···O2 vdW contact. The 9 key points, which can be considered as electron-topological "fingerprints" of the asynchronous concerted C·C* ↔ C*·C tautomerisation process via the DPT were revealed along the IRC and examined in detail. It was shown that the C·C*/C*·C base mispair is a thermodynamically and dynamically stable structure. Its lifetime is equal to 1.53 × 10(-7) s at the MP2/cc-pVQZ//B3LYP/6-311++G(d,p) level of theory in vacuum. All 6 low-frequency intermolecular vibrations are able to develop during this time span.

B-DNA characteristics are preserved in double stranded d(A)3·d(T)3 and d(G)3·d(C)3 mini-helixes: conclusions from DFT/M06-2X study

We report the results of the first comprehensive DFT study on the d(A)3·d(T)3 and d(G)3·d(C)3 nucleic acid duplexes. The ability of mini-helixes to preserve the conformation of B-DNA in the gas phase and under the influence of such factors as: solvent, uncompensated charge, and counter-ions was evaluated using M06-2X functional with 6-31G(d,p) basis set. The accuracy of the models was ascertained based on their ability to reproduce key structural features of natural B-DNA. Analysis of the helicity suggests that the helical conformations adopt geometrical parameters which are close to those of the B-DNA form. The torsion angles fall somewhere between the values observed for BI/BII conformational classes. The comparative analysis of parameters of isolated Watson-Crick base pairs versus B-DNA-like conformations indicates the same tendency of base-pair polarization and hydration. Specifically, effects of polarization of nucleobases in continuum type dielectric medium mimicking water are stronger than those caused by the presence of backbone. Polar environment as well as the presence of counterions stabilizes duplexes, facilitating helix formation. Substantial conformational changes of nucleotides upon duplex formation decrease the binding energy. In spite of structural and energetic changes, the placement of a mini-helix into the gas phase does not lead to significant disruption of the structure. On the contrary, the duplex preserves its helicity and the strands remain bound.

DPT tautomerisation of the wobble guanine·thymine DNA base mispair is not mutagenic: QM and QTAIM arguments

We have shown for the first time, connecting QM methods with QTAIM analysis and using the methodology of the sweeps of the energetical, electron-topological and geometrical parameters, that the tautomerisation of the wobble guanine·thymine (wG·T) DNA base mispair into the wG(*)·T(*) base mispair induced by the double proton transfer (DPT), which undergoes a concerted asynchronous pathway, is not mutagenic. The wG·T → wG(*)·T(*) DPT tautomerisation does not result in the transition of the G base into its mutagenic tautomeric form G(*) able to mispair with the T base within the Watson-Crick base pairing scheme. This observation is explained by the so-called quantum protection of the wG·T DNA base mispair from its mutagenic tautomerisation - the dynamical non-stability of the tautomerised wG(*)·T(*) base mispair and significantly negative value of the Gibbs free energy of activation for the reverse reaction of the wG·T → wG(*)·T(*) DPT tautomerisation.

Is the DPT tautomerization of the long A·G Watson-Crick DNA base mispair a source of the adenine and guanine mutagenic tautomers? A QM and QTAIM response to the biologically important question

Herein, we first address the question posed in the title by establishing the tautomerization trajectory via the double proton transfer of the adenine·guanine (A·G) DNA base mispair formed by the canonical tautomers of the A and G bases into the A*·G* DNA base mispair, involving mutagenic tautomers, with the use of the quantum-mechanical calculations and quantum theory of atoms in molecules (QTAIM). It was detected that the A·G ↔ A*·G* tautomerization proceeds through the asynchronous concerted mechanism. It was revealed that the A·G base mispair is stabilized by the N6H···O6 (5.68) and N1H···N1 (6.51) hydrogen bonds (H-bonds) and the N2H···HC2 dihydrogen bond (DH-bond) (0.68 kcal·mol(-1) ), whereas the A*·G* base mispair-by the O6H···N6 (10.88), N1H···N1 (7.01) and C2H···N2 H-bonds (0.42 kcal·mol(-1) ). The N2H···HC2 DH-bond smoothly and without bifurcation transforms into the C2H···N2 H-bond at the IRC = -10.07 Bohr in the course of the A·G ↔ A*·G* tautomerization. Using the sweeps of the energies of the intermolecular H-bonds, it was observed that the N6H···O6 H-bond is anticooperative to the two others-N1H···N1 and N2H···HC2 in the A·G base mispair, while the latters are significantly cooperative, mutually strengthening each other. In opposite, all three O6H···N6, N1H···N1, and C2H···N2 H-bonds are cooperative in the A*·G* base mispair. All in all, we established the dynamical instability of the А*·G* base mispair with a short lifetime (4.83·10(-14) s), enabling it not to be deemed feasible source of the A* and G* mutagenic tautomers of the DNA bases. The small lifetime of the А*·G* base mispair is predetermined by the negative value of the Gibbs free energy for the A*·G* → A·G transition. Moreover, all of the six low-frequency intermolecular vibrations cannot develop during this lifetime that additionally confirms the aforementioned results. Thus, the A*·G* base mispair cannot be considered as a source of the mutagenic tautomers of the DNA bases, as the A·G base mispair dissociates during DNA replication exceptionally into the A and G monomers in the canonical tautomeric form.