Correlated ab initio quantum chemical study of the interaction of the Na+, Mg2+, Ca2+, and Zn2+ ions with the tautomers of cytosine

The interaction of the metal ions Na(+), Mg(2+), Ca(2+), and Zn(2+) with cytosine have been reinvestigated at the density functional, M?ller-Plesset, and coupled cluster levels of theory, including hitherto unstudied tautomeric forms. It has been found that the interaction of the metal ion has a varying and often significant effect on the stabilities of the various tautomers, in some cases making most stable rare tautomeric forms. The results have been analyzed with respect to method and role of ion in binding, and confirm that, as has been found for the base cytosine tautomers, B3LYP does not give energetics consistent with highly accurate post-SCF methods for their interaction with these metal ions.     

Correlated ab initio study of nucleic acid bases and their tautomers in the gas phase,in a microhydrated environment and in aqueous solution. guanine: surprising stabilization of rare tautomers in aqueous solution

Altogether eight keto and enol tautomers of guanine were studied theoretically in the gas phase, in a microhydrated environment (1 and 2 water molecules) and in bulk water. The structures of isolated, as well as mono- and dihydrated tautomers were determined by means of the RI-MP2 method using the extended TZVPP (5s3p2d1f/3s2p1d) basis set. The relative energies of isolated tautomers included the correction to higher correlation energy terms evaluated at the CCSD(T)/aug-cc-pVDZ level. The relative enthalpies at 0 K and relative free energies at 298 K were based on the above-mentioned relative energies and zero-point vibration energies, temperature-dependent enthalpy terms and entropies evaluated at the MP2/6-31G level. The keto form having hydrogen atom at N7 is the global minimum while the canonical form having hydrogen atom at N9 represents the first local minimum at all theoretical levels in vacuo and in the presence of 1 and 2 water molecules. All three unusual rare tautomers having hydrogens at N3 and N7, at N3 and N9, and also at N9 and N7 are systematically considerably less stable and can be hardly detected in the gas phase. The theoretical predictions fully agree with existing theoretical as well as experimental results. The effect of bulk solvent on the relative stability of guanine tautomers was studied by self-consistent reaction field and molecular dynamics free energy calculations using the thermodynamic integration method. Bulk solvent, surprisingly, strongly favored these three rare tautomers over all remaining low-energy tautomers and probably only these forms can exist in water phase. The global minimum (tautomer with hydrogens at N3 and N7) is by 13 kcal/mol more stable than the canonical form (3rd local minimum). Addition of one or two water molecules does not change the relative stability order of isolated guanine tautomers but the respective trend clearly supports the surprising stabilization of three rare forms.     

Specific and nonspecific metal ion-nucleotide interactions at aqueous/solid interfaces functionalized with adenine, thymine, guanine, and cytosine oligomers

This article reports nonlinear optical measurements that quantify, for the first time directly and without labels, how many Mg(2+) cations are bound to DNA 21-mers covalently linked to fused silica/water interfaces maintained at pH 7 and 10 mM NaCl, and what the thermodynamics are of these interactions. The overall interaction of Mg(2+) with adenine, thymine, guanine, and cytosine is found to involve -10.0 ± 0.3, -11.2 ± 0.3, -14.0 ± 0.4, and -14.9 ± 0.4 kJ/mol, and nonspecific interactions with the phosphate and sugar backbone are found to contribute -21.0 ± 0.6 kJ/mol for each Mg(2+) ion bound. The specific and nonspecific contributions to the interaction energy of Mg(2+) with oligonucleotide single strands is found to be additive, which suggests that within the uncertainty of these surface-specific experiments, the Mg(2+) ions are evenly distributed over the oligomers and not isolated to the most strongly binding nucleobase. The nucleobases adenine and thymine are found to bind only three Mg(2+) ions per 21-mer oligonucleotide, while the bases cytosine and guanine are found to bind eleven Mg(2+) ions per 21-mer oligonucleotide.     

Photoelectron spectroscopy of adiabatically bound valence anions of rare tautomers of the nucleic acid bases

