Microhydration of the guanine-cytosine (GC) base pair in the neutral and anionic radical states: a density functional study

Secondary Electron Attachment-Induced Radiation Damage to Genetic Materials

Reactions of radiation-produced secondary electrons (SEs) with biomacromolecules (e.g., DNA) are considered one of the primary causes of radiation-induced cell death. In this Review, we summarize the latest developments in the modeling of SE attachment-induced radiation damage. The initial attachment of electrons to genetic materials has traditionally been attributed to the temporary bound or resonance states. Recent studies have, however, indicated an alternative possibility with two steps. First, the dipole-bound states act as a doorway for electron capture. Subsequently, the electron gets transferred to the valence-bound state, in which the electron is localized on the nucleobase. The transfer from the dipole-bound to valence-bound state happens through a mixing of electronic and nuclear degrees of freedom. In the presence of aqueous media, the water-bound states act as the doorway state, which is similar to that of the presolvated electron. Electron transfer from the initial doorway state to the nucleobase-bound state in the presence of bulk aqueous media happens on an ultrafast time scale, and it can account for the decrease in DNA strand breaks in aqueous environments. Analyses of the theoretically obtained results along with experimental data have also been discussed.

Attosecond charge migration following oxygen K-shell ionization in DNA bases and base pairs

Core ionization of DNA begins a cascade of events which could lead to cellular inactivation or death. The created core-hole following an impulse inner-shell ionization of molecules naturally decays in the auger timescale. We simulated charge migration (CM) phenomena following an impulsive core ionization of individual DNA bases at the oxygen K-edge which occurs before Auger decay of the oxygen. Our approach is based on real-time time dependent density functional theory (RT-TDDFT). It is shown that the pronounced hole fluctuation observed around bonds of the initial core-hole results in various valence orbital migrations. Also, the same photo-core-ionized dynamics is studied for the related base pairs. We investigate the role of base pairing and H-bonding interactions in the attosecond CM dynamics. In particular, the creation of a core-hole in the oxygen involved in H-bonding leads to an enhancement of charge migration relative to the respective single bases. Importantly, the hole oscillation of the adenine-thymine base pair upon creation of a core-hole at the oxygen, which does not contribute to the donor-acceptor interactions (not H-bonded), decreases compared to the single thymine base. Understanding the detailed dynamics of the localized core-hole initiating CM process would open the way for chemically controlling DNA damage/repair in the future.

Electron Attachment to Cytosine: The Role of Water

We present an EOM-CCSD-based quantum mechanical/molecular mechanical (QM/MM) study on the electron attachment process to solvated cytosine. The electron attachment in the bulk solvated cytosine occurs through a doorway mechanism, where the initial electron is localized on water. The electron is subsequently transferred to cytosine by the mixing of electronic and nuclear degrees of freedom, which occurs on an ultrafast time scale. The bulk water environment stabilizes the cytosine-bound anion by an extensive hydrogen-bond network and drastically enhances the electron transfer rate from that observed in the gas phase. Microhydration studies cannot reproduce the effect of the bulk water environment on the electron attachment process, and one needs to include a large number of water molecules in the calculation to obtain converged results. The predicted adiabatic electron affinity and electron transfer rate obtained from our QM/MM calculations are consistent with the available experimental results.

Electron-Induced Repair of 2'-Deoxyribose Sugar Radicals in DNA: A Density Functional Theory (DFT) Study

In this work, we used ωB97XD density functional and 6-31++G** basis set to study the structure, electron affinity, populations via Boltzmann distribution, and one-electron reduction potentials (E°) of 2′-deoxyribose sugar radicals in aqueous phase by considering 2′-deoxyguanosine and 2′-deoxythymidine as a model of DNA. The calculation predicted the relative stability of sugar radicals in the order C4′^• > C1′^• > C5′^• > C3′^• > C2′^•. The Boltzmann distribution populations based on the relative stability of the sugar radicals were not those found for ionizing radiation or OH-radical attack and are good evidence the kinetic mechanisms of the processes drive the products formed. The adiabatic electron affinities of these sugar radicals were in the range 2.6–3.3 eV which is higher than the canonical DNA bases. The sugar radicals reduction potentials (E°) without protonation (−1.8 to −1.2 V) were also significantly higher than the bases. Thus the sugar radicals will be far more readily reduced by solvated electrons than the DNA bases. In the aqueous phase, these one-electron reduced sugar radicals (anions) are protonated from solvent and thus are efficiently repaired via the “electron-induced proton transfer mechanism”. The calculation shows that, in comparison to efficient repair of sugar radicals by the electron-induced proton transfer mechanism, the repair of the cyclopurine lesion, 5′,8-cyclo-2′-dG, would involve a substantial barrier.

