The Adenine−Thymine Base Pair Radical Anion: Adding an Electron Results in a Major Structural Change

Understanding electron attachment to the DNA double helix: the thymidine monophosphate-adenine pair in the gas phase and aqueous solution

Electron attachment to the 2'-deoxythymidine-5'-monophosphate-adenine pairs (5'-dTMPH-A and 5'-dTMP(-)-A) has been investigated at a carefully calibrated level of theory (B3LYP/DZP++) to investigate the electron-accepting properties of thymine (T) in the DNA double helix under physiological conditions. All molecular structures have been fully optimized in vacuo and in solution. The adiabatic electron affinity of 5'-dTMPH-A in the gas phase has been predicted to be 0.67 eV. Solvent effects greatly increase the electron capture ability of 5'-dTMPH-A. In fact, the adiabatic electron affinity increases to 2.04 eV with solvation. The influence of the solvent environment on the electron-attracting properties of 5'-dTMPH-A arises not only from the stabilization of the corresponding radical anion through charge-dipole interactions, but also by changing the distribution of the unpaired electron in the molecular system. The unpaired electron is covalently bound even during vertical attachment, due to the solvent effects. Solvent effects also weaken the pairing interaction in the thymidine monophosphate-adenine complexes. The phosphate deprotonation is found to have a relatively minor influence on the capture of electrons by the 5'-dTMPH-A species in aqueous solution. The electron distributions, natural population analysis, and geometrical features of the models examined illustrate that the influence of the phosphate deprotonation is limited to the phosphate moiety in aqueous solution. Therefore, it is reasonable to expect that electron attachment to nucleotides will be independent of monovalent counterions in the vicinity of the phosphate group in aqueous solution.

Vibrational dynamics of DNA. I. Vibrational basis modes and couplings

Carrying out density functional theory calculations of four DNA bases, base derivatives, Watson-Crick (WC) base pairs, and multiple-layer base pair stacks, we studied vibrational dynamics of delocalized modes with frequency ranging from 1400 to 1800 cm(-1). These modes have been found to be highly sensitive to structure fluctuation and base pair conformation of DNA. By identifying eight fundamental basis modes, it is shown that the normal modes of base pairs and multilayer base pair stacks can be described by linear combinations of these vibrational basis modes. By using the Hessian matrix reconstruction method, vibrational coupling constants between the basis modes are determined for WC base pairs and multilayer systems and are found to be most strongly affected by the hydrogen bonding interaction between bases. It is also found that the propeller twist and buckle motions do not strongly affect vibrational couplings and basis mode frequencies. Numerically simulated IR spectra of guanine-cytosine and adenine-thymine bases pairs as well as of multilayer base pair stacks are presented and described in terms of coupled basis modes. It turns out that, due to the small interlayer base-base vibrational interactions, the IR absorption spectrum of multilayer base pair system does not strongly depend on the number of base pairs.

The 2'-deoxyadenosine-5'-phosphate anion, the analogous radical, and the different hydrogen-abstracted radical anions: molecular structures and effects on DNA damage

