Ring opening of the cyclobutane in a thymine dimer radical anion

DFT/MM Simulations for Cycloreversion Reaction of Cyclobutane Pyrimidine Dimer with Deprotonated and Protonated E283

DNA photolyase targets the primary ultraviolet (UV)-induced DNA lesion─cyclobutane pyrimidine dimer (CPD), attaches to it, and catalyzes its dissociation. The catalytic mechanism of DNA photolyase and the role of the conserved residue E283 remain subjects of debate. This study employs two-dimensional potential energy surface maps and minimum free energy paths calculated at the ωB97XD/6-31G/MM level to elucidate these mechanisms. Results suggest that the catalytic process follows a sequential, stepwise reaction in which the C5-C5 and C6-C6 bonds are cleaved in order, facilitated by a protonated E283. Activation free energies for these cleavages are calculated at 4.4 and 4.2 kcal·mol-1, respectively. Protonation of E283 reduces electrostatic repulsion with CPD and forms dual hydrogen bonds with it and provides better solvation, stabilizing the CPD radical anion, particularly during intermediate state. This stabilization renders the initial splitting step exergonic, slows reverse reactions of the C5-C5 bond cleavage and electron transfer, and ensures a high quantum yield. Furthermore, the protonation state of E283 significantly affects the type of bond cleavage. Other residues in the active site were also investigated for their roles in the mechanism.

Ab Initio Study of the Stepwise versus Concerted Fragmentation Pathways in Microhydrated Thymine Radical Cations

Thymine radiation-induced fragmentation is characterised by ring opening and the loss of HNCO/NCO. These pathways have been investigated using DFT calculations in the presence of zero, one and two water molecules. In addition to the already characterised stepwise fragmentation mechanism, we propose a novel concerted pathway reported here for the first time. We show that both the stepwise and concerted mechanisms are competitive with activation energies of 2.05 eV and 2.00 eV, respectively, in the gas phase. The effect of microhydration on these mechanisms are examined based on the most stable conformations found by an exploration of the potential energy surface performed by using DFT-based ab initio molecular dynamics. Microhydration is also accompanied by an increase in the activation energies, with respect to gas phase, amounting to 0.47 eV-an increase that is associated to a stabilising effect of water in agreement with recent experimental studies. However, we also point out that this effect is greatly dependent on the specific water arrangement around thymine and could be limited to only 0.13 eV for some configurations.

Exploring Cycloreversion Reaction of Cyclobutane Pyrimidine Dimers Quantum Mechanically

The cyclobutane pyrimidine dimer (CPD) is a major photoproduct of deoxyribonucleic acid (DNA) that is damaged by ultraviolet light. This DNA lesion can be repaired by DNA photolyase with the aid of UV light and two cofactors. To understand the repair mechanism of CPD and whether protonation of CPD participates in the DNA repair process, the cycloreversion reactions of four CPD models and proton transfers between the adjacent residue Glu283 and CPD models were explored through the quantum mechanical method. Two-dimensional maps of potential energy surface in a vacuum and in implicit water solution were calculated at the ωB97XD/6-311++G(2df,2pd) level. One-dimensional potential energy profiles were computed for proton transfer reactions. Among the models that have been considered, both in a vacuum and in water solution, the results indicate that the most likely repair mechanism involves CPD^•2- radical anion splitting in a stepwise manner. C5-C5' splits first, and C6-C6' splits later. The computed free energies of activation of the two splitting steps are 0.9 and 3.1 kcal/mol, respectively. The adjacent Glu283 may stabilize the CPD^•2- radical anion through hydrogen bond and increase the quantum yield; however, protonating the CPD radical anion by Glu283 cannot accelerate the rate of ring opening.

Metal nanoparticle catalyzed cyclobutane cleavage reaction

The pyridine substituted cyclobutane cleavage reaction can be catalyzed directly by metallic silver/gold nanoparticles.

Femtosecond stimulated Raman spectroscopy of the cyclobutane thymine dimer repair mechanism: a computational study

