The PHD1 finger of KDM5B recognizes unmodified H3K4 during the demethylation of histone H3K4me2/3 by KDM5B

KDM5B is a histone H3K4me2/3 demethylase. The PHD1 domain of KDM5B is critical for demethylation, but the mechanism underlying the action of this domain is unclear. In this paper, we observed that PHD1KDM5B interacts with unmethylated H3K4me0. Our NMR structure of PHD1KDM5B in complex with H3K4me0 revealed that the binding mode is slightly different from that of other reported PHD fingers. The disruption of this interaction by double mutations on the residues in the interface (L325A/D328A) decreases the H3K4me2/3 demethylation activity of KDM5B in cells by approximately 50% and increases the transcriptional repression of tumor suppressor genes by approximately twofold. These findings imply that PHD1KDM5B may help maintain KDM5B at target genes to mediate the demethylation activities of KDM5B.


INTRODUCTION
Covalent histone modifications, notably methylation, are reversible posttranslational modifications that play key roles in chromatin structure, gene transcription and the epigenetic state of the cell (Martin and Zhang, 2005;Mosammaparast and Shi, 2010). Six lysine residues, including histones H3 (i.e., H3K4, H3K9, H3K27, H3K36 and H3K79) and H4 (i.e., H4K20), can be mono-, di-, or trimethylated, and each modification is found in a stereotypical pattern with respect to the coding region of a gene and correlates with a different transcriptional outcome (Zhang and Reinberg, 2001;Margueron et al., 2005). As a general rule, the methylation of H3K9, H3K27 and H4K20 is associated with transcriptional repression, whereas the methylation of H3K4, H3K36 and H3K79 is related to transcriptional activation (Mosammaparast and Zhang and Reinberg, 2001). The methylation of H3K4 is a key regulator for RNA polymerase binding to an active gene (Sims et al., 2003) and of transcription factor binding within promoter elements (Szutorisz et al., 2005). One of the important aspects of H3K4 methylation is how this epigenetic mark is removed, thereby reducing the localization of RNA polymerase to the specific genes. The loss of H3K4 methylation appears to be a key step of differentiation (Mikkelsen et al., 2007). To date, two distinct classes of histone demethylases have been characterized. The first class includes two members, represented by LSD1 (lysine-specific demethylase 1, also known as KDM1A, which demethylates H3K4me1/2) and LSD2 (also called KDM1B or AOF1, which demethylates H3K4me2).
Both members of this class use FAD as a cofactor through an amine-oxidation reaction to remove the lysine methyl group of H3K4me1/2 (Shi et al., 2004;Ciccone et al., 2009). The other class contains a JmjC (i.e., Jumonji C) domain to catalyze histone lysine demethylation assisted by two cofactors: Fe 2+ ion and α-ketoglutarate (α-KG) Cloos et al., 2006;Tsukada et al., 2006;Whetstine et al., 2006, Shi and. Based on the sequence homology in the JmjC domain and the overall architecture of the associated motifs, JmjC domain-containing proteins have been classified into seven groups: JHDM1, PHF2/8, JARID, JHDM3/JMJD2, UTX/UTY, JHDM2 and JmjC domain only . Structural investigations of the members of this class, including JHDM1A, JHDM1D, JMJD2A and PHF8, have already been performed for their apo forms or in complex with the H3 peptide and α-KG substrates Chen et al., 2007;Couture et al., 2007;Ng et al., 2007;Horton et al., 2010;Yang et al., 2010). However, no structure has been published for the members of the JARID1 sub-group that can specifically remove methyl groups of di-or tri-methylated H3K4.
The members of the JARID1 subgroup are highly conserved from yeast to humans and contain a similar motif architecture, including JmjN, ARID (i.e., AT-rich interactive domain), JmjC, Zf-C 5 HC 2 (i.e., zinc finger) and two or three PHD domains (denoted PHD1, PHD2 and PHD3 from the N-terminus to the C-terminus). A total of four members are found in mammals (Fig. 1A): JARID1A (also called RBP2 or KDM5A), JARID1B (also named PLU-1 or KDM5B), JARID1C (i.e., SMCX or KDM5C) and JARID1D (also known as SMCY orKDM5D). These members have been identified to be H3K4me2/3 demethylases (Christensen et al., 2007;Iwase et al., 2007;Klose et al., 2007;Lee et al., 2007;Tahiliani et al., 2007;Yamane et al., 2007). All of these proteins are key transcriptional co-repressors because they can remove the transcription-activating marker H3K4me3. KDM5B is involved in transcriptional repression and breast cancer cell proliferation ; thus, mechanistic studies on KDM5B demethylation are necessary and useful to understand the development of breast cancer. Notably, the deletion of the N-terminal PHD1 finger (i.e., PHD1 KDM5B ) of KDM5B results in the loss of enzymatic demethylase activity, implying that PHD1 KDM5B is involved in H3K4me2/3 demethylation . This observation is consistent with the fact that the N-terminal PHD1 finger of Lid (i.e., PHD1 Lid ), a homologue of KDM5B in Drosophila, is also required for the demethylase activity of H3K4me3, whereas the PHD2 and PHD3 of Lid are not . However, the detailed mechanism underlying the function of PHD1 KDM5B in the demethylation process remains unclear.
Recently, the NAD + -dependent PARylation on KDM5B by poly (ADP-ribose) polymerase (PARP-1) was reported to be able to regulate chromatin structure and transcription through a KDM5B-dependent pathway (Krishnakumar and Kraus, 2010). The demethylation inhibited by PARylation on KDM5B was confirmed by an in vitro H3 binding assay. Previously, PHD1 Lid was reported to bind with unmodified H3K4me0 . These observations indicate that KDM5B may bind to the H3 peptide, most likely to unmodified H3K4. Thus, in this study, we first tested the interactions of the unmodified H3K4 peptide with full-length KDM5B or its truncated variants through a biotin-labeled peptide binding assay. By sequentially depleting different regions in KDM5B-N, we found that PHD1 in KDM5B (i.e., PHD1 KDM5B ) can specifically bind to the unmodified histone H3. To probe the structural basis for this interaction, we determined a solution structure of PHD1 KDM5B in complex with the unmodified H3K4 peptide. Through structural and biochemical data, we provide insights into the function of PHD1 KDM5B in KDM5B-regulated demethylation and tumor-suppressor gene transcription.

