Structural basis for effects of CpA modifications on C/EBPβ binding of DNA

Structural basis for specific DNA sequence recognition by the transcription factor NFIL3

The CCAAT/enhancer-binding proteins (C/EBPs) constitute a family of pivotal transcription factors involved in tissue development, cellular function, proliferation, and differentiation. NFIL3, as one of them, plays an important role in regulating immune cell differentiation, circadian clock system, and neural regeneration, yet its specific DNA recognition mechanism remains enigmatic. In this study, we showed by the ITC binding experiments that NFIL3 prefers to bind to the TTACGTAA DNA motif. Our structural studies revealed that the α-helical NFIL3 bZIP domain dimerizes through its leucine zipper region, and binds to DNA via its basic region. The two basic regions of the NFIL3 bZIP dimer were pushed apart upon binding to DNA, facilitating the snug accommodation of the two basic regions within the major grooves of the DNA. Remarkably, our binding and structural data also revealed that both NFIL3 and C/EBPα/β demonstrate a shared preference for the TTACGTAA sequence. Furthermore, our study revealed that disease-associated mutations within the NFIL3 bZIP domain result in either reduction or complete disruption of its DNA binding ability. These discoveries not only provide valuable insights into the DNA binding mechanisms of NFIL3 but also elucidate the causal role of NFIL3 mutations in disease pathogenesis.

The mitoepigenome responds to stress, suggesting novel mito-nuclear interactions in vertebrates

The mitochondria are central in the cellular response to changing environmental conditions resulting from disease states, environmental exposures or normal physiological processes. Although the influences of environmental stressors upon the nuclear epigenome are well characterized, the existence and role of the mitochondrial epigenome remains contentious. Here, by quantifying the mitochondrial epigenomic response of pineal gland cells to circadian stress, we confirm the presence of extensive cytosine methylation within the mitochondrial genome. Furthermore, we identify distinct epigenetically plastic regions (mtDMRs) which vary in cytosinic methylation, primarily in a non CpG context, in response to stress and in a sex-specific manner. Motifs enriched in mtDMRs contain recognition sites for nuclear-derived DNA-binding factors (ATF4, HNF4A) important in the cellular metabolic stress response, which we found to be conserved across diverse vertebrate taxa. Together, these findings suggest a new layer of mito-nuclear interaction in which the nuclear metabolic stress response could alter mitochondrial transcriptional dynamics through the binding of nuclear-derived transcription factors in a methylation-dependent context.

Structures of CTCF-DNA complexes including all 11 zinc fingers

The CCCTC-binding factor (CTCF) binds tens of thousands of enhancers and promoters on mammalian chromosomes by means of its 11 tandem zinc finger (ZF) DNA-binding domain. In addition to the 12-15-bp CORE sequence, some of the CTCF binding sites contain 5' upstream and/or 3' downstream motifs. Here, we describe two structures for overlapping portions of human CTCF, respectively, including ZF1-ZF7 and ZF3-ZF11 in complex with DNA that incorporates the CORE sequence together with either 3' downstream or 5' upstream motifs. Like conventional tandem ZF array proteins, ZF1-ZF7 follow the right-handed twist of the DNA, with each finger occupying and recognizing one triplet of three base pairs in the DNA major groove. ZF8 plays a unique role, acting as a spacer across the DNA minor groove and positioning ZF9-ZF11 to make cross-strand contacts with DNA. We ascribe the difference between the two subgroups of ZF1-ZF7 and ZF8-ZF11 to residues at the two positions -6 and -5 within each finger, with small residues for ZF1-ZF7 and bulkier and polar/charged residues for ZF8-ZF11. ZF8 is also uniquely rich in basic amino acids, which allows salt bridges to DNA phosphates in the minor groove. Highly specific arginine-guanine and glutamine-adenine interactions, used to recognize G:C or A:T base pairs at conventional base-interacting positions of ZFs, also apply to the cross-strand interactions adopted by ZF9-ZF11. The differences between ZF1-ZF7 and ZF8-ZF11 can be rationalized structurally and may contribute to recognition of high-affinity CTCF binding sites.

