Aqueous Ionization and Electron-Donating Properties of Dinucleotides: Sequence-Specific Electronic Effects on DNA Alkylation

Ionization Patterns and Chemical Reactivity of Cytosine-Guanine Watson-Crick Pairs

Gas-phase and aqueous vertical ionization potentials, vIPgas and vIPaq respectively and measurements of the molecular electrostatic and local ionization maps calculated at the DFT/B3LYP-D3/ 6-311+G** level of theory and the C-PCM reaction field model for single- and double-stranded CpG and 5MeCpG pairs show that the vIPaq for single- and double-stranded pairs of C-G and 5MeC-G are practically the same, in the range of 5.79 to 5.81 eV. The aqueous adiabatic ionization potentials for single-stranded CpG and 5MeCpG are 5.52 eV and 5.51 eV respectively and they reflect the nuclear reorganization that takes place after the abstraction of the electron. The aqueous adiabatic ionization energy values that correspond to the CpG+. radical cation and the hydrated electron, e-,, being at infinite distance, adIPaq+Vo, are 3.92 eV and 3.91 eV respectively with (Vo=-1.6 eV) Analysis of data suggest that the HOMO-LUMO energy gap in the hard/soft-acid/base (HSAB) concept cannot be used a priori to determine the effect of cytosine methylation on the guanine enhanced oxidative damage in DNA.

Using DNA origami nanostructures to determine absolute cross sections for UV photon-induced DNA strand breakage

We have characterized ultraviolet (UV) photon-induced DNA strand break processes by determination of absolute cross sections for photoabsorption and for sequence-specific DNA single strand breakage induced by photons in an energy range from 6.50 to 8.94 eV. These represent the lowest-energy photons able to induce DNA strand breaks. Oligonucleotide targets are immobilized on a UV transparent substrate in controlled quantities through attachment to DNA origami templates. Photon-induced dissociation of single DNA strands is visualized and quantified using atomic force microscopy. The obtained quantum yields for strand breakage vary between 0.06 and 0.5, indicating highly efficient DNA strand breakage by UV photons, which is clearly dependent on the photon energy. Above the ionization threshold strand breakage becomes clearly the dominant form of DNA radiation damage, which is then also dependent on the nucleotide sequence.

Methylation of zebularine investigated using density functional theory calculations

Deoxyribonucleic acid (DNA) methylation is an epigenetic phenomenon, which adds methyl groups into DNA. This study reveals methylation of a nucleoside antibiotic drug 1-(β-D-ribofuranosyl)-2-pyrimidinone (zebularine or zeb) with respect to its methylated analog, 1-(β-D-ribofuranosyl)-5-methyl-2-pyrimidinone (d5) using density functional theory calculations in valence electronic space. Very similar infrared spectra suggest that zeb and d5 do not differ by types of the chemical bonds, but distinctly different Raman spectra of the nucleoside pair reveal that the impact caused by methylation of zeb can be significant. Further valence orbital-based information details on valence electronic structural changes caused by methylation of zebularine. Frontier orbitals in momentum space and position space of the molecules respond differently to methylation. Based on the additional methyl electron density concentration in d5, orbitals affected by the methyl moiety are classified into primary and secondary contributors. Primary methyl contributions include MO8 (57a), MO18 (47a), and MO37 (28a) of d5, which concentrates on methyl and the base moieties, suggest certain connection to their Frontier orbitals. The primary and secondary methyl affected orbitals provide useful information on chemical bonding mechanism of the methylation in zebularine.

Solvent exposure associated with single abasic sites alters the base sequence dependence of oxidation of guanine in DNA in GG sequence contexts

The effect of exposure of guanine in double-stranded oligonucleotides to aqueous solvent during oxidation by one-electron oxidants was investigated by introducing single synthetic tetrahydrofuran-type abasic sites (Ab) either adjacent to or opposite tandem GG sequences. The selective oxidation of guanine was initiated by photoexcitation of the aromatic sensitizers riboflavin and a pyrene derivative, and by the relatively small negatively charged carbonate radical anion. The relative rates of oxidation of the 5'- and 3' side G in runs of 5'⋅⋅⋅GG⋅⋅⋅ (evaluated by standard hot alkali treatment of the damaged DNA strand followed by high resolution gel electrophoresis of the cleavage fragments) are markedly affected by adjacent abasic sites either on the same or opposite strand. For example, in fully double-stranded DNA or one with an Ab adjacent to the 5'-G, the 5'-G/3'-G damage ratio is ≥4, but is inverted (<1.0) with the Ab adjacent to the 3'-G. These striking effects of Ab are attributed to the preferential localization of the "hole" on the most solvent-exposed guanine regardless of the size, charge, or reduction potential of the oxidizing species.

