Structural basis of damage recognition by thymine DNA glycosylase: Key roles for N-terminal residues

Single-molecule analysis reveals TDG exhibits multiple modes of linear diffusion to process 5-formylcytosine

DNA methylation plays a key role in epigenetics, with 60-80% of CpG sites containing 5-methylcytosine. Base excision repair (BER) is suggested to be the main pathway involved in active DNA demethylation. 5-formylctyosine (5fC), an oxidized moiety of methylated cytosine, is recognized and removed by thymine DNA glycosylase (TDG) to generate an abasic site. TDG binds avidly to abasic sites and is product inhibited. Using single molecule fluorescence experiments, we saw TDG interact with DNA containing 5fC specifically and non-specifically with lifetimes of 72.9 and 7.5 seconds, respectively. These results indicate that TDG cleaves the 5fC and stays bound for an extended time at the generated abasic site. Mean squared displacement analysis and a two color TDG experiment indicate that TDG exhibits multiple modes of linear diffusion, including hopping and sliding, in search of a lesion. The catalytically crippled variants, N140A and R275A/L, have a reduced binding lifetime compared to wild type and Mean Squared Displacement (MSD) analysis indicates that R275L/A moves on the DNA with a faster diffusivity. These results indicate that mutating R275, but not N140 interferes with damage recognition by TDG. Our findings give insight into how TDG searches for its lesions in long stretches of undamaged DNA.

Backbone Conformational Equilibrium in Mismatched DNA Correlates with Enzyme Activity

T:G mismatches in mammals arise primarily from the deamination of methylated CpG sites or the incorporation of improper nucleotides. The process by which repair enzymes such as thymine DNA glycosylase (TDG) identify a canonical DNA base in the incorrect pairing context remains a mystery. However, the abundant contacts of the repair enzymes with the DNA backbone suggest a role for protein-phosphate interaction in the recognition and repair processes, where conformational properties may facilitate the proper interactions. We have previously used 31P NMR to investigate the energetics of DNA backbone BI-BII interconversion and the effect of a mismatch or lesion compared to canonical DNA and found stepwise differences in ΔG of 1-2 kcal/mol greater than equivalent steps in unmodified DNA. We have currently compared our results to substrate dependence for TDG, MBD4, M. HhaI, and CEBPβ, testing for correlations to sequence and base-pair dependence. We found strong correlations of our DNA phosphate backbone equilibrium (Keq) to different enzyme kinetics or binding parameters of these varied enzymes, suggesting that the backbone equilibrium may play an important role in mismatch recognition and/or conformational rearrangement and energetics during nucleotide flipping or other aspects of enzyme interrogation of the DNA substrate.

Chimeric d/l-DNA Probes of Base Excision Repair Enable Real-Time Monitoring of Thymine DNA Glycosylase Activity in Live Cells

The base excision repair (BER) pathway is a frontline defender of genomic integrity and plays a central role in epigenetic regulation through its involvement in the erasure of 5-methylcytosine. This biological and clinical significance has led to a demand for analytical methods capable of monitoring BER activities, especially in living cells. Unfortunately, prevailing methods, which are primarily derived from nucleic acids, are mostly incompatible with intracellular use due to their susceptibility to nuclease degradation and other off-target interactions. These limitations preclude important biological studies of BER enzymes and many clinical applications. Herein, we report a straightforward approach for constructing biostable BER probes using a unique chimeric d/l-DNA architecture that exploits the bioorthogonal properties of mirror-image l-DNA. We show that chimeric BER probes have excellent stability within living cells, where they were successfully employed to monitor relative BER activity, evaluate the efficiency of small molecule BER inhibitors, and study enzyme mutants. Notably, we report the first example of a fluorescent probe for real-time monitoring of thymine DNA glycosylase (TDG)-mediated BER of 5-formylcytosine and 5-carboxylcytosine in living cells, providing a much-needed tool for studying DNA (de)methylation biology. Chimeric probes offer a robust and highly generalizable approach for real-time monitoring of BER activity in living cells, which should enable a broad spectrum of basic research and clinical applications.

Characterization of a Novel Thermostable DNA Lyase Used To Prepare DNA for Next-Generation Sequencing

Recently, we constructed a hybrid thymine DNA glycosylase (hyTDG) by linking a 29-amino acid sequence from the human thymine DNA glycosylase with the catalytic domain of DNA mismatch glycosylase (MIG) from M. thermoautotrophicum, increasing the overall activity of the glycosylase. Previously, it was shown that a tyrosine to lysine (Y126K) mutation in the catalytic site of MIG could convert the glycosylase activity to a lyase activity. We made the corresponding mutation to our hyTDG to create a hyTDG-lyase (Y163K). Here, we report that the hybrid mutant has robust lyase activity, has activity over a broad temperature range, and is active under multiple buffer conditions. The hyTDG-lyase cleaves an abasic site similar to endonuclease III (Endo III). In the presence of β-mercaptoethanol (β-ME), the abasic site unsaturated aldehyde forms a β-ME adduct. The hyTDG-lyase maintains its preference for cleaving opposite G, as with the hyTDG glycosylase, and the hyTDG-lyase and hyTDG glycosylase can function in tandem to cleave T:G mismatches. The hyTDG-lyase described here should be a valuable tool in studies examining DNA damage and repair. Future studies will utilize these enzymes to quantify T:G mispairs in cells, tissues, and genomic DNA using next-generation sequencing.