Physicochemical models of protein-DNA binding with standard and modified base pairs

DNA-binding proteins play important roles in various cellular processes, but the mechanisms by which proteins recognize genomic target sites remain incompletely understood. Functional groups at the edges of the base pairs (bp) exposed in the DNA grooves represent physicochemical signatures. As these signatures enable proteins to form specific contacts between protein residues and bp, their study can provide mechanistic insights into protein-DNA binding. Existing experimental methods, such as X-ray crystallography, can reveal such mechanisms based on physicochemical interactions between proteins and their DNA target sites. However, the low throughput of structural biology methods limits mechanistic insights for selection of many genomic sites. High-throughput binding assays enable prediction of potential target sites by determining relative binding affinities of a protein to massive numbers of DNA sequences. Many currently available computational methods are based on the sequence of standard Watson-Crick bp. They assume that the contribution of overall binding affinity is independent for each base pair, or alternatively include dinucleotides or short k-mers. These methods cannot directly expand to physicochemical contacts, and they are not suitable to apply to DNA modifications or non-Watson-Crick bp. These variations include DNA methylation, and synthetic or mismatched bp. The proposed method, DeepRec, can predict relative binding affinities as function of physicochemical signatures and the effect of DNA methylation or other chemical modifications on binding. Sequence-based modeling methods are in comparison a coarse-grain description and cannot achieve such insights. Our chemistry-based modeling framework provides a path towards understanding genome function at a mechanistic level.

Structural basis for cell type specific DNA binding of C/EBPβ: The case of cell cycle inhibitor p15INK4b promoter

C/EBPβ is a key regulator of numerous cellular processes, but it can also contribute to tumorigenesis and viral diseases. It binds to specific DNA sequences (C/EBP sites) and interacts with other transcription factors to control expression of multiple eukaryotic genes in a tissue and cell-type dependent manner. A body of evidence has established that cell-type-specific regulatory information is contained in the local DNA sequence of the binding motif. In human epithelial cells, C/EBPβ is an essential cofactor for TGFβ signaling in the case of Smad2/3/4 and FoxO-dependent induction of the cell cycle inhibitor, p15INK4b. In the TGFβ-responsive region 2 of the p15INK4b promoter, the Smad binding site is flanked by a C/EBP site, CTTAA•GAAAG, which differs from the canonical, palindromic ATTGC•GCAAT motif. The X-ray crystal structure of C/EBPβ bound to the p15INK4b promoter fragment shows how GCGC-to-AAGA substitution generates changes in the intermolecular interactions in the protein-DNA interface that enhances C/EBPβ binding specificity, limits possible epigenetic regulation of the promoter, and generates a DNA element with a unique pattern of methyl groups in the major groove. Significantly, CT/GA dinucleotides located at the 5'ends of the double stranded element maintain local narrowing of the DNA minor groove width that is necessary for DNA recognition. Our results suggest that C/EBPβ would accept all forms of modified cytosine in the context of the CpT site. This contrasts with the effect on the consensus motif, where C/EBPβ binding is modestly increased by cytosine methylation, but substantially decreased by hydroxymethylation.

Interaction preferences between protein side chains and key epigenetic modifications 5-methylcytosine, 5-hydroxymethycytosine and N6-methyladenine

Covalent modifications of standard DNA/RNA nucleobases affect epigenetic regulation of gene expression by modulating interactions between nucleic acids and protein readers. We derive here the absolute binding free energies and analyze the binding modalities between key modified nucleobases 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC) and N⁶-methyladenine (m⁶A) and all non-prolyl/non-glycyl protein side chains using molecular dynamics simulations and umbrella sampling in both water and methanol, the latter mimicking the low dielectric environment at the dehydrated nucleic-acid/protein interfaces. We verify the derived affinities by comparing against a comprehensive set of high-resolution structures of nucleic-protein complexes involving 5mC. Our analysis identifies protein side chains that are highly tuned for detecting cytosine methylation as a function of the environment and can thus serve as microscopic readers of epigenetic marks. Conversely, we show that the relative ordering of sidechain affinities for 5hmC and m⁶A does not differ significantly from those for their precursor bases, cytosine and adenine, respectively, especially in the low dielectric environment. For those two modified bases, the effect is more nuanced and manifests itself primarily at the level of absolute changes in the binding free energy. Our results contribute towards establishing a quantitative foundation for understanding, predicting and modulating the interactions between modified nucleic acids and proteins at the atomistic level.

