A triplet mechanism for the formation of cyclobutane pyrimidine dimers in UV-irradiated DNA

A Theoretical Exploration of the Photoinduced Breaking Mechanism of the Glycosidic Bond in Thymine Nucleotide

DNA glycosidic bond cleavage may induce cancer under the ultraviolet (UV) effect. Yet, the mechanism of glycosidic bond cleavage remains unclear and requires more detailed clarification. Herein, quantum chemical studies on its photoinduced mechanism are performed using a 5'-thymidine monophosphate (5'-dTMPH) model. In this study, four possible paths were examined to study the glycosidic bond cleavage. The results showed that, upon excitation, the electronic transition from the π bonding to π antibonding orbitals of the thymine ring leads to the damage of the thymine ring. Afterwards, the glycosidic bond is cleaved. At first, the doublet ground state (GS) path of glycosidic bond cleavage widely studied by other groups is caused by free electron generated by photoirradiation, with a kinetically feasible energy barrier of ~23 kcal/mol. Additionally, then, the other three paths were proposed that also might cause the glycosidic bond cleavage. The first one is the doublet excited state (ES) path, triggered by free electron along with UV excitation, which can result in a very-high-energy barrier ~49 kcal/mol that is kinetically unfavorable. The second one is the singlet ES path, induced by direct UV excitation, which assumes DNA is directly excited by UV light, which features a very low-energy barrier ~16 kcal/mol that is favored in kinetics. The third one is the triplet ES path, from the singlet state via intersystem crossing (ISC), which refers to a feasible ~27 kcal/mol energy barrier. This study emphasizes the pivotal role of the DNA glycosidic bond cleavage by our proposed direct UV excitation (especially singlet ES path) in addition to the authorized indirect free-electron-induced path, which should provide essential insights to future mechanistic comprehension and novel anti-cancer drug design.

Triplet Excimer Formation in a DNA Duplex with Silver Ion-Mediated Base Pairs

The dynamics of excited electronic states in self-assembled structures formed between silver(I) ions and cytosine-containing DNA strands or monomeric cytosine derivatives were investigated by time-resolved infrared (TRIR) spectroscopy and quantum mechanical calculations. The steady-state and time-resolved spectra depend sensitively on the underlying structures, which change with pH and the nucleobase and silver ion concentrations. At pH ∼ 4 and low dC20 strand concentration, an intramolecularly folded i-motif is observed, in which protons, and not silver ions, mediate C-C base pairing. However, at the higher strand concentrations used in the TRIR measurements, dC20 strands associate pairwise to yield duplex structures containing C-Ag+-C base pairs with a high degree of propeller twisting. UV excitation of the silver ion-mediated duplex produces a long-lived excited state, which we assign to a triplet excimer state localized on a pair of stacked cytosines. The computational results indicate that the propeller-twisted motifs induced by metal-ion binding are responsible for the enhanced intersystem crossing that populates the triplet state and not a generic heavy atom effect. Although triplet excimer states have been discussed frequently as intermediates in the formation of cyclobutane pyrimidine dimers, we find neither computational nor experimental evidence for cytosine-cytosine photoproduct formation in the systems studied. These findings provide a rare demonstration of a long-lived triplet excited state that is formed in a significant yield in a DNA duplex, demonstrating that supramolecular structural changes induced by metal ion binding profoundly affect DNA photophysics.

Significance of Singlet Oxygen Molecule in Pathologies

