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

Structural underpinnings of mutation rate variations in the human genome

Single nucleotide mutation rates have critical implications for human evolution and genetic diseases. Importantly, the rates vary substantially across the genome and the principles underlying such variations remain poorly understood. A recent model explained much of this variation by considering higher-order nucleotide interactions in the 7-mer sequence context around mutated nucleotides. This model's success implicates a connection between DNA shape and mutation rates. DNA shape, i.e. structural properties like helical twist and tilt, is known to capture interactions between nucleotides within a local context. Thus, we hypothesized that changes in DNA shape features at and around mutated positions can explain mutation rate variations in the human genome. Indeed, DNA shape-based models of mutation rates showed similar or improved performance over current nucleotide sequence-based models. These models accurately characterized mutation hotspots in the human genome and revealed the shape features whose interactions underlie mutation rate variations. DNA shape also impacts mutation rates within putative functional regions like transcription factor binding sites where we find a strong association between DNA shape and position-specific mutation rates. This work demonstrates the structural underpinnings of nucleotide mutations in the human genome and lays the groundwork for future models of genetic variations to incorporate DNA shape.

DNA Methylation: Genomewide Distribution, Regulatory Mechanism and Therapy Target

DNA methylation is the most important epigenetic modification involved in the regulation of transcription, imprinting, establishment of X-inactivation, and the formation of a chromatin structure. DNA methylation in the genome is often associated with transcriptional repression and the formation of closed heterochromatin. However, the results of genome-wide studies of the DNA methylation pattern and transcriptional activity of genes have nudged us toward reconsidering this paradigm, since the promoters of many genes remain active despite their methylation. The differences in the DNA methylation distribution in normal and pathological conditions allow us to consider methylation as a diagnostic marker or a therapy target. In this regard, the need to investigate the factors affecting DNA methylation and those involved in its interpretation becomes pressing. Recently, a large number of protein factors have been uncovered, whose ability to bind to DNA depends on their methylation. Many of these proteins act not only as transcriptional activators or repressors, but also affect the level of DNA methylation. These factors are considered potential therapeutic targets for the treatment of diseases resulting from either a change in DNA methylation or a change in the interpretation of its methylation level. In addition to protein factors, a secondary DNA structure can also affect its methylation and can be considered as a therapy target. In this review, the latest research into the DNA methylation landscape in the genome has been summarized to discuss why some DNA regions avoid methylation and what factors can affect its level or interpretation and, therefore, can be considered a therapy target.

DNA Deformation Exerted by Regulatory DNA-Binding Motifs in Human Alkyladenine DNA Glycosylase Promotes Base Flipping

Human alkyladenine DNA glycosylase (AAG) is a key enzyme that corrects a broad range of alkylated and deaminated nucleobases to maintain genomic integrity. When encountering the lesions, AAG adopts a base-flipping strategy to extrude the target base from the DNA duplex to its active site, thereby cleaving the glycosidic bond. Despite its functional importance, the detailed mechanism of such base extrusion and how AAG distinguishes the lesions from an excess of normal bases both remain elusive. Here, through the Markov state model constructed on extensive all-atom molecular dynamics simulations, we find that the alkylated nucleobase (N3-methyladenine, 3MeA) everts through the DNA major groove. Two key AAG motifs, the intercalation and E131-N146 motifs, play active roles in bending/pressing the DNA backbone and widening the DNA minor groove during 3MeA eversion. In particular, the intercalated residue Y162 is involved in buckling the target site at the early stage of 3MeA eversion. Our traveling-salesman based automated path searching algorithm further revealed that a non-target normal adenine tends to be trapped in an exo site near the active site, which however barely exists for a target base 3MeA. Collectively, these results suggest that the Markov state model combined with traveling-salesman based automated path searching acts as a promising approach for studying complex conformational changes of biomolecules and dissecting the elaborate mechanism of target recognition by this unique enzyme.

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.

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.

Oxidation Enhances Binding of Extrahelical 5-Methyl-Cytosines by Thymine DNA Glycosylase

