Water-mediated tautomerization of cytosine to the rare imino form: An ab initio dynamics study

Tautomeric equilibrium and spectroscopic properties of 8-azaguanine revealed by quantum chemistry methods

8-azaguanine is a triazolopyrimidine nucleobase analog possessing potent antibacterial and antitumor activities, and it has been implicated as a lead molecule in cancer and malaria therapy. Its intrinsic fluorescence properties can be utilized for monitoring its interactions with biological polymers like proteins or nucleic acids. In order to better understand these interactions, it is important to know the tautomeric equilibrium of this compound. In this work, the tautomeric equilibrium of all natural neutral and anionic compound forms (except highly improbable imino-enol tautomers) as well as their methyl derivatives and ribosides was revealed by quantum chemistry methods. It was shown that, as expected, tautomers protonated at positions 1 and 9 dominate neutral forms both in gas phase and in aqueous solution. 8-azaguanines methylated at any position of the triazole ring are protonated at position 1. The computed vertical absorption and emission energies are in very good agreement with the experimental data. They confirm the validity of the assumption that replacing the proton with the methyl group does not significantly change the positions of absorption and fluorescence peaks.

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.

Non-poissonian Distribution of Point Mutations in DNA

In general, for chemical reactions occurring in systems, where fluctuations are not negligibly small, it is necessary to introduce a master equation for distribution of probability of fluctuations. It has been established that the monomolecular reactions of a type as A ↔ X are described by the master equation, which leads to a Poisson distribution with the variance equal to the average value N 0. However, the consideration of the Löwdin mechanism as autocatalytic non-linear chemical reactions such as A + X ↔ 2X and the corresponding master equation lead to a non-Poissonian probability distribution of fluctuations. In the presented work, first-order autocatalysis has been applied to the Löwdin's mechanism of spontaneous mutations in DNA. Describing double proton transfers between complimentary nucleotide bases along the chain by first-order autocatalytic reactions, the corresponding master equation for protons in tautomeric states becomes non-linear, and at non-equilibrium conditions this leads to the non-Poissonian distribution of spontaneous mutations in DNA. It is also suggested that the accumulation of large fluctuations of successive cooperative concerted protons along the chain may produce higher non-linearities which could have a significant impact on some biochemical processes, occurring in DNA.

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

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

Formic acid catalyzed isomerization of protonated cytosine: a lower barrier reaction for tautomer production of potential biological importance

Tautomerism in nucleotide bases is one of the possible mechanisms of DNA mutation. In spite of numerous studies on the structure and energy of protonated cytosine tautomers, little information is available on the process of their intra- and intermolecular tautomerizations. The catalytic ability of H₂O, HCOOH, and the HCOOHH₂O group to facilitate the tautomerism of the Cyt2t⁺ to CytN3⁺ isomer has been studied. It is shown that the activation free energies of tautomerism in the gas phase are 161.17, 58.96, 26.06, and 15.69 kJ mol^-1, respectively, when the reaction is carried out in the absence and presence of H₂O, HCOOH, or the HCOOHH₂O group. The formation of a doubly hydrogen bonded transition state is central to lowering the activation free energy and facilitating the intramolecular hydrogen atom transfer that is required for isomerization. In the aqueous phase, although the solvent effects of water significantly decrease the activation free energy of intramolecular tautomerization, the isomerization of the Cyt2t⁺ to CytN3⁺ isomer remains unfavorable, and the HCOOH and HCOOHH₂O group mediated mechanisms are still more favorable. Meanwhile, conventional transition state theory (CTST) followed by Wigner tunneling correction is then applied to estimate the rate constants. The rate constant with Wigner tunneling correction for direct tautomerization is obviously smaller than that of HCOOH-mediated tautomerization, which is the most plausible mechanism. Finally, another important finding is that the product complex (CytN3⁺HCOOH) is in the rapid tautomeric equilibrium with the reaction complex (Cyt2t⁺HCOOH) (τ_99.9% = 3.84 × 10^-12 s), which is implemented by the mechanism of the concerted synchronous double proton transfer. Its lifetime of the formed CytN3⁺HCOOH complex (τ = 8.33 × 10^-9 s) is almost one order of magnitude larger than the time required for the replication machinery to forcibly dissociate a base pair into the monomers during DNA replication (several ns), which is further dissociated into the CytN3⁺ and HCOOH monomers. The results of the present study demonstrate the feasibility of acid catalysis for DNA base isomerization reactions that would otherwise be forbidden.

