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

Quantitative and systematic NMR measurements of sequence-dependent A-T Hoogsteen dynamics uncovers unique conformational specificity in the DNA double helix

The propensities to form lowly-populated short-lived conformations of DNA could vary with sequence, providing an important source of sequence-specificity in biochemical reactions. However, comprehensively measuring how these dynamics vary with sequence is challenging. Using 1H CEST and 13C R 1 ρ NMR, we measured Watson-Crick to Hoogsteen dynamics for an A-T base pair in thirteen trinucleotide sequence contexts. The Hoogsteen population and exchange rate varied 4-fold and 16-fold, respectively, and were dependent on both the 3'- and 5'-neighbors but only weakly dependent on monovalent ion concentration (25 versus 100 mM NaCl) and pH (6.8 versus 8.0). Flexible TA and CA dinucleotide steps exhibited the highest Hoogsteen populations, and their kinetics rates strongly depended on the 3'-neighbor. In contrast, the stiffer AA and GA steps had the lowest Hoogsteen population, and their kinetics were weakly dependent on the 3'-neighbor. The Hoogsteen lifetime was especially short when G-C neighbors flanked the A-T base pair. The Hoogsteen dynamics had a distinct sequence-dependence compared to duplex stability and minor groove width. Thus, our results uncover a unique source of sequence-specificity hidden within the DNA double helix in the form of A-T Hoogsteen dynamics and establish the utility of 1H CEST to quantitively measure sequence-dependent DNA dynamics.

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.

All-or-None Selectivity in Probing Polarity-Determined Trinucleotide Repeat Foldings with a Parity Resolution by a Beyond-Size-Matching Ligand

Abnormal amplification of trinucleotide repeats (TNRs) is associated with neurodegenerative diseases by forming a particular hairpin bulge. It is well known that the polarity and parity of TNRs can regulate the formed hairpin structures. Therefore, there is a great challenge to efficiently discriminate the hairpin structures of TNRs with substantial selectivity. Herein, we developed a fluorescent ligand of pseudohypericin (Pse) with a beyond-size-matching (BSM) geometry to selectively sense hairpin structures of GTC and CTG TNRs. The GTC hairpin structures can bind with Pse dominantly at extreme T-T mismatches by the virtue of their most extrahelical conformations, while there is no binding event to occur with the polarity-inverted counterpart CTG hairpin structures because of the limited space provided by their intrahelical T-T mismatches. In addition, this all-or-none response with the polarity-dependent folding (PoDF) is independent of the length of these TNRs. Interestingly, the parity-dependent folding (PaDF) of GTC hairpin structures can also be resolved. Besides pure TNRs, the competency of this BSM ligand to sense the PoDF and PaDF effects was also generalized to DNAs with TNRs occurring at loop and stem end regions. To our knowledge, this is the first experimental observation with the state-of-the-art performance over the fluorescence measurement of PoDF and PaDF in TNRs. Our work provides an expedient way to elucidate the TNR folding by designing ligands having BSM features.

Rapid excision of oxidized adenine by human thymine DNA glycosylase

Oxidation of DNA bases generates mutagenic and cytotoxic lesions that are implicated in cancer and other diseases. Oxidative base lesions, including 7,8-dihydro-8-oxoguanine, are typically removed through base excision repair. In addition, oxidized deoxynucleotides such as 8-oxo-dGTP are depleted by sanitizing enzymes to preclude DNA incorporation. While pathways that counter threats posed by 7,8-dihydro-8-oxoguanine are well characterized, mechanisms protecting against the major adenine oxidation product, 7,8-dihydro-8-oxoadenine (oxoA), are poorly understood. Human DNA polymerases incorporate dGTP or dCTP opposite oxoA, producing mispairs that can cause A→C or A→G mutations. oxoA also perturbs the activity of enzymes acting on DNA and causes interstrand crosslinks. To inform mechanisms for oxoA repair, we characterized oxoA excision by human thymine DNA glycosylase (TDG), an enzyme known to remove modified pyrimidines, including deaminated and oxidized forms of cytosine and 5-methylcystosine. Strikingly, TDG excises oxoA from G⋅oxoA, A⋅oxoA, or C⋅oxoA pairs much more rapidly than it acts on the established pyrimidine substrates, whereas it exhibits comparable activity for T⋅oxoA and pyrimidine substrates. The oxoA activity depends strongly on base pairing and is 370-fold higher for G⋅oxoA versus T⋅oxoA pairs. The intrinsically disordered regions of TDG contribute minimally to oxoA excision, whereas two conserved residues (N140 and N191) are catalytically essential. Escherichia coli mismatch-specific uracil DNA-glycosylase lacks significant oxoA activity, exhibiting excision rates 4 to 5 orders of magnitude below that of its ortholog, TDG. Our results reveal oxoA as an unexpectedly efficient purine substrate for TDG and underscore the large evolutionary divergence of TDG and mismatch-specific uracil DNA-glycosylase.