The DNA repair protein thymine DNA glycosylase (TDG) removes mispaired or damaged bases, such as oxidized methyl-cytosine, from DNA by cleavage of the glycosidic bond between the sugar and the target base flipped into the enzyme's active site. The enzyme is active against formyl-cytosine and carboxyl-cytosine, whereas the lower oxidized hydroxymethyl-cytosine and methyl-cytosine itself are not processed by the enzyme. Molecular dynamics simulations with thermodynamic integration of TDG complexed to DNA carrying one of four different (oxidized) methyl-cytosine bases in extrahelcial conformation, methyl-cytosine (mC), hydroxymethyl-cytosine (hmC), formyl-cytosine (fC), or carboxyl-cytosine (caC), show a more favorable binding affinity of the higher oxidized forms, fC and caC, than the nonsubstrate bases hmC and mC. Despite rather comparable, reaction-competent conformations of the flipped bases in the active site of the enzyme, more and stronger interactions with active site residues account for the preferred binding of the higher oxidized bases. Binding of the negatively charged caC and the neutral fC are strengthened by interactions with positively charged His151. Our calculated proton affinities find this protonation state of His151 the preferred one in the presence of caC and conceivable in the presence of fC as well as increasing the binding affinity toward the two bases. Discrimination of the substrate bases is further achieved by the backbone of Tyr152 that forms a strong hydrogen bond to the carboxyl and formyl oxygen atoms of caC and fC, respectively, a contact that is completely lacking in mC and much weaker in hmC. Overall, our computational results indicate that the enzyme discriminates the different oxidation forms of methyl-cytosine already at the formation of the extrahelical complexes.

The base excision repair process: comparison between higher and lower eukaryotes

The base excision repair (BER) pathway is essential for maintaining the stability of DNA in all organisms and defects in this process are associated with life-threatening diseases. It is involved in removing specific types of DNA lesions that are induced by both exogenous and endogenous genotoxic substances. BER is a multi-step mechanism that is often initiated by the removal of a damaged base leading to a genotoxic intermediate that is further processed before the reinsertion of the correct nucleotide and the restoration of the genome to a stable structure. Studies in human and yeast cells, as well as fruit fly and nematode worms, have played important roles in identifying the components of this conserved DNA repair pathway that maintains the integrity of the eukaryotic genome. This review will focus on the components of base excision repair, namely, the DNA glycosylases, the apurinic/apyrimidinic endonucleases, the DNA polymerase, and the ligases, as well as other protein cofactors. Functional insights into these conserved proteins will be provided from humans, Saccharomyces cerevisiae, Drosophila melanogaster, and Caenorhabditis elegans, and the implications of genetic polymorphisms and knockouts of the corresponding genes.

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.

Polymerization and editing modes of a high-fidelity DNA polymerase are linked by a well-defined path

Proofreading by replicative DNA polymerases is a fundamental mechanism ensuring DNA replication fidelity. In proofreading, mis-incorporated nucleotides are excised through the 3'-5' exonuclease activity of the DNA polymerase holoenzyme. The exonuclease site is distal from the polymerization site, imposing stringent structural and kinetic requirements for efficient primer strand transfer. Yet, the molecular mechanism of this transfer is not known. Here we employ molecular simulations using recent cryo-EM structures and biochemical analyses to delineate an optimal free energy path connecting the polymerization and exonuclease states of E. coli replicative DNA polymerase Pol III. We identify structures for all intermediates, in which the transitioning primer strand is stabilized by conserved Pol III residues along the fingers, thumb and exonuclease domains. We demonstrate switching kinetics on a tens of milliseconds timescale and unveil a complete pol-to-exo switching mechanism, validated by targeted mutational experiments.

HMGB1 Recruits TET2/AID/TDG to Induce DNA Demethylation in STAT3 Promoter in CD4+ T Cells from aGVHD Patients