Reactive oxygen species, including singlet oxygen, play an important role in the onset and progression of disease, as well as in aging. Singlet oxygen can be formed non-enzymatically by chemical, photochemical, and electron transfer reactions, or as a byproduct of endogenous enzymatic reactions in phagocytosis during inflammation. The imbalance of antioxidant enzymes and antioxidant networks with the generation of singlet oxygen increases oxidative stress, resulting in the undesirable oxidation and modification of biomolecules, such as proteins, DNA, and lipids. This review describes the molecular mechanisms of singlet oxygen production in vivo and methods for the evaluation of damage induced by singlet oxygen. The involvement of singlet oxygen in the pathogenesis of skin and eye diseases is also discussed from the biomolecular perspective. We also present our findings on lipid oxidation products derived from singlet oxygen-mediated oxidation in glaucoma, early diabetes patients, and a mouse model of bronchial asthma. Even in these diseases, oxidation products due to singlet oxygen have not been measured clinically. This review discusses their potential as biomarkers for diagnosis. Recent developments in singlet oxygen scavengers such as carotenoids, which can be utilized to prevent the onset and progression of disease, are also described.

Role of charge transfer states into the formation of cyclobutane pyrimidine dimers in DNA

Photoinduced charge transfer between neighboring bases plays an important role in DNA. One of its important effects is shown in its ability to affect the photochemical yields of the formation of cyclobutane pyrimidine dimer (CPD) products between adjacent pyrimidine bases. In this work we examine how the energies of charge transfer states depend on the sequences of oligonucleotides using a hybrid quantum and molecular mechanics (QM/MM) methodology combined with the algebraic diagrammatic construction through a second order electronic structure method for excited states. Specifically, we examine 10 sequences with guanine being on the 5' or 3' position of two pyrimidine bases. The results show that the energies of charge transfer states are affected by the nature of the donor acceptor pair, by the distance between them, and by other electrostatic effects created by the surrounding environment.

Triplet-Induced Lesion Formation at CpT and TpC Sites in DNA

UV irradiation induces DNA lesions particularly at dipyrimidine sites. Using time-resolved UV pump (250 nm) and mid-IR probe spectroscopy the triplet pathway of cyclobutane pyrimidine dimer (CPD) formation within TpC and CpT sequences was studied. The triplet state is initially localized at the thymine base but decays with 30 ns under formation of a biradical state extending over both bases of the dipyrimidine. Subsequently this state either decays back to the electronic ground state on the 100 ns time scale or forms a cyclobutane pyrimidine dimer lesion (CPD). Stationary IR spectroscopy and triplet sensitization via 2'-methoxyacetophenone (2-M) in the UVA range shows that the lesions are formed with an efficiency of approximately 1.5 %. Deamination converts the cytosine moiety of the CPD lesions on the time scale of 10 hours into uracil which gives CPD(UpT) and CPD(TpU) lesions in which the coding potential of the initial cytosine base is vanished.

UV photostability of three 2-aminoazoles with key roles in prebiotic chemistry on the early earth

Three related molecules in the 2-aminoazole family are potentially important for prebiotic chemistry: 2-aminooxazole, 2-aminoimidazole, and 2-aminothiazole, which can provide critical functions as an intermediate in nucleotide synthesis, a nucleotide activating agent, and a selective agent, respectively. Here, we examine the wavelength-dependent photodegradation of these three molecules under mid-range UV light (210-290 nm). We then assess the implications of the observed degradation rates for the proposed prebiotic roles of these compounds. We find that all three 2-aminoazoles degrade under UV light, with half lives ranging from ≈7-100 hours under a solar-like spectrum. 2-Aminooxazole is the least photostable, while 2-aminoimidazole is the most photostable. The relative photostabilities are consistent with the order in which these molecules would be used prebiotically: AO is used first to build nucleotides and AI is used last to activate them.

Mechanistic insight into photocrosslinking reaction between triplet state 4-thiopyrimidine and thymine

Multiple nonadiabatic pathways greatly facilitate the proceeding of photocrosslinking reactions between 4-thiopyrimidine and thymine.

The effect of DNA backbone on the triplet mechanism of UV-induced thymine-thymine (6-4) dimer formation

