Benign decay vs. photolysis in the photophysics and photochemistry of 5-bromouracil. A computational study

Excited States of Bromopyrimidines Probed by VUV Photoabsorption Spectroscopy and Theoretical Calculations

We report absolute photoabsorption cross sections for gas-phase 2- and 5-bromopyrimidine in the 3.7-10.8 eV energy range, in a joint theoretical and experimental study. The measurements were carried out using high-resolution vacuum ultraviolet synchrotron radiation, with quantum chemical calculations performed through the nuclear ensemble approach in combination with time-dependent density functional theory, along with additional Franck-Condon Herzberg-Teller calculations for the first absorption band (3.7-4.6 eV). The cross sections of both bromopyrimidines are very similar below 7.3 eV, deviating more substantially from each other at higher energies. In the 7.3-9.0 eV range where the maximum cross-section is found, a single and broad band is observed for 5-bromopyrimidine, while more discernible features appear in the case of 2-bromopyrimidine. Several π* ← π transitions account for the most intense bands, while weaker ones are assigned to transitions involving the nitrogen and bromine lone pairs, the antibonding σ*Br orbital, and the lower-lying Rydberg states. A detailed comparison with the available photo-absorption data of bromobenzene is also reported. We have found significant differences regarding the main absorption band, which is more peaked in bromobenzene, becoming broader and shifting to higher energies in both bromopyrimidines. In addition, there is a significant suppression of vibrational structures and of Rydberg states in the pair of isomers, most noticeably for 2-bromopyrimidine.

Photoinduced electron transfer in 5-bromouracil labeled DNA. A contrathermodynamic mechanism revisited by electron transfer theories

The understanding of the 5-bromouracil (BrU) based photosensitization mechanism of DNA damage is of large interest due to the potential applications in photodynamic therapy. Photoinduced electron transfer (ET) in BrU labeled duplexes comprising the 5'-GBrU or 5'-ABrU sequence showed that a much lower reactivity was found for the 5'-GBrU pattern. Since the ionization potential of G is lower than that of A, this sequence selectivity has been dubbed a contrathermodynamic one. In the current work, we employ the Marcus and Marcus-Levich-Jortner theory of ET in order to shed light on the observed effect. By using a combination of Density Functional Theory (DFT) and solvation continuum models, we calculated the electronic couplings, reorganization energies, and thermodynamic stimuli for electron transfer which enabled the rates of forward and back ET to be estimated for the two considered sequences. The calculated rates show that the photoreaction could not be efficient if the ET process proceeded within the considered dimers. Only after introducing additional adenines between G and BrU, which accelerates the forward and slows down the back ET, is a significant amount of photodamage expected.

Vacuum-UV induced DNA strand breaks - influence of the radiosensitizers 5-bromouracil and 8-bromoadenine

Radiation therapy is a basic part of cancer treatment. To increase the DNA damage in carcinogenic cells and preserve healthy tissue at the same time, radiosensitizing molecules such as halogenated nucleobase analogs can be incorporated into the DNA during the cell reproduction cycle. In the present study 8.44 eV photon irradiation induced single strand breaks (SSB) in DNA sequences modified with the radiosensitizer 5-bromouracil (5BrU) and 8-bromoadenine (8BrA) are investigated. 5BrU was incorporated in the 13mer oligonucleotide flanked by different nucleobases. It was demonstrated that the highest SSB cross sections were reached, when cytosine and thymine were adjacent to 5BrU, whereas guanine as a neighboring nucleobase decreases the activity of 5BrU indicating that competing reaction mechanisms are active. This was further investigated with respect to the distance of guanine to 5BrU separated by an increasing number of adenine nucleotides. It was observed that the SSB cross sections were decreasing with an increasing number of adenine spacers between guanine and 5BrU until the SSB cross sections almost reached the level of a non-modified DNA sequence, which demonstrates the high sequence dependence of the sensitizing effect of 5BrU. 8BrA was incorporated in a 13mer oligonucleotide as well and the strand breaks were quantified upon 8.44 eV photon irradiation in direct comparison to a non-modified DNA sequence of the same composition. No clear enhancement of the SSB yield of the modified in comparison to the non-modified DNA sequence could be observed. Additionally, secondary electrons with a maximum energy of 3.6 eV were generated when using Si as a substrate giving rise to further DNA damage. A clear enhancement in the SSB yield can be ascertained, but to the same degree for both the non-modified DNA sequence and the DNA sequence modified with 8BrA.

