Base Release in Nucleosides Induced by Low-Energy Electrons: A DFT Study

Prebiotic chemical origin of biomolecular complementarity

The early Earth, devoid of the protective stratospheric ozone layer, must have sustained an ambient prebiotic physicochemical medium intensified by the co-existence of shortwave UV photons and very low energy electrons (vLEEs). Consequently, only intrinsically stable molecules against these two co-existing molecular destructors must have proliferated and thereby chemically evolved into the advanced molecules of life. Based on this view, we examined the stability inherent in nucleobases and their complementary pairs as resistance to the molecular damaging effects of shortwave UV photons and vLEEs. This leads to the conclusion that nucleobases could only proliferated as their complementary pairs under the unfavorable prebiotic conditions on early Earth. The complementary base pairing not only enhances but consolidates the intrinsic stability of nucleobases against short-range UV photons, vLEEs, and possibly many as-yet-unknown deleterious agents co-existed in the prebiotic conditions of the early Earth. In short, complementary base pairing is a manifestation of chemical evolution in the unfavorable prebiotic medium created by the absence of the stratospheric ozone layer.

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.

Why 6-Iodouridine Cannot Be Used as a Radiosensitizer of DNA Damage? Computational and Experimental Studies

Previous density functional theory (DFT) studies on 6-brominated pyrimidine nucleosides suggest that 6-iodo-2'-deoxyuridine (6IdU) should act as a better radiosensitizer than its 5-iodosubstituted 2'-deoxyuridine analogue. In this work, we show that 6IdU is unstable in an aqueous solution. Indeed, a complete disappearance of the 6IdU signal was observed during its isolation by reversed-phase high-performance liquid chromatography (RP-HPLC). As indicated by the thermodynamic characteristics for the S_N1-type hydrolysis of 6IdU obtained at the CAM-B3LYP/DGDZVP++ level and the polarizable continuum model (PCM) of water, 6-iodouracil (6IU) was already released quantitatively at ambient temperatures. The simulation of the hydrolysis kinetics demonstrated that a thermodynamic equilibrium was reached within seconds for the title compound. To assess the reliability of the calculations carried out, we synthesized 6-iodouridine (6IUrd), which was, unlike 6IdU, sufficiently stable in an aqueous solution at room temperature. The activation barrier for the N-glycosidic bond dissociation in 6IUrd was estimated experimentally using an Arrhenius plot. The stabilities in water calculated for 6IdU, 6IUrd, and 5-iodo-2'-deoxyuridine (5IdU) could be explained by the electronic and steric effects of the 2'-hydroxy group present in the ribose moiety. Our studies highlight the issue of the hydrolytic stability of potentially radiosensitizing nucleotides which, besides having favorable dissociative electron attachment (DEA) characteristics, must be stable in water to have any practical application.

Electrophilic Properties of 2'-Deoxyadenosine···Thymine Dimer: Photoelectron Spectroscopy and DFT Studies

The anion radical of the 2'-deoxyadenosine···thymine (dAT^•-) pair has been investigated experimentally and theoretically in the gas phase. By employing negative-ion photoelectron spectroscopy (PES), we have registered a spectrum typical for the valence-bound anion, featuring a broad peak at the electron-binding energy (EBE) between ∼1.5 and 2.2 eV with the maximum at ∼1.9 eV. The measured value of the adiabatic electron affinity (AEA) for dAT was estimated to be ∼1.1 eV. Calculations performed at the M06-2X/6-31++G(d,p) level revealed that the structure, where thymine is coordinated to the sugar of dA by two hydrogen bonds, is responsible for the observed PES signal. The AEA_G and the vertical detachment energy of 0.91 and 1.68 eV, respectively, calculated for this structure reproduce the experimental values well. The role of the possible proton transfer in the stabilization of anionic radical complexes is discussed.

First-Principles Simulations of Biological Molecules Subjected to Ionizing Radiation

Ionizing rays cause damage to genomes, proteins, and signaling pathways that normally regulate cell activity, with harmful consequences such as accelerated aging, tumors, and cancers but also with beneficial effects in the context of radiotherapies. While the great pace of research in the twentieth century led to the identification of the molecular mechanisms for chemical lesions on the building blocks of biomacromolecules, the last two decades have brought renewed questions, for example, regarding the formation of clustered damage or the rich chemistry involving the secondary electrons produced by radiolysis. Radiation chemistry is now meeting attosecond science, providing extraordinary opportunities to unravel the very first stages of biological matter radiolysis. This review provides an overview of the recent progress made in this direction, focusing mainly on the atto- to femto- to picosecond timescales. We review promising applications of time-dependent density functional theory in this context.