Density functional theory calculations were carried out to investigate the formation mechanism of the thymine-thymine (6-4) dimer ((6-4)TT), which is one of the main DNA lesions induced by ultraviolet radiation and is closely related to skin cancers. The DNA backbone was found to have nonnegligible effects on the triplet reaction pathway, particularly the reaction steps involving substantial base rotations. The mechanism for the isomerization from (6-4)TT to its Dewar valence isomer (DewarTT) was also explored, confirming the necessity of absorbing a second photon. In addition, the solvation effects were examined and showed considerable influence on the potential energy surface. Graphical Abstract DFT calculations on the influence of DNA backbone on the mechanism of UV-induced thymine-thymine (6-4) dimer formation.

A role for β,β-xanthophylls in Arabidopsis UV-B photoprotection

Plastidial isoprenoids, such as carotenoids and tocopherols, are important anti-oxidant metabolites synthesized in plastids from precursors generated by the methylerythritol 4-phosphate (MEP) pathway. In this study, we found that irradiation of Arabidopsis thaliana plants with UV-B caused a strong increase in the accumulation of the photoprotective xanthophyll zeaxanthin but also resulted in slightly higher levels of γ-tocopherol. Plants deficient in the MEP enzymes 1-deoxy-D-xylulose 5-phosphate synthase and 1-hydroxy-2-methyl-2-butenyl 4-diphosphate synthase showed a general reduction in both carotenoids and tocopherols and this was associated with increased DNA damage and decreased photosynthesis after exposure to UV-B. Genetic blockage of tocopherol biosynthesis did not affect DNA damage accumulation. In contrast, lut2 mutants that accumulate β,β-xanthophylls showed decreased DNA damage when irradiated with UV-B. Analysis of aba2 mutants showed that UV-B protection was not mediated by ABA (a hormone derived from β,β-xanthophylls). Plants accumulating β,β-xanthophylls also showed decreased oxidative damage and increased expression of DNA-repair enzymes, suggesting that this may be a mechanism for these plants to decrease DNA damage. In addition, in vitro experiments also provided evidence that β,β-xanthophylls can directly protect against DNA damage by absorbing radiation. Together, our results suggest that xanthophyll-cycle carotenoids that protect against excess illumination may also contribute to protection against UV-B.

Photochemical Formation of Cyclobutane Pyrimidine Dimers in DNA through Electron Transfer from a Flanking Base

Electron transfer (ET) to a pyrimidine base from external moieties is a common step involved in the quenching or repair of the cyclobutane pyrimidine dimer (CPD). In contrast, we present a pathway that is initiated by an ET from a flanking guanine base to a pyrimidine base, leading to the formation of a CPD. We studied a T5mCG sequence with a methylated cytosine and our results demonstrate that the pathway involves the formation of an exciplex and intersystem crossings. This pathway also provides an explanation for why the mutational hot spots are correlated with the methylated CpG sequences, which has been a significant issue in cancer research.

Drug-DNA complexation as the key factor in photosensitized thymine dimerization

The crucial role of photosensitizer@DNA complexation in the formation of cyclobutane pyrimidine dimers (CPDs) has been demonstrated using femtosecond and nanosecond transient absorption and emission measurements in combination with in vitro DNA damage assays. This finding opens the door to re-evaluate the mechanisms involved in CPDs photosensitized by other chemicals.

How Does Thymine DNA Survive Ultrafast Dimerization Damage?

The photodimerization reaction between the two adjacent thymine bases within a single strand has been the subject of numerous studies due to its potential to induce DNA mutagenesis and possible tumorigenesis in human skin cells. It is well established that the cycloaddition photoreaction takes place on a picosecond time scale along barrierless or low barrier singlet/triplet pathways. However, the observed dimerization quantum yield in different thymine multimer is considerable lower than might be expected. A reasonable explanation is required to understand why thymine in DNA is able to survive ultrafast dimerization damage. In this work, accurate quantum calculations based on the combined CASPT2//CASSCF/AMBER method were conducted to map the excited state relaxation pathways of the thymine monomer in aqueous solution and of the thymine oligomer in DNA. A monomer-like decay pathway, induced by the twisting of the methyl group, is found to provide a bypass channel to ensure the photostability of thymine in single-stranded oligomers. This fast relaxation path is regulated by the conical intersection between the bright S_CT(¹ππ*) state with the intra-base charge transfer character and the ground state to remove the excess excitation energy, thereby achieving the ground-state recovery with high efficiency.

