Where Quantum Biochemistry Meets Structural Bioinformatics: Excited Conformationally- Tautomeric States of the Classical A T DNA Base Pair

This Chapter summarizes recent quantum-chemical (QM) investigations of the novel conformational and tautomeric states on the potential energy hypersurface of the classical A T/A U nucleobase pairs. For the first time, it was observed 28 local minima for each base pair excluding enantiomers planar, non-planar base pairs and structures with wobble geometry. Considered excited conformationallytautomeric states of the classical A T DNA base pair have been revealed in the Nucleic Acid Database by structural bioinformatics. These data shed light on the biological significance of the unusual A T/A U nucleobase pairs for the functioning of the nucleic acids at the quantum level.


Introduction
Since the discovery of the spatial organization of the DNA molecule by James Watson & Francis Crick [1,2], it is traditionally believed that canonical Watson-Crick AÁT and GÁC DNA base pairs are quite conservative structures. These classical DNA base pairs almost do not have tautomeric variability and essential conformational mobility at the dynamical behavior of DNA molecule [2]. Generally, it is suggested that bases in the anti-conformation able to form a pair according to the so-called Watson-Crick (WC) scheme joined through the three intermolecular hydrogen (H) bonds [3]. At the same time, many biologists-contemporaries questioned the proposed Watson-Crick conformation, since X-ray resolution does not allow to establish for sure the precise conformation of the base pairs constituting to the DNA double helix. Exactly by this reason Maurice Wilkinsthe third author of the discovery of the DNA structureexplained the reason why Rosalind Franklin doubted in modeling the structure of DNA [4]. And even forefathers of the Discussed DNA base pairs are not static structures in the composition of DNA [27,28]. Thus, the spontaneous AÁT(WC)$AÁT(Н) conformational transition has been experimentally registered by the NMR spectroscopy on the DNA regions enriched by the classical AÁT nucleobase pairs [22]. Despite numerous theoretical investigations, microstructural nature of these transitions still remains incomprehensible [20,29].
These novel excited conformational and tautomeric states occur due to the quantum effects, e.g. amino group pyramidalization because of electron conjugation of the lone electron pair of nitrogen amino atom with π-electron system of the ring [54][55][56]. This data enables us to suggest the potential energy surface of the classical АÁТ base pairs and also to predict pathways of their interconversions. Moreover, this modeling could be used for the understanding and description in details of the physico-chemical mechanisms of the DNA functioning, in particular DNA "breathing", which has significant biological role [27,28].
Also, obtained data would enable to make new insights into the understanding the DNA and RNA structural biology, which are based on their conformational and tautomeric variety. By the methods of structural bioinformatics it was revealed unusual conformationally-tautomeric states of the AÁT DNA base pair in the Nucleic Acid Database, confirming their existence in biological systems. Altogether, further this could be extended to the area of epigenetics and experimental verification.

Computational methods
Equilibrium geometries of the investigated DNA base pairs, as well as their harmonic vibrational frequencies have been calculated at the B3LYP/6-311++G(d, p) level of theory [54][55][56][57][58], using Gaussian'09 package [59]. Applied level of theory has proved itself successful for the calculations of the similar systems [60][61][62]. A scaling factor that is equal to 0.9668 has been applied in the present work for the correction of the harmonic frequencies of all complexes [63,64].
All calculations have been carried out in the continuum with ε = 1, that adequately reflects the processes occurring in real biological systems without deprivation of the structurally functional properties of the bases in the composition of DNA, and in the continuum with ε = 4, which satisfactorily models the substantially hydrophobic recognition pocket of the DNA-polymerase machinery as a part of the replisome [65][66][67][68][69][70].
The Gibbs free energy G for all structures was obtained in the following way: where E el -electronic energy, while E corr -thermal correction. Electronic interaction energies ∆E int have been calculated at the MP2/6-311++G (2df,pd) level of theory as the difference between the total energy of the base pair and energies of the monomers and corrected for the basis set superposition error (BSSE) [73,74] through the counterpoise procedure [75,76].

Bioinformatical analysis
It was created original author's algorithm in order to reveal the unusual АÁТ base pairs in the Nucleic Acid Database by Rutgers University [93,94]. This algorithm is based on the comparison of the calculated structure of the АÁТ base pairs at the ε =4 with structure of the analogical base pairs in the Nucleic Acid Database.

Obtained results
It was found out novel tautomerization pathways for the formation of the rare tautomers of the A or T DNA bases: • AÁT(WC)$A*ÁT(w)/AÁT* O2 (w)/AÁT*(w) via the sequential proton transfer and shifting of the bases relatively each other [47]; • AÁT(rWC)/AÁT(H)/AÁT(rH)$AÁT*(rw WC )/AÁT*(w H )/AÁT*(rw H ) mutagenic tautomerization via the sequential proton transfer [48]; [49]; So, on the potential energy surface of the classical AÁT/AÁU base pair it was received 28 various conformationally-tautomeric states (Figure 1 Notably, that Gibbs free and electronic energies of the AÁT/AÁU base pairs are in the wide range of values, which insignificantly decrease at the transition from the continuum with ε = 1 to the continuum with ε = 4, while dipole moment increases at this ( Table 1).
We have carefully scanned all 28 unusual conformationally-tautomeric states of the АÁТ DNA base pairs in the Nucleic Acid Database by Rutgers University using original author's methodology for structural bioinformatics analysis. It was identified most part of the theoretically investigated by us excited conformationallytautomeric states of the classical АÁТ DNA base pair (Figure 1, Table 1).

Discussion of the obtained results
Let us start the discussion and more detailed analysis of the obtained results from the consideration of the traditional area of the biological applications of the prototropic tautomerism of the DNA bases [54], as well as their role in the origin of the spontaneous point mutationstransitions and transversions at the DNA biosynthesisso-called replication errors [58][59][60][61][62][63]. This physico-chemical model should satisfy strict conditions. Saying shortly, in order to point on the most important things, from one sidebarriers of the mutagenic tautomerization of the base pairs should not be quite high in view of the quite rigid kinetic requirements for the incorporation into the double strand of DNA by the DNA-polymerase during the one act of replication ($10 À4 s) [54]. At this, the lifetime of the tautomerized states of the pairs should exceed characteristic time of the inertial DNA-polymerase machinery ($10 À9 s). Only at this condition the inertial replicational DNApolymerase machinery would successfully dissociate tautomerized base pairs into the monomers, in particular into the rare tautomeric forms.
From the other side, these barriers should be quite high in order to overcome resistance of the stacking interactions and sugar-phosphate backbone of DNA on the way of the incorporation of the tautomerizing base pair into the double structure of DNA [54].   Table 1. Figure 1).

