Proton-coupled electron transfer in DNA on formation of radiation-produced ion radicals

A review on point mutations via proton transfer in DNA base pairs in the absence and presence of electric fields

This is a comprehensive review that focuses on spontaneous mutations that may occur during DNA replication, the fundamental process responsible for transferring genetic information. In 1963, Löwdin postulated that these mutations are primarily a result of proton transfer reactions within the hydrogen-bonded DNA base pairs. The single and double proton transfer reactions within the base pairs result in zwitterions and rare tautomers, respectively. For persistent mutations, these products must be generated at high rates and should be thermodynamically stable. This review covers the proton transfer reactions studied experimentally and computationally. The review also examines the influence of externally applied electric fields on the thermodynamics and kinetics of proton transfer reactions within DNA base pairs, and their biological implications.

The Excited State Dynamics of a Mutagenic Guanosine Etheno Adduct Investigated by Femtosecond Fluorescence Spectroscopy and Quantum Mechanical Calculations

Femtosecond fluorescence upconversion experiments were combined with CASPT2 and time dependent DFT calculations to characterize the excited state dynamics of the mutagenic etheno adduct 1,N2-etheno-2'-deoxyguanosine (ϵdG). This endogenously formed lesion is attracting great interest because of its ubiquity in human tissues and its highly mutagenic properties. The ϵdG fluorescence is strongly modified with respect to that of the canonical nucleoside dG, notably by an about 6-fold increase in fluorescence lifetime and quantum yield at neutral pH. In addition, femtosecond fluorescence upconversion experiments reveal the presence of two emission bands with maxima at 335 nm for the shorter-lived and 425 nm for the longer-lived. Quantum mechanical calculations rationalize these findings and provide absorption and fluorescence spectral shapes similar to the experimental ones. Two different bright minima are located on the potential energy surface of the lowest energy singlet excited state. One planar minimum, slightly more stable, is associated with the emission at 335 nm, whereas the other one, with a bent etheno ring, is associated with the red-shifted emission.

Simulation of Radiation-Induced DNA Damage and Protection by Histones Using the Code RITRACKS

(1) Background: DNA damage is of great importance in the understanding of the effects of ionizing radiation. Various types of DNA damage can result from exposure to ionizing radiation, with clustered types considered the most important for radiobiological effects. (2) Methods: The code RITRACKS (Relativistic Ion Tracks), a program that simulates stochastic radiation track structures, was used to simulate DNA damage by photons and ions spanning a broad range of linear energy transfer (LET) values. To perform these simulations, the transport code was modified to include cross sections for the interactions of ions or electrons with DNA and amino acids for ionizations, dissociative electron attachment, and elastic collisions. The radiochemistry simulations were performed using a step-by-step algorithm that follows the evolution of all particles in time, including reactions between radicals and DNA structures and amino acids. Furthermore, detailed DNA damage events, such as base pair positions, DNA fragment lengths, and fragment yields, were recorded. (3) Results: We report simulation results using photons and the ions 1H+, 4He2+, 12C6+, 16O8+, and 56Fe26+ at various energies, covering LET values from 0.3 to 164 keV/µm, and performed a comparison with other codes and experimental results. The results show evidence of DNA protection from damage at its points of contacts with histone proteins. (4) Conclusions: RITRACKS can provide a framework for studying DNA damage from a variety of ionizing radiation sources with detailed representations of DNA at the atomic scale, DNA-associated proteins, and resulting DNA damage events and statistics, enabling a broader range of future comparisons with experiments such as those based on DNA sequencing.

XMECP: Reaching State-of-the-Art MECP Optimization in Multiscale Complex Systems

The Python-based program, XMECP, is developed for realizing robust, efficient, and state-of-the-art minimum energy crossing point (MECP) optimization in multiscale complex systems. This article introduces the basic capabilities of the XMECP program by theoretically investigating the MECP mechanism of several example systems including (1) the photosensitization mechanism of benzophenone, (2) photoinduced proton-coupled electron transfer in the cytosine-guanine base pair in DNA, (3) the spin-flip process in oxygen activation catalyzed by an iron-containing 2-oxoglutarate-dependent oxygenase (Fe/2OGX), and (4) the photochemical pathway of flavoprotein adjusted by the intensity of an external electric field. MECPs related to multistate reaction and multistate reactivity in large-scale complex biochemical systems can be well-treated by workflows suggested by the XMECP program. The branching plane updating the MECP optimization algorithm is strongly recommended as it provides derivative coupling vector (DCV) with explicit calculation and can equivalently evaluate contributions from non-QM residues to DCV, which can be nonadiabatic coupling or spin-orbit coupling in different cases. In the discussed QM/MM examples, we also found that the influence on the QM region by DCV can occur through noncovalent interactions and decay with distance. In the example of DNA base pairs, the nonadiabatic coupling occurs across the π-π stacking structure formed in the double-helix system. In contrast to general intuition, in the example of Fe/2OGX, the central ferrous and oxygen part contribute little to the spin-orbit coupling; however, a nearby arginine residue, which is treated by molecular mechanics in the QM/MM method, contributes significantly via two hydrogen bonds formed with α-ketoglutarate (α-KG). This indicates that the arginine residue plays a significant role in oxygen activation, driving the initial triplet state toward the productive quintet state, which is more than the previous knowledge that the arginine residue can bind α-KG at the reaction site by hydrogen bonds.