Quantum Mechanics/Molecular Mechanics Free Energy Maps and Nonadiabatic Simulations for a Photochemical Reaction in DNA: Cyclobutane Thymine Dimer

The absorption of ultraviolet radiation by DNA may result in harmful genetic lesions that affect DNA replication and transcription, ultimately causing mutations, cancer, and/or cell death. We analyze the most abundant photochemical reaction in DNA, the cyclobutane thymine dimer, using hybrid quantum mechanics/molecular mechanics (QM/MM) techniques and QM/MM nonadiabatic molecular dynamics. We find that, due to its double helix structure, DNA presents a free energy barrier between nonreactive and reactive conformations leading to the photolesion. Moreover, our nonadiabatic simulations show that most of the photoexcited reactive conformations return to standard B-DNA conformations after an ultrafast nonradiative decay to the ground state. This work highlights the importance of dynamical effects (free energy, excited-state dynamics) for the study of photochemical reactions in biological systems.

Cyclobutane Thymine Photodimerization Mechanism Revealed by Nonadiabatic Molecular Dynamics

The formation of cyclobutane thymine dimers is one of the most important DNA carcinogenic photolesions induced by ultraviolet irradiation. The long debated question whether thymine dimerization after direct light excitation involves singlet or triplet states is investigated here for the first time using nonadiabatic molecular dynamics simulations. We find that the precursor of this [2 + 2] cycloaddition reaction is the singlet doubly π²π*² excited state, which is spectroscopically rather dark. Excitation to the bright ¹ππ* or dark ¹nπ* excited states does not lead to thymine dimer formation. In all cases, intersystem crossing to the triplet states is not observed during the simulated time, indicating that ultrafast dimerization occurs in the singlet manifold. The dynamics simulations also show that dimerization takes place only when conformational control happens in the doubly excited state.

Quantum Mechanical Studies on the Photophysics and the Photochemistry of Nucleic Acids and Nucleobases

The photophysics and photochemistry of DNA is of great importance due to the potential damage of the genetic code by UV light. Quantum mechanical studies have played a key role in interpretating the results of modern time-resolved pump-probe spectroscopy, and in elucidating the main photoactivated reactive paths. This review provides a concise, complete picture of the computational studies carried out, approximately, in the past decade. We start with an overview of the photophysics of the nucleobases in the gas phase and in solution. We discuss the proposed mechanisms for ultrafast decay to the ground state, that involve conical intersections, consider the role of triplet states, and analyze how the solvent modulates the photophysics. Then we move to larger systems, from dinucleotides to single- and double-stranded oligonucleotides. We focus on the possible role of charge transfer and delocalized or excitonic states in the photophysics of these systems and discuss the main photochemical paths. We finish with an outlook on the current challenges in the field and future directions of research.

Early events of DNA photodamage

Ultraviolet (UV) radiation is a leading external hazard to the integrity of DNA. Exposure to UV radiation triggers a cascade of chemical reactions, and many molecular products (photolesions) have been isolated that are potentially dangerous for the cellular system. The early steps that take place after UV absorption by DNA have been studied by ultrafast spectroscopy. The review focuses on the evolution of excited electronic states, the formation of photolesions, and processes suppressing their formation. Emphasis is placed on lesions involving two thymine bases, such as the cyclobutane pyrimidine dimer, the (6-4) lesion, and its Dewar valence isomer.

Computational reference data for the photochemistry of cyclobutane pyrimidine dimers

