Gloria
Ortiz
a,
James E.
Longbotham
a,
Sophia L.
Qin
a,
Meng Yao
Zhang
a,
Gregory M.
Lee
bc,
R. Jeffrey
Neitz
bc,
Mark J. S.
Kelly
c,
Michelle R.
Arkin
bc and
Danica Galonić
Fujimori
*ac
aDepartment of Cellular and Molecular Pharmacology, University of California San Francisco, San Francisco, CA 94158, USA. E-mail: danica.fujimori@ucsf.edu
bSmall Molecule Discovery Center (SMDC), University of California San Francisco, San Francisco, CA 94158, USA
cDepartment of Pharmaceutical Chemistry, University of California San Francisco, San Francisco, CA 94158, USA
First published on 18th December 2023
PHD fingers are a type of chromatin reader that primarily recognize chromatin as a function of lysine methylation state. Dysregulated PHD fingers are implicated in various human diseases, including acute myeloid leukemia. Targeting PHD fingers with small molecules is considered challenging as their histone tail binding pockets are often shallow and surface-exposed. The KDM5A PHD1 finger regulates the catalytic activity of KDM5A, an epigenetic enzyme often misregulated in cancers. To identify ligands that disrupt the PHD1-histone peptide interaction, we conducted a high-throughput screen and validated hits by orthogonal methods. We further elucidated structure–activity relationships in two classes of compounds to identify features important for binding. Our investigation offers a starting point for further optimization of small molecule PHD1 ligands.
PHD fingers bind N-terminal tails of histone proteins to stabilize interactions of chromatin-acting proteins with chromatin and can also serve as allosteric modulators of the associated enzymes. Aberrant PHD–histone interactions, and mutations in PHD fingers, are implicated in human pathologies, including cancer and neurological and immunological disorders.6 PHD fingers in KDM5A,7 PHF23,8,9 and BPTF10,11 are fusion partners of nucleoporin-98 (NUP98) in acute myeloid leukemia. Although PHD fingers share low primary sequence similarity, they adopt a conserved globular structure that binds its histone target within a shallow binding groove. The surface-exposed binding pockets and shallow binding grooves recognize the histone ligands, which can adopt a β-sheet confirmation or helical structure to accommodate extended interactions.7,12–14 The nature of the binding pocket makes PHD fingers challenging targets, and difficult to engage with small molecules.1,3,15 This is particularly the case for PHD fingers that engage H3K4me0, which primarily rely on polar interactions for ligand binding.
There are several literature reports describing the identification of small-molecule ligands for PHD fingers. The first report utilized a HaloTag pull-down screen to identify amiodarone, a KDM5A-PHD3 inhibitor, but subsequent structure-to-activity relationship (SAR) studies revealed this class of compounds displayed weak inhibition and poor selectivity.16,17 A 15N HSQC-based fragment screen focusing on identifying ligands for the PHD finger of Pygo identified a benzimidazole derivative that binds to a cleft adjacent to the histone binding pocket.18 Another report utilizing an NMR-based fragment screen identified millimolar-binding fragments for the PHD fingers of BAZ2A/BAZ2B that inhibited the PHD–histone interaction.19 More recently, ligands that displayed low micromolar potencies for the PHD finger in UHRF1 were discovered through an AlphaScreen®-based HTS.20 While these examples are encouraging, the reports of small molecule ligands for PHD fingers remain rare and with low affinities for the target of interest.
The paucity of small molecule ligands and their low affinity for PHD fingers prompted us to investigate large-scale high-throughput screens as a strategy to discover small molecule modulators for this class of chromatin readers. We selected the PHD1 finger of histone demethylase KDM5A as the target for our investigations due to its importance in the regulation of the catalytic activity.21,22 KDM5A is a multi-domain enzyme with an iron-containing active site consisting of jumonji-C and jumonji-N (JmjC/N) domains, an AT-rich interactive domain (ARID) that binds DNA, a zinc-finger (Znf) domain, and three PHD fingers. Overexpression of KDM5A is associated with breast cancer proliferation and drug resistance,23,24 lung tumorigenesis and drug resistance,25,26 and gastric cancer development and progression.27,28 Active site inhibitors of KDM5A, such as selective, nanomolar inhibitors CPI-455,29,30 KDOAM-25,31 and KDM-C49 and its ester form C-7032 have found success in cancer cell models. In particular, CPI-455 was shown to reduce drug-tolerant subpopulation of cancer cells in multiple cancer cell lines, and KDOAM-25 impaired the proliferation of human myeloma cells.30,31 However, the competition with α-KG hampers the cellular efficacy of active site inhibitors. Our recent investigations into the regulation of catalysis in KDM5A identified PHD1 as an allosteric regulatory site in this enzyme. Specifically, binding of the product of demethylation to the PHD1 finger stimulates catalysis, enabling feed-forward regulation.21,22 The NMR structure of apo and bound PHD1 construct revealed binding to the N-terminus of the H3 peptide in a helical conformation, with extended hydrogen bonding and electrostatic interactions (Fig. 1a).12 Here, we describe the identification of PHD1 ligands through a fluorescence polarization-based high-throughput screen and several rounds of structural derivatization. 1H–15N HSQC NMR was used to validate binding and identify regions of PHD1 that interact with the ligands. The nascent SAR of the benzofuran series of ligands paves the way for further optimization.
Fig. 1 Targeting the PHD1 finger of KDM5A. (a) NMR solution structure of the PHD1 finger of KDM5A (blue) bound to histone H3 10-mer peptide (orange) (PDB: 7klr). (b) Fluorescence polarization-based high-throughput screening plan to identify ligands for PHD1. (c) Relative FP of C-terminal 5-FAM H3 10-mer peptide binding to GST-PHD1. (d) Competitive displacement of 5-FAM H3 10-mer (10 nM) binding to GST-PHD1 (2 μM) by the unlabeled H3 10-mer peptide. Error bars are ± SD for triplicate measurements. |
The FP assay performance in a high-throughput, automated format was first assessed by investigating the Z prime (Z′) factor, a measure reflecting the dynamic range of the assay signal and the data variation associated with measurements.35 Varying the concentration of GST-PHD1 (0.5, 1, 1.5, and 2 μM) in the presence or absence of 0.1% BGG revealed that the FP signal arising from the negative and positive control buffers was not affected by the presence of BGG and that 1 μM GST-PHD1 was sufficient to maintain a good dynamic range, as indicated by Z′ values greater than 0.5 in buffer without BGG (Z′ = 0.56) and with 0.1% BGG (Z′ = 0.59) (Fig. S1, ESI†). The pilot screen was then carried out using a fragment library of 3,095 compounds to ensure optimized assay conditions. The fragments were screened at 50 μM due to low molecular weight (median MW = 207 Da) and typical low affinities. Parallel screens were run with or without bovine gamma globulin (BGG) in the assay buffer to determine whether the hit rate would be affected by the presence of a carrier protein. The pilot screen was robust, producing an average Z′ = 0.72 ± 0.03 without BGG and an average Z′ = 0.61 ± 0.03 with BGG. Compounds that displaced the 5-FAM-labeled H3 peptide and displayed inhibition greater than or equal to 3 standard deviations above the mean were considered as hits. The pilot fragment screen without BGG produced 61 hits (1.98% hit rate). For the pilot screen using BGG, 15 hits were produced (0.5% hit rate), indicating that BGG aids in the removal of false-positive hits (Fig. S2a and b, ESI†). Pan-assay interference (PAINS), compounds that often give false-positive results due to reactivity with biological nucleophiles, redox-reactivity, photo-reactivity, metal chelators, or interfering photochromic properties,36 were excluded. From the 15 hits obtained through the pilot screen using BGG, one was a PAINS compound and was excluded (Fig. S2c, ESI†). The commercially available compounds (13) were repurchased, assessed for purity using LCMS, and evaluated in dose–response using the FP assay. Out of the 12 soluble compounds tested, only 7 fragments showed dose-dependent inhibition while 5 compounds displayed no measurable activity (Fig. S2d–j, ESI†).
The binding of the ChemBridge screen hits to PHD1 was further evaluated using 2D 1H–15N HSQC NMR spectroscopy. This direct binding method allowed us to monitor binding by observing any ligand-induced changes in the previously assigned 1H–15N PHD1 amino acid resonances.12 The HSQC spectra of 50 μM 15N-labeled PHD1 with DMSO were overlaid with the spectra of the 15N-labeled PHD1 incubated with 1 mM compound. Compounds 714997 and 727150 did not induce any significant chemical shift perturbations (CSP) on any PHD1 resonances (data not shown). Compound 715958 was insoluble in the NMR buffer at this concentration and was not evaluated.
The lack of validated hits from the ChemBridge Diversity collection screen motivated us to expand our efforts to libraries better suited for the identification of protein–protein interaction inhibitors. Specifically, the screen was expanded to include 7027 compounds from the National Cancer Institute (NCI) collection library which encompassed protein–protein interaction inhibitors, macrocycles, and natural products (Fig. S4, ESI†). The libraries were screened at 50 μM, resulting in the identification of 38 hits with percent inhibition 3 standard deviations above the mean (0.54% hit rate; Z′ = 0.66 ± 0.04). Among these, the eighteen exact compounds and four structurally similar compounds that were commercially available were repurchased and evaluated by dose–response format using the FP assay. Six compounds displayed measurable dose-dependent inhibition of the PHD1-H3K4me0 interaction (Fig. 2).
The six NCI screen hits confirmed by dose response in the FP assay were further evaluated by 2D 1H–15N HSQC NMR at 1 mM to assess their direct binding to PHD1 (Fig. S5, ESI†). Compound 990670 was insoluble in the NMR buffer and was excluded from further consideration. Out of the five compounds evaluated, 991830, 990856, and 988385 caused very minor CSP on PHD1 (average shifts ≤ 0.002 ppm; Fig. S5a–c, ESI†). Conversely, 516969, and 990099 displayed greater average perturbations (0.0025–0.0033 ppm), while some residues were perturbed greater than 0.008 ppm (Fig. S5d and e, ESI†). The observed perturbations are low, but in the range of those induced by competitive binding fragment ligands evaluated at 5 mM with the PHD finger of BAZ2B, for which the perturbations observed were up to 0.035 ppm.19 Residue Cys322, one of the Zn(II) ligands in the first zinc finger in PHD1, is the most affected by the presence of oxazepine 516969 (Δδ of 0.009 ppm; Fig. S5d, ESI†). The most prominent perturbations caused by benzofuran-containing compound 990099 were those to residues N303, N304, D316, C322, and I324 (Fig. S5e, ESI†). Interestingly, N304 was also perturbed by aryl sulfonamide 988385 (Δδ > 0.01 ppm; Fig. S5c, ESI†). This residue is located near the first Zinc finger and is solvent-exposed.
The two molecules that caused the most prominent chemical shift perturbations, oxazepine 516969 and benzofuran 990099, were further evaluated. Given the complexity of the scaffold, we identified six commercially available close analogs of 516969 to evaluate whether potency could be improved. Derivatives of 516969 were evaluated by FP and their concentration-response curves are shown in Fig. S6 (ESI†). Compound 5624800 is structurally similar to NCI screen hit 516969, but contains the pyridine nitrogen atom at an ortho position. It displayed an approximate IC50 of 0.65 mM (0.43 to 1.1 mM 95% CI), showing a modest improvement compared to 516969 (Fig. S6b, ESI†). The phenyl analog 5487083 and substituted pyrazole analog 6378665 led to a decrease in potency (Fig. S6a and c, ESI†). Substituting the methyl group of 516969 with ethyl in 5620291 also resulted in a drop in potency (Fig. S6d, ESI†).
The structural similarity between benzofuran 990099 (Fig. 2e) and fragment hit 302387 (Fig. S2g, ESI†) motivated us to further investigate this series. Removal of the carboxylate from 990099 (8018–7366), drastically reduced potency, suggesting the carboxylic acid is important for PHD1 binding (Table 1 and Fig. S7a, b, ESI†). To systematically investigate how the R3 group on the sulfonamide substituent impacts binding, we prepared a series of derivatives where the C3 and C4 positions of the aryl ring carried various substitutions. To facilitate the rapid preparation of derivatives, the benzyl ester was replaced by a methyl ketone functionality. The synthesis of the benzofuran derivatives is outlined in Scheme S1 (ESI†) and follows a published route to obtain amino-substituted benzofuran 1.37,38 The route to access the benzofuran derivatives begins with a nucleophilic substitution reaction with commercially available sulfonyl chlorides and compound 1 to give substituted benzofurans B1–B7 (Scheme S1, ESI†).
Entry | Compound | R1 | R2 | R3 | %I at 1 mMa |
---|---|---|---|---|---|
a Data were collected in quadruplicate and display the % inhibition at 1 mM (n = 1 for R1 = OBn; n = 2 for R1 = Me). | |||||
1 | 990099 | OBn | H | 45% | |
2 | 8018-7366 | OBn | H | 6% | |
3 | B1 | Me | H | 38% | |
4 | 7917645 | Me | Ac | 30% | |
5 | B2 | Me | H | 45% | |
6 | B3 | Me | H | 12% | |
7 | B4 | Me | H | 9% | |
8 | B5 | Me | H | 11% | |
9 | B6 | Me | H | 4% | |
10 | B7 | Me | H | 13% |
The methyl ketone analog showed comparable binding affinity to the benzyl ester (Table 1, entry 1 vs. entry 3; Fig. S7a and c, ESI†). Capping the sulfonamide nitrogen with an acyl group in 7917645 (Table 1, entry 4; Fig. S7d, ESI†) slightly lowered binding affinity when compared to B1. Relative to p-carboxylate B1, placing the carboxylate substituent into the meta position on the aryl ring of B2 improved inhibition (Fig. S7c and e, ESI†). B2 displayed 45% inhibition at 1 mM, while B1 displayed 38% inhibition at 1 mM (Table 1, entries 3 vs. 5). Compound B2 displayed increased inhibition with an approximate IC50 of 1.3 mM (1.1 to 1.5 mM 95% CI) relative to 990099, which does not reach 50% inhibition when tested up to 1.5 mM (Fig. S7a and e, ESI†). Replacing the acid with an amide (B3 and B4) or an ester (B5 and B6) also decreases the percent inhibition observed (Table 1, entries 6–9; Fig. S7f–I, ESI†). Finally, replacing the carboxylic acid in B2 with a chloro group in B7 also weakens inhibition (Table 1, entry 10; Fig. S7j, ESI†).
Direct binding of the identified hits to PHD1 was further probed by NMR. The binding of the best-performing synthesized compound, B2, was assessed using HSQC NMR with 15N PHD1. Compound B2 was tested at 0.5, 1, 2, and 4 mM with 50 μM 15N PHD1. As shown in Fig. 3 and Fig. S8 (ESI†), many PHD1 residue signals broadened with increasing concentration of B2, indicating intermediate exchange. While the chemical shift perturbations are small, the residues perturbed greater than one standard deviation above the mean, shown in Fig. 3, include G313, C314, D316, F321, C322, I324, and K339. Some of the residues perturbed by B2 are near the second Zn finger, including G313, C314, D316, and K339 (Fig. 1a and 3). These findings suggest that compound B2 may be interacting with residues in this region or interacting with Zn directly.
Fig. 3 1H–15N PHD1 chemical shift perturbations induced by compound B2. (a) CSPs of PHD1 residues with 4 mM B2 (red) relative to 50 μM PHD1 with DMSO (blue). (b) CSPs induced by compound B2 that are greater than one standard deviation above the mean, indicated in panel c (PDB: 7klo). Residues colored in red display Δδ + 2σ, while residues in salmon display Δδ + σ. (c) CSP map of PHD1 residues affected by 4 mM compound B2 relative to PHD1 with DMSO. CSP bars highlighted by a magenta rectangle indicate residues that display broadening. |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3cb00214d |
This journal is © The Royal Society of Chemistry 2024 |