Electron Attachment to DNA Base Pairs: An Interplay of Dipole- and Valence-Bound States

How Accurately Can DFT Describe Non-valence Anions?

Weakly bound non-valence anions are molecular systems where the excess electron stabilizes in a very diffuse orbital whose size, shape, and binding energy (∼1-100 meV) are governed by the long-range electrostatic potential of the molecule. Its binding energy comes mainly from charge-dipole or charge-multipole interactions or dispersion forces. While highly correlated methods, like coupled cluster methods, are considered to be the state of the art for describing anionic systems, especially when the electron lies in a very diffuse orbital, we consider here the possibility to use DFT-based calculations. In such molecular anions, the outer electron experiences long-range exchange and correlation interactions. We show that DFT can describe long-range bound states provided that a correct asymptotic exchange and correlation potential is used, namely, that from a range-separated hybrid functional. This opens an alternative to the computationally demanding highly correlated method calculations. It is also suggested that the study of weakly bound anions could help in the construction of new DFT potentials to study systems where nonlocal effects are significant.

Secondary Electron Attachment-Induced Radiation Damage to Genetic Materials

Reactions of radiation-produced secondary electrons (SEs) with biomacromolecules (e.g., DNA) are considered one of the primary causes of radiation-induced cell death. In this Review, we summarize the latest developments in the modeling of SE attachment-induced radiation damage. The initial attachment of electrons to genetic materials has traditionally been attributed to the temporary bound or resonance states. Recent studies have, however, indicated an alternative possibility with two steps. First, the dipole-bound states act as a doorway for electron capture. Subsequently, the electron gets transferred to the valence-bound state, in which the electron is localized on the nucleobase. The transfer from the dipole-bound to valence-bound state happens through a mixing of electronic and nuclear degrees of freedom. In the presence of aqueous media, the water-bound states act as the doorway state, which is similar to that of the presolvated electron. Electron transfer from the initial doorway state to the nucleobase-bound state in the presence of bulk aqueous media happens on an ultrafast time scale, and it can account for the decrease in DNA strand breaks in aqueous environments. Analyses of the theoretically obtained results along with experimental data have also been discussed.

Electron Attachment to DNA: The Protective Role of Amino Acids

We have studied the effect of amino acids on the electron attachment properties of a DNA nucleobase, with cytosine as a model system. The equation of motion coupled cluster theory with an extended basis set has been used to simulate the electron-attached state of the DNA model system. Arginine, alanine, lysine, and glycine are the four amino acids considered to investigate their role in electron attachment to a DNA nucleobase. The electron attachment to cytosine in all the four cytosine-amino acid gas-phase dimer complexes follows a doorway mechanism, where the electron gets transferred from the initial dipole-bound doorway state to the final nucleobase-bound state through the mixing of electronic and nuclear degrees of freedom. When cytosine is bulk-solvated with glycine, the glycine-bound state acts as the doorway state, where the initial electron density is localized on the bulk amino acid and away from the nucleobase, thus leading to the physical shielding of the nucleobase from the incoming electron. At the same time, the presence of amino acids can increase the stability of the nucleobase-bound anionic state, which can suppress the sugar-phosphate bond rupture caused by dissociative electron attachment to DNA.

Recent Advances in the Study of Negative-Ion Resonances Using Multiconfigurational Propagator and a Complex Absorbing Potential

The transient resonances are a challenge to bound state quantum mechanics. These states lie in the continuum part of the spectrum of the Hamiltonian. For this, one has to treat a continuum problem due to electron-molecule scattering and the many-electron correlation problem simultaneously. Moreover, the description of a resonance requires a wavefunction that bridges the part that resembles a bound state with another that resembles a continuum state such that the continuity of the wavefunction and its first derivative with respect to the distance between the incoming projectile and the target is maintained. A review of the recent advances in the theoretical investigation of the negative-ion resonances (NIR) is presented. The NIRs are ubiquitous in nature. They result from the scattering of electrons off of an atomic or molecular target. They are important for numerous chemical processes in upper atmosphere, space and even biological systems. A contextual background of the existing theoretical methods as well as the newly-developed multiconfigurational propagator tools based on a complex absorbing potential are discussed.

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.

Experimental Observation of the Resonant Doorways to Anion Chemistry: Dynamic Role of Dipole-Bound Feshbach Resonances in Dissociative Electron Attachment

