Molecular insight into how the position of an abasic site modifies DNA duplex stability and dynamics.

Local perturbations to DNA base-pairing stability from lesions and chemical modifications can alter the stability and dynamics of an entire oligonucleotide. End effects may cause the position of a disruption within a short duplex to influence duplex stability and structural dynamics, yet this aspect of nucleic acid modifications is often overlooked. We investigate how the position of an abasic site (AP site) impacts the stability and dynamics of short DNA duplexes. Using a combination of steady-state and time-resolved spectroscopy and molecular dynamics simulations, we unravel an interplay between AP-site position and nucleobase sequence that controls energetic and dynamic disruption to the duplex. The duplex is disrupted into two segments by an entropic barrier for base-pairing on each side of the AP site. The barrier induces fraying of the short segment when an AP site is near the termini. Shifting the AP site inward promotes a transition from short-segment fraying to fully encompassing the barrier into the thermodynamics of hybridization, leading to further destabilization of the duplex. Nucleobase sequence determines the length scale for this transition by tuning the barrier height and base-pair stability of the short segment, and certain sequences enable out-of-register base-pairing to minimize the barrier height.

Molecular insight into how the position of an abasic site and its sequence environment influence DNA duplex stability and dynamics.

Local perturbations to DNA base-pairing stability from lesions and chemical modifications can alter the stability and dynamics of an entire oligonucleotide. End effects may cause the position of a disruption within a short duplex to influence duplex stability and structural dynamics, yet this aspect of nucleic acid modifications is often overlooked. We investigate how the position of an abasic site (AP site) impacts the stability and dynamics of short DNA duplexes. Using a combination of steady-state and time-resolved spectroscopy and molecular dynamics simulations, we unravel an interplay between AP-site position and nucleobase sequence that controls energetic and dynamic disruption to the duplex. The duplex is disrupted into two segments by an entropic barrier for base pairing on each side of the AP site. The barrier induces fraying of the short segment when an AP site is near the termini. Shifting the AP site inward promotes a transition from short-segment fraying to fully encompassing the barrier into the thermodynamics of hybridization, leading to further destabilization the duplex. Nucleobase sequence determines the length scale for this transition by tuning the barrier height and base-pair stability of the short segment, and certain sequences enable out-of-register base pairing to minimize the barrier height.

Thermodynamics and kinetics of DNA and RNA dinucleotide hybridization to gaps and overhangs.

Hybridization of short nucleic acid segments (<4 nt) to single-strand templates occurs as a critical intermediate in processes such as nonenzymatic nucleic acid replication and toehold-mediated strand displacement. These templates often contain adjacent duplex segments that stabilize base pairing with single-strand gaps or overhangs, but the thermodynamics and kinetics of hybridization in such contexts are poorly understood because of the experimental challenges of probing weak binding and rapid structural dynamics. Here we develop an approach to directly measure the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization using steady-state and temperature-jump infrared spectroscopy. Our results suggest that dinucleotide binding is stabilized through coaxial stacking interactions with the adjacent duplex segments as well as from potential noncanonical base-pairing configurations and structural dynamics of gap and overhang templates revealed using molecular dynamics simulations. We measure timescales for dissociation ranging from 0.2-40 µs depending on the template and temperature. Dinucleotide hybridization and dehybridization involve a significant free energy barrier with characteristics resembling that of canonical oligonucleotides. Together, our work provides an initial step for predicting the stability and kinetics of hybridization between short nucleic acid segments and various templates.

Direct monitoring of the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization from gaps and overhangs.

Hybridization of short nucleic acid segments (<4 nucleotides) to single-strand templates occurs as a critical intermediate in processes such as non-enzymatic nucleic acid replication and toehold-mediated strand displacement. These templates often contain adjacent duplex segments that stabilize base pairing with single-strand gaps or overhangs, but the thermodynamics and kinetics of hybridization in such contexts are poorly understood due to experimental challenges of probing weak binding and rapid structural dynamics. Here we develop an approach to directly measure the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization using steady-state and temperature-jump infrared spectroscopy. Our results suggest that dinucleotide binding is stabilized through coaxial stacking interactions with the adjacent duplex segments as well as from potential non-canonical base pairing configurations and structural dynamics of gap and overhang templates revealed using molecular dynamics simulations. We measure timescales for dissociation ranging from 0.2 to 40 µs depending on the template and temperature. Dinucleotide hybridization and dehybridization involves a significant free energy barrier with characteristics resembling that of canonical oligonucleotides. Together, our work provides an initial step for predicting the stability and kinetics of hybridization between short nucleic acid segments and various templates.

