On the Mechanism of the Mutagenic Action of 5-Bromouracil: A DFT Study of Uracil and 5-Bromouracil in a Water Cluster

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.

Water-Controlled Keto-Enol Tautomerization of a Prebiotic Nucleobase

Barbituric acid is believed to be a proto-RNA nucleobase that was used for biological information transfer on prebiotic earth before DNA and RNA in their present forms evolved. Nucleobases have various tautomeric forms and the relative stability of these forms is critical to their biological function. It has been shown that barbituric acid has a tri-keto form in the gas phase and an enol form in the solid state. However, its dominant tautomeric form in aqueous medium that is most relevant for biology has been investigated only to a limited extent and the findings are inconclusive. We have used multiple approaches, namely, molecular dynamics, quantum chemistry, NMR, and IR spectroscopy to determine the most stable tautomer of barbituric acid in solution. We find a delicate balance in the stability of the two tautomers, tri-keto and enol, which is tipped toward the enol as the extent of solvation by water increases.

Biomolecules of 2-Thiouracil, 4-Thiouracil and 2,4-Dithiouracil: A DFT Study of the Hydration, Molecular Docking and Effect in DNA:RNAMicrohelixes

The molecular structure of 2-thiouracil, 4-thiouracil and 2,4-dithiouracil was analyzed under the effect of the first and second hydration shell by using the B3LYP density functional (DFT) method, and the results were compared to those obtained for the uracil molecule. A slight difference in the water distribution appears in these molecules. On the hydration of these molecules several trends in bond lengths and atomic charges were established. The ring in uracil molecule appears easier to be deformed and adapted to different environments as compared to that when it is thio-substituted. Molecular docking calculations of 2-thiouracil against three different pathogens: Bacillus subtilis, Escherichia coli and Candida albicans were carried out. Docking calculations of 2,4-dithiouracil ligand with various targeted proteins were also performed. Different DNA: RNA hybrid microhelixes with uridine, 2-thiouridine, 4-thiouridine and 2,4-dithiouridine nucleosides were optimized in a simple model with three nucleotide base pairs. Two main types of microhelixes were analyzed in detail depending on the intramolecular H-bond of the 2'-OH group. The weaker Watson-Crick (WC) base pair formed with thio-substituted uracil than with unsubstituted ones slightly deforms the helical and backbone parameters, especially with 2,4-dithiouridine. However, the thio-substitution significantly increases the dipole moment of the A-type microhelixes, as well as the rise and propeller twist parameters.

Effect of the sulfur atom on S2 and S4 positions of the uracil ring in different DNA:RNA hybrid microhelixes with three nucleotide base pairs

The effect of the sulphur atom on the uracil ring was analyzed in different DNA:RNA microhelixes with three nucleotide base-pairs, including uridine, 2-thiouridine, 4-thiouridine, 2,4-dithiouridine, cytidine, adenosine and guanosine. Distinct backbone and helical parameters were optimized at different density functional (DFT) levels. The Watson-Crick pair with 2-thiouridine appears weaker than with uridine, but its interaction with water molecules appears easier. Two types of microhelixes were found, depending on the H-bond of H2' hydroxyl atom: A-type appears with the ribose ring in ³ E-envelope C_3' -endo, and B-type in ² E-envelope C_2' -endo. B-type is less common but it is more stable and with higher dipole-moment. The sulphur atoms significantly increase the dipole-moment of the microhelix, as well as the rise and propeller twist parameters. Simulations with four Na atoms H-bonded to the phosphate groups, and further hydration with explicit water molecules were carried out. A re-definition of the numerical value calculation of several base-pair and base-stacking parameters is suggested.

Energetic, Topological and Electric Field Analyses of Cation-Cation Nucleic Acid Interactions in Watson-Crick Disposition