The 2'-deoxyadenosine-5'-phosphate (5'-dAMP) anion and its related radicals have been studied by reliably calibrated theoretical approaches. This study reveals important physical characteristics of 5'-dAMP radical related processes. One-electron oxidation of the 5'-dAMP anion is found on both the phosphoryl group and the adenine base with electron detachment energies close to that of phosphate. Partial removal of electron density from the adenine fragment leads to an extended pi system which includes the amine group of the adenine. Although the radical-centered carbon increases the extent of bonding with its adjacent atoms, it usually weakens the chemical bonds between the atoms at the alpha- and beta-positions. This tendency should be important in predicting the reactivity of the sugar-based radicals. The overall stability sequence of the H-abstracted 5'-dAMP anionic radicals is consistent with the analogous results for the H-abstracted neutral radicals of the adenosine nucleoside: aliphatic radicals > aromatic radicals. The negatively charged phosphoryl group attached to atom C(5)' of the ribose does not change this energetic sequence. All the H-abstraction produced 5'-dAMP radical anions are distonic radical anions. Studies have shown that the charge-radical-separating feature of the distonic radical anions is biologically relevant. This result should be important in understanding the reactive properties of these H-abstraction-produced anion radicals.

Non-Watson-Crick base pairing in RNA. quantum chemical analysis of the cis Watson-Crick/sugar edge base pair family

Large RNA molecules exhibit an astonishing variability of base-pairing patterns, while many of the RNA base-pairing families have no counterparts in DNA. The cis Watson-Crick/sugar edge (cis WC/SE) RNA base pairing is investigated by ab initio quantum chemical calculations. A detailed structural and energetic characterization of all 13 crystallographically detected members of this family is provided by means of B3LYP/6-31G and RIMP2/aug-cc-pVDZ calculations. Further, a prediction is made for the remaining 3 cis WC/SE base pairs which are yet to be seen in the experiments. The interaction energy calculations point at the key role of the 2'-OH group in stabilizing the sugar-base contact and predict all 16 cis WC/SE base-pairing patterns to be nearly isoenergetic. The perfect correlation of the main geometrical parameters in the gas-phase optimized and X-ray structures shows that the principle of isosteric substitutions in RNA is rooted from the intrinsic structural similarity of the isolated base pairs. The present quantum chemical calculations for the first time analyze base pairs involving the ribose 2'-OH group and unambiguously correlate the structural information known from experiments with the energetics of interactions. The calculations further show that the relative importance and absolute value of the dispersion energy in the cis WC/SE base pairs are enhanced compared to the standard base pairs. This may by an important factor contributing to the strength of such interactions when RNA folds in its polar environment. The calculations further demonstrate that the Cornell et al. force field commonly used in molecular modeling and simulations provides satisfactory performance for this type of RNA interactions.

Hydrogen-bonding interactions in peptide nucleic acid and deoxyribonucleic acid: a comparative study

Peptide nucleic acid (PNA) is a synthetic analogue of deoxyribonucleic acid (DNA) capable of tightly binding to itself and DNA with high specificity. Using hybrid density functional methods, hydrogen-bond (H-bond) strengths have been evaluated for isolated Watson-Crick base pairs, PNA base pairs, and charged as well as neutral DNA base pairs. Heterogeneous base pairs of PNA with charged and neutral DNA have also been investigated. The competing effects of short-range H-bonding and long-range Coulombic repulsions in charged DNA base pairs have been analyzed. Polarizable continuum models have been employed to evaluate solvation effects on the binding energies.

Development of an ONIOM-G3B3 method to accurately predict C-H and N-H bond dissociation enthalpies of ribonucleosides and deoxyribonucleosides

The roles of ribonucleoside and deoxyribonucleoside radicals in DNA and RNA damage cannot be properly understood in the absence of knowledge of the C-H and N-H bond dissociation enthalpies (BDEs) depicting the energy cost to generate each of these radicals. However, because the nucleoside radicals tend to be extremely short-lived and it is very difficult to separate and identify different nucleoside radicals, experimental BDEs for nucleosides have remained elusive. Herein, we developed an ONIOM-G3B3 method in order to reliably predict the BDEs of nucleosides and we carefully benchmarked this new method against over 60 experimental BDEs of diverse sizable molecules. It was found that the accuracy of the ONIOM-G3B3 method was about 1.4 kcal/mol for BDE calculations. Using the ONIOM-G3B3 method, a full scale of C-H and N-H BDEs were obtained for the first time for ribonucleosides and deoxyribonucleosides with an estimated error bar of +/-1.4 kcal/mol. Discussions were then made about the interesting connections between these BDE values and previously reported experimental observations concerning radical-mediated DNA and RNA lesions. The significance of the work is twofold: (i) Nucleosides represent one of the most important groups of compounds in science. A full scale of reliable bond dissociation enthalpies for nucleosides is of fundamental importance. (ii) This work demonstrates the feasibility to accurately predict the bond strength of various sizable molecules ranging from nanosize molecular devices to biologically significant compounds.

MS-CASPT2 calculation of excess electron transfer in stacked DNA nucleobases

Calculations using the complete active space self-consistent field (CASSCF) and complete active space second-order perturbation (CASPT2) methods, and the multistate formulation of CASPT2 (MS-CASPT2), are performed for the ground and excited states of radical anions consisting of two pi-stacked nucleobases. The electronic couplings for excess electron transfer (EET) in the pi-stacks are estimated by using the generalized Mulliken-Hush approach. We compare results obtained within the different methods with data derived using Koopmans' theorem approximation at the Hartree-Fock level. The results suggest that although the one-electron scheme cannot be applied to calculate electron affinities of nucleobases, it provides reasonable estimates for EET energies. The electronic couplings calculated with KTA lie between the CASPT2 and the MS-CASPT2 based values in almost all cases.

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.

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.

Vibrational dynamics of DNA: IV. Vibrational spectroscopic characteristics of A-, B-, and Z-form DNA’s

Linear and nonlinear IR spectroscopic studies of nucleic acids can provide crucial information on solution conformations of DNA double helix and its complex with other molecules. Carrying out density functional theory calculations of A-, B-, and Z-form DNA's, the authors obtained vibrational spectroscopic properties as well as coupling constants between different basis modes. The vibrational couplings that determine the extent of exciton delocalization are strongly dependent on DNA conformation mainly because the interlayer distance between two neighboring base pairs changes with respect to the DNA conformation. The Z-DNA has comparatively small interlayer vibrational coupling constants so that its vibrational spectrum depends little on the number of base pairs, whereas the A-DNA shows a notable dependency on the size. Furthermore, it is shown that a few distinctively different line shape changes in both IR and two-dimensional IR spectra as the DNA conformation changes from B to A or from B to Z can be used as marker bands and characteristic features distinguishing different DNA conformations.

Estimates of electronic coupling for excess electron transfer in DNA

Electronic coupling V(da) is one of the key parameters that determine the rate of charge transfer through DNA. While there have been several computational studies of V(da) for hole transfer, estimates of electronic couplings for excess electron transfer (ET) in DNA remain unavailable. In the paper, an efficient strategy is established for calculating the ET matrix elements between base pairs in a pi stack. Two approaches are considered. First, we employ the diabatic-state (DS) method in which donor and acceptor are represented with radical anions of the canonical base pairs adenine-thymine (AT) and guanine-cytosine (GC). In this approach, similar values of V(da) are obtained with the standard 6-31G(*) and extended 6-31+ +G(**) basis sets. Second, the electronic couplings are derived from lowest unoccupied molecular orbitals (LUMOs) of neutral systems by using the generalized Mulliken-Hush or fragment charge methods. Because the radical-anion states of AT and GC are well reproduced by LUMOs of the neutral base pairs calculated without diffuse functions, the estimated values of V(da) are in good agreement with the couplings obtained for radical-anion states using the DS method. However, when the calculation of a neutral stack is carried out with diffuse functions, LUMOs of the system exhibit the dipole-bound character and cannot be used for estimating electronic couplings. Our calculations suggest that the ET matrix elements V(da) for models containing intrastrand thymine and cytosine bases are essentially larger than the couplings in complexes with interstrand pyrimidine bases. The matrix elements for excess electron transfer are found to be considerably smaller than the corresponding values for hole transfer and to be very responsive to structural changes in a DNA stack.

Does the A·T or G·C Base-Pair Possess Enhanced Stability? Quantifying the Effects of CH···O Interactions and Secondary Interactions on Base-Pair Stability Using a Phenomenological Analysis and ab Initio Calculations

An empirically based relationship between overall complex stability (-DeltaG degrees ) and various possible component interactions is developed to probe the question of whether the A.T/U and G.C base-pairs exhibit enhanced stability relative to similarly hydrogen-bonded complexes. This phenomenological approach suggests ca. 2-2.5 kcal mol-1 in additional stability for A.T owing to a group interaction containing a CH...O contact. Pairing geometry and the role of the CH...O interaction in the A.T base-pair were also probed using MP2/6-31+G(d,p) calculations and a double mutant cycle. The ab initio studies indicated that Hoogsteen geometry is preferred over Watson-Crick geometry in A.T by ca. 1 kcal mol-1. Factors that might contribute to the preference for Hoogsteen geometry are a shorter CH...O contact, a favorable alignment of dipoles, and greater distances between secondary repulsive sites. The CH...O interaction was also investigated in model complexes of adenine with ketene and isocyanic acid. The ab initio calculations support the result of the phenomenological approach that the A.T base-pair does have enhanced stability relative to hydrogen-bonded complexes with just N-H...N and N-H...O hydrogen bonds.

Hydrogen loss from nucleobase nitrogens upon electron attachment to isolated DNA and RNA nucleotide anions

Electron transfer to isolated nucleotide monoanions in collisions with Na vapor induces hydrogen loss from nitrogen of the transient nucleobase anion. The cross section for this process is linearly correlated with the number of N-H hydrogens and is highest for guanine. The process is much faster than microseconds since only dehydrogenated dianions survived for mass spectrometric detection. The lifetime of the adenosine 5(')-monophospate dianions was measured to be 0.2 ms in an electrostatic ion storage ring but also a longer-lived component with a lifetime of at least 10 ms was identified. Implications of dissociation along the N-H coordinate for a nucleotide in DNA are briefly discussed in terms of Watson-Crick base pairs.

Double proton transfer in the isolated and DNA-embedded guanine-cytosine base pair

The energetics and dynamics of double proton transfer (DPT) is investigated theoretically for the Watson-Crick conformation of the guanine-cytosine (GC) base pair. Using semiempirical density functional theory the isolated and DNA-embedded GC pair is considered. Differences in the energetics and dynamics of DPT thus addresses the question of how relevant studies of isolated base pairs are for the understanding of processes occurring in DNA. Two-dimensional potential energy surfaces involving the transferring hydrogen atoms and the proton donors and acceptors are presented for both systems. The DPT reaction is accompanied by a contraction of the distance between the two bases with virtually identical energetic barriers being 18.8 and 18.7 kcal/mol for the isolated and DNA-embedded system, respectively. However, the transition state for DPT in the DNA-embedded GC pair is offset by 0.1 A to larger N-H separation compared to the isolated GC pair. Using activated ab initio molecular dynamics, DPT is readily observed for the isolated base pair with a minimal amount of 21.4 kcal/mol of initial average kinetic energy along the DPT normal mode vector. On a time scale of approximately 100 fs DPT has occurred and the excess energy is redistributed. For the DNA-embedded GC pair considerably more kinetic energy is required (30.0 kcal/mol) for DPT and the process is completed within one hydrogen vibration. The relevance of studies of isolated base pairs and base pair analogs in regard of reactions or properties involving DNA is discussed.

Localized Orbital Corrections for the Calculation of Ionization Potentials and Electron Affinities in Density Functional Theory

This paper describes the extension of a previously reported empirical localized orbital correction model to the correction of ionization potential energies (IP) and electron affinities (EA) for atoms and molecules of first and second row elements. The B3LYP localized orbital correction version of the model (B3LYP-LOC) uses 22 heuristically determined parameters that improve B3LYP DFT IP and EA energy calculations on the G2 data set of 134 molecules from a mean absolute deviation (MAD) from experiment of 0.137 to 0.039 eV. The method significantly reduces the number of outliers and overall MAD to error levels below that achieved with G2 wave function based theory; furthermore, the new model has zero additional computational cost beyond standard DFT calculations. Although the model is heuristic and is based on a multiple linear regression to experimental errors, each of the parameters is justified on physical grounds, and each provides insight into the fundamental limitations of DFT, most importantly the failure of current DFT methods to accurately account for nondynamical electron correlation.

Deprotonated Cytosine Anions: A Theoretical Prediction of Photoelectron Spectra

Predictions on the photoelectron spectra of deprotonated cytosine anions (cytosinate, Cye(-)) have been made with ab initio electron propagator methods. Two imino-oxo forms are most stable, but four other isomers have energies within 10 kcal/mol. The first vertical electron detachment energies (VEDEs) for the three most stable Cye(-) isomers are approximately 3.4 eV. Imino-oxy VEDEs are about 0.3 eV smaller. For each anion, the lowest VEDE corresponds to a pi Dyson orbital. The order of higher final states is changed when relaxation and correlation effects are considered. Considerable mixing between lone-pair and bonding lobes occurs in the sigma Dyson orbitals.

Hydrogen Atom and Hydride Anion Addition to Adenine: Structures and Energetics

The radicals and anions derived from the 9H tautomer of adenine by adding a hydrogen atom to one of the four double bonds of the adenine framework have been studied. Computations were carried out using a carefully calibrated density functional (B3LYP) method and basis set (DZP++). Optimized geometries, energies, and vibrational frequencies are predicted for eight radicals and anions. The radicals are found to lie in a range of 22 kcal mol(-1), with the radical derived by addition to the C(8) carbon atom being the lowest lying energetically. The anions are predicted to be bound species in the gas phase with an energetic range of 43 kcal mol(-1). Anions produced by addition of a hydride ion to adenine carbon atoms are found to be the most favorable. Six of the anions are predicted to be stable species with respect to electron detachment. The adiabatic electron affinities, vertical electron affinities, and vertical detachment energies are computed for the first time. Electron affinities for these radicals range from 0.0 to 2.0 eV. Radicals produced by addition to a nitrogen atom have near-zero adiabatic electron affinities, while radicals produced by addition at carbon atoms have considerably higher electron affinities.

Microsolvation effects on the electron capturing ability of thymine: Thymine-water clusters

The effects of solvation on the stability of thymine and its negative ion have been investigated by explicitly considering the structures of complexes of thymine with up to five water molecules and the respective anions at the B3LYP/DZP++ level of theory. The vertical detachment energy of thymine was predicted to increase gradually with the hydration number, consistent with experimental observations from a photodetachment-photoelectron spectroscopy study J. Schiedt et al., [Chem. Phys. 239, 511 (1998)]. The adiabatic electron affinity of thymine was also found to increase with the hydration number, which implies that while the conventional valence anion of thymine is only marginally bound in the gas phase, it may form a stable anion in aqueous solution.

Binding of Gold Clusters with DNA Base Pairs: A Density Functional Study of Neutral and Anionic GC−Aun and AT−Aun (n = 4, 8) Complexes

Binding of clusters of gold atoms (Au) with the guanine-cytosine (GC) and adenine-thymine (AT) Watson-Crick DNA base pairs was studied using the density functional theory (DFT). Geometries of the neutral GC-Au(n) and AT-Au(n) and the corresponding anionic (GC-Au(n))(-1) and (AT-Au(n))(-1) (n = 4, 8) complexes were fully optimized in different electronic states, that is, singlet and triplet states for the neutral complexes and doublet and quartet states for the anionic complexes, using the B3LYP density functional method. The 6-31+G basis set was used for all atoms except gold. For gold atoms, the Los Alamos effective core potential (ECP) basis set LanL2DZ was employed. Vibrational frequency calculations were performed to ensure that the optimized structures corresponded to potential energy surface minima. The gold clusters around the neutral GC and AT base pairs have a T-shaped structure, which satisfactorily resemble those observed experimentally and in other theoretical studies. However, in anionic GC and AT base pairs, the gold clusters have extended zigzag and T-shaped structures. We found that guanine and adenine have high affinity for Au clusters, with their N3 and N7 sites being preferentially involved in binding with the same. The calculated adiabatic electron affinities (AEAs) of the GC-Au(n)complexes (n = 4, 8) were found to be much larger than those of the isolated base pairs.

The deprotonated guanine-cytosine base pair

Awareness of the harmful effects of radiation has increased interest in finding the mechanisms of DNA damage. Radical and anion formation among the DNA base pairs are thought to be important steps in such damage [Collins, G. P. (2003) Sci. Am. 289 (3), 26-27]. Energetic properties and optimized geometries of 10 radicals and their respective anions derived through hydrogen abstraction from the Watson-Crick guanine-cytosine (G-C) base pair have been studied using reliable theoretical methods. The most favorable deprotonated structure (dissociation energy 42 kcal x mol(-1), vertical detachment energy 3.79 eV) ejects the proton analogous to the cytosine glycosidic bond in DNA. This structure is a surprisingly large 12 kcal x mol(-1) lower in energy than any of the other nine deprotonated G-C structures. This system retains the qualitative G-C structure but with the H...O2 distance dramatically reduced from 1.88 to 1.58 A, an extremely short hydrogen bond. The most interesting deprotonated G-C structure is a "reverse wobble" incorporating two N-H...N hydrogen bonds. Three different types of relaxation energies (4.3-54 kcal x mol(-1)) are defined and reported to evaluate the energy released via different mechanisms for the preparation of the deprotonated species. Relative energies, adiabatic electron affinities (ranging from 1.93 to 3.65 eV), and pairing energies are determined to discern which radical will most alter the G-C properties. The most stable deprotonated base pair corresponds to the radical with the largest adiabatic electron affinity, 3.65 eV. This value is an enormous increase over the electron affinity (0.60 eV) of the closed-shell G-C base pair.

Microsolvation Effect, Hydrogen-Bonding Pattern, and Electron Affinity of the Uracil−Water Complexes U−(H2O)n (n = 1, 2, 3)

To achieve a systematic understanding of the influence of microsolvation on the electron accepting behaviors of nucleobases, the reliable theoretical method (B3LYP/DZP++) has been applied to a comprehensive conformational investigation on the uracil-water complexes U-(H(2)O)(n) (n = 1, 2, 3) in both neutral and anionic forms. For the neutral complexes, the conformers of hydration on the O2 of uracil are energetically favored. However, hydration on the O4 atom of uracil is more stable for the radical anions. The electron structure analysis for the H-bonding patterns reveal that the CH...OH(2) type H-bond exists only for di- and trihydrated uracil complexes in which a water dimer or trimer is involved. The electron density structure analysis and the atoms-in-molecules (AIM) analysis for U-(H(2)O)(n) suggest a threshold value of the bond critical point (BCP) density to justify the CH...OH(2) type H-bond; that is, CH...OH(2) could be considered to be a H-bond only when its BCP density value is equal to or larger than 0.010 au. The positive adiabatic electron affinity (AEA) and vertical detachment energy (VDE) values for the uracil-water complexes suggest that these hydrated uracil anions are stable. Moreover, the average AEA and VDE of U-(H(2)O)(n) increase as the number of the hydration waters increases.

Structural and Energetic Characterization of a DNA Nucleoside Pair and Its Anion: Deoxyriboadenosine (dA) − Deoxyribothymidine (dT)

The geometries of the DNA nucleoside pairs between 2'-deoxyriboadenosine (dA) and 2'-deoxyribothymidine (dT) and its anion (dAdT-) were fully optimized using carefully calibrated density functional methods. The addition of an electron to dAdT results in remarkable changes to the two hydrogen bonding distances, the H...O distance decreasing by 0.303 angstroms and the N...H distance increasing by 0.229 angstroms. The electron affinity of the dAdT pair was studied to reveal the correct trends of adiabatic electron affinity (EA(ad)) under the influence of the additional components to the individual bases. The consequence of negative charge in terms of structural variations, energetic changes, and charge distribution were explored. The EA(ad) of dAdT is predicted to be positive (0.60 eV), and it exhibits a substantial increase compared with those of the corresponding bases A and T and the nucleic acid base pair AT. The effects of pairing and the addition of the sugar moiety on the EA(ad) are well described as the summation of the individual influences. The influence of the pairing on the EA is comparable to that of the addition of 2-deoxyribose. The excess charge is mainly located on the thyminyl moiety in the anionic dAdT pair. The positive vertical electron affinity (VEA = 0.20 eV) for dAdT suggests that it is able to form a stable anion through electron attachment. A large vertical detachment energy (VDE = 1.14 eV) has been determined for the anionic dAdT nucleoside pair. Therefore, one may expect that the stable anionic dAdT nucleoside pair should be able to undergo the subsequent glycosidic bond cleavage process.

Marked Influences on the Adenine−Cytosine Base Pairs by Electron Attachment and Ionization

The reverse wobble and the reverse Hoogsteen adenine-cytosine mispairs regarding their radical cations and anions are studied with the hybrid three-parameter B3LYP density functional method and 6-31+G(d), 6-311+G(2df,2p) basis sets. Hydrogen bonding mispairs are remarkably influenced by electron attachment and ionization. Only one stronger hydrogen bond N6-H (in adenine)...N3 (in cytosine) exists in the radical pair, while the strengths of two N-H...N hydrogen bonds in the neutral pair are comparable. Geometrical coplanarity is found for the neutral and cationic pairs, in contrast to the anionic pairs in which the cytosine moiety exhibits significant deformation due to electron attachment. Dissociation energies for the neutral and radical pairs are slightly higher than those of the adenine-thymine pairs but much smaller than those of the guanine-cytosine pairs. Valence-bound anions of these two adenine-cytosine pairs are thermodynamically stable by 0.1-0.2 eV with respect to the neutral pairs. On the basis of the comparison between the experimental data of the solvated clusters and the calculated values, these two pairs can be quantitatively equivalent to the clusters in which each base is solvated by five water molecules.

Adiabatic Electron Affinities of the Polyhydrated Adenine−Thymine Base Pair: A Density Functional Study

Adiabatic electron affinities (AEAs) of the adenine-thymine (AT) base pair surrounded by 5 and 13 water molecules have been studied by density functional theory (DFT). Geometries of neutral AT x nH2O and anionic (AT x nH2O)- complexes (n = 5 and 13) were fully optimized, and vibrational frequency analysis was performed at the B3LYP/6-31+G** level of theory. The optimized structures of the neutral (AT x nH2O) and (AT x nH2O)- complexes were found to be somewhat nonplanar. Some of the water molecules are displaced away from the AT ring plane and linked with one another by hydrogen bonds. The optimized structures of the complexes are found to be in a satisfactory agreement with the observed experimental and molecular dynamics simulation results. In the optimized anionic complexes, the thymine (T) moiety was found to be puckered, whereas the adenine (A) moiety remained almost planar. Natural population analysis (NPA) performed using the B3LYP/6-31+G** method shows that the thymine moiety in the anionic (AT x nH2O)- complexes (n = 5 and 13) has most of the excess electronic charge, i.e., approximately -0.87 and approximately -0.81 (in the unit of magnitude of the electronic charge), respectively. The zero-point energy corrected adiabatic electron affinities of the hydrated AT base pair were found to be positive both for n = 5 and 13 and have the values of 0.97 and 0.92 eV, respectively, which are almost three times the AEA of the AT base pair. The results show that the presence of water molecules appreciably enhances the EA of the base pair.

AT Base Pair Anions versus (9-Methyl-A)(1-Methyl-T) Base Pair Anions

The anionic base pairs of adenine and thymine, (AT)(-), and 9-methyladenine and 1-methylthymine, (MAMT)(-), have been investigated both theoretically and experimentally in a complementary, synergistic study. Calculations on (AT)(-) found that it had undergone a barrier-free proton transfer (BFPT) similar to that seen in other dimer anion systems and that its structural configuration was neither Watson-Crick (WC) nor Hoogsteen (HS). The vertical detachment energy (VDE) of (AT)(-) was determined by anion photoelectron spectroscopy and found to be in agreement with the VDE value predicted by theory for the BFPT mechanism. An AT pair in DNA is structurally immobilized into the WC configuration, in part, by being bonded to the sugars of the double helix. This circumstance was mimicked by methylating the sites on both A and T where these sugars would have been tied, viz., 9-methyladenine and 1-methylthymine. Calculations found no BFPT in (MAMT)(-) and a resulting (MAMT)(-) configuration that was either HS or WC, with the configurations differing in stability by ca. 2 kcal/mol. The photoelectron spectrum of (MAMT)(-) occurred at a completely different electron binding energy than had (AT)(-). Moreover, the VDE value of (MAMT)(-) was in agreement with that predicted by theory. The configuration of (MAMT)(-) and its lack of electron-induced proton transfer are inter-related. While there may be other pathways for electron-induced DNA alterations, BFPT in the WC/HS configurations of (AT)(-) is not feasible.

(G-H)*-C and G-(C-H)* radicals derived from the guanine.cytosine base pair cause DNA subunit lesions

The radicals generated by the homolytic cleavage of an X-H bond from the guanine.cytosine (G.C) base pair were studied by using carefully calibrated theoretical methods. The gradient-corrected density functional B3LYP was applied in conjunction with double-zeta plus polarization and diffuse function basis sets. Optimized geometries, energies, and vibrational frequencies were obtained for all of the radicals considered. Structural perturbations along with energy relaxation due to radical formation were investigated. Dissociation energies of the G.C base pair and all of the radicals are predicted and compared with the dissociation energy of neutral G.C. The three lowest-energy base pair radicals all involve removal of an H atom from one of the N atoms in G.C. The lowest-energy base pair radical has the hydrogen atom removed from the guanine nitrogen atom used for the sugar phosphate linkage in DNA. This (G-H)(*)-C radical has a dissociation energy (to G-H(*) + C) of 30 kcal/mol, which may be compared with 27 kcal/mol for G.C. All of the radicals that are possible outcomes of direct ionizing radiation or oxidizing species were investigated for the presence of local minima with significant structural changes. Major structural deformations cause strain in the interstrand hydrogen bonding in the DNA double helix. Severe geometry changes were observed when the hydrogen was abstracted from interstrand hydrogen bonding sites, along with sizeable energy changes, indicating the potentially serious consequences to the G.C base pair.

Accurate interaction energies of hydrogen-bonded nucleic acid base pairs

Hydrogen-bonded nucleic acids base pairs substantially contribute to the structure and stability of nucleic acids. The study presents reference ab initio structures and interaction energies of selected base pairs with binding energies ranging from -5 to -47 kcal/mol. The molecular structures are obtained using the RI-MP2 (resolution of identity MP2) method with extended cc-pVTZ basis set of atomic orbitals. The RI-MP2 method provides results essentially identical with the standard MP2 method. The interaction energies are calculated using the Complete Basis Set (CBS) extrapolation at the RI-MP2 level. For some base pairs, Coupled-Cluster corrections with inclusion of noniterative triple contributions (CCSD(T)) are given. The calculations are compared with selected medium quality methods. The PW91 DFT functional with the 6-31G basis set matches well the RI-MP2/CBS absolute interaction energies and reproduces the relative values of base pairing energies with a maximum relative error of 2.6 kcal/mol when applied with Becke3LYP-optimized geometries. The Becke3LYP DFT functional underestimates the interaction energies by few kcal/mol with relative error of 2.2 kcal/mol. Very good performance of nonpolarizable Cornell et al. force field is confirmed and this indirectly supports the view that H-bonded base pairs are primarily stabilized by electrostatic interactions.

DNA Nucleosides and Their Radical Anions: Molecular Structures and Electron Affinities

The deoxyribonucleosides have been studied to determine the properties of combinations of 2-deoxyribose with each of the isolated DNA bases for both neutral and anionic species. We have used a carefully calibrated theoretical method [Chem. Rev. 2002, 102, 231], employing the B3LYP hybrid Hartree-Fock/DFT functional with the DZP++ basis set. Predictions are made of the geometric parameters, adiabatic electron affinities, charge distributions based on natural population analysis, and decomposition enthalpy for the neutral and anionic forms of the four 2'-deoxyribonucleosides in DNA: 2'-deoxyriboadenosine (dA), 2'-deoxyribocytidine (dC), 2'-deoxyriboguanosine (dG), and 2'-deoxyribothymidine (dT). Geometric changes in the anions show that the glycosidic bond exhibits little change with excess charge for the guanosine but significant shortening for the adenosine and for the pyrimidines. The zero-point corrected adiabatic electron affinities in eV for each of the 2'-deoxyribonucleosides are as follows: 0.06, dA; 0.09, dG; 0.33, dC; and 0.44, dT. These values are uniformly greater than those of the corresponding isolated bases (-0.28, A; -0.07, G; 0.03, C; 0.20, T) and the isolated 2-deoxyribose (-0.38) at the same level of theory. The vertical detachment energies of dT and dC are substantial, 0.72 and 0.94 eV, and these anions should be observable. A high VDE, 0.91 eV, is also found for dA but its anion is unlikely to be stable due to the small AEA of 0.06 eV. The high VDE reflects the fact that the molecular structures of the anions and the corresponding neutral species are quite different. Valence character is displayed for the SOMOs of dA, dC, and dT, while some component of diffuse character is visible in the SOMO of dG. Further analysis of electronic changes upon electron attachment include an examination of the NPA charges, which show that in the neutral 2'-deoxyribonucleosides the sum of NPA charges for every base is the same, -0.28 with the sum of 2-deoxyribose charges being positive, +0.28. In the anions, the trend in charge division varies based on the nature of the excess electron in the anions. Thermodynamically, the overall enthalpy change for the reaction of water with the neutral nucleosides to give bases and ribose is approximately zero. The analogous decomposition is exothermic by 8 to 11 kcal mol-1 for the anions, indicating possible challenges for anionic gas-phase nucleoside exploration in the presence of water.

Electron Binding to Nucleic Acid Bases. Experimental and Theoretical Studies. A Review

An in-depth knowledge of an excess electron binding mechanism to DNA and RNA nucleobases is important for our understanding of radiation damage influence on the biological functions of nucleic acids, as well as for the possible use of DNA molecules as wires in molecular electronic circuits. The of anions created by electron attachment to individual nucleic acid bases is discussed in detail. The principles of the experimental and theoretical approaches to the description of these anions are outlined, and the available results concerning valence- and dipole-bound anions of nucleic acid bases are reviewed. A review with 167 references.

A theoretical study of structures and electron affinities of radical anions of guanine-cytosine, adenine-thymine, and hypoxanthine-cytosine base pairs

Adiabatic electron affinities (AEA) and structural perturbations due to addition of an excess electron to each of the neutral guanine-cytosine (G-C), adenine-thymine (A-T), and hypoxanthine-cytosine (HX-C) base pairs were studied using the self-consistent charge, density functional tight-binding (SCC-DFTB-D) method, augmented by the empirical London dispersion energy term. Performance of the SCC-DFTB-D method was examined by comparing the calculated results using it with those obtained from experiment as well as ab initio and other different density functional theoretical studies. An excellent agreement between the SCC-DFTB-D results and those obtained by the other calculations regarding the structural modifications, hydrogen bonding, and dissociation energies of the neutral and radical anion base pairs was found. It is shown that adiabatic electron affinity can be better predicted by considering reaction enthalpies of formation of the respective neutral and anionic base pairs from their respective molecular components instead of taking the difference between their total energies. The calculated AEAs of the base pairs were compared with those obtained by the bracketing method from Schaefer and coworkers, where a satisfactory agreement was found. It shows applicability of the SCC-DFTB-D method to study charged DNA models at a highly economical computational cost.

DFT investigation of dehalogenation of adenine-halouracil base pairs upon low-energy electron attachment

The energetics of the dehalogenation of adenine-halouracil base pairs (A5XU), upon attachment of low-energy electrons, was investigated by use of density functional theory. These results are compared to those of single halouracils reported previously [J. Phys. Chem. A 2002, 106, 11248-11253]. Using the B3LYP functionals it was found that the gas phase adiabatic electron affinities (EA) of halogenated base pairs (A5BrU 0.59, A5ClU 0.56, A5FU 0.47 eV) are higher than that of AU (0.32 eV) and are slightly higher or comparable to the other DNA abundant base pair, guanine-cytosine (0.49 eV). Base pairing with adenine slightly decreases the EA of the halouracils, in contrast to the substantial increase in EA on base pairing of natural bases; as a result, the probability of electron capture by halouracils when in double-stranded DNA is suggested to be substantially reduced relative to that in single-stranded DNA. Even though the activation barriers for dehalogenation are small for both BrU-A and ClU-A, only the former has negative values of both DeltaH (-0.95 kcal/mol) and DeltaG (-1.52), while the latter has negative DeltaG (-0.28) but positive DeltaH (1.27). Infinite separations into halogen anions plus the remaining A-U-5-yl neutral radical are energetically unfavorable owing to sizable halide ion, radical interactions as reported earlier for non base paired halouracils. It is found that base pairing does not change the reactive nature of the uracil-5-yl radical. The results suggest that the radiosensitization properties of halouracils should be less effective in double-stranded DNA than in single-stranded DNA.