Site-specific dissociation of DNA bases by slow electrons at early stages of irradiation

Role of excited states in low-energy electron (LEE) induced strand breaks in DNA model systems: influence of aqueous environment

Low-energy electrons (LEEs) interact with a DNA model system to create a variety of excited states. In the gas phase, dissociative (sigma*) states are accessible by LEEs with energy <4 eV (see picture) and cause facile strand breaks through a dissociative electron attachment mechanism. However, under solvation these dissociative (sigma*) states are blue-shifted and are less accessible to LEEs.

Distinct hydroxy-radical-induced damage of 3'-uridine monophosphate in RNA: a theoretical study

RNA strand scission and base release in 3'-uridine monophosphate (UMP), induced by OH radical addition to uracil, is studied at the DFT B3LYP/6-31+G(d,p) level in the gas phase and in solution. In particular, the mechanism of hydrogen-atom transfer subsequent to radical formation, from C2' on the sugar to the C6 site on the base, is explored. The barriers of (C2'-)H2'(a) abstraction by the C6 radical site range from 11.2 to 20.0 kcal mol(-1) in the gas phase and 14.1 to 21.0 kcal mol(-1) in aqueous solution, indicating that the local surrounding governs the hydrogen-abstraction reaction in a stereoselective way. The calculated N1-C1' (N1-glycosidic bond) and beta-phosphate bond strengths show that homolytic and heterolytic bond-breaking processes are largely favored in each case, respectively. The barrier for beta-phosphate bond rupture is approximately 3.2-4.0 kcal mol(-1) and is preferred by 8-12 kcal mol(-1) over N1-glycosidic bond cleavage in both the gas phase and solution. The beta-phosphate bond-rupture reactions are exothermal in the gas phase and solution, whereas N1-C1' bond-rupture reactions require both solvation and thermal corrections at 298 K to be energetically favored. The presence of the ribose 2'-OH group and its formation of low-barrier hydrogen bonds with oxygen atoms of the 3'-phosphate linkage are highly important for hydrogen transfer and the subsequent bond-breakage reactions.

Electronic structure of DNA--unique properties of 8-oxoguanosine

8-Oxo-7,8-dihydroguanosine (8-oxoG) is among the most common forms of oxidative DNA damage found in human cells. The question of damage recognition by the repair machinery is a long standing one, and it is intriguing to suggest that the mechanism of efficiently locating damage within the entire genome might be related to modulations in the electronic properties of lesions compared to regular bases. Using laser-based methods combined with organizing various oligomers self-assembled monolayers on gold substrates, we show that indeed 8-oxoG has special electronic properties. By using oligomers containing 8-oxoG and guanine bases which were inserted in an all thymine sequences, we were able to determine the energy of the HOMO and LUMO states and the relative density of electronic states below the vacuum level. Specifically, it was found that when 8-oxoG is placed in the oligomer, the HOMO state is at higher energy than in the other oligomers studied. In contrast, the weakly mutagenic 8-oxo-7,8-dihydroadenosine (8-oxoA) has little or no effect on the electronic properties of DNA.

Progress in ab initio QM/MM free-energy simulations of electrostatic energies in proteins: accelerated QM/MM studies of pKa, redox reactions and solvation free energies

Hybrid quantum mechanical/molecular mechanical (QM/MM) approaches have been used to provide a general scheme for chemical reactions in proteins. However, such approaches still present a major challenge to computational chemists, not only because of the need for very large computer time in order to evaluate the QM energy but also because of the need for proper computational sampling. This review focuses on the sampling issue in QM/MM evaluations of electrostatic energies in proteins. We chose this example since electrostatic energies play a major role in controlling the function of proteins and are key to the structure-function correlation of biological molecules. Thus, the correct treatment of electrostatics is essential for the accurate simulation of biological systems. Although we will be presenting different types of QM/MM calculations of electrostatic energies (and related properties) here, our focus will be on pKa calculations. This reflects the fact that pKa's of ionizable groups in proteins provide one of the most direct benchmarks for the accuracy of electrostatic models of macromolecules. While pKa calculations by semimacroscopic models have given reasonable results in many cases, existing attempts to perform pKa calculations using QM/MM-FEP have led to discrepancies between calculated and experimental values. In this work, we accelerate our QM/MM calculations using an updated mean charge distribution and a classical reference potential. We examine both a surface residue (Asp3) of the bovine pancreatic trypsin inhibitor and a residue buried in a hydrophobic pocket (Lys102) of the T4-lysozyme mutant. We demonstrate that, by using this approach, we are able to reproduce the relevant side chain pKa's with an accuracy of 3 kcal/mol. This is well within the 7 kcal/mol energy difference observed in studies of enzymatic catalysis, and is thus sufficient accuracy to determine the main contributions to the catalytic energies of enzymes. We also provide an overall perspective of the potential of QM/MM calculations in general evaluations of electrostatic free energies, pointing out that our approach should provide a very powerful and accurate tool to predict the electrostatics of not only solution but also enzymatic reactions, as well as the solvation free energies of even larger systems, such as nucleic acid bases incorporated into DNA.

On the hydrogen loss from protonated nucleobases after electronic excitation or collisional electron capture

In this work, we have subjected protonated nucleobases MH(+) (M = guanine, adenine, thymine, uracil and cytosine) to a range of experiments that involve high-energy (50 keV) collision induced dissociation and electron capture induced dissociation. In the latter case, both neutralisation reionisation and charge reversal were done. For the collision induced dissociation experiments, the ions interacted with O(2). In neutral reionisation, caesium atoms were used as the target gas and the protonated nucleobases captured electrons to give neutrals. These were reionised to cations a microsecond later in collisions with O(2). In choosing Cs as the target gas, we have assured that the first electron transfer process is favourable (by about 0.1-0.8 eV depending on the base). In the case of protonated adenine, charge reversal experiments (two Cs collisions) were also carried out, with the results corroborating those from the neutralisation reionisation experiments. We find that while collisional excitation of protonated nucleobases in O(2) may lead to hydrogen loss with limited probabilities, this channel becomes dominant for electron capture events. Indeed, when sampling reionised neutrals on a microsecond timescale, we see that the ratio between MH(+) and M(+) is 0.2-0.4 when one electron is captured from Cs. There are differences in these ratios between the bases but no obvious correlation with recombination energies was found.

Valence anion of thymine in the DNA pi-stack

Most of theoretical data on the stability of radical anions supported by nucleic acid bases have been obtained for anions of isolated nucleobases, their nucleosides, or nucleotides. This approach ignores the hallmark forces of DNA, namely, hydrogen bonding and pi-stacking interactions. Since these interactions might be crucial for the electron affinities of nucleobases bound in DNA, we report for the first time on the stability of the thymine valence anion in trimers of complementary bases possessing the regular B-DNA geometry but differing in base sequence. In order to estimate the energetics of electron attachment to a trimer, we developed 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 thymine. All calculations were carried out at the MP2 level of theory with the aug-cc-pVDZ basis set. The two-body interaction energies were corrected for the basis set superposition error, and in benchmark systems, they were extrapolated to the basis set limit and supplemented with correction for higher order correlation terms calculated at the CCSD(T) level. We have demonstrated that the sequence of nucleic bases has a profound effect on the stability of the thymine valence anion: the anionic 5'-CTC-3' (6.0 kcal/mol) sequence is the most stable configuration, and the 5'-GTG-3' (-8.0 kcal/mol) trimer anion is the most unstable species. On the basis of obtained results, one can propose DNA sequences that are different in their vulnerability to damage by low energy electron.

The electron affinities of deprotonated adenine, guanine, cytosine, uracil, and thymine

Electron attachment rates and gas phase acidities for the canonical tautomers of the nucleobases and electron affinities for thymine, deprotonated thymine, and cytosine are reported The latter are from a new analysis of published photoelectron spectra. The values for deprotonated thymine are (all in eV) keto-N1-H, 3.327(5); enol-N3-H, 3.250(5), enol-C2OH, 3.120(5) enol-N1-H, 3.013(5), and enol-C4OH,3.123(5). The values for deprotonated cytosine, keto-N1-H, 3.184(5); trans-NH-H, 3.008(5); cis-NH-H, 3.039(5); and enol-N1-H, 2.750(5) and enol-O-H, 2.950(5). The gas phase acidities from these values are obtained from these values using experimental or theoretical calculations of bond dissociation energies. Kinetic and thermodynamic properties for thermal electron attachment to thymine are obtained from mass spectrometric data. We report an activation energy of 0.60 eV and electron affinity of thymine, 1.0(1) eV.

Solvation free energies of molecules. The most stable anionic tautomers of uracil

Anionic states of nucleic acid bases are suspected to play a role in the radiation damage processes of DNA. Our recent studies suggested that the excess electron attachment to the nucleic acid bases can stabilize some rare tautomers, i.e. imine-enamine tautomers and other tautomers with a proton being transferred from nitrogen sites to carbon sites (with respect to the canonical tautomer). So far, these new anionic tautomers have been characterized by the gas-phase electronic structure calculations and photoelectron spectroscopy experiments. In the current contribution we explore the effect of water solvation on the stability of the new anionic tautomers of uracil. The accurate free energies of solvation are calculated in a two step approach. The major contribution was calculated using the classical free-energy perturbation adiabatic-charging approach, where it is assumed that the solvated molecule has the charge distribution given by the polarizable continuum model. In the second step the free energy of solvation is refined by taking into account the real, average solvent charge distribution. This is done using our accelerated QM/MM simulations, where the QM energy of the solute is calculated in the mean potential averaged over many MD steps. We found that in water solution three of the recently identified anionic tautomers are 6.5-3.6 kcal mol(-1) more stable than the anion of the canonical tautomer.

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.

Structure and singly occupied molecular orbital analysis of anionic tautomers of guanine

Recently, we reported the discovery of adiabatically bound anions of guanine that might be involved in the processes of DNA damage by low-energy electrons and in charge transfer through DNA. These anions correspond to some tautomers that have been ignored thus far. They were identified using a hybrid quantum mechanical-combinatorial approach in which an energy-based screening was performed on the library of 499 tautomers with their relative energies calculated with quantum chemistry methods. In the current study, we analyze the adiabatically bound anions of guanine in two aspects: (1) the geometries and excess electron distributions are analyzed and compared with anions of the most stable neutrals to identify the sources of stability; (2) the chemical space of guanine tautomers is explored to verify if these new tautomers are contained in a particular subspace of the tautomeric space. The first task involves the development of novel approaches-the quantum chemical data like electron density, orbital, and information on its bonding/antibonding character are coded into holograms and analyzed using chemoinformatics techniques. The second task is completed using substructure analysis and clustering techniques performed on molecules represented by 2D fingerprints. The major conclusion is that the high stability of adiabatically bound anions originates from the bonding character of the pi orbital occupied by the excess electron. This compensates for the antibonding character that usually causes significant buckling of the ring. Also, the excess electron is more homogenously distributed over both rings than in the case of anions of the most stable neutral species. In terms of 2D substructure, the most stable anionic tautomers generally have additional hydrogen atoms at C8 and/or C2 and they do not have hydrogen atoms attached to C4, C5, and C6. They also form an "island of stability" in the tautomeric space of guanine.

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.

Dissociative electron attachment to gas phase glycine: exploring the decomposition pathways by mass separation of isobaric fragment anions

Dissociative electron attachment to gas phase glycine generates a number of fragment ions, among them ions observed at the mass numbers 15, 16 and 26 amu. From stoichiometry they can be assigned to the chemically rather different species NH(-)/CH(3)(-)(15 amu), O(-)/NH(2)(-)(16 amu) and CN(-)/C(2)H(2)(-)(26 amu). Here we use a high resolution double focusing two sector mass spectrometer to separate these isobaric ions. It is thereby possible to unravel the decomposition reactions of the different transient negative ions formed upon resonant electron attachment to neutral glycine in the energy range 0-15 eV. We find that within the isobaric ion pairs, the individual components generally arise from resonances located at substantial different energies. The corresponding unimolecular decompositions involve complex reaction sequences including multiple bond cleavages and substantial rearrangement in the precursor ion. To support the interpretation and assignments we also use (13)C labelling of glycine at the carboxylic group.

Interaction of low-energy electrons with the pyrimidine bases and nucleosides of DNA

We report computed cross sections for the elastic scattering of slow electrons by the pyrimidine bases of DNA, thymine and cytosine, and by the associated nucleosides, deoxythymidine and deoxycytidine. For the isolated bases, we carried out calculations both with and without the inclusion of polarization effects. For the nucleosides, we neglect polarization effects but estimate their influence on resonance positions by comparison with the results for the corresponding bases. Where possible, we compare our results with experiment and previous calculations.

Dissociative electron attachment to furan, tetrahydrofuran, and fructose

We study dissociative electron attachment to furan (FN) (C(4)H(4)O), tetrahydrofuran (THF) (C(4)H(8)O), and fructose (FRU) (C(6)H(12)O(6)) using crossed electron/molecular beams experiments with mass spectrometric detection of the anions. We find that FN and THF are weak electron scavengers and subjected to dissociative electron attachment essentially in the energy range above 5.5 eV via core excited resonances. In striking contrast to that, FRU is very sensitive towards low energy electrons generating a variety of fragment ions via a pronounced low energy feature close to 0 eV. These reactions are associated with the degradation of the ring structure and demonstrate that THF cannot be used as surrogate to model deoxyribose in DNA with respect to the attack of electrons at subexcitation energies (<3 eV). The results support the picture that in DNA the sugar moiety itself is an active part in the initial molecular processes leading to single strand breaks.

Phosphodiester and N-glycosidic bond cleavage in DNA induced by 4-15 eV electrons

Thin molecular films of the short single strand of DNA, GCAT, were bombarded under vacuum by electrons with energies between 4 and 15 eV. Ex vacuo analysis by high-pressure liquid chromatography of the samples exposed to the electron beam revealed the formation of a multitude of products. Among these, 12 fragments of GCAT were identified by comparison with reference compounds and their yields were measured as a function of electron energy. For all energies, scission of the backbone gave nonmodified fragments containing a terminal phosphate, with negligible amounts of fragments without the phosphate group. This indicates that phosphodiester bond cleavage by 4-15 eV electrons involves cleavage of the C-O bond rather than the P-O bond. The yield functions exhibit maxima at 6 and 10-12 eV, which are interpreted as due to the formation of transient anions leading to fragmentation. Below 15 eV, these resonances dominate bond dissociation processes. All four nonmodified bases are released from the tetramer, by cleavage of the N-glycosidic bond, which occurs principally via the formation of core-excited resonances located around 6 and 10 eV. The formation of the other nonmodified products leading to cleavage of the phosphodiester bond is suggested to occur principally via two different mechanisms: (1) the formation of a core-excited resonance on the phosphate unit followed by dissociation of the transient anion and (2) dissociation of the CO bond of the phosphate group formed by resonance electron transfer from the bases. In each case, phosphodiester bond cleavage leads chiefly to the formation of stable phosphate anions and sugar radicals with minimal amounts of alkoxyl anions and phosphoryl radicals.

Identification of Isomers of Nitrotoluene via Free Electron Attachment

Free electron attachment to the three different isomers of mononitrotoluene molecules in the gas phase is studied using a crossed electron-molecule beams technique. In contrast to previous studies for a large number of negative ions, the presently measured relative cross section curves are recorded with an electron energy resolution of better than 100 meV. For several product anions including the nitro anion NO(2)-, remarkable differences for the three isomers are observed. In almost all fragment anion efficiency curves, the 2-nitrotoluene exhibits pronounced differences from the two other isomers. In contrast, 3- and 4-nitrotoluene disagree only slightly in a few fragment anions from each other.

Vibrational Feshbach resonances in uracil and thymine

Sharp peaks in the dissociative electron attachment (DEA) cross sections of uracil and thymine at energies below 3 eV are assigned to vibrational Feshbach resonances (VFRs) arising from coupling between the dipole bound state and the temporary anion state associated with occupation of the lowest sigma* orbital. Three distinct vibrational modes are identified, and their presence as VFRs is consistent with the amplitudes and bonding characteristics of the sigma* orbital wave function. A deconvolution method is also employed to yield higher effective energy resolution in the DEA spectra. The site dependence of DEA cross sections is evaluated using methyl substituted uracil and thymine to block H atom loss selectively. Implications for the broader issue of DNA damage are briefly discussed.

Ion-induced biomolecular radiation damage: from isolated nucleobases to nucleobase clusters

A large number of studies are devoted to the investigation of the biomolecular ionization and fragmentation dynamics underlying biological radiation damage. Most of these studies have been based on gas-phase collisions with isolated DNA building blocks. The radiobiological significance of these studies is often questioned because of the lack of a chemical environment. To clarify this aspect, we studied interactions of keV ions with isolated nucleobases and with nucleobase clusters by means of coincidence time-of-flight spectrometry. Significant changes already show up in the molecular fragmentation patterns of very small clusters.

Bond selective dissociative electron attachment to thymine

Free-electron attachment to thymine and partially deuterated thymine, where D replaces H at all carbon atoms, is studied in the electron energy range from about 0 to 15 eV. The formation of fragment anions that are formed by the loss of one or two H (D) atoms is analyzed as a function of the incident electron energy using a crossed electron/neutral beam apparatus in combination with a quadrupole mass spectrometer. By using partially deuterated thymine and quantum-chemical calculation a bond selectivity for the loss of one and two hydrogen atoms is observed that is determined only by the kinetic energy of the incident electron.

Hydrogen-Abstracted Adenine−Thymine Radicals with Interesting Transferable Properties

The formation of radicals on DNA bases through various pathways can lead to harmful structural alterations. Such processes are of interest for preventing alteration of healthy DNA and, conversely, to develop more refined methods for inhibiting the replication of unwanted mutagenic DNA. In the present work, we explore theoretically the energetic and structural properties of the nine possible neutral radicals formed via hydrogen abstraction from the adenine-thymine base pair. The lowest energy radical is formed by loss of a hydrogen atom from the methyl group of thymine. The next lowest energy radicals, lying 8 and 9 kcal mol-1 higher than the global minimum, are those in which hydrogens are removed from the two nitrogens that would join the base pair to 2-deoxyribose in double-stranded DNA. The other six radicals lie between 16 and 32 kcal mol-1 higher in energy. Unlike the guanine-cytosine base pair, adenine-thymine (A-T) exhibits only minor structural changes upon hydrogen abstraction, with all A-T derived radicals maintaining planarity. Moreover, the energetic ordering for the radicals of the two isolated bases (adenine and thymine) is preserved upon formation of the base pair, though with a wider spread of energies. Even more significantly, the energetic interleaving of the (A-H)*-T and A-(T-H)* radicals is correctly predicted from the X-H bond dissociation energies of the isolated adenine and thymine. This suggests that the addition of the hydrogen-bonded complement base only marginally affects the bond energies.