Anion chemical dynamics of autodetachment and fragmentation mediated by the dipole-bound state (DBS) have been thoroughly investigated in a state-specific way by employing the picosecond time-resolved or the nanosecond frequency-resolved spectroscopy combined with the cryogenically cooled ion trap and velocity-map imaging techniques. For the ortho-, meta-, or para-iodophenoxide anion (o-, m-, or p-IPhO^-), the C-I bond rupture occurs via the nonadiabatic transition from the DBS to the nearby valence-bound states (VBS) of the anion where the vibronic coupling into the S₁ (πσ*) state (repulsive along the C-I bond extension coordinate) should be largely responsible. Dynamic details are governed by the isomer-specific nature of the potential energy surfaces in the vicinity of the DBS-VBS curve crossings, as manifested in the huge different chemical reactivity of o-, m-, or p-IPhO^-. It is confirmed here that the C-I bond dissociation is mediated by DBS resonances, providing the foremost evidence that the metastable DBS plays the critical role as the doorway into the anion chemistry especially of the dissociative electron attachment (DEA). The fragmentation channel is dominant when it is mediated by the DBS resonances located below the electron-affinity (EA) threshold, whereas it is kinetically adjusted by the competitive autodetachment when the DBS resonances above EA convey the electron to the valence orbitals. The product yield of the C-I bond cleavage is strongly mode-dependent as the rate of the concomitant autodetachment is much influenced by the characteristics of the individual vibrational modes, paving a new way of the reaction control of the anion chemistry.

Pair Natural Orbital Equation-of-Motion Coupled-Cluster Method for Core Binding Energies: Theory, Implementation, and Benchmark

We present the theory and implementation of a lower scaling core-valence separated equation-of-motion coupled-cluster approach based on domain-based local pair natural orbitals for core binding energies. The accuracy of the new method has been compared with that of the standard equation-of-motion coupled-cluster method and experimentally measured results. The use of pair natural orbitals significantly reduces the computation cost and can be applied to large molecules.

Dynamic role of the correlation effect revealed in the exceptionally slow autodetachment rates of the vibrational Feshbach resonances in the dipole-bound state

Real-time autodetachment dynamics of the loosely bound excess electron from the vibrational Feshbach resonances of the dipole-bound states (DBS) of 4-bromophonoxide (4-BrPhO^-) and 4-chlorophenoxide (4-ClPhO^-) anions have been thoroughly investigated. The state-specific autodetachment rate measurements obtained by the picosecond time-resolved pump-probe method on the cryogenically cooled anions exhibit an exceptionally long lifetime (τ) of ∼823 ± 156 ps for the 11'¹ vibrational mode of the 4-BrPhO^- DBS. Strong mode-dependency in the wide dynamic range has also been found, giving τ ∼ 5.3 ps for the 10'¹ mode, for instance. Though it is nontrivial to get the state-specific rates for the 4-ClPhO^- DBS, the average autodetachment lifetime of the 19'¹20'¹/11'¹ mode has been estimated to be ∼548 ± 108 ps. Observation of these exceptionally slow autodetachment rates of vibrational Feshbach resonances strongly indicates that the correlation effect may play a significant role in the DBS photodetachment dynamics. Fermi's golden rule has been invoked so that the correlation effect is taken into account in the form of the interaction between the charge and the induced dipole where the latter is given by the polarizable counterparts of the electron-rich halogenated compound and the diffuse non-valence electron. This report suggests that one may measure, from the real-time autodetachment dynamics, the extent of the correlation effect contribution to the stabilization and/or dynamics of the excess non-valence electron among many different types of long-range interactions of the DBS.

Doorway Mechanism for Electron Attachment Induced DNA Strand Breaks

We report a new doorway mechanism for the dissociative electron attachment to genetic materials. The dipole-bound state of the nucleotide anion acts as the doorway for electron capture in the genetic material. The electron gets subsequently transferred to a dissociative σ*-type anionic state localized on a sugar-phosphate or a sugar-nucleobase bond, leading to their cleavage. The electron transfer is mediated by the mixing of electronic and nuclear degrees of freedom. The cleavage rate of the sugar-phosphate bond predicted by this new mechanism is higher than that of the sugar-nucleobase bond breaking, and both processes are considerably slower than the formation of a stable valence-bound anion. The new mechanism can explain the relative rates of electron attachment induced bond cleavages in genetic materials.

