Chemical Methods to Identify Epigenetic Modifications in Cytosine Bases

Chemical modifications to Cytosine bases are among the most studied epigenetic markers and their detection in the human genome plays a crucial role in gaining more insights about gene regulation, prognosis of genetic disorders and unraveling genetic inheritance patterns. The Cytosine methylated at the 5th position and oxidized derivatives thereof generated in the demethylation pathways, perform separate and unique epigenetic functions in an organism. As the presence of various Cytosine modifications is associated with diverse diseases, including cancer, there has been a strong focus on developing methods, both chemical and alternative approaches, capable of detecting these modifications at a single-base resolution across the entire genome. In this comprehensive review, we aim to consolidate the various chemical methods and understanding their chemistry that have been established to date for the detection of various Cytosine modifications.

Is DNA methylation in the brain a mechanism of alcohol use disorder?

Alcohol use disorder (AUD) is a worldwide problem. Unfortunately, the molecular mechanisms of alcohol misuse are still poorly understood, therefore successful therapeutic approaches are limited. Accumulating data indicate that the tendency for compulsive alcohol use is inherited, suggesting a genetic background as an important factor. However, the probability to develop AUD is also affected by life experience and environmental factors. Therefore, the epigenetic modifications that are altered over lifetime likely contribute to increased risk of alcohol misuse. Here, we review the literature looking for the link between DNA methylation in the brain, a common epigenetic modification, and AUD-related behaviors in humans, mice and rats. We sum up the main findings, identify the existing gaps in our knowledge and indicate future directions of the research.

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.

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.

Structural basis of DNA methylation-dependent site selectivity of the Epstein-Barr virus lytic switch protein ZEBRA/Zta/BZLF1

In infected cells, Epstein-Barr virus (EBV) alternates between latency and lytic replication. The viral bZIP transcription factor ZEBRA (Zta, BZLF1) regulates this cycle by binding to two classes of ZEBRA response elements (ZREs): CpG-free motifs resembling the consensus AP-1 site recognized by cellular bZIP proteins and CpG-containing motifs that are selectively bound by ZEBRA upon cytosine methylation. We report structural and mutational analysis of ZEBRA bound to a CpG-methylated ZRE (meZRE) from a viral lytic promoter. ZEBRA recognizes the CpG methylation marks through a ZEBRA-specific serine and a methylcytosine-arginine-guanine triad resembling that found in canonical methyl-CpG binding proteins. ZEBRA preferentially binds the meZRE over the AP-1 site but mutating the ZEBRA-specific serine to alanine inverts this selectivity and abrogates viral replication. Our findings elucidate a DNA methylation-dependent switch in ZEBRA's transactivation function that enables ZEBRA to bind AP-1 sites and promote viral latency early during infection and subsequently, under appropriate conditions, to trigger EBV lytic replication by binding meZREs.

Towards a Better Understanding of Computational Models for Predicting DNA Methylation Effects at the Molecular Level

Human DNA is a very sensitive macromolecule and slight changes in the structure of DNA can have disastrous effects on the organism. When nucleotides are modified, or changed, the resulting DNA sequence can lose its information, if it is part of a gene, or it can become a problem for replication and repair. Human cells can regulate themselves by using a process known as DNA methylation. This methylation is vitally important in cell differentiation and expression of genes. When the methylation is uncontrolled, however, or does not occur in the right place, serious pathophysiological consequences may result. Excess methylation causes changes in the conformation of the DNA double helix. The secondary structure of DNA is highly dependent upon the sequence. Therefore, if the sequence changes slightly the secondary structure can change as well. These slight changes will then cause the doublestranded DNA to be more open and available in some places where large adductions can come in and react with the DNA base pairs. Computer models have been used to simulate a variety of biological processes including protein function and binding, and there is a growing body of evidence that in silico methods can shed light on DNA methylation. Understanding the anomeric effect that contributes to the structural and conformational flexibility of furanose rings through a combination of quantum mechanical and experimental studies is critical for successful molecular dynamic simulations.

Aqueous Wittig reaction-mediated fast fluorogenic identification and single-base resolution analysis of 5-formylcytosine in DNA

A highly reactive ylide tBuA was introduced, which could react rapidly with the 5-formyl and 4-amino groups of 5-formylcytosine (5fC) under mild conditions without any co-solvent or catalyst in a manner of Wittig olefination and intramolecular nucleophilic substitution to yield a cyclized fluorescent adduct, which meets the demands of both single-base resolution sequencing and fluorescence switch-on detection of 5fC in DNA.

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.

2'-Alkynyl spin-labelling is a minimally perturbing tool for DNA structural analysis

The determination of distances between specific points in nucleic acids is essential to understanding their behaviour at the molecular level. The ability to measure distances of 2-10 nm is particularly important: deformations arising from protein binding commonly fall within this range, but the reliable measurement of such distances for a conformational ensemble remains a significant challenge. Using several techniques, we show that electron paramagnetic resonance (EPR) spectroscopy of oligonucleotides spin-labelled with triazole-appended nitroxides at the 2' position offers a robust and minimally perturbing tool for obtaining such measurements. For two nitroxides, we present results from EPR spectroscopy, X-ray crystal structures of B-form spin-labelled DNA duplexes, molecular dynamics simulations and nuclear magnetic resonance spectroscopy. These four methods are mutually supportive, and pinpoint the locations of the spin labels on the duplexes. In doing so, this work establishes 2'-alkynyl nitroxide spin-labelling as a minimally perturbing method for probing DNA conformation.

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

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

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

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

