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 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.  more » « less
Award ID(s):
1920110
NSF-PAR ID:
10286737
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
Volume:
22
Issue:
9
ISSN:
1463-9076
Page Range / eLocation ID:
5078 to 5089
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    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 macromolecule 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. 
    more » « less
  2. Abstract

    We report relatively persistent, open‐shell thiophene‐based double helices, radical cations 1•+‐TMS12and 2•+‐TMS8. Closed‐shell neutral double helices, 1‐TMS12and 2‐TMS8, have nearly identical first oxidation potentials,E+/0 ≈ +1.33 V, corresponding to reversible oxidation to their radical cations. The radical cations are generated, using tungsten hexachloride in dichloromethane (DCM) as an oxidant,E+/0 ≈ +1.56 V. EPR spectra consist of a relatively sharp singlet peak with an unusually lowg‐value of 2.001–2.002, thus suggesting exclusive delocalization of spin density over π‐conjugated system consisting of carbon atoms only. DFT computations confirm these findings, as only negligible fraction of spin density is found on sulfur and silicon atoms and the spin density is delocalized over a single tetrathiophene moiety. For radical cation, 1•+‐TMS12, energy level of the singly occupied molecular orbital (SOMO) lies below the four highest occupied molecular orbitals (HOMOs), thus indicating the SOMO–HOMO inversion (SHI) and therefore, violating the Aufbau principle. 1•+‐TMS12has a half‐life of the order of only 5 min at room temperature. EPR peak intensity of 2•+‐TMS8, which does not show SHI, is practically unchanged over at least 2 h.

     
    more » « less
  3. Abstract

    The replacement of a CH group of benzene by a triel (Tr) atom places a positive region of electrostatic potential near the Tr atom in the plane of the aromatic ring. This σ‐hole can interact with an X lone pair of XCCH (X=F, Cl, Br, and I) to form a triel bond (TrB). The interaction energy between C5H5Tr and FCCH lies in the range between 2.2 and 4.4 kcal/mol, in the order Tr=B+cation above the ring pulls density toward itself and thus magnifies the Tr σ‐hole. The TrB to the XCCH nucleophile is thereby magnified as is the strength of the TrB. This positive cooperativity is particularly large for Tr=B.

     
    more » « less
  4. 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 is 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. 
    more » « less
  5. 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 splittings 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. 
    more » « less