A theoretical study of the effect of the diprotonation on the nucleic acid bases (A : U, A : T and G : C) in Watson-Crick conformation has been carried out by means of DFT computational methods in vacuum. In addition, the corresponding neutral and monoprotonated binary complexes have been considered. Most of the diprotonated species studied are stable, even though the binding energy is positive due to the overall repulsive electrostatic term. Local electrostatic attractive forces in the regions of hydrogen bonds (HBs) are responsible for equilibrium geometries, as shown by the electric field lines connecting the electrophilic and nucleophilic sites involved in the HB interactions. Secondary electrostatic effects also affect the assembling of the nucleic acid complexes in either neutral or cationic form. In particular, the electric field lines flowing from electrophilic sites in one base to nucleophilic sites in the other reinforce the linking between them. Hence, when the nucleophilic site concerns the free lone pair of the heteroatom involved in the HB interaction as acceptor, the HB distance shortens. However, if the free lone pair of the HB acceptor interacts with an electrophilic site in the same molecule, the HB distance elongates, weakening the HB interaction. The topological analysis of the electron density distribution in HB regions indicates that neutral, monoprotonated and diprotonated complexes show no differences in the nature of their HB's.

Relative Populations of Some Tautomeric Forms of 2'-Deoxyguanosine-5-Fluorouridine Mismatch

The importance of the 2'-deoxyguanosine-uridine mispair as the most occurring mismatch in transcriptional studies of RNAs from DNAs is multiplied when 5-halo-substituted uridine species cause a serious increase in the probability of its occurrence. Many studies relate this higher probability to the existence of possible tautomeric and ionic forms of its constituent bases. According to these statements, relative populations of mismatches between 5-fluorouridine and both keto and enol forms of 2'-deoxyguanosine are computed by using a conformational search. In order to have a complete scan of all of the highly probable conformers in a moderate computational time, an extensive conformational search methodology is employed here, which benefits from the advantages of both the molecular dynamics simulations and quantum mechanics calculations. The population of an enolic tautomer of normal wobble orientation is about 0.057% of that of its keto tautomer, whereas the population of an enolic tautomer of reverse wobble orientation is about 0.0054% of that of its keto tautomer. Totally, the reverse wobble orientation is about six times more populated than the normal wobble orientation. Calculated populations are in good agreement with experimental populations of closely related compounds. The reliability of the applied methodology is certified, in part, by a good agreement obtained between some experimental data and corresponding Boltzmann-weighted average data of most probable conformers such as NMR parameters. The validation of this methodology is certified with high accuracy by applying it on the substituted diuridine pairs, where experimental populations are available. Not only are the calculated populations and NMR parameters of this test in very good agreement with the experimental data, but also they are free of the ambiguities mentioned by experimentalists.

Tautomerism of uracil and related compounds: A mass spectrometry study

It has been demonstrated that uracil has a preponderant tautomeric form, but it is also known that different tautomers co-exist in this equilibrium. In this work, mass spectrometry is used as a helpful tool to analyse the equilibria, using derivative compounds to forbid the presence of some tautomers and ion trap mass spectrometry to follow relevant fragmentation pathways. Theoretical calculations were performed to confirm tautomers abundance by energy minimization in gas phase. Analysis of mass spectra of uracil, three methyl-substituted uracils, 2-thiouracil and three benzouracils suggest that uracil exists mainly as three tautomers in gas phase: one major structure that corresponds to the classical structure of uracil (pyrimidine-2,4(1H,3H)-dione) bearing two carbonyls and two NH moieties, and two minor enolic forms (4-hydroxypyrimidin-2(1H)-one and 2-hydroxypyrimidin-4(1H)-one). Such tautomeric distribution is supported by theoretical calculations, which show that they are the three most stable tautomers.

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.

Clustering of Uracil Molecules on Ice Nanoparticles

