Effects of Microsolvation on the Adenine−Uracil Base Pair and Its Radical Anion: Adenine−Uracil Mono- and Dihydrates

Intramolecular vibrational energy redistribution in nucleobases: Excitation of NH stretching vibrations in adenine–uracil + H2O

Redistribution of vibrational energy in the adenine–uracil base pair is studied when the base pair undergoes an intermolecular interaction with an overtone-bending vibration excited H2O(2[Formula: see text]bend) molecule. Energy transfer is calculated using the structural information obtained from density functional theory in the solution of the equations of motion. Intermolecular vibrational energy transfer (VET) from H2O(2[Formula: see text]bend) to the uracil–NH stretching mode is efficient and rapidly followed by intramolecular vibrational energy redistribution (IVR) resulting from coupling between vibrational modes. An important pathway is IVR carrying energy to the NH-stretching mode of the adenine moiety in a subpicosecond scale, the energy build-up being sigmoidal, when H2O interacts with the uracil–NH bond. The majority of intermolecular hydrogen bonds between the base pair and H2O are weakened but unbroken during the ultrafast energy redistribution period. Lifetimes of intermolecular HB are on the order of 0.5 ps. The efficiency of IVR in the base pair is due to near-resonance between coupled CC and CN vibrations. The resonance also exists between the frequencies of H2O bend and NH stretch, thus facilitating VET. When H2O interacts with the NH bond at the adenine end of the base pair, energy flow in the reverse direction to the uracil–NH stretch is negligible, the unidirectionality discussed in terms of the effects of uracil CH stretches. The energy distributed in the CH bonds is found to be significant. The IVR process is found to be nearly temperature independent between 200 and 400 K.

Excited-state calculations with TD-DFT: from benchmarks to simulations in complex environments

In this perspective, we present an overview of recent progress on Time-Dependent Density Functional Theory (TD-DFT) with a specific focus on its accuracy and on models able to take into account environmental effects, including complex media. To this end, we first summarise recent benchmarks and define an average TD-DFT accuracy in reproducing excitation energies when a conventional approach is used. Next, coupling of TD-DFT with models able to account for different kinds of interactions between a central chromophore and nearby chemical objects (solvent, organic cage, metal as well as semi-conducting surface) is investigated. Examples of application to excitation properties are presented, allowing to briefly describe several recent computational strategies. In addition, an extension of TD-DFT to describe a phenomenon involving interacting chromophores, e.g. the electronic energy transfer (EET), is presented to illustrate that this methodology can be applied to processes beyond the vertical excitation. This perspective therefore aims to provide to non-specialists a flavour of recent trends in the field of simulations of excited states in "realistic" situations.

DFT study on the radical anions formed by primaquine and its derivatives

The electron affinities (EA) of the 8-aminoquinoline antimalarial drug primaquine and several of its metabolites were studied using the density functional theory method. We first considered six substituents at the 5-position, -CH(3), -OH, -OCH(3), -Ph, -OPh, and -CHO. We found that in the gas phase the adiabatic EAs are similar to that of the parent primaquine for the -CH(3), -OH, and -OCH(3) substituents. In contrast, the -Ph, -OPh, and -CHO substituents all markedly increase the adiabatic EA. However, only the -CHO substituted compound is predicted to form a stable covalently bound radical anion in the gas phase due to its significant positive vertical EA relative to that of the parent primaquine. In addition, when the 8-position is substituted by the N-hydroxyl group or a quinone-imine structure is formed, the electron capture ability is significantly increased. In aqueous solution, all these molecules have significantly larger adiabatic EAs than in the gas phase. In addition, all of the vertical EAs are positive in aqueous solution. The implications of these findings for contributing to our mechanistic understanding of the red cell toxicity of 8-aminoquinoline compounds are further discussed.

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.

Excess electron localization in solvated DNA bases

We present a first-principles molecular dynamics study of an excess electron in condensed phase models of solvated DNA bases. Calculations on increasingly large microsolvated clusters taken from liquid phase simulations show that adiabatic electron affinities increase systematically upon solvation, as for optimized gas-phase geometries. Dynamical simulations after vertical attachment indicate that the excess electron, which is initially found delocalized, localizes around the nucleobases within a 15 fs time scale. This transition requires small rearrangements in the geometry of the bases.