Aberrant DNA hydroxymethylation reshapes transcription factor binding in myeloid neoplasms

Epigenetic abnormalities in DNA hydroxymethylation (5hmC) have been detected in patients with myeloid neoplasms, suggesting that 5hmC might act as a valuable epigenetic mark to reflect the disease status of myeloid neoplasms. Here, we report systematic genome-wide mapping of the DNA hydroxymethylomes in over 70 patients with myeloid neoplasms. Our integrative analysis leads to the identification of distinct 5hmC signatures that can sensitively discriminate patients from healthy individuals. At the molecular level, we unveiled dynamic 5hmC changes within key transcription factor (e.g., the CEBP family) binding motifs that are essential for hematopoiesis and myeloid lineage specification. 5hmC redistribution was found to alter the genome-wide binding of CEBP-α, thereby reprogramming transcriptional outputs to affect leukemia cell survival and stemness. Taken together, we provide a comprehensive 5hmC atlas representative of myeloid neoplasms, which sets the stage for future exploration on the epigenetic etiology of hematological malignancies. Mechanistically, our study further furnishes important insights into how abnormal 5hmC distribution in patients directly interrupts the binding of transcription factors to reshape transcriptional landscapes and aggravate leukemogenesis.

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.

Structural insights into glucocorticoid receptor function

The glucocorticoid receptor (GR) is a steroid hormone-activated transcription factor that binds to various glucocorticoid response elements to up- or down- regulate the transcription of thousands of genes involved in metabolism, development, stress and inflammatory responses. GR consists of two domains enabling interaction with glucocorticoids, DNA response elements and coregulators, as well as a large intrinsically disordered region that mediates condensate formation. A growing body of structural studies during the past decade have shed new light on GR interactions, providing a new understanding of the mechanisms driving context-specific GR activity. Here, we summarize the established and emerging mechanisms of action of GR, primarily from a structural perspective. This minireview also discusses how the current state of knowledge of GR function may guide future glucocorticoid design with an improved therapeutic index for different inflammatory disorders.

Structural basis for glucocorticoid receptor recognition of both unmodified and methylated binding sites, precursors of a modern recognition element

The most common form of DNA methylation involves the addition of a methyl group to a cytosine base in the context of a cytosine–phosphate–guanine (CpG) dinucleotide. Genomes from more primitive organisms are more abundant in CpG sites that, through the process of methylation, deamination and subsequent mutation to thymine–phosphate–guanine (TpG) sites, can produce new transcription factor binding sites. Here, we examined the evolutionary history of the over 36 000 glucocorticoid receptor (GR) consensus binding motifs in the human genome and identified a subset of them in regulatory regions that arose via a deamination and subsequent mutation event. GR can bind to both unmodified and methylated pre-GR binding sequences (GBSs) that contain a CpG site. Our structural analyses show that CpG methylation in a pre-GBS generates a favorable interaction with Arg447 mimicking that made with a TpG in a GBS. This methyl-specific recognition arose 420 million years ago and was conserved during the evolution of GR and likely helps fix the methylation on the relevant cytosines. Our study provides the first genetic, biochemical and structural evidence of high-affinity binding for the likely evolutionary precursor of extant TpG-containing GBS.

Molecular mechanism of methyl-dependent and spatial-specific DNA recognition of c-Jun homodimer