PHD1 KDM5B specifically interacts with the unmodified H3K4me0 peptide
To investigate the function of PHD1 KDM5B in the demethylation of H3K4me2/3 by KDM5B (as shown in Fig. 1B), the binding affinities of five recombinant KDM5B variants, namely KDM5B-N (N-terminal KDM5B), KDM5B-PC (only containing PHD1 and JmjC domains in KDM5B-N), KDM5B-N-△P (without the PHD1 domain in KDM5B-N), KDM5B-N-△C (without the JmjC domain in KDM5B), and KDM5B-N-△ZF (without the Zf-C 5 HC 2 domain in KDM5B-N), were tested with a biotinlabeled unmodified H3 peptide (Fig. 1C). The results indicate that KDM5B-N interacts with the unmodified H3 peptide. Similar to the KDM5B-N variant, the KDM5B-PC, KDM5B-N-△C and KDM5B-N-△ZF variants bind to the unmodified H3K4 peptide, implying that the JmjC domain, the Zf-C 5 HC 2 domain, and the JmjN-ARID domain are not involved in the interaction with the unmodified H3K4 peptide. However, the deletion of the N-terminal PHD1 KDM5B domain in KDM5B-N (i.e., KDM5B-N-△P) significantly impaired the interaction between KDM5B-N and the H3K4me0 peptide. This result suggested that PHD1 KDM5B primarily contributes to the binding of KDM5B to the H3K4me0 peptide. The glutathione S-transferase (GST) tag-fused PHD1 KDM5B (306-360 aa, similarly hereinafter) protein binds to the unmodified H3 tail (Fig. 1D), further supporting the function of PHD1 KDM5B in the specific binding to the unmodified H3. The methylation of H3R2 and H3K4 abolishes or weakens the binding affinity, suggesting that these two amino acids may be involved in the binding. In addition, the methylation of H3K9 does not notably inhibit this binding, which indicates that H3K9 does not participate in the binding. Moreover, an ITC analysis obtained the dissociation constants (K D ) of 6.4 ± 0.6 μmol/L for PHD1 KDM5B interacting with the unmodified H3K4 peptide, 25.6 ± 3.8 μmol/L for PHD1 KDM5B interacting with mono-methylated H3K4me1, 80.0 ± 7.9 μmol/L for PHD1 KDM5B interacting with dimethylated H3K4me2 and 103.7 ± 11.2 μmol/L for PHD1 KDM5B interacting with tri-methylated H3K4me3 (Fig. 1E). These K D s revealed that H3K4 is highly involved in the binding and that the methylation of H3K4 inhibits this binding.
Solution structure of PHD1 KDM5B in complex with H3K4me0 To understand the interaction of PHD1 KDM5B with unmodified H3K4, we initially attempted to crystallize PHD1 KDM5B in its free form and in complex with the H3K4me0 peptide. However, we only obtained the X-ray structure of the free PHD1 KDM5B with a high resolution of 1.65 Å (Fig. S1) . Thus, we determined the solution structures of free PHD1 KDM5B and of PHD1 KDM5B in complex with the H3K4 peptide (1-10 aa) using multidimensional heteronuclear NMR spectroscopy (Fig. 2). To probe whether the electronic properties of histidines in solution are similar to those in the crystal state (two histidines, H335 and H344, are in the amino acid sequence of PHD1 KDM5B ; the H335 residue was suggested to ligate with a zinc ion in the X-ray crystal structure of free PHD1 KDM5B ) before structural determination, we analyzed the 1 H-15 N LR-HSQC spectra of PHD1 KDM5B in the free and bound states (Fig. S2). This LR-HSQC experiment correlates the carbon-bound protons of the histidine rings with the imidazole nitrogen atoms and can unambiguously establish the tautomeric and protonation states of histidines in proteins (Pelton et al., 1993;Drohat et al., 1999 H344 in PHD1 KDM5B , revealing the strong chelation of H335 to an electropositive Zn 2+ . This observation is identical to that observed in the X-ray structure of the free PHD1 KDM5B . The NMR solution structure of free PHD1 KDM5B (Fig. 2B) was then determined using the program XPLOR-NIH (Kuszewski and Clore, 2000) with 670 NOEs, 10 hydrogen bonds and 106 dihedral angle constraints that were meaningful and acceptable (Table 1). The resulting structure was similar to that obtained by X-ray techniques with a backboneatom RMSD of 0.48 Å, which was determined by superimposing the backbone C α atoms in the secondary structural region (Fig. S1). Thus, in this paper, we only use the NMR structure of the free PHD1 KDM5B for comparison with that of the complex. In contrast, the solution structure of PHD1 KDM5B in complex with an unmodified H3K4me0 peptide (1-10 aa) was obtained through 1263 NOEs (including 150 observed intermolecular NOEs), 17 hydrogen bonds and 106 dihedral angle constraints. The PHD1 KDM5B structures in both free and bound states adopt a 'cross-braced' topology of zinc ion coordinated residues that was identical to that of all structurally characterized PHD fingers. These states are well defined by NMR data; the RMSD of the free PHD1 KDM5B and the PHD1 KDM5B in complex with the H3 peptide were 1.18 ± 0.17 Å and 0.63 ± 0.14 Å for the heavy atoms in the secondary structural regions of the 20 lowest energy structures, respectively. By superimposing the backbone C α atoms in the secondary structural region, the free and bound solution PHD1 KDM5B structures have a backbone-atom RMSD value of 1.13 Å (Fig. S1). This RMSD indicates that the H3K4me0 peptide binding does not induce major conformational changes to the backbone of PHD1 KDM5B .
In the structure of the PHD1 KDM5B -H3K4me0 complex, the H3 peptide binds to the surface of the PHD1 KDM5B as an anti-parallel β-sheet (residues 2-3 aa). This H3 β-sheet is shorter in this complex than in the complex of the PHD of human autoimmune regulator, (PHD AIRE ) with unmethylated H3K4 (PDB code: 2KFT) (Chakravarty et al., 2009) and in the complex of the PHD finger of the BHC80 protein in the LSD1 co-repressor (PHD BHC80 ) with unmethylated H3K4 (PDB code: 2PUY)   (Fig. 2D). In the current structure of the complex, the backbone atoms of residues H3R2 and H3K4 in the unmethylated H3 peptide form hydrogen bonds with the backbone carbonyl oxygen and nitrogen atoms of L324 and L326 in PHD1 KDM5B (Fig. 2F), respectively. The cognate PHD1 KDM5B only contacts the first six residues of the H3 peptide, whereas H3K9 is completely exposed to solvent by extending its side chain away from the complex (Fig. 2E). This finding supports the concept that the methylation of H3K9 has no effect on the binding. The conformation of the side chain of H3R8 is also flexible. In the 20 final NMR structures for the complex, H3R8 does not contact the protein. This coincides with the fact that no intermolecular NOEs were observed between the H3R8 side chain and the protein.

Analysis of the interaction between PHD1 KDM5B and the unmodified H3K4me0 peptide
According to the structural information provided by the structure of the complex of PHD1 KDM5B with the unmodified H3K4 peptide, the H3K4me0 specificity of PHD1 KDM5B is determined through the recognition of the residues in the H3 amino terminus, including H3A1, H3R2, H3T3, H3K4 and H3T6 ( Fig. 2F-K). The H3A1 methyl group is anchored by intermolecular hydrogen bonds with the backbone carbonyl oxygen atoms of residues P347 and G349 (Fig. 2F, in which the H3A1 backbone nitrogen supplies two hydrogen bonds) and by nonpolar hydrophobic interactions between the H3A1 methyl group and the side chains of L325, P347 and W351 in PHD1 KDM5B (Fig. 2G). These side chains are similar to those observed in the PHD AIRE -H3K4me0 complex and other complexes (Chakravarty et al., 2009;Lan et al., 2007, Pena et al., 2006. These observations support our site-directed mutation studies for both the PHD1 KDM5B protein and the H3K4me0 peptide (Table S1). On the one hand, both the W351A and L325A mutants of PHD1 KDM5B have nondetectable binding affinities with the H3K4me0 peptide. On the other hand, removing the methyl group from H3A1 by changing alanine to glycine results in an approximately 50-fold decrease in the binding affinity of the unmodified H3 s peptide to PHD1 KDM5B (K D H3K4-PHD1 = 6.4 ± 0.6 μmol/L and K D H3 A1G-PHD1 = 304.9 ± 11.7 μmol/L). In the 1 H-15 N HSQC spectra, compared with that of free PHD1 KDM5B , the binding to wild-type H3A1 (i.e., H3K4me0) produces a large shift in most of the cross-peaks of PHD1 KDM5B , whereas the binding of the H3G1 variant does not induce this shift. This result indicates the important roles of the H3A1 methyl group in the interaction. One of the positively charged side-chain nitrogen atoms (Nη) of H3R2 forms a salt bridge with one of the negatively charged oxygen atoms in the side-chain of D328 (Fig. 2H). D328 mutations in PHD1 KDM5B (the D328A mutant) or R2 mutations in the H3K4me0 peptide (H3 R2A or R2E mutants) remove the negative or positive charges in their side chains; therefore, the binding affinities were significantly reduced (K D H3A1-PHD1 D328A = 182.1 ± 21.9 μmol/L, the binding affinity decreased by 30-fold; K D H3 R2A-PHD1 = 370.4 ± 52.1 μmol/L, the binding affinity decreased by 60-fold; and the K D value for the interaction of PHD1 KDM5B to the H3 R2E mutant was too small to detect) ( Table S1) compared with that of the wild-type PHD1 KDM5B or the unmodified H3K4 peptide. These results were consistent with the observations mentioned for the above-described structures. In addition, the distance between the Cβ atom of L326 and the Cβ atom of H3R2 is 4.0 Å, indicating a hydrophobic interaction between the side chains of L326 and H3R2 (Fig. 2H).
The side-chain -OH group of H3T3 forms two additional hydrogen bonds with one of the side-chain carbonyl oxygen atoms and the backbone oxygen of E321 (Fig. 2I), further stabilizing the interaction between the H3 peptide and the ± 2.2 μmol/L). The distances between the H3T6-OH group and the backbone carboxyl oxygen of D308 and Y310 are less than 4 Å (Fig. 2K), indicating that weak hydrogen bonds might form between these atoms. The presence of these weak hydrogen bonds is supported by the measured K D values of the mutants (K D H3 T6V-PHD1 = 61.3 ± 1.6 μmol/L, the binding affinity decreased by approximately tenfold compared with that of the wild-type protein, which exhibits a K D H3K4-PHD1 of 6.4 ± 0.6 μmol/L).
Moreover, the side-chain NH 3 + group of H3K4 forms a rigid hydrogen bond with the side-chain carbonyl oxygen of D308 (Fig. 2J), which supports the results of the analysis of the structure of the complex and of the mutation studies (K D H3 K4E-PHD1 is non-detectable, whereas K D H3K4-PHD1 D308A was measured as 41.7 ± 2.2 μmol/L). In addition to interacting with D308, the side chain of H3K4 also displays weak hydrophobic interactions with the aromatic ring of Y310 and the methyl groups of L326 (Fig. 2J). The distances between the Cγ or Cδ atom of the H3K4 side chain and the aromatic Cγ atom of Y310 or the methyl group of L326 in PHD1 KDM5B are approximately 4.5 Å. The mutation of Y310 to F310 does not change the binding affinities of PHD1 KDM5B to the unmodified H3 peptide (K D H3K4-PHD1 = 6.4 ± 0.6 μmol/L and K D H3K4-PHD1 Y310F = 4.3 ± 0.2 μmol/L), suggesting that the-OH group in the side chain of residue Y310 may not be involved in the interaction. When Y310 is replaced by A310, the binding affinity between PHD1 KDM5B and the H3 peptide decreased by approximately fivefold (K D H3K4-PHD1 Y310A = 28.1 ± 1.3 μmol/L) compared with the binding affinity of the wild-type PHD1 KDM5B to the H3 peptide. The PHD1 KDM5B L326A mutant has an approximately two-fold weaker binding affinity (K D H3K4-PHD1 L326A = 14.1 ± 0.8 μmol/L) to the unmodified H3 peptide than the wild-type. These observations imply that the hydrophobic interactions between the side chains of L325 and Y310 in PHD1 KDM5B and H3K4 contribute less to the binding than the hydrogen bond between H3K4 and D308 in PHD1 KDM5B . The methylation of H3K4 may weaken these hydrophobic interactions because minimal space exists among residues L325 and Y310 in PHD1 KDM5B and H3K4. This may explain the reduction in the binding affinities of PHD1 KDM5B to H3K4me1 (by ∼4-fold), H3K4me2 (by >10-fold) and H3K4me3 (by >15fold) compared with the H3K4me0 peptide. Therefore, molecular recognition of the unmodified lysine 4 primarily occurs through hydrogen bonding to the unmodified epsilon amino group and steric elusion of the methyl groups on H3K4me3/2. In addition, residues D308, Y310 and L326 are conserved among the members in the KDM5 PHD1 family (Fig. 3B) with the exception of residues E323 in KDM5C and F331 in KDM5D. Therefore, our structure may reveal that the PHD1 domain in other members of the KDM5 family functions as a specific reader of unmodified H3K4, although structural models of the PHD1 KDM5A, KDM5C and KDM5D in complex with unmodified H3K4 are not available.

Recognition of H3K4me0 by PHD1 KDM5B affects KDM5B demethylase activity
The binding sites of unmodified H3K4 on PHD1 KDM5B were first determined by an NMR titration binding assay ( Fig. 2A  and 2B). The majority of residues with chemical shifts larger than 0.3 ppm are involved in the interaction; this result is consistent with the results from the ITC assay (Table S1) and with the structural information of the complex. In the solution structure of PHD1 KDM5B in complex with H3K4, residues D308, L325, D328 and W351 are involved in the interaction between PHD1 KDM5B and H3K4 through hydrogen bonds, hydrophobic interactions, or salt bridges (Fig. 2). The singlesite mutation of D308, L325, D328 or W351 to alanine disrupts the interaction between PHD1 KDM5B and the H3K4me0 peptide. Based on these structural and biochemical analyses, we then designed three full-length KDM5B mutants (D308A, L325A/D328A and W351A) and tested whether the in vivo H3K4me2/3 demethylation by KDM5B is affected by disrupting the interaction between PHD1 KDM5B and H3K4me0 through an immunofluorescence staining assay.
Compared with the wild-type (WT) KDM5B (99% of H3K4me3 was demethylated), the D308A mutant exhibited 79% demethylase activity on H3K4me3 (decreased by approximately 20%) (Fig. 4A). The mutants that completely lose their ability to binding to the unmodified H3K4 exhibited decreases in their demethylase activities to 72% (W351A) and 56% (L325A/D328A). For H3K4me2 demethylation (Fig. 4B), the wild-type KDM5B displays 93% demethylase activity, whereas the D308A, W351A and L325A/D328A mutants demonstrate 76%, 66% and 55% demethylase activities, respectively. In both cases, the L325A/D328A mutant has a higher effect on the demethylase activities of KDM5B. These data indicate that the KDM5B demethylase activity was partially affected by (but not dependent on) the interaction between the N-terminal PHD1 finger and the unmodified N-terminal H3K4 tail.
Recognition H3K4me0 by PHD1 KDM5B affects repression of tumor suppressor genes KDM5B has been reported to function as a transcriptional repressor, and the knockdown of KDM5B increases the expression level of several tumor suppressor genes, including 14-3-3s, BRCA1, CAV1 and HOXA5 (Tan et al., 2003;Yamane et al., 2007). To illustrate the functional importance of the binding between PHD1 KDM5B and the unmodified H3K4, we tested whether the mutants could affect the expression of tumor suppressor genes compared with wild-type KDM5B. As shown in Fig. 4C, the overexpression of WT KDM5B decreased the expression of tumor suppressor genes. Compared with WT KDM5B, the overexpression of the D308A, W351A and L325A/D328A mutants up-regulated the expression of the four genes to a lesser but still significant extent. Among the mutants, the L325A/D328A mutant increased the gene expression by more than twofold. Therefore, the gene expression results agree with the binding affinities of PHD1 KDM5B and its mutants with the H3 peptide. The overexpression of the mutants is expected to produce fewer unmodified H3K4, thereby increasing the gene repression level. In combination, these findings support the model in which PHD1 KDM5B specifically binds to the unmodified histone H3. Additionally, the H3 histone remains unmodified. This process is correlated with KDM5B demethylase activity and KDM5B-mediated gene repression.

Comparison between PHD1 KDM5B and other PHD fingers
The reported structure of PHD BHC80 in complex with unmodified H3K4 (PDB code: 2PUY) demonstrated that H3K4 residue forms hydrogen bonds with the side chains of D489 and E488 (Fig. 3A)  . Moreover, the Cβ atom of H487 restricts the interaction of PHD BHC80 with methylated H3K4, which means that the H3K4me1 and H3K4me2 peptides cannot bind to PHD BHC80 . However, in the PHD1 KDM5B -H3K4me0 complex, H3K4me0 recognition is enhanced by hydrophobic interactions of the side-chains of H3K4 with the aromatic ring of Y310 and the methyl group of L326. Therefore, the binding mode of PHD1 KDM5B -H3K4me0 differs slightly from that of PHD BHC80 . Recognition of H3K4me0  in PHD1 KDM5B , PHD BHC80 , PHD ARIE (Chakravarty et al., 2009), the PHD finger in UHRF1 (i.e., PHD UHRF1 , PDB code: 2LGG) Wang et al., 2011;Lallous et al., 2011;Rajakumara et al., 2011)    the PHD finger and the K4 side chain (Fig. 3). These differences result from the action of the conserved residues of glycine (G) in PHD BPTF and PHD ING2 or cysteine (C) in PHD UHRF1 and PHD1 KDM5B . In PHD BHC80 , this space is occupied by residue M502. In comparison, methylated H3K4me3 recognition was found to occur primarily through the interaction between the positively charged H3K4me3 and the conserved aromatic residues, such as W32 in PHD BPTF , W238 in PHD ING2 and W891 in PHD TAF3 , as shown in Fig. 3. Other residues, including Y10 and Y17 in PHD BPTF , Y315 and M226 in PHD ING2 , and W868 and M882 in PHD ING2 , associate through hydrophobic interactions or aromatic Π-cation ion interactions to form the binding cage. Thus, similar to PHD BHC80 and other PHD fingers (Fig. 3A), PHD1 KDM5B does not use the aromatic cage to specifically identify H3K4, which is present in structurally characterized methylated lysine PHD fingers, such as the PHD fingers in BPTF (i.e., PHD BPTF , PDB code: 2F6J), ING2 (i.e., PHD ING2 , PDB code: 2G6Q) and TAF3 (i.e., PHD TAF3 , PDB code: 2K17) (Fig. 3A) (Pena et al., 2006;Li et al., 2006). These PHD fingers adopt similar folds, engage the H3 peptide as an anti-parallel β-sheet on the surface and recognize the H3 N-amine and H3A1 side chain. H3R2 is buried in a pocket in PHD BPTF , PHD ING2 , PHD UHRF1 Wang et al., 2011;Lallous et al., 2011;Rajakumara et al., 2011;Xie et al., 2012) and PHD2 DPF3b (Zeng et al., 2010), bound in tiny surface grooves in PHD AIRE , PHD TAF3 and PHD1 KDM5B , and not contacted in PHD BHC80 .
The binding mode of PHD1 KDM5B to H3K4me0 is also distinguished from the caging of the di-and tri-methyl lysine by aromatic residues, as identified in the polycomb (Pc) and heterochromatin protein I (HP1) chromodomain (Fischle et al., 2003;Nielsen et al., 2002;Jacobs and Khorasanizadeh, 2002;Min et al., 2003;Tan et al., 2003). The methylation of lysine 9 in histone H3 is recognized by HP1; this methylation directs the binding of other proteins to control chromatin structure and gene expression. The structures of the complex between the Drosophila HP1 chromodomain and the histone H3 tail with a di-or tri-methylated K9 display histone tail inserts as a β strand, completing the β-sandwich architecture of the chromodomain. The methylated lysine is caged by the side chains of the aromatic residues Y21, W42 and F45, whereas adjacent residues form discerning contacts with one face of the chromodomain. The structure of the Pc chromodomain in complex with a H3 peptide bearing trimethylated K27 (Fischle et al., 2003) demonstrates that the methylated H3K27 is caged by four aromatic residues (i.e., Y26, W47, W50 and Y54) preceding the ARKS motif.