Disruption of energetic and dynamic base pairing cooperativity in DNA duplexes by an abasic site.

DNA duplex stability arises from cooperative interactions between multiple adjacent nucleotides that favor base pairing and stacking when formed as a continuous stretch rather than individually. Lesions and nucleobase modifications alter this stability in complex manners that remain challenging to understand despite their centrality to biology. Here, we investigate how an abasic site destabilizes small DNA duplexes and reshapes base pairing dynamics and hybridization pathways using temperature-jump infrared spectroscopy and coarse-grained molecular dynamics simulations. We show how an abasic site splits the cooperativity in a short duplex into two segments, which destabilizes small duplexes as a whole and enables metastable half-dissociated configurations. Dynamically, it introduces an additional barrier to hybridization by constraining the hybridization mechanism to a step-wise process of nucleating and zipping a stretch on one side of the abasic site and then the other.

Determining Sequence-Dependent DNA Oligonucleotide Hybridization and Dehybridization Mechanisms Using Coarse-Grained Molecular Simulation, Markov State Models, and Infrared Spectroscopy.

A robust understanding of the sequence-dependent thermodynamics of DNA hybridization has enabled rapid advances in DNA nanotechnology. A fundamental understanding of the sequence-dependent kinetics and mechanisms of hybridization and dehybridization remains comparatively underdeveloped. In this work, we establish new understanding of the sequence-dependent hybridization/dehybridization kinetics and mechanism within a family of self-complementary pairs of 10-mer DNA oligomers by integrating coarse-grained molecular simulation, machine learning of the slow dynamical modes, data-driven inference of long-time kinetic models, and experimental temperature-jump infrared spectroscopy. For a repetitive ATATATATAT sequence, we resolve a rugged dynamical landscape comprising multiple metastable states, numerous competing hybridization/dehybridization pathways, and a spectrum of dynamical relaxations. Introduction of a G:C pair at the terminus (GATATATATC) or center (ATATGCATAT) of the sequence reduces the ruggedness of the dynamics landscape by eliminating a number of metastable states and reducing the number of competing dynamical pathways. Only by introducing a G:C pair midway between the terminus and the center to maximally disrupt the repetitive nature of the sequence (ATGATATCAT) do we recover a canonical "all-or-nothing" two-state model of hybridization/dehybridization with no intermediate metastable states. Our results establish new understanding of the dynamical richness of sequence-dependent kinetics and mechanisms of DNA hybridization/dehybridization by furnishing quantitative and predictive kinetic models of the dynamical transition network between metastable states, present a molecular basis with which to understand experimental temperature jump data, and furnish foundational design rules by which to rationally engineer the kinetics and pathways of DNA association and dissociation for DNA nanotechnology applications.

Temperature-Jump 2D IR Spectroscopy with Intensity-Modulated CW Optical Heating.

Pulsed temperature-jump (T-jump) spectroscopy with infrared (IR) detection has been widely used to study biophysical processes occurring from nanoseconds to ∼1 ms with structural sensitivity. However, many systems exhibit structural dynamics on time scales longer than the millisecond barrier that is set by the time scale for thermal relaxation of the sample. We developed a linear and nonlinear infrared spectrometer coupled to an intensity-modulated continuous wave (CW) laser to probe T-jump-initiated chemical reactions from <1 ms to seconds. Time-dependent modulation of the CW laser leads to a <1 ms heating time as well as a constant final temperature (±3%) for the duration of the heating time. Temperature changes of up to 75 °C in D2O are demonstrated, allowing for nonequilibrium measurements inaccessible to standard pulsed optical T-jump setups. T-jump linear absorption, pump-probe, and two-dimensional IR (2D IR) spectroscopy are applied to the unfolding and refolding of ubiquitin and a model intercalated motif (i-motif) DNA sequence, and analysis of the observed signals is used to demonstrate the limits and utility of each method. Overall, the ability to probe temperature-induced chemical processes from <1 ms to many seconds with 2D IR spectroscopy provides multiple new avenues for time-dependent spectroscopy in chemistry and biophysics.

Oxidized Derivatives of 5-Methylcytosine Alter the Stability and Dehybridization Dynamics of Duplex DNA.