We generate a molecular beam of ice nanoparticles (H₂O)_N, N̅ ≈ 130-220, which picks up several individual gas phase uracil (U) or 5-bromouracil (BrU) molecules. The mass spectra of the doped nanoparticles prove that the uracil and bromouracil molecules coagulate to clusters on the ice nanoparticles. Calculations of U and BrU monomers and dimers on the ice nanoparticles provide theoretical support for the cluster formation. The (U)_mH⁺ and (BrU)_mH⁺ intensity dependencies on m extracted from the mass spectra suggest a smaller tendency of BrU to coagulate compared to U, which is substantiated by a lower mobility of bromouracil on the ice surface. The hydrated U_m·(H₂O)_nH⁺ series are also reported and discussed. On the basis of comparison with the previous experiments, we suggest that the observed propensity for aggregation on ice nanoparticles is a more general trend for biomolecules forming strong hydrogen bonds. This, together with their mobility, leads to their coagulation on ice nanoparticles which is an important aspect for astrochemistry.

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.

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.

Competition between hydrogen and halogen bonding in halogenated 1-methyluracil: Water systems

The competition between hydrogen- and halogen-bonding interactions in complexes of 5-halogenated 1-methyluracil (XmU; X = F, Cl, Br, I, or At) with one or two water molecules in the binding region between C5-X and C4=O4 is investigated with M06-2X/6-31+G(d). In the singly-hydrated systems, the water molecule forms a hydrogen bond with C4=O4 for all halogens, whereas structures with a halogen bond between the water oxygen and C5-X exist only for X = Br, I, and At. Structures with two waters forming a bridge between C4=O and C5-X (through hydrogen- and halogen-bonding interactions) exist for all halogens except F. The absence of a halogen-bonded structure in singly-hydrated ClmU is therefore attributed to the competing hydrogen-bonding interaction with C4=O4. The halogen-bond angle in the doubly-hydrated structures (150-160°) is far from the expected linearity of halogen bonds, indicating that significantly non-linear halogen bonds may exist in complex environments with competing interactions.

Structural grounds for the 2-aminopurine mutagenicity: a novel insight into the old problem of the replication errors

Mutagenic pressure of the 2AP molecule on DNA during its replication is realized via the more intensive generation of the T* mutagenic tautomers through the reaction 2AP·T(WC) → 2AP·T*(w).

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.

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.

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.

[Structural and energetic properties of the four configurations of the A.T and G.C DNA base pairs]

Using the methods of non-empirical quantum chemistry at the MP2/6-311++G(2df,pd)// B3LYP/6-311++G(d,p) level of theory it was established for the first time, that Hoogsteen, reverse Hoogsteen, Watson-Crick and reverse Watson-Crick configurations of the A.T and G.C DNA base pairs are isoelectronic and isomorphic structures with similar dynamic properties. Based on these results, non-ionisation mechanism of the Hoogsteen <"breathing" of the G*.C* DNA base pair, namely transformation of the tautomerised (Lowdin's) G-C base pair with Watson-Crick geometry into the Hoogsteen electroneutral G*.C* H base pair stabilized by the three O6H...N4, N3H...N7 and C8H...02 H-bonds, was postulated. It is suggested that such scenario activates only in those cases, when DNA is not located in aqueous solution, but works together with proteins and cytosine protonation at the N3 atom is precluded.

[Under what conditions does G.C Watson-Crick DNA base pair acquire all four configurations characteristic for A.T Watson-Crick DNA base pair?]

At the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of theory it was established for the first time, that the Löwdin's G*.C* DNA base pair formed by the mutagenic tautomers can acquire, as the A-T Watson-Crick DNA base pair, four biologically important configurations, namely: Watson-Crick, reverse Watson-Crick, Hoogsteen and reverse Hoogsteen. This fact demonstrates rather unexpected role of the tautomerisation of the one of the Watson-Crick DNA base pairs, in particular, via double proton transfer: exactly the G.C-->G*.C* tautomerisation allows to overcome steric hindrances for the implementation of the above mentioned configurations. Geometric, electron-topological and energetic properties of the H-bonds that stabilise the studied pairs, as well as the energetic characteristics of the latters are presented.

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.