Vertical ionization potentials of nucleobases in a fully solvated DNA environment

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

Solvent effects on ionization potentials of guanine runs and chemically modified guanine in duplex DNA: effect of electrostatic interaction and its reduction due to solvent

We examined the ionization potential (IP) corresponding to the free energy of a hole on duplex DNA by semiempirical molecular orbital theory with a continuum solvent model. As for the contiguous guanines (a guanine run), we found that the IP in the gas phase significantly decreases with the increasing number of nucleotide pairs of the guanine run, whereas the IP in water (OP, oxidation potential) only slightly does. The latter result is consistent with the experimental result for DNA oligomers in water. This decrease in the IP is mainly due to the attractive electrostatic interaction between the hole and a nucleotide pair in the duplex DNA. This interaction is reduced in water, which results in the small decrease in the IP in water. This mechanism explains the discrepancy between the experimental result and the previous computational results obtained by neglecting the solvent. As for the chemically modified guanine, the previous work showed that the removal of some solvent (water) molecules due to the attachment of a neutral functional group to a guanine in a duplex DNA stabilizes the hole on the guanine. One might naively have expected the opposite case, since a polar solvent usually stabilizes ions. This mechanism also explains this unexpected stabilization of a hole as follows. When some water molecules are removed, the attractive electrostatic interaction stabilizing the hole increases, and thus, the hole is stabilized. In order to design the hole energetics by a chemical modification of DNA, this mechanism has to be taken into account and can be used.

Characterizing radiation-induced oxidation of DNA by way of the monohydrated guanine-cytosine radical cation

The interaction of one water molecule with the guanine-cytosine radical cation has been studied with ab initio and density functional methods in order to help elucidate the nature of oxidized aqueous DNA. The theoretical spin density of [GC]*(+) reveals that the radical center is localized on guanine. The adiabatic ionization potential lowers from 7.63 to 6.71 eV in concurrence with the formation of the Watson-Crick base pair and hydration by one water molecule. A natural bond orbital analysis of partial charges shows that approximately 80% of the positive charge persists on guanine upon hydration and formation of the Watson-Crick base pair with cytosine. Hydration energies were computed with second-order Z-averaged perturbation theory (ZAPT2) using the aug-cc-pVDZ basis set at 11 stationary points on the B3LYP/DZP++ potential energy surface. The hydration energy at the global minimum is 14.2 kcal mol(-1). The lowest energy structures correspond to hydration near the glycosidic bond sites. Structural changes in the Watson-Crick base pair are predominantly seen for monohydration in the groove regions of double-helix DNA.

Base and phosphate electron detachment energies of deoxyribonucleotide anions

Photoelectron spectra of deoxyribonucleotide anions are interpreted with ab initio, electron propagator calculations. Ground-state structures display hydrogen bonds which are not present in less stable minima that resemble Watson-Crick fragment geometries. For the adenosine and thymidine anions, there are two vertical electron detachment energies (VEDEs) within 0.1 eV of each other that correspond to phosphate- and base-centered Dyson orbitals (DOs). The first VEDE of the cytidine anion belongs to a phosphate-centered DO. The anomalously low VEDE of the guanosine anion is assigned to a base-centered, pi DO. Higher VEDEs of all four anions also are assigned.

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.

Determination of Redox Potentials for the Watson−Crick Base Pairs, DNA Nucleosides, and Relevant Nucleoside Analogues

Redox potentials for the DNA nucleobases and nucleosides, various relevant nucleoside analogues, Watson-Crick base pairs, and seven organic dyes are presented based on DFT/B3LYP/6-31++G(d,p) and B3YLP/6-311+G(2df,p)//B3LYP/6-31+G* levels of calculations. The values are determined from an experimentally calibrated set of equations that correlate the vertical ionization (electron affinity) energy of 20 organic molecules with their experimental reversible oxidation (reduction) potential. Our results are in good agreement with those estimated experimentally for the DNA nucleosides in acetonitrile solutions (Seidel et al. J. Phys. Chem. 1996, 100, 5541). We have found that nucleosides with anti conformation exhibit lower oxidation potentials than the corresponding syn conformers. The lowering in the oxidation potential is due to the formation of an intramolecular hydrogen bonding interaction between the 5'-OH group of the sugar and the N3 of the purine bases or C2=O of the pyrimidine bases in the syn conformation. Pairing of adenine or guanine with its complementary pyrimidine base decreases its oxidation potential by 0.15 or 0.28 V, respectively. The calculated energy difference between the oxidation potential for the G.C base pair and that of the guanine base is in good agreement with the experimental value estimated recently (0.34 V: Caruso, T.; et al. J. Am. Chem. Soc. 2005, 127, 15040). The complete and consistent set of reversible redox values determined in this work for the DNA constituents is expected to be of considerable value to those studying charge and electronic energy transfer in DNA.