Effects of Positive and Negative Ionization on Prototropy in Pyrimidine Bases: An Unusual Case of Isocytosine

Intramolecular proton-transfers (prototropic conversions) have been studied for the guanine building block isocytosine (iC), and effects of positive ionization, called one-electron oxidation (iC - e → iC^+•), and negative ionization, called one-electron reduction (iC + e → iC^-•), on tautomeric conversions when proceeding from neutral to ionized isocytosine have been discussed. Although radical cations and radical anions are very short-lived species, the ionization effects could be investigated by quantum-chemical methods. Such kind of studies gives some information about the labile protons and the most basic positions in the neutral and radical forms of the tautomeric system. For investigations, the complete isomeric mixture of isocytosine has been considered and calculations performed in two extreme environments, apolar {DFT(B3LYP)/6-311+G(d,p)} and polar {PCM(water)//DFT(B3LYP)/6-311+G(d,p)}. For selected isomers, the G4 theory has also been applied. There are no good relations for energetic parameters of neutral and ionized forms. Ionization energies depend on localization of labile protons. Tautomeric equilibria for neutral and ionized isocytosine, favored sites of protonation and deprotonation, and favored structures of protonated and deprotonated forms strongly depend on environment. Acidity of iC^+• is close to that of the iC conjugate acid, and basicity of iC^-• is close to that of the iC conjugate base. This increase of acid-base properties of charged radicals explains the proton-transfer in ionized pairs of nucleobases. When compared to other pyrimidine bases such as uracil (U) and cytosine (C), which exhibit analogous tautomeric equilibria between nine prototropic tautomers as isocytosine, the tautomeric preferences for iC, iC^+•, iC^-•, U, U^+•, U^-•, C, C^+•, and C^-• are completely different. The differences suggest that acid-base properties of functional groups, their stabilities, and ionization energies play a principal role in proton-transfers for pyrimidine bases and influence compositions of tautomeric mixtures.

Microsolvation within the Systematic Molecular Fragmentation by Annihilation Approach

We have applied the systematic molecular fragmentation by annihilation (SMFA) fragmentation technique to glycine and DNA base pairs in water clusters, systems for which explicit solvation is believed to be important. The SMFA method was found to be capable of describing the structures, especially in handling the complexity of hydrogen bonding, with energies produced being comparable with those from full molecule results. Thus, the ability to break down large calculations into a manageable time without loss of accuracy shows promise for application to real biological systems for which these effects are relevant.

Electron interaction with a DNA duplex: dCpdC:dGpdG

Electron attachment to double-stranded cytosine-rich DNA, dCpdC:dGpdG, has been studied by density functional theory. This system represents a minimal descriptive unit of a cytosine-rich double-stranded DNA helix. A significant electron affinity for the formation of a cytosine-centered radical anion is revealed to be about 2.2 eV. The excess electron may reside on the nucleobase at the 5' position (dC˙(-)pdC:dGpdG) or at the 3' position (dCpdC˙(-):dGpdG). The inter-strand proton transfer between the radical anion centered cytosine (N3) and the paired guanine (HN1) results in the formation of radical anion center separated complexes dC1H˙pdC:dG2-H(-)pdG and dCpdC2H˙:dGpdG1-H(-). These distonic radical anions are found to be approximately 1 to 4 kcal mol(-1) more stable than the normal radical anions. Intra-strand cytosine π→π transition energies are below the electron detachment energy. Inter-strand π→π transitions of the excess electron from C to G are predicted to be less than 2.79 eV. Electron transfer might also be possible through the inter-strand base-jumping mode. An analysis of absorption visible spectra reveals the absorption bands ranging from 500 nm to 700 nm for the cytosine-rich radical anions of the DNA duplex. Electron attachment to cytidine oligomers might add color to the DNA duplex.

The influence of phosphorothioate on charge migration in single and double stranded DNA: a theoretical approach

In this study the influence of the phosphorothioate internucleotide bond on the electronic properties of single and double-stranded short nucleotides has been investigated at the M06-2X/6-31+G** level of theory in the gaseous phase. Due to the chirality of the phosphorus atom in a phosphorothioate (PT) internucleotide diester bond, the adiabatic/vertical mode of electron affinity/ionization potential, spin density and molecular orbital distribution, as well as structural analysis were taken under consideration for the single stranded (ss) R(P) and S(P) diastereomers of d[G(PS)G] and for double stranded (ds) d[G(PS)G]*d[C(PO)C], in comparison with the corresponding parent phosphate compounds. Moreover, the excitation states, HOMO and LUMO energies were calculated using a TD-DFT methodology at the M06-2X/6-31+G**//M06-2X/6-31++G** level of theory in the aqueous phase. The obtained results show that the PT plays a significant role in the case of ss-oligonucleotides, and to a much smaller extent in ds-oligomers.

New size-expanded RNA nucleobase analogs: a detailed theoretical study

Fluorescent nucleobase analogs have attracted much attention in recent years due to their potential applications in nucleic acids research. In this work, four new size-expanded RNA base analogs were computationally designed and their structural, electronic, and optical properties are investigated by means of DFT calculations. The results indicate that these analogs can form stable Watson-Crick base pairs with natural counterparts and they have smaller ionization potentials and HOMO-LUMO gaps than natural ones. Particularly, the electronic absorption spectra and fluorescent emission spectra are calculated. The calculated excitation maxima are greatly red-shifted compared with their parental and natural bases, allowing them to be selectively excited. In gas phase, the fluorescence from them would be expected to occur around 526, 489, 510, and 462 nm, respectively. The influences of water solution and base pairing on the relevant absorption spectra of these base analogs are also examined.

Mutagenic effects induced by the attack of NO2 radical to the guanine-cytosine base pair

We investigate the attack of the nitrogen dioxide radical (NO^•₂) to the guanine—cytosine (GC) base pair and the subsequent tautomeric reactions able to induce mutations, by means of density functional theory (DFT) calculations. The conducted simulations allow us to identify the most reactive sites of the GC base pair. Indeed, the computed relative energies demonstrate that the addition of the NO^•₂ radical to the C8 position of the guanine base forms to the most stable adduct. Although the initial adducts might evolve to non-canonical structures via inter-base hydrogen bonds rearrangements, the probability for the proton exchange to occur lies in the same range as that observed for undamaged DNA. As a result, tautomeric errors in NO₂-attacked DNA arises at the same rate as in canonical DNA, with no macroscopic impact on the overall stability of DNA. The potential mutagenic effects of the GC–NO^•₂ radical adducts likely involve side reactions, e.g., the GC deprotonation to the solvent, rather than proton exchange between guanine and cytosine basis.

DNA spontaneous mutation and its role in the evolution of GC-content: assessing the impact of the genetic sequence

We use theoretical tools to investigate the possible role played by a DNA sequence in the base pair tautomerization phenomena.

Hydrogen-bonded proton transfer in the hydrated adenine–thymine anion

The proton transfer processes of microhydrated adenine–thymine anions are studied using density functional theory with the B3LYP method and DZP++ basis set. The microhydration effects on the geometrical structures, adsorption site, and the proton transfer reaction of the adenine–thymine anion are investigated. The site N10 atom of the adenine moiety has a larger proton affinity than the site O24 atom of thymine, which facilitates the proton H26 transfers from the N25 site of thymine to the N10 site of adenine. Therefore, the first single-proton transfer pathway (SPT1) is found for the all of the monohydrated adenine–thymine anions (AN4T)−·H2O, (AN13T)−·H2O, (ATO24)−·H2O, and (ATO28)−·H2O and tetrahydrated adenine–thymine anions (AT)−·4H2O. The proton H9 at the N7 site of adenine is also found to transfer to the O24 site of thymine for (AN4T)−·H2O and (AN13T)−·H2O in the gas phase. The double-proton transferred pathway is found when one water molecule interacts with the O28 atom of thymine. The reactant structures before the proton transfer are more stable than the product structures, and the structural changes mainly occur in thymine. The reaction energies of the microhydrated adenine–thymine anion in the gas phase and in the aqueous environment predict that the proton transfer process of the microhydrated adenine–thymine anion are more favorable in the gas phase than in aqueous solution.

π- vs σ-radical states of one-electron-oxidized DNA/RNA bases: a density functional theory study

As a result of their inherent planarity, DNA base radicals generated by one-electron oxidation/reduction or bond cleavage form π- or σ-radicals. While most DNA base systems form π-radicals, there are a number of nucleobase analogues such as one-electron-oxidized 6-azauraci1, 6-azacytosine, and 2-thiothymine or one-electron reduced 5-bromouracil that form more reactive σ-radicals. Elucidating the availability of these states within DNA, base radical electronic structure is important to the understanding of the reactivity of DNA base radicals in different environments. In this work, we address this question by the calculation of the relative energies of π- and σ-radical states in DNA/RNA bases and their analogues. We used density functional theory B3LYP/6-31++G** method to optimize the geometries of π- and σ-radicals in Cs symmetry (i.e., planar) in the gas phase and in solution using the polarized continuum model (PCM). The calculations predict that σ- and π-radical states in one-electron-oxidized bases of thymine, T(N3-H)(•), and uracil, U(N3-H)(•), are very close in energy; i.e., the π-radical is only ca. 4 kcal/mol more stable than the σ-radical. For the one-electron-oxidized radicals of cytosine, C(•+), C(N4-H)(•), adenine, A(•+), A(N6-H)(•), and guanine, G(•+), G(N2-H)(•), G(N1-H)(•), the π-radicals are ca. 16-41 kcal/mol more stable than their corresponding σ-radicals. Inclusion of solvent (PCM) is found to stabilize the π- over σ-radical of each of the systems. U(N3-H)(•) with three discrete water molecules in the gas phase is found to form a three-electron σ bond between the N3 atom of uracil and the O atom of a water molecule, but on inclusion of full solvation and discrete hydration, the π-radical remains most stable.

One-electron oxidation of neutral sugar radicals of 2'-deoxyguanosine and 2'-deoxythymidine: a density functional theory (DFT) study

One electron oxidation of neutral sugar radicals has recently been suggested to lead to important intermediates in the DNA damage process culminating in DNA strand breaks. In this work, we investigate sugar radicals in a DNA model system to understand the energetics of sugar radical formation and oxidation. The geometries of neutral sugar radicals C(1')(•), C(2')(•), C(3')(•), C(4')(•), and C(5')(•) of 2'-deoxyguanosine (dG) and 2'-deoxythymidine (dT) were optimized in the gas phase and in solution using the B3LYP and ωB97x functionals and 6-31++G(d) basis set. Their corresponding cations (C(1')(+), C(2')(+), C(3')(+), C(4')(+), and C(5')(+)) were generated by removing an electron (one-electron oxidation) from the neutral sugar radicals, and their geometries were also optimized using the same methods and basis set. The calculation predicts the relative stabilities of the neutral sugar radicals in the order C(1')(•) > C(4')(•) > C(5')(•) > C(3')(•) > C(2')(•), respectively. Of the neutral sugar radicals, C(1')(•) has the lowest vertical ionization potential (IP(vert)), ca. 6.33 eV in the gas phase and 4.71 eV in solution. C(2')(•) has the highest IP(vert), ca. 8.02 eV, in the gas phase, and the resultant C(2') cation is predicted to undergo a barrierless hydride transfer from the C(1') site to produce the C(1') cation. One electron oxidation of C(2')(•) in dG is predicted to result in a low lying triplet state consisting of G(+•) and C(2')(•). The 5',8-cyclo-2'-deoxyguanosin-7-yl radical formed by intramolecular bonding between C(5')(•) and C(8) of guanine transfers spin density from C(5') site to guanine, and this structure has IP(vert) 6.25 and 5.48 eV in the gas phase and in solution.