Anionic states of nucleic acid bases (NABs) are involved in DNA damage by low-energy electrons and in charge transfer through DNA. Previous gas phase studies of free, unsolvated NAB parent anions probed mostly dipole-bound states, which are not present in condensed phase environments. Recently, we demonstrated that very rare tautomers of uracil (U), cytosine (C), adenine (A), and guanine (G), which are obtained from canonical tautomers through N-to-C proton transfers, support valence anionic states. Here we report the photoelectron spectrum of the final member of the NABs series: the valence state of the thymine (T) anion. Additionally, we summarized the work of all five NABs. All of the newfound anionic tautomers of the NABs may be formed via dissociative electron attachment followed by hydrogen atom reattachment to a carbon atom. Furthermore, these unusual tautomers may affect the structure and properties of DNA and RNA exposed to low-energy electrons. The new valence states observed here, unlike dipole bound states, could exist in condensed phases and may be relevant to radiobiological damage.     

Effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) and water coordination on the structure of glycine and zwitterionic glycine

Interactions between metal ions and amino acids are common both in solution and in the gas phase. Here, the effect of metal ions and water on the structure of glycine is examined. The effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) and water on structures of Gly.Mn+(H2O)m and GlyZwitt.Mn+(H2O)m (m = 0, 2, 5) complexes have been determined theoretically by employing the hybrid B3LYP exchange-correlation functional and using extended basis sets. Selected calculations were carried out also by means of CBS-QB3 model chemistry. The interaction enthalpies, entropies, and Gibbs energies of eight complexes Gly.Mn+ (Mn+ = Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) were determined at the B3LYP density functional level of theory. The computed Gibbs energies DeltaG degrees are negative and span a rather broad energy interval (from -90 to -1100 kJ mol(-1)), meaning that the ions studied form strong complexes. The largest interaction Gibbs energy (-1076 kJ mol(-1)) was computed for the NiGly2+ complex. Calculations of the molecular structure and relative stability of the Gly.Mn+(H2O)m and GlyZwitt.Mn+(H2O)m (Mn+ = Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+; m = 0, 2, and 5) systems indicate that in the complexes with monovalent metal cations the most stable species are the NO coordinated metal cations in non-zwitterionic glycine. Divalent cations Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+ prefer coordination via the OO bifurcated bonds of the zwitterionic glycine. Stepwise addition of two and five water molecules leads to considerable changes in the relative stability of the hydrated species. Addition of two water molecules at the metal ion in both Gly.Mn+ and GlyZwitt.Mn+ complexes reduces the relative stability of metallic complexes of glycine. For Mn+ = Li+ or Na+, the addition of five water molecules does not change the relative order of stability. In the Gly.K+ complex, the solvation shell of water molecules around K+ ion has, because of the larger size of the potassium cation, a different structure with a reduced number of hydrogen-bonded contacts. This results in a net preference (by 10.3 kJ mol(-1)) of the GlyZwitt.K+H2O5 system. Addition of five water molecules to the glycine complexes containing divalent cations Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+ results in a net preference for non-zwitterionic glycine species. The computed relative Gibbs energies are quite high (-10 to -38 kJ mol(-1)), and the NO coordination is preferred in the Gly.Mn+(H2O)5 (Mn+ = Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) complexes over the OO coordination.     

Properties of the optimal environment inducing tautomerization of complementary base pairs

The optimal environment charge configurations are predicted for the tautomerization of complementary base pairs into their corresponding rare forms, and vice versa. Results indicate that cations approaching the N₃ guanine site may induce tautomerization of the normal guanine—cytosine (G---C) base pair into its rare form. The reverse process requires that the cation approach the O₂ thymine site of the rare adenine*—thymine* pair (A*---31T*) or the O₆ guanine site of rare guanine*—cytosine* base pair (G*---C*). Possible mutagenic and antimutagenic roles of metal cations approaching base pairs are also discussed.     

Comprehensive study of the effects of methylation on tautomeric equilibria of nucleic acid bases