Energetic and polar characteristics of the conformers and tautomers of the AÁT/AÁU nucleobase pairs obtained at the MP2/6-311++G(2df,pd) // B3LYP/6-311++G(d,p) level of QM/PCM theory in the isolated state (ε = 1) and in the continuum with ε = 4 under normal conditions (see
Nowadays, just one single model satisfies these strict conditions [47]. According to this model (Figure 1), mutagenic tautomerization of the bases in the AÁT(WC) base pair is controlled by the transition states, which represent itself tight ion pairs A + ÁT À , and is realized through the step-by-step proton transfer along the intermolecular H-bonds and is assisted by the lateral changing of the configuration of the pairits transition from the Watson-Crick configuration to the wobble or shifted [47]. In fact, complementary A base plays a role of catalysator of the intramolecular mutagenic tautomerization of the T base within the AÁT(WC) base pair. Below it would be outlined experimental confirmations that wobble structures of the AÁT base pair, containing mutagenic tautomeric forms of the T base, are real objects of the structural biology. This fact, in our opinion, experimentally confirms reality of the tautomeric mechanisms of the origin of the replication errors [47].
We have demonstrated for the first time, that others three biologically important configurations of the AÁT base pair -AÁT(rWC), AÁT (H) and AÁT(rH) [47] tautomerises by the abovementioned and described mechanism of the tautomerization, forming wobble pairs by the participation of the mutagenic tautomers (Figure 1). Moreover, we have arrived to the conclusion by the comparison of their energetical characteristics, that Nature quite consciously choose evolutionary the most remote AÁT(WC) base pair for the building of the carrier of the genetic information in the form of the right-handed DNA [47].
In this regard, it arises quite logical question -"Whether Nature uses prototropic tautomerization of the DNA bases beyond the borders of classical tautomeric hypothesis?" Let us sayfor the supporting of the unusual DNA structures. Principle of economy of thinking (Entia non sunt multiplicanda praeter necessitate), which is quite often applied by the living nature, enables in principle, affirmatively answer on the quite interesting question. Below we would provide number of examples of the application in the structural biology of all without exception wobble configurations of the AÁT pair by the participation of the mutagenic tautomers.
Biological role of the prototropic tautomerism of the DNA bases is not limited by the presented here examples. It is quite more complex and wider. Let us attract readers' attention to the one more so-called unusual role of the tautomericconformational transformations in the DNA structural transitions. However, their mechanism of action could be explained only at the macroscopical level.
In the work [61] at the example of the hypoxanthine dimer it was revealed novel way of the conformationally-tautomeric transformations of the structures, which are joined by the neighboring antiparallel H-bonds, through the quasi-orthogonal transition state with the changing of the mutual orientation of the dimmers on 180 degree. Conformationally-tautomeric transitions of such a nature have been fixed in all without exception four configurations of the classical AÁT DNA base pair [53]. Combining these data with previous, concerning the WC/H$w WC /w H conformationally-tautomeric transitions [50], we have obtained joined picture of the WC/H$rWC/rH at the quantum level: Bioinformatical analysis. This data convincingly evidence on the real occurrence of these base pairs in the real biological systems [93,94] and thuson their biological importance. This situation remains for a long time the hidden side of the classical АÁТ DNA base pair. However, it became successfully resolved in the current work. 9 Where Quantum Biochemistry Meets Structural Bioinformatics: Excited Conformationally… DOI: http://dx.doi.org/10.5772/intechopen.94565

Conclusions
Concluding, we can state that it was received the most complete up to now quantum map of the biologically-important conformationally-tautomeric transitions of the classical AÁT/AÁU nucleobase pair, which enable to classify it as a quantum choreography with all further going consequences. But it is not a pursuit of a new term, but rather an attempt to realize the molecular logic of the quantum evolution at its initial stages, when it was formed its behavior, which is evolutionary programmed in its electronic structure.
For the first time, it was shown for the classical AÁT DNA base pair that prototropic tautomerism of the DNA bases is responsible both for the origin, as well as for the supporting of the unusual local structures in the constitution of DNA and in complexes with proteins and small biomolecules. Moreover, prototropic tautomerism of the classical AÁT DNA base pair significantly expands its conformational possibilities and its impact on the biological importance.
It is connected with the fact that presented mechanisms of the tautomerization are assisted by the significant changing of the geometry of the tautomerizing base pair. This means that they are conformationally-tautomeric transitions by their essence.
This conclusion is confirmed by the structural bioinformatics. Thus, it was identified hundredth of the structures containing tautomers of the DNA bases. This fact points that all described exited conformationally-tautomeric states of the AÁT and AÁU nucleobase pairs, corresponding to local minima, are real structures.