Noncatalytic Domains in DNA Glycosylases

Many proteins consist of two or more structural domains: separate parts that have a defined structure and function. For example, in enzymes, the catalytic activity is often localized in a core fragment, while other domains or disordered parts of the same protein participate in a number of regulatory processes. This situation is often observed in many DNA glycosylases, the proteins that remove damaged nucleobases thus initiating base excision DNA repair. This review covers the present knowledge about the functions and evolution of such noncatalytic parts in DNA glycosylases, mostly concerned with the human enzymes but also considering some unique members of this group coming from plants and prokaryotes.

1H NMR Chemical Exchange Techniques Reveal Local and Global Effects of Oxidized Cytosine Derivatives

5-Carboxycytosine (5caC) is a rare epigenetic modification found in nucleic acids of all domains of life. Despite its sparse genomic abundance, 5caC is presumed to play essential regulatory roles in transcription, maintenance and base-excision processes in DNA. In this work, we utilize nuclear magnetic resonance (NMR) spectroscopy to address the effects of 5caC incorporation into canonical DNA strands at multiple pH and temperature conditions. Our results demonstrate that 5caC has a pH-dependent global destabilizing and a base-pair mobility enhancing local impact on dsDNA, albeit without any detectable influence on the ground-state B-DNA structure. Measurement of hybridization thermodynamics and kinetics of 5caC-bearing DNA duplexes highlighted how acidic environment (pH 5.8 and 4.7) destabilizes the double-stranded structure by ∼10–20 kJ mol^–1 at 37 °C when compared to the same sample at neutral pH. Protonation of 5caC results in a lower activation energy for the dissociation process and a higher barrier for annealing. Studies on conformational exchange on the microsecond time scale regime revealed a sharply localized base-pair motion involving exclusively the modified site and its immediate surroundings. By direct comparison with canonical and 5-formylcytosine (5fC)-edited strands, we were able to address the impact of the two most oxidized naturally occurring cytosine derivatives in the genome. These insights on 5caC’s subtle sensitivity to acidic pH contribute to the long-standing questions of its capacity as a substrate in base excision repair processes and its purpose as an independent, stable epigenetic mark.

Computational investigations on target-site searching and recognition mechanisms by thymine DNA glycosylase during DNA repair process

DNA glycosylase, as one member of DNA repair machineries, plays an essential role in correcting mismatched/damaged DNA nucleotides by cleaving the N-glycosidic bond between the sugar and target nucleobase through the base excision repair (BER) pathways. Efficient corrections of these DNA lesions are critical for maintaining genome integrity and preventing premature aging and cancers. The target-site searching/recognition mechanisms and the subsequent conformational dynamics of DNA glycosylase, however, remain challenging to be characterized using experimental techniques. In this review, we summarize our recent studies of sequential structural changes of thymine DNA glycosylase (TDG) during the DNA repair process, achieved mostly by molecular dynamics (MD) simulations. Computational simulations allow us to reveal atomic-level structural dynamics of TDG as it approaches the target-site, and pinpoint the key structural elements responsible for regulating the translocation of TDG along DNA. Subsequently, upon locating the lesions, TDG adopts a base-flipping mechanism to extrude the mispaired nucleobase into the enzyme active-site. The constructed kinetic network model elucidates six metastable states during the base-extrusion process and suggests an active role of TDG in flipping the intrahelical nucleobase. Finally, the molecular mechanism of product release dynamics after catalysis is also summarized. Taken together, we highlight to what extent the computational simulations advance our knowledge and understanding of the molecular mechanism underlying the conformational dynamics of TDG, as well as the limitations of current theoretical work.

Kinetic Analysis of the Effect of N-Terminal Acetylation on Thymine DNA Glycosylase

Thymine DNA glycosylase (TDG) is tasked with initiating DNA base excision repair by recognizing and removing T, U, the chemotherapeutic 5-fluorouracil (5-FU), and many other oxidized and halogenated pyrimidine bases. TDG contains a long, unstructured N-terminus that contains four known sites of acetylation: lysine (K) residues 59, 83, 84, and 87. Here, K to glutamine (Q) mutants are used as acetyl-lysine (AcK) analogues to probe the effect of N-terminal acetylation on the kinetics of TDG. We find that mimicking acetylation affects neither the maximal single-turnover rate k_max nor the turnover rate k_TO, indicating that the steps after initial binding, through chemistry and product release, are not affected. Under subsaturating conditions, however, acetylation changes the processing of U substrates. Subtle differences among AcK analogues are revealed with 5-FU in single-stranded DNA. We propose that the subtleties observed among the AcK analogues may be amplified on the genomic scale, leading to regulation of TDG activity. N-terminal acetylation, though, may also play a structural, rather than kinetic role in vivo.

A systematic genome-wide account of binding sites for the model transcription factor Gcn4

