Atomistic mechanisms of the tautomerization of the G·C base pairs through the proton transfer: quantum-chemical survey.

This study is devoted to the investigation of the G·C*tO2(WC)↔G*NH3·C*t(WC), G·C*O2(WC)↔G*NH3·C*(WC) and G*·C*O2(WC)↔G*NH3·C(wWC)↓ tautomerization reactions occurring through the proton transfer, obtained at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of theory in gas phase under normal conditions ('WC' means base pair in Watson-Crick configuration, T=298.15 K). These reactions lead to the formation of the G*NH3·C*t(WC), G*NH3·C*(WC) and G*NH3·C(wWC)↓ base pairs by the participation of the G*NH3 base with NH3 group. Gibbs free energies of activation for these reactions are 6.43, 11.00 and 1.63 kcal·mol-1, respectively. All of these tautomerization reactions are dipole active. Finally, we believe that these non-dissociative processes, which are tightly connected with the tautomeric transformations of the G·C base pairs, play an outstanding role in supporting of the spatial structure of the DNA and RNA molecules with various functional purposes.

Novel mechanisms of the conformational transformations of the biologically important G·C nucleobase pairs in Watson-Crick, Hoogsteen and wobble configurations via the mutual rotations of the bases around the intermolecular H-bonds: a QM/QTAIM study.

At the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of quantum-mechanical theory, we provide for the first time a comprehensive investigation of the physico-chemical mechanisms of the 55 conformational transformations of the biologically-important G·C nucleobase pairs - Watson-Crick (WC), reverse Watson-Crick (rWC), Hoogsteen (H), reverse Hoogsteen (rH), wobble (w) and reverse wobble (rw) base pairs by the participation of the G and C bases in the canonical and rare tautomeric forms ("r" - means reverse configuration of the base pair). It was established that all these G·C nucleobase pairs can conformationally transform into each other without the changing of the tautomeric status of the G and C bases. These transitions occur through significantly non-planar transition states via the mutual rotation of the G and C bases relative to each other within the G·C nucleobase pair around the upper, middle or lower intermolecular H-bonds: WC ↔ rWC, WC ↔ rwWC, rWC ↔ WC, rWC ↔ wWC, wWC ↔ rwWC, H ↔ rH, H ↔ rwH, rH ↔ H, rH ↔ wH, wH ↔ rwH. Gibbs free energies ΔG of activation for these conformational transformations are ΔG = 2.96-19.04/3.58-13.36 kcal mol-1 (in vacuum under normal conditions (T = 298.15 K)), which means that these reactions proceed quite fast. Obtained conformational transformations are accompanied by the disruption and further formation of the intermolecular specific contacts in the G·C nucleobase pairs (H-bonds and attractive van der Waals contacts). As a result, 76 conformers of the G·C nucleobase pairs were established - 48 base pairs in WC, rWC, wWC and rwWC configurations and 28 base pairs in H, rH, wH and rwH configurations with relative Gibbs free ΔG/electronic ΔE energies in the range ΔG/ΔE = 0.00-44.73/0.00-46.99 and ΔG/ΔE = 0.00-37.52/0.00-38.54 kcal mol-1, respectively (in vacuum under normal conditions). Experimental investigation and verification of the novel G·C nucleobase pairs are promising tasks for the future research. Based on the obtained data, biologically important conclusions were made about the importance of the conformational mobility of the G·C nucleobase pairs for the understanding of the functioning of the DNA and RNA molecules and their transition from the parallel into the anti-parallel duplexes and vice versa.

A Quantum-Mechanical Looking Behind the Scene of the Classic G·C Nucleobase Pairs Tautomerization.