Minor tautomers of nucleic acid bases can result by intramolecular proton transfer. These rare tautomers could be stabilized through the addition of methyl groups to DNA bases. A comprehensive theoretical study of tautomers of methylated derivatives of guanine, adenine, cytosine, thymine, and uracil was performed. Molecular geometries of all tautomers were obtained at the density functional theory and MP2 levels with the 6-31G(d,p) basis set, and single-point calculations were performed at the CCSD(T)/6-311G(d,p) level. Tautomers obtained by protonation at the preferred protonation site for methylated isolated bases were compared to their nonmethylated counterparts. The effects of methylation on the relative stabilities of nucleic acid base tautomers are also studied and discussed in this work. The results suggest that some sites on the bases may not be mutagenic and may even stabilize the canonical Watson-Crick form. The results also indicate that a number of methylation sites can stabilize the tautomers, suggesting possible mechanisms for mutagenic changes.     

Electronic transitions of guanine tautomers, their stacked dimers, trimers and sodium complexes

Planar and nonplanar geometries of the keto-N9H and keto-N7H tautomers of the guanine base of DNA as well as the hydrogen bonded complexes of these species with three water molecules each were optimized using the density functional theory at the B3LYP/6-31G** level. Geometries of the isolated bases were also optimized using the ab initio approach at the MP2/6-31G** level. The isolated keto-N9H and keto-N7H tautomers as well as their hydrogen bonded complexes with three water molecules each were solvated in bulk water employing the polarized continuum model (PCM) of the self-consistent reaction field theory (SCRF). Stacked dimers and trimers of both the tautomers of guanine were generated by placing the planar forms of the species at a fixed distance of 3.5 A from the neighboring one and rotating one molecule with respect to the other by 110 degrees for the keto-N9H form and 90 degrees for the keto-N7H form which corresponded to total energy minima at the B3LYP/6-31G** level. Geometry optimization for the cation of the monomer of guanine was performed at the same level of theory, and its solvation in bulk water was treated using the PCM model of the SCRF theory. The geometries of complexes of the two tautomers of guanine with a Na+ ion each were optimized at the B3LYP/6-31G** level, and the Na+ ion is predicted to bind with the keto-N9H tautomer preferentially. While the complex of the keto-N7H form of guanine with three water molecules in gas phase is slightly more stable than the corresponding complex of the keto-N9H form of guanine, the reverse is true in bulk water. Stacking interactions enhance the relative stability of the keto-N9H tautomer over that of the keto-N7H tautomer, suggesting that in bulk solutions, the former would be dominant. Electronic spectra of the isolated tautomers of guanine, those of their complexes with three water molecules each, the (keto-n9h and keto-n7h) cation of guanine, the complexes of the tautomers with a Na+ ion each, the stacked dimers and trimers of the two tautomers were calculated using configuration interaction involving single electron excitations (CIS). The relative absorption intensities of the two tautomers of guanine near 275 and 248 nm in the monomer, dimer, and trimer are predicated to be in the opposite order. Thus the absorption intensity oscillation observed using a guanine aqueous solution can be explained in terms of oscillation of relative populations of the two tautomers of the molecule. The 248 nm absorption peak would be appreciably red-shifted on formation of the cation of guanine. Binding of the Na+ ion with the two tautomers of guanine reduces intensities of their transitions appreciably and also it causes large red-shifts in the same.     

Disentangling intrinsic ultrafast excited-state dynamics of cytosine tautomers

Gas-phase ultrafast excited-state dynamics of cytosine, 1-methylcytosine, and 5-fluorocytosine were investigated in molecular beams using femtosecond pump-probe photoionization spectroscopy to identify the intrinsic dynamics of the major cytosine tautomers. The results indicate that, upon photoexcitation in the first absorption band, the cytosine enol tautomer exhibits a significantly longer excited-state lifetime than its keto and imino counterparts. The initially excited states of the cytosine keto and imino tautomers decay with sub-picosecond dynamics for excitation wavelengths shorter than 300 nm, whereas that of the cytosine enol tautomer decays with time constants ranging from 3 to 45 ps for excitation between 260 and 285 nm.     