Atomistic understanding of the C·T mismatched DNA base pair tautomerization via the DPT: QM and QTAIM computational approaches

It was established that the cytosine·thymine (C·T) mismatched DNA base pair with cis-oriented N1H glycosidic bonds has propeller-like structure (|N3C4C4N3| = 38.4°), which is stabilized by three specific intermolecular interactions-two antiparallel N4H…O4 (5.19 kcal mol(-1)) and N3H…N3 (6.33 kcal mol(-1)) H-bonds and a van der Waals (vdW) contact O2…O2 (0.32 kcal mol(-1)). The C·T base mispair is thermodynamically stable structure (ΔG(int) = -1.54 kcal mol(-1) ) and even slightly more stable than the A·T Watson-Crick DNA base pair (ΔG(int) = -1.43 kcal mol(-1)) at the room temperature. It was shown that the C·T ↔ C*·T* tautomerization via the double proton transfer (DPT) is assisted by the O2…O2 vdW contact along the entire range of the intrinsic reaction coordinate (IRC). The positive value of the Grunenberg's compliance constants (31.186, 30.265, and 22.166 Å/mdyn for the C·T, C*·T*, and TS(C·T ↔ C*·T*), respectively) proves that the O2…O2 vdW contact is a stabilizing interaction. Based on the sweeps of the H-bond energies, it was found that the N4H…O4/O4H…N4, and N3H…N3 H-bonds in the C·T and C*·T* base pairs are anticooperative and weaken each other, whereas the middle N3H…N3 H-bond and the O2…O2 vdW contact are cooperative and mutually reinforce each other. It was found that the tautomerization of the C·T base mispair through the DPT is concerted and asynchronous reaction that proceeds via the TS(C·T ↔ C*·T*) stabilized by the loosened N4-H-O4 covalent bridge, N3H…N3 H-bond (9.67 kcal mol(-1) ) and O2…O2 vdW contact (0.41 kcal mol(-1) ). The nine key points, describing the evolution of the C·T ↔ C*·T* tautomerization via the DPT, were detected and completely investigated along the IRC. The C*·T* mispair was revealed to be the dynamically unstable structure with a lifetime 2.13·× 10(-13) s. In this case, as for the A·T Watson-Crick DNA base pair, activates the mechanism of the quantum protection of the C·T DNA base mispair from its spontaneous mutagenic tautomerization through the DPT.

Can DNA-binding proteins of replisome tautomerize nucleotide bases? Ab initio model study

Ab initio quantum-chemical study of specific point contacts of replisome proteins with DNA modeled by acetic acid with canonical and mutagenic tautomers of DNA bases methylated at the glycosidic nitrogen atoms was performed in vacuo and continuum with a low dielectric constant (ϵ ∼ 4) corresponding to a hydrophobic interface of protein-nucleic acid interaction. All tautomerized complexes were found to be dynamically unstable, because the electronic energies of their back-reaction barriers do not exceed zero-point vibrational energies associated with the vibrational modes whose harmonic vibrational frequencies become imaginary in the transition states of the tautomerization reaction. Additionally, based on the physicochemical arguments, it was demonstrated that the effects of biomolecular environment cannot ensure dynamic stabilization. This result allows suggesting that hypothetically generated by DNA-binding proteins of replisome rare tautomers will have no impact on the total spontaneous mutation due to the low reverse barrier allowing a quick return to the canonical form.

Vibrational spectra, tautomerism and thermodynamics of anticarcinogenic drug: 5-fluorouracil

