Sequential Electron Transfer Governs the UV-Induced Self-Repair of DNA Photolesions

Total Page:16

File Type:pdf, Size:1020Kb

Sequential Electron Transfer Governs the UV-Induced Self-Repair of DNA Photolesions Edinburgh Research Explorer Sequential electron transfer governs the UV-induced self-repair of DNA photolesions Citation for published version: Szabla, R, Kruse, H, Stadlbauer, P, Šponer, J & Sobolewski, AL 2018, 'Sequential electron transfer governs the UV-induced self-repair of DNA photolesions', Chemical Science, vol. 9, no. 12, pp. 3131-3140. https://doi.org/10.1039/c8sc00024g Digital Object Identifier (DOI): 10.1039/c8sc00024g Link: Link to publication record in Edinburgh Research Explorer Document Version: Publisher's PDF, also known as Version of record Published In: Chemical Science General rights Copyright for the publications made accessible via the Edinburgh Research Explorer is retained by the author(s) and / or other copyright owners and it is a condition of accessing these publications that users recognise and abide by the legal requirements associated with these rights. Take down policy The University of Edinburgh has made every reasonable effort to ensure that Edinburgh Research Explorer content complies with UK legislation. If you believe that the public display of this file breaches copyright please contact [email protected] providing details, and we will remove access to the work immediately and investigate your claim. Download date: 28. Sep. 2021 Chemical Science View Article Online EDGE ARTICLE View Journal | View Issue Sequential electron transfer governs the UV- induced self-repair of DNA photolesions† Cite this: Chem. Sci.,2018,9,3131 ab b bc b Rafał Szabla, * Holger Kruse, Petr Stadlbauer, Jiˇr´ı Sponerˇ and Andrzej L. Sobolewskia Cyclobutane pyrimidine dimers (CpDs) are among the most common DNA lesions occurring due to the interaction with ultraviolet light. While photolyases have been well known as external factors repairing CpDs, the intrinsic self-repairing capabilities of the GAT]T DNA sequence were discovered only recently and are still largely obscure. Here, we elucidate the mechanistic details of this self-repair process by means of MD simulations and QM/MM computations involving the algebraic diagrammatic construction to the second order [ADC(2)] method. We show that local UV-excitation of guanine may be followed by up to three subsequent electron transfers, which may eventually enable efficient CpD ring opening when the negative charge resides on the T]T dimer. Consequently, the molecular mechanism of GAT]T self- repair can be envisaged as sequential electron transfer (SET) occurring downhill along the slope of the S1 Creative Commons Attribution-NonCommercial 3.0 Unported Licence. potential energy surface. Even though the general features of the SET mechanism are retained in both of the studied stacked conformers, our optimizations of different S1/S0 state crossings revealed minor differences which could influence their self-repair efficiencies. We expect that such assessment of the Received 3rd January 2018 availability and efficiency of the SET process in other DNA oligomers could hint towards other sequences Accepted 22nd February 2018 exhibiting similar photochemical properties. Such explorations will be particularly fascinating in the DOI: 10.1039/c8sc00024g context of the origins of biomolecules on Earth, owing to the lack of external repairing factors in the rsc.li/chemical-science Archean age. This article is licensed under a Introduction formation and repair have been attracting signicant attention involving both experimental and theoretical studies.8–17 Dimerization of DNA bases is one of the most detrimental The vulnerability of DNA and RNA to the harmful effects of Open Access Article. Published on 22 February 2018. Downloaded 12/18/2020 10:02:34 AM. phenomena occurring during the exposure of nucleic acid UV-radiation is particularly intriguing in the context of the strands to UV-light.1 Cyclobutane pyrimidine dimers (CpDs) origins of life. First of all, it is reasonable to assume that the and (6-4) lesions can be classied among the most frequent highly sophisticated photolyases were absent during the emer- photodimers and are responsible for mutagenic processes gence of the rst oligonucleotides. Furthermore, recent theo- driven by hindered transcription and DNA replication.2 To retical estimates showed that much higher amounts of UV-light protect living organisms from such UV-activated stress, nature were reaching the surface of the Archean Earth, owing to the developed sophisticated enzymes, i.e. photolyases, which lack of oxygen in the atmosphere and higher Sun activity in the – selectively repair CpDs and (6-4) lesions via electron transfer or ultraviolet spectral range.18 22 Interestingly, a recent work by proton-coupled electron transfer processes.3,4 Despite the exis- Bucher et al.23 revealed that specic DNA sequences may tence of repairing machinery, DNA photolesions are a common promote very efficient self-repair driven by photoinduced elec- cause of skin cancer1,5–7 and the mechanistic details of their tron transfer to the CpD lesion containing two dimerized thymine bases (T]T). This remarkable property was discovered for the GAT]T sequence in single and double strands, while no 23 aInstitute of Physics, Polish Academy of Sciences, Al. Lotnikow´ 32/46, PL-02668 repair was reported for the TAT]T and AT]T oligomers. The Warsaw, Poland signicant selectivity of this process implies that UV-light bInstitute of Biophysics of the Czech Academy of Sciences, Kralovopolsk´ a´ 135, 61265 played an important role in the prebiotic selection of the Brno, Czech Republic. E-mail: [email protected] most photochemically stable nucleic acid sequences.24 cRegional Centre of Advanced Technologies and Materials, Department of Physical Many fundamental aspects of the photochemistry of nucleic ´ Chemistry, Faculty of Science, Palacky University, 17. Listopadu 1192/12, 77146 ] Olomouc, Czech Republic acid oligomers (such as GAT T) are still obscure, including the † Electronic supplementary information (ESI) available: Computational details of relative contribution of locally-excited (LE) and delocalized the MD simulations, results of calculations performed for the GA-syn conformer, excitations in DNA strands, the extent of the charge-transfer diabatic couplings and Cartesian coordinates of the stationary points. See DOI: (CT) character of the latter, and the availability of different 10.1039/c8sc00024g This journal is © The Royal Society of Chemistry 2018 Chem. Sci.,2018,9,3131–3140 | 3131 View Article Online Chemical Science Edge Article photorelaxation pathways in such assemblies.25,26 While some Computational methods valuable insights into these processes were provided by time- resolved (TR) spectroscopic (e.g. transient absorption UV) MD simulations techniques,27 it is generally difficult to study the clear origin of The initial geometry of the GAT]T tetranucleotide was the long-lived states due to overlapping absorption features of prepared in the B-DNA helical form. The T]T residue was different nucleobases.25 In addition, selective synthesis of prepared in the same way as the standard AMBER nucleotides.40 specic photodamaged sequences is challenging. On the other The torsion parameters for the cyclobutane ring in the T]T hand, considerable system sizes oen restrict the applications dimer given by the standard atomic types of the parmOL15 force of theoretical simulations to methods such as Time-Dependent eld41–44 were found satisfactory according to our test simula- Density Functional Theory (TDDFT). Even though TDDFT was tion. MD simulations were carried out under the parmOL15 employed with some success in investigations of photochemical force eld41–44 in an octahedral box of explicit water molecules and photophysical properties of oligonucleotides,28,29 the (either SPC/E,45 or OPC model46) with 0.15 M excess KCl,47 using applicability of this methodology to similar problems has been the AMBER soware package.48 Each of these simulations was criticized because of oen signicant overstabilization of CT propagated for 10 ms. The trajectories were clustered by states.26,30,31 While range-separated functionals can alleviate this a custom modied algorithm49 of Rodriguez et al.,50 using the problem to some extent for vertical excitations, nonadiabatic eRMSD metric.51 More details about the MD simulations are molecular dynamics simulations of adenine with TDDFT failed included in the ESI† to this article. to correctly describe the excited-state lifetimes and photo- deactivation pathways.32 QM/MM simulations In the light of the above discussion, characterization of photoinduced processes in oligonucleotides based on highly The geometries of selected (GA-anti and GA-syn) conformers accurate quantum-chemical methods is urgently needed. In obtained from the clustering procedure were further utilized in the QM/MM calculations. These geometries were solvated in Creative Commons Attribution-NonCommercial 3.0 Unported Licence. particular, a much more balanced description of the LE and CT – a spherical droplet of SPC/E explicit water molecules, with states, even outside the Franck Condon region, can be obtained ˚ with the algebraic diagrammatic construction to the second a radius of 25 A. The solvent was then equilibrated for 10 ps, in order [ADC(2)] method.33,34 The ADC(2) method was recently order to obtain a reasonable distribution of the water molecules ff proved to provide reliable results in the investigations of the within the sphere. We considered two di erent QM/MM setups: photochemistry of nucleic acids fragments within the QM/MM QMDNA/MM setup contained the whole tetranucleotide treated framework (referred
Recommended publications
  • Multidimensional Dynamics Involving a Conical Intersection: Wavepacket Calculations Using the MCTDH Method
    Multidimensional dynamics involving a conical intersection: Wavepacket calculations using the MCTDH method G. A. Worth†, H.-D. Meyer∗, and L. S. Cederbaum∗ ∗ Theoretical Chemistry, University of Heidelberg, INF 229, 69120 Heidelberg, Germany † Dept. of Chemistry, King’s College London, The Strand, London, WC2R 2LS, U.K. This article will appear in Conical Intersections, W. Domcke, D. R. Yarkony, and H. K¨oppel, Eds., World Scientific Co., (Singapore, 2004), pp. 573–607. June 1, 2004 1 Introduction The time-evolution of a molecular system through a conical intersection is multi- dimensional in nature, i.e. its description cannot be reduced to a single reaction coordinate. In addition to this, such processes are inherently quantum mechan- ical due to the mixing of electronic states as the intersection is traversed, and a powerful quantum dynamics method is required to treat these, often large, systems accurately. The most intuitive way of obtaining a numerically exact solution of the time- dependent Schr¨odinger equation is to use standard wavepacket propagation meth- ods. These represent the evolving system wavefunction using basis set expansion techniques or grid representations [1–3]. Unfortunately, the standard methods are restricted by an exponential increase in the computational resources required with the number of degrees of freedom treated, usually no more than 4-6 nuclear modes can be included. Systems in which vibronic coupling plays a role are often much larger than this. To treat these systems using the standard methods it is then necessary to turn to reduced dimensionality models in which only a few degrees of freedom are explicitly included, and the remainder either ignored or included in a phenomenological manner.
    [Show full text]
  • European Conference on the Dynamics of Molecular Systems Book of Abstracts
    European conference on the dynamics of molecular systems Book of abstracts Department of Chemistry Jesus College University of Oxford 1 2 Plenary talks 3 PLENARY TALK P1 Imaging studies of inelastic collisions David H. Parker Dept of Molecular and Laser Physics, IMM, Radboud University, Nijmegen, the Netherlands. E‐mail address: [email protected] State-to-state imaging studies of the interaction of small molecules including O2,H2O, OH, CO and ammonia with H2 and with He will be described. These molecules are prominent components of molecular interstellar matter. In order to extract quantitative information on these molecules from telescope spectra, models are employed that depend critically on the rates of rotational energy exchange due to collisions with molecular hydrogen and helium. Collision rates are currently determined by theory from the multidimensional Potential Energy Surface (PES) describing the interaction of H2O and H2 [1] or He [2] and OH with H2 [3], or He. Our velocity map imaging [4] measurements of state-to-state differential and relative integral cross sections of rotational inelastic collisions, also as a function of collision energy, are used to test these PESs. For the H2O/H2, He system experiment is compared with state-of-the-art theoretical calculations by the group of L. Wiesenfeld (Grenoble) [5]. Our studies of astrochemistry relevant small molecules with the collision partners He and H2 at collision energy relevant to that of the interstellar media should place the theoretically determined PESs and the collision rates extracted from the PES on a firmer basis. References 1. P. Valiron M. Wernli, A.
    [Show full text]
  • The Dissipative Photochemical Origin of Life: UVC Abiogenesis of Adenine
    Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 25 January 2021 doi:10.20944/preprints202101.0500.v1 Article The Dissipative Photochemical Origin of Life: UVC Abiogenesis of Adenine Karo Michaelian Department of Nuclear Physics and Application of Radiation, Instituto de Física, Universidad Nacional Autónoma de México, Circuito Interior de la Investigación Científica, Cuidad Universitaria, Cuidad de México, C.P. 04510.; karo@fisica.unam.mx Academic Editor: Karo Michaelian Version January 21, 2021 submitted to Entropy 1 Abstract: I describe the non-equilibrium thermodynamics and the photochemical mechanisms which 2 may have been involved in the dissipative structuring, proliferation and evolution of the fundamental 3 molecules at the origin of life from simpler and more common precursor molecules under the 4 impressed UVC photon flux of the Archean. Dissipative structuring of the fundamental molecules is 5 evidenced by their strong and broad wavelength absorption bands and rapid radiationless dexcitation 6 in this wavelength region. Molecular configurations with large photon dissipative efficacy become 7 dominant through dissipative selection. Proliferation arises from the auto- and cross-catalytic nature 8 of the intermediate products. This inherent non-linearity gives rise to numerous stationary states. 9 Amplification of a molecular concentration fluctuation near a bifurcation allows evolution of the 10 concentration profile towards states of generally greater photon disspative efficacy. An example is 11 given of photochemical dissipative abiogenesis of adenine from the precursors HCN and H2O within 12 a fatty acid vesicle on a hot ocean surface, driven far from equilibrium by the impressed UVC light. 13 The kinetic equations for the photochemical reactions with diffusion are resolved under different 14 environmental conditions and the results analyzed within the framework of Classical Irreversible 15 Thermodynamic theory.
    [Show full text]
  • The Dissipative Photochemical Origin of Life: UVC Abiogenesis of Adenine
    entropy Article The Dissipative Photochemical Origin of Life: UVC Abiogenesis of Adenine Karo Michaelian Department of Nuclear Physics and Applications of Radiation, Instituto de Física, Universidad Nacional Autónoma de México, Circuito Interior de la Investigación Científica, Cuidad Universitaria, Mexico City, C.P. 04510, Mexico; karo@fisica.unam.mx Abstract: The non-equilibrium thermodynamics and the photochemical reaction mechanisms are described which may have been involved in the dissipative structuring, proliferation and complex- ation of the fundamental molecules of life from simpler and more common precursors under the UVC photon flux prevalent at the Earth’s surface at the origin of life. Dissipative structuring of the fundamental molecules is evidenced by their strong and broad wavelength absorption bands in the UVC and rapid radiationless deexcitation. Proliferation arises from the auto- and cross-catalytic nature of the intermediate products. Inherent non-linearity gives rise to numerous stationary states permitting the system to evolve, on amplification of a fluctuation, towards concentration profiles providing generally greater photon dissipation through a thermodynamic selection of dissipative efficacy. An example is given of photochemical dissipative abiogenesis of adenine from the precursor HCN in water solvent within a fatty acid vesicle floating on a hot ocean surface and driven far from equilibrium by the incident UVC light. The kinetic equations for the photochemical reactions with diffusion are resolved under different environmental conditions and the results analyzed within the framework of non-linear Classical Irreversible Thermodynamic theory. Keywords: origin of life; dissipative structuring; prebiotic chemistry; abiogenesis; adenine; organic molecules; non-equilibrium thermodynamics; photochemical reactions Citation: Michaelian, K. The Dissipative Photochemical Origin of MSC: 92-10; 92C05; 92C15; 92C40; 92C45; 80Axx; 82Cxx Life: UVC Abiogenesis of Adenine.
    [Show full text]
  • Hydrogen-Adduction to Open-Shell Graphene Fragments: Spectroscopy, Thermochemistry and Cite This: Chem
    Chemical Science View Article Online EDGE ARTICLE View Journal | View Issue Hydrogen-adduction to open-shell graphene fragments: spectroscopy, thermochemistry and Cite this: Chem. Sci.,2017,8,1186 astrochemistry Gerard D. O'Connor,a Bun Chan,bc Julian A. Sanelli,d Katie M. Cergol,b Viktoras Dryza,d Richard J. Payne,b Evan J. Bieske,d Leo Radomb and Timothy W. Schmidt*a We apply a combination of state-of-the-art experimental and quantum-chemical methods to elucidate the electronic and chemical energetics of hydrogen adduction to a model open-shell graphene fragment. The lowest-energy adduct, 1H-phenalene, is determined to have a bond dissociation energy of 258.1 kJ molÀ1, while other isomers exhibit reduced or in some cases negative bond dissociation energies, the metastable species being bound by the emergence of a conical intersection along the high-symmetry dissociation coordinate. The gas-phase excitation spectrum of 1H-phenalene and its radical cation are recorded Received 24th August 2016 using laser spectroscopy coupled to mass-spectrometry. Several electronically excited states of both Accepted 24th September 2016 Creative Commons Attribution 3.0 Unported Licence. species are observed, allowing the determination of the excited-state bond dissociation energy. The DOI: 10.1039/c6sc03787a ionization energy of 1H-phenalene is determined to be 7.449(17) eV, consistent with high-level W1X-2 www.rsc.org/chemicalscience calculations. 1 Introduction Using density functional theory it has been calculated that semihydrogenated graphene (graphone) becomes a ferromag- Graphene is a material consisting of a single layer of carbon netic semiconductor with a small indirect gap.7 It has also been 1,2 atoms bonded in a hexagon lattice.
    [Show full text]
  • Femtosecond Intersystem Crossing in the DNA Nucleobase Cytosine
    Femtosecond Intersystem Crossing in the DNA Nucleobase Cytosine Martin Richter,y Philipp Marquetand,z Jes´usGonz´alez-Vazquez,{ Ignacio Sola,{ and Leticia Gonz´alez∗,z yInstitute of Physical Chemistry, Friedrich Schiller University Jena, Helmholtzweg 4, 07743 Jena, Germany zInstitute of Theoretical Chemistry, University of Vienna, W¨ahringerStr. 17, 1090 Vienna, Austria {Departamento de Qu´ımica F´ısica I, Universidad Complutense, 28040 Madrid, Spain E-mail: [email protected] Keywords: DNA photostability, excited Graphical TOC Entry state dynamics, intersystem crossing, spin-orbit coupling, conical intersection Abstract Ab initio molecular dynamics including non- adiabatic and spin-orbit couplings on equal footing is used to unravel the deactivation of cytosine after UV light absorption. Intersys- tem crossing (ISC) is found to compete directly with internal conversion in tens of femtosec- onds, thus making cytosine the organic com- pound with the fastest triplet population cal- culated so far. It is found that close degen- eracy between singlet and triplet states can more than compensate for very small spin-orbit couplings, leading to efficient ISC. The fem- tosecond nature of the intersystem crossing pro- cess highlights its importance in photochem- arXiv:2103.13070v1 [physics.chem-ph] 24 Mar 2021 istry and challenges the conventional view that large singlet-triplet couplings are required for an efficient population flow into triplet states. These findings are important to understand DNA photostability and the photochemistry and dynamics of organic molecules in general. 1 The interaction of DNA and RNA with ra- goes beyond what can be learned from quan- diation, from mobile-phone emissions1 to UV tum chemical calculations alone.
    [Show full text]
  • Metal-Organic Charge Transfer Can Produce Biradical States and Is Mediated by Conical Intersections
    Metal-organic charge transfer can produce biradical INAUGURAL ARTICLE states and is mediated by conical intersections Oksana Tishchenkoa,1, Ruifang Lia,b, and Donald G. Truhlara,1 aDepartment of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, MN 55455-0431; and bDepartment of Chemistry, Nankai University, Tianjin 300071, People’s Republic of China This article is part of the special series of Inaugural Articles by members of the National Academy of Sciences elected in 2009. Edited by Nicholas J. Turro, Columbia University, New York, NY, and approved September 20, 2010 (received for review July 14, 2010) The present paper illustrates key features of charge transfer be- theory revealed the presence of two types of equilibrium struc- tween calcium atoms and prototype conjugated hydrocarbons tures: one at a distance R between Ca and the closest carbon (ethylene, benzene, and coronene) as elucidated by electronic of about 3.8 Å with slightly negative charge on the calcium atom, structure calculations. One- and two-electron charge transfer is and another state with a smaller R and with considerable shift of controlled by two sequential conical intersections. The two lowest electrons from Ca to coronene; the latter is metastable. To get a electronic states that undergo a conical intersection have closed- better understanding of the structure of potential energy surfaces shell and open-shell dominant configurations correlating with that govern the charge transfer process and of the nature of the 4s2 and 4s13d1 states of Ca, respectively. Unlike the neutral- electronic states involved, we performed a detailed electronic ionic state crossing in, for example, hydrogen halides or alkali structure study using multiconfigurational wave functions for halides, the path from separated reactants to the conical intersec- the interaction of Ca with the smaller prototype molecules ben- tion region is uphill and the charge-transferred state is a biradical.
    [Show full text]
  • Photoprotective Properties of Eumelanin: Computational Insights Into the Photophysics of a Catechol:Quinone Heterodimer Model System
    Article Photoprotective Properties of Eumelanin: Computational Insights into the Photophysics of a Catechol:Quinone Heterodimer Model System Victoria C. Frederick 1, Thomas A. Ashy 1, Barbara Marchetti 1, Michael N. R. Ashfold 2 and Tolga N. V. Karsili 1,* 1 Department of Chemistry, University of Louisiana at Lafayette, Lafayette, LA 70503, USA; [email protected] (V.C.F.); [email protected] (T.A.A.); [email protected] (B.M.) 2 School of Chemistry, University of Bristol, Bristol BS8 1TS, UK; [email protected] * Correspondence: [email protected] Abstract: Melanins are skin-centered molecular structures that block harmful UV radiation from the sun and help protect chromosomal DNA from UV damage. Understanding the photodynamics of the chromophores that make up eumelanin is therefore paramount. This manuscript presents a multi-reference computational study of the mechanisms responsible for the experimentally observed photostability of a melanin-relevant model heterodimer comprising a catechol (C)–benzoquinone (Q) pair. The present results validate a recently proposed photoinduced intermolecular transfer of an H atom from an OH moiety of C to a carbonyl-oxygen atom of the Q. Photoexcitation of the ground state C:Q heterodimer (which has a π-stacked “sandwich” structure) results in population of a locally excited ππ* state (on Q), which develops increasing charge-transfer (biradical) character as it evolves Citation: Frederick, V.C.; Ashy, T.A.; to a “hinged” minimum energy geometry and drives proton transfer (i.e., net H atom transfer) from Marchetti, B.; Ashfold, M.N.R.; C to Q. The study provides further insights into excited state decay mechanisms that could contribute Karsili, T.N.V.
    [Show full text]
  • Evidence for Conical Intersection Dynamics Mediating Ultrafast Singlet Exciton Fission
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Apollo Evidence for Conical Intersection Dynamics Mediating Ultrafast Singlet Exciton Fission Andrew J. Musser1*, Matz Liebel2*, Christoph Schnedermann2, Torsten Wende2, Tom B. 1 1 2 Kehoe , Akshay Rao , Philipp Kukura 1 Cavendish Laboratory, University of Cambridge J J Thomson Avenue, Cambridge, UK 2 Physical and Theoretical Chemistry Laboratory, University of Oxford South Parks Road, Oxford, UK Email: [email protected], [email protected] *These authors contributed equally to this work. Introductory Paragraph Singlet exciton fission is the process in organic semiconductors through which a spin-singlet exciton converts into a pair of spin-triplet excitons residing on different chromophores, entangled in an overall spin-zero state. For some systems, singlet fission has been shown to occur on the 100 fs timescale and with a 200% yield, but the mechanism of this process remains uncertain. Here we study a model singlet fission system, TIPS-pentacene, using ultrafast vibronic spectroscopy. We observe that vibrational coherence in the initially photogenerated singlet state is transferred to the triplet state and show that this behaviour is effectively identical to that observed in ultrafast internal conversion for polyenes in solution. This similarity in vibronic dynamics suggest that both multi-molecular singlet fission and single-molecular internal conversion are mediated by the same underlying relaxation processes, based on strong coupling between nuclear and electronic degrees of freedom. In its most efficient form this leads to a conical intersection between the coupled electronic states. 1 Main Text Singlet exciton fission has long commanded interest as an exceptionally fast channel to generate triplet excitons in organic materials.
    [Show full text]
  • Ultrafast Deactivation of an Excited Cytosine-Guanine Base Pair in DNA Gerrit Groenhof,† Lars V
    Published on Web 05/08/2007 Ultrafast Deactivation of an Excited Cytosine-Guanine Base Pair in DNA Gerrit Groenhof,† Lars V. Scha¨fer,† Martial Boggio-Pasqua,‡ Maik Goette,† Helmut Grubmu¨ller,† and Michael A. Robb*,‡ Contribution from the Department of Theoretical and Computational Biophysics, Max-Planck-Institute for Biophysical Chemistry, Am Fassberg 11, D-37077 Go¨ttingen, Germany, and Department of Chemistry, Imperial College, London SW 2AZ, U.K. Received December 31, 2006; E-mail: [email protected] Abstract: Multiconfigurational ab initio calculations and QM/MM molecular dynamics simulations of a photoexcited cytosine-guanine base pair in both gas phase and embedded in the DNA provide detailed structural and dynamical insights into the ultrafast radiationless deactivation mechanism. Photon absorption promotes transfer of a proton from the guanine to the cytosine. This proton transfer is followed by an efficient radiationless decay of the excited state via an extended conical intersection seam. The optimization of the conical intersection revealed that it has an unusual topology, in that there is only one degeneracy- lifting coordinate. This is the central mechanistic feature for the decay both in vacuo and in the DNA. Radiationless decay occurs along an extended hyperline nearly parallel to the proton-transfer coordinate, indicating the proton transfer itself is not directly responsible for the deactivation. The seam is displaced from the minimum energy proton-transfer path along a skeletal deformation of the bases. Decay can thus occur anywhere along the single proton-transfer coordinate, accounting for the remarkably short excited- state lifetime of the Watson-Crick base pair. In vacuo, decay occurs after a complete proton transfer, whereas in DNA, decay can also occur much earlier.
    [Show full text]
  • Ultrafast 25-Fs Relaxation in Highly Excited States of Methyl Azide Mediated by Strong Nonadiabatic Coupling
    Ultrafast 25-fs relaxation in highly excited states of PNAS PLUS methyl azide mediated by strong nonadiabatic coupling William K. Petersa,b,1, David E. Coucha,b, Benoit Mignoletc, Xuetao Shid, Quynh L. Nguyena,b, Ryan C. Fortenberrye, H. Bernhard Schlegeld, Françoise Remaclec, Henry C. Kapteyna,b, Margaret M. Murnanea,b,1, and Wen Lid aJILA, University of Colorado, Boulder, CO 80309; bDepartment of Physics, University of Colorado, Boulder, CO 80309; cTheoretical Physical Chemistry, Unité de Recherche Molecular Systems (UR MOLSYS) University of Liège, B4000 Liège, Belgium; dDepartment of Chemistry, Wayne State University, Detroit, MI 48202; and eDepartment of Chemistry and Biochemistry, Georgia Southern University, Statesboro, GA 30460 Contributed by Margaret M. Murnane, October 14, 2017 (sent for review July 16, 2017; reviewed by Elliot Bernstein and Arthur Suits) Highly excited electronic states are challenging to explore exper- When molecules are excited by vacuum-UV (VUV) light in imentally and theoretically—due to the large density of states and the 7- to 10-eV range, just below the ionization potential, one or the fact that small structural changes lead to large changes in several excited states can be accessed depending on the band- electronic character with associated strong nonadiabatic dynamics. width of the pulse. In this energy range, there is typically a dense They can play a key role in astrophysical and ionospheric chemis- manifold of valence and Rydberg states that are often strongly try, as well as the detonation chemistry of high-energy density coupled, each with an associated vibrational structure (9). As a materials. Here, we implement ultrafast vacuum-UV (VUV)-driven result, VUV molecular absorption spectra typically contain electron–ion coincidence imaging spectroscopy to directly probe quasicontinuous but structured bands with overlapping sharp the reaction pathways of highly excited states of energetic mole- peaks, indicating the wide variety of timescales and states in- — cules in this case, methyl azide.
    [Show full text]
  • Ultrafast Radiationless Decay Mechanisms Through Conical Intersections in Cytosine. a New Semi-Planar Conical Intersection
    Indian Journal of Chemistry Vol. 53A, February 2014, pp. 143-151 Ultrafast radiationless decay mechanisms through conical intersections in cytosine. A new semi-planar conical intersection Saadullah G Aziza, Walid I Hassanb, Shaaban K Alroubya, c, Abdulrahman Alyoubia & Rifaat Hilala, b, * aChemistry Department, Faculty of Science, King Abdulaziz University, Jeddah, SA bChemistry Department, Faculty of Science Cairo University, Giza, Egypt cChemistry Department, Faculty of Science Beni-Suief University, Beni-Suief, Egypt Email: [email protected] Received 26 November 2013; revised and accepted 27 January 2014 Geometry, energetics and dipole moment of all possible conformers of cytosine in the ground state are calculated using density functional theory B3LYP method and the 6-311++G(3df,3pd) basis set. The most stable conformer is keto-amino conformer. The amino group hydrogen atoms are slightly out of plane by about 6.3° and 9.9° each. Ultrafast radiationless decay mechanism has been theoretically investigated using Complete Active Space Multiconfiguration SCF calculations. Effective pathways of ultrafast radiationless transitions from the optically allowed π π* state to the ground state S0 of cytosine are explored. The nπ*, nσ*, and the ππ* states have been taken into account as states involved in the radiationless process. Optimized geometry and conical intersections have been searched in the full dimensional space for the vibrational degrees of freedom. Three competing direct decay mechanisms through three possible conical intersections have been found to exist. The first pathway is through the bending of molecule in a sofa-like structure leading to conical intersection with ground state at 4.23 eV. The second pathway occurs through a twisted structure that has hydrogen twisted and the cytosine ring slightly deformed leading to conical intersection at 4.08 eV.
    [Show full text]