Cyclobutane thymine dimer, one of the major lesions in DNA formed by exposure to UV sunlight, is repaired in a photoreactivation process, which is essential to maintain life. The molecular mechanism of the central step, i.e., intradimer C-C bond splitting, still remains an open question. In a simulation study, we demonstrate how the time evolution of characteristic marker bands (C═O and C═C/C-C stretch vibrations) of cyclobutane thymine dimer and thymine dinucleotide radical anion, thymidylyl(3'→5')thymidine, can be directly probed with femtosecond stimulated Raman spectroscopy (FSRS). We construct a DFT(M05-2X) potential energy surface with two minor barriers for the intradimer C₅-C₅' splitting and a main barrier for the C₆-C₆' splitting, and identify the appearance of two C₅═C₆ stretch vibrations due to the C₆-C₆' splitting as a spectroscopic signature of the underlying bond splitting mechanism. The sequential mechanism shows only absorptive features in the simulated FSRS signals, whereas the fast concerted mechanism shows characteristic dispersive line shapes.

Charge transfer in E. coli DNA photolyase: understanding polarization and stabilization effects via QM/MM simulations

We study fast hole transfer events in E. coli DNA photolyase, a key step in the photoactivation process, using a multiscale computational method that combines nonadiabatic propagation schemes and linear-scaling quantum chemical methods with molecular mechanics force fields. This scheme allows us to follow the time-dependent evolution of the electron hole in an unbiased fashion; that is, no assumptions about hole wave function localization, time scale separation, or adiabaticity of the process have to be made beforehand. DNA photolyase facilitates an efficient long-range charge transport between its flavin adenine dinucleotide (FAD) cofactor and the protein surface via a chain of evolutionary conserved Trp residues on the sub-nanosecond time scale despite the existence of multiple potential trap states. By including a large number of aromatic residues along the charge transfer pathway into the quantum description, we are able to identify the main pathway among alternative possible routes. The simulations show that charge transfer, which is extremely fast in this protein, occurs on the same time scale as the protein response to the electrostatic changes; that is, time-scale separation as often presupposed in charge transfer studies seems to be inappropriate for this system. Therefore, coupled equations of motion, which propagate electrons and nuclei simultaneously, appear to be necessary. The applied computational model is shown to capture the essentials of the reaction kinetics and thermodynamics while allowing direct simulations of charge transfer events on their natural time scale.

Thymine dimer photoreversal in purine-containing trinucleotides

Cyclobutane-pyrimidine dimer yields in UV-irradiated DNA are controlled by the equilibrium between forward and reverse photoreactions. Past studies have shown that dimer yields are suppressed at sites adjacent to a purine base, but the underlying causes are unclear. In order to investigate whether this suppression is the result of repair by electron transfer from a neighboring nucleobase, the yields and dynamics of the reverse reaction were studied using trinucleotides containing a cis-syn dimer (T<>T) flanked on the 5' or the 3' side by adenine or guanine. The probability of forming an excited state on T<>T or on the purine base was varied by tuning the irradiation wavelength between 240 and 280 nm. Cleavage quantum yields decrease by an order of magnitude over this wavelength range and are less than 1% at 280 nm, a wavelength that excites the purine base with more than 95% probability. Conditional quantum yields of cleavage for the trinucleotides given excitation of T<>T are similar in magnitude to the quantum yield of cleavage of unmodified T<>T. These results indicate that within experimental uncertainty all photoreversal in these single-stranded substrates is the result of direct electronic excitation of T<>T. Photolyase-like repair of T<>T due to electron transfer from an adjacent purine is negligible in these substrates. Instead, the observed variation in photoreversal quantum yields for adenine- versus guanine-flanked cis-syn dimer could be due to uncertainties in absorption cross sections or to a modest quenching effect by the purine on the excited state of T<>T. Pump-probe measurements reveal that the excited-state lifetimes of A or G in the dimer-containing trinucleotides are unperturbed by the neighboring dimer, indicating that electron transfer from purine base to T<>T is not competitive with rapid excited-state deactivation. Pump-probe measurements on unmodified T<>T in aqueous solution indicate that cleavage is most likely complete on a picosecond or subpicosecond time scale.

