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

Photo-induced radical transformations of tosyl cyanide

Tosyl cyanide is a commonly used reagent for cyanation and sulfonylation in organic synthesis and pharmaceutical chemistry. The photocatalytic transformations of tosyl cyanide are generally conducted under mild conditions. This minireview summarizes the recent progress of radical-involved transformations of tosyl cyanide via photo-induced cyanation or sulfonylcyanation.

Theoretical Modeling of Redox Potentials of Biomolecules

The estimation of the redox potentials of biologically relevant systems by means of theoretical-computational approaches still represents a challenge. In fact, the size of these systems typically does not allow a full quantum-mechanical treatment needed to describe electron loss/gain in such a complex environment, where the redox process takes place. Therefore, a number of different theoretical strategies have been developed so far to make the calculation of the redox free energy feasible with current computational resources. In this review, we provide a survey of such theoretical-computational approaches used in this context, highlighting their physical principles and discussing their advantages and limitations. Several examples of these approaches applied to the estimation of the redox potentials of both proteins and nucleic acids are described and critically discussed. Finally, general considerations on the most promising strategies are reported.

Direct formation of the C5'-radical in the sugar-phosphate backbone of DNA by high-energy radiation

Neutral sugar radicals formed in DNA sugar-phosphate backbone are well-established as precursors of biologically important damage such as DNA strand scission and cross-linking. In this work, we present electron spin resonance (ESR) evidence showing that the sugar radical at C5' (C5'(•)) is one of the most abundant (ca. 30%) sugar radicals formed by γ- and Ar ion-beam irradiated hydrated DNA samples. Taking dimethyl phosphate as a model of sugar-phosphate backbone, ESR and theoretical (DFT) studies of γ-irradiated dimethyl phosphate were carried out. CH(3)OP(O(2)(-))OCH(2)(•) is formed via deprotonation from the methyl group of directly ionized dimethyl phosphate at 77 K. The formation of CH(3)OP(O(2)(-))OCH(2)(•) is independent of dimethyl phosphate concentration (neat or in aqueous solution) or pH. ESR spectra of C5'(•) found in DNA and of CH(3)OP(O(2)(-))OCH(2)(•) do not show an observable β-phosphorus hyperfine coupling (HFC). Furthermore, C5'(•) found in DNA does not show a significant C4'-H β-proton HFC. Applying the DFT/B3LYP/6-31G(d) method, a study of conformational dependence of the phosphorus HFC in CH(3)OP(O(2)(-))OCH(2)(•) shows that in its minimum energy conformation, CH(3)OP(O(2)(-))OCH(2)(•), has a negligible β-phosphorus HFC. On the basis of these results, the formation of radiation-induced C5'(•) is proposed to occur via a very rapid deprotonation from the directly ionized sugar-phosphate backbone, and the rate of this deprotonation must be faster than that of energetically downhill transfer of the unpaired spin (hole) from ionized sugar-phosphate backbone to the DNA bases. Moreover, C5'(•) in irradiated DNA is found to be in a conformation that does not exhibit β-proton or β-phosphorus HFCs.