How a Single 5 eV Electron Can Induce Double-Strand Breaks in DNA: A Time-Dependent Density Functional Theory Study

Low-energy (<20 eV) electrons (LEEs) can resonantly interact with DNA to form transient anions (TAs) of fundamental units, inducing single-strand breaks (SSBs), and cluster damage, such as double-strand breaks (DSBs). Shape resonances, which arise from electron capture in a previously unfilled orbital, can induce only a SSB, whereas a single core-excited resonance (i.e., two electrons in excited orbitals of the field of a hole) has been shown experimentally to cause cluster lesions. Herein, we show from time-dependent density functional theory (TDDFT) that a core-excited resonance can produce a DSB, i.e., a single 5 eV electron can induce two close lesions in DNA. We considered the nucleotide with the G-C base pair (ds[5'-G-3']) as a model for electron localization in the DNA double helix and calculated the potential energy surfaces (PESs) of excited states of the ground-state TA of ds[5'-G-3'], which correspond to shape and core-excited resonances. The calculations show that shape TAs start at ca. 1 eV, while core-excited TAs occur only above 4 eV. The energy profile of each excited state and the corresponding PES are obtained by simultaneously stretching both C5'-O5' bonds of ds[5'-G-3']. From the nature of the PES, we find two dissociative (σ*) states localized on the PO4 groups at the C5' sites of ds[5'-G-3']. The first σ* state at 1 eV is due to a shape resonance, while the second σ* state is induced by a core-excited resonance at 5.4 eV. As the bond of the latter state stretches and arrives close to the dissociation limit, the added electron on C transfers to C5' phosphate, thus demonstrating the possibility of producing a DSB with only one electron of ca. 5 eV.

Interaction of low-energy electrons with radiosensitizers

We provide an experimentalist's perspective on the present state-of-the-art in the studies of low-energy electron interactions with common radiosensitizers, including compounds used in combined chemo-radiation therapy and their model systems. Low-energy electrons are important secondary species formed during the interaction of ionizing radiation with matter. Their role in the radiation chemistry of living organisms has become an important topic for more than 20 years. With the increasing number of works and reviews in the field, we would like to focus here on a very narrow area of compounds that have been shown to have radio-sensitizing properties on the one hand, and high reactivity towards low-energy electrons on the other hand. Gas phase experiments studying electron attachment to isolated molecules and environmental effects on reaction dynamics are reviewed for modified DNA components, nitroimidazoles, and organometallics. In the end, we provide a perspective on the future directions that may be important for transferring the fundamental knowledge about the processes induced by low-energy electrons into practice in the field of rational design of agents for concomitant chemo-radiation therapy.

Triplet Photochemistry of Effluent Organic Matter in Degradation of Extracellular Antibiotic Resistance Genes

Wastewater effluent is a major source of extracellular antibiotic resistance genes (eArGs) in the aquatic environment, a threat to human health and biosecurity. However, little is known about the extent to which organic matter in the wastewater effluent (EfOM) might contribute to photosensitized oxidation of eArGs. Triplet states of EfOM were found to dominate the degradation of eArGs (accounting for up to 85%). Photo-oxidation proceeded mainly via proton-coupled electron transfer reactions. They broke plasmid strands and damaged bases. O₂^•- was also involved, and it coupled with the reactions' intermediate radicals of eArGs. The second-order reaction rates of bla_TEM-1 and tet-A segments (209-216 bps) with the triplet state of 4-carboxybenzophenone were calculated to be (2.61-2.75) × 10⁸ M^-1 s^-1. Besides as photosensitizers, the antioxidant moieties in EfOM also acted as quenchers to revert intermediate radicals back to their original forms, reducing the rate of photodegradation. However, the terrestrial origin natural organic matter was unable to photosensitize because it formed less triplets, especially high-energy triplets, so its inhibitory effects predominated. This study advances our understanding of the role of EfOM in the photo-oxidation of eArGs and the difference between EfOM and terrestrial-origin natural organic matter.

Electron Attachment to Wobble Base Pairs

