1. Introduction
Rapid improvements in next-generation sequencing techniques have contributed to the routine application of whole genome and tissue-specific sequencing for the identification of recurrent somatic mutations, many of which are driver mutations of tumorigenesis (Chi et al. Reference Chi, Allis and Wang2010; Dawson & Kouzarides, Reference Dawson and Kouzarides2012; Shih et al. Reference Shih, Abdel-Wahab, Patel and Levine2012), thereby providing genomic cues of human disease. Complementary chromatin immunoprecipitation-based sequencing techniques have provided a comprehensive map of the epigenome, in the process revealing that many cancers are accompanied by changes in the epigenome (Dawson & Kouzarides, Reference Dawson and Kouzarides2012). Further, accumulating evidence has established that disruption of the epigenome is a fundamental mechanism in cancer, which has prompted the development of therapeutic approaches to target the epigenetic machinery. DNA methyltransferases (DNMTs) and histone deacetylaseses (HDACs) were among the first epigenetic regulatory enzymes to be targeted for drug design (Minucci & Pelicci, Reference Minucci and Pelicci2006; Szyf, Reference Szyf2009), with several of these inhibitors approved by the Food and Drug Administration (FDA) for clinical use in hematological malignancies (Arrowsmith et al. Reference Arrowsmith, Bountra, Fish, Lee and Schapira2012). These advances in turn led to a concerted effort toward identification and characterization of novel inhibitors for other histone-modifying enzymes. Histone methylation has long been considered one of the most prevalent chromatin modifications, one that is written by histone lysine methyltransferases (KMTs) and erased by histone lysine demethylases (KDMs). The establishment and maintenance of an appropriate histone methylation landscape is not only critical for normal cell development, but when altered, closely related to human disease. In addition, an increasing amount of biological data has indicated that aberrant enzymatic activities of a majority of these histone KMTs and KDMs are closely related to human cancer. Within the past decade, tremendous progress has been made toward identification initially of inhibitors that target histone KMTs (Copeland et al. Reference Copeland, Solomon and Richon2009), and soon thereafter, their corresponding KDMs, with potent clinical potential. In this review, we focus on continuing progress toward the design of inhibitors that target human KMTs and KDMs. We start with a short introduction on the targeted enzymes, followed by an emphasis on the contributions of structural biology toward the development of potent inhibitors against specific KMTs and KDMs. We anticipate that a compendium of available successful examples of epigenetic inhibitor discovery against histone KMTs and KDMs will propel the field forward toward the identification and application of new approaches for drug discovery targeted to these epigenetic enzymes in the coming decade.
2. Small molecule inhibitors of histone lysine methyltransferases
Except for DOT1L, all identified histone KMTs are composed of a conserved approx. 130 amino acid (aa)-long SET domain, which adopts a unique topology composed of a series of β-strands folded into three β-sheets that surround a knot-like structure (Cheng et al. Reference Cheng, Collins and Zhang2005). In addition, functional SET domain folds are usually flanked by two closely packed cysteine-rich modules called pre-SET and post-SET domains that are keys to enzymatic activity. Functionally, all histone KMTs catalyze methyl transfer from donor S-adenosyl methionine (SAM) (Fig. 1a) to substrate lysine, resulting in formation of methylated lysine product and S-adenosyl homocysteine (SAH) (Fig. 1 b).
Illustration of the methyl-transfer reaction and inhibitors of H3K79 methyltransferase. (a, b) Chemical formula-based illustration of the methyl-transfer reaction. The active methyl group is colored in red. Chemical formulas of SAM and SAH are shown in panels (a) and (b), respectively. (c) Crystal structure of DOT1L with bound SAM, shown in a stick representation (PDB: 1NW3). The activation and substrate-binding loops are colored in purple. The N-terminal domain and the open α/β structure are colored in cyan and blue, respectively. (d) Chemical formula of EPZ004777. (e) Structure of complex of DOT1L with bound EPZ004777 (PDB: 4ER3). Color representations are the same as in panel (c).
2.1 Inhibitors of H3K79 methyltransferases
Methylation of H3K79 was first identified as a bona fide histone modification in the nucleosome core by mass spectroscopy, for which DOT1L was identified as the catalytic enzyme (Feng et al. Reference Feng, Wang, Ng, Erdjument-Bromage, Tempst, Struhl and Zhang2002). Although not a translocation partner of MLL1, DOT1L is recruited by the super elongation complex formed by MLL1 fusion partners (Chi et al. Reference Chi, Allis and Wang2010), with aberrant enzymatic activity of DOT1L shown to be a driving force of leukemogenesis (Bernt et al. Reference Bernt, Zhu, Sinha, Vempati, Faber, Krivtsov, Feng, Punt, Daigle, Bullinger, Pollock, Richon, Kung and Armstrong2011). Thus, DOT1L has long been considered to be an excellent therapeutic target for MLL1-related leukemia. As the only non-SET domain containing KMT identified in humans, the catalytic core of DOT1L adopts an elongated structure composed of an N-terminal domain and an open α/β-fold for the C-terminal domain (Fig. 1 c), a folding topology more similar to histone arginine methyltransferases (RMTs) (Min et al. Reference Min, Feng, Li, Zhang and Xu2003; Sawada et al. Reference Sawada, Yang, Horton, Collins, Zhang and Cheng2004). The bound SAM cofactor exhibited excellent shape complementary within the pocket formed between an activation loop and the open α/β-fold in the structure of human DOT1L, whereas the potential histone substrate is expected to occupy the cleft between the activation loop and the substrate-binding loop (Fig. 1 c). Guided by the above structural information, epizyme synthesized a SAM mimic named EPZ004777 (Fig. 1 d), which inhibits DOT1L enzymatic activity with extremely high potency (IC50=0.4 nM) (Daigle et al. Reference Daigle, Olhava, Therkelsen, Majer, Sneeringer, Song, Johnston, Scott, Smith, Xiao, Jin, Kuntz, Chesworth, Moyer, Bernt, Tseng, Kung, Armstrong, Copeland, Richon and Pollock2011). The chemical structure of EPZ004777 retains the nucleoside core of the SAM substrate, while the aa moiety was coupled with a bulky para-tert-butylphenyl group through a urea linkage, which should be unable to fit into the SAM-binding pocket. The structure of DOT1L bound to EPZ004777 validated the concept that the deazaadenosine moiety exhibits similar interactions as SAM with DOT1L, while the tert-butylphenyl end induced a large structural rearrangement within the substrate-binding and activation loops (Yu et al. Reference Yu, Chory, Wernimont, Tempel, Scopton, Federation, Marineau, Qi, Barsyte-Lovejoy, Yi, Marcellus, Iacob, Engen, Griffin, Aman, Wienholds, Li, Pineda, Estiu, Shatseva, Hajian, Al-Awar, Dick, Vedadi, Brown, Arrowsmith, Bradner and Schapira2012). Notably, Tyr312, a key residue within the substrate-binding loop known to be important for KMT activity, was rotated away from this cleft (Fig. 1 e), thereby indicating that EPZ004777 may also affect substrate binding. Despite the broad utilization of SAM by all histone KMTs, the SAM mimic EPZ004777 displayed remarkable inhibition selectivity toward DOT1L over all other methyltransferases tested in this study. EPZ004777 could selectively inhibit H3K79 methylation and blocked expression of leukemogenic genes in MLL cells, with little effect on non-MLL-translocated cells, indicative of promising potential as a drug candidate. However, poor pharmacokinetic properties of EPZ004777 limited its utility in an in vivo setting (Daigle et al. Reference Daigle, Olhava, Therkelsen, Majer, Sneeringer, Song, Johnston, Scott, Smith, Xiao, Jin, Kuntz, Chesworth, Moyer, Bernt, Tseng, Kung, Armstrong, Copeland, Richon and Pollock2011). Taking advantage of the hydrophobic environment surrounding the adenine ring of SAM in DOT1L, further optimization on EPZ004777 led to the identification of a more potent compound SGC0946, which added a bromine atom at position 7 of the adenine ring (Yu et al. Reference Yu, Chory, Wernimont, Tempel, Scopton, Federation, Marineau, Qi, Barsyte-Lovejoy, Yi, Marcellus, Iacob, Engen, Griffin, Aman, Wienholds, Li, Pineda, Estiu, Shatseva, Hajian, Al-Awar, Dick, Vedadi, Brown, Arrowsmith, Bradner and Schapira2012). SGC0946 showed improved potency in in vitro studies (IC50=0.3 nm), as well as in cellular assays. Of interest, just addition of a single bromine atom at the N7 position of SAH (BrSAH) resulted in an 8-fold increase in potency against DOT1L (IC50=77 nM) (Yu et al. Reference Yu, Smil, Li, Tempel, Fedorov, Nguyen, Bolshan, Al-Awar, Knapp, Arrowsmith, Vedadi, Brown and Schapira2013), which verified that the N7 position on the adenine ring is available for further optimization. In a parallel study, the Song laboratory showed that addition of a methyl or benzyl group to the 6-NH2 group of SAM resulted in more selective inhibitors toward DOT1L, consistent with DOT1L having a more spacious hydrophobic binding pocket able to accommodate extra groups (Yao et al. Reference Yao, Chen, Diao, Cheng, Deng, Anglin, Prasad and Song2011). Both examples provide valuable insights into the potential impact of beneficial structure-based modifications positioned either on the adenine ring or the aa moiety toward SAM-based inhibitor design.
2.2 Inhibitors of H3K4 and H3K27 methyltransferases
Histone H3K4 methylation is a hallmark of transcriptional activation, which in humans is established by the Trithorax group of proteins. H3K4 KMTs include SET1A and SET1B, as well as the MLL family proteins MLL1 to MLL4, which exhibit full activity only when complexed with the proteins WDR5, RbBP5, Ash2L and other components (Ruthenburg et al. Reference Ruthenburg, Allis and Wysocka2007). Although aberrant expression of these histone KMTs are closely related to human cancer, as shown for MLL1 (Krivtsov & Armstrong, Reference Krivtsov and Armstrong2007) and MLL2 (Dawson & Kouzarides, Reference Dawson and Kouzarides2012), the loss of enzymatic activity is not the causal reason for cancer. In addition, some of these H3K4 methyltransferases exhibit redundant function in cellular processes. In this regard, inhibitors that target the catalytic activity of H3K4 methyltransferases may not be promising candidates for cancer therapy, partially accounting for the slow progress in inhibitor design toward these enzymes.
In strong contrast to H3K4 methylation, H3K27 methylation is a hallmark of gene repression established by the Polycomb group of proteins, which function in opposition to the human Trithorax group of proteins. In practice, both Trithorax and Polycomb group of proteins function together in maintaining ‘cellular memory’. The core enzymes that add H3K27 methylation marks are EZH2 and EZH1, which exhibit full activity only when complexed with the proteins EED, Suz12, RbBP4 and other components required to form the PRC2 complex (Margueron & Reinberg, Reference Margueron and Reinberg2011). EZH2 overexpression is implicated in several human cancers including both leukemia and solid tumors (Chi et al. Reference Chi, Allis and Wang2010; Shih et al. Reference Shih, Abdel-Wahab, Patel and Levine2012). Especially, somatic mutations of a specific tyrosine (Tyr641), which result in altered catalytic efficiency, were reported to be associated with follicular lymphoma and diffuse large B-cell lymphoma (Shih et al. Reference Shih, Abdel-Wahab, Patel and Levine2012), indicating that pharmacological inhibition of EZH2 activity may provide a promising treatment for EZH2-associated cancer. An initial effort aimed toward identification of HDAC inhibitors found that SAH hydrolase inhibitor 3-deazaneplanocin A (DZNep, Fig. 2 a) could also function as an EZH2 inhibitor through depletion of the cellular levels of PRC2 components (Tan et al. Reference Tan, Yang, Zhuang, Jiang, Chen, Lee, Karuturi, Tan, Liu and Yu2007). Inhibition of SAH hydrolase resulted in the cellular accumulation of SAH, which in turn caused by-product inhibition of SAM-dependent histone KMTs (Chiang & Cantoni, Reference Chiang and Cantoni1979). However, as a SAH hydrolase inhibitor, DZNep could also non-specifically inhibit other KMTs (Miranda et al. Reference Miranda, Cortez, Yoo, Liang, Abe, Kelly, Marquez and Jones2009), indicating that DZNep was a selective but non-specific EZH2 inhibitor.
Inhibitors of H3K27 and H3K9 methyltransferases. (a, b, c, d) Chemical formulas of inhibitors DZNep (panel a), EPZ005687 (panel b), GSK126 (panel c) and EI1 (panel d) that target EZH2. (e, f, g) Chemical formulas of Chaetocin (panel e), PS-ETP-1 (panel f) and BX-01294 (panel g). (h) Structure of the complex of GLP with bound SAH and inhibitor BIX-01294 (PDB: 3FPD). SAH and BIX-01294 are presented in stick representations. Pre-SET and SET domains are colored in blue, while Post-SET is colored in purple. (i) Chemical formula of UNC0638. (j) Structure of the complex of G9a with bound SAH and inhibitor UNC0638 (PDB: 3RJW). Color representations are the same as in panel (h). (k) Chemical formula of BRD4770.
Recently, three groups reported on the development of a new generation of potent EZH2 specific inhibitors through a combination of high throughput screening and chemistry-based optimization. All three inhibitors exhibited high potency against EZH2-containing PRC2 (EPZ005687 has an IC50 of 54 nM, GSK126 has an IC50 of 9.9 nM and EI1 showed an IC50 of 15 nM; chemical formulas in Figs 2 b–2 d), while not inhibiting other tested KMTs, indicative of enhanced selectivity (Knutson et al. Reference Knutson, Wigle, Warholic, Sneeringer, Allain, Klaus, Sacks, Raimondi, Majer, Song, Scott, Jin, Smith, Olhava, Chesworth, Moyer, Richon, Copeland, Keilhack, Pollock and Kuntz2012; McCabe et al. Reference McCabe, Ott, Ganji, Korenchuk, Thompson, Van Aller, Liu, Graves, Della Pietra, Diaz, LaFrance, Mellinger, Duquenne, Tian, Kruger, McHugh, Brandt, Miller, Dhanak, Verma, Tummino and Creasy2012; Qi et al. Reference Qi, Chan, Teng, Li, Chuai, Zhang, Zeng, Li, Fan, Lin, Gu, Ardayfio, Zhang, Yan, Fang, Mi, Zhang, Zhou, Feng, Chen, Li, Yang, Zhao, Liu, Yu, Lu, Atadja and Li2012). Additionally, the inhibition was 50- to 150-fold less potent for closely related EZH1-containing PRC2. The chemical composition of these three compounds are very similar, as reflected in the central indazole group in EPZ005687 and the indole rings in GSK126 and EI1, the invariant pyridone-containing group attached to the 4 position of the indole/indazole rings, as well as the bulky lipophilic groups at the N1 position of the indole/indazole scaffolds. Although structural information is not yet available for the catalytic domain of EZH2, homologous modeling studies indicate that the indole rings could occupy the sugar ring position of SAM, while the pyridone-containing group could occupy the adenine ring position of SAM (McCabe et al. Reference McCabe, Ott, Ganji, Korenchuk, Thompson, Van Aller, Liu, Graves, Della Pietra, Diaz, LaFrance, Mellinger, Duquenne, Tian, Kruger, McHugh, Brandt, Miller, Dhanak, Verma, Tummino and Creasy2012). In addition, it has been proposed that four of the six residue differences between EZH2 and EZH1 within the post-SET domain could account for the observed selectivity, although this claim remains to be validated following eventual structure determination of these key enzymes. All three compounds could decrease global H3K27me3 levels and also effectively inhibit the proliferation of cell lines that harbor EZH2 mutations, indicative of the potential clinical usage of EZH2 inhibitors in cancer therapy. Clearly, the structure of EZH2 and its interacting components would provide a critical platform for future inhibitor design, followed by clinical testing and implementation of epigenetic drugs in cancer treatment.
2.3 Inhibitors of H3K9 methyltransferases
Similar to H3K27 methylation, histone H3K9 methylation is also closely related to repressive chromatin states. The K9me3 and K9me2 marks on H3 serve as docking sites for HP1 (heterochromatin protein 1), a process critical for heterochromatin formation (Grewal & Jia, Reference Grewal and Jia2007). Methylation of H3K9 in humans is dependent on the KMTs SUV39H1, SUV39H2, G9a, GLP, SETDB1, SETDB2 and RIZ1 (Krauss, Reference Krauss2008). Most of these H3K9 methyltransferases are reported to be related to various kinds of diseases, with G9a showing increased expression in lung cancer cell lines (Chen et al. Reference Chen, Hua, Kao, Chi, Wei, Johansson, Shiah, Chen, Jeng, Cheng, Lai, Chang, Jan, Chien, Yang, Huang, Hsiao and Kuo2010), and SUV39H1 associated with colon cancer (Kang et al. Reference Kang, Lee, Kim, Chang, Kyu Park, Chun, Song, Park and Kim2007). Although it has been shown that the KMT product SAH and the natural product analog Sinefungin could serve as inhibitors through SAM-competitive inhibition (Couture et al. Reference Couture, Hauk, Thompson, Blackburn and Trievel2006), these natural products lack selectivity. Chaetocin was also shown to efficiently inhibit SUV39H1 (IC50=0.8 μM, Fig. 2 e) in a SAM-competitive manner (Greiner et al. Reference Greiner, Bonaldi, Eskeland, Roemer and Imhof2005). However, as it also inhibited G9a, several other methyltransferases, and some functionally unrelated proteins, further research indicated that Chaetocin may inhibit histone KMTs in a non-specific manner through chemical modification of enzymes mediated by the disulfide groups of this inhibitor (Cherblanc et al. Reference Cherblanc, Chapman, Brown and Fuchter2013). A simplified Chaetocin derivative PS-ETP-1 (Fig. 2 f) was shown to exhibit a comparable potency against G9a (IC50=5.2 μM) with significantly reduced cytotoxicity (Fujishiro et al. Reference Fujishiro, Dodo, Iwasa, Teng, Sohtome, Hamashima, Ito, Yoshida and Sodeoka2013), which may help to resolve this controversy. In addition, structural information will be required as part of an effort to resolve this ambiguity.
An extensive screening targeted to G9a led to the discovery of the first potent and selective histone KMT inhibitor, BIX-01294 (Kubicek et al. Reference Kubicek, O'Sullivan, August, Hickey, Zhang, Teodoro, Rea, Mechtler, Kowalski, Homon, Kelly and Jenuwein2007). The drug BIX-01294 is a diazepin–quinazolin–amine derivative (Fig. 2 g), which specifically inhibited G9a (IC50=1.7 μM) and the closely related GLP enzyme. The structure of BIX-01294 bound to the SET domain of GLP established that this inhibitor occupied a region of the histone substrate-binding groove that overlapped with the region bound by the K4 to R8 segment of histone H3 (Fig. 2 h) (Chang et al. Reference Chang, Zhang, Horton, Upadhyay, Spannhoff, Liu, Snyder, Bedford and Cheng2009), indicative of a substrate-competitive inhibition mechanism. However, BIX-01294 exhibited cytotoxicity at concentrations higher than 4.1 μM, which limited its usage in cell-based assays. Subsequent structure-guided chemical optimization resulted in the development of two BIX-01294 derivatives, named E72 (Chang et al. Reference Chang, Ganesh, Horton, Spannhoff, Liu, Sun, Zhang, Bedford, Shinkai, Snyder and Cheng2010) and UNC0321 (Liu et al. Reference Liu, Chen, Allali-Hassani, Quinn, Wigle, Wasney, Dong, Senisterra, Chau, Siarheyeva, Norris, Kireev, Jadhav, Herold, Janzen, Arrowsmith, Frye, Brown, Simeonov, Vedadi and Jin2010), which exhibited improved inhibitory potency through addition of a 7-alkoxyamine tethered to the quinazoline core so as to mimic the K9 side chain of H3 substrate. Especially, UNC0321 had a Morrision Ki of 63 pM, which was 250-fold more potent than BIX-01294. However, UNC0321 was less potent in cellular assays, possibly due to poor cell membrane permeability. Further optimization aimed at increasing lipophilicity led to the discovery of UNC0638 (Fig. 2 i), which was also a substrate-competitive inhibitor, but displayed superior selectivity, as well as improved cellular potency, together with low toxicity (Vedadi et al. Reference Vedadi, Barsyte-Lovejoy, Liu, Rival-Gervier, Allali-Hassani, Labrie, Wigle, Dimaggio, Wasney, Siarheyeva, Dong, Tempel, Wang, Chen, Chau, Mangano, Huang, Simpson, Pattenden, Norris, Kireev, Tripathy, Edwards, Roth, Janzen, Garcia, Petronis, Ellis, Brown, Frye, Arrowsmith and Jin2011). The crystal structure of the G9a-UNC0638-SAH complex established that the core segment of UNC0638 occupied the histone peptide-binding channel, while the 7-(3-pyrrolidin-1-yl)-propoxy side chain was inserted into the lysine-binding pocket (Fig. 2 j), thereby explaining the improved potency over BIX-01294 (Vedadi et al. Reference Vedadi, Barsyte-Lovejoy, Liu, Rival-Gervier, Allali-Hassani, Labrie, Wigle, Dimaggio, Wasney, Siarheyeva, Dong, Tempel, Wang, Chen, Chau, Mangano, Huang, Simpson, Pattenden, Norris, Kireev, Tripathy, Edwards, Roth, Janzen, Garcia, Petronis, Ellis, Brown, Frye, Arrowsmith and Jin2011). Recently, a SAM mimic inhibitor BRD4770 (Fig. 2 k) was synthesized and tested for efficacy against G9a (Yuan et al. Reference Yuan, Wang, Paulk, Kubicek, Kemp, Adams, Shamji, Wagner and Schreiber2012). Although BRD4770 exhibited lower binding affinity (EC50=∼5 μM), it displayed no apparent cellular toxicity, indicative of its potential usage to decipher the role of G9a in cancer cell biology.
2.4 Inhibitors of H3K36 methyltransferases
Histone H3K36 methylation is a well-established histone mark primarily associated with transcriptional activation, although it has also been implicated in some instances in transcriptional repression. H3K36 methylation has been identified in processes such as alternative splicing, dosage compensation, as well as DNA repair and recombination (Wagner & Carpenter, Reference Wagner and Carpenter2012). The well-established enzymes that catalyze H3K36 methylation include the NSD family of methyltransferases (NSD1, NSD2 and NSD3), SETD2, ASH1L and SMYD2 (Wagner & Carpenter, Reference Wagner and Carpenter2012). Many of these enzymes are also closely correlated with human cancer, with the SMYD2 gene frequently amplified in esophageal squamous cell carcinoma and other solid tumors (Komatsu et al. Reference Komatsu, Imoto, Tsuda, Kozaki, Muramatsu, Shimada, Aiko, Yoshizumi, Ichikawa, Otsuji and Inazawa2009; Rooney et al. Reference Rooney, Murray, Stevenson, Haites, Cassidy and McLeod1999), whereas defects in SETD2 are correlated with sporadic clear renal cell carcinoma (Wagner & Carpenter, Reference Wagner and Carpenter2012).
Compared with other members of this family of methyltransferases, SMYD2 may not be a bona fide H3K36 methyltransferase, given that it also targets non-histone substrates, such as p53 and RB (retinoblastoma) proteins. Through an application of AlphaScreen technology (Eglen et al. Reference Eglen, Reisine, Roby, Rouleau, Illy, Bosse and Bielefeld2008), Ferguson and colleagues identified a small molecule compound AZ505 (Fig. 3 a) as a SMYD2-specific inhibitor (Ferguson et al. Reference Ferguson, Larsen, Howard, Pollard, Green, Grande, Cheung, Garcia-Arenas, Cowen, Wu, Godin, Chen and Keen2011). AZ505 has an IC50 of 0.12 μM against SMYD2, while showing much weaker activity against other histone KMTs. Structures of complexes of SMYD2 with bound p53 (366–374) peptide (Fig. 3 b) or AZ505 (Fig. 3 c) revealed that AZ505 occupied the p53-binding pocket within the concaved cleft of SMYD2, which is formed by the SET domain together with imbedded MYND domain, post-SET and C-terminal lobe of the enzyme (Ferguson et al. Reference Ferguson, Larsen, Howard, Pollard, Green, Grande, Cheung, Garcia-Arenas, Cowen, Wu, Godin, Chen and Keen2011), suggestive of AZ505 being a peptide-competitive inhibitor. The benzooxazinone group of AZ505 was positioned within the lysine-binding channel of the substrate, where it formed extensive electrostatic and hydrophobic interactions with the donor methyl group of SAM and residues from SMYD2, consistent with the observation that AZ505 binding was dependent on bound SAM. The cyclohexyl group of AZ505 was positioned in a hydrophobic pocket normally occupied by Leu369 of p53 peptide, whereas the dichlorophenethyl moiety extended across the peptide-binding groove of SMYD2 to reach another hydrophobic pocket adjacent to the cofactor-binding site, thereby most likely accounting for the observed selectivity. AZ505 binding was mainly mediated by hydrophobic interactions, as shown by limited formation of hydrogen bonds. Given that the inhibitory effect of AZ505 required the presence of SAM, its potency in cellular-based assays needs to be tested.
Inhibitors of H3K36 and H4K20 methyltransferase. (a) Chemical formula of AZ505. (b) Structure of complex of SMYD2 with bound SAM and peptide substrate p53 (366–378) (PDB: 3S7F). Pre-SET and SET domain are colored in blue. The imbedded MYDN domain is colored in green, the Post-SET is colored in purple, whereas the C-terminal domain is colored in cyan. Bound SAM and peptide are shown in stick representations. (c) A blow-up view of the structure of the complex of SMYD2 with bound SAM and AZ505 (PDB: 3S7B). Color schemes are the same as in panel (b). (d, e) Chemical formulas of Sinefungin (panel d) and Pr-SNF (panel e). (f) Structure of complex of SETD2 with bound SAH (PDB: 4H12). N-terminal region and the SET domain are colored in blue; the imbedded AWS domain is colored in green, while the Post-SET is colored in purple. Bound zinc ions are shown as white spheres. SAH is shown a stick representation. (g) A blow-up view of the structure of the complex of SETD2 with bound Pr-SNF (PDB: 4FMU), with the overlapped post-SET loop shown in panel f, colored in silver. Note the conformational change of the post-SET loop in the two structures, and the movement of the Arg within the loop.
In an attempt to identify SAM-based inhibitors through modification of the known SAM analog Sinefungin (Fig. 3 d), the Luo laboratory synthesized several N-alkyl sinefungin derivatives and identified N-propyl-modified sinefungin (Pr-SNF, Fig. 3 e), which displayed 2- to 200-fold preference for SETD2 relative to other KMTs (Zheng et al. Reference Zheng, Ibanez, Wu, Blum, Zeng, Dong, Li, Hajian, Allali-Hassani, Amaya, Siarheyeva, Yu, Brown, Schapira, Vedadi, Min and Luo2012). In addition, Pr-SNF also showed more than 10-fold potency against SETD2 compared with unmodified Sinefungin. A structure-activity relationship study revealed that Pr-SNF occupied the same binding pocket as SAH (Fig. 3 f), with the propyl group of Pr-SNF partially extending into the lysine-binding pocket to induce a conformational change of the post-SET loop (Fig. 3 g), which otherwise adopted an auto-inhibitory conformation with an Arg from the post-SET loop forming a hydrogen bond with the cofactor SAH (Fig. 3 f). In this substrate accessible conformation, the secondary amine of Pr-SNF is ideally positioned to form two hydrogen bonds with the backbone carbonyl oxygens from SETD2 (Fig. 3 g), thereby potentially accounting for the observed selectivity. This approach has identified a promising route for inhibitor design through optimization of previously identified pan-inhibitors.
2.5 Inhibitors for H4K20 methyltransferases
Histone H4K20 methylation is the only well-known methylation mark on H4, with SETD8, Suv4-20h1 and Suv4-20h2 representing major KMTs responsible for methylation of this site (Beck et al. Reference Beck, Oda, Shen and Reinberg2012). H4K20 methylation is related to biological processes ranging from DNA damage response, to mitotic condensation, DNA replication and gene regulation. Generally speaking, H4K20 methylation is considered as a mark for transcriptional repression. Recently, the first selective SAM-competitive inhibitor of SETD8, Nahuoic acid A (Fig. 3 h), was identified through a screen using a library of marine organism extracts and pure marine natural products (Williams et al. Reference Williams, Dalisay, Li, Amphlett, Maneerat, Chavez, Wang, Matainaho, Yu, Brown, Arrowsmith, Vedadi and Andersen2013). Nahuoic acid A specifically inhibited SETD8 (IC50=6.5 μM) but not other tested KMTs. Although cellular based analysis has not yet been reported for Nahuoic acid A, it has been anticipated that this first natural product inhibitor of an H4K20 methyltransfrase would encourage follow-up studies directed toward this family of KMTs.
3. Small molecule inhibitors of histone lysine demethylases
Histone lysine methylation used to be considered as a permanent modification, until the seminal discovery of the first lysine-specific demethylase 1 (LSD1) (Shi et al. Reference Shi, Lan, Matson, Mulligan, Whetstine, Cole and Casero2004). LSD1, and its closely related family member LSD2, are the only KDMs discovered to date that remove methyl groups from methylated-lysines through an flavin adenine dinucleotide (FAD dependent oxidation mechanism. LSD1 and LSD2 cannot remove the methyl group from trimethylated-lysines, since the demethylation mechanism requires a protonated amine (shown schematically in Fig. 4 a). In a subsequent groudbreaking discovery, the large Jumonji C-domain-containing family of enzymes were also shown to be histone KDMs (Tsukada et al. Reference Tsukada, Fang, Erdjument-Bromage, Warren, Borchers, Tempst and Zhang2006), which catalyze the reaction through a hydroxylation pathway using Fe2+ and α-ketoglutarate (α-KG) as cofactors (shown schematically in Fig. 4 b). Since they do not require the existence of a lone pair of electrons on the nitrogen atom of the methylated-lysines, the Jumonji family of KDMs can also demethylate trimethylated-lysines.
Illustration of two identified mechanisms for lysine demethylation. (a) LSD1-mediated lysine demethylation mechanism. (b) Jumonji demethylase family-mediated lysine demethylation mechanism.

3.1 Inhibitors of H3K4 demethylases
Histone H3K4 demethylases include LSD1, LSD2 and the JARID1 family of proteins, with most of these enzymes impacting on cancer. LSD1 is aberrantly expressed in breast and prostate cancers (Chi et al. Reference Chi, Allis and Wang2010), making it an attractive drug target candidate. As the first known histone KDM, extensive efforts have been made using different approaches toward LSD1 inhibitor design (Suzuki & Miyata, Reference Suzuki and Miyata2011). The first generation of LSD1 inhibitors that were discovered was based on its homology with monoamine oxidases (MAOases), which also use FAD as the cofactor. Known MAOase inhibitors such as phenelzine, pargyline and trans-2-phenylcyclopropylamine (PCPA, IC50=243 μM) (chemical formulas in Figs 5 a–5 c) were tested and found to possess inhibitory activity toward LSD1 (Lee et al. Reference Lee, Wynder, Schmidt, McCafferty and Shiekhattar2006; Metzger et al. Reference Metzger, Wissmann, Yin, Muller, Schneider, Peters, Gunther, Buettner and Schule2005; Schmidt & McCafferty, Reference Schmidt and McCafferty2007). These inhibitors possibly impact on the demethylation activity by forming a covalent adduct with the FAD, as shown in the structure of LSD1 with bound PCPA (Yang et al. Reference Yang, Culhane, Szewczuk, Jalili, Ball, Machius, Cole and Yu2007b), wherein the cyclopropyl ring of PCPA adducted with FAD to form a five-membered ring through covalent bonding with the N(5) and C(4a) atoms (Fig. 5 d). In turn, structure-guided generation of PCPA derivatives showed improved potency (Benelkebir et al. Reference Benelkebir, Hodgkinson, Duriez, Hayden, Bulleid, Crabb, Packham and Ganesan2011; Binda et al. Reference Binda, Valente, Romanenghi, Pilotto, Cirilli, Karytinos, Ciossani, Botrugno, Forneris, Tardugno, Edmondson, Minucci, Mattevi and Mai2010; Gooden et al. Reference Gooden, Schmidt, Pollock, Kabadi and McCafferty2008; Mimasu et al. Reference Mimasu, Umezawa, Sato, Higuchi, Umehara and Yokoyama2010), but most of them have relatively poor selectivity, as they also inhibit MAOases. Recently, a novel PCPA derivative named OG-L002 (Fig. 5 e) has been identified to be a highly specific LSD1 inhibitor (IC50≈0.02 μM), showing a >36-fold selectivity over MAOases (Liang et al. Reference Liang, Quenelle, Vogel, Mascaro, Ortega and Kristie2013). In search of more selective inhibitors, a series of peptide-based inhibitors, such as compounds with the propargylamine or hydrazine functionality were synthesized in the context of a histone H3(1–21) peptide (Culhane et al. Reference Culhane, Szewczuk, Liu, Da, Marmorstein and Cole2006, Reference Culhane, Wang, Yen and Cole2010; Dancy et al. Reference Dancy, Ming, Papazyan, Jelinek, Majumdar, Sun, Dancy, Drury, Cotter, Taverna and Cole2012; Szewczuk et al. Reference Szewczuk, Culhane, Yang, Majumdar, Yu and Cole2007; Yang et al. Reference Yang, Culhane, Szewczuk, Gocke, Brautigam, Tomchick, Machius, Cole and Yu2007a). These peptide analogs are more selective and potent than the MAOase inhibitors, with hydrazine-H3 (Fig. 5 f) exhibiting the most potent LSD1 inhibitory activity (Ki=4.35 nM). A structural characterization revealed that the N-methylpropargyl-K4-modified H3 peptide formed a stable covalent complex with the cofactor FAD through the N(5) position (Fig. 5 g), which provided a structural basis toward understanding the peptide-based H3K4 demethylation mechanism by LSD1 (Yang et al. Reference Yang, Culhane, Szewczuk, Gocke, Brautigam, Tomchick, Machius, Cole and Yu2007a). However, poor membrane permeability limited the usage of these peptide-based inhibitors in cellular assays. In search of a selective LSD1 inhibitor suitable for in vivo assays, a PCPA-lysine hybrid inhibitor NCL-1 was synthesized (Ueda et al. Reference Ueda, Suzuki, Mino, Tsumoto, Nakagawa, Hasegawa, Sasaki, Mizukami and Miyata2009). NCL-1 (IC50=2.5 μM) exhibited selectivity for LSD1 over MAOases of 400 to 11 000 times higher than that of PCPA, with the optically active (1S, 2R)-NCL-1 (Fig. 5 h) exhibiting 4-fold more potency than its (1R, 2S)-NCL-1 stereoisomeric counterpart (Ogasawara et al. Reference Ogasawara, Suzuki, Mino, Ueda, Khan, Matsubara, Koseki, Hasegawa, Sasaki, Nakagawa, Mizukami and Miyata2011). Several N-alkylated PCPA analogs, such as compound A (Fig. 5 i), which exhibit high selectivity and potency (Ki=5 nM) have been patented (Suzuki & Miyata, Reference Suzuki and Miyata2011), but a detailed description of their mode of action remains to be been disclosed.
Inhibitors of LSD1. (a–c) Chemical formulas of inhibitors that target LSD1. Pheneizine is shown in panel (a), Pargyline is shown in panel (b) and 2-PCPA is shown in panel (c). (d) Structural illustration of the FAD–PCPA adduct (PDB: 2UXX). (e, f) Chemical formulas of inhibitors that target LSD1. OG-L002 is shown in panel (e) and H3(1–21)K4(hydrazine) is shown in panel (f). (g) Structural illustration of the FAD-H3(1–21)K4(N-methylpropargyl) adduct (PDB: 2UXN). (h–m) Chemical formulas of inhibitors that target LSD1. (1S, 2R)-NCL-1 is shown in panel (h), Compound A is shown in panel (i), Polyamine 2d is shown in panel (j), PG-11144 is shown in panel (k), CBB1007 is shown in panel (l) and Amidoxine-22 is shown in panel (m).
Besides the above-mentioned family of MAOases, LSD1 was found to be homologous with two polyamine oxidases, SMO (spermine oxidase) and APAO (N1-acetylpolyamine oxidase). Subsequently, polyamine-based inhibitors were designed and proved to be potent and selective LSD1 inhibitors (Huang et al. Reference Huang, Greene, Murray Stewart, Goodwin, Baylin, Woster and Casero2007, Reference Huang, Stewart, Wu, Baylin, Marton, Perkins, Jones, Woster and Casero2009; Sharma et al. Reference Sharma, Wu, Steinbergs, Crowley, Hanson, Casero and Woster2010). The most potent inhibitor, polyamine-2d (Fig. 5 j, IC50≈1 μM), exhibited noncompetitive inhibition kinetics, whereas another compound, PG-11144 (Fig. 5 k, IC50=5 μM), functions through a substrate-competitive inhibition mechanism, with both these compounds producing an increase in H3K4me2 levels. Importantly, both compounds when combined with DNMT inhibitor 5-Aza, exhibited marked increase in the inhibition of tumor growth in a human colon cancer tumor model (Huang et al. Reference Huang, Stewart, Wu, Baylin, Marton, Perkins, Jones, Woster and Casero2009). These compounds represent a new branch of epigenetic modulators targeted to LSD1 with the potential for usage as antitumor agents.
On the basis of the structural features of the active site of LSD1, several inhibitors with a non-peptide chemical scaffold were also designed, with the most potent compound CBB1007 (Fig. 5 l), displaying an IC50 of 5.27 μM (Wang et al. Reference Wang, Lu, Ren, Sun, Xu, Lan, Liu, Ward, Quan, Ye and Zhang2011). These compounds selectively inhibited LSD1 but not LSD2 and JARID1A in vitro, and also inhibited the proliferation of pluripotent cancer cells. These compounds may function in a similar manner as the polyamine-derivatives, as they are structurally related to (bis)guanidine and (bis)urea analogs (Huang et al. Reference Huang, Greene, Murray Stewart, Goodwin, Baylin, Woster and Casero2007; Sharma et al. Reference Sharma, Wu, Steinbergs, Crowley, Hanson, Casero and Woster2010).
A series of small molecule amidoximes were recently reported as LSD1 inhibitors discovered through a virtual screening strategy followed by chemical synthesis and optimization (Hazeldine et al. Reference Hazeldine, Pachaiyappan, Steinbergs, Nowotarski, Hanson, Casero and Woster2012). The best compound, Amidoxime-22 (Fig. 5 m), despite exhibiting moderate LSD1 inhibitory activity in vitro (IC50=16.8 μM), promoted dramatic increases in H3K4me2 levels in Calu-6 lung carcinoma cells. In addition, Amidoxime-22 and its derivatives are more drug-like (Lipinski et al. Reference Lipinski, Lombardo, Dominy and Feeney2001) and can be easily synthesized, representing a promising scaffold for further optimization.
3.2 Inhibitors of H3K9 demethylases
Histone H3K9 demethylases include LSD1, the JMJD2 subfamily of enzymes and JHDM2A (Shi, Reference Shi2007). Since all members of the Jumonji family of demethylases use Fe2+ and α-KG (Fig. 6 a) as cofactors to catalyze the demethylation process, the α-KG analog N-oxalylglycine (NOG, Fig. 6 b) has generally been used as an inhibitor for in vitro studies (Cloos et al. Reference Cloos, Christensen, Agger, Maiolica, Rappsilber, Antal, Hansen and Helin2006). The structure of the catalytic domain of JMJD2A with bound cofactor α-KG (or NOG), in the presence of inactive Ni2+ cation (as a replacement for Fe2+) and substrate H3 (1–15)K9me3 peptide has provided definitive insights into the demethylation mechanism (Couture et al. Reference Couture, Collazo, Ortiz-Tello, Brunzelle and Trievel2007; Ng et al. Reference Ng, Kavanagh, McDonough, Butler, Pilka, Lienard, Bray, Savitsky, Gileadi, von Delft, Rose, Offer, Scheinost, Borowski, Sundstrom, Schofield and Oppermann2007). In the complex, the catalytic JmjC domain adopts a conserved β-barrel fold associated with the cupin superfamily, surrounded by the JmjN domain, the connecting mixed domain from the N-terminal segment and a stabilizing C-terminal segment (Fig. 6 c). The α-KG interacts with amino acids from the active site from one side and with the Ni2+ cation in a bidentate manner from another side, with the Ni2+ cation hexa-coordinated through additional interactions with the His–His–Glu triad and a water molecule (Fig. 6 d).
Inhibitors for H3K9 demethylases. (a, b) Chemical formulas of cofactor/inhibitors targeted to H3K9 demethylases. α-KG is shown in panel (a) and NOG is shown in panel (b). (c) Structure of complex of JMJD2A with bound α-KG, Ni2+ cation and its substrate peptide H3(1–15)K9me3 (PDB: 2Q8C). JmjN domain is colored in green, the mixed region is colored in salmon, JmjC domain is colored in blue and the C-terminal region is colored in cyan. Ni2+ cation is colored in green. (d) A blow-up view of the active site shown in panel (c), with the histone peptide shown in a stick representation and JMJD2A shown in a surface view. (e–i) Chemical formulas of cofactor/inhibitors targeted to H3K9 demethylases. α-HG is shown in panel (e), PDCA is shown in panel (f), 5-carboxy-8HQ is shown in panel (h), and Bipyridine A is shown in panle i. (g, j) Blow-up views of the catalytic sites of JMJD2A with bound inhibitor PDCA (panel g) (PDB: 2VD7) and Bipyridine A (panel j) (PDB: 2PDQ). Color schemes are the same as in panel (c). (k–p) Chemical formulas of cofactor/inhibitors targeted to H3K9 demethylases. Caffeic acid is shown in panel (k), 4-hydroxypyrazole-A is shown in panel (l), Hydroxamate analog A is shown in panel (m), Methylstat is shown in panel (n), Ebselen is shown in panel (o) and H3(7–14)-T11(DNOC) is shown in panel (p).
In this regard, α-KG analogs were extensively investigated as cofactor-competitive inhibitors. Besides NOG, another α-KG analog, the oncometabolite α-hydroxyglutarate (α-HG, Fig. 6 e) was found to inhibit JMJD2A, JMJD2C and JHDM1A with IC50 values ranging from 24 μM to 106 μM (Chowdhury et al. Reference Chowdhury, Yeoh, Tian, Hillringhaus, Bagg, Rose, Leung, Li, Woon, Yang, McDonough, King, Clifton, Klose, Claridge, Ratcliffe, Schofield and Kawamura2011). To enhance the selectivity, a range of NOG derivatives were synthesized (Hamada et al. Reference Hamada, Kim, Suzuki, Itoh, Tsumoto, Nakagawa, Janknecht and Miyata2009; Rose et al. Reference Rose, Ng, Mecinovic, Lienard, Bello, Sun, McDonough, Oppermann and Schofield2008, Reference Rose, Woon, Kingham, King, Mecinovic, Clifton, Ng, Talib-Hardy, Oppermann, McDonough and Schofield2010), most of which showed limited selectivity. A screen using known inhibitors from other Fe2+/α-KG dependent oxygenases identified 2,4-pyridinedicarboxylic acid (PDCA, Fig. 6 f), which showed potent inhibitory activity on JMJD2E (IC50=1.4 μM) (Rose et al. Reference Rose, Ng, Mecinovic, Lienard, Bello, Sun, McDonough, Oppermann and Schofield2008). The structure of JMJD2A with bound PDCA verified that PDCA functions in an α-KG competitive manner. PDCA binds the Ni2+ cation in a bidentate manner via its N-atom and 2-carboxylate, whereas the 4-carboxylate mimics α-KG binding by forming two hydrogen bonds with a lysine and a tyrosine in the active site (Fig. 6 g). Structural information showed that the C-3 position of PDCA was compatible for further modification, which led to the identification of C-3-substituted PDCA with enhanced selectivity for JMJD2E over PHD2 (Thalhammer et al. Reference Thalhammer, Mecinovic, Loenarz, Tumber, Rose, Heightman and Schofield2011). High-throughput screening also identified PDCA-related compound 5-carboxy-8-hydroxyquinoline (5-Carboxy-8HQ, Fig. 6 h), which has a greater potency for JMJD2E (IC50=200 nM) (King et al. Reference King, Li, Sakurai, Kawamura, Rose, Ng, Quinn, Rai, Mott, Beswick, Klose, Oppermann, Jadhav, Heightman, Maloney, Schofield and Simeonov2010). Another family of pyridine-based compounds, with a 4-carboxy-2,2′-bipyridyl scaffold (abbreviated here as Bipyridine A, Fig. 6 i), have also been reported to show enhanced inhibitory activity toward JMJD2 subfamily of demethylases (Chang et al. Reference Chang, King, Tumber, Woon, Heightman, McDonough, Schofield and Rose2011). Structural characterization validated that this bipyridine compound also functions through α-KG competition. The carboxylate group of the compound is positioned to interact with the same Tyr and Lys in a manner analogous to that observed for α-KG and PDCA, whereas the active site Ni2+ cation was chelated through both pyridinyl nitrogen atoms (Fig. 6 j). This bipyridine scaffold may represent an important element for active site metal chelation, given that it has also been shown to be a key scaffold in the JMJD3 complex with bound inhibitor GSK-J1 (to be discussed below) (Kruidenier et al. Reference Kruidenier, Chung, Cheng, Liddle, Che, Joberty, Bantscheff, Bountra, Bridges, Diallo, Eberhard, Hutchinson, Jones, Katso, Leveridge, Mander, Mosley, Ramirez-Molina, Rowland, Schofield, Sheppard, Smith, Swales, Tanner, Thomas, Tumber, Drewes, Oppermann, Patel, Lee and Wilson2012).
More recently, screening efforts using a library of pharmacologically active compounds or natural products led to the identification of hydroxamic acids and catechols as JMJD2 demethylase inhibitors (Nielsen et al. Reference Nielsen, Kristensen, Stephansen, Kristensen, Helgstrand, Lees, Cloos, Helin, Gajhede and Olsen2012; Sakurai et al. Reference Sakurai, Rose, Schultz, Quinn, Jadhav, Ng, Oppermann, Schofield and Simeonov2010). One such compound, caffeic acid (Fig. 6 k), could inhibit JMJD2C with an IC50 of 13.7 μM. However, these compounds are more promiscuous as they also inhibit HDACs and DNMTs. A 4-hydroxypyrazole scaffold-containing compound named 4-hydroxypyrazole-A (Fig. 6 l) was also reported as an inhibitor of JMJD2C, with the best hit exhibiting an IC50 of 147 μM. Based on the crystal structure of JMJD2A and a homology model of JMJD2C, Hamada and colleagues designed a series of hydroxamate analogs bearing a tertiary amine and identified one such analog, named Hydroxamate analogue A (Fig. 6 m), to be a potent and selective JMJD2 inhibitor, showing 500-fold greater JMJD2C inhibitory activity over NOG (Hamada et al. Reference Hamada, Suzuki, Mino, Koseki, Oehme, Flamme, Ozasa, Itoh, Ogasawara, Komaarashi, Kato, Tsumoto, Nakagawa, Hasegawa, Sasaki, Mizukami and Miyata2010). They also identified prodrugs of this analog, which showed synergistic growth inhibition of cancer cells in combination with a LSD1 inhibitor. With the idea of combining the mimicry of both lysine and α-KG based on available structural information of the active site, Luo et al. synthesized a novel JMJD-selective inhibitor, which could selectively inhibit the JMJD2 family of KDMs (Luo et al. Reference Luo, Liu, Kubicek, Myllyharju, Tumber, Ng, Che, Podoll, Heightman, Oppermann, Schreiber and Wang2011). In addition, its prodrug methylstat (Fig. 6 n) selectively inhibited the Jumonji family of demethylases in cells and could inhibit cell growth of the esophageal carcinoma cell line KYSE150 (GI50=5.1 μM). Both methylstat and the hydroxamate analog showed cancer cell inhibitory activity, indicative of clinical promise in anticancer treatment.
In an additional interesting finding, Sekirnik and colleagues showed that Zn-ejecting compounds including disulfiram and ebselen (Fig. 6 o, IC50=10.6 μM) could also inhibit JMJD2A by removing the zinc ion from the enzyme (Sekirnik et al. Reference Sekirnik, Rose, Thalhammer, Seden, Mecinovic and Schofield2009). As zinc finger motifs are a ubiquitous scaffold among many enzymes including several members of the Jumonji KDM family, Zn ejecting compounds should be used with caution for inhibitor design.
Most of above ligand-based inhibitors showed limited selectivity between subfamilies/isoforms of these Jumonji demethylases. To enhance selectivity, a peptidomimetic approach facilitated by structural information of the catalytic site (Fig. 6 d) was implemented for inhibitor design (Lohse et al. Reference Lohse, Nielsen, Kristensen, Helgstrand, Cloos, Olsen, Gajhede, Clausen and Kristensen2011; Woon et al. Reference Woon, Tumber, Kawamura, Hillringhaus, Ge, Rose, Ma, Chan, Walport, Che, Ng, Marsden, Oppermann, McDonough and Schofield2012), similar to that used for the design of peptide-based LSD1 inhibitors (Szewczuk et al. Reference Szewczuk, Culhane, Yang, Majumdar, Yu and Cole2007). Lohse and colleagues combined a uracil-based group for metal chelation and a minimum H3 peptide required for recognition of H3K9 demethylase to obtain a 5-mer inhibitor, which although showing low potency, exhibited a 4-fold selectivity for JMJD2C (Ki=27 μM) over JMJD2A (Lohse et al. Reference Lohse, Nielsen, Kristensen, Helgstrand, Cloos, Olsen, Gajhede, Clausen and Kristensen2011). Woon and colleagues examined the structure of JMJD2A with a bound H3K9me3-containing peptide (Fig. 6 d), and designed several peptide-based inhibitors by linking the α-KG analog with the peptide through the adjacent position occupied by T11. Subsequently, histone peptide H3(7–14)K9me3 and DNOC (N-oxalyl-d-cysteine), a NOG analog, were linked at peptide position 11 to yield an inhibitor (Fig. 6 p) that could potently inhibit both JMJD2E (IC50 of 90 nM) and JMJD2A (IC50 of 270 nM), with little or no inhibition against other Jumonji subfamily members (Woon et al. Reference Woon, Tumber, Kawamura, Hillringhaus, Ge, Rose, Ma, Chan, Walport, Che, Ng, Marsden, Oppermann, McDonough and Schofield2012). They also proved that a change in peptide substrate could enable these inhibitors to discriminate between closely related Jumonji subfamily members. The combination of both cofactor and substrate structural information should provide a promising avenue for future inhibitor design with enhanced selectivity and potency toward this subfamily of Jumonji KDMs.
3.3 Inhibitors of H3K27 demethylases
UTX and JMJD3 are H3K27-specific demethylases that function in opposition to EZH2 in regulating H3K27 methylation levels, with both KDMs possessing tumor-suppressive characteristics (Pedersen & Helin, Reference Pedersen and Helin2010). UTX is a frequent target of somatic mutations in human cancers, whereas JMJD3 has been demonstrated to activate the INK4A-ARF locus (Pedersen & Helin, Reference Pedersen and Helin2010). An effort championed by the GlaxoSmithKline (GSK) epigenetics group and participation by our group solved crystal structures of JMJD3 with Ni2+ cation, NOG, both in the absence and presence of bound H3K27me3 peptide, as well as potent inhibitors (Kruidenier et al. Reference Kruidenier, Chung, Cheng, Liddle, Che, Joberty, Bantscheff, Bountra, Bridges, Diallo, Eberhard, Hutchinson, Jones, Katso, Leveridge, Mander, Mosley, Ramirez-Molina, Rowland, Schofield, Sheppard, Smith, Swales, Tanner, Thomas, Tumber, Drewes, Oppermann, Patel, Lee and Wilson2012). An extensive screening approach with subsequent optimization identified the lead compound, GSK-J1 (chemical formula in Fig. 7 a), which exhibited an IC50 of 60 nM for JMJD3. In the structure of JMJD3 bound to GSK-J1, the propanoic acid of GSK-J1 mimics α-KG by forming hydrogen bonds with residues within the catalytic site (Figs 7 b and 7 c), whereas the pyridyl-pyrimidine biaryl ring system makes a bidentate interaction with the catalytic metal (Co2+ was used to mimic Fe2+), which induced a 2.34 Å shift of the cation from its original position. Such cation movement resulted in a disruption of the conserved chelating bond between the cation and conserved His1470 (Fig. 7 b). This distinctive dynamic metal shift was unexpected, and introduced a new concept associated with mediating inhibition. Overlaying the structures of inhibitor and peptide bound JMJD3 complexes established that the tetrahydrobenzazepine ring system of GSK-J1 clashes with the Ala29-Pro30 segment of the bound histone peptide (Fig. 7 c), indicating that GSK-J1 may also perturb substrate binding. In addition, the prodrug counterpart, GSK-J4, could efficiently reduce lipopolysaccharide induced proinflammatory cytokine production by human primary macrophages (Kruidenier et al. Reference Kruidenier, Chung, Cheng, Liddle, Che, Joberty, Bantscheff, Bountra, Bridges, Diallo, Eberhard, Hutchinson, Jones, Katso, Leveridge, Mander, Mosley, Ramirez-Molina, Rowland, Schofield, Sheppard, Smith, Swales, Tanner, Thomas, Tumber, Drewes, Oppermann, Patel, Lee and Wilson2012). Such a structure-guided small-molecule and chemoproteomics approach holds promise for the design of inhibitors targeted to other Jumonji family KMTs. Since JMJD3 and UTX function in different pathways, further optimization could facilitate identification of additional inhibitors with the potential for selectively discriminating between these two Jumonji KDM family members.
Inhibitors of H3K27, H3K36/H4K20 and H3K4 demethylases. (a) Chemical formulas of inhibitor GSK-J1. (b) A blow-up view of the catalytic site of JMJD3 with bound GSK-J1 (PDB: 4ASK). Co2+ cation is colored in pink and the JmjC domain is colored in blue. Bound GSK-J1 is shown in a stick representation. (c) Structure of complex of JMJD3 with bound H3(20–34)K27me3 substrate (PDB: 4EZH). Bound histone peptide and α-KG are colored in orange and Ni2+ cation is colored in green. Overlayed GSK-J1 and Co2+ are colored as shown in panel (b). (d) Chemical formula of Daminozide. (e) A blow-up view of the catalytic site of PHF8 with bound Daminozide (PDB: 4DO0). Zn2+ is colored in grey. Daminozide is shown in a stick representation. (f) Chemical formula of PBIT.
3.4 Inhibitors of other histone lysine demethylases
To date, there have been no reports of an H3K79 demethylase, and hence the principles underlying removal of this methylation mark remains a mystery. H3K36-specific demethylases include JMJD2 family members JMJD2A, JMJD2B and JMJD2C, which also possess activity for H3K9 demethylation. Their specific inhibitors have been discussed above in the context of H3K9 KDMs. JHDM1A and JHMD1B (FBXL10/KDM2) are also H3K36 histone KDMs. PHF8 (KDM7) and PHF2 have recently been identified as H4K20-specific demethylases, with aberrant PHF8 closely correlated with prostate cancer (Crea et al. Reference Crea, Sun, Mai, Chiang, Farrar, Danesi and Helgason2012).
Recently, plant growth regulator Daminozide (Fig. 7 d) was found to be a selective α-KG competitive inhibitor against JHMD1B and PHF8 (IC50=2 μM for both), exhibiting >60-fold selectivity over other families of demethylases tested (Rose et al. Reference Rose, Woon, Tumber, Walport, Chowdhury, Li, King, Lejeune, Ng, Krojer, Chan, Rydzik, Hopkinson, Che, Daniel, Strain-Damerell, Gileadi, Kochan, Leung, Dunford, Yeoh, Ratcliffe, Burgess-Brown, von Delft, Muller, Marsden, Brennan, McDonough, Oppermann, Klose, Schofield and Kawamura2012). The structure of PHF8 with bound Daminozide revealed the general mode of inhibition, whereby Daminozide occupied the α-KG binding pocket and recognized the catalytic metal through its acylhydrazide carbonyl and dimethylamino groups (Fig. 7 e). Importantly, both the nitrogen atoms of Daminozide are within the hydrogen bond distance of the hydroxyl group of Tyr257 within the catalytic site. This tyrosine is not present in the catalytic sites of JMJD2A and JMJD3, which may account for the observed selectivity. Although Daminozide itself is not suitable for medicinal use due to its potential human toxicity, its biological impact on animal and other organisms would be of interest and worth further investigation given the high degree of selectivity.
The JARID1 family of proteins belongs to the Jumonji family of H3K4 demethylases. Progress toward development of potent inhibitors for this family of Jumonji KDMs has lagged behind, most likely due to a lack of structural information. Recently, the Yan laboratory identified several novel selective inhibitors of JARID1B, amongst which N-phenyl-benzisothiazolinone (PBIT, Fig. 7 f) exhibited an IC50 of 3 μM (Sayegh et al. Reference Sayegh, Cao, Zou, Morales, Blair, Norcia, Hoyer, Tackett, Merkel and Yan2013). PBIT inhibited several JARID1 subfamily members but not JMJD3 and UTX. Given that its analog ebselen (Fig. 6 o), a zinc-rejecting compound (Sekirnik et al. Reference Sekirnik, Rose, Thalhammer, Seden, Mecinovic and Schofield2009), also inhibited JARID1B (contains several zinc-fingers) with an IC50 of 6 μM, additional research is needed to elucidate the inhibitory mechanism of PBIT.
4. Inhibitors that target non-catalytic sites
All the inhibitors mentioned above target catalytic active sites of KMTs and KDMs through either cofactor-competitive or substrate-competitive mechanisms or both, while allosteric inhibitors that target non-catalytic sites have not been reported for histone KMTs and KDMs. Allosteric inhibitors could potentially trap alternate states through binding to sites distal from the active site, which could constitute a potentially powerful approach toward regulating these histone-modifying enzymes (Hardy & Wells, Reference Hardy and Wells2004). There is however available structure-based information on allosteric inhibition of other epigenetic enzymes that is discussed below.
A case in point is the discovery of an allosteric inhibitor for PRMT3 (Siarheyeva et al. Reference Siarheyeva, Senisterra, Allali-Hassani, Dong, Dobrovetsky, Wasney, Chau, Marcellus, Hajian, Liu, Korboukh, Smil, Bolshan, Min, Wu, Zeng, Loppnau, Poda, Griffin, Aman, Brown, Jin, Al-Awar, Arrowsmith, Schapira and Vedadi2012), a protein arginine methyltransferase (RMT) that catalyzes arginine methylation using SAM as the cofactor, for both histone and non-histone proteins. A screen of 16 000 drug-like compounds led to the discovery and identification of Compound 1 (Fig. 8 a), a potent PRMT3 inhibitor. This small molecule inhibited full-length PRMT3 with an IC50 of 1.6 μM, while showing either weak or no inhibition for several other KMTs and RMTs tested. The structure of PRMT3 bound to Compound 1 revealed that this small molecule bound within a newly created pocket positioned at the dimeric interface that is distal from both the SAM and substrate binding pockets (Fig. 8 b). Complex formation induced conformational side-chain rearrangements at the junction of the α-Y and α-X helices of the dimeric pair, which in turn destabilized the α-X segment of the activation helix. The α-X helix, which is critical for the catalytic activity of PRMT3, was disordered in this complex, accounting for the loss of both cofactor binding and activity. This example demonstrates the potential of allosteric inhibitors for targeting epigenetic enzymes, and highlights a potential alternate approach for targeting histone KMTs and KDMs.
An allosteric inhibitor of PRMT3 and two inhibitors that target non-catalytic sites of MLL1. (a) Chemical formula of Compound 1. (b) Structure of complex of PRMT3 with bound inhibitor Compound 1 (PDB: 3SMQ). The dimeric pairs are colored in blue and salmon respectively. Bound Compound 1 is shown in a space-filling representation. Overlayed structure of PRMT3 with bound SAH is shown in silver. The related α-X and α-Y helices are labeled, with α-Y, which is close to Compound 1, colored in purple. (c) Chemical formula of MI-2-2. (d) A blow-up view of the structure of the Menin binding pocket with bound MI-2-2 (PDB: 4GQ4). Menin is colored in blue and MI-2-2 is shown in a stick representation. (e) Chemical formula of WDR5-0103. (f) Structure of complex of WDR5 with bound compound WDR5-0103 (PDB: 3UR4). WDR5 is colored in blue. Compound WDR5-0103 is shown in a space-filling representation.
In addition to the promise for targeting allosteric sites within histone KMTs and KDMs, additional potential target sites for inhibitor design include regulatory proteins, as shown by the example of MLL1-binding proteins. MLL1 has been extensively studied due to its close association with human leukemia (Krivtsov & Armstrong, Reference Krivtsov and Armstrong2007). However, as discussed earlier, the catalytic SET domain of MLL1 does not constitute an optimal drug target. By contrast, the MLL1 N-terminal region, which interacts with Menin, has been shown to be indispensable for leukemogenic transformation (Yokoyama et al. Reference Yokoyama, Somervaille, Smith, Rozenblatt-Rosen, Meyerson and Cleary2005), and could serve as a potential drug target. Systematic screening for inhibitors that target the Menin-MLL1 interaction resulted in the identification of a specific inhibitor named MI-1. Its optimized derivative MI-2-2 (Fig. 8 c), displayed excellent potency in disrupting the interaction between MLL1 (4–43) and Menin (IC50=520 nM). High resolution structures of complexes of Menin with bound MI-2-2 or MLL1 (4–15) peptide revealed that both MI-2-2 and MLL1 (4–15) peptide occupied the same binding pocket within Menin (Grembecka et al. Reference Grembecka, He, Shi, Purohit, Muntean, Sorenson, Showalter, Murai, Belcher, Hartley, Hess and Cierpicki2012; Shi et al. Reference Shi, Murai, He, Lund, Hartley, Purohit, Reddy, Chruszcz, Grembecka and Cierpicki2012), with MI-2-2 overlapping with the Phe9-Pro13 segment of the MLL1 peptide, indicative a substrate-competitive inhibition mechanism. The bound MI-2-2 is mainly stabilized by hydrophobic interactions, with only a pair of hydrogen bonds formed between N1 and N3 nitrogen atoms of the thienopyrimidine ring and the side chains of a Tyr and an Asn, respectively (Fig. 8 d). MI-2-2 displayed specific and pronounced activity on MLL1 leukemia cells, thereby demonstrating proof of concept that chemical inhibitors can be designed to target and inhibit protein-protein interactions.
Similarly, a screen for inhibitors that targeted the interaction between MLL1 WIN-peptide and WDR5 (Patel et al. Reference Patel, Vought, Dharmarajan and Cosgrove2008) identified small molecule WDR5-0101 (Senisterra et al. Reference Senisterra, Wu, Allali-Hassani, Wasney, Barsyte-Lovejoy, Dombrovski, Dong, Nguyen, Smil, Bolshan, Hajian, He, Seitova, Chau, Li, Poda, Couture, Brown, Al-Awar, Schapira, Arrowsmith and Vedadi2013). Its optimized analog, WDR5-0103 (Fig. 8 e), exhibited a peptide-displacement constant K dis of 3.0 μM. WDR5-0103 bound within the central cavity associated with the arginine-binding pocket of WDR5 mainly through a shape-complementary mechanism (Fig. 8 f), thereby displaying direct competition with the peptide substrate. The inhibitory potency of WDR5-0103 on MLL1 core complexes was dependent on substrate concentration, consistent with the stabilization contribution of WDR5 within the MLL1 complex. WDR5-0103 does not inhibit other histone KMTs tested, indicative of good selectivity. However, since these MLL1 binding partners also interact with other regulatory proteins, the side effects of these compounds require further evaluation, due to their potential for competitive targeting of other binding partners. Examples of such side effects include the interaction between Menin and JUND through the same binding pocket within the MLL1 N-terminus (Huang et al. Reference Huang, Gurung, Wan, Matkar, Veniaminova, Wan, Merchant, Hua and Lei2012), and the interaction between WDR5 with both histone H3 tails and several other MLL family members through targeting of the same binding pocket within the MLL1 WIN motif (Zhang et al. Reference Zhang, Lee, Brunzelle and Couture2012).
5. Future challenges and prospects
One problem that impedes inhibitor design is the low activity of some of the enzymes, given that quite a few enzymes only show detectable activity in the presence of their associated factors, as shown for MLL1 and EZH2, and reconstitution of the functional complexes would be a critical step for inhibitor screening. Another factor that has not yet been fully appreciated reflects the requirement of the nucleosome as the substrate, as shown for the enzymes DOT1L, EZH2 complex, and NSD2, for which histone peptide-based substrates result in undetectable activity or misleading results. In some cases, modified nucleosomes also exhibit enhanced activity, as in the case of H2B120 ubiqutination, which greatly enhanced DOT1L activity (Nguyen & Zhang, Reference Nguyen and Zhang2011). In such circumstances, designer nucleosomes should be used as substrates (Fierz & Muir, Reference Fierz and Muir2012), thereby facilitating the process of inhibitor design at the biologically relevant nucleosomal level.
Recent developments in the successful design of inhibitors that target histone reader modules (Dawson et al. Reference Dawson, Kouzarides and Huntly2012) should draw additional attention to inhibitor design for histone methylated-lysine writers and erasers discussed in this review, given that most histone KMTs and KDMs also possess one or more histone reader modules. Examples include MLL1, which possesses a PHD finger for readout of the H3K4 methylation mark (Wang et al. Reference Wang, Song, Milne, Wang, Li, Allis and Patel2010), and the JMJD2A double Tudor domain, which recognizes both H3K4 and H4K20 methylation marks (Huang et al. Reference Huang, Fang, Bedford, Zhang and Xu2006; Lee et al. Reference Lee, Thompson, Botuyan and Mer2008). The histone reader modules of these enzymes are critical for their in vivo enzymatic activity, given that recruitment of these enzymes to the proper loci constitutes an important step toward their in vivo function. Targeting such reader modules within histone KMTs and KDMs offer additional opportunities for inhibitor design.
A major bottleneck in inhibitor design has involved the search for and identification of target-specific hits from within an unlimited choice of candidates. As shown above, the majority of successful studies have made extensive use of structural biology, which although dispensable at the initial hit stage, has been of invaluable assistance in the small molecule optimization stage, and has also been considered as an essential contributor for the fragment-based drug design approach (Murray & Blundell, Reference Murray and Blundell2010). Structural biology-based knowledge of enzyme binding pockets has also contributed to in silico approaches to inhibitor design and optimization (Ghosh et al. Reference Ghosh, Nie, An and Huang2006). We anticipate that the combined improvements in screening technology, medicinal chemistry, structural and computational biology, and cellular-based analysis should greatly expedite the process of pharmacologically potent inhibitor identification and optimization targeted toward these disease-impacted epigenetic modifying enzymes, which in turn should ultimately lead to the application of efficacious and safe drugs beneficial for human health.
6. Acknowledgements
This research was supported by the Leukemia and Lymphoma Society and the STARR Foundation to D. J. P., as well as the ‘Thousand Young Talents’ Program of China and start-up funds from Beijing Normal University to Z. W.
7. Competing financial interests
The authors declare no competing financial interests.