Sequence-specific DNA-binding transcription factors are central to gene regulation. They are often associated with consensus binding sites that predict far more genomic sites than are bound in vivo. One explanation is that most sites are blocked by nucleosomes, such that only sites in nucleosome-depleted regulatory regions are bound. We compared the binding of the yeast transcription factor Gcn4 in vivo using published ChIP-seq data (546 sites) and in vitro, using a modified SELEX method (“G-SELEX”), which utilizes short genomic DNA fragments to quantify binding at all sites. We confirm that Gcn4 binds strongly to an AP-1-like sequence (TGACTCA) and weakly to half-sites. However, Gcn4 binds only some of the 1078 exact matches to this sequence, even in vitro. We show that there are only 166 copies of the high-affinity RTGACTCAY site (exact match) in the yeast genome, all occupied in vivo, largely independently of whether they are located in nucleosome-depleted or nucleosomal regions. Generally, RTGACTCAR/YTGACTCAY sites are bound much more weakly and YTGACTCAR sites are unbound, with biological implications for determining induction levels. We conclude that, to a first approximation, Gcn4 binding can be predicted using the high-affinity site, without reference to chromatin structure. We propose that transcription factor binding sites should be defined more precisely using quantitative data, allowing more accurate genome-wide prediction of binding sites and greater insight into gene regulation.

Interaction of Thymine DNA Glycosylase with Oxidised 5-Methyl-cytosines in Their Amino- and Imino-Forms

Thymine DNA Glycosylase (TDG) is an enzyme of the base excision repair mechanism and removes damaged or mispaired bases from DNA via hydrolysis of the glycosidic bond. Specificity is of high importance for such a glycosylase, so as to avoid the damage of intact DNA. Among the substrates reported for TDG are mispaired uracil and thymine but also formyl-cytosine and carboxyl-cytosine. Methyl-cytosine and hydroxylmethyl-cytosine are, in contrast, not processed by the TDG enzyme. We have in this work employed molecular dynamics simulations to explore the conformational dynamics of DNA carrying a formyl-cytosine or carboxyl-cytosine and compared those to DNA with the non-cognate bases methyl-cytosine and hydroxylmethyl-cytosine, as amino and imino tautomers. Whereas for the mispairs a wobble conformation is likely decisive for recognition, all amino tautomers of formyl-cytosine and carboxyl-cytosine exhibit the same Watson–Crick conformation, but all imino tautomers indeed form wobble pairs. The conformational dynamics of the amino tautomers in free DNA do not exhibit differences that could be exploited for recognition, and also complexation to the TDG enzyme does not induce any alteration that would indicate preferable binding to one or the other oxidised methyl-cytosine. The imino tautomers, in contrast, undergo a shift in the equilibrium between a closed and a more open, partially flipped state, towards the more open form upon complexation to the TDG enzyme. This stabilisation of the more open conformation is most pronounced for the non-cognate bases methyl-cytosine and hydroxyl-cytosine and is thus not a likely mode for recognition. Moreover, calculated binding affinities for the different forms indicate the imino forms to be less likely in the complexed DNA. These findings, together with the low probability of imino tautomers in free DNA and the indifference of the complexed amino tautomers, suggest that discrimination of the oxidised methyl-cytosines does not take place in the initial complex formation.

Structural Insights into the Mechanism of Base Excision by MBD4

DNA glycosylases remove damaged or modified nucleobases by cleaving the N-glycosyl bond and the correct nucleotide is restored through subsequent base excision repair. In addition to excising threatening lesions, DNA glycosylases contribute to epigenetic regulation by mediating DNA demethylation and perform other important functions. However, the catalytic mechanism remains poorly defined for many glycosylases, including MBD4 (methyl-CpG binding domain IV), a member of the helix-hairpin-helix (HhH) superfamily. MBD4 excises thymine from G·T mispairs, suppressing mutations caused by deamination of 5-methylcytosine, and it removes uracil and modified uracils (e.g., 5-hydroxymethyluracil) mispaired with guanine. To investigate the mechanism of MBD4 we solved high-resolution structures of enzyme-DNA complexes at three stages of catalysis. Using a non-cleavable substrate analog, 2'-deoxy-pseudouridine, we determined the first structure of an enzyme-substrate complex for wild-type MBD4, which confirms interactions that mediate lesion recognition and suggests that a catalytic Asp, highly conserved in HhH enzymes, binds the putative nucleophilic water molecule and stabilizes the transition state. Observation that mutating the Asp (to Gly) reduces activity by 2700-fold indicates an important role in catalysis, but probably not one as the nucleophile in a double-displacement reaction, as previously suggested. Consistent with direct-displacement hydrolysis, a structure of the enzyme-product complex indicates a reaction leading to inversion of configuration. A structure with DNA containing 1-azadeoxyribose models a potential oxacarbenium-ion intermediate and suggests the Asp could facilitate migration of the electrophile towards the nucleophilic water. Finally, the structures provide detailed snapshots of the HhH motif, informing how these ubiquitous metal-binding elements mediate DNA binding.

Gas-Phase Experimental and Computational Studies of 5-Halouracils: Intrinsic Properties and Biological Implications

