Manifestations of intramolecular H-bonds of CH O and OH C type in quercetin molecule: Analysis of IR spectra by mean of density functional theory.

The IR spectra of 48 conformers of quercetin which represent full conformation space of its tautomers have been modeled at B3LYP/6-311++G(d,p) level of the density functional theory. The presence of intramolecular H-bonds C2'H/C6'H…O3 and O3H…C2'/C6' was characterized by their spectral manifestations. The C2'H/C6'H…O3 contacts were found to have a spectral blue-shift. The O3H…C2'/C6' contacts were mostly red-shifted. The stretching vibrations of H-bonds C2'H/C6'H…O3 demonstrate an increase in the intensity of the modes of stretching vibrations ν(C2'H)/ν(C6'H) and an increase in the frequency of their out-of-plane vibrations γ(C2'H)/γ(C6'H). Most of the spectral parameters correlate a little with the energy of the H-bonds.

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.

Effect of Microenvironment on the Geometrical Structure of d(A)5 d(T)5 and d(G)5 d(C)5 DNA Mini-Helixes and the Dickerson Dodecamer: A Density Functional Theory Study.

We report a comprehensive quantum-chemical study on d(A)5·d(T)5 and d(G)5·d(C)5 DNA mini-helixes and the Dickerson dodecamer d[CGCGAATTCGCG]. The research was performed to model the evolution of the spatial structure of d(A)5·d(T)5 and d(G)5 d(C)5 DNA mini-helixes all the way from vacuum to water bulk. The influence of external factors such as the presence of counterions and the extent of hydration was included. Also, for comparison, limited calculations have been carried out on the Dickerson dodecamer. The study has been performed at the density functional theory level using B97D3 and ωB97XD exchange-correlation functionals augmented by the Def2SVP basis set. We found that the (dA)5·(dT)5 anion when placed in vacuum forms a DNA duplex, which possesses an intermediate form between a helix and a ladder. The presence of compensating Na+ counterions or explicit microhydration of minor and major grooves stabilizes a DNA mini-helix of B-shape. Factors such as water bulk play a minor role. Somewhat different behavior has been found in the case of the (dG)5·(dC)5 duplex. In this case, we observe the formation of B-type mini-helixes even for the (dG)5·(dC)5 anion placed in vacuum. This is due to an additional stabilization originated from the appearance of an extra hydrogen bond, compared to an AT base pair. To assess whether the obtained results are transferable to different sizes of mini-helixes, similar calculations have been performed for the duplex formed by the Dickerson dodecamer which contains a total of 12 dG·dC and dA·dT base pairs. It has been found that in vacuum, analogous to the d(A)5·d(T)5 duplex, this system possesses a shape which is also quite close to a ladder. However, the presence of factors such as hydration restores the B-type geometry. Also, our results completely in line with the results of electrospray-ionization experiments suggest that uncompensated by counterions the DNA backbone preserves the duplex geometry in vacuum. We present arguments that this state is kinetically unstable.

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.

Structural transitions in poly(A), poly(C), poly(U), and poly(G) and their possible biological roles.

The homopolynucleotide (homo-oligonucleotide) tracts function as regulatory elements at various stages of mRNAs life cycle. Numerous cellular proteins specifically bind to these tracts. Among them are the different poly(A)-binding proteins, poly(C)-binding proteins, multifunctional fragile X mental retardation protein which binds specifically both to poly(G) and poly(U) and others. Molecular mechanisms of regulation of gene expression mediated by homopolynucleotide tracts in RNAs are not fully understood and the structural diversity of these tracts can contribute substantially to this regulation. This review summarizes current knowledge on different forms of homoribopolynucleotides, in particular, neutral and acidic forms of poly(A) and poly(C), and also biological relevance of homoribopolynucleotide (homoribo-oligonucleotide) tracts is discussed. Under physiological conditions, the acidic forms of poly(A) and poly(C) can be induced by proton transfer from acidic amino acids of proteins to adenine and cytosine bases. Finally, we present potential mechanisms for the regulation of some biological processes through the formation of intramolecular poly(A) duplexes.

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.