Kinetics of cyclobutane thymine dimer splitting by DNA photolyase directly monitored in the UV

CPD photolyase uses light to repair cyclobutane pyrimidine dimers (CPDs) formed between adjacent pyrimidines in UV-irradiated DNA. The enzyme harbors an FAD cofactor in fully reduced state (FADH(-)). The CPD repair mechanism involves electron transfer from photoexcited FADH(-) to the CPD, splitting of its intradimer bonds, and electron return to restore catalytically active FADH(-). The two electron transfer processes occur on time scales of 10(-10) and 10(-9) s, respectively. Until now, CPD splitting itself has only been poorly characterized by experiments. Using a previously unreported transient absorption setup, we succeeded in monitoring cyclobutane thymine dimer repair in the main UV absorption band of intact thymine at 266 nm. Flavin transitions that overlay DNA-based absorption changes at 266 nm were monitored independently in the visible and subtracted to obtain the true repair kinetics. Restoration of intact thymine showed a short lag and a biexponential rise with time constants of 0.2 and 1.5 ns. We assign these two time constants to splitting of the intradimer bonds (creating one intact thymine and one thymine anion radical T(∘-)) and electron return from T(∘-) to the FAD cofactor with recovery of the second thymine, respectively. Previous model studies and computer simulations yielded various CPD splitting times between < 1 ps and < 100 ns. Our experimental results should serve as a benchmark for future efforts to model enzymatic photorepair. The technique and methods developed here may be applied to monitor other photoreactions involving DNA.

Exploration of photochemical reactions of N-trimethylsilylmethyl-substituted uracil, pyridone, and pyrrolidone derivatives

Photochemical reactions of N-trimethylsilylmethyl-substituted uracil, pyridone and pyrrolidone derivatives were carried out to determine if silicone containing substituents have an impact on excited state reaction profiles. The results show that ultraviolet irradiation of N-trimethylsilylmethyl substituted uracils in the presence of substituted alkenes leads to efficient formation of both dimeric and cross [2+2]-cycloaddition products. Qualitatively similar observations were made in a study of the photochemistry of N-trimethylsilylmethyl-2-pyridone. The combined results demonstrate that [2+2]-photocycloaddition is a more efficient excited state reaction pathway for the uracil and pyridone substrates as compared to other processes, such as ylide-forming trimethylsilyl group C-to-O migration. Finally, photoreactions of N-trimethylsilylmethyl-2-pyrrolidone in solutions containing dipolarophiles, such as methyl acrylate, lead to the formation of the desilylation product, N-methyl-2-pyrrolidone by way of a simple, non-ylide generating, protodesilylation process. In addition, observations were made which show that orbital symmetry allowed photocycloreversion reactions of dimeric uracil derivatives, involving cyclobutane ring splitting, to take place. These processes, which lead to the formation of monomeric uracils, appear to be stimulated by the presence of electron donor groups on the cyclobutane ring, a likely result of a new SET promoted cyclobutane ring cleavage pathway. In the cases of N-trimethylsilylmethyl-substituted cyclobutane derivatives that possess phthalimide groups, highly efficient excited state cleavage of the cyclobutane moiety occurs to produce uracil derivatives and corresponding vinyl phthalimide.

An AIMD study of the CPD repair mechanism in water: reaction free energy surface and mechanistic implications

