The Charge-Transfer Band of an Oxidized Watson-Crick Guanosine-Cytidine Complex

Hole Transfer and the Resulting DNA Damage

In this review, we focus on the one-electron oxidation of DNA, which is a multipart event controlled by several competing factors. We will discuss the oxidation free energies of the four nucleobases and the electron detachment from DNA, influenced by specific interactions like hydrogen bonding and stacking interactions with neighboring sites in the double strand. The formation of a radical cation (hole) which can migrate through DNA (hole transport), depending on the sequence-specific effects and the allocation of the final oxidative damage, is also addressed. Particular attention is given to the one-electron oxidation of ds-ODN containing G:C pairs, including the complex mechanism of the deprotonation vs. hydration steps of a G:C•+ pair, as well as to the modes of formation of the two guanyl radical tautomers after deprotonation. Among the reactive oxygen species (ROS) generated in aerobic organisms by cellular metabolisms, several oxidants react with DNA. The mechanism of stable product formation and their use as biomarkers of guanine oxidation in DNA damage are also addressed.

Quantum Chemical Insights into DNA Nucleobase Oxidation: Bridging Theory and Experiment

The oxidation free energies of DNA nucleobases in aqueous solution are still matter of extensive discussion because of the contrasting results reported so far. With the aim of settling a longstanding debate about the oxidation potentials of DNA constituents, herein we report the results of state-of-the-art DFT-based molecular dynamics simulations, in which the whole solvent environment is modeled at the atomistic level, by using DFT supercell calculations, with periodic boundary conditions. Calculated vertical ionization energies are very close to those observed by photoelectron spectroscopy both in the gas phase and in solution. One-electron oxidation free energies in aqueous solution agree well with the results of differential pulse voltammetry measurements and with those inferred by photoelectron spectroscopy with the aid of theoretical computations to estimate vibrational relaxation.

Hydroxyl Radical vs. One-Electron Oxidation Reactivities in an Alternating GC Double-Stranded Oligonucleotide: A New Type Electron Hole Stabilization

We examined the reaction of hydroxyl radicals (HO•) and sulfate radical anions (SO4•-), which is generated by ionizing radiation in aqueous solutions under anoxic conditions, with an alternating GC doubled-stranded oligodeoxynucleotide (ds-ODN), i.e., the palindromic 5'-d(GCGCGC)-3'. In particular, the optical spectra of the intermediate species and associated kinetic data in the range of ns to ms were obtained via pulse radiolysis. Computational studies by means of density functional theory (DFT) for structural and time-dependent DFT for spectroscopic features were performed on 5'-d(GCGC)-3'. Comprehensively, our results suggest the addition of HO• to the G:C pair moiety, affording the [8-HO-G:C]• detectable adduct. The previous reported spectra of one-electron oxidation of a variety of ds-ODN were assigned to [G(-H+):C]• after deprotonation. Regarding 5'-d(GCGCGC)-3' ds-ODN, the spectrum at 800 ns has a completely different spectral shape and kinetic behavior. By means of calculations, we assigned the species to [G:C/C:G]•+, in which the electron hole is predicted to be delocalized on the two stacked base pairs. This transient species was further hydrated to afford the [8-HO-G:C]• detectable adduct. These remarkable findings suggest that the double-stranded alternating GC sequences allow for a new type of electron hole stabilization via delocalization over the whole sequence or part of it.

Electron Transfer Rates in Solution: Toward a Predictive First Principle Approach

Using a very recently proposed theoretical model, electron transfer rates in solution are calculated from first principles for different donor-acceptor pairs in tetrahydrofuran. We show that this approach, which integrates tunneling effects into a classical treatment of solvent motion, is able to provide reliable rate constants and their temperature dependence, even in the case of highly exergonic reactions, where Marcus’ theory usually fails.

Duplex DNA Retains the Conformational Features of Single Strands: Perspectives from MD Simulations and Quantum Chemical Computations

Molecular dynamics simulations and geometry optimizations carried out at the quantum level as well as by quantum mechanical/molecular mechanics methods predict that short, single-stranded DNA oligonucleotides adopt conformations very similar to those observed in crystallographic double-stranded B-DNA, with rise coordinates close to ≈3.3 Å. In agreement with the experimental evidence, the computational results show that DNA single strands rich in adjacent purine nucleobases assume more regular arrangements than poly-thymine. The preliminary results suggest that single-stranded poly-cytosine DNA should also retain a substantial helical order in solution. A comparison of the structures of single and double helices confirms that the B-DNA motif is a favorable arrangement also for single strands. Indeed, the optimal geometry of the complementary single helices is changed to a very small extent in the formation of the duplex.

