Excited state proton-coupled electron transfer in 8-oxoG-C and 8-oxoG-A base pairs: a time dependent density functional theory (TD-DFT) study

An Update on Molecularly Imprinted Polymer Design through a Computational Approach to Produce Molecular Recognition Material with Enhanced Analytical Performance

Molecularly imprinted polymer (MIP) computational design is expected to become a routine technique prior to synthesis to produce polymers with high affinity and selectivity towards target molecules. Furthermore, using these simulations reduces the cost of optimizing polymerization composition. There are several computational methods used in MIP fabrication and each requires a comprehensive study in order to select a process with results that are most similar to properties exhibited by polymers synthesized through laboratory experiments. Until now, no review has linked computational strategies with experimental results, which are needed to determine the method that is most appropriate for use in designing MIP with high molecular recognition. This review will present an update of the computational approaches started from 2016 until now on quantum mechanics, molecular mechanics and molecular dynamics that have been widely used. It will also discuss the linear correlation between computational results and the polymer performance tests through laboratory experiments to examine to what extent these methods can be relied upon to obtain polymers with high molecular recognition. Based on the literature search, density functional theory (DFT) with various hybrid functions and basis sets is most often used as a theoretical method to provide a shorter MIP manufacturing process as well as good analytical performance as recognition material.

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.

Role of Electron-Driven Proton-Transfer Processes in the Ultrafast Deactivation of Photoexcited Anionic 8-oxoGuanine-Adenine and 8-oxoGuanine-Cytosine Base Pairs

It has been reported that 8-oxo-7,8-dihydro-guanosine (8-oxo-G), which is the main product of oxidative damage of DNA, can repair cyclobutane pyrimidine dimer (CPD) lesions when incorporated into DNA or RNA strands in proximity to such lesions. It has therefore been suggested that the 8-oxo-G nucleoside may have been a primordial precursor of present-day flavins in DNA or RNA repair. Because the electron transfer leading to the splitting of a thymine-thymine pair in a CPD lesion occurs in the photoexcited state, a reasonably long excited-state lifetime of 8-oxo-G is required. The neutral (protonated) form of 8-oxo-G exhibits a very short (sub-picosecond) intrinsic excited-state lifetime which is unfavorable for repair. It has therefore been argued that the anionic (deprotonated) form of 8-oxo-G, which exhibits a much longer excited-state lifetime, is more likely to be a suitable cofactor for DNA repair. Herein, we have investigated the exited-state quenching mechanisms in the hydrogen-bonded complexes of deprotonated 8-oxo-G^- with adenine (A) and cytosine (C) using ab initio wave-function-based electronic-structure calculations. The calculated reaction paths and potential-energy profiles reveal the existence of barrierless electron-driven inter-base proton-transfer reactions which lead to low-lying S₁/S₀ conical intersections. The latter can promote ultrafast excited-state deactivation of the anionic base pairs. While the isolated deprotonated 8-oxo-G^- nucleoside may have been an efficient primordial repair cofactor, the excited states of the 8-oxo-G^--A and 8-oxo-G^--C base pairs are likely too short-lived to be efficient electron-transfer repair agents.

Do Solvated Electrons (e(aq)⁻) Reduce DNA Bases? A Gaussian 4 and Density Functional Theory-Molecular Dynamics Study

The solvated electron (e(aq)⁻) is a primary intermediate after an ionization event that produces reductive DNA damage. Accurate standard redox potentials (E(o)) of nucleobases and of e(aq)⁻ determine the extent of reaction of e(aq)⁻ with nucleobases. In this work, E(o) values of e(aq)⁻ and of nucleobases have been calculated employing the accurate ab initio Gaussian 4 theory including the polarizable continuum model (PCM). The Gaussian 4-calculated E(o) of e(aq)⁻ (-2.86 V) is in excellent agreement with the experimental one (-2.87 V). The Gaussian 4-calculated E(o) of nucleobases in dimethylformamide (DMF) lie in the range (-2.36 V to -2.86 V); they are in reasonable agreement with the experimental E(o) in DMF and have a mean unsigned error (MUE) = 0.22 V. However, inclusion of specific water molecules reduces this error significantly (MUE = 0.07). With the use of a model of e(aq)⁻ nucleobase complex with six water molecules, the reaction of e(aq)⁻ with the adjacent nucleobase is investigated using approximate ab initio molecular dynamics (MD) simulations including PCM. Our MD simulations show that e(aq)⁻ transfers to uracil, thymine, cytosine, and adenine, within 10 to 120 fs and e(aq)⁻ reacts with guanine only when a water molecule forms a hydrogen bond to O6 of guanine which stabilizes the anion radical.