The naturally occurring nucleobase 5-methylcytosine (mC) and its oxidized derivatives 5-hydroxymethylcytosine (hmC), 5-formylcytosine (fC), and 5-carboxylcytosine (caC) play important roles in epigenetic regulation and, along with cytosine (C), represent nucleobases currently implicated in the active cytosine demethylation pathway. Despite considerable interest in these modified bases, their impact on the thermodynamic stability of double-stranded DNA (dsDNA) remains ambiguous and their influence on hybridization kinetics and dynamics is even less well-understood. To address these unknowns, we employ steady-state and time-resolved infrared spectroscopy to measure the influence of cytosine modification on the thermodynamics and kinetics of hybridization by assessing the impact on local base pairing dynamics, shifts in the stability of the duplex state, and changes to the hybridization transition state. Modification with mC leads to more tightly bound base pairing below the melting transition and stabilizes the duplex relative to canonical DNA, but the free energy barrier to dehybridization at physiological temperature is nevertheless reduced slightly. Both hmC and fC lead to an increase in local base pair fluctuations, a reduction in the cooperativity of duplex melting, and a lowering of the dissociation barrier, but these effects are most pronounced when the 5-position is formylated. The caC nucleobase demonstrates little impact on dsDNA under neutral conditions, but we find that this modification can dynamically switch between C-like and fC-like behavior depending on the protonation state of the 5-position carboxyl group. Our results provide a consistent thermodynamic and kinetic framework with which to describe the modulation of the physical properties of double-stranded DNA containing these modified nucleobases.

5-Carboxylcytosine and Cytosine Protonation Distinctly Alter the Stability and Dehybridization Dynamics of the DNA Duplex.

Applications associated with nucleobase protonation events are grounded in their fundamental impact on DNA thermodynamics, structure, and hybridization dynamics. Of the canonical nucleobases, N3 protonation of cytosine (C) is the most widely utilized in both biology and nanotechnology. Naturally occurring C derivatives that shift the N3 pKa introduce an additional level of tunability. The epigenetic nucleobase 5-carboxylcytosine (caC) presents a particularly interesting example since this derivative forms Watson-Crick base pairs of similar stability and displays pH-dependent behavior over the same range as the canonical nucleobase. However, the titratable group in caC corresponds to the exocyclic carboxyl group rather than N3, and the implications of these divergent protonation events toward DNA hybridization thermodynamics, kinetics, and base pairing dynamics remain poorly understood. Here, we study the pH dependence of these physical properties using model oligonucleotides containing C and caC with FTIR and temperature-jump IR spectroscopy. We demonstrate that N3 protonation of C completely disrupts duplex stability, leading to large shifts in the duplex/single-strand equilibrium, a reduction in the cooperativity of melting, and an acceleration in the rate of duplex dissociation. In contrast, while increasing 5-carboxyl protonation in caC-containing duplexes induces an increase in base pair fluctuations, the DNA duplex can tolerate substantial protonation without significant perturbation to the duplex/single-strand equilibrium. However, 5-carboxyl protonation has a large impact on hybridization kinetics by reducing the transition state free energy. Our thermodynamic and kinetic analysis provides new insight on the impact of two divergent protonation mechanisms in naturally occurring nucleobases on the biophysical properties of DNA.

Photochemical and Photodynamical Properties of Sulfur-Substituted Nucleic Acid Bases.

Sulfur-substituted nucleobases (a.k.a., thiobases) are among the world's leading prescriptions for chemotherapy and immunosuppression. Long-term treatment with azathioprine, 6-mercaptopurine and 6-thioguanine has been correlated with the photoinduced formation of carcinomas. Establishing an in-depth understanding of the photochemical properties of these prodrugs may provide a route to overcoming these carcinogenic side effects, or, alternatively, a basis for developing effective compounds for targeted phototherapy. In this review, a broad examination is undertaken, surveying the basic photochemical properties and excited-state dynamics of sulfur-substituted analogs of the canonical DNA and RNA nucleobases. A molecular-level understanding of how sulfur substitution so remarkably perturbs the photochemical properties of the nucleobases is presented by combining experimental results with quantum-chemical calculations. Structure-property relationships demonstrate the impact of site-specific sulfur substitution on the photochemical properties, particularly on the population of the reactive triplet state. The value of fundamental photochemical investigations for driving the development of ultraviolet-A chemotherapeutics is showcased. The most promising photodynamic agents identified thus far have been investigated in various carcinoma cell lines and shown to decrease cell proliferation upon exposure to ultraviolet-A radiation. Overarching principles have been elucidated for the impact that sulfur substitution of the carbonyl oxygen has on the photochemical properties of the nucleobases.

Photochemical relaxation pathways of S6-methylthioinosine and O6-methylguanosine in solution.