Discovery of conical intersection mediated photochemistry with growing string methods

Conical intersections (CIs) are important features of photochemistry that determine yields and selectivity. Traditional CI optimizers require significant human effort and chemical intuition, which typically restricts searching to only a small region of the CI space. Herein, a systematic approach utilizing the growing string method is introduced to locate multiple CIs. Unintuitive MECI are found using driving coordinates that can be generated using a combinatorial search, and subsequent optimization allows reaction pathways, transition states, products, and seam-space pathways to be located. These capabilities are demonstrated by application to two prototypical photoisomerization reactions and the dimerization of butadiene. In total, many reaction pathways were uncovered, including the elusive stilbene hula-twist mechanism, and a previously unidentified product in butadiene dimerization. Overall, these results suggest that growing string methods provide a predictive strategy for exploring photochemistry.

UV-induced electron transfer between triethylamine and 5-bromo-2'-deoxyuridine. A puzzle concerning the photochemical debromination of labeled DNA

5-Bromo-2'-deoxyuridine (BrdU) photosensitizes DNA to strand break formation. However, this type of photodamage is completely quenched by the presence of triethylamine (TEA) which originates from RP-HPLC purification commonly employed by oligonucleotide providers. While the presence of TEA in oligonucleotide samples does not interfere with PCR or other molecular biology applications, the mechanism of photochemical reaction proceeding in the labeled DNA is dramatically changed due to the photoinduced electron transfer (PET) between the photoexcited BrdU and the ground state TEA. For the first time, we demonstrated that the latter process produces 2'-deoxyuridne2'-deoxyuridine (debromination) in the labeled DNA instead of the expected strand break. PET between TEA and BrdU was additionally confirmed by the UV irradiations of aqueous solutions containing both species. Indeed, the efficient formation of 2'-deoxyuridine was observed in the studied photolytes. Moreover, we showed the formation of an additional product in these binary mixtures, i.e. imidazole derivative, that is not formed in DNA and was reported in the literature in the context of dark rather than photochemical processes. Using mass spectrometry we demonstrated that the amount of TEA impurity in the commercial samples of oligos exceeds up to 3 orders of magnitude that of the purchased DNA.

Insights into the deactivation of 5-bromouracil after ultraviolet excitation

5-Bromouracil is a nucleobase analogue that can replace thymine in DNA strands and acts as a strong radiosensitizer, with potential applications in molecular biology and cancer therapy. Here, the deactivation of 5-bromouracil after ultraviolet irradiation is investigated in the singlet and triplet manifold by accurate quantum chemistry calculations and non-adiabatic dynamics simulations. It is found that, after irradiation to the bright ππ* state, three main relaxation pathways are, in principle, possible: relaxation back to the ground state, intersystem crossing (ISC) and C-Br photodissociation. Based on accurate MS-CASPT2 optimizations, we propose that ground-state relaxation should be the predominant deactivation pathway in the gas phase. We then employ different electronic structure methods to assess their suitability to carry out excited-state dynamics simulations. MRCIS (multi-reference configuration interaction including single excitations) was used in surface hopping simulations to compute the ultrafast ISC dynamics, which mostly involves the 1nOπ* and 3ππ* states.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.

Early events in the photochemistry of 5-diazo Meldrum's acid: formation of a product manifold in C–N bound and pre-dissociated intersection seam regions

Photolysis of 5-diazo Meldrum's acid yields three different primary products in a sub-ps scale. High-level computations show that this is due to the presence of an extended seam of intersection associated to different reaction paths.

Photochemistry and photophysics at extended seams of conical intersection