DNA methylation is important in regulation of gene expression and normal development because it alters the interplay between protein and DNA. Experiments have shown that a single 5-methylcytosine at different CpG sites (mCpG) might have different effects on specific recognition, but the atomistic origin and dynamic details are largely unclear. In this work, we investigated the mechanism of monomethylation at different CpG sites in the cognate motif and the cooperativity of full methylation. By constructing four models of c-Jun/Jun protein binding to the 5[Formula: see text]-XGAGTCA-3[Formula: see text] (X represents C or methylated C) motif, we characterized the dynamics of the contact interface using the all-atom molecular dynamics method. Free energy analysis of MM/GBSA suggests that regardless of whether the C12pG13 site of the bottom strand is methylated, the effects from mC25 of the top strand are dominant and can moderately enhance the binding by [Formula: see text] 31 kcal/mol, whereas mC12 showed a relatively small contribution, in agreement with the experimental data. Remarkably, we found that this spatial-specific influence was induced by different regulatory rules. The influence of the mC25 site is mainly mediated by steric hindrance. The additional methyl group leads to the conformational changes in nearby residues and triggers an obvious structural bending in the protein, which results in the formation of a new T-Asn-C triad that enhances the specific recognition of TCA half-sites. The substitution of the methyl group at the mC12 site of the bottom strand breaks the original H-bonds directly. Such changes in electrostatic interactions also lead to the remote allosteric effects of protein by multifaceted interactions but have negligible contributions to binding. Although these two influence modes are different, they can both fine-tune the local environment, which might produce remote allosteric effects through protein-protein interactions. Further analysis reveals that the discrepancies in these two modes are primarily due to their location. Moreover, when both sites are methylated, the major determinant of binding specificity depends on the context and the location of the methylation site, which is the result of crosstalk and cooperativity.

Enhanced C/EBP binding to G·T mismatches facilitates fixation of CpG mutations in cancer and adult stem cells

Somatic mutations in regulatory sites of human stem cells affect cell identity or cause malignant transformation. By mining the human genome for co-occurrence of mutations and transcription factor binding sites, we show that C/EBP binding sites are strongly enriched with [C > T]G mutations in cancer and adult stem cells, which is of special interest because C/EBPs regulate cell fate and differentiation. In vitro protein-DNA binding assay and structural modeling of the CEBPB-DNA complex show that the G·T mismatch in the core CG dinucleotide strongly enhances affinity of the binding site. We conclude that enhanced binding of C/EBPs shields CpG·TpG mismatches from DNA repair, leading to selective accumulation of [C > T]G mutations and consequent deterioration of the binding sites. This mechanism of targeted mutagenesis highlights the effect of a mutational process on certain regulatory sites and reveals the molecular basis of putative regulatory alterations in stem cells.

Preferential CEBP binding to T:G mismatches and increased C-to-T human somatic mutations

DNA cytosine methylation in mammals modulates gene expression and chromatin accessibility. It also impacts mutation rates, via spontaneous oxidative deamination of 5-methylcytosine (5mC) to thymine. In most cases the resulting T:G mismatches are repaired, following T excision by one of the thymine DNA glycosylases, TDG or MBD4. We found that C-to-T mutations are enriched in the binding sites of CCAAT/enhancer binding proteins (CEBP). Within a CEBP site, the presence of a T:G mismatch increased CEBPβ binding affinity by a factor of >60 relative to the normal C:G base pair. This enhanced binding to a mismatch inhibits its repair by both TDG and MBD4 in vitro. Furthermore, repair of the deamination product of unmethylated cytosine, which yields a U:G DNA mismatch that is normally repaired via uracil DNA glycosylase, is also inhibited by CEBPβ binding. Passage of a replication fork over either a T:G or U:G mismatch, before repair can occur, results in a C-to-T mutation in one of the daughter duplexes. Our study thus provides a plausible mechanism for accumulation of C-to-T human somatic mutations.

Tet2 Inactivation Enhances the Antitumor Activity of Tumor-Infiltrating Lymphocytes