We have analyzed the low-energy electron attachment to wobble base pairs using the equation of motion coupled cluster method and extended basis sets. A doorway mechanism exists for the attachment of the additional electron to the base pairs, where the initially formed dipole-bound anion captures the incoming electron. The doorway dipole-bound anionic state subsequently leads to the formation of a valence-bound state, and the transfer of extra electron occurs by mixing of electronic and nuclear degrees of freedom. The formation of the valence-bound anion is associated with proton transfer in hypoxanthine-cytosine and hypoxanthine-adenine base pairs, which happens through a concerted electron-proton transfer process. The calculated rate constant for the dipole-bound to valence-bound transition in wobble base pairs is slower than that observed in the Watson-Crick guanine-cytosine base pair.

Proton transfer during DNA strand separation as a source of mutagenic guanine-cytosine tautomers

Proton transfer between the DNA bases can lead to mutagenic Guanine-Cytosine tautomers. Over the past several decades, a heated debate has emerged over the biological impact of tautomeric forms. Here, we determine that the energy required for generating tautomers radically changes during the separation of double-stranded DNA. Density Functional Theory calculations indicate that the double proton transfer in Guanine-Cytosine follows a sequential, step-like mechanism where the reaction barrier increases quasi-linearly with strand separation. These results point to increased stability of the tautomer when the DNA strands unzip as they enter the helicase, effectively trapping the tautomer population. In addition, molecular dynamics simulations indicate that the relevant strand separation time is two orders of magnitude quicker than previously thought. Our results demonstrate that the unwinding of DNA by the helicase could simultaneously slow the formation but significantly enhance the stability of tautomeric base pairs and provide a feasible pathway for spontaneous DNA mutations.

Electromagnetic bioeffects: a multiscale molecular simulation perspective

Electromagnetic bioeffects remain an enigma from both the experimental and theoretical perspectives despite the ubiquitous presence of related technologies in contemporary life. Multiscale computational modelling can provide valuable insights into biochemical systems and predict how they will be perturbed by external stimuli. At a microscopic level, it can be used to determine what (sub)molecular scale reactions various stimuli might induce; at a macroscopic level, it can be used to examine how these changes affect dynamic behaviour of essential molecules within the crowded biomolecular milieu in living tissues. In this review, we summarise and evaluate recent computational studies that examined the impact of externally applied electric and electromagnetic fields on biologically relevant molecular systems. First, we briefly outline the various methodological approaches that have been employed to study static and oscillating field effects across different time and length scales. The practical value of such modelling is then illustrated through representative case-studies that showcase the diverse effects of electric and electromagnetic field on the main physiological solvent - water, and the essential biomolecules - DNA, proteins, lipids, as well as some novel biomedically relevant nanomaterials. The implications and relevance of the theoretical multiscale modelling to practical applications in therapeutic medicine are also discussed. Finally, we summarise ongoing challenges and potential opportunities for theoretical modelling to advance the current understanding of electromagnetic bioeffects for their modulation and/or beneficial exploitation in biomedicine and industry.

Proton-Transfer Reactions in One-Electron-Oxidized G-Quadruplexes: A Density Functional Theory Study