One- and Two-Photon Ionization of DNA Single and Double Helices Studied by Laser Flash Photolysis at 266 nm

The ionization of the DNA single and double helices (dA)20, (dT)20, (dAdT)10(dAdT)10 and (dA)20(dT)20, induced by nanosecond pulses at 266 nm, is studied by time-resolved absorption spectroscopy. The variation of the hydrated electron concentration with the absorbed laser intensity shows that, in addition to two-photon ionization, one-photon ionization takes place for (dAdT)10(dAdT)10, (dA)20(dT)20 and (dA)20 but not for (dT)20. The spectra of all adenine-containing oligomers at the microsecond time-scale correspond to the adenine deprotonated radical formed in concentrations comparable to that of the hydrated electron. The quantum yield for one-photon ionization of the oligomers (ca. 10(-3)) is higher by at least 1 order of magnitude than that of dAMP, showing clearly that organization of the bases in single and double helices leads to an important lowering of the ionization potential. The propensity of (dAdT)10(dAdT)10, containing alternating adenine-thymine sequences, to undergo one-photon ionization is lower than that of (dA)20(dT)20 and (dA)20, containing adenine runs. Pairing of the (dA)20 with the complementary strand leads to a decrease of quantum yield for one photon ionization by about a factor of 2.

Activation barriers for DNA alkylation by carcinogenic methane diazonium ions

Methylation reactions of the DNA bases with the methane diazonium ion, which is the reactive intermediate formed from several carcinogenic methylating agents, were examined. The SN2 transition states of the methylation reactions at N7, N3, and O6 of guanine; N7, N3, and N1 of adenine; N3 and O2 of cytosine; and O2 and O4 of thymine were calculated using the B3LYP density functional method. Solvation effects were examined using the conductor-like polarizable continuum method and the combined discrete/SCRF method. The transition states for reactions at guanine N3, adenine N7, and adenine N1 are influenced by steric interactions between the methane diazonium ion and exocyclic amino groups. Both in the gas phase and in aqueous solution, the methylation reactions at N atoms have transition states that are looser, and generally occur earlier along the reaction pathways than reactions at O atoms. The forming bonds in the transition states in water are 0.03 to 0.13 A shorter than those observed in the gas phase, and the activation energies are 13 to 35 kcal/mol higher. The combined discrete/SCRF solvation energy calculations using base-water complexes with three water molecules yield base solvation energies that are larger than those obtained from the CPCM continuum method, especially for cytosine. Reactivities calculated using barriers obtained with the discrete/SCRF method are consistent with the experimentally observed high reactivity at N7 of guanine.

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

At the very early time of irradiation, ballistic secondary electrons are produced as the most abundant of the radiolytic species directly within DNA or its environment. Here, we demonstrate the propensity of such low-energy (<3 eV) electrons to damage DNA bases via an effective loss of hydrogen located at the specific nitrogen positions. Since this site is directly implicated in the bonding of nucleobases within DNA and since dehydrogenation of the nucleic acid bases has been observed to be the predominant dissociative channel, the present findings foreshadow significant implications for the initial molecular processes leading to genotoxicity in living cells following unwanted or intended exposure to ionizing radiation (e.g., sunbathing, air travel, radiotherapy, etc.).

Theoretical ab initio study of the effects of methylation on structure and stability of G:C Watson-Crick base pair

Methylation of DNA occurs most readily at N(3), N(7), and O(6) of purine bases and N(3) and O(2) of pyrimidines. Methylated bases are continuously formed through endogenous and exogenous mechanisms. The results of a theoretical ab initio study on the methylation of G:C base pair components are reported. The geometries of the local minima were optimized without symmetry restrictions by the gradient procedure at DFT level of theory and were verified by energy second derivative calculations. The standard 6-31G(d) basis set was used. The single-point calculations have been performed at the MP2/6-31G(d,p), MP2/6-31++G(d,p), and MP2/6-311++G(2d,2p) levels of theory. The geometrical parameters, relative stability and counterpoise corrected interaction energies are reported. Also, using a variation-perturbation energy decomposition scheme we have found the vital contributions to the total interaction energy.

Biomolecular electronic devices based on self-organized deoxyguanosine nanocrystals

We report on a new class of hybrid electronic devices based on a DNA nucleoside (deoxyguanosine lipophilic derivative) whose assembled polymeric ribbons interconnect a submicron metallic gate. The device exhibits large conductivity at room temperature, rectifying behavior and strong current-voltage hysteresis. The transport mechanism through the molecules is investigated by comparing films with different self-assembling morphology. We found that the main transport mechanism is connected to pi-pi interactions between guanosine molecules and to the formation of a strong dipole along ribbons, consistently with the results of our first-principles calculations.