Proposed biological function of PHD1 KDM5B in demethylation by KDM5B
The binding of the PHD1 KDM5B domain to the H3K4me0 peptide, which is a demethylation product of the full-length KDM5B, indicates that PHD1 KDM5B may function downstream of KDM5B demethylase activity. This observation is similar to that found for BHC80 in LSD1-mediated H3K4me2 demethylation and repression . The regulation of histone methylation is highly dynamic and involves the actions of both a methyltransferase and demethylase on identical target promoters. Therefore, we suggest that PHD1 KDM5B may be important in maintaining KDM5B at the target promoters and preventing the re-methylation of H3K4. The downstream effector of PHD1 KDM5B is therefore required for the KDM5Binduced demethylation of H3K4me2/3. This was confirmed by our binding assay results, which showed that the removal of PHD1 KDM5B resulted in decreased binding of KDM5B-N to the K4 unmethylated histone H3 peptide in vitro (Fig. 1C). Our findings indicate that PHD1 KDM5B is important for KDM5B association with its reaction product H3K4me0 after demethylation. Thus, PHD1 KDM5B operates both as a reader and a protector of unmethylated H3K4.
Similar to PHD1 KDM5B , PHD1 Lid has also been suggested to interact with H3K4me0 . PHD1 Lid was proposed to bind to non-DNA elements, such as local chromatin environments, during H3K4me2/3 demethylation by Lid. The C-terminal PHD3 region of Lid (i.e., PHD3 Lid ) was observed to specifically bind to H3K4me2/3 through the interaction of aromatic residues in the PHD finger with the positively charged methylated H3K4. By aligning the amino acid sequences of PHD3 KDM5B , PHD3 Lid and other PHD fingers that recognize the H3K4me2/3 peptide (Fig. 5), we found that the residues most likely interacting with the H3K4me2/3 site are highly conserved. Residues W1781 in PHD3 Lid and W1512 in PHD3 KDM5B are conserved corresponding to W32 in PHD BPTF , W238 in PHD ING2 , and W891 in PHD TAF3 . We thus suggested that PHD3 KDM5B might also specifically bind to methylated H3K4me2/3 and that the W1502 residue in PHD3 KDM5B may have a biological function similar to that of residues W1771 in PHD3 Lid , Y17 in PHD BPTF , M226 in PHD ING2 and M882 in PHD TAF3 . In mammalian cells, c-Myc prefers to bind to E-boxes located within a chromatin context that contain highly di-and tri-methylated nucleosomal histone H3K4 (Guccione et al., 2006). However, the mechanism through which Myc recognizes the chromatin landscape remains unclear. Here, we propose that KDM5B may utilize its H3K4me2/3-binding C-terminal PHD3 finger to tether Myc to its preferred chromatin context. This process may be enhanced by the interaction between PHD1 KDM5B and the unmethylated H3K4me0 N-terminal tail, thereby permitting the selection of biologically important E boxes. Further experiments are required to more precisely define the roles of KDM5B PHD fingers in cell growth.
In summary, we identified a specific interaction between PHD1 KDM5B and the unmodified H3K4 peptide. We further provided structural insights into the binding. The specific recognition of unmodified H3K4 by the PHD1 domain of KDM5B is important for the KDM5B histone demethylase activity in cells and for the transcriptional repression of tumor suppressor genes.