Deuteration of Aromatic Rings under Very Mild Conditions through Keto-Enamine Tautomeric Amplification

We have discovered a surprising, mild method for deuteration of select aromatic compounds that is facilitated by a keto-enamine tautomeric intermediate. The mechanism of the reaction has been studied using kinetics experiments and detailed computational analysis. It was found that a chain of water molecules has a substantial role in lowering the activation barrier to the tautomerization-enhanced deuteration reaction. Our results demonstrate that tautomeric forms of aromatic molecules can be exploited to bring about enhanced reactivity.

Gas-phase hydration thermochemistry of sodiated and potassiated nucleic acid bases

Hydration reactions of sodiated and potassiated nucleic acid bases (uracil, thymine, cytosine, and adenine) produced by electrospray have been studied in a gas phase using the pulsed ion-beam high-pressure mass spectrometer. The thermochemical properties, ΔH(o)(n), ΔS(o)(n), and ΔG(o)(n), for the hydrated systems were obtained from hydration equilibrium measurement. The structural aspects of the hydrated complexes are discussed in conjunction with available literature data. The correlation between water binding energies in the hydrated complexes and the corresponding metal ion affinities of nucleobases suggests that a significant (if not dominant) amount of the canonical structure of cytosine undergoes tautomerization during electrospray ionization, and the thermochemical values for cationized cytosine probably correspond to a mixture of tautomeric complexes.

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.

Structural evidence for the rare tautomer hypothesis of spontaneous mutagenesis

Even though high-fidelity polymerases copy DNA with remarkable accuracy, some base-pair mismatches are incorporated at low frequency, leading to spontaneous mutagenesis. Using high-resolution X-ray crystallographic analysis of a DNA polymerase that catalyzes replication in crystals, we observe that a C • A mismatch can mimic the shape of cognate base pairs at the site of incorporation. This shape mimicry enables the mismatch to evade the error detection mechanisms of the polymerase, which would normally either prevent mismatch incorporation or promote its nucleolytic excision. Movement of a single proton on one of the mismatched bases alters the hydrogen-bonding pattern such that a base pair forms with an overall shape that is virtually indistinguishable from a canonical, Watson-Crick base pair in double-stranded DNA. These observations provide structural evidence for the rare tautomer hypothesis of spontaneous mutagenesis, a long-standing concept that has been difficult to demonstrate directly.

Modeling the enolization of succinimide derivatives, a key step of racemization of aspartic acid residues: importance of a two-H2O mechanism

Racemization of aspartic acid residues in peptides and proteins is assumed to proceed via succinimide intermediates. An enolization of the succinimide intermediate is required for the racemization to occur. In this study, we modeled the enolization step by density-functional theory (DFT) calculations (B3LYP/6-31+G**), using two model compounds, N-methylsuccinimide (1) and its formylamino derivative 2. Three mechanisms were investigated for 1, i.e., the direct mechanism without active participation of H(2)O molecules, and one-H(2)O and two-H(2)O mechanisms, in which one or two H(2)O molecules actively participate in the reaction. We found that the two-H(2)O mechanism was the most favorable with an activation barrier of 37 kcal mol(-1). In the two-H(2)O mechanism, a concerted bond reorganization involving a triple H-atom transfer occurred in an eight-membered cyclic structure formed between the imide and two H(2)O molecules. For 2, we investigated only the two-H(2)O mechanism and found that the activation barrier was lowered to 31 kcal mol(-1) due to an H-bond between the CO O-atom of the formylamino group ('the neighboring residue') and one of the H(2)O molecules. Our results suggest that, in proteins, the Asp racemization is severely controlled by the accessibility of H(2)O molecules to the reaction site of the succinimide intermediate.