Intermediate state representation approach to physical properties of molecular electron-attached states: Theory, implementation, and benchmarking

Computational schemes for comprehensive studies of molecular electron-attached states and the calculation of electron affinities (EAs) are formulated and implemented employing the intermediate state representation (ISR) formalism and the algebraic-diagrammatic construction approximation for the electron propagator (EA-ADC). These EA-ADC(n)/ISR(m) schemes allow for a consistent treatment of not only electron affinities and pole strengths up to third-order of perturbation theory (n = 3) but also one-electron properties of electron-attached states up to second order (m = 2). The EA-ADC/ISR equations were implemented in the Q-Chem program for Ŝ_z-adapted intermediate states, allowing also open-shell systems to be studied using unrestricted Hartree-Fock references. For benchmarking of the EA-(U)ADC/ISR schemes, EAs and dipole moments of various electron-attached states of small closed- and open-shell molecules were computed and compared to full configuration interaction data. As an illustrative example, EA-ADC(3)/ISR(2) has been applied to the thymine-thymine (6-4) DNA photolesion.

Electron Attachment to Cytosine: The Role of Water

We present an EOM-CCSD-based quantum mechanical/molecular mechanical (QM/MM) study on the electron attachment process to solvated cytosine. The electron attachment in the bulk solvated cytosine occurs through a doorway mechanism, where the initial electron is localized on water. The electron is subsequently transferred to cytosine by the mixing of electronic and nuclear degrees of freedom, which occurs on an ultrafast time scale. The bulk water environment stabilizes the cytosine-bound anion by an extensive hydrogen-bond network and drastically enhances the electron transfer rate from that observed in the gas phase. Microhydration studies cannot reproduce the effect of the bulk water environment on the electron attachment process, and one needs to include a large number of water molecules in the calculation to obtain converged results. The predicted adiabatic electron affinity and electron transfer rate obtained from our QM/MM calculations are consistent with the available experimental results.

Efficient implementation of the single-reference algebraic diagrammatic construction theory for charged excitations: Applications to the TEMPO radical and DNA base pairs

We present an efficient implementation of the second- and third-order single-reference algebraic diagrammatic construction (ADC) theory for electron attachment and ionization energies and spectra [EA/IP-ADC(n), n = 2, 3]. Our new EA/IP-ADC program features spin adaptation for closed-shell systems, density fitting for efficient handling of the two-electron integral tensors, and vectorized and parallel implementation of tensor contractions. We demonstrate capabilities of our efficient implementation by applying the EA/IP-ADC(n) (n = 2, 3) methods to compute the photoelectron spectrum of the (2,2,6,6-tetramethylpiperidin-1-yl)oxyl (TEMPO) radical, as well as the vertical and adiabatic electron affinities of TEMPO and two DNA base pairs (guanine-cytosine and adenine-thymine). The spectra and electron affinities computed using large diffuse basis sets with up to 1028 molecular orbitals are found to be in good agreement with the best available results from the experiment and theoretical simulations.

Fate of Dipole-Bound Anion States when Hydrated

Many strongly polar molecules can form an anion by attaching an electron to either an empty or half-filled valence-bound (VB) orbital or a so-called dipole-bound (DB) orbital. These two families of orbitals can be very different in their radial extent (the former are usually more compact, while the latter are quite diffuse) and in the degree to which they are affected by surrounding solvent molecules. In this study, the effects of hydration (representative of strong solvation) on the DB state of a model polar species are investigated with an eye toward determining whether this state is stabilized or even persists when a few to 100 water molecules surround the polar molecule. It is found that in the presence of up to ca. 10-12 water molecules, the excess electron can remain in a DB orbital. However, once there are enough water molecules to form a complete first hydration shell (or more), the excess electron migrates into an orbital localized on the outer surface of the water solvent cage. These findings have implications on the possible role of DB states as doorways to facilitating electron attachment and subsequent electron transfer to VB states. It is shown that even when the electron is bound to the surface of the surrounding solvent, the dipole potential of the solute molecule can influence where on the surface the electron binds. It is also illustrated that using continuum dielectric methods to describe the hydration of DB states is fraught with danger because much of the outermost electron density in such states penetrates outside the boundary of the cavity used in these methods.