Recently, G-quadruplexes (Gq) formed in B-DNA as secondary structures are found to be important therapeutic targets and material for developing nanodevices. Gq are guanine-rich and thus susceptible to oxidative damage by producing short-lived intermediate radicals via proton-transfer reactions. Understanding the mechanisms of radical formation in Gq is of fundamental interest to understand the early stages of DNA damage. Herein, we used density functional theory including aqueous phase (ωB97XD-PCM/6-31++G**) and considered single layer of Gq [G-quartets (G4): association of four guanines in a cyclic Hoogsteen hydrogen-bonded arrangement (Scheme 1)] to unravel the mechanisms of formation of intermediates by calculating the relative Gibbs free energies and spin density distributions of one-electron-oxidized G4 and its various proton-transfer states: G^•+, G(N₁-H)^•, G(N₂-H')^•, G(N₂-H″)^•, G(N₁-H)^•-(H⁺O₆)G, and G(N₂-H)^•-(H⁺N₇)G. The present calculation predicts the formation of G(N₂-H)^•-(H⁺N₇)G, which is only ca. 0.8 kcal/mol higher in energy than the initially formed G^•+. The formation of G(N₂-H)^•-(H⁺N₇)G plays a key role in explaining the formation of 8-OG along with G(N₁-H)^• formation via tautomerization from G(N₂-H)^•, as proposed recently.

Modulation of the Directionality of Hole Transfer between the Base and the Sugar-Phosphate Backbone in DNA with the Number of Sulfur Atoms in the Phosphate Group

This work shows that S atom substitution in phosphate controls the directionality of hole transfer processes between the base and sugar-phosphate backbone in DNA systems. The investigation combines synthesis, electron spin resonance (ESR) studies in supercooled homogeneous solution, pulse radiolysis in aqueous solution at ambient temperature, and density functional theory (DFT) calculations of in-house synthesized model compound dimethylphosphorothioate (DMTP(O^-)═S) and nucleotide (5'-O-methoxyphosphorothioyl-2'-deoxyguanosine (G-P(O^-)═S)). ESR investigations show that DMTP(O^-)═S reacts with Cl₂^•- to form the σ²σ*¹ adduct radical -P-S[Formula: see text]Cl, which subsequently reacts with DMTP(O^-)═S to produce [-P-S[Formula: see text]S-P-]^-. -P-S[Formula: see text]Cl in G-P(O^-)═S undergoes hole transfer to Gua, forming the cation radical (G^•+) via thermally activated hopping. However, pulse radiolysis measurements show that DMTP(O^-)═S forms the thiyl radical (-P-S^•) by one-electron oxidation, which did not produce [-P-S[Formula: see text]S-P-]^-. Gua in G-P(O^-)═S is oxidized unimolecularly by the -P-S^• intermediate in the sub-picosecond range. DFT thermochemical calculations explain the differences in ESR and pulse radiolysis results obtained at different temperatures.

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.

Electrophilic Properties of 2'-Deoxyadenosine···Thymine Dimer: Photoelectron Spectroscopy and DFT Studies

The anion radical of the 2'-deoxyadenosine···thymine (dAT^•-) pair has been investigated experimentally and theoretically in the gas phase. By employing negative-ion photoelectron spectroscopy (PES), we have registered a spectrum typical for the valence-bound anion, featuring a broad peak at the electron-binding energy (EBE) between ∼1.5 and 2.2 eV with the maximum at ∼1.9 eV. The measured value of the adiabatic electron affinity (AEA) for dAT was estimated to be ∼1.1 eV. Calculations performed at the M06-2X/6-31++G(d,p) level revealed that the structure, where thymine is coordinated to the sugar of dA by two hydrogen bonds, is responsible for the observed PES signal. The AEA_G and the vertical detachment energy of 0.91 and 1.68 eV, respectively, calculated for this structure reproduce the experimental values well. The role of the possible proton transfer in the stabilization of anionic radical complexes is discussed.

Hydrogen Bonding in Natural and Unnatural Base Pairs-A Local Vibrational Mode Study

In this work hydrogen bonding in a diverse set of 36 unnatural and the three natural Watson Crick base pairs adenine (A)-thymine (T), adenine (A)-uracil (U) and guanine (G)-cytosine (C) was assessed utilizing local vibrational force constants derived from the local mode analysis, originally introduced by Konkoli and Cremer as a unique bond strength measure based on vibrational spectroscopy. The local mode analysis was complemented by the topological analysis of the electronic density and the natural bond orbital analysis. The most interesting findings of our study are that (i) hydrogen bonding in Watson Crick base pairs is not exceptionally strong and (ii) the N-H⋯N is the most favorable hydrogen bond in both unnatural and natural base pairs while O-H⋯N/O bonds are the less favorable in unnatural base pairs and not found at all in natural base pairs. In addition, the important role of non-classical C-H⋯N/O bonds for the stabilization of base pairs was revealed, especially the role of C-H⋯O bonds in Watson Crick base pairs. Hydrogen bonding in Watson Crick base pairs modeled in the DNA via a QM/MM approach showed that the DNA environment increases the strength of the central N-H⋯N bond and the C-H⋯O bonds, and at the same time decreases the strength of the N-H⋯O bond. However, the general trends observed in the gas phase calculations remain unchanged. The new methodology presented and tested in this work provides the bioengineering community with an efficient design tool to assess and predict the type and strength of hydrogen bonding in artificial base pairs.

Electron-Induced Repair of 2'-Deoxyribose Sugar Radicals in DNA: A Density Functional Theory (DFT) Study

In this work, we used ωB97XD density functional and 6-31++G** basis set to study the structure, electron affinity, populations via Boltzmann distribution, and one-electron reduction potentials (E°) of 2′-deoxyribose sugar radicals in aqueous phase by considering 2′-deoxyguanosine and 2′-deoxythymidine as a model of DNA. The calculation predicted the relative stability of sugar radicals in the order C4′^• > C1′^• > C5′^• > C3′^• > C2′^•. The Boltzmann distribution populations based on the relative stability of the sugar radicals were not those found for ionizing radiation or OH-radical attack and are good evidence the kinetic mechanisms of the processes drive the products formed. The adiabatic electron affinities of these sugar radicals were in the range 2.6–3.3 eV which is higher than the canonical DNA bases. The sugar radicals reduction potentials (E°) without protonation (−1.8 to −1.2 V) were also significantly higher than the bases. Thus the sugar radicals will be far more readily reduced by solvated electrons than the DNA bases. In the aqueous phase, these one-electron reduced sugar radicals (anions) are protonated from solvent and thus are efficiently repaired via the “electron-induced proton transfer mechanism”. The calculation shows that, in comparison to efficient repair of sugar radicals by the electron-induced proton transfer mechanism, the repair of the cyclopurine lesion, 5′,8-cyclo-2′-dG, would involve a substantial barrier.

Recent progress in approximate quantum dynamics methods for the study of proton-coupled electron transfer reactions

Proton-coupled electron transfer (PCET) reactions are ubiquitous natural processes at the heart of energy conversion reactions in photosynthesis and respiration, DNA repair, and diverse enzymatic reactions. Theoretical formulation and computational method developments have eyed modeling of thermal and photoinduced PCET for the last three decades. The accumulation of these studies, collected in dozens of reviews, accounts, and perspectives, has firmly established the influence of quantum effects, including non-adiabatic electronic transitions, vibrational relaxation, zero-point energy, and proton tunneling, on the rate and mechanism of PCET reactions. Here, we focus on some recently-developed methods, spanning the last eight years, that can quantitatively capture these effects in the PCET context and provide efficient means for their qualitative description in complex systems. The theoretical background of each method and their accuracy with respect to exact results are discussed and the results of relevant PCET simulations based on each method are presented.

Theoretical Insight on the Biosensing Applications of 2D Materials

The research on the design of efficient, reliable, and cost-effective biosensors is expanding given its high demand in various fields such as health care, environmental surveillance, agriculture, diagnostics, industries, and so forth. In the last decade, various fascinating and interesting 2D materials with extraordinary properties have been experimentally synthesized and theoretically predicted. 2D materials have been explored for the sensing of different biomolecules because of their large surface area and strong interaction with different biomolecules. Theoretical simulations can bring important insight on the interaction of biomolecules on 2D materials, charge transfer, orbital interactions, and so forth and may play an important role in the development of efficient biosensors. Quantum simulation techniques, such as density functional theory (DFT), are very powerful and are gaining popularity especially with the advent of high-speed computing facilities. This review article provides theoretical insight regarding the interaction of various biomolecules on different 2D materials and the charge transfer between the biomolecules and 2D materials leading to electrochemical signals, which can then provide experimentalists the useful design parameters for fabrication of biosensors. It also includes an overview of quantum simulations, use of the DFT method for simulating biomolecules on 2D materials, parameters obtained from theoretical simulations and sensitivity, and limitations of computational techniques for sensing biomolecules on 2D materials. Furthermore, this review summarizes the recent work in first-principles investigation of 2D materials for the purpose of biomolecule sensing. Beyond the traditional graphene or 2D transition-metal dichalcogenides, some novel and recently proposed 2D materials such as pentagraphene, haeckelite, MXenes, and so forth which have exhibited good sensing applications have also been highlighted.

Analysis of protein-DNA interactions in chromatin by UV induced cross-linking and mass spectrometry

Protein–DNA interactions are key to the functionality and stability of the genome. Identification and mapping of protein–DNA interaction interfaces and sites is crucial for understanding DNA-dependent processes. Here, we present a workflow that allows mass spectrometric (MS) identification of proteins in direct contact with DNA in reconstituted and native chromatin after cross-linking by ultraviolet (UV) light. Our approach enables the determination of contact interfaces at amino-acid level. With the example of chromatin-associated protein SCML2 we show that our technique allows differentiation of nucleosome-binding interfaces in distinct states. By UV cross-linking of isolated nuclei we determined the cross-linking sites of several factors including chromatin-modifying enzymes, demonstrating that our workflow is not restricted to reconstituted materials. As our approach can distinguish between protein–RNA and DNA interactions in one single experiment, we project that it will be possible to obtain insights into chromatin and its regulation in the future.

Cross-linking mass spectrometry (XLMS) allows mapping of protein-protein and protein-RNA interactions, but the analysis of protein-DNA complexes remains challenging. Here, the authors develop a UV light-based XLMS workflow to determine protein-DNA interfaces in reconstituted chromatin and isolated nuclei.

Tautomers of homophthalic anhydride in the ground and excited electronic states: analysis through energy, hardness and vibrational signatures

The keto-enol tautomerisation in homophthalic anhydride (HA) is investigated in the ground (S₀) and excited (S₁) electronic states. The keto form with a dicarbonyl structure is found to be the most stable form in S₀ and enol form with a monocarbonyl structure in S₁ indicating an excited state intramolecular proton transfer (ESIPT) process. The computed results show consistency with the change in basis sets and methods of calculations. Apart from the two tautomers, transition states are also identified. The barrier to interconversion is found to reduce substantially in S₁. Internal reaction coordinate (IRC) calculations confirm the pathway of interconversion between the two forms in S₀ and S₁. The observed FT-IR spectra corroborate well with our computed spectra. The appearance of two strong lines around 1800 cm^-1 confirms the lowest energy structure to be the keto tautomer with a dicarbonyl form in S₀. Our computations corroborate well with the crystal structure data for an analogous molecule. Electron distribution in HOMO and LUMO indicate the excitation process as π → π^* in nature. The qualitative chemical concepts like hardness and electrophilicity are calculated to estimate the stability of the tautomers. The energy and hardness profiles with the variation of IRC are opposite to each other, verifying the principle of maximum hardness.

Vibrations of the guanine-cytosine pair in chloroform: an anharmonic computational study

We compute at the anharmonic level the vibrational spectra of the Watson-Crick dimer formed by guanosine (G) and cytidine (C) in chloroform, together with those of G, C and the most populated GG dimer. The spectra for deuterated and partially deuterated GC are also computed. We use DFT calculations, with B3LYP and CAM-B3LYP as reference functionals. Solvent effects from chloroform are included via the Polarizable Continuum Model (PCM), and by performing tests on models including up two chloroform molecules. Both B3LYP and CAM-B3LYP calculations reproduce the shape of the experimental spectra well in the fingerprint region (1500-1700 cm-1) and in the N-H stretching region (2800-3600 cm-1), with B3LYP providing better quantitative agreement with experiments. According to our calculations, the N-H amido streching mode of G falls at ∼2900 cm-1, while the N-H amino of G and C falls at ∼3100 cm-1 when hydrogen-bonded, or ∼3500 cm-1 when free. Overtone and combination bands strongly contribute to the absorption band at ∼3300 cm-1. Inclusion of bulk solvent effects significantly increases the accuracy of the computed spectra, while solute-solvent interactions have a smaller, though still noticeable, effect. Some key aspects of the anharmonic treatment of strongly vibrationally coupled supermolecular systems and the related methodological issues are also discussed.

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 (E0) 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 E0 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 E0 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.

Excited state hydrogen transfer dynamics in phenol-(NH3)2 studied by picosecond UV-near IR-UV time-resolved spectroscopy

Time-evolutions of excited state hydrogen transfer (ESHT) in phenol (PhOH)-(NH₃)₂ clusters have been measured by three-color picosecond (ps) ultraviolet (UV)-near infrared (NIR)-UV pump-probe ion dip spectroscopy. The formation of a reaction product, ˙NH₄NH₃, is detected by its NIR absorption due to a 3p-3s Rydberg transition. The ESHT reactions from all of the vibronic levels show biexponential time-evolutions, even from the S₁ origin. Based on the biexponential time-evolution, it is suggested that there is a second reaction path via the triplet πσ* state, which gives the slow component. The fast time-evolution of the ESHT reaction from the S₁ origin is measured to be 268 ps, which is 10-times slower than that in PhOH-(NH₃)₃, and a higher barrier between the ππ* and reactive πσ* states is suggested. The size dependence of the ESHT reaction rates is discussed based on a potential distortion due to the proton transferred state in the ππ* potential surface.

Alternative Approach for the Determination of Mean Free Paths of Electron Scattering in Liquid Water Based on Experimental Data

Mean free paths of low-energy electrons in liquid water are of importance for modeling many physicochemical processes, but neither theoretical predictions nor experimental results have converged for these parameters. We therefore introduce an approach to determine elastic and inelastic mean free paths (EMFP, IMFP) based on experimental data. We show that ab initio calculations of electron scattering with water clusters converge with cluster size, thus providing access to condensed-phase scattering. The results are used in Monte Carlo simulations to extract EMFP and IMFP from recent liquid-microjet experiments that determined the effective attenuation length (EAL) and the photoelectron angular distribution (PAD) following oxygen 1s-ionization of liquid water. For electron kinetic energies from 10 to 300 eV, we find that the IMFP is noticeably larger than the EAL. The EMFP is longer than that of gas-phase water and the IMFP is longer compared to latest theoretical estimations, but both EMFP and IMFP are much shorter than suggested by experimental measurements of integral cross sections for amorphous ice.

One-electron oxidation of TAT-motif triplex DNA and the ensuing Hoogsteen hydrogen-bonding dissociation

One-electron oxidation of adenine (A) leads initially to the formation of adenine radical cation (A^•+). Subsequent deprotonation of A^•+ can provoke deoxyribonucleic acid (DNA) damage, which further causes senescence, cancer formation, and even cell death. However, compared with considerable reports on A^•+ reactions in free deoxyadenosine (dA) and duplex DNA, studies in non-B-form DNA that play critical biological roles are rare at present. It is thus of vital importance to explore non-B-form DNA, among which the triplex is an emerging topic. Herein, we investigate the deprotonation behavior of A^•+ in the TAT triplex with continuous A bases by time-resolved laser flash photolysis. The rate constants for the one-oxidation of triplex 8.4 × 10⁸ M^-1 s^-1 and A^•+ deprotonation 1.3 × 10⁷ s^-1 are obtained. The kinetic isotope effect of A^•+ deprotonation in the TAT triplex is 1.8, which is characteristic of a direct release of the proton into the solvent similar to free base dA. It is thus elucidated that the A^•+ proton bound with the third strand is most likely to be released into the solvent because of the weaker Hoogsteen H-bonding interaction and the presence of the highly mobile hydration waters within the third strand. Additionally, it is confirmed through Fourier transform infrared spectroscopy that the deprotonation of A^•+ results in the dissociation of the third strand and disruption of the secondary structure of the triplex. These results provide valuable kinetic data and in-depth mechanistic insights for understanding the adenine oxidative DNA damage in the triplex.

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.

The Role of Proton Transfer on Mutations

Hydrogen bonds play a critical role in nucleobase studies as they encode genes, map protein structures, provide stability to the base pairs, and are involved in spontaneous and induced mutations. Proton transfer mechanism is a critical phenomenon that is related to the acid-base characteristics of the nucleobases in Watson-Crick base pairs. The energetic and dynamical behavior of the proton can be depicted from these characteristics and their adjustment to the water molecules or the surrounding ions. Further, new pathways open up in which protonated nucleobases are generated by proton transfer from the ionized water molecules and elimination of a hydroxyl radical in this review, the analysis will be focused on understanding the mechanism of untargeted mutations in canonical, wobble, Hoogsteen pairs, and mutagenic tautomers through the non-covalent interactions. Further, rare tautomer formation through the single proton transfer (SPT) and the double proton transfer (DPT), quantum tunneling in nucleobases, radiation-induced bystander effects, role of water in proton transfer (PT) reactions, PT in anticancer drugs-DNA interaction, displacement and oriental polarization, possible models for mutations in DNA, genome instability, and role of proton transfer using kinetic parameters for RNA will be discussed.

Reaction of Electrons with DNA: Radiation Damage to Radiosensitization

This review article provides a concise overview of electron involvement in DNA radiation damage. The review begins with the various states of radiation-produced electrons: Secondary electrons (SE), low energy electrons (LEE), electrons at near zero kinetic energy in water (quasi-free electrons, (e-qf)) electrons in the process of solvation in water (presolvated electrons, e-pre), and fully solvated electrons (e-aq). A current summary of the structure of e-aq, and its reactions with DNA-model systems is presented. Theoretical works on reduction potentials of DNA-bases were found to be in agreement with experiments. This review points out the proposed role of LEE-induced frank DNA-strand breaks in ion-beam irradiated DNA. The final section presents radiation-produced electron-mediated site-specific formation of oxidative neutral aminyl radicals from azidonucleosides and the evidence of radiosensitization provided by these aminyl radicals in azidonucleoside-incorporated breast cancer cells.

Modeling the action of environment on proton tunneling in the adenine-thymine base pair

Effect of environment coupling on the quantum-biological phenomenon of proton tunneling in the hydrogen bonds of the adenine-thymine base pair in DNA was modeled within the framework of quantum statistics and perturbation theory. A number of important thermodynamic indicators including partition function, free energy, and entropy were then calculated and examined. The proton was then assumed to be subject to an attraction represented by a double-well potential energy surface with a small asymmetry, which was considered as the perturbation introduced to the system. The action of environment manifested itself in the form of a global minimum in the free energy curve, as an implicit implication of the tendency of the system toward randomness and disorder, at which no spontaneous change such as quantum tunneling will accordingly occur. Furthermore, assuming the free energy to be in a close neighborhood of its minimum truly explained the smallness of the contribution of environment coupling to the tunneling probability reported in the literature based on the fact that the closer the free energy to its minimum, the less the transition probability to this point.

Excited States of One-Electron Oxidized Guanine-Cytosine Base Pair Radicals: A Time Dependent Density Functional Theory Study

One-electron oxidized guanine (G^•+) in DNA generates several short-lived intermediate radicals via proton transfer reactions resulting in the formation of neutral guanine radicals. The identification of these radicals in DNA is of fundamental interest to understand the early stages of DNA damage. Herein, we used time-dependent density functional theory (TD-ωB97XD-PCM/6-31G(3df,p)) to calculate the vertical excitation energies of one-electron oxidized G and G-cytosine (C) base pair in various protonation states: G^•+, G(N1-H)^•, and G(N2-H)^•, as well as G^•+-C, G(N1-H)^•-(H⁺)C, G(N1-H)^•-(N4-H⁺)C), G(N1-H)^•-C, and G(N2-H)^•-C in aqueous phase. The calculated UV-vis spectra of these radicals are in good agreement with the experiment for the G radical species when the calculated values are red-shifted by 40-70 nm. The present calculations show that the lowest energy transitions of proton transfer species (G(N1-H)^•-(H⁺)C, G(N1-H)^•-(N4-H⁺)C, and G(N1-H)^•-C) are substantially red-shifted in comparison to the spectrum of G^•+-C. The calculated spectrum of G(N2-H)^•-C shows intense absorption (high oscillator strength), which matches the strong absorption in the experimental spectra of G(N2-H)^• at 600 nm. The present calculations predict the lowest charge transfer transition of C → G^•+ is π → π* in nature and lies in the UV region (3.4-4.3 eV) with small oscillator strength.