4'-SCF3 -Labeling Constitutes a Sensitive 19 F NMR Probe for Characterization of Interactions in the Minor Groove of DNA

Fluorinated nucleotides are invaluable for ¹⁹ F NMR studies of nucleic acid structure and function. Here, we synthesized 4'-SCF₃ -thymidine (T 4 ' - SCF 3 ${{^{4{^\prime}\hbox{-}{\rm SCF}{_{3}}}}}$ ) and incorporated it into DNA by means of solid-phase DNA synthesis. NMR studies showed that the 4'-SCF₃ group exhibited a flexible orientation in the minor groove of DNA duplexes and was well accommodated by various higher order DNA structures. The three magnetically equivalent fluorine atoms in 4'-SCF₃ -DNA constitute an isolated spin system, offering high ¹⁹ F NMR sensitivity and excellent resolution of the positioning of T 4 ' - SCF 3 ${{^{4{^\prime}\hbox{-}{\rm SCF}{_{3}}}}}$ within various secondary and tertiary DNA structures. The high structural adaptability and high sensitivity of T 4 ' - SCF 3 ${{^{4{^\prime}\hbox{-}{\rm SCF}{_{3}}}}}$ make it a valuable ¹⁹ F NMR probe for quantitatively distinguishing diverse DNA structures with single-nucleotide resolution and for monitoring the dynamics of interactions in the minor groove of double-stranded DNA.

Measuring thermodynamic preferences to form non-native conformations in nucleic acids using ultraviolet melting

Thermodynamic preferences to form non-native conformations are crucial for understanding how nucleic acids fold and function. However, they are difficult to measure experimentally because this requires accurately determining the population of minor low-abundance (<10%) conformations in a sea of other conformations. Here, we show that melting experiments enable facile measurements of thermodynamic preferences to adopt nonnative conformations in DNA and RNA. The key to this "delta-melt" approach is to use chemical modifications to render specific minor non-native conformations the major state. The validity and robustness of delta-melt is established for four different non-native conformations under various physiological conditions and sequence contexts through independent measurements of thermodynamic preferences using NMR. Delta-melt is faster relative to NMR, simple, and cost-effective and enables thermodynamic preferences to be measured for exceptionally low-populated conformations. Using delta-melt, we obtained rare insights into conformational cooperativity, obtaining evidence for significant cooperativity (1.0 to 2.5 kcal/mol) when simultaneously forming two adjacent Hoogsteen base pairs. We also measured the thermodynamic preferences to form G-C⁺ and A-T Hoogsteen and A-T base open states for nearly all 16 trinucleotide sequence contexts and found distinct sequence-specific variations on the order of 2 to 3 kcal/mol. This rich landscape of sequence-specific non-native minor conformations in the DNA double helix may help shape the sequence specificity of DNA biochemistry. Thus, melting experiments can now be used to access thermodynamic information regarding regions of the free energy landscape of biomolecules beyond the native folded and unfolded conformations.

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.

Kinetics and thermodynamics of BI-BII interconversion altered by T:G mismatches in DNA

T:G mismatches in DNA result in humans primarily from deamination of methylated CpG sites. They are repaired by redundant systems, such as thymine DNA glycosylase (TDG) and methyl-binding domain enzyme (MBD4), and maintenance of these sites has been implicated in epigenetic processes. The process by which these enzymes identify a canonical DNA base in the incorrect basepairing context remains a mystery. However, the conserved contacts of the repair enzymes with the DNA backbone suggests a role for protein-phosphate interaction in the recognition and repair processes. We have used 31P NMR to investigate the energetics of DNA backbone BI-BII interconversion, and for this work have focused on alterations to the activation barriers to interconversion and the effect of a mismatch compared with canonical DNA. We have found that alterations to the ΔG of interconversion for T:G basepairs are remarkably similar to U:G basepairs in the form of stepwise differences in ΔG of 1-2 kcal/mol greater than equivalent steps in unmodified DNA, suggesting a universality of this result for TDG substrates. Likewise, we see perturbations to the free energy (∼1 kcal/mol) and enthalpy (2-5 kcal/mol) of activation for the BI-BII interconversion localized to the phosphates flanking the mismatch. Overall our results strongly suggest that the perturbed backbone energetics in T:G basepairs play a significant role in the recognition process of DNA repair enzymes.

