|Home | About | Journals | Submit | Contact Us | Français|
The JmjC histone demethylases (KDMs) are linked to tumour cell proliferation and are current cancer targets; however, very few highly selective inhibitors for these are available. Here we report cyclic peptide inhibitors of the KDM4A-C with selectivity over other KDMs/2OG oxygenases, including closely related KDM4D/E isoforms. Crystal structures and biochemical analyses of one of the inhibitors (CP2) with KDM4A reveals that CP2 binds differently to, but competes with, histone substrates in the active site. Substitution of the active site binding arginine of CP2 to N--trimethyl-lysine or methylated arginine results in cyclic peptide substrates, indicating that KDM4s may act on non-histone substrates. Targeted modifications to CP2 based on crystallographic and mass spectrometry analyses results in variants with greater proteolytic robustness. Peptide dosing in cells manifests KDM4A target stabilization. Although further development is required to optimize cellular activity, the results reveal the feasibility of highly selective non-metal chelating, substrate-competitive inhibitors of the JmjC KDMs.
Eukaryotic gene expression is substantially regulated by a complex set of covalent modifications to DNA and histones1,2. Aberrant histone modifications are associated with diseases, including cancer, inflammatory disorders, neurological and cardiovascular diseases3,4. The multi-domain proteins that read and write posttranslational histone modifications are often mutated/dysregulated in such diseases. Structural similarities between individual domains of these ‘epigenetic' proteins make the identification of potent and selective inhibitors of these proteins challenging.
Histone lysine methylation is dynamically regulated by methyltransferases and demethylases (KDMs). Two families of human KDMs catalyse the removal of methyl groups from methylated lysines on histone tails in a sequence- and methylation-state-dependent manner5; the flavin-dependent lysine-specific demethylases (KDM1/LSD)6 and the JmjC-domain containing KDMs (JmjC-KDMs), which are Fe(II) and 2-oxoglutarate (2OG)-dependent oxygenases (KDM2–7)7,8. Only a few JmjC-KDM inhibitors are reported, the majority of which operate via active-site metal chelation9. Achieving both selectivity and potency for specific JmjC-KDMs remains a major challenge9.
Natural product screening is an attractive method for identifying inhibitor hits and leads. However, natural product-based structure activity studies are often limited by difficulties in synthesis and/or systematically manipulating the output of biosynthesis. An alternative approach, Random nonstandard Peptides Integrated Discovery (RaPID) system has been developed, which involves efficient selection using ribosomally expressed ‘natural product'-inspired macrocycles10,11. De novo macrocycles that tightly bind to target proteins can be efficiently selected from the >1012 members of the library; derivatives of the initial hits are then chemically synthesized for structural optimizaion11,12,13.
We describe the use of the RaPID methodology for discovery of highly selective and potent cyclic peptide inhibitors of KDMs, which, after structure- and activity-guided modifications, show evidence of on-target engagement in cells. We targeted the KDM4 subfamily, which represent biomedically attractive but challenging targets. Although the catalytic domains (JmjC-domain) and active sites are highly conserved, all KDM4s remove the repressive H3K9me3 mark, but only KDM4A-C are additionally capable of demethylating the activating H3K36me3 mark14,15. Intra-subfamily selective inhibitors will be useful tools to dissect the roles of the opposing histone modifications and of the KDM4 isoforms in disease.
A messenger RNA template library was designed with the general form AUG-(NNK)4–12-UGC, where the AUG start codon was reassigned from Met to either N-chloroacetyl-L-Tyr or N-chloroacetyl-D-Tyr, which induce cyclization of the translated peptides through reaction with an appropriately positioned cysteine (UGC) downstream of the random sequence to yield lantibiotic-like thioether rings16,17. The selection used two thioether-macrocycle libraries (LTyr and DTyr libraries)12 with screening against the catalytic domain of His-tagged KDM4A (residues 1–359) (Fig. 1 and Supplementary Fig. 1a). Enriched complementary DNA pools from both libraries were cloned and analysed. Out of 21 (DTyr-library) and 23 (LTyr-library) clones sequenced (Supplementary Fig. 1b,c), five cyclic sequences (CP1–5, 11–14 residues in length) were selected for further analysis based on hit frequency and similarity; these were synthesized by solid-phase peptide synthesis. Catalytic inhibition of KDM4A/C by the selected peptides was tested using a luminescence-based AlphaScreen activity assay18. From the KDM4A screen, two cyclic peptides (CP2 and CP4) demonstrated potencies of half-maximal inhibitory concentration (IC50)<50nM against KDM4A (Table 1); CP5 was less potent (IC50<500nM). The other two hits (CP1 and CP3) were much less effective (IC50>100μM).
Clear selectivity (that is, >160-fold by IC50) of all inhibitory CPs (CP2, CP4 and CP5) was observed for KDM4 over representatives from four other human JmjC-demethylase families and two related 2OG oxygenases, prolyl hydroxylase domain 2 and factor inhibiting HIF, as well as a flavin-dependent KDM (KDM1A), indicating that these CPs are exceptionally selective (Table 1). Intra-subfamily selectivity was then investigated by determining IC50 values across all KDM4 subfamily members. CP2 and CP4 displayed high potency (IC50<50nM) against KDM4A–C; importantly, they were much less active against KDM4D–E, revealing intra-subfamily selectivity of >100-fold and that the CPs are therefore even capable of distinguishing the subtle differences between KDM4A–C and KDM4D/E. Binding studies demonstrated that CP2, CP4 and CP5 are tight binders of KDM4A with slow on and off rates; in contrast, the histone H3K9 peptide (KDM4A substrate) exhibits substantially faster on/off rates and weaker binding affinity (Kd>1,000-fold; Table 1 and Supplementary Fig. 2a,b). Kinetic analyses with CP2 revealed it is a competitive inhibitor with respect to H3K9me3 (Supplementary Fig. 2c). CP2 also inhibits the demethylation activity of full-length KDM4A (residues 1–1,064) purified from HEK293T cells (Supplementary Fig. 3), indicating CP2 can inhibit KDM4A when other domains (for example, PHD/tudor domains) are present.
To investigate the structural basis of the unprecedented selectivity of the cyclic peptides, we crystallized CP2 with KDM4A (2.7Å resolution, Fig. 2a). Interestingly, the structure reveals that CP2 occupies the histone-binding groove of KDM4A, but not the 2OG-binding pocket, the latter of which has been exploited for most reported 2OG oxygenase, including JmjC KDM, inhibitors. Moreover, CP2 exhibits a binding mode, which partially overlaps, but which is clearly distinct from those observed for H3K9me(n) and H3K36me(n) histone substrates (Fig. 2b). Analysis of the inhibitor conformation revealed that CP2 adopts a distorted β-sheet fold with two turns; a type 1 β-turn (Thr5-Gly8) is closest of the two β-turns to the active site metal (Fig. 2c). CP2 engages in multiple hydrogen bonds with KDM4A, many of which are not observed in the reported structures of KDM4A-H3K9me3 (PDB 2OQ6) and KDM4A-H3K36me3 enzyme substrate complexes (PDB 2P5B)19. Strikingly, the structure reveals Arg6-CP2 bound in the sub-pocket usually occupied by the extended trimethyl-lysine side chain of the H3K9 and K36 substrate residues. In addition to mimicking the positive charge of the lysine, the guanidino group of Arg6-CP2 is positioned to make hydrogen bonds with KDM4A residues (Tyr177, Ser288 and Asn290; Supplementary Fig. 4a,b). Two KDM4A residues in the CP2-complexed structure (Tyr175 and Arg309) move, relative to uncomplexed KDM4A, to open pockets that accommodate CP2, suggesting flexibility in these regions. The Arg309 side chain rotates by approximately 90° relative to the H3K9me3/H3K36me3 bound KDM4A structures to accommodate Trp9-CP2 (Fig. 2d,e). Tyr175 also shifts by approximately 90° to accommodate Arg6-CP2 (Fig. 2f). This movement relocates an adjacent loop (between α6 and DSBH βI (β7), residues 161–173), which resembles a conformation found in KDM4C(apo) structure (Fig. 2f). Other notable conformational changes in the KDM4A residues induced by CP2 binding include at Lys241 and Asp311 (Supplementary Fig. 5).
Although the linearized sequence of CP2 maintains modest activity against KDM4A, its IC50 is four fold higher than that of CP2, probably due to the increased conformational freedom of the linear sequence (Table 2, peptide 2). However, changing the stereochemistry of the ‘cyclizing' amino-terminal tyrosine (D/L-Tyr) does not affect the potency of CP2 (Table 2, peptide 3); this may reflect the lack of interactions that CP2 makes with KDM4A in the regions surrounding the thioether bond (Tyr1, Val2 and Tyr12-Cys14) as observed crystallographically (Fig. 2c and Supplementary Fig. 4a,b).
Analysis of the polar contacts in the KDM4A-CP2 structure provides a structural basis for the selectivity of the cyclic peptides for KDM4A-C over D/E (Supplementary Fig. 4c). Despite the sequence and structural similarity in the catalytic domains between KDM4A and KDM4D (77%)20, 5 of the 14 residues directly interacting with CP2 in KDM4A–C (Asn86, Gln88, Ser288, Arg309 and Asp311) are different in KDM4D/E. Interestingly, three of the five residues (Asn86, Gln88 and Arg309) coincide with those identified as important for the lack of activity of KDM4D/E on H3K36me3 (ref. 20).
To investigate the role of Arg6 of CP2, which binds in a similar manner to the H3Kme(n) residue, we substituted it with other residues (Table 2 and Supplementary Fig. 6); substitution with alanine (peptide 4) or with hydrophobic/uncharged side chains including N--acetyl lysine and citrulline (peptides 5, 6 and 7) leads to loss of potency (IC50>2,400nM for KDM4A). Inhibition of activity is regained upon substitution with positively charged residues (lysine (peptide 8)/N--trimethyllysine (peptide 9), IC50=24 and 12nM, respectively, for KDM4A). The relevance of the crystallographically observed orientation of Arg6 in solution was explored using the trimethyl-lysine CP2 variant (peptide 9). Despite the lack of similarity between CP2 and the histone substrates, N--trimethyl-lysine-CP2 (R6Kme3) was efficiently demethylated by KDM4A to Kme2/Kme1/Kme0 (Fig. 3b). A co-crystal structure of CP2(R6Kme3) in complex with N-oxalylglycine (NOG, an inactive form of 2OG) and KDM4A supports the productive binding conformation of the CP2(R6Kme3) (Fig. 3a); the acceptance of CP2(R6Kme3) as a KDM4A substrate unequivocally demonstrates that the CP2 series does not compete with 2OG. Notably, KDM4A was also able to demethylate methylated Arg6 (CP2(R6me2a), peptide 10) to the non-methylated, inhibitory CP2 sequence (Fig. 3c). The different sequences and binding modes between the well-established histone substrates (H3K9me3/K36me3) and CP2 suggests that KDM4s might be more promiscuous in their substrate selectivity and binding mode than presently perceived.
We then investigated the utility of CP2 as a cellular probe. Strain-promoted azide-alkyne cycloaddition was used to generate CP2 derivatives with a fluorescein dye (CP2(Fl)) (Supplementary Fig. 6a)21,22. Time-dependent increases in intracellular fluorescence were observed in HeLa cells dosed with CP2(Fl), indicating cellular uptake (Supplementary Fig. 7). The fluorescence signal was observed to accumulate outside the nucleus (>10h), possibly indicating peptide degradation. Cellular thermal shift assays (CETSA)23 were then used to investigate intracellular target engagement of CP2 in cells. U2OS cells ectopically expressing FLAG-tagged KDM4A1–1,064 were dosed with CPs for 20h, harvested, washed and the intact cells were subjected to a temperature gradient. CETSA melt curves were determined for KDM4A and actin by western blot analysis. The melting temperature (Tm) for FLAG-KDM4A was shifted from 47°C (dimethyl sulfoxide, DMSO) to 59°C with CP2(T13Z) (the precursor to the labeled CP2(Fl), peptide 25) (Fig. 5b and Supplementary Fig. 6b). Isothermal dose–response function (ITDRF) at 55°C revealed FLAG-KDM4A stabilization by CP2 in a dose-dependent manner (EC50~1μM), supporting intracellular target engagement of the inhibitor (Fig. 5c). Correlating with the selective inhibition of CP2 for isolated KDM4A-C, FLAG-KDM4E was not stabilized under the same CETSA conditions (Supplementary Fig. 8). However, despite apparent intracellular KDM4A engagement, CP2-treated cells did not have any apparent effect on the global histone methylation levels at H3K9me3 and H3K36me3 mark (Supplementary Fig. 12).
Addition of a poly-arginine sequence to CP2 aimed at enhancing the cellular uptake (CP2(polyR), peptide 11) unexpectedly increased the potency against KDM4s in vitro (KDM4A IC50=1.8nM, KDM4C IC50=0.8nM; Table 2). Interestingly, polyR alone is a potent KDM4A inhibitor (IC50=40nM); thus, the increased potency of CP2(polyR) is likely to be a combined effect of the two inhibitory elements. However, although cytotoxicity was observed at high concentrations (>3μM) with significant reduction in cell numbers, no inhibition of cellular KDM4A demethylase activity by CP2(polyR) was detected (Supplementary Figs 11 and 13). An analogous phenomenon has been previously reported with disulphide linked cyclic peptide generated against KDM4C using phage display24; the potency of a proposed allosteric binding cyclic peptide inhibitor (IC50=52μM) was improved to IC50=0.6μM on addition of a poly arginine/lysine (TAT) tag, but no cell activity was observed.25
We then modified CP2 by backbone amide N-methylation, which is reported to improve peptide stability and cellular uptake.25 Guided by the CP2-KDM4A crystal structure (Fig. 2), we synthesized mono N-methylated variants at backbone amide positions not engaging in critical hydrogen bonding interactions with KDM4A and within the CP2 β-sheet. The selected residues, Cys14, Thr13, Val2 and Tyr1, were mainly in the non-interacting, thioether-linkage region of CP2 (Table 2 and Fig. 4a; peptides 13–16). Mass spectrometric (MS) analysis of CP2 degradation patterns after incubation in cell lysates confirmed that these sites are prone to proteolysis (Fig. 4b and Supplementary Fig. 9). D-Ala was also substituted for Gly8 of CP2 with the aim of improving proteolytic stability, as the crystal structure predicted this substitution would be tolerated and hydrolysis at this site was observed by MS (Fig. 4b). 4-Fluorophenylalanine was further incorporated with the aim of improving cellular uptake by increasing hydrophobicity (peptides 17 and 18). Although some of the CP2 modifications reduced activity, in general the structure-guided modifications were well tolerated. Combinations of tolerated modifications were prepared (named as CP2.1, CP2.2 and CP2.3, peptides 19–21, respectively) (Table 2, Fig. 5a and Supplementary Fig. 6); Potency and selectivity for KDM4A over other KDM families were maintained for these modified CPs with isolated enzymes. (Supplementary Table 1). Interestingly, the inhibitory effect of CP2.3 against KDM4B was much weaker than for CP2, but was retained for KDM4A/C. CP2.3 demonstrated greater stability in HeLa cell extracts (t1/2~5h) compared with CP2 (t1/2~1h), as determined by MS (Supplementary Fig. 10). Furthermore, the modified CP2.3 was able to stabilize FLAG-KDM4A as observed in CETSA, at concentrations lower than the unmodified CP2 (EC50<100nM) (Fig. 5c), suggesting that the target engagement was maintained on modification and that peptide cellular stability was increased.
We tested CP2.1, CP2.2 and CP2.3 for KDM inhibition in HeLa cells by immunofluorescence analysis of histone methylation marks (72h treatment). Although application of the CP2 parent gave no discernible changes, treatment with CP2.1–3 manifested increases in nuclear H3K9me3/H3K36me3 staining over the 10–20μM inhibitor range (Supplementary Fig. 12a–g). No changes were observed in the apparent H3K9me3 levels when cells were treated with the inactive CP2.3(R6A) peptide (peptide 22; Table 2 and Supplementary Fig. 12f). A significant reduction in cell numbers was observed in HeLa cells treated with CP2.1–3 (>10μM) (Supplementary Fig. 12e). As hypermethylation can correlate with cellular stress/toxicity, it is unclear whether the observed hypermethylation is due to on-target KDM4 inhibition or cytotoxic effects, or a combination of these. However, dosing HeLa cells ectopically overexpressing FLAG-KDM4A (wild-type or catalytic inactive mutant H188A) with CP2.3 for 24h showed no significant changes in H3K9me3 levels (Supplementary Fig. 13) at non-cytotoxic concentration range (<10μM). Thus, the biological relevance of the cellular observations of changes in global histone methylation status following treatment with the KDM inhibitors should be treated with caution.
In conclusion, we have used in vitro selection from a ribosomally synthesized library of cyclic peptides to identify natural product-like inhibitors of KDM4A-C, which act via a previously unidentified binding mode and which have unprecedented selectivity and potency. The RaPID display approach is substantially more efficient than traditional medicinal chemistry and is likely to be of widespread utility in target-based probe discovery. The method is well-suited to identify new inhibitor binding modes, as revealed by the structures of KDM4A complexed with CP2 and CP(R6Kme3), and associated biochemical results. The binding mode of CP2 is distinct from reported KDM4C peptide inhibitors (with IC50 values in the μM range) based on the outputs of a phage display library screen, which probably do not bind at the active site (structural studies are not available)24. The sequence of CP2 is clearly distinct from that of well-characterized histone substrates for KDM4A–C. The importance of the anchoring residue Arg6 within the CP2 sequence for potent KDM4A inhibition, suggests that arginine residues can compete with methylated lysines binding to KDM4A. This is significant, given the recent findings that some, but not all, JmjC-KDMs, including some KDM4 subfamily members, can also act as N-methyl arginine demethylases26. KDM4A can both catalyse the demethylation of CP2(R6Kme3) and CP2(R6Rme2a), revealing its potential for modulating both lysine and arginine methylation status in non-histone sequences. Structure- and MS-guided derivatization of an initial hit led to proteolytically more stable macrocycles that stabilize the intracellular target. Changes in the global histone H3K9/K36 methylation status were observed in cells treated with CP2.3 for 72h, although further studies are needed to link this observation directly to KDM4 inhibition. The biological effects of substrate competing cyclic peptides targeting KDM4A may arise both from simple inhibition of KDM4A catalysis, as well as disruption of its binding to chromatin complex; the latter mechanism is likely to be different to the effect of catalytically inactivating small-molecule inhibitors, which compete with the 2OG cosubstrate, but not with histone substrates. We also cannot rule out the possibility that the observed intracellular KDM4A stabilization may in part be due to CP2/CP2.3 proteolytically cleaved fragments, which may, or may not, exert inhibitory effects. Although further studies are necessary to investigate the cell permeability/subcellular distribution, stability and activities of the cyclic peptides, our work highlights the utility of the approach for the generation of potent, substrate competitive and selective KDM inhibitors, and more generally, for developing new approaches to the inhibition of challenging targets.
Flexizyme RNA and tRNAfMetCAU were prepared as reported16. ClAcLTyr-tRNAfMetCAU and ClAcDTyr-tRNAfMetCAU were synthesized by using the reported flexizyme technology11. Translation factors, enzymes and ribosome were prepared and mixed as previously described27,28,29 to reconstitute an in vitro translation system used for reprogramming of translation initiation11,17. The translation reaction mixture contained final concentrations of 50mM Hepes-KOH (pH 7.6), 100mM potassium acetate, 2mM GTP, 2mM ATP, 1mM CTP, 1mM UTP, 20mM creatine phosphate, 12mM Mg(OAc)2, 2mM spermidine, 2mM dithiothreitol, 1.5ml−1 Escherichia coli total transfer RNA (Roche), 1.2μM ribosome, 0.6μM MTF, 2.7μM IF1, 0.4μM IF2, 1.5μM IF3, 30μM EF-Tu, 30μM EF-Ts, 0.26μM EF-G, 0.25μM RF2, 0.17μM RF3, 0.5μM RRF, 4μgml−1 creatine kinase, 3μgml−1 myokinase, 0.1μM pyrophosphatase, 0.1μM nucleotide-diphosphatase kinase, 0.1μM T7 RNA polymerase, 0.73μM AlaRS, 0.03μM ArgRS, 0.38μM AsnRS, 0.13μM AspRS, 0.02μM CysRS, 0.06μM GlnRS, 0.23μM GluRS, 0.09μM GlyRS, 0.02μM HisRS, 0.4μM IleRS, 0.04μM LeuRS, 0.11μM LysRS, 0.03μM MetRS, 0.68μM PheRS, 0.16μM ProRS, 0.04μM SerRS, 0.09μM ThrRS, 0.03μM TrpRS, 0.02μM TyrRS, 0.02μM ValRS and 200μM each proteinogenic amino acids, except for methionine, and 50μM ClAcLTyr-tRNAfMetCAU or ClAcDTyr-tRNAfMetCAU.
RNAs consisting of 4−12 repeated NNK random sequences (5′-GGGUU, AACUU UAAGA AGGAG AUAUA CAU AUG (NNK)m UGC GGC AGC GGC AGC GGC AGC UAG GACGG GGGGC GGAAA-3′, m=4−12) were prepared by in vitro transcription according to the reported method12. The resulting RNAs were mixed in the following ratio—(NNK)4:(NNK)5:(NNK)6:(NNK)7:(NNK)8:(NNK)9:(NNK)10:(NNK)11:(NNK)12=20−3:20−2:20−1:1:10:10:10:10:10. The mRNA library was ligated with a puromycin linker (5′-CTCCC GCCCC CCGTC C-(SPC18)5-CC-puromycin-3′) by T4 RNA ligase. The ligated product was purified by phenol–chloroform extraction and ethanol precipitation.
Translation of the first round selection was performed using 156pmol mRNA-puromycin and 150μl of translation mixture at 37°C for 30min. Subsequently, the translation mixture was incubated at room temperature for 12min to conjugate the translated peptide with the corresponding mRNA–puromycin. To this solution was added 15μl of 200mM EDTA (pH 7.5) then further incubated at 37°C for 30min, to proceed peptide macrocyclization. After exchange with the selection buffer (50mM Tris-HCl (pH 7.6), 150mM NaCl, 0.05% Tween 20) by gel filtration, 165μl of blocking solution (100mM Tris-HCl (pH7.6), 300mM NaCl, 0.2% acetyl BSA) was added. The resulting peptide library was treated with Dynabeads magnetic beads (His-tag isolation and pulldown, M-280, Invitrogen) at 4°C for 30min. The supernatant was then collected, mixed with KDM4A-immobilized magnetic beads (final concentration of KDM4A was adjusted to 1μM) and incubated at 4°C for 30min. After the supernatant was removed, the bead was washed three times with 800μl of cold selection buffer. To the bead was added 40μl of room temperature reaction buffer I (50mM Tris-HCl (pH 8.3), 75mM KCl, 3mM MgCl2, 10mM dithiothreitol, 0.5mM dNTPs, 2μM CGS3an13.R39 (5′-TTTCC GCCCC CCGTC CTAGC TGCCG CTGCC GCTGC CGCA-3′)) containing 200 units of M-MLV reverse transcriptase (Promega) and 8 units of RNase inhibitor (Promega), and reverse transcribed at 42°C for 60min. The collected cDNA was eluted with 800μl of PCR buffer (10mM Tris-HCl (pH 7.5), 50mM KCl, 0.1% Triton X-100, 2.5mM MgCl2, 0.25mM dNTPs, 0.25μM T7g10M.F48 (5′-TAATA CGACT CACTA TAGGG TTAAC TTTAA GAAGG AGATA TACAT ATG-3′), 0.25μM CGS3an13.R39) at 95°C for 5min. After addition of Taq DNA polymerase to the eluate, the mixture was used for PCR amplification. The amplified DNA was purified by phenol/chloroform extraction followed by ethanol purification. The resulting DNA library was transcribed in vitro, ligated with the puromycin linker and used for the next round of selection.
For the 2nd–6th rounds of selection, the translation reaction (5μl scale) was carried out (37°C for 30min, then room temperature for 12min). After incubation with 1μl of 100mM EDTA (pH 7.5) at 37°C for 30min, 3.43μl of room temperature reaction buffer II (74mM Tris-HCl (pH 8.3), 44mM Mg(OAc)2, 0.74mM dNTPs, 29mM KOH, 4.9μM CGS3an13.R39, 15Uμl−1) M-MLV reverse transcriptase without RNase H activity (Promega)) was added. Reverse transcription of the mRNAs fused on peptides was then performed at 42°C for 60min. After addition of 21μl of selection buffer to the library solution and subsequent buffer exchange with selection buffer by gel filteration, 30μl of blocking solution was added. The resulting peptide library was applied to free Dynabeads magnetic beads at 4°C for 30min and then the supernatant was recovered (negative selection). After repeating the negative selection step 3 times, half of the recovered peptide solution was mixed with KDM4A-immobilized magnetic beads (final concentration of KDM4A was adjusted to 500nM) and incubated at 4°C for 30min. After the supernatant was removed, the beads were washed three times with 60μl of cooled selection buffer. The collected cDNA–mRNA–peptide fusion was recovered with 100μl of PCR buffer by incubation at 95°C for 5min and amplified by Taq DNA polymerase. The amplified DNA was purified by phenol/chloroform extraction followed by ethanol purification. The resulting DNA library was transcribed in vitro, ligated with the puromycin linker, then used for the next round of selection. After the fifth and sixth rounds of the selection, the resulting cDNA libraries were cloned into the pGEM-T Easy Vector (Promega) and the individual clones were sequenced.
PolyR (Arg)9 peptide was purchased from AnaSpec (AS-621204). The synthesis of peptides used in this study are as described:
Amino acids were from CS Bio, Novabiochem, Sigma, TCI, Alfa Aesar and AGTC Bioproducts. Azide-containing amino acids were from Iris Biotech, BCN-conjugates for copper-free click chemistry were purchased from Synaffix (Nijmwegen, The Netherlands), Azides for click chemistry were obtained Baseclick (Tutzing, Germany). All peptides were analysed using an Agilent 1200 series LC-MS system (6120 quadrupole MS) with a Waters Sunfire column. Preparative HPLC purifications were carried out using Shimazu prominence LC-20AP system with a Merck Chromolith Prep column (200–25mm) or Dionex Ultimate 3,000 system with a Grace Vydac 218TP101522 column.
Method A: peptides were synthesized with an amidated carboxy terminus on a 25μmol scale using NovaPEG Rink Amide resin (Merck) using standard Fmoc protection chemistry. After coupling of the N-terminal L-tyrosine or D-tyrosine residue, their Fmoc group was removed. The resulting N-terminal α-amino group was chloroacetylated by incubating with a solution of 0.5M chloroacetyl N-hydroxysuccinimide ester chloromethylcarbonyloxysuccinimide (ClAc-OSu) in N-methylpyrrolidone with rotation for 40min at room temperature. After washing the resin with dimethylformamide (2ml, three times) and dichloromethane (2ml, six times), the peptides were cleaved from the resin and deprotected by incubation with a solution of trifluoroacetic acid (TFA)/1,2-ethanedithiol/triisopropyl silane/water (92.5:2.5:2.5:2.5) with rotation at room temperature for 3h. The solutions were concentrated in vacuo and the cleaved peptides were then precipitated with diethyl ether. The resulting linear peptide pellets were dissolved in 20ml of water/DMSO (1:1) and triethylamine was added to the solution to give a pH of around 10. After incubation at room temperature for 30min to effect cyclization, the peptide solution was acidified using TFA then purified by reverse-phase HPLC (Shimazu prominence LC-20AP system with a Merck Chromolith Prep column) in 0.1% aqueous TFA/acetonitrile containing 0.1% TFA gradient. The purified peptides were lyophilized and dissolved in DMSO.
Method B: for the synthesis of cyclic peptides, the linear precursors were prepared by standard solid-phase synthesis using a CS Bio CS336X peptide synthesizer (100μmol scale) using DIC as coupling reagent. For methylated peptides coupling times were tripled for reaction of the amino acid following the N-methylated one.
After cleavage of the N-terminal Fmoc-protecting group a solution of 150mg of ClAc-OSu in 4ml dimethylformamide was added to the resin and the mixture was shaken for 3h. The resin was filtered off and subsequently treated with 4ml of deprotection solution TFA/triisopropyl silane/water (95:2.5:2.5). After 3h, the volume was reduced to 1ml under a nitrogen stream and the peptides were precipitated with cold Et2O. The mixture was centrifuged and the supernatant discarded. The solid was taken up in 1.5ml of triethylammoniumacetate buffer (1M, pH 8.5) and the pH readjusted to >8 if necessary. In a microwave (Biotage Initiator), the mixture was heated to 80°C for 10min and subsequently purified by HPLC (0–45% MeCN in 45min, 0.1% TFA, Dionex Ultimate 3,000 series, Grace Vydac 218TP101522 column).
ClAc-OSu was synthesized using a modified literature procedure30. A round-bottom flask was charged with 4.00ml chloroacetylchloride, 6.90g of K2CO3 and 80ml CH2Cl2, and cooled to 0°C. N-hydroxysuccinimide (5.78g) was dissolved in 80ml of MeCN and added slowly to the reaction vessel. The mixture was allowed to warm to room temperature overnight. Subsequently, the slurry was diluted with 500ml CH2Cl2, washed with 200ml HCl (1M) and the aqueous layer was then counter-extracted with 500ml CH2Cl2. The combined organic layers were dried over MgSO4, then concentrated to give 3.50g (36%) of chloromethylcarbonyloxysuccinimide as a colourless solid. 1H NMR (400MHz, CDCl3): δ (p.p.m.)=4.37 (s, 2 H), 2.85 (s, 4 H). 13C NMR (101MHz, CDCl3) δ (p.p.m.)=168.4, 163.2, 37.8, 25.5.
N-terminally Hexa-His-tagged KDM4A (His-KDM4A1–359) was produced and purified as described31. The KDM4A gene encoding for KDM4A protein (residues 1-359) was cloned into the pNIC-Bio3 plasmid, which incorporates a C-terminal AviTag and expressed with biotin supplement in BL21(DE3)Rosetta co-transfected with pCDF-BirA32. His-KDM4A with biotin (His-KDM4A-Bio) was purified using nickel affinity chromatography. Complete (>95%) biotinylation of KDM4A was confirmed using electrospray ionization–MS (Waters Micromass LCT Premier). The enzyme activity was comparable to that of His-KDM4A and the biotin-tag was demonstrated to be accessible to streptavidin binding as determined by AlphaScreen-binding assay33. The full-length FLAG-tagged KDM4A1–1,064 (ref. 15) was produced in HEK293T cells. HEK293T cells were grown to 80% confluency in DMEM medium (Sigma) supplemented with Glutamax (Life Technologies), 10% fetal bovine serum (FBS, Sigma). FLAG-tagged KDM4A was transfected using polyethyleneimine. Cells were harvested after 48h, washed with PBS (Sigma) and lysed (45min, 4°C) in 50mM HEPES (pH 7.5), 150mM NaCl buffer containing protease inhibitor cocktail (Sigma), 10μM DNase (Sigma) and 0.1% NP40 (Roche). FLAG-tagged KDM4A was purified by immunoprecipitation using anti-FLAG M2 magnetic beads (Sigma) and washed in in 50mM HEPES (pH 7.5) 150mM NaCl. The final wash was performed in 50mM HEPES (pH 7.5) and the resin was used as an enzyme source for the full-length KDM4A.
IC50 values for KDM2–6 were measured using an AlphaScreen assay as described18. The appropriate KDM (5μl) was pre-incubated with an inhibitor (0.1μl in DMSO) or DMSO for 15min before the addition of 5μl mixture of biotinylated peptide substrates, 2OG, ferrous ammonium sulfate (1–10μM final) and ascorbate (100μM final) to initiate the demethylation reaction. The reaction was quenched with 5μl of 30mM EDTA at predetermined reaction times for each enzyme (initial rate). Demethylated product was detected by addition of 5μl of AlphaScreen donor and acceptor beads (AlphaScreen General IgG Detection Kit, 0.02mgml−1 final each) pre-incubated with antibody. The plates were read after 1h incubation in a BMG Labtech Pherastar FS plate reader. Data were normalized to a no enzyme control and IC50 values were calculated using nonlinear regression curve fits using GraphPad Prism. Assay buffer (50mM HEPES pH 7.5, BSA (0.1% w/v) and Tween-20 (0.01% v/v)) was used throughout, except for dilution of stock ferrous ammonium sulfate solutions where water was used. The enzyme, substrate and antibody concentrations used are described in Supplementary Table 3. A fluorescence-based LSD1 inhibitor screening assay kit (Cayman, 700120) was used for KDM1/LSD1 assay. Assays for prolyl hydroxylase domain 2 and factor inhibiting HIF using matrix-assisted laser desorption/ionization–time of flight MS were performed as previously described31. Prism 5.0 (GraphPad) was used to generate and to fit dose–response curves. IC50 values were determined by plotting percentage inhibition versus compound concentration relative to positive (vehicle) and negative (no enzyme) controls, using normalized dose–response (variable slope) fit. The reported values are the mean of multiple experiments (n≥3). The formaldehyde dehydrogenase coupled KDM assay was used to determine the mode of inhibition34, using 0.5μM His-KDM4A1–359. Prism 5.0 (GraphPad) was used to fit in mix-mode inhibition, to determine the kinetic parameters. Binding constants for CPs and H3(1–15)K9me3 were measured using biolayer interferometry (OctetRed 384, ForteBio) in buffer (50mM HEPES (pH 7.5), 150mM NaCl, 10μM ammonium iron(II) sulphate). His-KDM4A1–359-Bio (70nM) was immobilized onto super streptavidin biosensor and subjected to different peptide concentrations in buffer. Data analysis was performed using ForteBio Data Aalysis (v6.4).
The KDM4A.CP2 complex was prepared by adding CP2 to a final concentration of 0.2mM in a solution of 0.25mM (11mgml−1) His-KDM4A1–359 in gel filtration buffer (10mM HEPES pH 7.5, 500mM NaCl and 5% glycerol). Crystals were grown by vapour diffusion at 4°C in 200nl sitting drops with a 1:1 ratio of sample to well solution (0.1M propionate-cacodylate-bis-tris propane pH 9.0 and 25% w/v polyethylene glycol 1,500). Crystals of KDM4A.CP2(R6Kme3) complex (0.25mM His-KDM4A1–359, 1mM NiCl2, 2.5mM NOG and 5mM CP2(R6Kme3) were grown by vapour diffusion at 4°C in 300nl sitting drops with a 1:2 ratio of sample to well solution (0.1M bicine pH 9.0 and 10% w/v polyethylene glycol 6,000). Crystals were cryo-protected by transferring to a solution of mother liquor supplemented with 25% v/v glycerol before being flash frozen in liquid N2. An additional 0.5mM CP2 was added to the cryoprotectant for KDM4A.CP2 crystals.
Data were collected from single crystals at 100K at the Diamond Light Source, beamline I04-1. Data for the KDM4A.CP2 complex were processed using XDS35 and SCALA36; data for the CP2(R6Kme3) complex were processed using HKL2000 (ref. 37). The structures were solved by molecular replacement using PHASER38 and PDB 2OX0 as the initial model. The majority of the KDM4A.CP2 structure refinement was carried out with Phenix39, with iterative rebuilding of the model using COOT40. Final rounds of slowcool-simulated annealing refinement using the maximum-likelihood function and bulk-solvent modelling in CNS41 proceeded until the decreasing R/Rfree no longer converged. KDM4A.CP2(R6Kme3) complex was refined by alternative cycles of CNS41 and PHENIX39 using the ML function and bulk-solvent modelling. All residues were in acceptable regions of the Ramachandran plot as calculated by PROCHECK42 and/or MOLPROBITY43. Data collection and refinement statistics are in Supplementary Table 2.
HEK293T, HeLa and U2OS cells were obtained from the ATCC (LGC Standards, UK). FBS, L-glutamine and DMEM medium were purchased from Sigma; PBS and OptiMEM were from ThermoFisher Scientific. All cell lines were cultured in DMEM supplemented with 10% FBS, 2mM L-glutamine and 1 × Pen-Strep at 37°C and 5% CO2 unless otherwise stated.
HEK 293T cells (approximately 1 × 107 cells) in a 75cm2 flask were washed with PBS, harvested with trypsin, quenched with 10% FBS (Sigma)-containing DMEM (Sigma) and centrifuged at 200g for 3min at 4°C. The pelleted cells were re-suspended in 1ml buffer containing 20mM Tris-HCl (pH 7.6) and 140mM NaCl. The suspension was freeze–thawed five times using liquid nitrogen and a heating block set at 25°C, vortexing briefly between each thawing, followed by centrifugation (14,000g, 10min, 4°C).
CP2(T13Z) was added to 250μl lysate, yielding a final concentration of 200μM peptide. The contents were incubated at 37°C and 30μl aliquots withdrawn at various time points. Analyses were performed using a Thermo Exactive mass spectrometer equipped with Waters Acquity liquid chromatography system. Instrument control and data processing were performed using Thermo Xcalibur Software. The system was calibrated on the day of the analysis. The mass accuracy with external calibration (as used for these experiments) was better than 5ppm for 24 hours following calibration. Electrospray source conditions were adjusted to maximize sensitivity. Each sample (1μl) was injected into the system with a 2.1 × 50mm (3μm) Purospher STAR RP-18 endcapped column at a flow rate of 0.4mlmin−1. The oven was held at 40°C throughout the analysis. Solvent A was aqueous 0.1% formic acid and solvent B was acetonitrile. Peptides were eluted from the column using a gradient of 5–100% solvent B over 8min. All raw data files were analysed with MestReNova V.10.
For CETSA melting curve experiments, 80% confluent U2OS cells in a 75cm2 flask were transfected with a plasmid encoding full-length Flag-tagged KDM4A1–1,064 (ref. 15) or KDM4E1–506 (ref. 20) using Lipofectamine 2000 (Thermo Scientific) and incubated for 4h. The cells were washed with PBS buffer and media replaced with fresh DMEM with 10% FBS, supplemented with cyclic peptides in DMSO (0.1% DMSO final). Cell dosing was carried out for 20h at 37°C. Control cells were dosed only with DMSO. The cells were trypsinized and washed thrice with PBS, by centrifugation (300g, 3min, 4°C). Cells were re-suspended in PBS supplemented with Complete Protease inhibitor cocktail (Sigma) and divided into 8–10 aliquots (100μl) containing equal cell numbers in PCR tubes. For ITDRF CETSA experiments, approximately 6 × 106 cells were seeded evenly into nine 25cm2 flask and adhered before transfection. After 4h of incubation, the cells were washed with PBS buffer and the media replaced with fresh DMEM with 10% FBS supplemented with cyclic peptides at a final concentration between 0 and 4,000nM (0.1% DMSO final). Cells were trypsinized and washed thrice with PBS using centrifugation (300g, 3min, 4°C). Cells were re-suspended in PBS supplemented with protease inhibitor cocktail (Sigma), aliquots (100μl) containing equal cell numbers in PCR tubes were prepared. For CETSA melting curve experiments, each PCR tube was incubated for 3min in a heating block at different temperatures, then at room temperature for 3min. For ITDRF, tubes were incubated at approximately the Tm of the proteins of interest as determined by CETSA melting curve experiments, for example, 55°C (Tm for actin) for 3min, followed by room temperature for 3min. The tubes were centrifuged (300g, 3min, 4°C). The supernatant was removed and cells suspended in lysis buffer (100mM HEPES, 300mM NaCl, 2% NP-40 and 10mM EDTA, pH 7.4), supplemented with protease inhibitor cocktail. The tubes were incubated at 4°C for 1h, with vortexing every 20min. Samples were centrifuged for 30min at >16,000g at 4°C to pellet cell debris and precipitated proteins.
The samples (supernatant) were run on SDS–PAGE (NuPage Novex Bis-Tris 4–12% polyacrylamide gels with NuPAGE MES SDS running buffer (Thermo Scientific)) and transferred onto nitrocellulose membranes using either the iBlot blotting system. FLAG-KDM4A levels were determined by immunoblotting using a monoclonal mouse anti-FLAG primary antibody (F1804, Sigma), at a dilution of 1: 1,000 in PBS-T buffer containing Seablock Reagent Salmon Plasma (2% v/v, Calbiochem). An additional chicken anti-β-actin primary antibody (ab13822, Abcam) was also used. A goat anti-mouse fluorescent secondary antibody (800nm, 1:10,000, Thermo Scientific) and anti-chicken fluorescent secondary antibody (680nm, 1:10,000, Thermo Scientific) were used for imaging on a LI-COR Odyssey CLx imaging system. Western blotting band intensities were quantified using Image Studio Ver. 3.1 (LI-COR Biosciences) and plotted in Prism 6.07 software (GraphPad Software Inc.) to generate thermal stability curves.
The crystal structures KDM4A.Ni(II).CP2 and KDN4A.Ni(II).CP2(R6Kme3) have been deposited under PDB accession codes 5LY1 and 5LY2, respectively. All other data are available from the authors upon reasonable request.
How to cite this article: Kawamura, A. et al. Highly selective inhibition of histone demethylases by de novo macrocyclic peptides. Nat. Commun. 8, 14773 doi: 10.1038/ncomms14773 (2017).
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary figures, supplementary tables, supplementary methods and supplementary references.
We thank Dr Jumpei Morimoto for invaluable discussions and James Wickens for his assistance in MS analysis. This work was supported by the JSPS Grant-in-Aid for the Specially Promoted Research (21000005), the Core Research for Evolutional Science and Technology (CREST) Program of Japan Science and Technology Agency (JST), BHF Centre of Research Excellence Oxford (RE/08/004), Engineering and Physical Sciences Research Council (EP/L003376/1), Cancer Research UK (C8717/A18245), the People Programme (Marie Curie Actions) of the European Union's Seventh Framework Programme (FP7/2007–2013) under REA grant agreement number 298603 (Marie Curie IEF Fellowship to M.M.), the Wellcome Trust, and the European Research Council Starter Grant (679479). R.J.H. acknowledges a William R Miller Junior Research Fellowship from St Edmund Hall Oxford and A.K. is supported by the Royal Society Dorothy Hodgkin Research Fellowship. The SGC is a registered charity (number 1097737) that receives funds from AbbVie, Boehringer Ingelheim, the Canada Foundation for Innovation, the Canadian Institutes for Health Research, Genome Canada, GlaxoSmithKline, Janssen, Lilly Canada, the Novartis Research Foundation, the Ontario Ministry of Economic Development and Innovation, Pfizer, Takeda and the Wellcome Trust (092809/Z/10/Z).
The authors declare no competing financial interests.
Author contributions H.S., C.J.S. and A.K. conceived and designed the study experiments, and drafted the manuscript. T.K. and T. Kato performed the RaPID selections. T.K., Y.G., M.M. and R.J.H. synthesized the peptides. A.K., A.T., O.N.F.K. and B.B. conducted in vitro analysis and MS, and protein crystallography was carried out by R.C., O.N.F.K. and S.M. The cell-based assay was performed by A.K., C.Y., B.B., S.B.H., L.J.W. and T.P. All authors analysed the experimental data, discussed the results and were involved in preparation of the manuscript.