The gas-phase acidity and proton affinity (PA) of 5-halouracils (5-fluorouracil, 5-chlorouracil, 5-bromouracil, and 5-iodouracil) have been examined using both theoretical and experimental methods. This work represents a comprehensive study of the thermochemical properties of these nucleobases. Other than 5-fluorouracil acidity, the intrinsic acidity and PA of these halouracils have not been heretofore measured; these new experimental data provide a benchmark for the computational values. Furthermore, we examine these 5-halouracils in the context of the enzyme thymine DNA glycosylase (TDG), which is an enzyme that protects the genome by cleaving these substrates from DNA. Our gas-phase results are compared and contrasted to TDG excision rates to afford insights into the TDG mechanism.

Reversible chromatin condensation by the DNA repair and demethylation factor thymine DNA glycosylase

Chromatin structures (and modulators thereof) play a central role in genome organization and function. Herein, we report that thymine DNA glycosylase (TDG), an essential enzyme involved in DNA repair and demethylation, has the capacity to alter chromatin structure directly through its physical interactions with DNA. Using chemically defined nucleosome arrays, we demonstrate that TDG induces decompaction of individual chromatin fibers upon binding and promotes self-association of nucleosome arrays into higher-order oligomeric structures (i.e. condensation). Chromatin condensation is mediated by TDG’s disordered polycationic N-terminal domain, whereas its C-terminal domain antagonizes this process. Furthermore, we demonstrate that TDG-mediated chromatin condensation is reversible by growth arrest and DNA damage 45 alpha (GADD45a), implying that TDG cooperates with its binding partners to dynamically control chromatin architecture. Finally, we show that chromatin condensation by TDG is sensitive to the methylation status of the underlying DNA. This new paradigm for TDG has specific implications for associated processes, such as DNA repair, DNA demethylation, and transcription, and general implications for the role of DNA modification ‘readers’ in controlling chromatin organization.

Single-Base Lesions and Mismatches Alter the Backbone Conformational Dynamics in DNA

DNA damage has been implicated in numerous human diseases, particularly cancer, and the aging process. Single-base lesions and mismatches in DNA can be cytotoxic or mutagenic and are recognized by a DNA glycosylase during the process of base excision repair. Altered local dynamics and conformational properties in damaged DNAs have previously been suggested to assist in recognition and specificity. Herein, we use solution nuclear magnetic resonance to quantify changes in BI-BII backbone conformational dynamics due to the presence of single-base lesions in DNA, including uracil, dihydrouracil, 1,N6-ethenoadenine, and T:G mismatches. Stepwise changes to the %BII and ΔG of the BI-BII dynamic equilibrium compared to those of unmodified sequences were observed. Additionally, the equilibrium skews toward endothermicity for the phosphates nearest the lesion/mismatched base pair. Finally, the phosphates with the greatest alterations correlate with those most relevant to the repair of enzyme binding. All of these results suggest local conformational rearrangement of the DNA backbone may play a role in lesion recognition by repair enzymes.

Atomic resolution of short-range sliding dynamics of thymine DNA glycosylase along DNA minor-groove for lesion recognition

Thymine DNA glycosylase (TDG), as a repair enzyme, plays essential roles in maintaining the genome integrity by correcting several mismatched/damaged nucleobases. TDG acquires an efficient strategy to search for the lesions among a vast number of cognate base pairs. Currently, atomic-level details of how TDG translocates along DNA as it approaches the lesion site and the molecular mechanisms of the interplay between TDG and DNA are still elusive. Here, by constructing the Markov state model based on hundreds of molecular dynamics simulations with an integrated simulation time of ∼25 μs, we reveal the rotation-coupled sliding dynamics of TDG along a 9 bp DNA segment containing one G·T mispair. We find that TDG translocates along DNA at a relatively faster rate when distant from the lesion site, but slows down as it approaches the target, accompanied by deeply penetrating into the minor-groove, opening up the mismatched base pair and significantly sculpturing the DNA shape. Moreover, the electrostatic interactions between TDG and DNA are found to be critical for mediating the TDG translocation. Notably, several uncharacterized TDG residues are identified to take part in regulating the conformational switches of TDG occurred in the site-transfer process, which warrants further experimental validations.

TET-TDG Active DNA Demethylation at CpG and Non-CpG Sites

In mammalian genomes, cytosine methylation occurs predominantly at CG (or CpG) dinucleotide contexts. As part of dynamic epigenetic regulation, 5-methylcytosine (mC) can be erased by active DNA demethylation, whereby ten-eleven translocation (TET) enzymes catalyze the stepwise oxidation of mC to 5-hydroxymethylcytosine (hmC), 5-formylcytosine (fC), and 5-carboxycytosine (caC), thymine DNA glycosylase (TDG) excises fC or caC, and base excision repair yields unmodified cytosine. In certain cell types, mC is also enriched at some non-CG (or CH) dinucleotides, however hmC is not. To provide biochemical context for the distribution of modified cytosines observed in biological systems, we systematically analyzed the activity of human TET2 and TDG for substrates in CG and CH contexts. We find that while TET2 oxidizes mC more efficiently in CG versus CH sites, this context preference can be diminished for hmC oxidation. Remarkably, TDG excision of fC and caC is only modestly dependent on CG context, contrasting its strong context dependence for thymine excision. We show that collaborative TET-TDG oxidation-excision activity is only marginally reduced for CA versus CG contexts. Our findings demonstrate that the TET-TDG-mediated demethylation pathway is not limited to CG sites and suggest a rationale for the depletion of hmCH in genomes rich in mCH.

