Ultraviolet light is strongly absorbed by DNA, producing excited electronic states that sometimes initiate damaging photochemical reactions. Fully mapping the reactive and nonreactive decay pathways available to excited electronic states in DNA is a decades-old quest. Progress toward this goal has accelerated rapidly in recent years, in large measure because of ultrafast laser experiments. Here we review recent discoveries and controversies concerning the nature and dynamics of excited states in DNA model systems in solution. Nonradiative decay by single, solvated nucleotides occurs primarily on the subpicosecond timescale. Surprisingly, excess electronic energy relaxes one or two orders of magnitude more slowly in DNA oligo- and polynucleotides. Highly efficient nonradiative decay pathways guarantee that most excited states do not lead to deleterious reactions but instead relax back to the electronic ground state. Understanding how the spatial organization of the bases controls the relaxation of excess electronic energy in the double helix and in alternative structures is currently one of the most exciting challenges in the field.

It was recently shown that thymine dimers in the all-thymine oligonucleotide (dT)(18) are fully formed in <1 ps after ultraviolet excitation. The speed and low quantum yield of this reaction suggest that only a small fraction of the conformers of this structurally disordered oligonucleotide are in a position to react at the instant of photon absorption. In this work, we explore the hypothesis that conventional molecular dynamics simulations can be used to predict the yield of cyclobutane pyrimidine dimers in DNA. Conformations obtained from simulations of thymidylyl-(3'-5')-thymidine in various cosolvents were classified as dimerizable or nondimerizable depending on the distance between the C5-C6 double bonds of the adjacent thymine bases and the torsion angle between them. The quantum yield of cyclobutane pyrimidine dimer formation was calculated as the number of dimerizable conformations divided by the total number of conformations. The experimental quantum yields measured in the different solvents were satisfactorily reproduced using physically reasonable values for the two parameters. The mean dimerizable structure computed by averaging all of the dimerizable cis-syn conformations is structurally similar to the actual cis-syn dimer. Compared to the canonical B-form TT step, the most important structural property of a dimerizable conformation is its reduced helical twist angle of 22 degrees.

The kinetics of thymine-thymine cyclobutane pyrimidine dimer (TT-CPD) formation was studied at 23 thymine-thymine base steps in two 247-base pair DNA sequences irradiated at 254 nm. Damage was assayed site-specifically and simultaneously on both the forward and reverse strands by detecting emission from distinguishable fluorescent labels at the 5’-termini of fragments cleaved at CPD sites by T4 pyrimidine dimer glycosylase and separated by gel electrophoresis. The total DNA strand length of nearly 1000 bases made it possible to monitor damage at all 9 tetrads of the type XTTY, where X and Y are non-thymine bases. TT-CPD yields for different tetrads were found to vary by as much as an order of magnitude, but similar yields were observed at all instances of a given tetrad. Kinetic analysis of CPD formation at 23 distinct sites reveals that both the formation and reversal photoreactions depend sensitively on the identity of the nearest-neighbour bases on the 5’ and 3’ side of a photoreactive TT base step. The lowest formation and reversal rates occur when two purine bases flank a TT step, while the highest formation and reversal rates are observed for tetrads with at least one flanking C. Overall, the results show that the probabilities of CPD formation and photoreversal depend principally on interactions with nearest-neighbour bases.

In this work, we examine the importance of nuclear quantum effects on capturing the line broadening and vibronic structure of optical spectra. We determine the absorption spectra of three aromatic molecules indole, pyridine, and benzene using time dependent density functional theory with several molecular dynamics sampling protocols: force-field based empirical potentials, ab initio simulations, and finally path-integrals for the inclusion of nuclear quantum effects. We show that the absorption spectrum for all these chromophores are similarly broadened in the presence of nuclear quantum effects regardless of the presence of hydrogen bond donor or acceptor groups. We also show that simulations incorporating nuclear quantum effects are able to reproduce the heterogeneous broadening of the absorption spectra even with empirical force fields. The spectral broadening associated with nuclear quantum effects can be accounted for by the broadened distribution of chromophore size as revealed by a particle in the box model. We also highlight the role that nuclear quantum effects have on the underlying electronic structure of aromatic molecules as probed by various electrostatic properties.