Energy Landscapes for Base-Flipping in a Model DNA Duplex

We explore the process of base-flipping for four central bases, adenine, guanine, cytosine, and thymine, in a deoxyribonucleic acid (DNA) duplex using the energy landscape perspective. NMR imino-proton exchange and fluorescence correlation spectroscopy studies have been used in previous experiments to obtain lifetimes for bases in paired and extrahelical states. However, the difference of almost 4 orders of magnitude in the base-flipping rates obtained by the two methods implies that they are exploring different pathways and possibly different open states. Our results support the previous suggestion that minor groove opening may be favored by distortions in the DNA backbone and reveal links between sequence effects and the direction of opening, i.e., whether the base flips toward the major or the minor groove side. In particular, base flipping along the minor groove pathway was found to align toward the 5′ side of the backbone. We find that bases align toward the 3′ side of the backbone when flipping along the major groove pathway. However, in some cases for cytosine and thymine, the base flipping along the major groove pathway also aligns toward the 5′ side. The sequence effect may be caused by the polar interactions between the flipping-base and its neighboring bases on either of the strands. For guanine flipping toward the minor groove side, we find that the equilibrium constant for opening is large compared to flipping via the major groove. We find that the estimated rates of base opening, and hence the lifetimes of the closed state, obtained for thymine flipping through small and large angles along the major groove differ by 6 orders of magnitude, whereas for thymine flipping through small angles along the minor groove and large angles along the major groove, the rates differ by 3 orders of magnitude.

Resolution of chemical shift anisotropy in 19F ENDOR spectroscopy at 263 GHz/9.4 T

Pulsed ¹⁹F ENDOR spectroscopy provides a selective method for measuring angstrom to nanometer distances in structural biology. Here, the performance of ¹⁹F ENDOR at fields of 3.4 T and 9.4 T is compared using model compounds containing one to three ¹⁹F atoms. CF₃ groups are included in two compounds, for which the possible occurrence of uniaxial rotation might affect the distance distribution. At 9.4 T, pronounced asymmetric features are observed in many of the presented ¹⁹F ENDOR spectra. Data analysis by spectral simulations shows that these features arise from the chemical shift anisotropy (CSA) of the ¹⁹F nuclei. This asymmetry is also observed at 3.4 T, albeit to a much smaller extent, confirming the physical origin of the effect. The CSA parameters are well consistent with DFT predicted values and can be extracted from simulation of the experimental data in favourable cases, thereby providing additional information about the geometrical and electronic structure of the spin system. The feasibility of resolving the CSA at 9.4 T provides important information for the interpretation of line broadening in ENDOR spectra also at lower fields, which is relevant for developing methods to extract distance distributions from ¹⁹F ENDOR spectra.

Natural deep eutectic solvents as green and biocompatible reaction medium for carbonic anhydrase catalysis

Easy deactivation of free enzymes under non-native condition has become a stumbling block to the industrial application of biocatalysis. Natural deep eutectic solvent (NADES) has been exploited as a novel reaction medium for improving enzyme stability. The present work focused on preserving and enhancing the activity of carbonic anhydrase (CA) in a more economical and biocompatible NADES system. We synthesized six choline chloride/betaine-based NADES and analyzed the effects of compositions and concentrations of NADES on their physicochemical properties. The Bet-Gly (1: 2) NADES (55%) was proved to be more suitable as reaction medium for CA by analyzing enzyme activity in the presence of NADES. The enhancement in the stability of CA was found to be as a result of a three-dimensional hydrogen bonding network, rather than the individual or the synergistic effect of betaine and glyceride. The conformational change of CA to become more compact was confirmed both by fluorescence spectrum analysis and circular dichroism analysis. It is worth mentioning that a remarkable thermal stability was maintained when CA was incubated at temperature below 60 °C, and about 96% of activity was still restored in 55% NADES at 60 °C for 12 h.

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.

Developments in solution-state NMR yield broader and deeper views of the dynamic ensembles of nucleic acids