The cis-syn cyclobutane pyrimidine dimer is one of the major classes of carcinogenic UV-induced DNA photoproducts. In this work, diverse high-level quantum-chemical methods were used to determine the spectroscopic properties of neutral (singlet and triplet) and charged (cation and anion) species of thymine dimers. Maps of potential energy, charge distribution, electron affinity, and ionization potential of the thymidine dimers were computed along the two dimerization coordinates for neutral and charged species, as well as for the singlet excited state. This set of data aims at providing consistent results computed with the same methods as for photodamage and repair. Based on these results, several different photo-, heat-, and charge-induced mechanisms of dimerization and repair are characterized and discussed. Additionally, a new stable dimer with methylmethylidene-hexahydropyrimidine structure was found in the S0 state.

Physical quenching in competition with the formation of cyclobutane pyrimidine dimers in DNA photolesion

The potential energy profiles toward formation of cyclobutane pyrimidine dimers CPD and the physical quenching after UV excitation were explored for the dinucleotide thymine dinucleoside monophosphate (TpT) using density functional theory (ωB97XD) and the time-dependent density functional theory (TD-ωB97XD). The ωB97XD functional that includes empirical dispersion correction is shown to be an appropriate method to obtain rational results for the current large reaction system of TpT. Photophysical quenching is shown to be predominant over the photochemical CPD formation. Following the initial excitation to the (1)ππ* state, the underlying dark (1)nπ* state bifurcates the excited population to the prevailing IC to S0 and the small ISC to the long-lived triplet state T1 via T4 ((3)ππ*) state that has negligible energy gap with (1)nπ* state. Even for the reactive T1 state, two physical quenching pathways resulting in the conversion back to ground-state reactant via the T1/S0 crossing points are newly located, which are in strong competition with CPD formation. These results provide rationale for the recently observed nanosecond triplet decay rates in the single-stranded (dT)18 and inefficiency of deleterious CPD formation, which allow for a deeper understanding of DNA photostability.

Mechanism of the Decay of Thymine Triplets in DNA Single Strands

The decay of triplet states and the formation of cyclobutane pyrimidine dimers (CPDs) after UV excitation of the all-thymine oligomer (dT)18 and the locked dinucleotide TLpTL were studied by nanosecond IR spectroscopy. IR marker bands characteristic for the CPD lesion and the triplet state were observed from ∼1 ns (time resolution of the setup) onward. The amplitudes of the CPD marker bands remain constant throughout the time range covered (up to 10 μs). The triplet decays with a time constant of ∼10 ns presumably via a biradical intermediate (lifetime ∼60 ns). This biradical has often been invoked as an intermediate for CPD formation via the triplet channel. The present results lend strong support to the existence of this intermediate, yet there is no indication that its decay contributes significantly to CPD formation.

Contribution of cytosine-containing cyclobutane dimers to DNA damage produced by photosensitized triplet-triplet energy transfer

Mutagenic cyclobutane pyrimidine dimers (CPDs) can be induced in DNA through either direct excitation or photosensitized triplet-triplet energy transfer (TTET). In the latter pathway, thymines are expected to receive the excitation energy from the photosensitizer and react with adjacent pyrimidines. By using state-of-the art analytical tools, we provide herein additional information on the formation of cytosine-containing CPDs. We thus determined the yield of all possible CPDs upon TTET in a series of natural DNAs with various base compositions. We show that the distribution of CPDs cannot be explained only by excitation of individual thymines. We propose that the mechanism for TTET involves at least dinucleotides as the minimal targets. The observation of the formation of cytosine-cytosine CPDs also suggests that additional pathways are involved in this photosensitized reaction.

Thymine dimer splitting in the T<>T-G trinucleotide model system: a semiclassical dynamics and TD-DFT study

The mechanism leading to bond cleavage of a thymine-thymine cyclobutane dimer (T<>T) in a model system consisting of the dimer flanked by guanine trinucleotide was studied using semiclassical dynamics simulation. Pulsed laser excitation of the guanine molecule is found to cause electron transfer from the guanine molecule to the dimer, which then dissociates via sequential cleavage of the C5C5' and C6C6' bonds. Subsequently, electrons transfer back to the guanine molecule as the dimer splits into two monomers. The splitting of the cyclobutane dimer was found to be in the femtosecond time scale.