Structure and bonding of Ag(I)-DNA base complexes and Ag(I)-adenine-cytosine mispairs: an ab initio study

High level ab initio calculations have been carried out to characterize the structure, bonding and energetics of Ag(I)-DNA base complexes, including adenine or cytosine, as well as Ag(I)-adenine-cytosine mispairs. The interactions of the Ag cation in all binding sites of all adenine and cytosine tautomers have been considered. The calculations show that in gas phase the canonical form of cytosine is stabilized upon metalation, whereas the lowest energy structure of Ag-adenine correspond to a rare tautomer. Interestingly, the theoretical inspection of metalated adenine-cytosine mispair reveals that the most stable structures are formed from the canonical cytosine and adenine tautomers. The lowest energy structure is planar with adenine and cytosine hydrogen-bonded. Within few kcal/mol nonplanar, conformationally very flexible structures are found, in which the Ag(I) crosslinks an endocyclic nitrogen of adenine and the oxygen of cytosine. Metalated reverse-Wobble type of structures, on the contrary, are predicted much higher in energy.     

Effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) and water coordination on the structure and properties of L-arginine and zwitterionic L-arginine

Interactions between metal ions and amino acids are common both in solution and in the gas phase. The effect of metal ions and water on the structure of L-arginine is examined. The effects of metal ions (Li(+), Na(+), K(+), Mg(2+), Ca(2+), Ni(2+), Cu(2+), and Zn(2+)) and water on structures of Arg x M(H2O)m , m = 0, 1 complexes have been determined theoretically by employing the density functional theories (DFT) and using extended basis sets. Of the three stable complexes investigated, the relative stability of the gas-phase complexes computed with DFT methods (with the exception of K(+) systems) suggests metallic complexes of the neutral L-arginine to be the most stable species. The calculations of monohydrated systems show that even one water molecule has a profound effect on the relative stability of individual complexes. Proton dissociation enthalpies and Gibbs energies of arginine in the presence of the metal cations Li(+), Na(+), K(+), Mg(2+), Ca(2+), Ni(2+), Cu(2+), and Zn(2+) were also computed. Its gas-phase acidity considerably increases upon chelation. Of the Lewis acids investigated, the strongest affinity to arginine is exhibited by the Cu(2+) cation. The computed Gibbs energies DeltaG(o) are negative, span a rather broad energy interval (from -150 to -1500 kJ/mol), and are appreciably lowered upon hydration.     

Shape resonance states of the low-energy electron attachments to DNA base tautomers

Two different types of shape resonance states, π* and σ*, formed in the low-energy electron attachments to the low-lying tautomers of DNA bases (adenine, guanine, thymine, and cytosine) in the gas phase are investigated using the quantum scattering method with the non-empirical model potentials in a symmetry-adapted, single-center expansion. Four and three π* states are predicted for purines and pyrimidines, respectively. Comparing the different tautomers of a certain DNA base, we find distinct differences both in the resonance energies and the resonant wave function patterns of π* states. As for the lowest three π* states, the energetic values predicted in this work are also compared with the theoretical and experimental results available in the literature.     

Interactions of nucleobases with alkali earth metal cations--electrospray ionization mass spectrometric study

Interactions of nucleobases with alkali earth metal cations have been studied by electrospray ionization mass spectrometry (ESI-MS). Nucleobases containing at least one oxygen atom form stable complexes with alkali earth metal cations. This phenomenon can be explained on the grounds of the well known theory of hard and soft acids and bases. Uracil and thymine make complexes only when in their deprotonoted forms. The cations of great radii (Sr(2+), Ba(2+)) are more prone to form complexes of stoichiometry 1:1 with uracil and thymine than the cations of small radii (Mg(2+), Ca(2+)). On the other hand, Mg(2+) forms complexes of stoichiometry 2:1 and 3:2 with uracil and thymine. Gas-phase stabilities of the 1:1 complexes are higher for the cations of small radii, in contrast to the solution stabilities. For cytosine and 9- methylhypoxantine the 1:1 complexes of their deprotonated forms are observed at higher cone voltage as a result of HCl molecule loss from the complexes containing the counter ion (Cl(-)). In solution, more stable complexes are formed with metal cations of low radii. Gas-phase stability of the complexes formed by deprotonated 9- methyl-hypoxantine increases with increasing metal cation radius.     