The role of extended seams of conical intersection in excited-state mechanisms is reviewed. Seams are crossings of the potential energy surface in many dimensions where the decay from the excited to the ground state can occur, and the extended seam is composed of different segments lying along a reaction coordinate. Every segment is associated with a different primary photoproduct, which gives rise to competing pathways. This idea is first illustrated for fulvene and ethylene, and then it is used to explain more complex cases such as the dependence of the isomerisation of retinal chromophore isomers on the protein environment, the dependence of the efficiency of the azobenzene photochemical switch on the wavelength of irradiation and the direction of the isomerisation, and the coexistence of different mechanisms in the photo-induced Wolff rearrangement of diazonaphthoquinone. The role of extended seams in the photophysics of the DNA nucleobases and the relationship between two-state seams and three-state crossings is also discussed. As an outlook, the design of optical control strategies based on the passage of the excited molecule through the seam is considered, and it is shown how the excited-state lifetime of fulvene can be modulated by shaping the energy of the seam.

Photophysics of fulvene under the non-resonant stark effect. Shaping the conical intersection seam

We introduce a mechanistic strategy to control the excited state lifetime of fulvene based on shaping the topography of an extended seam of intersection with the non-resonant dynamic Stark effect. Fulvene has a very short excited state lifetime due to an energetically accessible seam of intersection which lies along the methylene torsion coordinate, and the initial decay occurs at the seam segment around the planar conical intersection structure. We have followed a three-step approach to simulate the control. First, we have calculated the effect of a non-resonant electric field on the potential energy surface at the ab initio level, including the field in a self-consistent way. The relative energy of the planar segment of the seam is increased by the non-resonant field. In the second step we simulate the control carrying out MCTDH quantum dynamics propagations under a static non-resonant field to derive the main control mechanisms. At moderately intense fields (epsilon < or = 0.03 a.u.) the decay is faster as compared to the field free case because the vibrational overlap between the excited and ground state vibrational functions is increased. However, at more intense fields (epsilon = 0.04 a.u.) the planar conical intersection is energetically inaccessible and the decay occurs at a slower time scale, at the segment of the seam with more twisted geometries. In the third step, the control over the dynamics is exerted with a non-resonant dynamic field. The acceleration of the decay due to the improved vibrational overlap does not occur, but the decay can be made slower with a dynamic field of 0.08 a.u. The results show the viability of our approach to control the photophysics shaping the topology of the conical intersection seam, and they prove that the extended nature of the seam is crucial for simulating and understanding the control.

Wave Packet Dynamics at an Extended Seam of Conical Intersection: Mechanism of the Light-Induced Wolff Rearrangement

Quantum dynamics calculations on a model surface based on CASPT2//CASSCF calculations are carried out to probe the traversal of a wave packet through an extended seam of conical intersection during the light-induced Wolff rearrangement of diazonaphtoquinone. The reaction is applied in the fabrication of integrated circuits. It consists of nitrogen elimination and ring rearrangement to yield a ketene. After excitation, the wave packet relaxes and reaches the extended seam. A fraction of the wave packet decays to the ground state at a region of the seam connected to a carbene intermediate, while the remaining part decays at a region leading to the ketene. The passage of the wave packet through the extended seam explains the competition between concerted ketene formation and a stepwise mechanism involving a carbene. The two primary photoproducts are formed in the first 100 fs of the simulation, in agreement with recent ultrafast spectroscopy measurements.

Quantum dynamics study of fulvene double bond photoisomerization: the role of intramolecular vibrational energy redistribution and excitation energy

The double bond photoisomerization of fulvene has been studied with quantum dynamics calculations using the multi-configuration time-dependent Hartree method. Fulvene is a test case to develop optical control strategies based on the knowledge of the excited state decay mechanism. The decay takes place on a time scale of several hundred femtoseconds, and the potential energy surface is centered around a conical intersection seam between the ground and excited state. The competition between unreactive decay and photoisomerization depends on the region of the seam accessed during the decay. The dynamics are carried out on a four-dimensional model surface, parametrized from complete active space self-consistent field calculations, that captures the main features of the seam (energy and locus of the seam and associated branching space vectors). Wave packet propagations initiated by single laser pulses of 5-25 fs duration and 1.85-4 eV excitation energy show the principal characteristics of the first 150 fs of the photodynamics. Initially, the excitation energy is transferred to a bond stretching mode that leads the wave packet to the seam, inducing the regeneration of the reactant. The photoisomerization starts after the vibrational energy has flowed from the bond stretching to the torsional mode. In our propagations, intramolecular energy redistribution (IVR) is accelerated for higher excess energies along the bond stretch mode. Thus, the competition between unreactive decay and isomerization depends on the rate of IVR between the bond stretch and torsion coordinates, which in turn depends on the excitation energy. These results set the ground for the development of future optical control strategies.

