Oliver D.
Coleman
abc,
Jessica
Macdonald
b,
Ben
Thomson
b,
Jennifer A.
Ward
de,
Christopher J.
Stubbs
f,
Tom E.
McAllister
abc,
Shane
Clark
b,
Siddique
Amin
a,
Yimang
Cao
b,
Martine I.
Abboud
b,
Yijia
Zhang
b,
Hitesh
Sanganee
g,
Kilian V. M.
Huber
de,
Tim D. W.
Claridge
*b and
Akane
Kawamura
*abc
aSchool of Natural and Environmental Sciences - Chemistry, Newcastle University, Newcastle, NE1 7RU, UK. E-mail: akane.kawamura@ncl.ac.uk
bDepartment of Chemistry, University of Oxford, Chemistry Research Laboratory, Mansfield Road, Oxford, OX1 3TA, UK. E-mail: tim.claridge@chem.ox.ac.uk
cRadcliffe Department of Medicine, Wellcome Centre for Human Genetics, University of Oxford, Roosevelt Drive, Old Road Campus, Oxford, OX3 7BN, UK
dTarget Discovery Institute, Nuffield Department of Medicine, University of Oxford, Old Road Campus, Oxford OX3 7FZ, UK
eCentre for Medicines Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, OX3 7FZ, UK
fMechanistic and Structural Biology, Discovery Sciences, R&D, AstraZeneca, Cambridge, CB4 0WG, UK
gEmerging Innovations Unit, Discovery Sciences, R&D, AstraZeneca, Cambridge, UK
First published on 17th April 2023
Plant homeodomain fingers (PHD-fingers) are a family of reader domains that can recruit epigenetic proteins to specific histone modification sites. Many PHD-fingers recognise methylated lysines on histone tails and play crucial roles in transcriptional regulation, with their dysregulation linked to various human diseases. Despite their biological importance, chemical inhibitors for targeting PHD-fingers are very limited. Here we report a potent and selective de novo cyclic peptide inhibitor (OC9) targeting the Nε-trimethyllysine-binding PHD-fingers of the KDM7 histone demethylases, developed using mRNA display. OC9 disrupts PHD-finger interaction with histone H3K4me3 by engaging the Nε-methyllysine-binding aromatic cage through a valine, revealing a new non-lysine recognition motif for the PHD-fingers that does not require cation-π interaction. PHD-finger inhibition by OC9 impacted JmjC-domain mediated demethylase activity at H3K9me2, leading to inhibition of KDM7B (PHF8) but stimulation of KDM7A (KIAA1718), representing a new approach for selective allosteric modulation of demethylase activity. Chemoproteomic analysis showed selective engagement of OC9 with KDM7s in T cell lymphoblastic lymphoma SUP T1 cells. Our results highlight the utility of mRNA-display derived cyclic peptides for targeting challenging epigenetic reader proteins to probe their biology, and the broader potential of this approach for targeting protein–protein interactions.
PHD-fingers are prevalent among chromatin associated proteins, including >50 found in lysine methyltransferases (KMTs) and Jumonji-C (JmjC) histone demethylases (KDMs),13,14 which dynamically alter histone lysine methylation as part of epigenetic regulation in health and disease.2,3,15–18 The importance of PHD-fingers in such context is exemplified by KDM7B (PHF8, JHDM1F) which reads Nε-trimethyllysine at H3K4 (H3K4me3), generally considered an active histone mark at transcriptional start sites, to influence the linked JmjC-domain catalysed demethylation of H3K9me2.19 KDM7B is upregulated in melanoma to transcriptionally promote the TGFβ pathway through its demethylase activity and increase cell invasiveness.20 Inactivating mutations to either the catalytic JmjC-domain or PHD-finger greatly reduces metastasis in cell and animal studies, demonstrating that the PHD-finger function is crucial to its oncogenic activity.20 KDM7B is also involved in other cancer progressions, as well as neural and developmental processes,21–26 with other sub-family members KDM7A (KIAA1718) and KDM7C (PHF2) linked to similar conditions.27–30
Despite the biological importance, there are currently no potent and selective inhibitors for the PHD-fingers,14,31,32 which is in stark contrast to other families of chromatin reader domains (e.g. bromodomains, MBT, PWWP).33–35 The low druggability of the PHD-fingers with small molecules has been attributed to their shallow and open lysine binding pocket,36 and their extended histone binding interactions across the protein surface. H3K4me3 binding PHD-fingers are believed to interact with positively charged Nε-trimethyllysine primarily through cation–π interaction with their aromatic cage, typically consisted of two to four aromatic and hydrophobic amino acids, as well as via additional salt-bridging interactions.37 These features can allow some PHD-fingers to bind histones with remarkable methylation-state selectivity (KD can differ 10–1000 fold).
Some histone peptide-derived PHD-finger inhibitors containing K4me3 as a targeting motif have been reported which show promise,38,39 however, achieving selectivity over other epigenetic proteins remains a major challenge.
Here we report the development of a de novo natural product-like cyclic peptide (CP) ligand OC9 targeting the PHD-finger of histone demethylase KDM7 sub-family. We used cyclic peptide mRNA-display to screen a library with >1012 diversity of structurally constrained 3D scaffolds to identify OC9, which exhibits nanomolar affinity and high selectivity against isolated proteins and in cellular lysate context. OC9 disrupts the PHD-finger–H3K4me3 interaction through binding at the aromatic cage via an unprecedented valine motif. Moreover, PHD-finger inhibition gives allosteric control of JmjC demethylase activity at H3K9me2, inhibiting KDM7B but stimulating KDM7A, representing specificity of functional modulation within the KDM sub-family members.
Our work provides new molecular insight into PHD-finger recognition and their role in fine-tuned epigenetic control, and a new PHD-finger targeting scaffold for further inhibitor development.
Fig. 1 mRNA-display identifies high affinity cyclic peptides for KDM7B. (a) Crystal structure of KDM7B PHD finger (pink) and JmjC (grey) domains in complex with H31–24K4me3K9me2 (green) and 2-oxoglutarate co-factor mimic N-oxalyl glycine (NOG, yellow) (PDB ID 3KV4). The ‘YYMW’ aromatic cage of the PHD-finger binds K4me3 and the JmjC demethylates K9me2. (b) KDM7 sub-family domain architectures, and the main KDM7 constructs used in this study. (c) Structure of OC9, discovered by mRNA-display. (d) BLI traces for association and dissociation (separated by vertical dotted line) of OC9 to KDM7B domains. Each peptide concentration (2-fold serial dilution) response is a different colour, with curve fittings (1:1 binding model) overlaid in black. See Fig. S5† for further info. (e) ITC binding curves of OC9 with KDM7 domains, fitted with a 1:1 binding model. See Fig. S10† for further info. |
Six CPs (Fig. S1†) were initially characterised by bio-layer interferometry (BLI) to examine binding to KDM7B. Different CPs exhibited distinctive binding profiles to KDM7B protein constructs (dual domain, PHD-finger, JmjC-domain, Fig. 1b and S2†). CP hits were clustered into those that are: (i) specific to the dual domain only (OC3), (ii) specific to a single domain (OC4 – JmjC, OC9 – PHD) or (iii) non-binders. OC9 (Fig. 1c) exhibited near-equivalent nanomolar affinities and dissociative half-lives in the order of minutes for both the dual domain (KDBLI = 6 nM, t1/2 = 2.8 min) and PHD-finger (KDBLI = 4 nM, t1/2 = 8.1 min) of KDM7B, with no affinity for the isolated JmjC, indicating that it likely binds the PHD-finger alone (Fig. 1d). Thus, OC9 was prioritised for further development. Isothermal titration calorimetry (ITC) confirmed OC9 affinity for KDM7B with KD values of 11 nM (dual domain) and 12 nM (PHD-finger) (Fig. 1e). OC9 also showed similar affinity and binding profiles against the related KDM7A and KDM7C dual domains and PHD-fingers, consistent with high KDM7 PHD-finger sequence homology (>80% identity to KDM7B, Fig. S3†). The affinities of OC9 against the KDM7 PHD-fingers were at least 10-fold greater than the natural histone peptide ligand (H31–21K4me3, ITC KD = 115 nM to KDM7B PHD-finger). For further BLI and ITC data, see Fig. S4–S10.†
To further evaluate the interactions of OC9 with KDM7B we used hydrogen–deuterium exchange mass spectrometry (HDX-MS), which informs on protein dynamics by following deuterium incorporation into the protein amide backbone (Fig. 2). Briefly, protein was incubated with or without ligand (CP or H31–21K4me3) in deuterated buffer over time and quenched by acidification and cooling. Proteolysis and LC-MS analysis of the nascent peptides and comparison with a DMSO-control identified protein regions affected by ligand binding. For the dual domain, OC9 and H31–21K4me3 only reduced deuterium uptake in peptide fragments within the PHD-finger (Fig. 2, S11–S14, HDX-MS ESI†). The lack of observable effect on the JmjC by H31–21K4me3 could be due to a very weak or non-substantial interaction of its unmethylated H3K9 with the JmjC. However, the dual domain specific CP OC3 (ITC KD = 163 nM, Fig. S9†) affected deuterium uptake in both domains, confirming that dual domain binding can be mapped using HDX-MS.
Fig. 2 OC9 specifically targets the PHD-finger. (a) Regions of KDM7B dual domain and PHD-finger proteins affected in HDX-MS after 5 minutes of incubation with OC9, OC3 or H31–21K4me3, assessed in triplicate relative to DMSO. The residue level view was generated by averaging mass changes across overlapping peptide fragments. Supporting data available in ESI.† K4me3 binding aromatic cage residues (‘YYMW’ Fig. 1a) are indicated by red asterix. (b) Combined HDX-MS data for all time points (1, 5, 15, 30, 60 minutes) for dual domain and PHD-finger protein mapped to a structure (PDB ID 3KV4) as protected (blue), not significantly affected (white) and without sufficient coverage (black). |
Further HDX-MS on KDM7B PHD-finger alone showed that H31–21K4me3 reduced deuterium uptake across a substantial proportion of the domain over a 60 minute time-course, whilst OC9 affected an even greater proportion including the regions that contain the function-critical aromatic cage residues ‘YYMW’ (Fig. 2a and S11A†). Thus, OC9 exclusively binds at the PHD-finger in agreement with BLI and ITC.
Other positions (L3, Y7) were also important as they showed 10 to 100-fold reductions in potency when replaced with alanine.45–47 Residues T4 and H5 tolerated L-alanine and D-alanine substitutions, whilst P8 only tolerated L-alanine. All tyrosines could be exchanged for phenylalanine, individually and in combination, with minimal impact. Bioinformatic analysis of NGS data from mRNA-display supported these empirical SAR findings, with a Logo plot for the CP family cluster related to OC9 (902 distinct sequences in total) highlighting the enrichment of sequences with the ‘SVXR’ motif (where X = F/Y/W), aromatic residues favoured at positions 2, 7, 11 and greater variability found at positions 4 and 5 (Fig. S20†). Notably, V10 was the most highly conserved residue. Characterisation of an OC9 variant with the second-most abundant residue (leucine) at position 10 (OC9 V10L) resulted in >80-fold reduction in KD and >75-fold reduction in t1/2 in BLI, and 10-fold reduction in IC50AS compared to OC9 (Fig. S16, S18B and S19†), in line with its significantly lower enrichment relative to OC9 by NGS (0.003% vs. 3.2% abundance). Thus, future SAR work may be accelerated by further exploiting such NGS-based peptide family analysis.
Fig. 4 Protein NMR reveals valine contacting of the aromatic cage. (a) 1H–15N HSQC of KDM7B PHD-finger in the apo (black) or bound forms with OC9 (purple, 1:1) or H31–21K4me3 (green, 3:1). (b) 1H and 15N chemical shift changes for apo vs. + OC9 were combined and weighted to give chemical shift perturbation (CSP) (Δδ) values. Standard deviation (σ) of Δδ are indicated (>1 pink, >2 purple). (c) TALOS+ predicts similar secondary structure elements in apo and OC9-bound solution protein complexes (light grey) to those found in PDB ID 3KV4 (dark grey). (d) Mapping of CSP σ from (b) onto PDB ID 3KV4 shows OC9 binding surfaces primarily along one face, which presents the cage (annotated residues). (e) Shielded 1H-NMR signals centred around −0.5 ppm belonging to the OC9 V10 γ-methyl groups are observed in some OC9–PHD complexes, but not in apo protein, or complexes with OC9 V10A and OC9 V10DV. The V10OC9 γ-methyl signals have NOE interactions (black semi-circles) with the W29PHD indole, which itself is shielded upon binding of OC9 (indicated by ΔδH values). |
OC9 affected the PHD-finger differently to H31–21K4me3 based on spectra overlay, and was in slow exchange consistent with a slow off-rate on the NMR timescale and high affinity. Among the residues affected by OC9 binding, all four aromatic cage residues (Y7PHD, Y14PHD, M20PHD, W29PHD) were shifted, with changes in the W29PHD indole NH also observed (ΔδH = +0.86 ppm). The Y14PHD, M20PHD and W29PHD cage residues also exhibited changes in their 1H amide temperature coefficients upon OC9 complexation (Fig. S22†) to further suggest the aromatic cage is directly affected, despite absence of any lysine residue in OC9. Other affected residues included F19PHD and Y48PHD which, alongside I21PHD, E22PHD and I45PHD, have been found to interact with H3 residues A1, R2 and T3 in crystallographic studies19 and so their disruption could also contribute to observed OC9 histone displacement.
The PHD-finger remained well folded despite the extensive number of affected residues, as overall 1H–15N HSQC signals were well dispersed with similar secondary structure elements in solution predicted by TALOS+49 for apo and bound protein using chemical shift data (Fig. 4c and S23†). These elements were in broad agreement with previous X-ray crystallography (PDB ID 3KV4),19 likely due to the stabilising influence of two zinc coordination centres. Mapping of all affected residues showed they cumulatively formed a probable binding interface involving the β1, β2, β3, and β4 strands and α2 region, with negligible effect on α1 and α3 on the reverse face of the PHD-finger (Fig. 4d), in reasonable agreement with HDX-MS affected regions.
1H-NMR data also revealed a pair of distinct signals in the OC9 complex around −0.5 ppm (Fig. 4e, S24A and S25†). Titration of OC9 L3A and F11A saw these signals progressively increase, with their overall 1H–15N HSQC shift patterns indicating binding surfaces similar to OC9 (Fig. S26–S31†). However, the −0.5 ppm signals were not observed with OC9 V10A, lacking γ-methyl groups (Fig. 4e). Subsequent TOCSY and NOESY experiments did attribute the −0.5 ppm resonances to the two V10OC9 side chain γ-methyl groups, and also, via NOEs, found them to be in close proximity to the W29PHD indole ring (Fig. 4e and S32†). A 2D-TOCSY investigation showed cage residue W29PHD exhibited shielded indole proton 1H chemical shifts (5–6 ppm) in the OC9 complex relative to the apo PHD-finger (Fig. 4e and S33†), which could result from shielding induced by V10OC9 or adjacent residues.
Notably, inversion of the valine stereocentre in OC9 V10DV severely weakened affinity (NMR KD = 503 ± 1 μM, pKD < 3.3, supported by BLI KD > 10 μM, Fig. S34†) and did not present the shielded methyl signals. Protein 1H-NMR also showed clear changes in the aromatic region (∼6 ppm) upon OC9 complexation with new signals appearing, but weak-binding OC9 V10DV instead elicited broadening of peaks (Fig. S24A†). Together, this suggests the presence and spatial orientation of the OC9 valine, imposed by conformational constraint of cyclisation, is critical to OC9 potency by binding at the aromatic cage to disrupt the methyllysine host site. This contrasts with the ability of H3K4me3 reader domains (PHD-fingers and tandem-Tudor) to accommodate D-K4me3 H3 derivatives (albeit with weaker affinity than L-K4me3).50
Interestingly, when V10 was replaced with Kme3 in OC9 V10Kme3, PHD-finger binding was moderately retained (KDBLI = 533 nM, IC50AS = 316 nM) whereas no binding was observed when replaced with K (OC9 V10K) (KDBLI, IC50AS > 10 μM) (Fig. S16, S18B and S19†). The higher affinity for the methylated over unmethylated lysine variant of OC9 reflects the pattern of histone PTM recognition for KDM7 PHD-fingers (Fig. S8†) and provides further evidence for the engagement of OC9 position 10 with the aromatic cage. The binding of Kme3 in KDM7 and other PHD-fingers is proposed to be driven by a combination of favourable cation–π interaction between the positively charged quaternary ammonium group and the electron-rich aromatic residues, and the displacement of high-energy water molecules within the cage.52 However, binding to an aromatic cage without a cation–π interaction, as proposed for OC9 valine insertion, is precedented by the all-carbon, neutral, non-natural tert-butyl norleucine (t-BuNle) substituted H3K4 derived peptides which retain affinity against Kme3-binding PHD-fingers.52,53 In such cases, dispersion forces and hydrophobic effects are proposed to drive the aromatic cage engagement, as exemplified for CBX5 chromodomain, where affinity increases with increasing size of alkyl substituent on lysine ε-N.54
Further characterisation of OC9 binding was carried out from partial 1H chemical shift assignments derived for bound-OC9 obtained from isotope filtered 2D experiments (Fig. S24B†). While multiple interconverting conformers were observed for free OC9 (Fig. S24A†), only a single conformer population was observed when OC9 was bound to PHD-finger (Fig. S24B†). The conformer switching in the free form is likely driven by isomerism of the two prolines present in OC9. Substitution of OC9 with 4-fluoroproline at P6 gave higher affinity (>10-fold) for the trans-fluoro exo-puckered ring favouring the trans-amide bond45–47 over the cis-fluoro (Fig. S18 and S19†), suggesting trans-to be the more active PHD-finger binding conformer. In addition to the very substantially shielded 1H resonances of the OC9 V10 sidechain, the bound-OC9 NMR data also showed that the ortho aromatic protons of F11OC9 were similarly shielded by ∼1.3 ppm (relative to random coil values51), appearing at 6.0 ppm. Although no clear NOEs could be observed between the shielded F11OC9 ortho protons and W29PHD, NOEs could be identified to the spatially nearby I21PHD methyl groups. The beta protons of R12OC9 also experienced appreciable shielding (∼1 ppm versus random coil). The H3R2 residue is known to anchor to a cage-adjacent polar/acidic patch of the KDM7 PHD-finger through a guanidinium-carboxylate salt bridge to E22PHD to assist binding,45–47 with similar concomitant recognition of unmodified H3R2 and H3K4me3 also characteristic of many other PHD-fingers including BPTF, ING2 and KDM5A PHD3.6 Thus, since variant OC9 R12A also exhibited a considerably faster off-rate (∼100×) by BLI, R12OC9 could plausibly mimic the anchoring role of the histone H3R2 residue to control OC9 residence time. Collectively, these data suggest the OC9 V10-F11-R12 residues to be interacting significantly with the β2–β3–α2 face of the PHD finger, although there is insufficient data to definitively position OC9.
We tested the effect of OC9 on KDM7B demethylation of H31–21K4me3K9me2 and the less-efficient H31–21K9me2 substrate by MALDI-TOF MS based assay.14,19 OC9 inhibited KDM7B demethylation of both substrates with similar potencies, with IC50 values of 1.19 μM and 0.51 μM respectively (Fig. 5a, S39 and S40†). Given the lack of affinity of unmethylated H3K4 for the PHD-finger (Fig. S8†),19 the inhibitory effect of OC9 on H31–21K9me2 demethylation was unexpected, but highlights the likely important allosteric role of the PHD-finger beyond PTM recognition. The potencies of OC9 were comparable to that of the KDM2/7 JmjC-active site small molecule inhibitor daminozide (DAM),58 and an equimolar mix of OC9 and DAM showed an additive inhibitory effect. Inhibition by OC9 was further confirmed by an orthogonal AS activity assay58 (Fig. S41†). OC9 also reduced the demethylation activity of KDM7B on H3K9me2 in calf histone extracts, which possess naturally occurring combinatorial PTMs (Fig. 5a and S42†). This suggests KDM7B demethylase activity against H3K9me2 is largely dependent on PHD-finger driven positioning of the N-terminal histone H3 tail regardless of H3K4 methylation state. As KDM7B aromatic cage mutation results in the same phenotypic outcome as JmjC catalytic site inactivation (e.g. KDM7B demethylase-dependent melanoma metastasis20), our results emphasise the crucial role of PHD-fingers in regulating KDM7 demethylase activity and the need to explore PHD-finger viability as therapeutic targets.
Fig. 5 KDM7 PHD-finger inhibition by OC9 modulates dual domain JmjC activity. (a) (i–iii) OC9 inhibits demethylation of H3K9me2, with or without K4me3, by KDM7B (0.4 μM) on H3 peptide (2.5 μM). (iv) OC9 (20 μM) also inhibits KDM7B (1 μM) demethylation of calf histones (50 μg). (b) (i–iii) OC9 enhances demethylation of H3K4me3K9me2 by KDM7A (0.55 μM) on H3 peptide (2.5 μM), and (iv) enhances KDM7A (0.5 μM) demethylation of calf histone extracts (50 μg). Peptide demethylation assays analysed by MALDI-TOF-MS, with potency comparisons of N = 2 experiments in technical triplicate. Inset MALDI-TOF-MS overlays are from respective highest OC9 concentration points compared to DMSO alone. KDM7 PHD-finger only proteins had no demethylase activity (Fig. S40†). Calf histone H3K9me2 was quantitated relative to H4 band intensities at each time point, normalised to DMSO at t = 0. Full Western blots available in Fig. S42.† Relative rates (Krel min−1) from a one-phase decay model, average from N = 3 independent experiments. |
In contrast to KDM7B, KDM7A lacks activity against H3K4me3K9me2, but demethylates H31–21K9me2.57 OC9, however, restored the activity of KDM7A against H31–21K4me3K9me2 from <5% to a level (∼30%) that approached demethylation on H31–21K9me2 with DMSO alone (∼40%) (Fig. 5b and S39†). Daminozide remained inhibitory against KDM7A and prevented restoration of activity in the equimolar mix with OC9. Activity enhancement by OC9 translated to an increased rate of demethylation against calf histone extracts (Fig. 5b and S42†). These combined results indicated OC9 blocks PHD-finger binding to K4me3 to allow demethylation of K9me2 by KDM7A, and further emphasises that OC9 binding is confined to the PHD-finger. The effect of PHD-finger inhibition by OC9 for KDM7C is unclear, due to its lack of demethylase activity in vitro.59 Due to the negligible permeability of OC9 in Caco-2 cell assay, the activity was not tested in cells. Further optimisations are likely required to improve the cellular permeability of OC9.
Screening against other isolated KDM sub-families (KDM2/3/4/5/6) with histone peptide substrates saw no significant inhibition by OC9. The panel included full-length recombinant KDM2A and KDM2B which have the closest JmjC-domains to KDM7s, but with unrelated PHD-fingers (Fig. S43 and S44†). Hence, OC9 not only represents the first selective and potent inhibitor for the PHD-fingers, but also of the demethylase activity for the KDM7 sub-family.60
While the use of de novo CPs to selectively inhibit epigenetic eraser and reader domains have been reported (e.g. KDM4 JmjC-domain,61 SIRT2,62 TET63 and bromodomain (BRD)64 sub-families), the PHD-finger targeting CPs provide unique context-dependent allosteric control via modulation of the tandem JmjC-domain. Thus, OC9 also exemplifies a PHD-finger targeting strategy for KDM inhibitor development, which can overcome selectivity challenges encountered by small molecules binding at the highly conserved active site of JmjC-domains.58,65,66
A further chemoproteomic profiling analysis of endogenous protein capture by OC9-Bt from KDM7B-expressing T-cell lymphoblastic lymphoma SUP-T1 nuclear lysates67 identified 454 proteins. Only the relative abundance of KDM7B and KDM7C were significantly affected in a dose-dependent manner by competitive titration of unlabelled OC9 during pull-down (Fig. 6b, Proteomics ESI, and Western blot Fig. S47†). This proteomics data suggests remarkable selectivity and affinity of OC9 for KDM7 PHD-fingers, even in complex nuclear lysates abundant in chromatin associated proteins including H3K4me3 binders.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d2sc05944d |
This journal is © The Royal Society of Chemistry 2023 |