skip to main content

Title: One-electron oxidation of ds(5′-GGG-3′) and ds(5′-G(8OG)G-3′) and the nature of hole distribution: a density functional theory (DFT) study
Of particular interest in radiation-induced charge transfer processes in DNA is the extent of hole localization immediately after ionization and subsequent relaxation. To address this, we considered double stranded oligomers containing guanine (G) and 8-oxoguanine (8OG), i.e. , ds(5′-GGG-3′) and ds(5′-G8OGG-3′) in B-DNA conformation. Using DFT, we calculated a variety of properties, viz. , vertical and adiabatic ionization potentials, spin density distributions in oxidized stacks, solvent and solute reorganization energies and one-electron oxidation potential ( E 0 ) in the aqueous phase. Calculations for the vertical state of the -GGG- cation radical showed that the spin was found mainly (67%) on the middle G. However, upon relaxation to the adiabatic -GGG- cation radical, the spin localized (96%) on the 5′-G, as observed in experiments. Hole localizations on the middle G and 3′-G were higher in energy by 0.5 kcal mol −1 and 0.4 kcal mol −1 , respectively, than that of 5′-G. In the -G8OGG- cation radical, the spin localized only on the 8OG in both vertical and adiabatic states. The calculated vertical ionization potentials of -GGG- and -G8OGG- stacks were found to be lower than that of the vertical ionization potential of a single G in DNA. The calculated more » E 0 values of -GGG- and -G8OGG- stacks are 1.15 and 0.90 V, respectively, which owing to stacking effects are substantially lower than the corresponding experimental E 0 values of their monomers (1.49 and 1.18 V, respectively). SOMO to HOMO level switching is observed in these oxidized stacks. Consequently, our calculations predict that local double oxidations in DNA will form triplet diradical states, which are especially significant for high LET radiations. « less
; ; ;
Award ID(s):
Publication Date:
Journal Name:
Physical Chemistry Chemical Physics
Page Range or eLocation-ID:
5078 to 5089
Sponsoring Org:
National Science Foundation
More Like this
  1. Among the radicals (hydroxyl radical (•OH), hydrogen atom (H•), and solvated electron (esol−)) that are generated via water radiolysis, •OH has been shown to be the main transient species responsible for radiation damage to DNA via the indirect effect. Reactions of these radicals with DNA-model systems (bases, nucleosides, nucleotides, polynucleotides of defined sequences, single stranded (ss) and double stranded (ds) highly polymeric DNA, nucleohistones) were extensively investigated. The timescale of the reactions of these radicals with DNA-models range from nanoseconds (ns) to microseconds (µs) at ambient temperature and are controlled by diffusion or activation. However, those studies carried out in dilute solutions that model radiation damage to DNA via indirect action do not turn out to be valid in dense biological medium, where solute and water molecules are in close contact (e.g., in cellular environment). In that case, the initial species formed from water radiolysis are two radicals that are ultrashort-lived and charged: the water cation radical (H2O•+) and prethermalized electron. These species are captured by target biomolecules (e.g., DNA, proteins, etc.) in competition with their inherent pathways of proton transfer and relaxation occurring in less than 1 picosecond. In addition, the direct-type effects of radiation, i.e., ionization of macromoleculemore »plus excitations proximate to ionizations, become important. The holes (i.e., unpaired spin or cation radical sites) created by ionization undergo fast spin transfer across DNA subunits. The exploration of the above-mentioned ultrafast processes is crucial to elucidate our understanding of the mechanisms that are involved in causing DNA damage via direct-type effects of radiation. Only recently, investigations of these ultrafast processes have been attempted by studying concentrated solutions of nucleosides/tides under ambient conditions. Recent advancements of laser-driven picosecond electron accelerators have provided an opportunity to address some long-term puzzling questions in the context of direct-type and indirect effects of DNA damage. In this review, we have presented key findings that are important to elucidate mechanisms of complex processes including excess electron-mediated bond breakage and hole transfer, occurring at the single nucleoside/tide level.« less
  2. The spectroscopic identification of Bi 4 has been very elusive. Two constitutional Bi 4 isomers of T d and C 2v symmetry are investigated and each is found to be a local energetic minimum. The optimized geometries and vibrational frequencies of these two isomers are obtained at the CCSD(T)/cc-pVQZ-PP level of theory, utilizing the Stoll, Metz, and Dolg 60-electron effective core potential. The fundamental frequencies of the T d isomer are obtained at the same level of theory. The focal point analysis method, from a maximum basis set of cc-pV5Z-PP, and proceeding to a maximum correlation method of CCSDTQ, was employed to determine the dissociation energy of Bi 4 ( T d ) into two Bi 2 and the adiabatic energy difference between the C 2v and T d isomers of Bi 4 . These quantities are predicted to be +65 kcal mol −1 and +39 kcal mol −1 , respectively. Two electron vertical excitation energies between the T d and C 2v electronic configurations are computed to be 156 kcal mol −1 for the T d isomer and 9 kcal mol −1 for the C 2v isomer. The most probable approach to laboratory spectroscopic identification of Bi 4 ismore »via an infrared spectrum. The predicted fundamentals (cm −1 ) with harmonic IR intensities in parentheses (km mol −1 ) are 94(0), 123(0.23), and 167(0) for the T d isomer. The moderate IR intensity for the only allowed fundamental may explain why Bi 4 has yet to be observed. Through natural bond orbital analysis, the C 2v isomer of Bi 4 was discovered to exhibit “long-bonding” between the furthest apart ‘wing’ atoms. This long-bonding is postulated to be facilitated by the σ-bonding orbital between the ‘spine’ atoms of the C 2v isomer.« less
  3. First, high-resolution sub-Doppler infrared spectroscopic results for cyclopentyl radical (C 5 H 9 ) are reported on the α-CH stretch fundamental with suppression of spectral congestion achieved by adiabatic cooling to T rot ≈ 19(4) K in a slit jet expansion. Surprisingly, cyclopentyl radical exhibits a rotationally assignable infrared spectrum, despite 3N − 6 = 36 vibrational modes and an upper vibrational state density (ρ ≈ 40–90 #/cm −1 ) in the critical regime (ρ ≈ 100 #/cm −1 ) necessary for onset of intramolecular vibrational relaxation (IVR) dynamics. Such high-resolution data for cyclopentyl radical permit detailed fits to a rigid-rotor asymmetric top Hamiltonian, initial structural information for ground and vibrationally excited states, and opportunities for detailed comparison with theoretical predictions. Specifically, high level ab initio calculations at the coupled-cluster singles, doubles, and perturbative triples (CCSD(T))/ANO0, 1 level are used to calculate an out-of-plane bending potential, which reveals a C 2 symmetry double minimum 1D energy surface over a C 2v transition state. The inversion barrier [V barrier ≈ 3.7(1) kcal/mol] is much larger than the effective moment of inertia for out-of-plane bending, resulting in localization of the cyclopentyl wavefunction near its C 2 symmetry equilibrium geometry and tunneling splittingsmore »for the ground state too small (<1 MHz) to be resolved under sub-Doppler slit jet conditions. The persistence of fully resolved high-resolution infrared spectroscopy for such large cyclic polyatomic radicals at high vibrational state densities suggests a “deceleration” of IVR for a cycloalkane ring topology, much as low frequency torsion/methyl rotation degrees of freedom have demonstrated a corresponding “acceleration” of IVR processes in linear hydrocarbons.« less
  4. Hydrogen-rich cation radicals (GATT + 2H) + ˙ and (AGTT + 2H) + ˙ represent oligonucleotide models of charged hydrogen atom adducts to DNA. These tetranucleotide cation radicals were generated in the gas phase by one-electron reduction of the respective (GATT + 2H) 2+ and (AGTT + 2H) 2+ dications in which the charging protons were placed on the guanine and adenine nucleobases. We used wavelength-dependent UV/Vis photodissociation in the valence-electron excitation region of 210–700 nm to produce action spectra of (GATT + 2H) + ˙ and (AGTT + 2H) + ˙ that showed radical-associated absorption bands in the near-UV (330 nm) and visible (400–440 nm) regions. Born–Oppenheimer molecular dynamics and density-functional theory calculations were used to obtain and rank by energy multiple (GATT + 2H) dication and cation-radical structures. Time-dependent density functional theory (TD-DFT) calculations of excited-state energies and electronic transitions in (GATT + 2H) + ˙ were augmented by vibronic spectra calculations at 310 K for selected low-energy cation radicals to provide a match with the action spectrum. The stable product of one-electron reduction was identified as having a 7,8-dihydroguanine cation radical moiety, formed by intramolecular hydrogen atom migration from adenine N-1–H. The hydrogen migration was calculated tomore »have a transition state with a low activation energy, E a = 96.5 kJ mol −1 , and positive activation entropy, Δ S ‡ = 75 J mol −1 K −1 . This allowed for a fast isomerization of the primary reduction products on the ion-trap time scale of 150 ms that was substantially accelerated by highly exothermic electron transfer.« less
  5. Intermolecular interactions between radicals and closed-shell molecules are ubiquitous in chemical processes, ranging from the benchtop to the atmosphere and extraterrestrial space. While energy decomposition analysis (EDA) schemes for closed-shell molecules can be generalized for studying radical–molecule interactions, they face challenges arising from the unique characteristics of the electronic structure of open-shell species. In this work, we introduce additional steps that are necessary for the proper treatment of radical–molecule interactions to our previously developed unrestricted Absolutely Localized Molecular Orbital (uALMO)-EDA based on density functional theory calculations. A “polarize-then-depolarize” (PtD) scheme is used to remove arbitrariness in the definition of the frozen wavefunction, rendering the ALMO-EDA results independent of the orientation of the unpaired electron obtained from isolated fragment calculations. The contribution of radical rehybridization to polarization energies is evaluated. It is also valuable to monitor the wavefunction stability of each intermediate state, as well as their associated spin density profiles, to ensure the EDA results correspond to a desired electronic state. These radical extensions are incorporated into the “vertical” and “adiabatic” variants of uALMO-EDA for studies of energy changes and property shifts upon complexation. The EDA is validated on two model complexes, H 2 O⋯˙F and FH⋯˙OH. It is thenmore »applied to several chemically interesting radical–molecule complexes, including the sandwiched and T-shaped benzene dimer radical cation, complexes of pyridine with benzene and naphthalene radical cations, binary and ternary complexes of the hydroxyl radical with water (˙OH(H 2 O) and ˙OH(H 2 O) 2 ), and the pre-reactive complexes and transition states in the ˙OH + HCHO and ˙OH + CH 3 CHO reactions. These examples suggest that this second generation uALMO-EDA is a useful tool for furthering one's understanding of both energetic and property changes associated with radical–molecule interactions.« less