Photoinduced processes in nucleic acids

Photoinduced processes in nucleic acids are phenomena of fundamental interest in diverse fields, from prebiotic studies, through medical research on carcinogenesis, to the development of bioorganic photodevices. In this contribution we survey many aspects of the research across the boundaries. Starting from a historical background, where the main milestones are identified, we review the main findings of the physical-chemical research of photoinduced processes on several types of nucleic-acid fragments, from monomers to duplexes. We also discuss a number of different issues which are still under debate.

On the Formation of Thymine Photodimers in Thymine Single Strands and Calf Thymus DNA

Solar light leads to thymine dimers that are mutagenic and primary cause of skin cancer. Here, we report absorption and synchrotron radiation circular dichroism (CD) spectra of Tn single strands with different number n of bases (n = 2-7, 10, 11) recorded after various 254 nm irradiation times. From a principal component analysis of the CD spectra, we extract fingerprint spectra of both the cyclobutane pyrimidine dimer (CPD) and the pyrimidine (6-4) pyrimidone photoadduct (64PP). Extending the CD measurements to the vacuum ultraviolet region in combination with systematic examinations of size effects is a new approach to gain insight on the dimeric photoproducts. We find a simple linear correlation between n and average number of dimers formed after 1 h of irradiation. The probability for a thymine to engage in a dimer increases from 32% for n = 2 to 41% for n = 11, which implies limited effects of terminal thymines, i.e., the reaction does not occur preferentially at the extremities of the single strands as previously stated. It is even possible to form two dimers with only two bridging thymines. Finally, experiments conducted on calf thymus DNA provided a similar signature of the photodimer, but differences are also evident.

Photophysics and photochemistry of thymine deoxy-dinucleotide in water: a PCM/TD-DFT quantum mechanical study

We here report a fully quantum mechanical study of the main photochemical and photophysical decay routes in aqueous solution of thymine deoxy-dinucleotide (TpT(-) and TpTNa) and of its analogue locked in C3-endo puckering, characterizing five different representative backbone conformers and discussing the chemical physical effects modulating the yield of the different photoproducts. Our approach is based on time-dependent DFT calculations, using the last generation M052X functional, whereas solvent effects are included by means of the polarizable continuum model. Especially when at least one of the sugars adopts C3-endo puckering, a barrierless path on the bright ππ* excitons leads to the S(1)/S(0) crossing region corresponding to the formation of cyclobutane pyrimidine dimer. Charge transfer excited states involving the transfer of an electron from the 5' Thy toward the 3' Thy are involved in the formation of the oxetane intermediate in the path leading to 6-4 pyrimidine pyrimidinone adducts. A non-negligible energy barrier is associated with this latter pathway, which is possible only when one of the two nucleotides adopts C2-endo puckering. Monomer-like decay pathways, involving ππ* or nπ* excited states localized on a single base, are shown to be operative also for loosely stacked bases.

Electronic excited states responsible for dimer formation upon UV absorption directly by thymine strands: joint experimental and theoretical study

The study addresses interconnected issues related to two major types of cycloadditions between adjacent thymines in DNA leading to cyclobutane dimers (T<>Ts) and (6-4) adducts. Experimental results are obtained for the single strand (dT)(20) by steady-state and time-resolved optical spectroscopy, as well as by HPLC coupled to mass spectrometry. Calculations are carried out for the dinucleoside monophosphate in water using the TD-M052X method and including the polarizable continuum model; the reliability of TD-M052X is checked against CASPT2 calculations regarding the behavior of two stacked thymines in the gas phase. It is shown that irradiation at the main absorption band leads to cyclobutane dimers (T<>Ts) and (6-4) adducts via different electronic excited states. T<>Ts are formed via (1)ππ* excitons; [2 + 2] dimerization proceeds along a barrierless path, in line with the constant quantum yield (0.05) with the irradiation wavelength, the contribution of the (3)ππ* state to this reaction being less than 10%. The formation of oxetane, the reaction intermediate leading to (6-4) adducts, occurs via charge transfer excited states involving two stacked thymines, whose fingerprint is detected in the fluorescence spectra; it involves an energy barrier explaining the important decrease in the quantum yield of (6-4) adducts with the irradiation wavelength.

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.