Bond energies and attachments sites of sodium and potassium cations to DNA and RNA nucleic acid bases in the gas phase

Gas-phase metal affinities of DNA and RNA bases for the Na(+) and K(+) ions were determined at density functional level employing the hybrid B3LYP exchange correlation potential in connection with the 6-311+G(2df,2p) basis set. All the molecular complexes, obtained by the interaction between several low-lying tautomers of nucleic acid bases and the alkali ions on the different binding sites, were considered. Structural features of the sodium and potassium complexes were found to be similar except in some uracil and thymine compounds in which the tendency of potassium ion toward monocoordination appeared evident. B3LYP bond energies for both metal ions were in agreement with the available experimental results in the cases of uracil and thymine for which the most stable complex was obtained starting from the most stable tautomer of the free nucleic acid base. For adenine, although the interaction of the ions with the most stable free tautomer generated the least stable molecular complex, the best agreement with experiment was found in just this case. For the remaining cytosine and guanine bases, our calculations indicated that the metal ion affinity value closest to experiment should be determined taking into account the role played by the different tautomers of the free bases with similar energy and all the possible complexes obtained by them.     

Effect of Na+, Mg2+, and Zn2+ chlorides on the structural and thermodynamic properties of water/n-heptane interfaces

The effect on the structural and thermodynamic properties in water/n-heptane interfaces on addition of NaCl, MgCl(2), and ZnCl(2) has been examined through five independent 100-ns molecular dynamics simulations. Results indicate that the interfacial thickness within the framework of the capillary-wave model decreases on addition of electrolytes in the order Na(+) < Mg(2+) < Zn(2+), whereas the interfacial tension increases in the same order. Ionic density profiles and self-diffusion coefficients are strongly influenced by the strength of the first hydration shell, which varies in the order Na(+) < Mg(2+) < Zn(2+). On the other hand, the Cl(-) behavior, that is, diffusion and solvation sphere, is influenced by its counterion. Accordingly, cations are strongly expelled from the interface, which is especially remarkable for the small divalent cations. This fact alters the water geometry near the interface and in a lesser extent n-heptane order and number of hydrogen bonds per water molecule close to the interface.     

带电组氨酸侧链与DNA碱基间非键作用强度的理论研究

采用MP2方法和6-31+G(d,p)基组优化得到了带有一个正电荷的组氨酸侧链与4个DNA碱基间形成的18个氢键复合物的气相稳定结构, 从文献中获取了组氨酸侧链与DNA碱基间形成的12个堆积和T型复合物的气相稳定结构, 使用包含基组重叠误差(BSSE)校正的MP2方法和aug-cc-pVTZ基组及密度泛函理论M06-2X-D3方法和aug-cc-pVDZ基组计算了这些复合物的结合能. 研究结果表明, 包含BSSE校正的M06-2X-D3方法和aug-cc-pVDZ基组能够给出较准确的结合能; 气相条件下, 组氨酸侧链与同种DNA碱基间的离子氢键作用明显强于堆积作用和T型作用, 组氨酸侧链最易通过离子氢键与胞嘧啶C和鸟嘌呤G作用形成氢键复合物, 组氨酸与胞嘧啶C和鸟嘌呤G间的T型作用强于与腺嘌呤A和胸腺嘧啶T间的离子氢键作用; 水相条件下, 组氨酸侧链与同种DNA碱基间的离子氢键作用仍明显强于堆积作用和T型作用, 组氨酸侧链更易与胞嘧啶C和鸟嘌呤G相互作用形成氢键复合物, 但是最强的组氨酸侧链与胞嘧啶C间的T型作用明显弱于与腺嘌呤A和胸腺嘧啶T间的离子氢键作用, 说明水相条件下组氨酸侧链与DNA碱基间主要通过离子氢键作用形成氢键复合物.     