S6-Methylthioinosine and O6-methylguanosine are byproducts resulting from the enzymatic reactions of sulfur-substituted prodrugs in cells and from the interaction of alkylating agents with cellular DNA, respectively. Their photochemistry has not been investigated, and it is currently unknown whether light absorption by these byproducts may pose any threat to the cell. In this contribution, their photoinduced processes upon absorption of UVB radiation are reported using broadband transient absorption spectroscopy. Plausible electronic relaxation mechanisms are proposed for both biological molecules, which are supported by steady-state absorption and emission measurements, and by singlet and triplet vertical excitation energies performed on a large subset of ground-state optimized conformational isomers in solution. The results are compared to the body of knowledge gathered in the scientific literature about the light-induced processes in the sulfur-substituted and canonical purine monomers. In particular, it is shown that S6-methylation decreases the rate to populate the lowest-energy triplet state and blueshifts the ground-state absorption spectrum compared to those for the sulfur-substituted prodrugs and for the 6-thioguanosine metabolite. Similarly, O6-methylation decreases the rate of internal conversion to the ground state observed in the guanine monomers by more than 10-fold in acetonitrile and 40-fold in aqueous solution, while it redshifts the ground-state absorption spectrum. Collectively, this investigation provides relevant new insights about the relationship between structural modifications of the purine chromophore and the electronic relaxation mechanisms in this important group of biological molecules.

Excited-State Dynamics in O6-Methylguanosine: Impact of O6-Methylation on the Relaxation Mechanism of Guanine Monomers.

Absorption of ultraviolet radiation by DNA bases results in ultrafast internal conversion to the ground state, which minimizes photodamage. However, exogenous and endogenous alkylating agents present in the cellular environment can methylate the nucleobases in DNA. In particular, methylation of guanosine at the O6 position in DNA leads to the formation of the O6-methylguanosine adduct, which may alter the photostability of DNA. This contribution demonstrates that O6-methylation of guanosine red shifts its ground-state absorption spectrum and slows down the rate of internal conversion to the ground state by ∼40-fold in aqueous solution. The 40-fold decrease in the rate of excited-state decay increases the probability of photodamage within cellular DNA. It is proposed that the longer decay lifetime corresponds to relaxation of the excited-state population in O6-methylguanosine along a C6-puckered reaction coordinate in the 1ππ*(La) potential energy surface that runs parallel to an ultrafast internal conversion pathway along a C2-puckered coordinate.

2-Thiouracil intersystem crossing photodynamics studied by wavelength-dependent photoelectron and transient absorption spectroscopies.

Single-atom substitution within a natural nucleobase-such as replacing oxygen by sulfur in uracil-can result in drastic changes in the relaxation dynamics after UV excitation. While the photodynamics of natural nucleobases like uracil are dominated by pathways along singlet excited states, the photodynamics of thiobases like 2-thiouracil populate the triplet manifold with near unity quantum yield. In the present study, a synergistic approach based on time-resolved photoelectron spectroscopy (TRPES), time-resolved absorption spectroscopy (TRAS), and ab initio computations has been particularly successful at unraveling the underlying photophysical principles and describing the dissimilarities between the natural and substituted nucleobases. Specifically, we find that varying the excitation wavelength leads to differences between gas-phase and condensed-phase experimental results. Systematic trends are observed in the intersystem crossing time constants with varying excitation wavelength, which can be readily interpreted in the context of ab initio calculations performed both in vacuum and including solvent effects. Thus, the combination of TRPES and TRAS experiments with high-level computational techniques allows us to characterize the topology of the potential energy surfaces defining the relaxation dynamics of 2-thiouracil in both gas and condensed phases, as well as investigate the accessibility of conical intersections and crossings, and potential energy barriers along the associated relaxation coordinates.

Photochemical Reactivity of dTPT3: A Crucial Nucleobase Derivative in the Development of Semisynthetic Organisms.

In 2017, the Romesberg group successfully developed the dTPT3·dNaM unnatural base pair to create a semisynthetic organism with enhanced genetic fidelity and the ability to store additional genetic information indefinitely. It is also desirable that the newly developed genetic material remains stable upon exposure to radiation. However, the photochemical properties of dTPT3 are presently unknown. In this contribution, excitation of dTPT3 with near-visible radiation is shown to efficiently populate a reactive triplet state in high yield and on a sub-1 ps time scale; a state that is able to survive for up to a few microseconds in aqueous solution. The triplet state can also generate singlet oxygen in ca. 30% yield, suggesting that dTPT3 can act as a pervasive photosensitizer to accelerate oxidatively generated damage within DNA and to other biological molecules within cells.

Solvatochromic Effects on the Absorption Spectrum of 2-Thiocytosine.