Electron attachment to the guanine-cytosine nucleic acid base pair and the effects of monohydration and proton transfer

The guanine-cytosine (GC) radical anion and its interaction with a single water molecule is studied using ab initio and density functional methods. Z-averaged second-order perturbation theory (ZAPT2) was applied to GC radical anion for the first time. Predicted spin densities show that the radical character is localized on cytosine. The Watson-Crick monohydrated GC anion is compared to neutral GC·H2O, as well as to the proton-transferred analogue on the basis of structural and energetic properties. In all three systems, local minima are identified that correspond to water positioned in the major and minor grooves of macromolecular DNA. On the anionic surface, two novel structures have water positioned above or below the GC plane. On the neutral and anionic surfaces, the global minimum can be described as water interacting with the minor groove. These structures are predicted to have hydration energies of 9.7 and 11.8 kcal mol(-1), respectively. Upon interbase proton-transfer (PT), the anionic global minimum has water positioned in the major groove, and the hydration energy increases to 13.4 kcal mol(-1). PT GC·H2O(•-) has distonic character; the radical character resides on cytosine, while the negative charge is localized on guanine. The effects of proton transfer are further investigated through the computed adiabatic electron affinities (AEA) of GC and monohydrated GC, and the vertical detachment energies (VDE) of the corresponding anions. Monohydration increases the AEAs and VDEs by only 0.1 eV, while proton-transfer increases the VDEs substantially (0.8 eV). The molecular charge distribution of monohydrated guanine-cytosine radical anion depends heavily on interbase proton transfer.

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.

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.