What are the Potential Sites of DNA Attack by N-Acetyl-p-benzoquinone Imine (NAPQI)?

Metabolic bioactivation of small molecules can produce electrophilic metabolites that can covalently modify proteins and DNA. Paracetamol (APAP) is a commonly used over-the-counter analgesic, and its hepatotoxic side effects have been postulated to be due to the formation of the electrophilic metabolite N-acetylbenzoquinone imine (NAPQI). It has been established that NAPQI reacts to form covalent bonds to the side-chain functional groups of cysteine, methionine, tyrosine, and tryptophan residues. While there have been scattered reports that APAP can form adducts with DNA the nature of these adducts have not yet been fully characterised. Here the four deoxynucleosides, deoxyguanosine (dG), deoxyadenosine (dA), deoxycytidine (dC), and deoxythymidine (dT) were reacted with NAPQI and the formation of adducts was profiled using liquid chromatography–mass spectrometry with positive-ion mode electrospray ionisation and collision-induced dissociation. Covalent adducts were detected for dG, dA, and dC and tandem mass spectrometry (MS/MS) spectra revealed common neutral losses of deoxyribose (116 amu) arising from cleavage of the glyosidic bond with formation of the modified nucleobase. Of the four deoxynucleosides, dC proved to be the most reactive, followed by dG and dA. A pH dependence was found, with greater reactivity being observed at pH 5.5. The results of density functional theory calculations aimed at understanding the relative reactivities of the four deoxynucleosides towards NAPQI are described.

Strand Breaks Induced by Very Low Energy Electrons: Product Analysis and Mechanistic Insight into the Reaction with TpT