The FT-IR and FT-Raman spectra of 5-Fluorouracil were recorded in the solid phase in the regions 400-4000 cm(-1) and 50-4000 cm(-1), respectively. The vibrational spectra were analysed and the observed fundamentals were assigned to different normal modes of vibration. The experimental wavenumbers were compared with the scaled vibrational values using DFT methods: the Ar matrix data were related to gas phase calculations, while the values of the solid state spectra were compared to those with dimer simulations. The study indicates that some features that are characteristic of vibrational spectra of uracil and its derivatives are retained in the spectrum of 5-fluorouracil and it exists in ketonic form in the solid phase. The tautomerism was also studied and the spectra of the two most stable forms were simulated. The calculated wavenumbers have been employed to yield thermodynamic properties.

A study of nucleic acid base-stacking by the Monte Carlo method: Extended cluster approach

The adenine-thymine (AT), adenine-uracil (AU) and guanine-cytosine (GC) base associates in clusters containing 400 water molecules were studied using a newly implemented Metropolis Monte Carlo algorithm based on the extended cluster approach. Starting from the hydrogen-bonded Watson-Crick geometries, all three base pairs are transformed into more favorable stacked configurations during the simulation. The obtained results show, for the first time, the transition from planar base pairs to stacked base associates in the Monte Carlo framework. Analysis of the interaction energies shows that, in the water cluster, the stacked dimers are energetically preferable compared to the corresponding Watson-Crick base pairs. This is due to the larger base-water interaction in the stacked structures. The water-water interaction is one of the main factors promoting the formation of stacked dimers, and the obtained data confirm the crucial role of the water-water interactions in base stacking.

Tautomeric equilibrium of uracil and thymine in model protein-nucleic acid contacts. Spectroscopic and quantum chemical approach

This work deals with tautomeric transformations of uracil (Ura) and thymine (Thy) in their model complexes with the deprotonated carboxylic group. Essential changes in the UV spectra of the bases upon their interaction with NaAc, vanishing signals of both imino protons in (1)H NMR spectra, and a perceptible decrease in intensity of both IR bands, related to the stretching vibrations nu(C=O) of the carbonyl groups, imply involvement of enolic tautomers. Results of quantum chemical calculations of the double complexes of the Ura(Thy) tautomers with CH(3)COO(-) at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of theory proved to be incompatible with the spectral features: despite the fact that the complexes of the enolic tautomers are much closer in energy to the diketo ones as compared to isolated tautomers, the energy gap between them is such that in tautomeric equilibrium dominate diketo forms. Calculations of triple complexes of the type CH(3)COO(-):Ura(Thy) tautomer:Na(+), taking into account the effect of the Na(+) coordination with tautomers, show that three triple complexes formed by enolic tautomers appeared more stable than those formed by diketo ones. This makes the UV and (1)H NMR data understandable, but the high residual intensity of the nu(C=O) bands in the IR spectra remains unclear. At that ion, Na(+) itself was not able to disturb the tautomeric equilibrium in the coordination complexes of the type Ura(Thy) tautomer:Na(+). To evaluate the DMSO effect, the CPCM solvation model was applied to triple complexes of the Ura tautomers. It appeared that in the solution there is coexistence between the diketo and enolic tautomers in a ratio of 53%:47%. This makes possible reconciliation of our experimental data. The biological significance of high-energy tautomers of nucleotide bases is discussed.

Replication across template T/U by human DNA polymerase-iota

Human DNA polymerase-iota (Poliota) incorporates correct nucleotides opposite template purines with a much higher efficiency and fidelity than opposite template pyrimidines. In fact, the fidelity opposite template T is so poor that Poliota inserts an incorrect dGTP approximately 10 times better than it inserts the correct dATP. We determine here how a template T/U is accommodated in the Poliota active site and why a G is incorporated more efficiently than an A. We show that in the absence of incoming dATP or dGTP (binary complex), template T/U exists in both syn and anti conformations, but in the presence of dATP or dGTP (ternary complexes), template T/U is predominantly in the anti conformation. We also show that dATP and dGTP insert differently opposite template T/U, and that the basis of selection of dGTP over dATP is a hydrogen bond between the N2 amino group of dGTP and Gln59 of Poliota.