PHD1 KDM5B preparation for crystallization and NMR experiments
The details of the preparation, purification, crystallization and determination of the PHD KDM5B (residues 306-360 aa) X-ray structure were previously described . The modified pGEX-6p-1 vectors expressed proteins with the N-terminal GST tag, which is removable by cleavage with a 3C protease, enabling the use of non-tagged proteins in our studies. Only the X-ray structure of free PHD1 KDM5B was obtained.
A similar PHD KDM5B construct was used for the sample required in the NMR experiments. 13 C-and 15 N-labeled PHD KDM5B were prepared in M9 medium. Site-directed mutagenesis was performed using a QuikChange site-directed mutagenesis kit (Stratagene Inc. La Jolla, California, U.S.A). All of the DNA constructs were sequenced, and the molecular weights of the recombinant proteins were verified by mass spectrometry (MALDI).

GST pull-down and biotin pull down assays
For the biotin pull-down assay, biotin-labeled H3 peptides were used to pull down the PHD1 KDM5B protein in order to determine the binding affinities. Briefly, 0.5 μg of the peptide was mixed with the protein at an identical molar ratio, and the mixture was incubated with streptavidin beads at 4°C overnight. After washing five times, the beads were boiled in SDS loading buffer and separated on a SDS-PAGE gel.
To investigate whether the PHD1 KDM5B finger and its mutants interact with unmodified or methylated H3K4 or H3K4mutants, the binding affinities of PHD1 KDM5B or its mutants to H3 peptides were studied. An ITC-200 microcalorimeter (GE Healthcare) was used with a buffer containing 20 mmol/L Tris, 150 mmol/L NaCl, and 1% Triton X-100, pH 7.5 at 25°C. The reference titration of small molecules in the buffer was subtracted from the experimental data, and the data were fitted using the Origin 7.0 (OriginLab Corporation) software. The results are summarized in Table S1.