Electron Transfer Rates in Polar and Non-Polar Environments: a Generalization of Marcus' Theory to Include an Effective Treatment of Tunneling Effects

A multistep kinetic model in which solvent motion is treated in the framework of Marcus theory and the rates of the elementary electron transfer step are evaluated at full quantum mechanical level is proposed and applied to the calculation of the rates of intramolecular electron transfer reactions in rigidly spaced D-Br-A (D = 1,1'-biphenyl radical anion, Br = androstane) compounds, for five acceptors (A) in three organic solvents with different polarity. The calculated rates agree well with experimental ones, and their temperature dependence is almost quantitatively reproduced.

The Time Scale of Electronic Resonance in Oxidized DNA as Modulated by Solvent Response: An MD/QM-MM Study

The time needed to establish electronic resonant conditions for charge transfer in oxidized DNA has been evaluated by molecular dynamics simulations followed by QM/MM computations which include counterions and a realistic solvation shell. The solvent response is predicted to take ca. 800-1000 ps to bring two guanine sites into resonance, a range of values in reasonable agreement with the estimate previously obtained by a kinetic model able to correctly reproduce the observed yield ratios of oxidative damage for several sequences of oxidized DNA.

The Dynamics of Hole Transfer in DNA

High-energy radiation and oxidizing agents can ionize DNA. One electron oxidation gives rise to a radical cation whose charge (hole) can migrate through DNA covering several hundreds of Å, eventually leading to irreversible oxidative damage and consequent disease. Understanding the thermodynamic, kinetic and chemical aspects of the hole transport in DNA is important not only for its biological consequences, but also for assessing the properties of DNA in redox sensing or labeling. Furthermore, due to hole migration, DNA could potentially play an important role in nanoelectronics, by acting as both a template and active component. Herein, we review our work on the dynamics of hole transfer in DNA carried out in the last decade. After retrieving the thermodynamic parameters needed to address the dynamics of hole transfer by voltammetric and spectroscopic experiments and quantum chemical computations, we develop a theoretical methodology which allows for a faithful interpretation of the kinetics of the hole transport in DNA and is also capable of taking into account sequence-specific effects.

Transient and Enduring Electronic Resonances Drive Coherent Long Distance Charge Transport in Molecular Wires

It is shown that the yields of oxidative damage observed in double-stranded DNA oligomers consisting of two guanines separated by adenine-thymine (A:T) _n bridges of various lengths are reliably accounted for by a multistep mechanism, in which transient and nontransient electronic resonances induce charge transport and solvent relaxation stabilizes the hole transfer products. The proposed multistep mechanism leads to results in excellent agreement with the observed yield ratios for both the short and the long distance regime; the almost distance independence of yield ratios for longer bridges ( n ≥ 3) is the consequence of the significant energy decrease of the electronic levels of the bridge, which, as the bridge length increases, become quasi-degenerate with those of the acceptor and donor groups (enduring resonance). These results provide significant guidelines for the design of novel DNA sequences to be employed in organic electronics.

Single-Stranded DNA Oligonucleotides Retain Rise Coordinates Characteristic of Double Helices

The structures of single-stranded DNA oligonucleotides from dimeric to hexameric sequences have been thoroughly investigated. Computations performed at the density functional level of theory including dispersion forces and solvation show that single-stranded helices adopt conformations very close to crystallographic B-DNA, with rise coordinates amounting up to 3.3 Å. Previous results, suggesting that single strands should be shorter than double helices, largely originated from the incompleteness of the adopted basis set. Although sensible deviations with respect to standard B-DNA are predicted, computations indicate that sequences rich in stacked adenines are the most ordered ones, favoring the B-DNA pattern and inducing regular arrangements also on flanking nucleobases. Several structural properties of double helices rich in adenine are indeed already reflected by the corresponding single strands.

Modeling DNA oxidation in water

A novel set of hole-site energies and electronic coupling parameters to be used, in the framework of the simplest tight-binding approximation, for predicting DNA hole trapping efficiencies and rates of hole transport in oxidized DNA is proposed. The novel parameters, significantly different from those previously reported in the literature, have been inferred from reliable density functional calculations, including both the sugar-phosphate ionic backbone and the effects of the aqueous environment. It is shown that most of the experimental oxidation free energies of DNA tracts and of oligonucleotides available from photoelectron spectroscopy and voltammetric measurements are reproduced with great accuracy, without the need for introducing sequence dependent parameters.

Hole delocalization over adenine tracts in single stranded DNA oligonucleotides

Adiabatic ionization energies of single stranded DNA oligonucleotides containing adenine tracts of different sizes have been computed at the DFT level and compared with the oxidation potentials determined by differential pulse voltammetry. Geometry optimizations have been performed at the full quantum mechanical level, including the sugar phosphate backbone and solvent effects. The observed progressive lowering of the ionization energy upon increasing the number of consecutive adenines is well predicted, the computed ionization potential shifts being in very good agreement with the experimental outcomes, both by using pure and hybrid functionals. The spin density of the oligonucleotide radical cations is distributed almost over the whole adenine tract, forming delocalized polarons.

Delocalized hole domains in Guanine-rich DNA oligonucleotides

Differential pulse voltammetries of guanine-rich single- and double-stranded oligonucleotides containing up to six consecutive guanines are reported. The observed progressive lowering of the first voltammetric peak potential as the number of adjacent guanines increases unambiguously points toward the establishment of delocalized hole domains; the hole stabilization energy is ca. 0.1 eV per GG step, significantly lower than that observed for AA steps.

The association constant of 5',8-cyclo-2'-deoxyguanosine with cytidine

The association of 5′,8-cyclo-2′-deoxyguanosine (cdG), a DNA tandem lesion, with its complementary base cytosine has been studied by voltammetry and NMR in chloroform, using properly silylated derivatives of the two nucleobases for increasing their solubilities. Both voltammetric data and NMR titrations indicated that the Watson-Crick complex of cytidine with cdG is weaker than that with guanosine, the difference being approximately of one order of magnitude between the two association constants.

Vibronic couplings and coherent electron transfer in bridged systems

A computational strategy to analyze the dynamics of coherent electron transfer processes in bridged systems, involving three or more electronic states, is presented.

Stacking interactions between adenines in oxidized oligonucleotides

The effects of stacking interactions on the oxidation potentials of single strand oligonucleotides containing up to four consecutive adenines, alternated with thymines and cytosines in different sequences and ratios, have been determined by means of differential pulse voltammetry. Voltammetric measurements point toward the establishment in solution of structured oligonucleotide conformations, in which the nucleobases are well stacked altogether. Molecular dynamics simulations confirm that finding, indicating that single strands assume geometrical parameters characteristic of the B-DNA form. The analysis of the voltammetric signals in terms of a simple effective tight binding quantum model leads one to infer a robust set of parameters for treating hole transfer in one-electron-oxidized DNA containing adenines and thymines.

Franck–Condon factors—Computational approaches and recent developments

Algorithms and methodologies for the calculation of Franck–Condon integrals are reviewed. Starting from the standard approach based on the Cartesian representation of the normal modes and the use of Duschinsky's transformation and recursion formulas, methods for treating large displacements of the equilibrium positions and anharmonic and non-Condon effects are considered. Finally, focusing attention on problems arising from the use of recurrence relations, some of the proposed solutions are critically reviewed along with new alternative approaches.

Reactivity of nucleosides with a hydroxyl radical in non-aqueous medium

DNA damage: The reactivity of HO(.) with silylated 2'-deoxyribonucleosides was investigated in acetonitrile by means of a time-resolved technique. The obtained rate constants were in general slightly lower than those reported for the natural nucleosides in water. Analysis of the reaction mixture by UPLC-MS revealed that HO(.) attack occurred at the nucleobase (see scheme).

Excited state dependent electron transfer of a rhenium-dipyridophenazine complex intercalated between the base pairs of DNA: a time-resolved UV-visible and IR absorption investigation into the photophysics of fac-[Re(CO)3(F2dppz)(py)]+ bound to either [poly(dA-dT)]2 or [poly(dG-dC)]2

The transient species formed following excitation of fac-[Re(CO)(3)(F(2)dppz)(py)](+) (F(2)dppz = 11,12-difluorodipyrido[3,2-a:2',3'-c]phenazine) bound to double-stranded polynucleotides [poly(dA-dT)](2) or [poly(dG-dC)](2) have been studied by transient visible and infra-red spectroscopy in both the picosecond and nanosecond time domains. The latter technique has been used to monitor both the metal complex and the DNA by monitoring the regions 1900-2100 and 1500-1750 cm(-1) respectively. These data provide direct evidence for electron transfer from guanine to the excited state of the metal complex, which proceeds both on a sub-picosecond time scale and with a lifetime of 35 ps, possibly due to the involvement of two excited states. No electron transfer is found for the [poly(dA-dT)](2) complex, although characteristic changes are seen in the DNA-region TRIR consistent with changes in the binding of the bases in the intercalation site upon excitation of the dppz-complex.