Impact of 5-formylcytosine on the melting kinetics of DNA by 1H NMR chemical exchange

5-Formylcytosine (5fC) is a chemically edited, naturally occurring nucleobase which appears in the context of modified DNA strands. The understanding of the impact of 5fC on dsDNA physical properties is to date limited. In this work, we applied temperature-dependent 1H Chemical Exchange Saturation Transfer (CEST) NMR experiments to non-invasively and site-specifically measure the thermodynamic and kinetic influence of formylated cytosine nucleobase on the melting process involving dsDNA. Incorporation of 5fC within symmetrically positioned CpG sites destabilizes the whole dsDNA structure-as witnessed from the ∼2°C decrease in the melting temperature and 5-10 kJ mol-1 decrease in ΔG°-and affects the kinetic rates of association and dissociation. We observed an up to ∼5-fold enhancement of the dsDNA dissociation and an up to ∼3-fold reduction in ssDNA association rate constants, over multiple temperatures and for several proton reporters. Eyring and van't Hoff analysis proved that the destabilization is not localized, instead all base-pairs are affected and the transition states resembles the single-stranded conformation. These results advance our knowledge about the role of 5fC as a semi-permanent epigenetic modification and assist in the understanding of its interactions with reader proteins.

Excision of 5-Carboxylcytosine by Thymine DNA Glycosylase