Nucleic acids do not fold into a single conformation, and dynamic ensembles are needed to describe their propensities to cycle between different conformations when performing cellular functions. We review recent advances in solution-state nuclear magnetic resonance (NMR) methods and their integration with computational techniques that are improving the ability to probe the dynamic ensembles of DNA and RNA. These include computational approaches for predicting chemical shifts from structure and generating conformational libraries from sequence, measurements of exact nuclear Overhauser effects, development of new probes to study chemical exchange using relaxation dispersion, faster and more sensitive real-time NMR techniques, and new NMR approaches to tackle large nucleic acid assemblies. We discuss how these advances are leading to new mechanistic insights into gene expression and regulation.

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.

Inequality in the Frequency of the Open States Occurrence Depends on Single 2H/1H Replacement in DNA

In the present study, the effect of 2H/1H isotopic exchange in hydrogen bonds between nitrogenous base pairs on occurrence and open states zones dynamics is investigated. These processes are studied using mathematical modeling, taking into account the number of open states between base pairs. The calculations of the probability of occurrence of open states in different parts of the gene were done depending on the localization of the deuterium atom. The mathematical modeling study demonstrated significant inequality (dependent on single 2H/1H replacement in DNA) among three parts of the gene similar in length of the frequency of occurrence of the open states. In this paper, the new convenient approach of the analysis of the abnormal frequency of open states in different parts of the gene encoding interferon alpha 17 was presented, which took into account both rising and decreasing of them that allowed to make a prediction of the functional instability of the specific DNA regions. One advantage of the new algorithm is diminishing the number of both false positive and false negative results in data filtered by this approach compared to the pure fractile methods, such as deciles or quartiles.

2'-Fluorinated Hydantoins as Chemical Biology Tools for Base Excision Repair Glycosylases

The guanine oxidation products, 5-guanidinohydantoin (Gh) and spiroiminodihydantoin (Sp), are mutagenic and toxic base lesions that are removed by Fpg, Nei, and the Nei-like (NEIL) glycosylases as the first step in base excision repair (BER). The hydantoins are excellent substrates for the NEIL glycosylases in a variety of DNA contexts beyond canonical duplex DNA, implicating the potential impact of repair activity on a multitude of cellular processes. In order to prepare stable derivatives as chemical biology tools, oligonucleotides containing fluorine at the 2'-position of the sugar of 8-oxo-7,8-dihydro-2'-deoxyguanosine2'-F-OG) were synthesized in ribo and arabino configuration. Selective oxidation of 2'-F-OG within a DNA oligonucleotide provided the corresponding 2'-F-Gh or 2'-F-Sp containing DNA. The 2'-F-hydantoins in duplex DNA were found to be highly resistant to the glycosylase activity of Fpg and NEIL1 compared to the unmodified lesion substrates. Surprisingly, however, some glycosylase-mediated base removal from both the 2'-F-ribo- and 2'-F-arabinohydantoin duplex DNA was observed. Notably, the associated β-lyase strand scission reaction of the 2'-F-arabinohydantoins was inhibited such that the glycosylases were "stalled" at the Schiff-base intermediate. Fpg and NEIL1 showed high affinity for the 2'-F-Gh duplexes in both ribo and arabino configurations. However, binding affinity assessed using catalytically inactive variants of Fpg and NEIL1 indicated higher affinity for the 2'-F-riboGh-containing duplexes. The distinct features of glycosylase processing of 2'-F-ribohydantoins and 2'-F-arabinohydantoins illustrate their utility to reveal structural insight into damage recognition and excision by NEIL and related glycosylases and provide opportunities for delineating the impact of lesion formation and repair in cells.

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.

Formation and biological consequences of 5-Formylcytosine in genomic DNA

5-Formyl-2'-deoxycytidine (5fdC) is a naturally occurring nucleobase that is broadly distributed in genomic DNA. 5fdC is produced via the oxidation of 5-methylcytosine (5mdC) by ten-eleven translocation enzyme (TET) and can be further converted to 5-carboxylcytosine (5cadC) by TET. Both 5fdC and 5cadC can be restored to dC by TDG-mediated base excision repair and direct deformylation/decarboxylation. Thus, 5fdC is considered an intermediate in the TET-mediated DNA demethylation pathway. 5fdC also alters the structure and stability of genomic DNA and affects genetic expression. This review summarizes the recent research on 5fdC, detailing its formation, detection and distribution, biological functions and transformation in cells. The challenges and future prospects to further explore the function and metabolism of 5fdC are briefly discussed at the end.