For the first time, at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of theory, a comprehensive quantum-mechanical investigation of the physico-chemical mechanism of the tautomeric wobblization of the four biologically-important G·C nucleobase pairs by the participation of the monomers in rare, in particular mutagenic, tautomeric forms (marked with an asterisk) was provided. These novel tautomeric transformations (wobblization or shifting of the bases within the pair) are intrinsically inherent properties of the G·C nucleobase pairs. In this study, we have obtained intriguing results, lying far beyond the existing representations. Thus, it was shown that Löwdin's G*·C*(WC) base pair does not tautomerize according to the wobblization mechanism. Tautomeric wobblization of the G*·C*(rWC) (relative Gibbs free energy ΔG = 0.00/relative electronic energy ΔE = 0.00 kcal·mol-1) ("r"-means the configuration of the base pair in reverse position; "WC"-the classic Watson-Crick configuration) and G*t·C*(H) (ΔG = -0.19/ΔE = 0.29 kcal·mol-1) ("H"-Hoogsteen configuration;"t" denotes the O6H hydroxyl group in the trans position) base pairs are preceded by the stages of the base pairs tautomerization by the single proton transfer (SPT). It was established that the G*t·C*(rH) (ΔG = 2.21/ΔE = 2.81 kcal·mol-1) base pair can be wobbled through two different pathways via the traditional one-stage mechanism through the TSs, which are tight G+·C- ion pairs, stabilized by the participation of only two intermolecular H-bonds. It was found out that the G·C base pair is most likely incorporated into the DNA/RNA double helix with parallel strands in the G*·C*(rWC), G·C*(rwwc), and G*·C(rwwc) ("w"-wobble configuration of the pair) tautomeric forms, which are in rapid tautomeric equilibrium with each other. It was proven that the G*·C*(rWC) nucleobase pair is also in rapid tautomeric equilibrium with the eight tautomeric forms of the so-called Levitt base pair. It was revealed that a few cases of tautomerization via the DPT of the nucleobase pairs by the participation of the C8H group of the guanine had occurred. The biological role of the obtained results was also made apparent.

Conformational transitions of the quercetin molecule via the rotations of its rings: a comprehensive theoretical study.