Intramolecular electron-induced proton transfer and its correlation with excited-state intramolecular proton transfer

Electron-induced proton transfer depicts the proton motion coupled with the attachment of a low-energy electron to a molecule, which helps to understand copious fundamental chemical processes. Intramolecular electron-induced proton transfer is a similar process that occurs within a single molecule. To date, there is only one known intramolecular example, to the best of our knowledge. By studying the 10-hydroxybenzo[h]quinoline and 8-hydroxyquinoline molecules using anion photoelectron spectroscopy and density functional theory, and by theoretical screening of six other molecules, here we show the intramolecular electron-induced proton transfer capability of a long list of molecules that meanwhile have the excited-state intramolecular proton transfer property. Careful examination of the intrinsic electronic signatures of these molecules reveals that these two distinct processes should occur to the same category of molecules. Intramolecular electron-induced proton transfer could have potential applications such as molecular devices that are responsive to electrons or current.

Effects of Positive and Negative Ionization on Prototropy in Pyrimidine Bases: An Unusual Case of Isocytosine

Intramolecular proton-transfers (prototropic conversions) have been studied for the guanine building block isocytosine (iC), and effects of positive ionization, called one-electron oxidation (iC - e → iC^+•), and negative ionization, called one-electron reduction (iC + e → iC^-•), on tautomeric conversions when proceeding from neutral to ionized isocytosine have been discussed. Although radical cations and radical anions are very short-lived species, the ionization effects could be investigated by quantum-chemical methods. Such kind of studies gives some information about the labile protons and the most basic positions in the neutral and radical forms of the tautomeric system. For investigations, the complete isomeric mixture of isocytosine has been considered and calculations performed in two extreme environments, apolar {DFT(B3LYP)/6-311+G(d,p)} and polar {PCM(water)//DFT(B3LYP)/6-311+G(d,p)}. For selected isomers, the G4 theory has also been applied. There are no good relations for energetic parameters of neutral and ionized forms. Ionization energies depend on localization of labile protons. Tautomeric equilibria for neutral and ionized isocytosine, favored sites of protonation and deprotonation, and favored structures of protonated and deprotonated forms strongly depend on environment. Acidity of iC^+• is close to that of the iC conjugate acid, and basicity of iC^-• is close to that of the iC conjugate base. This increase of acid-base properties of charged radicals explains the proton-transfer in ionized pairs of nucleobases. When compared to other pyrimidine bases such as uracil (U) and cytosine (C), which exhibit analogous tautomeric equilibria between nine prototropic tautomers as isocytosine, the tautomeric preferences for iC, iC^+•, iC^-•, U, U^+•, U^-•, C, C^+•, and C^-• are completely different. The differences suggest that acid-base properties of functional groups, their stabilities, and ionization energies play a principal role in proton-transfers for pyrimidine bases and influence compositions of tautomeric mixtures.

Unveiling the complex vibronic structure of the canonical adenine cation

Adenine, a DNA base, exists as several tautomers and isomers that are closely lying in energy and that may form a mixture upon vaporization of solid adenine. Indeed, it is challenging to bring adenine into the gas phase, especially as a unique tautomer. The experimental conditions were tuned to prepare a jet-cooled canonical adenine (9H-adenine). This isolated DNA base was ionized by single VUV photons from a synchrotron beamline and the corresponding slow photoelectron spectrum was compared to ab initio computations of the neutral and ionic species. We report the vibronic structure of the X+ 2A'' (D0), A+ 2A' (D1) and B+ 2A'' (D2) electronic states of the 9H adenine cation, from the adiabatic ionization energy (AIE) up to AIE + 1.8 eV. Accurate AIEs are derived for the 9H-adenine (X[combining tilde] 1A') + hν → 9H-adenine+ (X+ 2A'', A+ 2A', B+ 2A'') + e- transitions. Close to the AIE, we fully assign the rich vibronic structure solely to the 9H-adenine (X 1A') + hν → 9H-adenine+ (X+ 2A'') transition. Importantly, we show that the lowest cationic electronic states of canonical adenine are coupled vibronically. The present findings are important for understanding the effects of ionizing radiation and the charge distribution on this elementary building block of life, at ultrafast, short, and long timescales.

Communication: Site-selective bond excision of adenine upon electron transfer

This work demonstrates that selective excision of hydrogen atoms at a particular site of the DNA base adenine can be achieved in collisions with electronegative atoms by controlling the impact energy. The result is based on analysing the time-of-flight mass spectra yields of potassium collisions with a series of labeled adenine derivatives. The production of dehydrogenated parent anions is consistent with neutral H loss either from selective breaking of C-H or N-H bonds. These unprecedented results open up a new methodology in charge transfer collisions that can initiate selective reactivity as a key process in chemical reactions that are dominant in different areas of science and technology.