金属离子与腺嘌呤异构体复合物的稳定性

在B3LYP/6-311++G**//6-31+G*基组水平上结合PCM模型系统的优化了气相和液相环境中14种腺嘌呤异构体与四种金属离子(Na+,K+,Mg2+,Ca2+)形成的132个稳定复合物。通过能量对比,得到了所有复合物异构体在气液两相中的稳定性顺序及不同金属离子与同一异构体形成的复合物的能量排序,首次给出了液相中这些不同金属离子复合物的最稳定结构。结果发现溶剂效应导致了液相中的复合物稳定性顺序与气相中的相比发生了很大变化;同一异构体与不同金属离子形成的复合物其稳定性在排序中也变化很大。对于这些变化,本文分别从金属离子与腺嘌呤的复合物在气相中的结合能(EBE)及在液相中的溶质溶剂相互作用能(Epol)能等方面进行了系统的阐述     

Ab initio determination of the ionization potentials of DNA and RNA nucleobases

Quantum chemical high level ab initio coupled-cluster and multiconfigurational perturbation methods have been used to compute vertical and adiabatic ionization potentials of the five canonical DNA and RNA nucleobases: uracil, thymine, cytosine, adenine, and guanine. Several states of their cations have been also calculated. The present results represent a systematic compendium of these magnitudes, establishing theoretical reference values at a level not reported before, calibrating computational strategies, and guiding the assignment of the features in the experimental photoelectron spectra.     

Electronic structure and spectroscopy of nucleic acid bases: ionization energies, ionization-induced structural changes, and photoelectron spectra

We report high-level ab initio calculations and single-photon ionization mass spectrometry study of ionization of adenine (A), thymine (T), cytosine (C), and guanine (G). For thymine and adenine, only the lowest-energy tautomers were considered, whereas for cytosine and guanine we characterized the five lowest-energy tautomeric forms. The first adiabatic and several vertical ionization energies were computed using the equation-of-motion coupled-cluster method for ionization potentials with single and double substitutions. Equilibrium structures of the cationic ground states were characterized by DFT with the ωB97X-D functional. The ionization-induced geometry changes of the bases are consistent with the shapes of the corresponding molecular orbitals. For the lowest-energy tautomers, the magnitude of the structural relaxation decreases in the following series, G > C > A > T, the respective relaxation energies being 0.41, 0.32, 0.25, and 0.20 eV. The computed adiabatic ionization energies (8.13, 8.89, 8.51-8.67, and 7.75-7.87 eV for A, T, C, and G, respectively) agree well with the onsets of the photoionization efficiency (PIE) curves (8.20 ± 0.05, 8.95 ± 0.05, 8.60 ± 0.05, and 7.75 ± 0.05 eV). Vibrational progressions for the S(0)-D(0) vibronic bands computed within double-harmonic approximation with Duschinsky rotations are compared with previously reported experimental photoelectron spectra and differentiated PIE curves.     

Influence of the environment on the specificity of the Mg(II) binding to uracil

Interactions of uracil with the Mg(2+) ion were studied theoretically in the gas phase and in solution. The bare Mg(2+) prefers bidentate N-C═O binding sites stabilizing rare keto-enol forms of the base. Hydration and/or phosphate binding of the Mg(2+) ion shield its positive charge, which leads to preference of monodentate binding to the oxygen keto atoms, shifting fully the equilibrium between the tautomers back toward the canonical diketo tautomer. In solution, a direct inner-sphere metal binding to uracil is not clearly advantageous compared to the outer-sphere metal binding. Similar trends were also obtained for the Ca(2+) ion. Results are supported by the natural bond orbital (NBO) and atoms in molecule (AIM) analyses and the combined extended transition-state energy decomposition analysis and natural orbitals for chemical valence (ETS-NOCV).