ZMYND 11 links histone H 3 . 3 K 36 trimethylation to transcription elongation and tumor suppression
Hong Wen,Yuanyuan Li,Yuanxin Xi,Shiming Jiang,Sabrina Stratton,Danni,Peng,Kaori Tanaka,Yongfeng Ren,Zheng Xia,Jun Wu,Bing Li,C. Michelle,Barton,Wei Li,Haitao Li,Xiaobing Shi
2014-01-01
Abstract:Recognition of modified histones by “reader” proteins plays a critical role in the regulation of chromatin1. H3K36 trimethylation (H3K36me3) is deposited onto the nucleosomes in the transcribed regions following RNA polymerase II (Pol II) elongation. In yeast, this mark in turn recruits epigenetic regulators to reset the chromatin to a relatively repressive state thus suppressing Correspondence and requests for materials should be addressed to xbshi@mdanderson.org, lht@tsinghua.edu.cn, WL1@bcm.edu, or hwen@mdanderson.org. *These authors contribute equally to this work. Author Contributions H.W., H.L. and X.S. conceived this study. H.W. performed the biochemical and cellular studies; Y.L. and H.L. performed structural and calorimetric studies; Y.X. and Z.X. performed bioinformatic analysis; S. J and S. S. performed xenograft study; D.P., K.T. and Y.R. provided technical assistance; J.W. and B. L. performed nucleosome binding experiments; X.S., H.L., H.W., W.L. and M.C.B analyzed data and wrote the paper. Supplementary Information is linked to the online version of the paper at www.nature.com/nature. The authors declare no competing financial interest. Structure data are deposited in PDB with the accession numbers of 4N4G (free BP), 4N4H (BP-H3.1K36me3 complex), and 4N4I (BP-H3.3K36me3 complex). The ChIP-seq and RNA-seq data are deposited in the GEO database with the accession number of GSE48423. NIH Public Access Author Manuscript Nature. Author manuscript; available in PMC 2014 October 10. Published in final edited form as: Nature. 2014 April 10; 508(7495): 263–268. doi:10.1038/nature13045. N IH -P A A uhor M anscript N IH -P A A uhor M anscript N IH -P A A uhor M anscript cryptic transcription2. However, much less is known about the role of H3K36me3 in transcription regulation in mammals. This is further complicated by the transcription-coupled incorporation of the histone variant H3.3 in gene bodies3. Here we show that the candidate tumor suppressor ZMYND11 specifically recognizes H3K36me3 on H3.3 (H3.3K36me3) and regulates Pol II elongation. Structural studies reveal that in addition to the trimethyl-lysine binding by an aromatic cage within the PWWP domain, the H3.3-dependent recognition is mediated by the encapsulation of the H3.3-specific “Ser31” residue in a composite pocket formed by the tandem bromo-PWWP domains of ZMYND11. ChIP-sequencing analyses reveal a genome-wide colocalization of ZMYND11 with H3K36me3 and H3.3 in gene bodies, and its occupancy requires the predeposition of H3.3K36me3. Although ZMYND11 is associated with highly expressed genes, it functions as an unconventional transcription corepressor via modulating Pol II at the elongation stage. ZMYND11 is critical for the repression of a transcriptional program that is essential for tumor cell growth; low expression level of ZMYND11 in breast cancer patients correlates with worse prognosis. Consistently, overexpression of ZMYND11 suppresses cancer cell growth in vitro and tumor formation in mice. Together, this study identifies ZMYND11 as an H3.3-specific reader of H3K36me3 that links the histone variant-mediated transcription elongation control to tumor suppression. ZMYND11 (a.k.a. BS69) is an adenovirus E1A-interacting protein that functions as a corepressor of E1A and cellular transcription factors including the c-Myb and ETS-2 oncoproteins4-6. It contains several histone “reader” modules, including a plant homeodomain (PHD), a bromodomain, and a PWWP domain (Fig. 1a), indicating a role in chromatin regulation. To determine whether these domains can recognize histone modification(s), we utilized a histone peptide array (Supplementary Table 1) and found that the recombinant PHD-bromo-PWWP (PBP) domains specifically bound H3K36me3 peptides (Fig. 1b), which was also confirmed by peptide pull-down assays (Fig. 1c). In addition, ZMYND11 PBP also bound H3K36me3 on bulk histones and the reconstituted H3K36-specific methyl-lysine analog-containing nucleosomes (Extended Data Fig. 1a, b). Several PWWP domains can recognize methylation on histone H3K36 or H3K797-9; therefore we asked whether ZMYND11 PWWP alone was able to bind H3K36me3. The results revealed that neither the isolated PWWP nor the PHD and bromodomain bound H3K36me3, whereas the tandem bromo-PWWP (BP) retained a weak interaction (Fig. 1d and Extended Data Fig. 1c). Isothermal titration calorimetry (ITC) revealed a Kd of 56 μM for PBP and 71 μM for BP to the H3K36me3 peptide, whereas weak or no binding to di-, mono-methylated and the un-methylated H3K36 peptides (Fig. 1e and Extended Data Fig. 1d). Furthermore, deletion of the PWWP, bromo, or PHD domain greatly diminished fulllength ZMYND11 binding to H3K36me3 peptide and chromatin, whereas deletion of the MYND domain had no effect (Extended Data Fig. 1e,f). To gain molecular insights into the reader function of ZMYND11, we solved the crystal structure of the tandem BP domains in its free state and in complex with H3K36me3 peptide at 1.95 Å and 2.0 Å, respectively (Extended Data Table 1 and Supplementary Discussion). Successful crystallization of the BP-H3K36me3 complex was facilitated by alanine mutation of D234 and E236 within the bromodomain to disrupt unfavorable packing contacts without affecting protein folding and histone binding (Extended Data Fig. 2a-e). The tandem BromoWen et al. Page 2 Nature. Author manuscript; available in PMC 2014 October 10. N IH -P A A uhor M anscript N IH -P A A uhor M anscript N IH -P A A uhor M anscript PWWP domains adopt a V-shaped structure imparted by a “kink” around a previously uncharacterized zinc finger (ZnF) motif (Fig. 1f, g). Several structural features imply that ZMYND11 bromodomain is unlikely a histone acetyllysine-binding module (Extended Data Fig. 2f-l). The PWWP domain adopts a five-bladed β-barrel fold with an extended Cterminal α-helix (αP), with a conserved H3K36me3-binding aromatic cage7-9 formed by F291 and W294 of the β1-β2 loop and F310 of the β3-β4 loop (Fig. 1h and Extended Data Fig. 3a-c). Alanine mutations abrogated ZMYND11 binding to the H3K36me3 peptide and chromatin without affecting the overall protein folding (Figs. 1i, 2i and Extended Data Fig. 3d-f). Despite extensive interactions observed in crystals for the histone peptide, the relatively weak binding affinity suggests that an alternative binding mechanism is plausible. In deed, PWWP domain was previously identified as a DNA-binding motif10. Consistent with a recent report11, we found that ZMYND11 BP bound to DNA at low micromolar level (Extended Data Fig. 3g). The higher binding affinity for DNA than that for histone peptides is in agreement with observations for some other PWWP domains12,13. Although not detected in our pulldowns using peptides and oligos, all these features suggest that a synergy between histone and DNA may contribute to ZMYND11 binding at the nucleosomal level (Supplementary Discussion). The H3K36me3 peptides used in the biochemical and structural studies were histone H3.3based peptides (Fig. 2a). The histone variant H3.3 possesses a signature motif “S31⋯A87AIG90” that is distinct from the “A31...S87AVM90” sequence signature of the canonical histone H3.1/23. Our structure revealed that the H3.3-specific “S31” residue and its neighboring T32 were encapsulated at the bottom of the “Bromo-ZnF-PWWP” valley (Fig. 2b), contributing to the ZMYND11-H3 interactions in addition to the binding of K36me3 by the PWWP aromatic cage. Indeed, the binding of ZMYND11 to H3.3K36me3 was much stronger than its binding to the H3.1 counterpart, with 7-8 fold difference in binding affinity (Fig. 2c,d and Extended Data Fig. 4a,b). Furthermore, Flag-tagged histone H3.3, but not H3.1 strongly co-immunoprecipitated ZMYND11 in cells (Extended Data Fig. 4c). Together, these results suggest that the H3.3-specific S31 residue is indispensable for ZMYND11-H3K36me3 interaction. We also solved the structure of ZMYND11 BP-H3.1K36me3 complex at 2.3 Å for comparison (Extended Data Table 1 and Extended Data Fig. 4d). A number of hydrogen bonding interactions around the “A29-G33” segment in the BP-H3.3 complex were lost in the BP-H3.1 complex following the “Ser to Ala” conversion at position 31 (Fig. 2e,f and Extended data Fig. 4e,f). Superimposition of the complex structures revealed a conformational adjustment of the H3.1 “A29-G33” segment (Fig. 2e), indicating a non-ideal encapsulation. As highlighted in the Ligplot diagram14, eight direct hydrogen bonds and two sets of water-mediated hydrogen bonding were detected upon H3.3 variant binding (Fig. 2g), whereas only three direct hydrogen bonds and two sets of water-mediated hydrogen bonding were observed for H3.1 recognition (Fig. 2h). Notably, residues from all three modules, bromo, ZnF and PWWP, contributed to the A29-G33 segment recognition through the formation of a composite pocket, calling attention to the theme of “module integration” for extended reader function15. Unique hydrogen bonding interactions for the H3.3-specific readout were particularly supplied by “R168-A29”, “E251-S31”, “E254-T32” and “N266S31” of the ZMYND11-H3.3 pair along with an adaptive flipping of the R168 side chain Wen et al. Page 3 Nature. Author manuscript; available in PMC 2014 October 10. N IH -P A A uhor M anscript N IH -P A A uhor M anscript N IH -P A A uhor M anscript (Fig. 2e, g), and mutations diminished the bindings (Fig. 2i and Extended Data Fig. 4g). Together, the results strongly support a role of ZMYND11 as an H3.3-specific reader of H3K36me3. We next sought to determine whether ZMYND11 colocalizes with H3K36me3 and the H3.3 variant in the genome by chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq) analysis. We identified 9,759 ZMYND11 peaks distributed over 3,353 genes in human osteosarcoma U2OS cells (Supplementary Table 2). The peaks were highly enriched in intron and exon (herein defined as gene body) but not promoter (Fig. 3a). Strikingly, the average d