Detailed mechanism for photoinduced cytosine dimerization: a semiclassical dynamics simulation

Semiclassical dynamics simulation is used to study dimerization of two stacked cytosine molecules following excitation by ultrashort laser pulses (25 fs fwhm, Gaussian, 4.1 eV photon energy). The initial excited state was found to form an ultrashort exciton state, which eventually leads to the formation of an excimer state by charge transfer. When the interbase distance, defined as an average value of C(5)-C(5)' and C(6)-C(6)', becomes less than 3 Å, charge recombination occurs due to strong intermolecular interaction, eventually leading to an avoided crossing within 20-30 fs. Geometries at the avoided crossing, with average intermolecular distance of about 2.1 Å, are in accord with CASSCF/CASPT2 calculations. Results indicate that the C(2)-N(1)-C(6)-C(5) and C(2)'-N(1)'-C(6)'-C(5)' dihedral angles' bending vibrations play a significant role in the vibronic coupling between the HOMO and LUMO, which leads to a nonadiabatic transition to the electronic ground state.

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.

Photoinduced dissociation of cyclobutane thymine dimer studied by semiclassical dynamics simulation

Photoinduced dissociation of the thymine dimer is studied in a semiclassical dynamics simulation. The simulation follows excitation of an isolated thymine dimer by a 25 fs fwhm laser pulse, and finds that dissociation proceeds via an asynchronously concerted mechanism, in which the C(5)-C(5)' bond breaks soon after application of the laser pulse, followed by cleavage of the C(6)-C(6)' bond. The dissociation results in two thymine monomers, one in an electronically excited state and the other in the ground state. The former decays to the electronic ground state through an avoided crossing induced by deformation of the pyrimidine ring at the C(5)' and C(6)' sites.

Stabilization of radical anion states of nucleobases in DNA

Trapping of an electron by DNA leads to the formation of radical anion states of pyrimidine bases. Because these states play an important role in biological and chemical processes, their computational treatment is of particular interest. We show that simple electrostatic and quantum chemical models can accurately reproduce the adiabatic electron affinities (EAs) of short DNA stacks recently derived from high-level ab initio calculations (M. Kobylecka, J. Leszczynski, and J. Rak, J. Am. Chem. Soc., 2008, 130, 15683). The electrostatic interaction of an excess electron localized on cytosine or thymine with intra- and inter-strand adjacent nucleobases is found to strongly affect the energy of the radical anions. This interaction is the main origin of the dependence of EA of nucleobases on the nature of neighboring base pairs. In particular, the states XT(-)Y and XC(-)Y, where X and Y = C, T, are, by ca. 0.7 eV, more stable than radical anions GT(-)G and GC(-)G. We find that second-neighbor effects can also significantly modulate EAs, although being smaller than the effects of adjacent bases. The strongest destabilizing effect is found for 5'-GC and 3'-GC, while the 5'-AT base pair stabilizes the radical anion states. Using a combined QM/MD approach, we consider how structural fluctuations of DNA influence the stability of the radical anion states. Despite large dispersions of the stabilization energies due to conformational dynamics of DNA, there are only few thermally accessible structures where GT(-)G and GC(-)G are energetically more favorable than the corresponding pyrimidine triplets. Although stabilization energies calculated for stacks of regular structure are in qualitative agreement with the QM/MD results, structural fluctuations of pi stacks should be taken into account for more accurate description of the excess electron trapped by DNA. The results obtained in this study suggest that simple electrostatic models, in combination with MD simulations, can be very helpful to explore the long time scale behavior of radical anions in DNA.