STAT3 is highly expressed in aGVHD CD4⁺ T cells and plays a critical role in inducing or worsening aGVHD. In our preceding studies, DNA hypomethylation in STAT3 promoter was shown to cause high expression of STAT3 in aGVHD CD4⁺ T cells, and the process could be modulated by HMGB1, but the underlying mechanism remains unclear. TET2, AID, and TDG are indispensable in DNA demethylation; meanwhile, TET2 and AID also serve extremely important roles in immune response. So, we speculated these enzymes involved in the STAT3 promoter hypomethylation induced by HMGB1 in aGVHD CD4⁺ T cells. In this study, we found that the binding levels of TET2/AID/TDG to STAT3 promoter were remarkably increased in CD4⁺T cells from aGVHD patients and were significantly negatively correlated with the STAT3 promoter methylation level. Simultaneously, we revealed that HMGB1 could recruit TET2, AID, and TDG to form a complex in the STAT3 promoter region. Interference with the expression of TET2/AID/TDG inhibited the overexpression of STAT3 caused by HMGB1 downregulation of the STAT3 promoter DNA methylation. These data demonstrated a new molecular mechanism of how HMGB1 promoted the expression of STAT3 in CD4⁺ T cells from aGVHD patients.

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.

Disclosing the Mechanism of Spontaneous Aggregation and Template-Induced Misfolding of the Key Hexapeptide (PHF6) of Tau Protein Based on Molecular Dynamics Simulation

The microtubule-associated protein tau is critical for the development and maintenance of the nervous system. Tau dysfunction is associated with a variety of neurodegenerative diseases called tauopathies, which are characterized by neurofibrillary tangles formed by abnormally aggregated tau protein. Studying the aggregation mechanism of tau protein is of great significance for elucidating the etiology of tauopathies. The hexapeptide ³⁰⁶VQIVYK³¹¹ (PHF6) of R3 has been shown to play a vital role in promoting tau aggregation. In this study, long-term all-atom molecular dynamics simulations in explicit solvent were performed to investigate the mechanisms of spontaneous aggregation and template-induced misfolding of PHF6, and the dimerization at the early stage of nucleation was further specifically analyzed by the Markov state model (MSM). Our results show that PHF6 can spontaneously aggregate to form multimers enriched with β-sheet structure and the β-sheets in multimers prefer to exist in a parallel way. It is observed that PHF6 monomer can be induced to form a β-sheet structure on either side of the template but in a different way. In detail, the β-sheet structure is easier to form on the left side but does not extend well, but on the right side, the monomer can form the extended β-sheet structure. Furthermore, MSM analysis shows that the formation of dimer mainly occurs in three steps. First, the separated monomers collide with each other at random orientations, and then a dimer with short β-sheet structure at the N-terminal forms; finally, β-sheets elongate to form an extended parallel β-sheet dimer. During these processes, multiple intermediate states are identified and multiple paths can form a parallel β-sheet dimer from the disordered coil structure. Moreover, the residues I308, V309, and Y310 play an essential role in the dimerization. In a word, our results uncover the aggregation and misfolding mechanism of PHF6 from the atomic level, which can provide useful theoretical guidance for rational design of effective therapeutic drugs against tauopathies.

Effect of the thymine-DNA glycosylase rs4135050 variant on Saudi smoker population

Background

Thymine‐DNA glycosylase (TDG) is an essential DNA‐repair enzyme which works in both epigenetic regulation and genome maintenance. It is also responsible for efficient correction of multiple endogenous DNA lesions which occur commonly in mammalian genomes. Research of genetic variants such as SNPs, resulting in disease, is predicted to yield clinical advancements through the identification of sensitive genetic markers and the development of disease prevention and therapy. To that end, the main objective of the present study is to identify the possible interactions between cigarette smoking and the rs4135050 variant of the TDG gene, situated in the intron position, among Saudi individuals.

Methods

TDG rs4135050 (A/T) was investigated by genotyping 239, and 235 blood specimens were obtained from nonsmokers and smokers of cigarette respectively.

Results

T allele frequency was found which showed a significant protective effect on Saudi male smokers (OR = 0.64, p = 0.0187) compared to nonsmoking subjects, but not in female smokers. Furthermore, smokers aged less than 29 years, the AT and AT+TT genotypes decreased more than four times the risk of initiation of smoking related‐diseases compare to the ancestral AA homozygous genotype. Paradoxically, the AT (OR = 3.88, p = 0.0169) and AT+TT (OR = 2.86, p = 0.0420) genotypes were present at a higher frequency in smoking patients aged more than 29 years as compared to nonsmokers at the same ages.

Conclusion

Depending on the gender and age of patients, TDG rs4135050 may provide a novel biomarker for the early diagnosis and prevention of several diseases caused by cigarette smoking.

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.