Oakley B.
Cox
abc,
Tobias
Krojer
a,
Patrick
Collins
c,
Octovia
Monteiro
ab,
Romain
Talon
a,
Anthony
Bradley
a,
Oleg
Fedorov
ab,
Jahangir
Amin
d,
Brian D.
Marsden
ae,
John
Spencer
d,
Frank
von Delft
*acf and
Paul E.
Brennan
*ab
aStructural Genomics Consortium (SGC), University of Oxford, Oxford OX3 7DQ, UK
bTarget Discovery Institute (TDI), Nuffield Department of Medicine, University of Oxford, Oxford OX3 7FZ, UK. E-mail: paul.brennan@sgc.ox.ac.uk
cDiamond Light Source (DLS), Harwell Science and Innovation Campus, Didcot, OX11 0DE, UK. E-mail: frank.von-delft@diamond.ac.uk
dDepartment of Chemistry, School of Life Sciences, University of Sussex, Brighton, BN1 9QJ, UK
eKennedy Institute of Rheumatology, Nuffield Department of Orthopaedics, Rheumatology and Musculoskeletal Sciences, University of Oxford, Roosevelt Drive, Headington, Oxford OX3 7FY, UK
fDepartment of Biochemistry, University of Johannesburg, Aukland Park 2006, South Africa
First published on 22nd December 2015
Research into the chemical biology of bromodomains has been driven by the development of acetyl-lysine mimetics. The ligands are typically anchored by binding to a highly conserved asparagine residue. Atypical bromodomains, for which the asparagine is mutated, have thus far proven elusive targets, including PHIP(2) whose parent protein, PHIP, has been linked to disease progression in diabetes and cancers. The PHIP(2) binding site contains a threonine in place of asparagine, and solution screening have yielded no convincing hits. We have overcome this hurdle by combining the sensitivity of X-ray crystallography, used as the primary fragment screen, with a strategy for rapid follow-up synthesis using a chemically-poised fragment library, which allows hits to be readily modified by parallel chemistry both peripherally and in the core. Our approach yielded the first reported hit compounds of PHIP(2) with measurable IC50 values by an AlphaScreen competition assay. The follow-up libraries of four poised fragment hits improved potency into the sub-mM range while showing good ligand efficiency and detailed structural data.
Fig. 1 Bromodomain family tree with sub-families (I–VIII) as defined by Filippakopoulos et al.1 Brds for which SGC Chemical Probes have been released are highlighted and the structures depicted: SGC-CBP30,11 I-CBP112, (+)-JQ1,3 PFI-1,7 PFI-3, BAZ2-ICR,15 GSK cmpd 3, OF-1, LP99,12 I-BRD9 and the pan-Brd inhibitor bromosporine.16 Atypical bromodomains are highlighted for which the conserved asparagine residue in the KAc binding site is replaced by tyrosine (blue), threonine (green) or aspartic acid (purple). |
Most Brd inhibitors are anchored by a highly conserved asparagine (Asn) residue and a network of four water molecules in the KAc binding pocket. Together the Asn and first water in the chain of four donate two H-bonds to a pair of H-bond acceptors on the inhibitor. The importance of the conserved Asn residue in inhibitor binding reflects its critical role in KAc recognition of the natural ligand (ESI Fig. S1A†).1 However, only forty-eight of the sixty-one known human Brds (79%) Brds have an Asn residue in the KAc binding pocket. The remaining thirteen Brds (21%) have a tyrosine, threonine or aspartic acid in place of the asparagine (Fig. 1 and ESI Fig. S1B†). To date, no inhibitors of these atypical, non-asparagine bromodomains have been reported.
Fig. 2 PHIP(2) (green sticks and ribbons) in complex with NMP (cyan sticks). The substitution of an Asn for the smaller Thr allows an extra water molecule, the PHIP water (red CPK, PW) to fill the space between T1396 and the usual four Brd waters (red CPK I–IV) (PDB ID: 3MB3).1 |
PHIP(2) belongs to the bromodomain subfamily III alongside CREBBP (Fig. 1).1 The atypical binding site of PHIP(2) is characterised by a threonine (Thr) residue (T1396) in place of the conserved asparagine (Fig. 2). Because the Thr side chain is one atom shorter than Asn, there is space for an additional water molecule, dubbed the ‘PHIP water’, to mediate an interaction between T1396 and the typical water network at the back of the binding pocket. The water-network is otherwise comparable to the network observed in most Brd binding sites, consisting of four water molecules and anchored by a tyrosine residue (Y1353). Despite being identified as a druggable Brd,19 to date no small molecule inhibitors of PHIP(2) have been reported. An unpublished screening campaign by the Structural Genomics Consortium (SGC) did not yield any hits from a library of approximately seven thousand Brd-focused compounds.
Fragment screening has been especially effective in discovering Brd inhibitors using both focused and diverse libraries.23,24 Fragments have been optimized into more potent chemical probes but in all cases the Brds targeted contained the key Asn residue to anchor an acetyl lysine mimetic fragment.11,25
A major challenge of fragment screening is that, although fragment sized compounds make high quality interactions with their targets and bind with high ligand efficiency, overall binding affinity is typically weak due to their small size.26 Therefore the techniques used to detect fragment binding must be sensitive on the μM to mM scale. Biophysical and biochemical solution-based techniques are convenient and widely used but drawbacks include insufficient sensitivity, high rates of false positives and screening concentrations limited to low mM.27 In contrast, NMR and X-ray crystallography are ideal for detecting weak binding by yielding direct structural information, although historically both techniques have been labour intensive and low-throughput. However, in recent years X-ray crystallography has seen an order-of-magnitude speed-up thanks to robotics, improved algorithms and detectors and technical advances at synchrotron facilities.28,29 Moreover, at beamline I04-1 at Diamond Light Source (DLS),30 other recent developments around soaking, harvesting and data analysis have further reduced the effort to such an extent that a 1000-compound fragment screen by crystallography can be completed within a week. As a consequence, X-ray crystallography is now a viable primary screen for academic and industrial researchers.31
In particular, this makes it realistic to screen by soaking single compounds per crystal, rather than cocktails, so that concentrations >100 mM can be achieved, one hundred to one thousand times higher than possible in solution-based techniques, with accordingly increased sensitivity.31
The very weak hits identified by X-ray screening pose a practical problem, in that they are in general too weak to be verified by orthogonal solution assays; yet such information is a common requirement in medicinal chemistry operations, especially if the crystallographic evidence is ambiguous. Optimising weak compounds by traditional SAR is equally undermined without reliable assay readout in the range of the compounds' binding affinity. In any case, progressing weak compounds to potency in general requires changes not only to the periphery, where conventional SAR operates, but also to the core of the compound.
To overcome the problem of optimising weak compounds discovered by high concentration fragment soaking, we report the design and use of a poised fragment screening library to identify hits and the parallel synthesis of analogues to deliver inhibitors with measurable activity. Furthermore, we demonstrate how the method is effective on the previously intractable target, PHIP(2).
In 2011, Roughley and Jordan published an analysis of the most commonly used reactions in drug discovery.32 The authors propose the reactions are regularly used because experienced medicinal chemists know they provide reliable methods to synthesise drug-like molecules. The top ten reactions identified by the analysis, including amide coupling, reductive amination and Suzuki-type aryl–aryl coupling, can be performed using standardised procedures for a wide range of substrates. The starting materials for each reaction contain common functionalities which are well represented in commercial supplier space and are compatible with this diverse range of commercial compounds. The ten reactions were selected to create poised scaffolds with which to perform substructure searches. The poised fragment chemical space was further augmented with twelve heterocycle forming reactions as defined by Hartenfeller et al.33 and an oxazole formation developed in our own lab.34,35 The complete list of poised scaffolds is depicted in Fig. 4A.
Fig. 4 Design of DSPL1 (A) scaffolds used to identify poised fragments. Core scaffolds based upon the most commonly used medicinal chemistry reactions as described by Roughley and Jordan.32 Heterocycle scaffolds based on reactions described by Hartenfeller et al.37 (B) Poised reaction motif vs. cluster. Circle size represents the number of compounds in each group. 407 compounds for DSPL1 were selected from all SGC poised fragments based on diversity and poised reaction motif. |
Encouraged by the success of soaking DSPL1 in PHIP(2) (see below), a second generation poised library was designed (DSPL2). DSPL2 was designed using similar methodology as DSPL1 but was assembled entirely from commercial fragments to ensure that anyone could purchase and use the library. The ZINC fragment-like library was downloaded36 and filtered based on vendor ID for suppliers from which we had routinely sourced compounds, yielding ∼192k compounds. Filters were applied to the compound library to ensure drug-like characteristics (ESI Fig. S2†). From the resulting library of 41271 compounds, 28438 (68%) were found to be poised by our definition. While initially surprising that such a large number of commercial fragments could be synthesised by only twenty-one different chemical reactions, we believe this reflects the very limited coverage of potential synthetic compounds by commercial vendors. Further filters were applied to ensure each compound is compatible with the reaction by which it was found to be poised. For example an amine substituted amide would be removed from the library as an asymmetric diamine synthons could give a mixture of amides if used to synthesise a poised fragment. Compounds without commercially available starting materials were also removed.
In theory, the resulting library of 10448 compounds covers the whole of commercially accessible poised fragment space using our chosen reactions. The library is dominated by amides, with 55.1% of fragments containing a poised amide bond; ether bonds (18.2%) and reductive amination products (10.2%) were also heavily represented. The least represented chemistries were Sonogashira products (0.2%) and sulfonamides (2.0%) (ESI Fig. S3†). Unfortunately, poised heterocyclic reactions were observed in just 2.2% of the compounds, despite thirteen reactions (only eight of which were found to be present) being used to identify such fragments and the importance of such scaffolds in drugs and chemical tools.38 The finding demonstrates the bias of commercially available compounds towards the most commonly used reactions identified by Roughley and Jordan.32
Finally, the USRCAT algorithm was used to ensure chemical diversity was maintained upon selecting 1000 compounds for DSPL2 from the total of >10000.39 USRCAT was selected because it was believed the additional conformation information makes it a superior measure of diversity for library design in comparison to traditional chemical diversity measures (MACCS, Morgan fingerprints, etc.). Although DSPL1 and 2 do not contain the exact same compounds, they were constructed using similar principles and occupy similar chemical space (ESI Fig. S4†). The full identities of DSPL1 and DSPL2 can be downloaded as sdf files from the ESI.†
Due to the relatively potent binding of DMSO to Brds,41 compounds were dissolved in ethylene glycol at a nominal concentration of 400 mM. The supernatant of compounds which did not fully dissolve was used with the assumption it was a saturated solution. The solutions were soaked into protein crystals in crystal buffer in a 1:1 volume ratio to give an approximate final compound concentration of 200 mM. All 406 poised fragments in DSPL1 were soaked into PHIP(2) Brd crystals and when the structures were solved, compounds 1–4 were found to have bound to the protein (1.2% hit rate).
All of the poised fragment hits 1–4 bound to PHIP(2) Brd in the KAc recognition pocket (Fig. 5). The thiourea sulfur of compound 1 displaces water I at the back of the binding pocket allowing the formation of hydrogen bonds with Y1353 and water II. The secondary thiourea nitrogen forms an H-bond to T1396 via the structural PHIP water (PW). The primary thiourea nitrogen lies in the hydrophilic pocket made by waters II, III and IV. The phenyl ring of compound 1 makes an edge–face aromatic interaction with Y1350, which is shifted significantly towards the ligand compared to the NMP-bound conformation (Fig. 2). The ortho-methyl substituent on the phenyl ring makes additional interactions with the hydrophobic residues I1403 and P1340.
Compound 2 binds like a classical acetyl-lysine mimetic (Fig. 5B). The acetyl group lies deep in the binding pocket allowing the methyl group to sit in the hydrophobic pocket created by the four structural waters in a manner analogous to inhibitors of typical Brds and acetyl lysine. The carbonyl oxygen interacts with T1396 through a hydrogen bond mediated by the PHIP water. Like compound 1, the phenyl ring of compound 2 appears to make an edge–face aromatic interaction with the shifted Y1350 side chain, allowing the 2,6-dichloro substituents to occupy the hydrophobic pockets formed by V1345 and I1403 above and below the ring.
Oxazoles 3 and 4 share the same binding mode. The primary amine and oxazole oxygen form a donor and acceptor interaction with T1396 (Fig. 5C, compound 4 not shown for clarity). Y1350 forms an H-bond to the oxazole nitrogen. The orientation of the oxazole rings of compounds 3 and 4 is striking, and allows for the respective lipophilic benzyl and isobutyl groups to lie in the methyl-binding pocket created by the four structural waters. Furthermore, the 4-cyano and 4-amino groups of compounds 3 and 4 interact with a neighbouring protein in the crystal lattice. The cyano group acts as an H-bond acceptor with the backbone N–H of D1352, while the amino group forms a salt bridge with the D1352 carboxylate group.
Scheme 1 Reagents and conditions: (a) amine 1, thio-CDI, DCM, rt, 30 min, then amine 2, rt, 12 h. (b) Triethylamine, DCM, rt, 3 h. (c) 4-Toluenesulfonic acid, NMP, MW irradiation, 120 °C. |
Cmpd | R1 | R3 | Yielda | X-tal soakb | PHIP(2) pIC50d | LE |
---|---|---|---|---|---|---|
a Synthesised using procedure described in Scheme 1. NI: not isolated. b NL: no ligand found in model. Error: experiment failure during soaking, at beamline or during data processing or model refinement. c Binding pose uncertain as a result of partial occupancy. d By AlphaScreen peptide displacement assay. | ||||||
1 | 2-Me | H | 66% | Success | 3.11 ± 0.06 (2) | 0.40 |
5 | 2-Me-3-OMe | H | 71% | Error | 2.97 ± 0.18 (2) | 0.32 |
6 | 2,3-Me2 | H | 12% | Successc | 3.78 ± 0.06 (2) | 0.45 |
7 | 2-Me-3-CF3 | H | 33% | NL | 3.85 ± 0.21 (2) | 0.36 |
8 | 2-Me-3-Cl | H | 26% | Successc | 3.59 ± 0.03 (2) | 0.42 |
9 | 2,6-Me2 | H | 8% | NL | 3.32 ± 0.08 (2) | 0.39 |
10 | 2-Me | Me | 14% | NL | 3.38 ± 0.30 (2) | 0.39 |
11 | 2-Cl | Me | 19% | NL | 3.89 ± 0.05 (2) | 0.45 |
The thiourea series was tested for inhibition of PHIP(2) ligand binding using an AlphaScreen assay based upon a previously reported method.41 Ligand efficiencies (LE) were calculated using the relationship described by Hopkins et al. in which LE is taken to be proportional to the IC50 contribution per non-hydrogen atom in the ligand.26
The original fragment hit 1 had pIC50 of 3.11 in the AlphaScreen assay. Addition of a substituent at the 3-position of the phenyl ring to give compounds 5–8 increased potency. The methyl derivative 6 and the trifluoromethyl derivative 7 had improved pIC50 ∼ 3.8. The addition of a methyl group at the 6-position of the phenyl ring to give compound 9 slightly increased potency to pIC50 3.32 as did methylation of the primary thiourea nitrogen to give compound 10 (pIC50 3.38). The SAR improvements were not cumulative (see ESI,† compounds s1–s20), but replacement of 2-methylphenyl with 2-chlorophenyl on the methylated thiourea, giving compound 11, improved PHIP(2) activity to pIC50 3.89. The LE of fragment hit 1 was calculated to be 0.40. The value was increases to 0.45 for the most potent thioureas 6 and 11.
When the thiourea compounds 5–11 were subjected to soaking experiments with PHIP(2) crystals, two of the most potent compounds, 6 and 8, appeared to bind in the PHIP(2) KAc recognition site (ESI Fig. S5A and B†). The apparent poses are similar to compound 1, displacing structural water.
Cmpd | R1 | R3a | Yieldb | X-tal soakf | PHIP(2) pIC50h | LE |
---|---|---|---|---|---|---|
a R3 = N–H unless stated. b Synthesised using procedure described in Scheme 1 unless stated. c Reagents and conditions: EDC, DMAP, triethylamine, DMF, rt. d Over three steps. See ESI for details. e Triethylamine, AcN, rt. f NL: no ligand found in model. g Binding pose uncertain as a result of partial occupancy. h By AlphaScreen peptide displacement assay. | ||||||
2 | 2,6-Cl2 | C–Me | 44% | Success | <2.30 (2) | <0.25 |
12 | 2,6-(OMe)2 | C–Me | 79% | Success | 3.72 ± 0.04 (2) | 0.35 |
13 | 2,6-Cl2 | C–CH2OH | 15%c | Success | <2.30 (2) | <0.23 |
14 | 2,6-Cl2 | N–Me, C–Me | 83% | Success | <2.30 (2) | <0.23 |
15 | 2,6-(OMe)2 | C–CH2OH | 23%d | Successg | <2.30 (2) | <0.20 |
16 | 2,6-(OMe)2 | 5-Lactam | 6%e | NL | 3.25 ± 0.04 (2) | 0.27 |
17 | 2,6-(OMe)2 | 6-Lactam | 18%e | NL | 3.51 ± 0.04 (2) | 0.27 |
Compounds 2 and 12–17 were tested for inhibition of PHIP(2) ligand binding using the AlphaScreen assay. The original fragment hit 2 was inactive against PHIP(2) at the highest possible assay concentration and would not have been detected as a hit if AlphaScreen had been used as the primary screen. Attempts to modify the acetyl amide group or the benzene ring gave no improvements in potency. The exception was the 2,6-dimethoxy analogue 12 with a pIC50 of 3.72 and LE of 0.35.
The benzyl amide compounds 12–17 were soaked into PHIP(2) crystals. Surprisingly the inactive compounds 13 and 14 were found to bind in the KAc binding site as well as the active compound 12. Compound 12 has the same binding mode as fragment hit 2, with the 2,6-dimethoxy substitution pattern on the benzene ring occupying the V1345 and I1403 lipophilic pockets above and below the ring (Fig. 6A).
Alcohol 13 binds in the same position as fragment hit 2, with the additional hydroxyl group making an H-bond interaction to the protein backbone at P1340 (Fig. 6B). In contrast, the N-methyl group of the tertiary amide 14 forms a hydrophobic interaction with the P1340 sidechain (Fig. 6C).
Despite binding in a similar way to the other amide compounds, tertiary amide 14 appears to displace structural water II.
The activity and crystallographic data resulted in three new compounds being proposed for synthesis. All three were synthesized in low yield within two weeks (Table 2, compounds 15–17). The 2,6-dimethoxy analogue of alcohol 13 (compound 15) is inactive against PHIP(2) but appears to bind to the protein when subjected to soaking experiments (ESI Fig. S5C†).
Overlay of fragment hit 2 with the previously reported NMP-bound structure (Fig. 2) showed good overlap between the amide binding modes for each compound. A scaffold merging approach was adopted and the lactam containing compounds 16 and 17 proposed and synthesised. Both lactams were found to inhibit the activity of PHIP(2). pIC50’s of 3.25 and 3.51 for compounds 16 and 17 respectively demonstrate good activity but LE of only 0.27. The best compound of the amide series, compound 12, has a calculated LE of 0.35.
Cmpd | R1 | Yielda | X-tal soakc | PHIP(2) pIC50d | LE |
---|---|---|---|---|---|
a Synthesised using procedure described in Scheme 1 unless stated. b Over three steps. See Spencer et al. for details.34 c NL: no ligand found in model. Error: experiment failure during soaking, at beamline or during data processing or model refinement. d By AlphaScreen peptide displacement assay. | |||||
3 | Bn | 62% | Success | 3.23 ± 0.10 (2) | 0.30 |
4 | iBu | 42% | Success | 3.57 ± 0.10 (2) | 0.42 |
18 | Prc | 93% | NL | 3.48 ± 0.03 (2) | 0.44 |
19 | iPr | 54% | NL | <2.30 (2) | <0.29 |
20 | 4-Cl-Bn | 46% | Error | 3.95 ± 0.03 (2) | 0.35 |
21 | Hexc | 71% | NL | 3.31 ± 0.06 (2) | 0.33 |
22 | Hexc | 30%b | NL | 3.71 ± 0.16 (2) | 0.29 |
In contrast, the assessment of the follow-up compounds requires different criteria. Our observation that some of the more potent compounds such as 11 are not visible in the crystal structure upon soaking is line with common anecdotal experience, that crystallographic clarity is not a useful ranking criterion for micromolar compounds. In the fragment field, this manifests as the widely-reported lack of hit overlap between techniques.42 A variety of factors influence whether soaking is successful, including compound solubility in the crystallization buffer.
Our in-house poised fragments were used in combination with a novel medium-throughput crystallographic screen to identify the first reported inhibitors of PHIP(2), an atypical bromodomain. By exclusively utilising poised fragments in our screen, we have rapidly generated follow-up leads and shown increases in compound activity alongside highly detailed structural data.
The thiourea 1 displaces one of the structural waters at the back of the PHIP(2) binding site, something rarely seen in other bromodomain ligand–protein complexes. The series was expanded to give highly ligand efficient inhibitors with IC50's of 100–200 μM and calculated LE values greater than 0.40 (compounds 6, 7, and 11).
Amide 2 was found to be inactive in an AlphaScreen assay against PHIP(2), yet one round of poised synthesis yielded compound 12 with an IC50 of 190 μM and LE 0.35. Structural data from AlphaScreen-inactive compounds show displacement of structural waters and allowed further development of the series. The 5- and 6-member lactams, 16 and 17, showed sub-millimolar activity and provide a number of additional vectors for elaboration into unexploited regions of the PHIP(2) binding pocket.
The binding mode of the aminooxazoles 3 and 4 appears to be driven by interactions with T1396, with no displacement of the structural waters or the PHIP water deeper in the binding pocket. The 2-isobutyl and 2-cyclopropyl oxazoles, 4 and 18, show excellent ligand efficiency, with IC50's 226 μM and 329 μM respectively and LE > 0.40.
Ongoing studies are focusing on improving the affinities of the PHIP(2) hit series reported here with the aim of developing a viable PHIP(2) chemical probe to help further explore the role of this Brd in disease, notably cancer. In parallel, through the SGC's network of synthetic chemistry collaborators, we are actively expanding our poised fragment library beyond the simple chemistry described here with compounds derived from sp3 rich, stereochemically controlled poised reactions.
We submit that the use of a crystallographic primary screen followed by rapid poised chemistry to generate a follow-up library offers a new, powerful method in hit discovery and lead series selection which can be utilised by both academic and industrial researchers. The success in finding efficient hits for the low druggability PHIP(2) bromodomain indicates the power of this method in addressing difficult protein–protein interaction targets.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c5sc03115j |
This journal is © The Royal Society of Chemistry 2016 |