Description of Conical Intersections with Density Functional Methods

Conical intersections are perhaps the most significant mechanistic features of chemical reactions occurring through excited states. By providing funnels for efficient non-adiabatic population transfer, conical intersections govern the branching ratio of products of such reactions, similar to what the transition states do for ground-state reactivity. In this regard, intersections between the ground and the lowest excited states play a special role, and the correct description of the potential energy surfaces in their vicinity is crucial for understanding the mechanism and dynamics of excited-state reactions. The methods of density functional theory, such as time-dependent density functional theory, are widely used to describe the excited states of large molecules. However, are these methods suitable for describing the conical intersections or do they lead to artifacts and, consequently, to erroneous description of reaction dynamics? Here we address the first part of this question and analyze the ability of several density functional approaches, including the linear-response time-dependent approach as well as the spin-flip and ensemble formalisms, to provide the correct description of conical intersections and the potential energy surfaces in their vicinity. It is demonstrated that the commonly used linear-response time-dependent theory does not yield a proper description of these features and that one should instead use alternative computational approaches.

Ultrafast excited-state dynamics and vibrational cooling of 8-oxo-7,8-dihydro-2'-deoxyguanosine in D2O

Nguyen and Burrows recently demonstrated that UV-B irradiation of 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodG), a signature product of oxidatively damaged DNA, can repair cyclobutane pyrimidine dimers in double-stranded DNA (J. Am. Chem. Soc. 2011, 133, 14586 - 14589). In order to test the hypothesis that repair occurs by photoinduced electron transfer, it is critical to determine basic photophysical parameters of 8-oxodG including the excited-state lifetime. Here, femtosecond transient absorption spectroscopy was used to study the ultrafast excited-state dynamics of 8-oxodG with excitation in the UV and probing at visible and mid-IR wavelengths. The excited-state lifetimes of both neutral and basic forms of 8-oxodG in D2O are reported for the first time by monitoring the disappearance of excited-state absorption at 570 nm. The lifetime of the first excited state of the neutral form is 0.9 ± 0.1 ps, or nearly twice as long as that of 2'-deoxyguanosine. The basic form of 8-oxodG exhibits a much longer excited-state lifetime of 43 ± 3 ps. Following ultrafast internal conversion by neutral 8-oxodG, a vibrationally hot ground state is created that dissipates its excess vibrational energy to the solvent on a time scale of 2.4 ± 0.4 ps. Femtosecond time-resolved IR experiments provide additional insights into excited-state dynamics and the vibrational relaxation of several modes in the fingerprint region.

π- vs σ-radical states of one-electron-oxidized DNA/RNA bases: a density functional theory study

As a result of their inherent planarity, DNA base radicals generated by one-electron oxidation/reduction or bond cleavage form π- or σ-radicals. While most DNA base systems form π-radicals, there are a number of nucleobase analogues such as one-electron-oxidized 6-azauraci1, 6-azacytosine, and 2-thiothymine or one-electron reduced 5-bromouracil that form more reactive σ-radicals. Elucidating the availability of these states within DNA, base radical electronic structure is important to the understanding of the reactivity of DNA base radicals in different environments. In this work, we address this question by the calculation of the relative energies of π- and σ-radical states in DNA/RNA bases and their analogues. We used density functional theory B3LYP/6-31++G** method to optimize the geometries of π- and σ-radicals in Cs symmetry (i.e., planar) in the gas phase and in solution using the polarized continuum model (PCM). The calculations predict that σ- and π-radical states in one-electron-oxidized bases of thymine, T(N3-H)(•), and uracil, U(N3-H)(•), are very close in energy; i.e., the π-radical is only ca. 4 kcal/mol more stable than the σ-radical. For the one-electron-oxidized radicals of cytosine, C(•+), C(N4-H)(•), adenine, A(•+), A(N6-H)(•), and guanine, G(•+), G(N2-H)(•), G(N1-H)(•), the π-radicals are ca. 16-41 kcal/mol more stable than their corresponding σ-radicals. Inclusion of solvent (PCM) is found to stabilize the π- over σ-radical of each of the systems. U(N3-H)(•) with three discrete water molecules in the gas phase is found to form a three-electron σ bond between the N3 atom of uracil and the O atom of a water molecule, but on inclusion of full solvation and discrete hydration, the π-radical remains most stable.