Formation mechanism and structure of a guanine-uracil DNA intrastrand cross-link

The formation and structure of the 5'-G[8-5]U-3' intrastrand cross-link are studied using density functional theory and molecular dynamics simulations due to the potential role of this lesion in the activity of 5-halouracils in antitumor therapies. Upon UV irradiation of 5-halouracil-containing DNA, a guanine radical cation reacts with the uracil radical to form the cross-link, which involves phosphorescence or an intersystem crossing and a rate-determining step of bond formation. Following ionizing radiation, guanine and the uracil radical react, with a rate-limiting step involving hydrogen atom removal. Although cross-link formation from UV radiation is favored, comparison of calculated reaction thermokinetics with that for related experimentally observed purine-pyrimidine cross-links suggests this lesion is also likely to form from ionizing radiation. For the first time, the structure of 5'-G[8-5]U-3' within DNA is identified by molecular dynamics simulations. Furthermore, three conformations of cross-linked DNA are revealed, which differ in the configuration of the complementary bases. Distortions, such as unwinding, are localized to the cross-linked dinucleotide and complementary nucleotides, with minimal changes to the flanking bases. Global changes to the helix, such as bending and groove alterations, parallel cisplatin-induced distortions, which indicate 5'-G[8-5]U-3', may contribute to the cytotoxicity of halouracils in tumor cell DNA using similar mechanisms.

Developing a computational model that accurately reproduces the structural features of a dinucleoside monophosphate unit within B-DNA

The ability of a dinucleoside monophosphate to mimic the conformation of B-DNA was tested using a combination of different phosphate models (anionic, neutral, counterion), environments (gas, water), and density functionals (B3LYP, MPWB1K, M06-2X) with the 6-31G(d,p) basis set. Three sequences (5'-GX(Py)-3', where X(Py) = T, U or (Br)U) were considered, which vary in the (natural or modified) 3' pyrimidine nucleobase (X(Py)). These bases were selected due to their presence in natural DNA, structural similarity to T and/or applications in radical-initiated anti-tumour therapies. The accuracy of each of the 54 model, method and sequence combinations was judged based on the ability to reproduce key structural features of natural B-DNA, including the stacked base-base orientation and important backbone torsion angles. B3LYP yields distorted or tilted relative base-base orientations, while many computational challenges were encountered for MPWB1K. Despite wide use in computational studies of DNA, the neutral (protonated) phosphate model could not consistently predict a stacked arrangement for all sequences. In contrast, stacked base-base arrangements were obtained for all sequences with M06-2X in conjunction with both the anionic and (sodium) counterion phosphate models. However, comparison of calculated and experimental backbone conformations reveals the charge-neutralized counterion phosphate model best mimics B-DNA. Structures optimized with implicit solvent (water) are comparable to gas-phase structures, which suggests similar results should be obtained in an environment of intermediate polarity. We recommend the use of either gas-phase or water M06-2X optimizations with the counterion phosphate model to study the structure and/or reactivity of natural or modified DNA with a dinucleoside monophosphate.

Unexpected photoproduct generated via the acetone-sensitized photolysis of 5-bromo-2'-deoxyuridine in a water/isopropanol solution: experimental and computational studies

The acetone-sensitized photolysis of 5-bromo-2'-deoxyuridine (5-BrdU) in a water/isopropanol solution with 300 nm photons leads to the formation of 2'-deoxyuridine (dU) and a comparable amount of another photoproduct that has not been reported in the literature so far. The negative and positive mass spectra recorded for this species indicate that they originate from the molecular mass of 286 Da, which corresponds to an adduct of 2'-deoxyuridine and 2-propanol. Quantum chemical calculations carried out at the DFT and TDDFT levels reveal both the structure and the UV spectrum of that adduct. The latter computational characteristic matches well the experimental UV spectrum of the new photoproduct. Our findings indicate that the acetone-sensitized photolysis of 5-BrdU is more complicated than has hitherto been assumed. Nevertheless, since electron transfer is one of the pathways responsible for 5-BrdU decay, acetone-sensitized photolysis of the halogen derivatives of nucleobases could be a convenient tool for studying their radiosensitivity in aqueous solutions.