Vertical ionization potentials of nucleobases in a fully solvated DNA environment

Developing a computational model that accurately reproduces the structural features of a dinucleoside monophosphate unit within B-DNA

The ability of a dinucleoside monophosphate to mimic the conformation of B-DNA was tested using a combination of different phosphate models (anionic, neutral, counterion), environments (gas, water), and density functionals (B3LYP, MPWB1K, M06-2X) with the 6-31G(d,p) basis set. Three sequences (5'-GX(Py)-3', where X(Py) = T, U or (Br)U) were considered, which vary in the (natural or modified) 3' pyrimidine nucleobase (X(Py)). These bases were selected due to their presence in natural DNA, structural similarity to T and/or applications in radical-initiated anti-tumour therapies. The accuracy of each of the 54 model, method and sequence combinations was judged based on the ability to reproduce key structural features of natural B-DNA, including the stacked base-base orientation and important backbone torsion angles. B3LYP yields distorted or tilted relative base-base orientations, while many computational challenges were encountered for MPWB1K. Despite wide use in computational studies of DNA, the neutral (protonated) phosphate model could not consistently predict a stacked arrangement for all sequences. In contrast, stacked base-base arrangements were obtained for all sequences with M06-2X in conjunction with both the anionic and (sodium) counterion phosphate models. However, comparison of calculated and experimental backbone conformations reveals the charge-neutralized counterion phosphate model best mimics B-DNA. Structures optimized with implicit solvent (water) are comparable to gas-phase structures, which suggests similar results should be obtained in an environment of intermediate polarity. We recommend the use of either gas-phase or water M06-2X optimizations with the counterion phosphate model to study the structure and/or reactivity of natural or modified DNA with a dinucleoside monophosphate.

Effect of solvation on the vertical ionization energy of thymine: from microhydration to bulk

The effect of hydration on the vertical ionization energy (VIE) of thymine was characterized using equation-of-motion ionization potential coupled-cluster (EOM-IP-CCSD) and effective fragment potential (EFP) methods. We considered several microsolvated clusters as well as thymine solvated in bulk water. The VIE in bulk water was computed by averaging over solvent-solute configurations obtained from equilibrium molecular dynamics trajectories at 300 K. The effect of microsolvation was analyzed and contrasted against the combined effect of the first solvation shell in bulk water. Microsolvation reduces the ionization energy (IE) by about 0.1 eV per water molecule, while the first solvation shell increases the IE by 0.1 eV. The subsequent solvation lowers the IE, and the bulk value of the solvent-induced shift of thymine's VIE is approximately -0.9 eV. The combined effect of the first solvation shell was explained in terms of specific solute-solvent interactions, which were investigated using model structures. The convergence of IE to the bulk value requires the hydration sphere of approximately 13.5 Å radius. The performance of the EOM-IP-CCSD/EFP scheme was benchmarked against full EOM-IP-CCSD using microhydrated structures. The errors were found to be less than 0.01-0.02 eV. The relative importance of the polarization and higher multipole moments in EFP model was also investigated.

Theoretical study of the tautomerism in the one-electron oxidized guanine-cytosine base pair

Ionizing radiation on DNA mainly generates one-electron oxidized guanine-cytosine base pair (G(+·):C), and in the present paper we study all possible tautomers of G(+·):C by using ab initio approaches. Our calculations reveal that the tautomeric equilibrium follows a peculiar path, characterized by a stepwise mechanism: first the proton in the central hydrogen bond N1(G)-H1-N3(C) migrates from guanine to cytosine, and then the cytosine cation releases one proton from its amino group. During this second step, water acts as a proton acceptor, localizing the positive charge on one of the water molecules interacting with the guanine radical. In agreement with experimental findings, the computed energy barriers show that the deprotonation of the cytosine cation is the speed-limiting step in the tautomeric equilibrium. The influence of the number of water molecules incorporated in the theoretical model is analyzed in detail. The evolution of electronic properties along the reaction path is also discussed on the basis of partial atomic charges and spin density distributions. This work demonstrates that water indeed plays a crucial role in the tautomeric equilibra of base pairs.

Ionization of purine tautomers in nucleobases, nucleosides, and nucleotides: from the gas phase to the aqueous environment

We have simulated ionization of purine nucleic acid components in the gas phase and in a water environment. The vertical and adiabatic ionization processes were calculated at the PMP2/aug-cc-pVDZ level with the TDDFT method applied to obtain ionization from the deeper lying orbitals. The water environment was modeled via microsolvation approach and using a nonequilibrium polarizable continuum model. We have characterized a set of guanine tautomers and investigated nucleosides and nucleotides in different conformations. The results for guanine, i.e., the nucleic acid base with the lowest vertical ionization potential, were also compared to those for the other purine base, adenine. The main findings of our study are the following: (i) Guanine remains clearly the base with the lowest ionization energy even upon aqueous solvation. (ii) Water solvent has a strong effect on the ionization energetics of guanine and adenine and their derivatives; the vertical ionization potential (VIP) is lowered by about 1 eV for guanine while it is ∼1.5 eV higher in the nucleotides, overall resulting in similar VIPs for GMP(-), guanosine and guanine in water. (iii) Water efficiently screens the electrostatic interactions between nucleic acid components. Consequently, ionization in water always originates from the base unit of the nucleic acid and all the information about conformational state is lost in the ionization energetics. (iv) The energy splitting between ionization of the two least bound electrons increases upon solvation. (v) Tautomerism does not contribute to the width of the photoelectron spectra in water. (vi) The effect of specific short-range interactions with individual solvent molecules is negligible for purine bases, compared to the long-range dielectric effects of the aqueous medium.