NMR spectroscopy and analysis
The NMR samples contained 1.5 mmol/L uniformly 13 C/ 15 N-labelled PHD KDM5B and the unlabeled H3 peptide H3K4me0 in complex at a PHD KDM5B -to-H3-peptide molar ratio of 1:6 in NMR buffer (20 mmol/L Na 2 HPO 3 , 100 mmol/L NaCl, 0.01% NaN 3 , pH 7.4 and 10% D 2 O). All of the NMR experiments were performed at 20°C on a Varian Unity Inova 600 NMR spectrometer equipped with a triple resonances cryoprobe and pulsed field gradients. The standard suite of experiments for assigning the 1 H, 13 C and 15 N backbone, determining the side-chain chemical shifts of PHD KDM5B in complex with the H3 peptide and collecting the NOE-based distance restraints were measured (Bax and Grzesiek, 1993;Clore and Gronenborn, 1998), and these included 2D 13 C-edited HSQC and 15 N-edited HSQC; 3D HNCA, HNCO, HN(CO)CA, HNCACB, CBCA(CO)NH, 15 N-resolved HSQC-TOCSY and HCCH-TOCSY in both aliphatic and aromatic regions; 15 N-resolved HSQC-NOESY; 13 C-resolved HSQC-NOESY for both aliphatic and aromatic resonances; and 2D hbcbcgcdceheA and hbcbcgcdhdA spectra for the correlation of Cβ and Hδ or Hε in the aromatic ring that is used for aromatic proton assignment (Yamazaki et al., 1993). The proton NMR signals of the bound H3 peptides were assigned by analyzing the 2D 13 C-filtered, 15 N-filtered and J-resolved NOESY and TOCSY spectra recorded for the 13 C-and 15 N-labeled protein with the unlabeled H3 peptide H3K4me0 and the 2D 1 H-1 H COSY, NOESY and TOCSY spectra  Figure 5. Sequence alignment of PHD fingers (binding to H3K4me3/2) of KDM5B (i.e., PHD3 KDM5B ), Lid (i.e., PHD3 Lid ), BPTF, ING2 and TAF3. The zinc-binding residues, H3A1-binding residues, H3R2-binding residues and H3K4me3/2-binding residues are highlighted in grey (and star on the top of KDM5B), blue, red and brown, respectively. recorded for the unlabeled free H3 peptides in the NMR buffer mentioned above, respectively. The intermolecular NOEs between the labeled protein and the unlabeled H3 peptides were obtained by analyzing the 3D 13 C-F1 edited and 13 C/ 15 N-F3 filtered NOESY spectra. The spectra were processed with the NMRPipe program (Delaglio et al., 1995) and analyzed using Sparky 3 (http://www.cgl. ucsf.edu/home/sparky/).

Determining the NMR structure
The calculations were performed using a standard simulated annealing protocol implemented in the XPLOR-2.19 program (NIH version) (Kuszewski and Clore, 2000). The inter-proton distance restraints derived from the NOE intensities were grouped into three distance ranges, namely 1.8-2.9 Å, 1.8-3.5 Å and 1.8-6.0 Å, which corresponds to strong, medium and weak NOEs, respectively. The dihedral angles phi and psi were derived from the backbone chemical shifts (HN, HA, CO and CA) using the program TALOS (Cornilescu et al., 1999). Slow-exchanging amide protons identified in the 2D 15 N-1 H HSQC spectra recorded after the H 2 O buffer was exchanged for a D 2 O buffer were used in the structure calculated with the NOE distance restraints to generate hydrogen bonds for the final structure calculation as previously described in the literature (Chakravarty et al., 2009). Constraints between the protein ligands and the zinc ion were added using a previously reported procedure (Neuhaus et al., 1992;Cao et al., 2003). A total of ten iterations were performed (50 structures in the initial eight iterations). In total, 100 structures were computed during the last two iterations, and the 20 conformers with the lowest energy were used to represent the 3D structures. The conformers of these two bundles (free PHD1 KDM5B and PHD1 KDM5B in complex with the H3K4me0 peptide) do not violate the following constraints: NOE > 0.3 Å and dihedral angle >3°. The entire structure statistics were evaluated with PROCHECK (Laskowski et al., 1993) and PROCHECK-NMR (Laskowski et al., 1996) and are summarized in Table 1. All of the structure figures were generated using the PyMOL (http://pymol.org/) and MOLMOL programs (Koradi et al., 1996).

Immunofluorescence staining
The 293T cells were transfected with pcDNA3-KDM5B-Myc-His WT and mutants using Lipofectamine 2000. After 48 h, the cells were fixed with 4% paraformaldehyde for 10 min and permeabilized with 0.5% Triton in PBS for 15 min. After blocking, the cells were incubated with a primary antibody (Millipore, anti-H3K4me2, 1:200 dilution), and washed three times with PBS, and incubated with fluorescence-labeled secondary antibody for 1 h (Molecular Probes, Alexa Fluor 555 goat anti-rabbit, 1:250 dilution). After extensive rinsing with PBS, cover slips were mounted with an antifade reagent and DAPI (Molecular Probes) and examined on an Olympus IX51 microscope.

Real-time qPCR assay
The pcDNA3-KDM5B-Myc-His WT and mutants were transfected into 293T cells with Lipofectamine 2000. Forty-eight hours after transfection, the total RNA was extracted using a standard protocol. Reverse transcription was conducted using the reverse transcriptase M-MLV from Promega. Real-time PCR was performed in triplicate using the SYBR Green PCR Mix (Promega) on an ABI 7500 sequence detection system (Applied Biosystems). Quantitative PCR reactions were performed under conditions that were standardized for each primer.