Vertical ionization potentials of nucleobases in a fully solvated DNA environment

Vertical ionization potentials (IPs) of nucleobases embedded in a fully solvated DNA fragment (12-mer B-DNA fragment + 22 sodium counterions + 5760 water molecules equilibrated to 298 K) have been calculated using a combined quantum mechanical molecular mechanics (QM/MM) approach. Calculations of the vertical IP of the anion Cl(-) are reported that support the accuracy of the application of a QM/MM method to this problem. It is shown that the pi nucleotide HOMO origin for the emitted electron is localized on the base by the hydration structure surrounding the DNA in a way similar to that recently observed for pyrimidine nucleotides in aqueous solutions (Slavicek, P.; et al. J. Am. Chem. Soc. 2009, 131, 6460). In a first step, a high level of theory, CCSD(T)/aug-cc-pVDZ, was used to calculate the vertical IP of each of the four single bases isolated in the QM region while the remaining DNA fragment, counterions, and water solvent molecules were included in the MM region. The calculated vertical IPs show a large positive shift of 3.2-3.3 eV compared to the corresponding gas-phase values. This shift is similar for all four DNA bases. The origin of the large increase in vertical IPs of nucleobases is found to be the long-range electrostatic interactions with the solvation structure outside the DNA helix. Thermal fluctuations in the fluid can result in IP changes of roughly 1 eV on a picosecond time scale. IPs of pi-stacked and H-bonded clusters of DNA bases were also calculated using the same QM/MM model but with a lower level of theory, B3LYP/6-31G(d=0.2). An IP shift of 4.02 eV relative to the gas phase is found for a four-base-pair B-DNA duplex configuration. The primary goal of this work was to estimate the influence of long-range solvation interactions on the ionization properties of DNA bases rather than provide highly precise IP evaluations. The QM/MM model presented in this work provides an attractive method to treat the difficult problem of incorporating a detailed long-range structural model of physiological conditions into investigations of the electronic processes in DNA.

Elementary lesions in DNA subunits: electron, hydrogen atom, proton, and hydride transfers

When DNA is damaged by ionizing radiation, the genes in a cell may acquire mutations or the cell could die. The smallest known DNA-damaging unit is an electron, often low-energy secondary electrons. Additional electrons and transfers involving hydrogen atoms, protons, and hydride anions can damage DNA subunits, including individual nucleobases and nucleoside pairs. Researchers would like to better understand the molecular mechanisms involved in DNA damage from ionizing radiation. In this Account, we highlight our theoretical investigations of the molecular mechanisms of DNA damage using quantum mechanical models. Our investigations use robust theoretical methods with computations conducted in the gas phase and with solution models. We calculate adiabatic electron affinities (AEAs), which describe the energetics of electronic attachment to closed-shell DNA subunits, for the free bases, nucleosides, nucleotides, base pairs, and single and double DNA strand units. Electron affinities for free nucleobases yield the order uracil > thymine > cytosine > guanine > adenine and the same order for the DNA nucleosides, mononucleotides, and nucleoside 3',5'-diphosphates. AEA values increase steadily with the size and complexity of the system because of greater hydration, glycosylation, nucleotide formation, and base pairing. We predict and experimental results partially confirm that most of the more complex and hydrated species are observable as radical anions. Our modeling studies indicate that depyrimidination reactions of radical anion nucleosides release cytosine more often than thymine. Recent experimental results support those findings. Our theoretical studies of DNA base-pair radical anions predict increases in electron affinity accompanying H bonding and solvation. Electron addition facilitates some proton transfers in these pairs, which results in strongly perturbed pairing configurations. Of all nucleobase moieties within the more complex radical anion systems, thymine is best able to retain a negative charge. Charge and spin are well-separated in some of these systems. Radical species derived via hydrogen abstraction from DNA subunits yield large AEA values because they form closed-shell anions. Our studies predict single-strand breaks following H abstraction from nucleotides. Some H-abstraction processes in the DNA base pairs lead to severe distortions in pairing configuration based on our calculations. This body of systematic theoretical studies provides realistic descriptions of some events that lead to elementary DNA lesions, while providing rationalizations for many observed phenomena. Such approaches can lead to the design of new experiments, which would contribute to our understanding of the chemical physics of nucleic acids.