UV photodissociation dynamics of deprotonated 2'-deoxyadenosine 5'-monophosphate [5'-dAMP-H]-

The UV photodissociation dynamics of deprotonated 2'-deoxyadenosine 5'-monophosphate ([5'-dAMP-H](-)) has been studied using a unique technique based on the coincident detection of the ion and the neutral fragments. The observed fragment ions are m/z 79 (PO(3)(-)), 97 (H(2)PO(4)(-)), 134 ([A-H](-)), 177 ([dAMP-H-A-H(2)O](-)), and 195 ([dAMP-H-A](-)), where "A" refers to a neutral adenine molecule. The relative abundances are comparable to that found in previous studies on [5'-dAMP-H](-) employing different excitation processes, i.e., collisions and UV photons. The fragmentation times of the major channels have been measured, and are all found to be on the microsecond time scale. The fragmentation mechanisms for all channels have been characterized using velocity correlation plots of the ion and neutral fragment(s). The findings show that none of the dissociation channels of [5'-dAMP-H](-) is UV specific and all proceed via statistical fragmentation on the ground state after internal conversion, a result similar to fragmentations induced by collisions.

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.

Determination of the electron-detachment energies of 2'-deoxyguanosine 5'-monophosphate anion: influence of the conformation

The vertical electron-detachment energies (VDEs) of the singly charged 2'-deoxyguanosine 5'-monophosphate anion (dGMP-) are determined by using the multiconfigurational second-order perturbation CASPT2 method at the MP2 ground-state equilibrium geometry of relevant conformers. The origin of the unique low-energy band in the gas phase photoelectron spectrum of dGMP-, with maximum at around 5.05 eV, is unambiguously assigned to electron detachment from the highest occupied molecular orbital of pi-character belonging to guanine fragment of a syn conformation. The presence of a short H-bond linking the 2-amino and phosphate groups, the guanine moiety acting as proton donor, is precisely responsible for the pronounced decrease of the computed VDE with respect to that obtained in other conformations. As a whole, the present research supports the nucleobase as the site with the lowest ionization potential in negatively charged (deprotonated) nucleotides at the most stable conformations as well as for B-DNA-like type arrangements, in agreement with experimental evidence.

Elementary lesions in DNA subunits: electron, hydrogen atom, proton, and hydride transfers

When DNA is damaged by ionizing radiation, the genes in a cell may acquire mutations or the cell could die. The smallest known DNA-damaging unit is an electron, often low-energy secondary electrons. Additional electrons and transfers involving hydrogen atoms, protons, and hydride anions can damage DNA subunits, including individual nucleobases and nucleoside pairs. Researchers would like to better understand the molecular mechanisms involved in DNA damage from ionizing radiation. In this Account, we highlight our theoretical investigations of the molecular mechanisms of DNA damage using quantum mechanical models. Our investigations use robust theoretical methods with computations conducted in the gas phase and with solution models. We calculate adiabatic electron affinities (AEAs), which describe the energetics of electronic attachment to closed-shell DNA subunits, for the free bases, nucleosides, nucleotides, base pairs, and single and double DNA strand units. Electron affinities for free nucleobases yield the order uracil > thymine > cytosine > guanine > adenine and the same order for the DNA nucleosides, mononucleotides, and nucleoside 3',5'-diphosphates. AEA values increase steadily with the size and complexity of the system because of greater hydration, glycosylation, nucleotide formation, and base pairing. We predict and experimental results partially confirm that most of the more complex and hydrated species are observable as radical anions. Our modeling studies indicate that depyrimidination reactions of radical anion nucleosides release cytosine more often than thymine. Recent experimental results support those findings. Our theoretical studies of DNA base-pair radical anions predict increases in electron affinity accompanying H bonding and solvation. Electron addition facilitates some proton transfers in these pairs, which results in strongly perturbed pairing configurations. Of all nucleobase moieties within the more complex radical anion systems, thymine is best able to retain a negative charge. Charge and spin are well-separated in some of these systems. Radical species derived via hydrogen abstraction from DNA subunits yield large AEA values because they form closed-shell anions. Our studies predict single-strand breaks following H abstraction from nucleotides. Some H-abstraction processes in the DNA base pairs lead to severe distortions in pairing configuration based on our calculations. This body of systematic theoretical studies provides realistic descriptions of some events that lead to elementary DNA lesions, while providing rationalizations for many observed phenomena. Such approaches can lead to the design of new experiments, which would contribute to our understanding of the chemical physics of nucleic acids.

Ionization energies of the nucleotides

The vertical ionization energies of the four nucleotides have been computed. Geometries have been chosen to mimic orientations as they appear in B-DNA. The negative charge on the phosphate was neutralized by protonation, and also by the inclusion of counterions. Calculations have been performed with electron propagator methods (P3), Møller-Plesset second-order perturbation theory, and density functional theory to determine the nature of the orbitals associated with the highest lying ionization energies. Calculations at the MP2/6-311G(d,p)//P3/6-311G(d,p) level of theory yield vertical ionization energies for 5'-dTMP 9.05 eV, for 5'-dCMP 8.40 eV, for 5'-dAMP 8.16 eV and for 5'-dGMP 7.96 eV. In all cases the highest occupied molecular orbital resides on the base moieties.

One-electron oxidation of 2'-deoxyadenosine-5'-phosphate: comparisons of theoretical calculations with experimental values

A recent paper by Hou et al. (Hou, R.; Gu, J.; Xie, Y.; Yi, X.; Schaefer, H. F. J. Phys. Chem. B 2005, 109, 22053) on 2'-deoxyadenosine-5'-phosphate (5'-dAMP) reports calculations on one-electron oxidation of the 5'-dAMP anion. The paper presents a very interesting observation that, for the radical produced by electron removal, the unpaired spin density resides on both the phosphate and the adenine base moieties. There are also indications that this radical has a weakened C5'-O5' bond, and it is said that this may be the origin of a single-strand break in DNA. New calculations have been performed to show that the spin density on the phosphate is dependent on the charge on the phosphate. The use of the B3LYP method with the 6-31G(d) basis set yields results very similar to those obtained with the much larger B3LYP/DZP++ basis set in computing the structures of one electron oxidized 5'-dAMP. New calculations on the isotropic hyperfine couplings in 5'-dAMP are presented to show under just what conditions one might expect to see small amounts of unpaired spin density on the phosphates. Results show that this may occur in gas-phase studies of nucleotides but, most likely, not in DNA.

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.