Quercetin is an important flavonoid compound, usually extracted from plants, vegetables and fruits such as blueberries, apples, green tea, wine, onions and possessing broad range of pharmacological properties, in particular, powerful antioxidant, antitoxic, antiinflammation and antimicrobial effects due to its various reactive sites. The structure of this phenolic compound consists of three (A + C) and B rings, bearing five hydroxyl groups. Primarily, the chemical structure of quercetin determines its physico-chemical properties. Earlier, it was established that isolated quercetin molecule can acquire 48 stable conformations (24 planar and 24 non-planar) due to the mobility of its hydroxyl groups and (A + C) and B rings with relative Gibbs free energies in the range of 0.0-25.3 kcal·mol-1 under normal conditions (Brovarets' et al., 2019c). In this work by quantum-mechanical calculations at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of theory and Bader's 'Quantum Theory of Atoms in Molecules', we have theoretically modeled the interconversions in the 24 pairs of the conformers of the quercetin molecule, occuring via the rotation of its non-deformable (A+С) and B rings around the С2-С1' bond through the quasi-orthogonal transition state with low values of the imaginary frequencies (28-33/29-36 cm-1) and Gibbs free energies of activation in the range of 2.17-5.68/1.86-4.90 kcal·mol-1 in the continuum with dielectric permittivity ε = 1/ε = 4 under normal conditions. Also, we studied the changes of the number of physico-chemical characteristics of all intramolecular-specific contacts - hydrogen bonds and attractive van der Waals contacts during these conformational rearrangements.Communicated by Ramaswamy H. Sarma.

Conformational diversity of the quercetin molecule: a quantum-chemical view.

This paper focuses on the comprehensive conformational analysis of the quercetin molecule with a broad range of the therapeutic and biological actions. All possible conformers of these molecule, corresponding to the local minima on the potential energy hypersurface, have been obtained by the sequential rotation of its five hydroxyl groups and also by the rotation of its (A + C) and B rings relatively each other. Altogether, it was established 48 stable conformers, among which 24 conformers possess planar structure and 24 conformers - nonplanar structure. Their structural, symmetrical, energetical and polar characteristics have been investigated in details. Quantum-mechanical calculations indicate that conformers of the quercetin molecule are polar structures with a dipole moment, which varies within the range from 0.35 to 9.87 Debay for different conformers. Relative Gibbs free energies of these conformers are located within the range from 0.0 to 25.3 kcal·mol-1 in vacuum under normal conditions. Impact of the continuum with ε = 4 leads to the decreasing of the Gibbs free energies (-0.19-18.15 kcal·mol-1) and increasing of the dipole moment (0.57-12.48 D). It was shown that conformers of the quercetin molecule differ from each other by the intramolecular specific contacts (two or three), stabilizing all possible conformers of the molecule - H-bonds (both classical ОН…Ðž and so-called unusual С'Н…Ðž and ОН…Ð¡') and attractive van-der-Waals contacts О…Ðž. Obtained conformational analysis for the quercetin molecule enables to provide deeper understanding of the 'structure-function' relationship and also to suggest its mechanisms of the therapeutic and biological actions.Communicated by Ramaswamy H. Sarma.

A new era of the prototropic tautomerism of the quercetin molecule: A QM/QTAIM computational advances.

In this study for the first time we have revealed and investigated in details 123 different prototropic tautomers of the most stable conformer of the quercetin molecule using quantum-mechanical calculations at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of QM theory. We have found that in the most energetically favorable prototropic tautomer mobile hydrogen atoms are localized at the О3, О3', О4', О5, and О7 exocyclic oxygen atoms. Molecular tautomers are in the range of the Gibbs free energies from 0.0 to 69.8 kcal·mol-1, while zwitterionic ones - from 30.1 до 172.8 kcal·mol-1 at normal conditions. It was also reliably established that the weakest point causing the decyclization of the molecule is its C ring - this reaction is launched by the transition of the proton from the C8H group to the endocyclic O1 oxygen atom. All prototropic tautomers, except two cases, are joined by the intramolecular cooperative specific interactions (from 1 to 5) - H-bonds and attractive van der Waals contacts, which have been revealed and characterized by QTAIM analysis. Communicated by Ramaswamy H. Sarma.

Intramolecular tautomerization of the quercetin molecule due to the proton transfer: QM computational study.

Quercetin molecule (3, 3', 4', 5, 7-pentahydroxyflavone, C15H10O7) is an important flavonoid compound of natural origin, consisting of two aromatic A and B rings linked through the C ring with endocyclic oxygen atom and five hydroxyl groups attached to the 3, 3', 4', 5 and 7 positions. This molecule is found in many foods and plants, and is known to have a wide range of therapeutic properties, like an anti-oxidant, anti-toxic, anti-inflammatory etc. In this study for the first time we have revealed and investigated the pathways of the tautomeric transformations for the most stable conformers of the isolated quercetin molecule (Brovarets' & Hovorun, 2019) via the intramolecular proton transfer. Energetic, structural, dynamical and polar characteristics of these transitions, in particular relative Gibbs free and electronic energies, characteristics of the intramolecular specific interactions-H-bonds and attractive van der Waals contacts, have been analysed in details. It was demonstrated that the most probable process among all investigated is the proton transfer from the O3H hydroxyl group of the C ring to the C2' carbon atom of the C2'H group of the B ring along the intramolecular O3H…C2' H-bond with the further formation of the C2'H2 group. It was established that the proton transfer from the hydroxyl groups to the carbon atoms of the neighboring CH groups is assisted at the transition states by the strong intramolecular HCH…O H-bond (~28.5 kcal∙mol-1). The least probable path of the proton transfer-from the C8H group to the endocyclic O1 oxygen atom-causes the decyclization of the C ring in some cases. It is shortly discussed the biological importance of the obtained results.

Novel Tautomerisation Mechanisms of the Biologically Important Conformers of the Reverse Löwdin, Hoogsteen, and Reverse Hoogsteen G*·C* DNA Base Pairs via Proton Transfer: A Quantum-Mechanical Survey.

For the first time, in this study with the use of QM/QTAIM methods we have exhaustively investigated the tautomerization of the biologically-important conformers of the G*·C* DNA base pair-reverse Löwdin G*·C*(rWC), Hoogsteen G*'·C*(H), and reverse Hoogsteen G*'·C*(rH) DNA base pairs-via the single (SPT) or double (DPT) proton transfer along the neighboring intermolecular H-bonds. These tautomeric reactions finally lead to the formation of the novel G· C O 2 * (rWC), G N 2 * · C(rWC), G*'N2·C(rWC), G N 7 * · C(H), and G*'N7·C(rH) DNA base mispairs. Gibbs free energies of activation for these reactions are within the range 3.64-31.65 kcal·mol-1 in vacuum under normal conditions. All TSs are planar structures (Cs symmetry) with a single exception-the essentially non-planar transition state TSG*·C*(rWC)↔G+·C-(rWC) (C1 symmetry). Analysis of the kinetic parameters of the considered tautomerization reactions indicates that in reality only the reverse Hoogsteen G*'·C*(rH) base pair undergoes tautomerization. However, the population of its tautomerised state G*'N7·C(rH) amounts to an insignificant value-2.3·10-17. So, the G*·C*(rWC), G*'·C*(H), and G*'·C*(rH) base pairs possess a permanent tautomeric status, which does not depend on proton mobility along the neighboring H-bonds. The investigated tautomerization processes were analyzed in details by applying the author's unique methodology-sweeps of the main physical and chemical parameters along the intrinsic reaction coordinate (IRC). In general, the obtained data demonstrate the tautomeric mobility and diversity of the G*·C* DNA base pair.

Atomistic mechanisms of the double proton transfer in the H-bonded nucleobase pairs: QM/QTAIM computational lessons.

In this Review, we have summarized and generalized the results of the investigation of the microstructural mechanisms of the tautomerization by the counter movement of the protons along the neighboring intermolecular H-bonds in 22 biologically important pairs of nucleotide bases in the framework of the original method, which allows to trace the evolution of the physicochemical parameters, that characterize these processes along the intrinsic reaction coordinate (IRC). It was demonstrated the performance of the introduction of the conception of the key points (KPs) (from nine to five, depending on the symmetry and nature of system), which exhaustively characterize the flow of the tautomerization processes. It was proved that for all tautomerizing base pairs the extrema of the first derivative of the electron energy of the complex by IRC coincide with the second and penultimate KPs, in which the Laplacian of the electron density equals zero at the corresponding (3,-1) bond critical points of the H-bonds. It was established the linear dependence of the width of the transition state zone of the DPT tautomerization on the degree of its asynchrony. Authors emphasize that the tautomerization reaction through the DPT of the H-bonded pairs of nucleotide bases can be considered successful in those and only in those case if the tautomerized complex is a dynamically stable system, during lifetime of which low-frequency intermolecular vibrations could develop. Perspectives of the application of the obtained approaches to the thorough study of the proton transfer processes in the biologically important objects have been briefly discussed.

Key microstructural mechanisms of the 2-aminopurine mutagenicity: Results of extensive quantum-chemical research.

As of today, a great amount of experimental and theoretical phenomenological data have been collected in the literature according the mutagenic action of the classical mutagen - 2-aminopurine (2AP). However, so far they have not received proper explanation and substantiation. In this Opinion Piece, we provide an overview of recent progress in computational design and modeling of the physico-chemical mechanisms of the mutagenic action of 2AP. Results of quantum-chemical studies, aimed at the elucidation of the key microstructural mechanisms of the mutagenicity of 2AP, have been summarized here. In this context, for the first time it was outlined the most important surveys: Why 2AP is incorporated into DNA in trace concentrations? Whether classical mechanisms presented in the literature according the formation of the rare tautomers of canonical DNA bases work also for base analogue - 2AP? In what way 2AP induces replication and incorporation errors? Whether the amino-imino tautomerisation of 2AP is related to its mutagenicity, that is whether the 2AP* rare tautomer is mutagenic? It is emphasized that the applied approach has a proper theoretical substantiation, since it is based on our microstructural theory of the spontaneous point mutagenesis in DNA, and at the same time it accumulates scenarios of the origin of the induced point errors - transitions and transversions, which the classical Watson-Crick tautomeric hypothesis permits. Moreover, using author's methodology, the profiles of the main physico-chemical characteristics for the tautomerisation reactions involving 2AP, which are integral parts of the biologically important tautomerically-conformational transformations, have been presented. Obtained results open new perspectives for prediction and design of the mutagenic derivatives of the nucleotide bases of any structure and origin before their synthesis and also for planning of new experiments and interpretation of the existing data. Abbreviations 2AP 2-Aminopurine A adenine C cytosine DPT double proton transfer G guanine IRC intrinsic reaction coordinate KP key point T thymine w wobble WC Watson-Crick vdW van der Waals H-bond hydrogen bond Communicated by Ramaswamy H. Sarma.

Unexpected Routes of the Mutagenic Tautomerization of the T Nucleobase in the Classical A·T DNA Base Pairs: A QM/QTAIM Comprehensive View.

In this paper using quantum-mechanical (QM) calculations in combination with Bader's quantum theory of "Atoms in Molecules" (QTAIM) in the continuum with ε = 1, we have theoretically demonstrated for the first time that revealed recently highly-energetic conformers of the classical A·T DNA base pairs - Watson-Crick [A·T(wWC)], reverse Watson-Crick [A·T(wrWC)], Hoogsteen [A·T(wH)] and reverse Hoogsteen [A·T(wrH)] - act as intermediates of the intrapair mutagenic tautomerization of the T nucleobase owing to the novel tautomerisation pathways: A·T(wWC)↔A·T*(w⊥ WC); A·T(wrWC)↔A· T O 2 * (w⊥ rWC); A·T(wH)↔A·T*(w⊥ H); A·T(wrH)↔A· T O 2 * (w⊥ rH). All of them occur via the transition states as tight ion pairs (A+, protonated by the N6H2 amino group)·(T-, deprotonated by the N3H group) with quasi-orthogonal geometry, which are stabilized by the participation of the strong (A)N6+H···O4-/O2-(T) and (A)N6+H···N3-(T) H-bonds. Established tautomerizations proceed through a two-step mechanism of the protons moving in the opposite directions along the intermolecular H-bonds. Initially, proton moves from the N3H imino group of T to the N6H2 amino group of A and then subsequently from the protonated N6+H3 amino group of A to the O4/O2 oxygen atom of T, leading to the products - A·T*(w⊥ WC), A· T O 2 * (w⊥ rWC), A·T*(w⊥ H), and A· T O 2 * (w⊥ rH), which are substantially non-planar, conformationally-labile complexes. These mispairs are stabilized by the participation of the (A)N6H/N6H'···N3(T) and (T)O2H/O4H···N6(A) H-bonds, for which the pyramidalized amino group of A is their donor and acceptor. The Gibbs free energy of activation of these mutagenic tautomerizations lies in the range of 27.8-29.8 kcal·mol-1 at T = 298.15 K in the continuum with ε = 1.

Non-dissociative structural transitions of the Watson-Crick and reverse Watson-Crick А·Ð¢ DNA base pairs into the Hoogsteen and reverse Hoogsteen forms.

In this study it was theoretically shown that discovered by us recently (Brovarets' et al., Frontiers in Chemistry, 2018, 6:8; doi: 10.3389/fchem.2018.00008) high-energetical, significantly non-planar (symmetry C1), short-lived wobbled conformers of the classical Watson-Crick А·Ð¢(WC), reverse Watson-Crick А·Ð¢(rWC), Hoogsteen А·Ð¢(Н) and reverse Hoogsteen А·Ð¢(rН) DNA base pairs are the intermediates of their pairwise А∙Т(WC)/А∙Т(rWC) ↔ А∙Т(H)/А∙Т(rH) conformational transformations. These transitions do not require for their realization the energy-consumable anisotropic rotation of the amino group of A around the exocyclic C6-N6 bond. They are controlled by the non-planar transition states with quasi-orthogonal geometry (symmetry C1) joined by the single intermolecular (Т)N3H···N6(А) H-bond (~4 kcal∙mol-1). The Gibbs free energies of activation for these non-dissociative, dipole-active conformational transitions consist 7.33 and 7.81 kcal∙mol-1, accordingly. Quantum-mechanical (QM) calculations in combination with Bader's quantum theory of "Atoms in Molecules" (QTAIM) have been performed at the MP2/aug-cc-pVDZ//B3LYP/6-311++G(d,p) level of QM theory in the continuum with ε = 4 under normal conditions.

Novel pathway for mutagenic tautomerization of classical А∙Т DNA base pairs via sequential proton transfer through quasi-orthogonal transition states: A QM/QTAIM investigation.

In this paper we have theoretically predicted a novel pathway for the mutagenic tautomerization of the classical A∙T DNA base pairs in the free state, the Watson-Crick A·Ð¢(WC), reverse Watson-Crick A·Ð¢(rWC), Hoogsteen A·Ð¢(H) and reverse Hoogsteen A·Ð¢(rH) pairs, via sequential proton transfer accompanied by a significant change in the mutual orientation of the bases. Quantum-mechanical (QM) calculations were performed at the MP2/aug-cc-pVDZ//B3LYP/6-311++G(d,p) level in vacuum phase, along with Bader's quantum theory of Atoms in Molecules (QTAIM). These processes involve transition states (TSs) with quasi-orthogonal structures (symmetry C1), which are highly polar, tight ion pairs (A-, N6H2-deprotonated)∙(T+, O4/O2-protonated). Gibbs free energies of activation for the A∙T(WC) / A∙T(rWC) ↔ A*∙Т(rwWC) / A*∙Т(wWC) tautomeric transitions (~43.5 kcal∙mol-1) are lower than for the A∙T(H) / A∙T(rH) ↔ A*N7∙Т(rwH) / A*N7∙Т(wH) tautomerisations (~53.0 kcal∙mol-1) (rare tautomers are marked by an asterisk; w-wobble configured tautomerisation products). The (T)N3+H⋯N1-(A), (T)O4+H⋯N1-(A) / (T)N3+H⋯N1-(A) and (T)O2+H⋯N1-(A) H-bonds are found in the transition states TSA-·T+A·T(WC)↔A*·T(rwWC) / TSA-·T+A·T(rWC)↔A*·T(wWC). However, in the transition state TSA-·T+A·Ð¢(H)↔A*N7·T(rwH) / TSA-·T+A·Ð¢(rH)↔A*N7·T(wH), the (T)N3+H⋯N7-(A), (T)O4+H⋯N7-(A) / (T)N3+H⋯N7-(A) and (T)O2+H⋯N7-(A) H-bonds are supplemented by the attractive (T)O4+/O2+⋯N6-(A) van der Waals contacts. It was demonstrated that the products of the tautomerization of the classical A∙T DNA base pairs-A*∙Т(rwWC), A*N7∙Т(rwH) and A*N7∙Т(wH) (symmetry Cs)-further transform via double proton transfer into the energetically favorable wobble A∙T*(rwWC), A∙T*(rwH) and A∙T*O2(wH) base mispairs (symmetry Cs).

Surprising Conformers of the Biologically Important A·T DNA Base Pairs: QM/QTAIM Proofs.

For the first time novel high-energy conformers-A·T(wWC) (5.36), A·T(wrWC) (5.97), A·T(wH) (5.78), and A·T(wrH) (ΔG = 5.82 kcal·mol-1) (See Graphical Abstract) were revealed for each of the four biologically important A·T DNA base pairs - Watson-Crick A·T(WC), reverse Watson-Crick A·T(rWC), Hoogsteen A·T(H) and reverse Hoogsteen A·T(rH) at the MP2/aug-cc-pVDZ//B3LYP/6-311++G(d,p) level of quantum-mechanical theory in the continuum with ε = 4 under normal conditions. Each of these conformers possesses substantially non-planar wobble (w) structure and is stabilized by the participation of the two anti-parallel N6H/N6H'…O4/O2 and N3H…N6 H-bonds, involving the pyramidalized amino group of the A DNA base as an acceptor and a donor of the H-bonding. The transition states - TSA·T(WC)↔A·T(wWC), TSA·T(rWC)↔A·T(wrWC), TSA·T(H)↔A·T(wH), and TSA·T(rH)↔A·T(wrH), controlling the dipole-active transformations of the conformers from the main plane-symmetric state into the high-energy, significantly non-planar state and vice versa, were localized. They also possess wobble structures similarly to the high-energy conformers and are stabilized by the participation of the N6H/N6H'…O4/O2 and N3H…N6 H-bonds. Discovered conformers of the A·T DNA base pairs are dynamically stable short-lived structures [lifetime τ = (1.4-3.9) ps]. Their possible biological significance and future perspectives have been briefly discussed.

A QM/QTAIM detailed look at the Watson-Crick↔wobble tautomeric transformations of the 2-aminopurine·pyrimidine mispairs.

This work is devoted to the careful QM/QTAIM analysis of the evolution of the basic physico-chemical parameters along the intrinsic reaction coordinate (IRC) of the biologically important 2AP·T(WC)↔2AP·T*(w) and 2AP·C*(WC)↔2AP·C(w) Watson-Crick(WC)↔wobble(w) tautomeric transformations obtained at each point of the IRC using original authors' methodology. Established profiles reflect the high similarity between the courses of these processes. Basing on the scrupulous analysis of the profiles of their geometric and electron-topological parameters, it was established that the dipole-active WC↔w tautomerizations of the Watson-Crick-like 2AP·T(WC)/2AP·C*(WC) mispairs, stabilized by the two classical N3H⋯N1, N2H⋯O2 and one weak C6H⋯O4/N4 H-bonds, into the wobble 2AP·T*(w)/2AP·C(w) base pairs, respectively, joined by the two classical N2H⋯N3 and O4/N4H⋯N1 H-bonds, proceed via the concerted stepwise mechanism through the sequential intrapair proton transfer and subsequent large-scale shifting of the bases relative each other, through the planar, highly stable, zwitterionic transition states stabilized by the participation of the four H-bonds - N1+H⋯O4-/N4-, N1+H⋯N3-, N2+H⋯N3-, and N2+H⋯O2-. Moreover, it was found out that the 2AP·T(WC)↔2AP·T*(w)/2AP·C*(WC)↔2AP·C(w) tautomerization reactions occur non-dissociatively and are accompanied by the consequent replacement of the 10 unique patterns of the specific intermolecular interactions along the IRC. Obtained data are of paramount importance in view of their possible application for the control and management of the proton transfer, e.g. by external electric or laser fields.

The A·T(rWC)/A·T(H)/A·T(rH) ↔ A·T*(rwWC)/A·T*(wH)/A·T*(rwH) mutagenic tautomerization via sequential proton transfer: a QM/QTAIM study.

In this study for the first time we have revealed by QM and QTAIM calculations at the MP2/aug-cc-pVDZ//B3LYP/6-311++G(d,p) level of QM theory the novel routes of the mutagenic tautomerization of three biologically important A·T DNA base pairs - reverse Watson-Crick A·T(rWC), Hoogsteen A·T(H) and reverse Hoogsteen A·T(rH) - followed by their rebuilding into the wobble (w) A·T*(rwWC), A·T*(wH) and A·T*(rwH) base mispairs by the participation of the mutagenic tautomers of the DNA bases (denoted by asterisk) and vice versa, thus complementing the physico-chemical property of the canonical A·T(WC) Watson-Crick DNA base pair reported earlier (Brovarets' et al., RSC Adv., 2015, 5, 99594-99605). These non-dissociative tautomeric transformations in the classical A·T(rWC), A·T(H) and A·T(rH) DNA base pairs proceed similarly to the canonical A·T(WC) DNA base pair via the intrapair sequential proton transfer with shifting towards major or minor grooves of DNA followed by further double proton transfer along the intermolecular H-bonds and are controlled by the plane symmetric and highly stable transition states - tight ion pairs formed by the A+ nucleobase, protonated by the N1/N7 nitrogen atoms, and T- nucleobase, deprotonated by the N3H imino group. Comparison of the estimated populations of the tautomerised states (10-21 to 10-14) with similar characteristics for the canonical A·T(WC) DNA base pair (10-8 to 10-7) leads authors to the conclusion, that only a base pair with WC architecture can be a building block of the DNA macromolecule as a genetic material, which is able for the evolutionary self-development. Among all four classical DNA base pairs, only A·T(WC) DNA base pair can ensure the proper rate of the spontaneous point errors of replication in DNA.

Physico-chemical profiles of the wobble ↔ Watson-Crick G*·2AP(w) ↔ G·2AP(WC) and A·2AP(w) ↔ A*·2AP(WC) tautomerisations: a QM/QTAIM comprehensive survey.

This study is intended to clarify in detail the tautomeric transformations of the wobble (w) G*·2AP(w) and A·2AP(w) nucleobase mispairs involving 2-aminopurine (2AP) into the Watson-Crick (WC) G·2AP(WC) and A*·2AP(WC) base mispairs (asterisks denote mutagenic tautomers of the DNA bases), respectively, by quantum-mechanical methods and Bader's Quantum Theory of Atoms in Molecules. Our previously reported methodology has been used, which allows the evolution of the physico-chemical parameters to be tracked along the entire internal reaction coordinate (IRC), not exclusively in the stationary states of these reactions. These biologically important G*·2AP(w) ↔ G·2AP(WC) and A·2AP(w) ↔ A*·2AP(WC) w ↔ WC tautomerisations, which are involved in mutagenic tautomerically-conformational pathways, determine the origin of the transitions and transversions induced by 2AP. In addition, it is established that they proceed through planar, highly stable, zwitterionic transition states and they exhibit similar physico-chemical profiles and stages of sequential intrapair proton transfer, followed by spatial rearrangement of the nucleobases relative to each other within the base pairs. These w ↔ WC tautomerisations occur non-dissociatively and are accompanied by a significant alteration in geometry (from wobble to Watson-Crick and vice versa) and redistribution of the specific intermolecular interactions, which can be divided into 10 patterns including AHB H-bonds and loosened A-H-B covalent bridges along the IRC of tautomerisation. Based on the redistribution of the geometrical and electron-topological parameters of the intrapair hydrogen bonds, exactly 9 key points have been allocated to characterize the evolution of these reactions.

Whether 2-aminopurine induces incorporation errors at the DNA replication? A quantum-mechanical answer on the actual biological issue.

In this paper, we consider the mutagenic properties of the 2-aminopurine (2AP), which has intrigued molecular biologists, biophysicists and physical chemists for a long time and been widely studied by both experimentalists and theorists. We have shown for the first time using QM calculations, that 2AP very effectively produces incorporation errors binding with cytosine (C) into the wobble (w) C·2AP(w) mispair, which is supported by the N4H⋯N1 and N2H⋯N3 H-bonds and is tautomerized into the Watson-Crick (WC)-like base mispair C*·2AP(WC) (asterisk denotes the mutagenic tautomer of the base), that quite easily in the process of the thermal fluctuations acquires enzymatically competent conformation. 2AP less effectively produces transversions forming the wobble mispair with A base - A·2AP(w), stabilized by the participation of the N6H⋯N1 and N2H⋯N1 H-bonds, followed by further tautomerization A·2AP(w) â†’ A*·2AP(WC) and subsequent conformational transition A*·2AP(WC) â†’ A*·2APsyn thus acquiring enzymatically competent structure. In this case, incorporation errors occur only in those case, when 2AP belongs to the incoming nucleotide. Thus, answering the question posed in the title of the article, we affirm for certain that 2AP induces incorporation errors at the DNA replication. Obtained results are consistent well with numerous experimental data.

Proton tunneling in the A∙T Watson-Crick DNA base pair: myth or reality?

The results and conclusions reached by Godbeer et al. in their recent work, that proton tunneling in the A∙T(WC) Watson-Crick (WC) DNA base pair occurs according to the Löwdin's (L) model, but with a small (~10(-9)) probability were critically analyzed. Here, it was shown that this finding overestimates the possibility of the proton tunneling at the A∙T(WC)↔A*∙T*(L) tautomerization, because this process cannot be implemented as a chemical reaction. Furthermore, it was outlined those biologically important nucleobase mispairs (A∙A*↔A*∙A, G∙G*↔G*∙G, T∙T*↔T*∙T, C∙C*↔C*∙C, H∙H*↔H*∙H (H - hypoxanthine)) - the players in the field of the spontaneous point mutagenesis - where the tunneling of protons is expected and for which the application of the model proposed by Godbeer et al. can be productive.