The solvatochromic effects of six different solvents on the UV absorption spectrum of 2-thiocytosine have been studied by a combination of experimental and theoretical techniques. The steady-state absorption spectra show significant shifts of the absorption bands, where in more polar solvents the first absorption maximum shifts to higher transition energies and the second maximum to lower energies. The observed solvatochromic shifts have been rationalized using three popular solvatochromic scales and with high-level multireference quantum chemistry calculations including implicit and explicit solvent effects. It has been found that the dipole moments of the excited states account for some general shifts in the excitation energies, whereas the explicit solvent interactions explain the differences in the spectra recorded in the different solvents.

Excited-State Dynamics of the Thiopurine Prodrug 6-Thioguanine: Can N9-Glycosylation Affect Its Phototoxic Activity?

6-Thioguanine, an immunosuppressant and anticancer prodrug, has been shown to induce DNA damage and cell death following exposure to UVA radiation. Its metabolite, 6-thioguanosine, plays a major role in the prodrug's overall photoreactivity. However, 6-thioguanine itself has proven to be cytotoxic following UVA irradiation, warranting further investigation into its excited-state dynamics. In this contribution, the excited-state dynamics and photochemical properties of 6-thioguanine are studied in aqueous solution following UVA excitation at 345 nm in order to provide mechanistic insight regarding its photochemical reactivity and to scrutinize whether N9-glycosylation modulates its phototoxicity in solution. The experimental results are complemented with time-dependent density functional calculations that include solvent dielectric effects by means of a reaction-field solvation model. UVA excitation results in the initial population of the S2(ππ*) state, which is followed by ultrafast internal conversion to the S1(nπ*) state and then intersystem crossing to the triplet manifold within 560 ± 60 fs. A small fraction (ca. 25%) of the population that reaches the S1(nπ*) state repopulates the ground state. The T1(ππ*) state decays to the ground state in 1.4 ± 0.2 µs under N2-purged conditions, using a 0.2 mM concentration of 6-thioguanine, or it can sensitize singlet oxygen in 0.21 ± 0.02 and 0.23 ± 0.02 yields in air- and O2-saturated solution, respectively. This demonstrates the efficacy of 6-thioguanine to act as a Type II photosensitizer. N9-glycosylation increases the rate of intersystem crossing from the singlet to triplet manifold, as well as from the T1(ππ*) state to the ground state, which lead to a ca. 40% decrease in the singlet oxygen yield under air-saturated conditions. Enhanced vibronic coupling between the singlet and triplet manifolds due to a higher density of vibrational states is proposed to be responsible for the observed increase in the rates of intersystem crossing in 6-thioguanine upon N9-glycosylation.

Can a Six-Letter Alphabet Increase the Likelihood of Photochemical Assault to the Genetic Code?

In 2014, two unnatural nucleosides, d5SICS and dNaM, were shown to selectively base pair and replicate with high fidelity in a modified strain of E. coli, thus effectively expanding its genetic alphabet from four to six letters. More recently, a significant reduction in cell proliferation was reported in cells cultured with d5SICS, and putatively with dNaM, upon exposure to brief periods of near-visible radiation. The photosensitizing properties of the lowest-energy excited triplet state of both d5SICS and dNaM were implicated in their cytotoxicity. Importantly, however, the excited-state mechanisms by which near-visible excitation populates the triplet states of d5SICS and dNaM are currently unknown. In this study, steady-state and time-resolved spectroscopies are combined with quantum-chemical calculations in order to reveal the excited-state relaxation mechanisms leading to efficient population of the triplet states in these unnatural nucleosides in solution. It is shown that excitation of d5SICS or dNaM with near-visible light leads overwhelmingly to ultrafast population of their triplet states on the femtosecond time scale. The results presented in this work lend strong support to the proposal that photoexcitation of these unnatural nucleosides can accelerate oxidatively generated damage to DNA and other biomolecules within the cellular environment.

Unintended Consequences of Expanding the Genetic Alphabet.

The base pair d5SICS·dNaM was recently reported to incorporate and replicate in the DNA of a modified strain of Escherichia coli, thus making the world's first stable semisynthetic organism. This newly expanded genetic alphabet may allow organisms to store considerably more information in order to translate proteins with unprecedented enzymatic activities. Importantly, however, there is currently no knowledge of the photochemical properties of d5SICS or dNaM-properties that are central to the chemical integrity of cellular DNA. In this contribution, it is shown that excitation of d5SICS or dNaM with near-visible light leads to efficient trapping of population in the nucleoside's excited triplet state in high yield. Photoactivation of these long-lived, reactive states is shown to photosensitize cells, leading to the generation of reactive oxygen species and to a marked decrease in cell proliferation, thus warning scientists of the potential phototoxic side effects of expanding the genetic alphabet.