Dear Editor,
Histone recognition by reader modules constitutes a major mechanism for epigenetic regulation (
Jenuwein and Allis 2001). BAHD1 (bromo adjacent homology domain containing protein 1) is a vertebrate-specific nuclear protein (Fig. S1) involved in gene silencing by promoting heterochromatin formation. BAHD1 is characteristic with an N-terminal proline-rich region, a nuclear localization signal motif, and a C-terminal bromo adjacent homology (BAH) domain (Fig. 1A). Previous study revealed that BAHD1 could act as a scaffold protein and tether diverse heterochromatin-associated factors including HP1, MBD1, SETDB1, HDAC5, and several transcriptional factors to trigger facultative heterochromatin formation (
Bierne, Tham et al. 2009). Consistent with a “repressive” role, BAHD1 binds to CpG-rich P3 promoter region of
IGF2 (insulin-like growth factor II) then represses
IGF2 and
IGF2 antisense transcription via the recruitment of MBD1 and HDAC5 (
Bierne, Tham et al. 2009). Intriguingly, BAHD1 is also involved in host-pathogen interplay. For example, at early
L. monocytogenes infection state, BAHD1 forms a complex with TRIM28 and HP1 to repress interferon-stimulated genes, including IFNL1, IFNL2, and IFNL3. At specific infection stages,
Listeria secretes a virulence factor, LntA, which could physically interact with BAHD1 to activate interferon (IFN)-stimulated genes (ISGs) (
Lebreton, Lakisic et al. 2011). Despite a repressive role of BAHD1, the molecular mechanism underlying BAHD1 heterochromatin targeting remains largely unexplored.
BAH domain is an evolutionarily conserved motif which is found in several chromatin-associated proteins such as Sir3, ORC1, Rsc2, ZMET2, and DNMT1. BAH is characteristic of a conserved β-sheet core (typically 9-bladed) flanked with function-specific N-, C-, and β6-β7 insertions. Recent structural and functional studies revealed multi-facet roles of BAH domain in chromatin regulation (
Yang and Xu 2013). For example, in yeast, ORC1 BAH could act as a scaffold to mediate ORC1-Sir1 interaction to form a silencing complex; Sir3 BAH could function as a nucleosome-targeting unit to induce heterochromatin formation; moreover, Rsc2 BAH could bind histone H3 and its interaction interface is conserved in a subset of Rsc-like BAH domains (
Chambers, Pearl et al. 2013). Interestingly, in metazoan species, ORC1 BAH domain has acquired a histone methylation reader activity and recognizes H4K20me2 to prompt DNA replication licensing (
Kuo, Song et al. 2012). BAH domains also exist in DNA methyltransferases of mammalian DNMT1 and plant ZMET2. Noteworthily, ZMET2 but not DNMT1 BAH domain displays histone H3K9me2 binding activity and thus directly mediates a cross-talk between histone and DNA methylations (
Du, Johnson et al. 2015). Previous study showed that deletion of the C-terminal BAH domain interfered with co-localization of BAHD1 with H3K27me3 at nuclear foci
in vivo (
Bierne, Tham et al. 2009), suggesting a role of BAH
BAHD1 in histone H3K27me3 recognition. In order to test this hypothesis, we recombinantly expressed BAH
BAHD1 (aa 589–780) with an N-terminal GST tag, and carried out modified histone peptide array screening (Fig. 1B). Many black dots of different intensities were detected in the grid, supporting a histone binding activity of BAH
BAHD1. In-depth data analysis revealed that all the positive hits can be classified into two major categories: the H3K9me2/3-containing and the H3K27me2/3-containing clusters. Among these hits, H3K27me3-containing peptides displayed strongest signal, consistent with a proposed role of BAH
BAHD1 in H3K27me3 recognition.
In order to quantitatively characterize the array data, we performed isothermal titration calorimetry (ITC) under an optimized buffer condition that gives a better melting temperature (Tm
optimized = 39° C vs. Tm
unoptimized = 35.5°C) of BAH
BAHD1 in thermal shift assays (TSA) (Fig. S2). ITC titration revealed a dissociation constant (
KD) of 15.9 μmol/L between BAH
BAHD1 and H3
15-42K27me3 peptide. By contrast, the binding affinity dropped to 2.6 mmol/L in the case of H3
1-15K9me3, and no bindings were observed for H3
1-10R8me2s, H3
1-15K14ac, and H3
28-41K36me3 peptides (Fig. 1C, left), suggesting K27 site-specificity. We next explored the methylation-state preference of H3K27 by BAH
BAHD1. Following the loss of methyl groups, the binding affinity dropped 8-fold for H3K27me2 (
KD = 116 μmol/L), 129-fold for H3K27me1 (
KD = 2.0 mmol/L) and lost completely for unmodified H3K27, suggesting BAH
BAHD1 is a histone H3K27 trimethyllysine-specific reader (Fig. 1C, right). Remarkably, phosphorylation of H3S28 (H3S28ph) dramatically reduced H3K27me3 binding by 42-fold (
KD = 680 μmol/L), suggesting a “methyl-phos” binary switch mechanism of BAH
BAHD1 (
Fischle, Wang et al. 2003).
Both H3K9me3 and H3K27me3 are hallmarks for heterochromatin and gene silencing (
Kim and Kim 2012). To evaluate the
in vivo functional distinction between the two marks in BAHD1 recruitment, we next performed co-localization analysis of BAHD1 with H3K9me3 or H3K27me3 in HeLa cells by immunofluorescence. As shown in Fig. 1D, ectopically expressed EGFP-BAHD1 overlapped nicely with the punctate staining pattern of H3K27me3 as evidenced by the yellow appearance of merged signals. By contrast, little overlapping between EGFP-BAHD1 and H3K9me3 was observed. This result confirms the functional connection of BAHD1 with H3K27me3 but not H3K9me3 at cellular level. H3K9me3 and H3K27me3 share a common “ARKS” sequence motif (Fig. 1E). Additional ITC titrations using H3K27me3 peptides in shorter H3 frames of 19–34 and 21–32 revealed 3.9- and 7.8-fold binding reduction (Fig. 1F), suggesting that distal sequence motifs other than “ARKS” contribute to H3K27me3 recognition and thus discriminate against H3K9me3.
We next performed hydrogen exchange mass spectrometry (HXMS) to map the responsible regions of BAH
BAHD1 for histone H3
15-42K27me3 peptide binding (
Wales and Engen 2006). Three peptide segments spanning “633–653”, “666–691”, and “674–698” of BAHD1 displayed obvious reduced exchange rates in the presence of the H3K27me3 peptide (Fig. 2A), suggesting their involvement in histone recognition. By contrast, other peptides show little or no change in deuterium uptake level upon binding to histone ligand. Structural modelling of BAH
BAHD1 revealed that the three peptide segments cluster together to form a surface that contains an aromatic cage for methyllysine binding (Figs. 2B and S3). This aromatic cage is formed by Y645, W667, and Y669, which are conserved among the BAH domains of mouse ORC1 and plant ZMET2 that are known to recognize H4K20me2 and H3K9me2, respectively (Fig. S4).
To test the importance of the aromatic residues in H3K27me3 readout by BAHBAHD1, we generated single point mutant of Y645A, W667A, Y669A, and performed ITC titration using H315-42K27me3 peptide. As expected, alanine mutation of the aromatic residues disrupted binding between BAHBAHD1 and H3K27me3 peptide (Fig. 2C), supporting a critical role of the aromatic cage in methyllysine recognition. The importance of these aromatic residues was further confirmed by immunofluorescence analysis in HeLa cells. As shown in Fig. 2D and quantification in Fig. 2E, Y645A, W667A, and Y669A mutant EGFP-BAHD1 failed to co-localize with H3K27me3 compared with the wild type protein, thus supporting the functional importance of the aromatic cage for heterochromatin targeting by BAHD1 in vivo.
In sum, combining peptide array screen, quantitative binding, hydrogen exchange MS, and cellular co-localization studies, we established that BAHD1 BAH domain is an H3K27me3-specific reader that discriminates against H3K9me3. H3K27me3 often marks facultative heterochromatin with important functional implications in gene regulation, cell differentiation and development (
Gaydos, Wang et al. 2014). By contrast, histone H3K9me3 represents a hallmark for constitutive heterochromatin and can be recognized by effector proteins such as HP1 to maintain structurally condensed chromatin conformation. Our quantitative ITC assays revealed that the binding affinity of BAH
BAHD1 to H3K27me3 (
KD = 15.9 μmol/L) is more than two orders of magnitude stronger than H3K9me3 (
KD = 2.6 mmol/L). This recognition preference suggests that distal sequences flanking the K9/K27 consensus “ARKS” motif contribute to H3K27me3 recognition. In support, an optimal binding between H3K27me3 and BAH
BAHD1 is achieved in a long frame of 15–42 but not in shorter segments of 19–34 or 21–32. Exact molecular basis underlying H3 distal sequence recognition calls for complex structure determination in the future. Utilizing hydrogen exchange MS, we were able to map key segments within BAH
BAHD1 that are responsible for H3K27me3 binding. Notably, spatial arrangement of these key segments in modelled BAH
BAHD1 structure underscored the role of an aromatic cage consisting of Y645, W667, and Y669 for H3K27me3 readout, a mechanism conserved in many other methyllysine readers (
Patel and Wang 2013). Subsequent ITC titration and immunofluorescence studies comparing wild type and mutant BAH
BAHD1 further confirmed the importance of the aromatic cage in histone H3K27me3 readout both
in vitro and
in vivo.Previously reported H3K27me3 readers include EED WD40 repeats of the PRC2 complex (
Margueron, Justin et al. 2009) and Pc family chromodomain of the PRC1 complex (
Kaustov, Ouyang et al. 2011). Here we characterized human BAHD1 BAH domain as a third class of histone H3K27me3 reader that functions to facilitate BAHD1 heterochromatin targeting and subsequent gene silencing. Moreover, our work revealed that the binding of BAH
BAHD1 to H3K27me3 was markedly disrupted by adjacent H3S28 phosphorylation—a hallmark of the transcriptional response to stress signaling (
Sawicka, Hartl et al. 2014). BAHD1 is overexpressed in peripheral blood mononuclear cells and pancreas, and excels critical function to maintain a repressive state of interferon-stimulated genes and insulin-like growth factors (
Bierne, Tham et al. 2009,
Lebreton, Lakisic et al. 2011). Given the critical role of H3S28ph in signal-induced transcription, the observed “binary switch” between H3K27me3 and H3S28ph may serve as an important mechanism to derepress BAHD1-mediated gene silencing.