Inactivation of tumor-infiltrating lymphocytes (TIL) is one of the mechanisms mitigating antitumor immunity during tumor onset and progression. Epigenetic abnormalities are regarded as a major culprit contributing to the dysfunction of TILs within tumor microenvironments. In this study, we used a murine model of melanoma to discover that Tet2 inactivation significantly enhances the antitumor activity of TILs with an efficacy comparable to immune checkpoint inhibition imposed by anti-PD-L1 treatment. Single-cell RNA-sequencing analysis suggested that Tet2-deficient TILs exhibit effector-like features. Transcriptomic and ATAC-sequencing analysis showed that Tet2 ablation reshapes chromatin accessibility and favors binding of transcription factors geared toward CD8⁺ T-cell activation. Furthermore, the ETS family of transcription factors contributed to augmented CD8⁺ T-cell function following Tet2 depletion. Overall, our study establishes that Tet2 constitutes one of the epigenetic barriers that account for dysfunction of TILs and that Tet2 inactivation could promote antitumor immunity to suppress tumor growth. SIGNIFICANCE: This study suggests that ablation of TET2⁺ from TILs could promote their antitumor function by reshaping chromatin accessibility for key transcription factors and enhancing the transcription of genes essential for antitumor activity.

Analysis of Il36a induction by C/EBPβ via a half-CRE•C/EBP element in murine macrophages in dependence of its CpG methylation level

Interleukin-36α is a novel member of the IL-1 cytokine family that is highly expressed in epithelial tissues and several myeloid-derived cell types after induction. The transcription factor (TF) C/EBPβ binds specifically to an essential half-CRE•C/EBP motif in the Il36a promoter to induce Il36a expression upon LPS stimulation. C/EBPs regulate gene expression by binding to recognition sequences that can contain 5'-cytosine-phosphate-guanine-3' dinucleotides (CpG), whose methylation can influence TF binding and gene expression. Herein we show that the half-CRE•C/EBP element in the Il36a promoter is differentially methylated in the murine RAW264.7 macrophage cell line and in primary murine macrophages. We demonstrate that C/EBPβ binding to the half-CRE•C/EBP element in the Il36a promoter following LPS stimulation is insensitive to CpG methylation and that methylation of the CpG in the half-CRE•C/EBP element does not alter LPS-induced Il36a promoter activity which correlated with similar Il36a mRNA copy numbers and pro-IL-36α protein amount in both cell types. Taken together, our data indicate that C/EBPβ binding to the half-CRE•C/EBP element and subsequent gene activation occurs independently of the CpG methylation status of the half-CRE•C/EBP motif and underlines the potential of C/EBPs to recognize methylated as well as unmethylated motifs.

Waterborne copper exposure up-regulated lipid deposition through the methylation of GRP78 and PGC1α of grass carp Ctenopharyngodon idella

Early molecular events after the exposure of heavy metals, such as aberrant DNA methylation, suggest that DNA methylation was important in regulating physiological processes for animals and accordingly could be used as environmental biomarkers. In the present study, we found that copper (Cu) exposure increased lipid content and induced the DNA hypermethylation at the whole genome level. Especially, Cu induced hypermethylation of glucose-regulated protein 78 (grp78) and peroxisome proliferator-activated receptor gamma coactivator-1α (pgc1α). CCAAT/enhancer binding protein α (C/EBPα) could bind to the methylated sequence of grp78, whereas C/EBPβ could not bind to the methylated sequence of grp78. These synergistically influenced grp78 expression and increased lipogenesis. In contrast, DNA methylation of PGC1α blocked the specific protein 1 (SP1) binding and interfered mitochondrial function. Moreover, Cu increased reactive oxygen species (ROS) production, activated endoplasmic reticulum (ER) stress and damaged mitochondrial function, and accordingly increased lipid deposition. Notably, we found a new toxicological mechanism for Cu-induced lipid deposition at DNA methylation level. The measurement of DNA methylation facilitated the use of these epigenetic biomarkers for the evaluation of environmental risk.

Sleep deprivation and diet affect human GH gene expression in transgenic mice in vivo