In a series of two papers, we report the detailed mechanism of cyclobutane pyrimidine dimer repair in aqueous solvent using ab initio molecular dynamics simulations (AIMD). Umbrella sampling is used to determine the free energy surface for dimer splitting. The two-dimensional free energy surface for splitting of the C5-C5' and C6-C6' bonds on the anion surface is reported. The splitting of the C5-C5' and C6-C6' bonds occurs on a picosecond time scale. The transition state along the splitting coordinate in the anion state coincides with a maximum in the free energy along the same coordinate on the neutral surface. The implication is that back electron transfer occurring before the anion reaches the transition state leads to reformation of the cyclobutane dimer, while back electron transfer after transit through the transition state, leads to successful repair. On the basis of our calculations for CPD splitting in water, we propose a framework for understanding how various factors, such as solvent polarity, can control repair efficiency. This framework explains why back electron transfer leads predominantly to unsuccessful repair in some situations, and successful repair in others. A key observation is that the same free energy surfaces that control dimer splitting also govern how the back electron transfer rate changes during the splitting process. Configurational changes of the dimer along the splitting coordinate are also documented.

Investigation of excess-electron transfer in DNA double-duplex systems allows estimation of absolute excess-electron transfer and CPD cleavage rates

To investigate the parameters and rates that determine excess-electron transfer processes in DNA duplexes, we developed a DNA double-duplex system containing a reduced and deprotonated flavin donor at the junction of two duplexes with either the same or different electron acceptors in the individual duplex substructures. This model system allows us to bring the two electron acceptors in the duplex substructures into direct competition for injected electrons and this enables us to decipher how the kind of acceptor influences the transfer data. Measurements with the electron acceptors 8-bromo-dA (BrdA), 8-bromo-dG (BrdG), 5-bromo-dU (BrdU), and a cyclobutane pyrimidine dimer, which is a UV-induced DNA lesion, allowed us to obtain directly the maximum overall reaction rates of these acceptors and especially of the T=T dimer with the injected electrons in the duplex. In line with previous observations, we detected that the overall dimer cleavage rate is about one order of magnitude slower than the debromination of BrdU. Furthermore, we present a more detailed explanation of why sequence dependence cannot be observed when a T=T dimer is used as the acceptor and we estimate the absolute excess-electron hopping rates.

Reaction mechanisms of DNA photolyase

DNA photolyase uses visible light and a fully reduced flavin cofactor FADH(-) to repair major UV-induced lesions in DNA, the cyclobutane pyrimidine dimers (CPDs). Electron transfer from photoexcited FADH(-) to CPD, splitting of the two intradimer bonds, and back electron transfer to the transiently formed flavin radical FADH° occur in overall 1ns. Whereas the kinetics of FADH° was resolved, the DNA-based intermediates escaped unambiguous detection yet. Another light reaction, named photoactivation, reduces catalytically inactive FADH° to FADH(-) without implication of DNA. It involves electron hopping along a chain of three tryptophan residues in 30ps, as elucidated in detail by transient absorption spectroscopy. The same triple tryptophan chain is found in cryptochrome blue-light photoreceptors and may be involved in their primary photoreaction.

Formation of pyrimidine dimer radical anions in the gas phase

Crossed-beam experiments revealed that attachment of a free electron to the cyclobutane pyrimidine dimers c,s-DMT<>DMT and c,a-DMT<>DMT leads to the formation of dimer radical anions with the lifetime of at least 80 micros, thus showing that the latter are much more stable than previously believed.

Computational study of thymine dimer radical anion splitting in the self-repair process of duplex DNA

Formation of the thymine dimer is one of the most important types of photochemical damage in DNA, responsible for several biological pathologies. Though specifically designed proteins (photolyases) can efficiently repair this type of damage in living cells, an autocatalytic activity of the DNA itself was recently discovered, allowing for a self-repair mechanism. In this paper, we provide the first molecular dynamics study of the splitting of thymine dimer radical anions, using a quantum mechanical/molecular mechanics (QM/MM) approach based on density functional theory (DFT) to describe the quantum region. A set of seven statistically representative molecular dynamics trajectories is analyzed. Our calculations predict an asynchronously concerted process in which C5-C5' bond breaking is barrierless while C6-C6' bond breaking is characterized by a small free energy barrier. An upper bound of 2.5 kcal/mol for this barrier is estimated. Moreover, the molecular dynamics study and the low free energy barrier involved in C6-C6' bond breaking characterize the full process as being an ultrafast reaction.