Numerous experimental studies show that 5-15 eV electrons induce strand breaks in DNA at energies below the ionization threshold of DNA components. In this energy range, DNA damage arises principally by the formation of transient negative ions, decaying into dissociative electron attachment (DEA) and electronic excitation of dissociative states. Here, we carried out LC-MS/MS analysis of the degradation products arising from bombardment of TpT, a DNA model compound, irradiated with very low energy electrons (vLEEs; ∼1.8 eV). The formation of thymidine 5'-monophosphate (TMP5') together with 2',3'-dideoxythymidine (ddT3') can be explained by cleavage of the C3'-O bond of TpT, whereas thymidine 3'-monophosphate (TMP3') and 2',5'-dideoxythymidine (ddT5') are formed by cleavage of the C5'-O bond. The formation of ddT3' and ddT5' decreased upon irradiation of either TMP5' or TMP3', and even further in the case of thymidine, underlining the critical role of the phosphate group. Interestingly, the yield of TMP5' and TMP3' was higher than that of the corresponding ddT3' and ddT5' products, suggesting alternative fates of C3' and C5'-centered sugar radicals. In contrast, the release of thymine from TpT was minor (<20%) and did not result in the formation of expected products from DEA-mediated cleavage at the N-glycosidic bond. Lastly, vLEE induced the conversion of thymine to 5,6-dihydrothymine (5,6-dhT) within TpT, a reaction likely involving thymine anion radicals. In summary, we show that a major pathway of vLEEs involves DEA-mediated cleavage of the C3'-O and C5'-O bonds of TpT, resulting in the formation of specific fragments, which represent a prompt single strand break in DNA.

Electron-Induced Dissociation of the Potential Radiosensitizer 5-Selenocyanato-2'-deoxyuridine

5-Selenocyanato-2'-deoxyuridine (SeCNdU) is a recently proposed radiosensitizer based on 2'-deoxyuridine (dU) with the electron-affinic selenocyanato (-SeCN) side group attached at the C5 position of uracil. Since electron interaction processes may be an important source of DNA damage by ionizing radiation, we have studied low-energy dissociative electron attachment to SeCNdU in the gas phase. Negative ion formation has been obtained by means of mass spectrometry, where a rich fragmentation pattern is observed even at ∼0 eV. The reaction pathways exhibiting the highest ion yields are C₄N₂O₂H₂Se^•- and CN^-, both involving a cleavage of the Se-CN bond. The heaviest fragment anion observed is C₉N₂O₅H₁₀Se^•-, where besides the charged species, the hydrogen and cyano radicals are also formed. Further decomposition channels also yield the highly reactive hydroxyl radical, which possesses a high DNA damage potential. All observed channels have experimentally determined onsets at 0 eV, which are supported by calculations performed at the M06-2X/aug-cc-pVTZ level. The calculations comprise the thermochemical thresholds at standard and experimental (428.15 K, 3 × 10^-11 atm) conditions together with the adiabatic electron affinities. The present study shows that low-energy electrons very effectively decompose SeCNdU upon attachment of thermal electrons, producing a large variety of charged fragments and radicals.

Notable effect of water on excess electron attachment to aqueous DNA deoxyribonucleosides

As excess electrons are vertically attached to aqueous deoxyribonucleosides, ∼50% of excess electrons would be delocalized over the water molecules.

Base Release and Modification in Solid-Phase DNA Exposed to Low-Energy Electrons

Ionization generates a large number of secondary low-energy electrons (LEEs) with a most probable energy of approximately 10 eV, which can break DNA bonds by dissociative electron attachment (DEA) and lead to DNA damage. In this study, we investigated radiation damage to dry DNA induced by X rays (1.5 keV) alone on a glass substrate or X rays combined with extra LEEs (average energy of 5.8 eV) emitted from a tantalum (Ta) substrate under an atmosphere of N₂ and standard ambient conditions of temperature and pressure. The targets included calf-thymus DNA and double-stranded synthetic oligonucleotides. We developed analytical methods to measure the release of non-modified DNA bases from DNA and the formation of several base modifications by LC-MS/MS with isotopic dilution for precise quantification. The results show that the yield of non-modified bases as well as base modifications increase by 20-30% when DNA is deposited on a Ta substrate compared to that on a glass substrate. The order of base release (Gua > Ade > Thy ∼ Cyt) agrees well with several theoretical studies indicating that Gua is the most susceptible site toward sugar-phosphate cleavage. The formation of DNA damage by LEEs is explained by DEA leading to the release of non-modified bases involving the initial cleavage of N1-C1', C3'-O3' or C5'-O5' bonds. The yield of base modifications was lower than the release of non-modified bases. The main LEE-induced base modifications include 5,6-dihydrothymine (5,6-dHT), 5,6-dihydrouracil (5-dHU), 5-hydroxymethyluracil (5-HmU) and 5-formyluracil (5-ForU). The formation of base modifications by LEEs can be explained by DEA and cleavage of the C-H bond of the methyl group of Thy (giving 5-HmU and 5-ForU) and by secondary reactions of H atoms and hydride anions that are generated by primary LEE reactions followed by subsequent reaction with Cyt and Thy (giving 5,6-dHU and 5,6-dHT).

Low energy electron induced cytosine base release in 2'-deoxycytidine-3'-monophosphate via glycosidic bond cleavage: a time-dependent wavepacket study

Low energy electron (LEE) induced cytosine base release in a selected pyrimidine nucleotide, viz., 2'-deoxycytidine-3'-monophosphate is investigated using ab initio electronic structure methods and time dependent quantum mechanical calculations. It has been noted that the cytosine base scission is comparatively difficult process than the 3' C-O bond cleavage from the lowest π* shape resonance in energy region <1 eV. This is mainly due to the high activation energy barrier associated with the electron transfer from the π* orbital of the base to the σ* orbital of the glycosidic N-C bond. In addition, the metastable state formed after impinging LEE (0-1 eV) has very short lifetime (10 fs) which may decay in either of the two competing auto-detachment or dissociation process simultaneously. On the other hand, the selected N-C mode may cleave to form the cytosine base anion at higher energy regions (>2 eV) via tunneling of the glycosidic bond. Resonance states generated within this energy regime will exist for a duration of ~35-55 fs. Comparison of salient features of the two dissociation events, i.e., 3' C-O single strand break and glycosidic N-C bond cleavage in 3'-dCMPH molecule are also provided.

DFT reinvestigation of DNA strand breaks induced by electron attachment

The benchmark study of DFT methods on the activation energies of phosphodiester C3'-O and C5'-O bond ruptures and glycosidic C1'-N bond ruptures induced by electron attachment was performed. While conventional pure and hybrid functionals provide a relatively reasonable description for the C1'-N bond rupture, they significantly underestimate the energy barriers of the C-O bond ruptures. This is because the transition states of the later reactions, which are characterized by an electron distribution delocalized from the nucleobase to sugar-phosphate backbone, suffer from a severe self-interaction error in common DFT methods. CAM-B3LYP, M06-2X, and ωB97XD are the top three methods that emerged from the benchmark study; the mean absolute errors relative to the CCSD(T) values are 1.7, 1.9, and 2.2 kcal/mol, respectively. The C-O bond cleavages of 3'- and 5'-dXMP(•-), where X represents four nucleobases, were then recalculated at the M06-2X/6-31++G**//M06-2X/6-31+G* level, and it turned out that the C-O bond cleavages do not proceed as easily as previously predicted by the B3LYP calculations. Our calculations revealed that the C-O bonds of purine nucleotides are more susceptible than pyrimidine nucleotides to the electron attachment. The energies of electron attachment to nucleotides were calculated and discussed as well.

Electron interaction with phosphate cytidine oligomer dCpdC: base-centered radical anions and their electronic spectra

Computational chemistry approach was applied to explore the nature of electron attachment to cytosine-rich DNA single strands. An oligomer dinucleoside phosphate deoxycytidylyl-3',5'-deoxycytidine (dCpdC) was selected as a model system for investigations by density functional theory. Electron distribution patterns for the radical anions of dCpdC in aqueous solution were explored. The excess electron may reside on the nucleobase at the 5' position (dC(•-)pdC) or at the 3' position (dCpdC(•-)). From comparison with electron attachment to the cytosine related DNA fragments, the electron affinity for the formation of the cytosine-centered radical anion in DNA is estimated to be around 2.2 eV. Electron attachment to cytosine sites in DNA single strands might cause perturbations of local structural characteristics. Visible absorption spectroscopy may be applied to validate computational results and determine experimentally the existence of the base-centered radical anion. The time-dependent DFT study shows the absorption around 550-600 nm for the cytosine-centered radical anions of DNA oligomers. This indicates that if such species are detected experimentally they would be characterized by a distinctive color.

Photoelectron spectroscopy and density functional theory studies on the uridine homodimer radical anions

We report the photoelectron spectrum (PES) of the homogeneous dimer anion radical of uridine, (rU)(2)(●-). It features a broad band consisting of an onset of ∼1.2 eV and a maximum at the electron binding energy (EBE) ranging from 2.0 to 2.5 eV. Calculations performed at the B3LYP∕6-31++G∗∗ level of theory suggest that the PES is dominated by dimeric radical anions in which one uridine nucleoside, hosting the excess charge on the base moiety, forms hydrogen bonds via its O8 atom with hydroxyl of the other neutral nucleoside's ribose. The calculated adiabatic electron affinities (AEAGs) and vertical detachment energies (VDEs) of the most stable homodimers show an excellent agreement with the experimental values. The anionic complexes consisting of two intermolecular uracil-uracil hydrogen bonds appeared to be substantially less stable than the uracil-ribose dimers. Despite the fact that uracil-uracil anionic homodimers are additionally stabilized by barrier-free electron-induced proton transfer, their relative thermodynamic stabilities and the calculated VDEs suggest that they do not contribute to the experimental PES spectrum of (rU)(2)(●-).

Low-energy electron-induced DNA damage: effect of base sequence in oligonucleotide trimers

DNA damage induced by low-energy electrons (LEEs) has attracted considerable attention in recent years because LEEs represent a large percentage of the total energy deposited by ionizing radiation and because LEEs have been shown to damage DNA components. In this article, we have studied the effect of base sequences in a series of oligonucleotide trimers by the analysis of damage remaining within the nonvolatile condensed phase after LEE irradiation. The model compounds include TXT, where X represents one of the four normal bases of DNA (thymine (T), cytosine (C), adenine (A), and guanine (G)). Using HPLC-UV analysis, several known fragments were quantified from the release of nonmodified nucleobases (T and X) as well as from phosphodiester C-O bond cleavage (pT, pXT, Tp, and TXp). The total damage was estimated by the disappearance of the parent peaks in the chromatogram of nonirradiated and irradiated samples. When trimers were irradiated with LEE (10 eV), the total damage decreased 2-fold in the following order: TTT > TCT > TAT > TGT. The release of nonmodified nuclobases (giving from 17 to 24% of the total products) mainly occurred from the terminal sites of trimers (i.e., T) whereas the release of central nucleobases was minor (C) or not at all detected (A and G). In comparison, the formation of products arising from phosphodiester bond cleavage accounted for 9 to 20% of the total damage and it partitioned to the four possible sites of cleavage present in trimers. This study indicates that the initial LEE capture and subsequent bond breaking within the intermediate anion depend on the sequence and electron affinity of the bases, with the most damage attributed to the most electronegative base, T.

Electron attachment-induced DNA single-strand breaks at the pyrimidine sites

To elucidate the contribution of pyrimidine in DNA strand breaks caused by low-energy electrons (LEEs), theoretical investigations of the LEE attachment-induced C_3′–O_3′, and C_5′–O_5′ σ bond as well as N-glycosidic bond breaking of 2′-deoxycytidine-3′,5′-diphosphate and 2′-deoxythymidine-3′,5′-diphosphate were performed using the B3LYP/DZP++ approach. The base-centered radical anions are electronically stable enough to assure that either the C–O or glycosidic bond breaking processes might compete with the electron detachment and yield corresponding radical fragments and anions. In the gas phase, the computed glycosidic bond breaking activation energy (24.1 kcal/mol) excludes the base release pathway. The low-energy barrier for the C_3′–O_3′ σ bond cleavage process (∼6.0 kcal/mol for both cytidine and thymidine) suggests that this reaction pathway is the most favorable one as compared to other possible pathways. On the other hand, the relatively low activation energy barrier (∼14 kcal/mol) for the C_5′–O_5′ σ bond cleavage process indicates that this bond breaking pathway could be possible, especially when the incident electrons have relatively high energy (a few electronvolts). The presence of the polarizable medium greatly increases the activation energies of either C–O σ bond cleavage processes or the N-glycosidic bond breaking process. The only possible pathway that dominates the LEE-induced DNA single strands in the presence of the polarizable surroundings (such as in an aqueous solution) is the C_3′–O_3′ σ bond cleavage (the relatively low activation energy barrier, ∼13.4 kcal/mol, has been predicted through a polarizable continuum model investigation). The qualitative agreement between the ratio for the bond breaks of C_5′–O_5′, C_3′–O_3′ and N-glycosidic bonds observed in the experiment of oligonucleotide tetramer CGAT and the theoretical sequence of the bond breaking reaction pathways have been found. This consistency between the theoretical predictions and the experimental observations provides strong supportive evidences for the base-centered radical anion mechanism of the LEE-induced single-strand bond breaking around the pyrimidine sites of the DNA single strands.

Stability of the valence anion of cytosine is governed by nucleobases sequence in the double stranded DNA pi-stack: A computational study

The stabilities of the valence anion of cytosine (C(-)) in model trimers of complementary base pairs that possess the B-DNA geometry but differ in base sequence are reported. In order to estimate the energetics of electron attachment to the middle cytosine incorporated in the trimer, a thermodynamic cycle employing all possible two-body interaction energies in the neutral and anionic duplex as well as the adiabatic electron affinity of isolated cytosine were developed. All calculations were carried out at the MP2 level of theory with the aug-cc-pVDZ basis set. We have demonstrated that contrary to the literature reports, concerning single stranded DNA, the sequence of nucleic bases has a profound effect on the stability of the cytosine valence anion. The anionic 3(')-CCC-5(') complex is the most stable configuration (EA=0.399 eV) and the 3(')-GCG-5(') trimer anion is the most unstable species (EA=-0.193 eV). Moreover, with the energetic correction for the presence of sugar-phosphate backbone all possible double stranded DNA sequences lead to the stable C(-). The predicted electron affinities of the cytosine anion have been compared to the results of analogous studies on the thymine anion published recently [M. Kobyłecka et al., J. Am. Chem. Soc. 130, 15683 (2008)]. The consequences of low-energy barrier proton transfer in the GC anion have been discussed in the context of induced by electrons DNA single strand breaks. The DNA sequences that should dramatically differ in their vulnerability to be damaged by low energy electrons have been proposed.

Mechanisms of strand break formation in DNA due to the direct effect of ionizing radiation: the dependency of free base release on the length of alternating CG oligodeoxynucleotides

The question of how NA base sequence influences the yield of DNA strand breaks produced by the direct effect of ionizing radiation was investigated in a series of oligodeoxynucleotides of the form (d(CG)(n))(2) and (d(GC)(n))(2). The yields of free base release from X-irradiated DNA films containing 2.5 waters/nucleotide were measured by HPLC as a function of oligomer length. For (d(CG)(n))(2), the ratio of the Gua yield to Cyt yield, R, was relatively constant at 2.4-2.5 for n = 2-4 and it decreased to 1.2 as n increased from 5 to 10. When Gua was moved to the 5' end, for example going from d(CG)(5) to d(GC)(5), R dropped from 1.9 +/- 0.1 to 1.1 +/- 0.1. These effects are poorly described if the chemistry at the oligomer ends is assumed to be independent of the remainder of the oligomer. A mathematical model incorporating charge transfer through the base stack was derived to explain these effects. In addition, EPR was used to measure the yield of trapped-deoxyribose radicals at 4 K following X-irradiation at 4 K. The yield of free base release was substantially greater, by 50-100 nmol/J, than the yield of trapped-deoxyribose radicals. Therefore, a large fraction of free base release stems from a nonradical intermediate. For this intermediate, a deoxyribose carbocation formed by two one-electron oxidations is proposed. This reaction pathway requires that the hole (electron loss site) transfers through the base stack and, upon encountering a deoxyribose hole, oxidizes that site to form a deoxyribose carbocation. This reaction mechanism provides a consistent way of explaining both the absence of trapped radical intermediates and the unusual dependence of free base release on oligomer length.

Effects of electron attachment on C5'-O5' and C1'-N1 bond cleavages of pyrimidine nucleotides: A theoretical study

Sugar-base C(1')-N(1) and phosphate-sugar C(5')-O(5') bond breakings of 2'-deoxycytidine-5'-monophosphates (dCMP) and 2'-deoxythymidine-5'- monophosphates (dTMP) and their radical anions have been explored theoretically at the B3LYP/DZP++ level of theory. Calculations show that the low-energy electrons attachment to the pyrimidine nucleotides results in remarkable structural and chemical bonding changes. Predicted Gibbs free energies of reaction DeltaG for the C(5')-O(5') bond dissociation process of the radical anions are -14.6 and -11.5 kcal mol(-1), respectively, and such dissociation processes may be intrinsically spontaneous in the gas phase. Furthermore, the C(5')-O(5') bond cleavage processes of the anionic dCMP and dTMP were predicted to have activation energies of 6.9 and 8.0 kcal mol(-1) in the gas phase, respectively, much lower than the barriers for the C(1')-N(1) bond breaking process, showing that the C-O bond dissociation in DNA single strand breaks is a dominant process as observed experimentally.

Electron attachment to the hydrogenated Watson-Crick guanine cytosine base pair (GC+H): conventional and proton-transferred structures

The anionic species resulting from hydride addition to the Watson-Crick guanine-cytosine (GC) DNA base pair are investigated theoretically. Proton-transferred structures of GC hydride, in which proton H1 of guanine or proton H4 of cytosine migrates to the complementary base-pair side, have been studied also. All optimized geometrical structures are confirmed to be minima via vibrational frequency analyses. The lowest energy structure places the additional hydride on the C6 position of cytosine coupled with proton transfer, resulting in the closed-shell anion designated 1T (G(-)C(C6)). Energetically, the major groove side of the GC pair has a greater propensity toward hydride/hydrogen addition than does the minor grove side. The pairing (dissociation) energy and electron-attracting ability of each anionic structure are predicted and compared with those of the neutral GC and the hydrogenated GC base pairs. Anion 8T (G(O6)C(-)) is a water-extracting complex and has the largest dissociation energy. Anion 2 (GC(C4)(-)) and the corresponding open-shell radical GC(C4) have the largest vertical electron detachment energy and adiabatic electron affinity, respectively. From the difference between the dissociation energy and electron-removal ability of the normal GC anion and the most favorable structure of GC hydride, it is clear that one may dissociate the GC anion and maintain the integrity of the GC hydride.

Theoretical studies of damage to 3'-uridine monophosphate induced by electron attachment

Low-energy electrons (LEE) are well known to induce nucleic acid damage. However, the damage mechanisms related to charge state and structural features remain to be explored in detail. In the present work, we have investigated the N1-glycosidic and C3'-O(P) bond ruptures of 3'-UMP (UMP=uridine monophosphate) and the protonated form 3'-UMPH with -1 and zero charge, respectively, based on hybrid density functional theory (DFT) B3 LYP together with the 6-31+G(d,p) basis set. The glycosidic bond breakage reactions of the 3'UMP and 3'UMPH electron adducts are exothermic in both cases, with barrier heights of 19-20 kcal mol(-1) upon inclusion of bulk solvation. The effects of the charge state on the phosphate group are marginal, but the C2'-OH group destabilizes the transition structure of glycosidic bond rupture of 3'-UMPH in the gas phase by approximately 5.0 kcal mol(-1). This is in contrast with the C3'-O(P) bond ruptures induced by LEE in which the charge state on the phosphate influences the barrier heights and reaction energies considerably. The barrier towards C3'-O(P) bond dissociation in the 3'UMP electron adduct is higher in the gas phase than the one corresponding to glycosidic bond rupture and is dramatically influenced by the C2'-OH group and bulk salvation, which decreases the barrier to 14.7 kcal mol(-1). For the C3'-O(P) bond rupture of the 3'UMPH electron adduct, the reaction is exothermic and the barrier is even lower, 8.2 kcal mol(-1), which is in agreement with recent results for 3'-dTMPH and 5'-dTMPH (dTMPH=deoxythymidine monophosphate). Both the Mulliken atomic charges and unpaired-spin distribution play significant roles in the reactions.

Photoelectron spectroscopy of the parent anions of the nucleotides, adenosine-5'-monophosphate and 2'deoxyadenosine-5'-monophosphate

The parent anions of the nucleotides, adenosine-5(')-monophosphate (AMPH) and 2(')deoxyadenosine-5(')-monophosphate (dAMPH) were generated in a novel source and their photoelectron spectra recorded with 3.49 eV photons. Vertical detachment energy (VDE) and the adiabatic electron affinity (EA(a)) values were extracted from each of the two spectra. Concurrently, Kobylecka et al. [J. Chem. Phys. 128, 044315 (2008)] conducted calculations which explored electron attachment to dAMPH. Based on the agreement between their calculated and our measured VDE and EA(a) values, we conclude that the dAMPH(-) anions studied in these experiments were formed by electron-induced, intramolecular, (barrier-free) proton-transfer as predicted by the calculations. Given the similarities between the photoelectron spectra of dAMPH(-) and AMPH(-), it is likely that AMPH(-) can be described in the same manner.

Barrier-free proton transfer in the valence anion of 2'-deoxyadenosine-5'-monophosphate. II. A computational study

The propensity of four representative conformations of 2(')-deoxyadenosine-5(')-monophosphate (5(')-dAMPH) to bind an excess electron has been studied at the B3LYP6-31++G(d,p) level. While isolated canonical adenine does not support stable valence anions in the gas phase, all considered neutral conformations of 5(')-dAMPH form adiabatically stable anions. The type of an anionic 5(')-dAMPH state, i.e., the valence, dipole bound, or mixed (valence/dipole bound), depends on the internal hydrogen bond(s) pattern exhibited by a particular tautomer. The most stable anion results from an electron attachment to the neutral syn-south conformer. The formation of this anion is associated with a barrier-free proton transfer triggered by electron attachment and the internal rotation around the C4(')-C5(') bond. The adiabatic electron affinity of the a_south-syn anion is 1.19 eV, while its vertical detachment energy is 1.89 eV. Our results are compared with the photoelectron spectrum (PES) of 5(')-dAMPH(-) measured recently by Stokes et al., [J. Chem. Phys. 128, 044314 (2008)]. The computational VDE obtained for the most stable anionic structure matches well with the experimental electron binding energy region of maximum intensity. A further understanding of DNA damage might require experimental and computational studies on the systems in which purine nucleotides are engaged in hydrogen bonding.

Which DNA damage is likely to be relevant in hormetic responses?

Working under the assumption that hormesis is triggered by specific types of DNA damage, this report focuses on the types of damage which form the signature of ionizing radiation. The key attribute of the signature is the clustering of damage, arising from clusters of energy deposition such that more than one site within a 10 base pair segment of DNA has been chemically altered. A brief overview is given on what is currently believed to be the primary components of clustered damage produced by the direct effect. The overview draws primarily on studies that utilize electron paramagnetic resonance to measure free radical intermediates and gel electrophoresis to measure clustered damage in plasmid DNA. Based on this information, the threshold for a radiation induced biological response is calculated.

Intrinsic electrophilic properties of nucleosides: Photoelectron spectroscopy of their parent anions

The nucleoside parent anions 2(')-deoxythymidine(-), 2(')-deoxycytidine(-), 2(')-deoxyadenosine(-), uridine(-), cytidine(-), adenosine(-), and guanosine(-) were generated in a novel source, employing a combination of infrared desorption, electron photoemission, and a gas jet expansion. Once mass selected, the anion photoelectron spectrum of each of these was recorded. In the three cases in which comparisons were possible, the vertical detachment energies and likely adiabatic electron affinities extracted from these spectra agreed well with the values calculated both by Richardson et al. [J. Am. Chem. Soc. 126, 4404 (2004)] and by Li et al. [Radiat. Res. 165, 721 (2006)]. Through the combination of our experimental results and their theoretical calculations, several implications emerge. (1) With the possible exception of dG(-), the parent anions of nucleosides exist, and they are stable. (2) These nucleoside anions are valence anions, and in most cases the negative charge is closely associated with the nucleobase moiety. (3) The nucleoside parent anions we have generated and studied are the negative ions of canonical, neutral nucleosides, similar to those found in DNA.

Density Functional Study toward Understanding Dehydrogenation of the Adenine−Thymine Base Pair and Its Anion

The dehydrogenated radicals and anions of Watson-Crick adenine-thymine (A-T) base pair have been investigated by the B3LYP/DZP++ approach. Calculations show that the dehydrogenated radicals and anions have relatively high stabilities compared with the single base adenine and thymine. The electron attachment to the A-T base pair and its derivatives significantly modifies the hydrogen bond interactions and results in remarkable structural changes. As for the dehydrogenated A-T radicals, they have relatively high electron affinities and different dehydrogenation properties with respect to their constituent elements. The relatively low-cost hydrogen eliminations correspond to the (N9)-H (adenine) and (N1)-H (thymine) bonds cleavage. Both dehydrogenation processes have Gibbs free energies of reaction DeltaG degrees of 13.4 and 17.2 kcal mol-1, respectively. The solvent water exhibits significant effect on electron attachment and dehydrogenation properties of the A-T base pair and its derivatives. In the dehydrogenating process, the anionic A-T fragment gradually changes its electronic configuration from pi* to sigma* state, like the single bases adenine and thymine.

Low-energy electron attachment to 5'-thymidine monophosphate: modeling single strand breaks through dissociative electron attachment

Mechanisms of low-energy electron (LEE) attachment and subsequent single-strand break (SSB) formation are investigated by density functional theory treatment of a simple model for DNA, i.e., the nucleotide, 5'-thymidine monophosphate (5'-dTMPH). In the present study, the C5'-O5' bond dissociation due to LEE attachment has been followed along the adiabatic as well as on the vertical (electron attached to the optimized geometry of the neutral molecule) anionic surfaces using B3LYP functional and 6-31G* and 6-31++G** basis sets. Surprisingly, it is found that the PES of C5'-O5' bond dissociation in the anion radicals have approximately the same barrier for both adiabatic and vertical pathways. These results provide support for the hypothesis that transiently bound electrons (shape resonances) to the virtual molecular orbitals of the neutral molecule likely play a key role in the cleavage of the sugar-phosphate C5'-O5' bond in DNA resulting in the direct formation of single strand breaks without significant molecular relaxation. To take into account the solvation effects, we considered the neutral and anion radical of 5'-dTMP surrounded by 5 or 11 water molecules with Na+ as a counterion. These structures were optimized using the B3LYP/6-31G** level of theory. We find the barrier height for adiabatic C5'-O5' bond dissociation of 5'-dTMP anion radical in aqueous environment is so substantially higher than in the gas phase that the adiabatic route will not contribute to DNA strand cleavage in aqueous systems. This result is in agreement with experiment.

Interaction of low-energy electrons with the purine bases, nucleosides, and nucleotides of DNA

The authors report results from computational studies of the interaction of low-energy electrons with the purine bases of DNA, adenine and guanine, as well as with the associated nucleosides, deoxyadenosine and deoxyguanosine, and the nucleotide deoxyadenosine monophosphate. Their calculations focus on the characterization of the pi* shape resonances associated with the bases and also provide general information on the scattering of slow electrons by these targets. Results are obtained for adenine and guanine both with and without inclusion of polarization effects, and the resonance energy shifts observed due to polarization are used to predict pi* resonance energies in associated nucleosides and nucleotides, for which static-exchange calculations were carried out. They observe slight shifts between the resonance energies in the isolated bases and those in the nucleosides.

Dissociative electron attachment to abasic DNA

Thin films of the short single DNA strand, GCAT, in which one of the bases has been removed were bombarded with 3 to 15 eV electrons. The yield functions of the H(-), O(-) and OH(-) ions desorbed from these films exhibit a broad peak near 9 eV, which is attributed to dissociative electron attachment to the basic molecules. Whereas removal of any one of the bases considerably decreases N-glycosidic and backbone C-O bond scission, the creation of basic sites does not appreciably modify bond rupture leading to anion electron stimulated desorption. These seemingly contradictory results make it possible to propose a detailed mechanism leading to the transfer of electrons in the range 5-13 eV within DNA.