Evidence for novel epigenetic marks within plants

Variation in patterns of gene expression can result from modifications in the genome that occur without a change in the sequence of the DNA; such modifications include methylation of cytosine to generate 5-methylcytosine (5mC) resulting in the generation of heritable epimutation and novel epialleles. This type of non-sequence variation is called epigenetics. The enzymes responsible for generation of such DNA modifications in mammals are named DNA methyltransferases (DNMT) including DNMT1, DNMT2 and DNMT3. The later stages of oxidations to these modifications are catalyzed by Ten Eleven Translocation (TET) proteins, which contain catalytic domains belonging to the 2-oxoglutarate dependent dioxygenase family. In various mammalian cells/tissues including embryonic stem cells, cancer cells and brain tissues, it has been confirmed that these proteins are able to induce the stepwise oxidization of 5-methyl cytosine to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and finally 5-carboxylcytosine (5caC). Each stage from initial methylation until the end of the DNA demethylation process is considered as a specific epigenetic mark that may regulate gene expression. This review discusses controversial evidence for the presence of such oxidative products, particularly 5hmC, in various plant species. Whereas some reports suggest no evidence for enzymatic DNA demethylation, other reports suggest that the presence of oxidative products is followed by the active demethylation and indicate the contribution of possible TET-like proteins in the regulation of gene expression in plants. The review also summarizes the results obtained by expressing the human TET conserved catalytic domain in transgenic plants.

A method for identifying allele-specific hydroxymethylation

We previously identified sequence-dependent allele-specific methylation (sd-ASM) in adult human peripheral blood leukocytes, in which ASM occurs in cis depending on adjacent polymorphic sequences. A number of groups have identified sd-ASM sites in the human and mouse genomes, illustrating the prevalence of sd-ASM in mammalian genomes. In addition, sd-ASM can lead to sequence-dependent allele-specific expression of neighbouring genes. Imprinted genes also often exhibit parent-of-origin-dependent allele-specific methylation (pd-ASM), which causes parent-of-origin-dependent allele-specific expression. However, whether most of the already known sd-ASM and pd-ASM sites are methylated or hydroxymethylated remains unclear due to technical restrictions. Accordingly, a novel method that enables examination of allelic methylation and hydroxymethylation status and also overcomes the drawbacks of conventional methods is needed. Such a method could also be used to elucidate the mechanisms underlying polymorphism-associated inter-individual differences in disease susceptibility and the mechanism of genomic imprinting. Here, we developed a simple method to determine allelic hydroxymethylation status and identified novel sequence- and parent-of-origin-dependent allele-specific hydroxymethylation sites. Correlation analyses of TF binding sequences and methylation or hydroxymethylation between three mouse strains revealed the involvement of Pax5 in strain-specific methylation and hydroxymethylation in exon 7 of Pdgfrb.

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.

Epigenetic control of the angiotensin-converting enzyme in endothelial cells during inflammation

The angiotensin-converting enzyme (ACE) plays a central role in the renin-angiotensin system, which is involved in the regulation of blood pressure. Alterations in ACE expression or activity are associated with various pathological phenotypes, particularly cardiovascular diseases. In human endothelial cells, ACE was shown to be negatively regulated by tumor necrosis factor (TNF) α. To examine, whether or not, epigenetic factors were involved in ACE expression regulation, methylated DNA immunoprecipitation and RNA interference experiments directed against regulators of DNA methylation homeostasis i.e., DNA methyltransferases (DNMTs) and ten-eleven translocation methylcytosine dioxygenases (TETs), were performed. TNFα stimulation enhanced DNA methylation in two distinct regions within the ACE promoter via a mechanism linked to DNMT3a and DNMT3b, but not to DNMT1. At the same time, TET1 protein expression was downregulated. In addition, DNA methylation decreased the binding affinity of the transcription factor MYC associated factor X to the ACE promoter. In conclusion, DNA methylation determines the TNFα-dependent regulation of ACE gene transcription and thus protein expression in human endothelial cells.

NMR and MS-based Stable Isotope-Resolved Metabolomics and Applications in Cancer Metabolism

There is considerable interest in defining metabolic reprogramming in human diseases, which is recognized as a hallmark of human cancer. Although radiotracers have a long history in specific metabolic studies, stable isotope-enriched precursors coupled with modern high resolution mass spectrometry and NMR spectroscopy have enabled systematic mapping of metabolic networks and fluxes in cells, tissues and living organisms including humans. These analytical platforms are high in information content, are complementary and cross-validating in terms of compound identification, quantification, and isotope labeling pattern analysis of a large number of metabolites simultaneously. Furthermore, new developments in chemoselective derivatization and in vivo spectroscopy enable tracking of labile/low abundance metabolites and metabolic kinetics in real-time. Here we review developments in Stable Isotope Resolved Metabolomics (SIRM) and recent applications in cancer metabolism using a wide variety of stable isotope tracers that probe both broad and specific aspects of cancer metabolism required for proliferation and survival.

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.

Fluorescent Wittig reagent as a novel ratiometric probe for the quantification of 5-formyluracil and its application in cell imaging

For the first time a Wittig reagent was introduced into the design of a fluorescent probe for the quantification of 5-formyluracil.