Human (h) growth hormone (GH) production studies are largely limited to effects on secretion. How pituitary hGH gene (hGH-N/GH1) expression is regulated is important in our understanding of the role hGH plays in physiology and disease. Here we assess for the first time the effect of sleep deprivation (SD) and high-fat diet (HFD) on hGH-N expression in vivo using partially humanized 171hGH/CS transgenic (TG) mice, and attempted to elucidate a role for DNA methylation. Activation of hGH-N expression requires interactions between promoter and upstream locus control region (LCR) sequences including pituitary-specific hypersensitive site (HS) I/II. Both SD and diet affect hGH secretion, but the effect of SD on hGH-N expression is unknown. Mice fed a HFD or regular chow diet for 3 days underwent SD (or no SD) for 6 h at Zeitgeber time (ZT) 3. Serum and pituitaries were assessed over 24 h at 6-h intervals beginning at ZT 14. SD and HFD caused significant changes in serum corticosterone and insulin, as well as hGH and circadian clock-related gene RNA levels. No clear association between DNA methylation and the negative effects of SD or diet on hGH RNA levels was observed. However, a correlation with increased methylation at a CpG (cytosine paired with a guanine) in a putative E-box within the hGH LCR HS II was suggested in situ. Methylation at this site also increased BMAL1/CLOCK-related nuclear protein binding in vitro. These observations support an effect of SD on hGH synthesis at the level of gene expression.

Toward a Mechanistic Understanding of DNA Methylation Readout by Transcription Factors

Epigenetic DNA modification impacts gene expression, but the underlying molecular mechanisms are only partly understood. Adding a methyl group to a cytosine base locally modifies the structural features of DNA in multiple ways, which may change the interaction with DNA-binding transcription factors (TFs) and trigger a cascade of downstream molecular events. Cells can be probed using various functional genomics assays, but it is difficult to disentangle the confounded effects of DNA modification on TF binding, chromatin accessibility, intranuclear variation in local TF concentration, and rate of transcription. Here we discuss how high-throughput in vitro profiling of protein-DNA interactions has enabled comprehensive characterization and quantification of the methylation sensitivity of TFs. Despite the limited structural data for DNA containing methylated cytosine, automated analysis of structural information in the Protein Data Bank (PDB) shows how 5-methylcytosine (5mC) can be recognized in various ways by amino acid side chains. We discuss how a context-dependent effect of methylation on DNA groove geometry can affect DNA binding by homeodomain proteins and how principled modeling of ChIP-seq data can overcome the confounding that makes the interpretation of in vivo data challenging. The emerging picture is that epigenetic modifications affect TF binding in a highly context-specific manner, with a direction and effect size that depend critically on their position within the TF binding site and the amino acid sequence of the TF. With this improved mechanistic knowledge, we have come closer to understanding how cells use DNA modification to acquire, retain, and change their identity.

Detection of DNA Modifications by Sequence-Specific Transcription Factors

The establishment, detection, and alteration or elimination of epigenetic DNA modifications are essential to controlling gene expression ranging from bacteria to mammals. The DNA methylations occurring at cytosine and adenine are carried out by SAM-dependent methyltransferases. Successive oxidations of 5-methylcytosine (5mC) by Tet dioxygenases generate 5-hydroxymethyl (5hmC), 5-formyl (5fC), and 5-carboxyl (5caC) derivatives; thus, DNA elements with multiple methylation sites can have a wide range of modification states. In contrast, oxidation of N6-methyladenine by homologs of Escherichia coli AlkB removes the methyl group directly. Both Tet and AlkB enzymes are 2-oxoglutarate- and Fe(II)-dependent dioxygenases. DNA-binding proteins decode the modification status of specific genomic regions. This article centers on two families of sequence-specific transcription factors: bZIP (basic leucine-zipper) proteins, exemplified by the AP-1 and CEBPβ recognition of 5mC; and bHLH (basic helix-loop-helix) proteins, exemplified by MAX and TCF4 recognition of 5caC. We discuss the impact of template strand DNA modification on the activities of DNA and RNA polymerases, and the varied tendencies of modifications to alter base pairing and their interactions with DNA repair enzymes.