Hydration and stability of nucleic acid bases and base pairs

Microhydration of protonated nucleic acid bases and protonated nucleosides in the gas phase

Thermochemical data, DeltaH(o)(n), DeltaS(o)(n), and DeltaG(o)(n), for the hydration of protonated nucleic acid bases and protonated nucleosides have been experimentally studied by equilibrium measurements using an electrospray high-pressure mass spectrometer equipped with a pulsed ion-beam reaction chamber. For protonated nucleobases the hydration enthalpies were found to be similar for all studied systems and varied between 12.4-13.1 kcal/mol for the first and 11.2-11.5 kcal/mol for the second water molecule. While for protonated nucleosides the water binding enthalpies (11.7-13.3 kcal/mol) are very close to those for protonated nucleobases, the entropy values are "more negative." The structural and energetic aspects of hydrated ions are discussed in conjunction with the available theoretical data.

Characterizing radiation-induced oxidation of DNA by way of the monohydrated guanine-cytosine radical cation

The interaction of one water molecule with the guanine-cytosine radical cation has been studied with ab initio and density functional methods in order to help elucidate the nature of oxidized aqueous DNA. The theoretical spin density of [GC]*(+) reveals that the radical center is localized on guanine. The adiabatic ionization potential lowers from 7.63 to 6.71 eV in concurrence with the formation of the Watson-Crick base pair and hydration by one water molecule. A natural bond orbital analysis of partial charges shows that approximately 80% of the positive charge persists on guanine upon hydration and formation of the Watson-Crick base pair with cytosine. Hydration energies were computed with second-order Z-averaged perturbation theory (ZAPT2) using the aug-cc-pVDZ basis set at 11 stationary points on the B3LYP/DZP++ potential energy surface. The hydration energy at the global minimum is 14.2 kcal mol(-1). The lowest energy structures correspond to hydration near the glycosidic bond sites. Structural changes in the Watson-Crick base pair are predominantly seen for monohydration in the groove regions of double-helix DNA.

Microhydration of guanine...cytosine base pairs, a theoretical Study on the role of water in stability, structure and tautomeric equilibrium

The potential energy surfaces of guanine...cytosine complexes and microhydrated guanine...cytosine (one and two water molecules) were investigated by the molecular dynamics/quenching method (MD/Q), using the empirical potential Parm94 force field, implemented in the Amber program package. The calculations were conducted for all the possible combinations of the four most stable tautomers of guanine and three of cytosine (covering the canonical forms in both cases). The obtained structures were sorted by their structural motifs into three main groups: planar hydrogen-bonded; stacked; and T-shaped structures. The most stable structures found at the empirical potential energy surfaces were fully reoptimised at the second-order Møller-Plesset perturbation theory as well as using the density functional method with an empirical dispersion term (DFT-D). A combination of the canonical form of guanine and cytosine and canonical cytosine with a guanine tautomer where the hydrogen is switched from position N9 to N7 are energetically preferred in microsolvated systems as well as those without the presence of a solvent. The rising number of water molecules leads to smaller differences between the stability of the various combinations of the tautomers of bases in the base pairs. For some of the tautomer combinations (mainly the enol-enol combination), two water molecules are sufficient for the preference of stacked structures over the H-bonded ones. The interaction energies and geometries obtained by the second-order Møller-Plesset perturbation theory method and the much less computationally demanding DFT-D method are comparable, except for stacked complexes, where the interaction energies are overestimated on average by 3 kcal mol(-1) at the MP2 level.

Nuclear magnetic resonance investigation of primer--template models: formation of a pyrimidine bulge upon misincorporation

Our previous studies have shown that misaligned structures can occur upon misincorporation of a dNTP opposite thymine templates. The formation of misaligned structures during DNA replication, if not repaired properly, can be bypassed and extended by low-fidelity polymerases and ultimately lead to mutations. In this study, the base pair structures at the replicating sites of a set of primer-template models which mimic the situation upon misincorporation of a dNTP opposite cytosine templates have been determined. High-resolution NMR structural results show that misaligned structures with a C-bulge can be formed upon incorporation of dCTP, dTTP, and dATP opposite 5'-GC, 5'-AC, and 5'-TC templates, respectively. The stabilities of misaligned structures depend on the types of terminal base pairs at the replicating sites. Together with the structural findings in thymine templates, we conclude that terminal G.C and C.G base pairs always contribute a larger stabilizing effect to the misaligned structures containing a pyrimidine bulge than terminal A.T and T.A base pairs. Misalignment and thus deletion mutation are more likely to occur if misincorporation of a nucleotide opposite a pyrimidine template can cause template slippage to form a terminal G.C or C.G base pair. Although misalignment also occurs when the newly formed terminal base pair is an A.T base pair or a T.A base pair, both misaligned and mismatched conformers coexist, which can lead to deletion and substitution mutations, respectively.

Determination of Redox Potentials for the Watson−Crick Base Pairs, DNA Nucleosides, and Relevant Nucleoside Analogues

Redox potentials for the DNA nucleobases and nucleosides, various relevant nucleoside analogues, Watson-Crick base pairs, and seven organic dyes are presented based on DFT/B3LYP/6-31++G(d,p) and B3YLP/6-311+G(2df,p)//B3LYP/6-31+G* levels of calculations. The values are determined from an experimentally calibrated set of equations that correlate the vertical ionization (electron affinity) energy of 20 organic molecules with their experimental reversible oxidation (reduction) potential. Our results are in good agreement with those estimated experimentally for the DNA nucleosides in acetonitrile solutions (Seidel et al. J. Phys. Chem. 1996, 100, 5541). We have found that nucleosides with anti conformation exhibit lower oxidation potentials than the corresponding syn conformers. The lowering in the oxidation potential is due to the formation of an intramolecular hydrogen bonding interaction between the 5'-OH group of the sugar and the N3 of the purine bases or C2=O of the pyrimidine bases in the syn conformation. Pairing of adenine or guanine with its complementary pyrimidine base decreases its oxidation potential by 0.15 or 0.28 V, respectively. The calculated energy difference between the oxidation potential for the G.C base pair and that of the guanine base is in good agreement with the experimental value estimated recently (0.34 V: Caruso, T.; et al. J. Am. Chem. Soc. 2005, 127, 15040). The complete and consistent set of reversible redox values determined in this work for the DNA constituents is expected to be of considerable value to those studying charge and electronic energy transfer in DNA.