5-Methylcytosine (mC) is an epigenetic mark that is written by methyltransferases, erased through passive and active mechanisms, and impacts transcription, development, diseases including cancer, and aging. Active DNA demethylation involves TET-mediated stepwise oxidation of mC to 5-hydroxymethylcytosine, 5-formylcytosine (fC), or 5-carboxylcytosine (caC), excision of fC or caC by thymine DNA glycosylase (TDG), and subsequent base excision repair. Many elements of this essential process are poorly defined, including TDG excision of caC. To address this problem, we solved high-resolution structures of human TDG bound to DNA with cadC (5-carboxyl-2'-deoxycytidine) flipped into its active site. The structures unveil detailed enzyme-substrate interactions that mediate recognition and removal of caC, many involving water molecules. Importantly, two water molecules contact a carboxylate oxygen of caC and are poised to facilitate acid-catalyzed caC excision. Moreover, a substrate-dependent conformational change in TDG modulates the hydrogen bond interactions for one of these waters, enabling productive interaction with caC. An Asn residue (N191) that is critical for caC excision is found to contact N3 and N4 of caC, suggesting a mechanism for acid-catalyzed base excision that features an N3-protonated form of caC but would be ineffective for C, mC, or hmC. We also investigated another Asn residue (N140) that is catalytically essential and strictly conserved in the TDG-MUG enzyme family. A structure of N140A-TDG bound to cadC DNA provides the first high-resolution insight into how enzyme-substrate interactions, including water molecules, are impacted by depleting the conserved Asn, informing its role in binding and addition of the nucleophilic water molecule.

Defining the impact of sumoylation on substrate binding and catalysis by thymine DNA glycosylase

Thymine DNA glycosylase (TDG) excises thymine from mutagenic G·T mispairs generated by deamination of 5-methylcytosine (mC) and it removes two mC derivatives, 5−formylcytosine (fC) and 5−carboxylcytosine (caC), in a multistep pathway for DNA demethylation. TDG is modified by small ubiquitin-like modifier (SUMO) proteins, but the impact of sumoylation on TDG activity is poorly defined and the functions of TDG sumoylation remain unclear. We determined the effect of TDG sumoylation, by SUMO-1 or SUMO-2, on substrate binding and catalytic parameters. Single turnover experiments reveal that sumoylation dramatically impairs TDG base-excision activity, such that G·T activity is reduced by ≥45-fold and fC and caC are excised slowly, with a reaction half-life of ≥9 min (37°C). Fluorescence anisotropy studies reveal that unmodified TDG binds tightly to G·fC and G·caC substrates, with dissociation constants in the low nanomolar range. While sumoylation of TDG weakens substrate binding, the residual affinity is substantial and is comparable to that of biochemically-characterized readers of fC and caC. Our findings raise the possibility that sumoylation enables TDG to function, at least transiently, as reader of fC and caC. Notably, sumoylation could potentially facilitate TDG recruitment of other proteins, including transcription factors or epigenetic regulators, to these sites in DNA.

N-terminal domain of human uracil DNA glycosylase (hUNG2) promotes targeting to uracil sites adjacent to ssDNA-dsDNA junctions

The N-terminal domain (NTD) of nuclear human uracil DNA glycosylase (hUNG2) assists in targeting hUNG2 to replication forks through specific interactions with replication protein A (RPA). Here, we explored hUNG2 activity in the presence and absence of RPA using substrates with ssDNA–dsDNA junctions that mimic structural features of the replication fork and transcriptional R-loops. We find that when RPA is tightly bound to the ssDNA overhang of junction DNA substrates, base excision by hUNG2 is strongly biased toward uracils located 21 bp or less from the ssDNA–dsDNA junction. In the absence of RPA, hUNG2 still showed an 8-fold excision bias for uracil located <10 bp from the junction, but only when the overhang had a 5′ end. Biased targeting required the NTD and was not observed with the hUNG2 catalytic domain alone. Consistent with this requirement, the isolated NTD was found to bind weakly to ssDNA. These findings indicate that the NTD of hUNG2 targets the enzyme to ssDNA–dsDNA junctions using RPA-dependent and RPA-independent mechanisms. This structure-based specificity may promote efficient removal of uracils that arise from dUTP incorporation during DNA replication, or additionally, uracils that arise from DNA cytidine deamination at transcriptional R-loops during immunoglobulin class-switch recombination.

Computational and Experimental Druggability Assessment of Human DNA Glycosylases

Due to a polar or even charged binding interface, DNA-binding proteins are considered extraordinarily difficult targets for development of small-molecule ligands and only a handful of proteins have been targeted successfully to date. Recently, however, it has been shown that development of selective and efficient inhibitors of 8-oxoguanine DNA glycosylase is possible. Here, we describe the initial druggability assessment of DNA glycosylases in a computational setting and experimentally investigate several methods to target endonuclease VIII-like 1 (NEIL1) with small-molecule inhibitors. We find that DNA glycosylases exhibit good predicted druggability in both DNA-bound and -unbound states. Furthermore, we find catalytic sites to be highly flexible, allowing for a range of interactions and binding partners. One flexible catalytic site was rationalized for NEIL1 and further investigated experimentally using both a biochemical assay in the presence of DNA and a thermal shift assay in the absence of DNA.

Base-flipping dynamics from an intrahelical to an extrahelical state exerted by thymine DNA glycosylase during DNA repair process

Thymine DNA glycosylase (TDG) is a DNA repair enzyme that excises a variety of mismatched or damaged nucleotides (nts), e.g. dU, dT, 5fC and 5caC. TDG is shown to play essential roles in maintaining genome integrity and correctly programming epigenetic modifications through DNA demethylation. After locating the lesions, TDG employs a base-flipping strategy to recognize the damaged nucleobases, whereby the interrogated nt is extruded from the DNA helical stack and binds into the TDG active site. The dynamic mechanism of the base-flipping process at an atomistic resolution, however, remains elusive. Here, we employ the Markov State Model (MSM) constructed from extensive all-atom molecular dynamics (MD) simulations to reveal the complete base-flipping process for a G.T mispair at a tens of microsecond timescale. Our studies identify critical intermediates of the mispaired dT during its extrusion process and reveal the key TDG residues involved in the inter-state transitions. Notably, we find an active role of TDG in promoting the intrahelical nt eversion, sculpturing the DNA backbone, and penetrating into the DNA minor groove. Three additional TDG substrates, namely dU, 5fC, and 5caC, are further tested to evaluate the substituent effects of various chemical modifications of the pyrimidine ring on base-flipping dynamics.

Dynamics of the excised base release in thymine DNA glycosylase during DNA repair process

Thymine DNA glycosylase (TDG) initiates base excision repair by cleaving the N-glycosidic bond between the sugar and target base. After catalysis, the release of excised base is a requisite step to terminate the catalytic cycle and liberate the TDG for the following enzymatic reactions. However, an atomistic-level understanding of the dynamics of the product release process in TDG remains unknown. Here, by employing molecular dynamics simulations combined with the Markov State Model, we reveal the dynamics of the thymine release after the excision at microseconds timescale and all-atom resolution. We identify several key metastable states of the thymine and its dominant releasing pathway. Notably, after replacing the TDG residue Gly142 with tyrosine, the thymine release is delayed compared to the wild-type (wt) TDG, as supported by our potential of mean force (PMF) calculations. These findings warrant further experimental tests to potentially trap the excised base in the active site of TDG after the catalysis, which had been unsuccessful by previous attempts. Finally, we extended our studies to other TDG products, including the uracil, 5hmU, 5fC and 5caC bases in order to compare the product release for different targeting bases in the TDG–DNA complex.

Reaction Mechanism for the N-Glycosidic Bond Cleavage of 5-Formylcytosine by Thymine DNA Glycosylase

Thymine DNA glycosylase (TDG) initiates the base excision repair mechanism for the deamination and oxidation products of cytosine and 5-methylcytosine. This enzyme has a key role in epigenetic regulation, and its catalytic inactivation results in, e.g., mice embryo lethality. Here, we employ molecular dynamics simulations and quantum mechanics/molecular mechanics calculations to investigate the reaction mechanism of the TDG-catalyzed N-glycosidic bond hydrolysis of the modified base 5-formylcytosine. Our results reveal a reaction pathway, which in its first step features a reorganization of the substrate that lowers the barrier height for the subsequent C1'-N1 bond dissociation. The suggested mechanism is consistent with the experimental data, as it is not acid-catalyzed and proceeds through an oxocarbenium-like transition state. It also provides insights into the catalytic roles of the Thr197 and Asn140 residues.

Defining the Role of Nucleotide Flipping in Enzyme Specificity Using 19F NMR

A broad range of proteins employ nucleotide flipping to recognize specific sites in nucleic acids, including DNA glycosylases, which remove modified nucleobases to initiate base excision repair. Deamination, a pervasive mode of damage, typically generates lesions that are recognized by glycosylases as being foreign to DNA. However, deamination of 5-methylcytosine (mC) generates thymine, a canonical DNA base, presenting a challenge for damage recognition. Nevertheless, repair of mC deamination is important because the resulting G·T mispairs cause C → T transition mutations, and mC is abundant in all three domains of life. Countering this threat are three types of glycosylases that excise thymine from G·T mispairs, including thymine DNA glycosylase (TDG). These enzymes must minimize excision of thymine that is not generated by mC deamination, in A·T pairs and in polymerase-generated G·T mispairs. TDG preferentially removes thymine from DNA contexts in which cytosine methylation is prevalent, including CG and one non-CG site. This remarkable context specificity could be attained through modulation of nucleotide flipping, a reversible step that precedes base excision. We tested this idea using fluorine NMR and DNA containing 2'-fluoro-substituted nucleotides. We find that dT nucleotide flipping depends on DNA context and is efficient only in contexts known to feature cytosine methylation. We also show that a conserved Ala residue limits thymine excision by hindering nucleotide flipping. A linear free energy correlation reveals that TDG attains context specificity for thymine excision through modulation of nucleotide flipping. Our results provide a framework for characterizing nucleotide flipping in nucleic acids using ¹⁹F NMR.

Facilitated Diffusion Mechanisms in DNA Base Excision Repair and Transcriptional Activation

Preservation of the coding potential of the genome and highly regulated gene expression over the life span of a human are two fundamental requirements of life. These processes require the action of repair enzymes or transcription factors that efficiently recognize specific sites of DNA damage or transcriptional regulation within a restricted time frame of the cell cycle or metabolism. A failure of these systems to act results in accumulated mutations, metabolic dysfunction, and disease. Despite the multifactorial complexity of cellular DNA repair and transcriptional regulation, both processes share a fundamental physical requirement that the proteins must rapidly diffuse to their specific DNA-binding sites that are embedded within the context of a vastly greater number of nonspecific DNA-binding sites. Superimposed on the needle-in-the-haystack problem is the complex nature of the cellular environment, which contains such high concentrations of macromolecules that the time frame for diffusion is expected to be severely extended as compared to dilute solution. Here we critically review the mechanisms for how these proteins solve the needle-in-the-haystack problem and how the effects of cellular macromolecular crowding can enhance facilitated diffusion processes. We restrict the review to human proteins that use stochastic, thermally driven site-recognition mechanisms, and we specifically exclude systems involving energy cofactors or circular DNA clamps. Our scope includes ensemble and single-molecule studies of the past decade or so, with an emphasis on connecting experimental observations to biological function.

Interactions of the DNA Repair Enzyme Human Thymine DNA Glycosylase with Cognate and Noncognate DNA

Glycosylases specifically recognize and flip their target base out of the DNA helix into the enzyme's active site. Our simulations show that a partially flipped state, already present in free DNA carrying a T:G mispair, becomes the more probable state compared to the closed state after binding of thymine DNA glycosylase (TDG). Paired thymine (T:A) or methyl-cytosine (mC:G) does not exhibit a partially flipped state in free or complexed DNA. Important enzyme-DNA interactions exhibit significant strength in the intrahelical and extrahelical TDG-DNA complexes. The computed binding free energy differences suggest these interactions account for the stabilization of the partially flipped state, thereby driving the T:G mispair toward base flip. In the fully flipped state, the cognate base thymine is significantly better accommodated in the enzyme's active site than noncognate bases are, suggesting the hydrolysis step as the last of several stages at which base recognition can be achieved.

Nucleosomes and the three glycosylases: High, medium, and low levels of excision by the uracil DNA glycosylase superfamily

Human cells express the UDG superfamily of glycosylases, which excise uracil (U) from the genome. The three members of this structural superfamily are uracil DNA glycosylase (UNG/UDG), single-strand selective monofunctional uracil DNA glycosylase (SMUG1), and thymine DNA glycosylase (TDG). We previously reported that UDG is efficient at removing U from DNA packaged into nucleosome core particles (NCP) and is minimally affected by the histone proteins when acting on an outward-facing U in the dyad region. In an effort to determine whether this high activity is a general property of the UDG superfamily of glycosylases, we compare the activity of UDG, SMUG1, and TDG on a U:G wobble base pair using NCP assembled from Xenopus laevis histones and the Widom 601 positioning sequence. We found that while UDG is highly active, SMUG1 is severely inhibited on NCP and this inhibition is independent of sequence context. Here we also provide the first report of TDG activity on an NCP, and found that TDG has an intermediate level of activity in excision of U and is severely inhibited in its excision of T. These results are discussed in the context of cellular roles for each of these enzymes.

Uncovering universal rules governing the selectivity of the archetypal DNA glycosylase TDG

Thymine DNA glycosylase (TDG) is a pivotal enzyme with dual roles in both genome maintenance and epigenetic regulation. TDG is involved in cytosine demethylation at CpG sites in DNA. Here we have used molecular modeling to delineate the lesion search and DNA base interrogation mechanisms of TDG. First, we examined the capacity of TDG to interrogate not only DNA substrates with 5-carboxyl cytosine modifications but also G:T mismatches and nonmismatched (A:T) base pairs using classical and accelerated molecular dynamics. To determine the kinetics, we constructed Markov state models. Base interrogation was found to be highly stochastic and proceeded through insertion of an arginine-containing loop into the DNA minor groove to transiently disrupt Watson-Crick pairing. Next, we employed chain-of-replicas path-sampling methodologies to compute minimum free energy paths for TDG base extrusion. We identified the key intermediates imparting selectivity and determined effective free energy profiles for the lesion search and base extrusion into the TDG active site. Our results show that DNA sculpting, dynamic glycosylase interactions, and stabilizing contacts collectively provide a powerful mechanism for the detection and discrimination of modified bases and epigenetic marks in DNA.

A low-complexity region in the YTH domain protein Mmi1 enhances RNA binding

Mmi1 is an essential RNA-binding protein in the fission yeast Schizosaccharomyces pombe that eliminates meiotic transcripts during normal vegetative growth. Mmi1 contains a YTH domain that binds specific RNA sequences, targeting mRNAs for degradation. The YTH domain of Mmi1 uses a noncanonical RNA-binding surface that includes contacts outside the conserved fold. Here, we report that an N-terminal extension that is proximal to the YTH domain enhances RNA binding. Using X-ray crystallography, NMR, and biophysical methods, we show that this low-complexity region becomes more ordered upon RNA binding. This enhances the affinity of the interaction of the Mmi1 YTH domain with specific RNAs by reducing the dissociation rate of the Mmi1-RNA complex. We propose that the low-complexity region influences RNA binding indirectly by reducing dynamic motions of the RNA-binding groove and stabilizing a conformation of the YTH domain that binds to RNA with high affinity. Taken together, our work reveals how a low-complexity region proximal to a conserved folded domain can adopt an ordered structure to aid nucleic acid binding.

Nei-like 1 (NEIL1) excises 5-carboxylcytosine directly and stimulates TDG-mediated 5-formyl and 5-carboxylcytosine excision

Thymine DNA glycosylase (TDG) and Nei-like 1 (NEIL1) have both been implicated in the base excision repair step of active DNA demethylation. The robust glycosylase activity of TDG on DNA substrates containing 5-formylcytosine (5fC) or 5-carboxylcytosine (5caC) is universally accepted, but the mode of action of NEIL1 is still debated. Based on genetic experiments, it has been suggested that NEIL1 acts redundantly with TDG and excises 5fC and 5caC directly. However, this result has been disputed, and it was suggested instead that NEIL1 is recruited by the monofunctional TDG for the 2′-deoxyribose excision step. Using purified human NEIL1 and its catalytically impaired P2T and E3Q variants as controls, we detect NEIL1 activity on 5caC, but not a 5fC containing dsDNA substrate. We confirm direct NEIL1 TDG binding and NEIL1 mediated 2′-deoxyribose excision downstream of TDG glycosylase activity. NEIL1 acts not only downstream of TDG, but also enhances TDG activity on 5fC or 5caC containing DNA. NEIL1 mediated enhancement of the TDG glycosylase activity is substrate specific and does not occur for dsDNA with a T/G mismatch.

Base excision repair of oxidative DNA damage: from mechanism to disease

Reactive oxygen species continuously assault the structure of DNA resulting in oxidation and fragmentation of the nucleobases. Both oxidative DNA damage itself and its repair mediate the progression of many prevalent human maladies. The major pathway tasked with removal of oxidative DNA damage, and hence maintaining genomic integrity, is base excision repair (BER). The aphorism that structure often dictates function has proven true, as numerous recent structural biology studies have aided in clarifying the molecular mechanisms used by key BER enzymes during the repair of damaged DNA. This review focuses on the mechanistic details of the individual BER enzymes and the association of these enzymes during the development and progression of human diseases, including cancer and neurological diseases. Expanding on these structural and biochemical studies to further clarify still elusive BER mechanisms, and focusing our efforts toward gaining an improved appreciation of how these enzymes form co-complexes to facilitate DNA repair is a crucial next step toward understanding how BER contributes to human maladies and how it can be manipulated to alter patient outcomes.

Structural Basis for Excision of 5-Formylcytosine by Thymine DNA Glycosylase

Thymine DNA glycosylase (TDG) is a base excision repair enzyme with key functions in epigenetic regulation. Performing a critical step in a pathway for active DNA demethylation, TDG removes 5-formylcytosine and 5-carboxylcytosine, oxidized derivatives of 5-methylcytosine that are generated by TET (ten-eleven translocation) enzymes. We determined a crystal structure of TDG bound to DNA with a noncleavable (2'-fluoroarabino) analogue of 5-formyldeoxycytidine flipped into its active site, revealing how it recognizes and hydrolytically excises fC. Together with previous structural and biochemical findings, the results illustrate how TDG employs an adaptable active site to excise a broad variety of nucleobases from DNA.