Christoph
Selg
a,
Andrea
Schöler
a,
Julian
Schliehe-Diecks
b,
Maria
Hanl
c,
Laura
Sinatra
a,
Arndt
Borkhardt
b,
Menyhárt B.
Sárosi
d,
Sanil
Bhatia
b,
Evamarie
Hey-Hawkins
d and
Finn K.
Hansen
*c
aInstitute for Drug Discovery, Medical Faculty, Leipzig University, Brüderstraße 34, 04103 Leipzig, Germany
bDepartment of Pediatric Oncology, Hematology and Clinical Immunology, Medical Faculty, Heinrich-Heine University Düsseldorf, Düsseldorf, Germany
cPharmaceutical Institute, Department of Pharmaceutical and Cell Biological Chemistry, University of Bonn, An der Immenburg 4, 53121 Bonn, Germany. E-mail: finn.hansen@uni-bonn.de
dInstitute of Inorganic Chemistry, Faculty of Chemistry and Mineralogy, Leipzig University, Johannisallee 29, 04103 Leipzig, Germany
First published on 4th August 2021
The elevated expression of histone deacetylases (HDACs) in various tumor types renders their inhibition an attractive strategy for epigenetic therapeutics. One key issue in the development of improved HDAC inhibitors (HDACis) is the selectivity for single HDAC isoforms over unspecific pan inhibition to minimize off-target toxicity. Utilizing the carborane moiety as a fine-tuning pharmacophore, we herein present a robust solid phase synthetic approach towards tailor-made HDACis meeting both ends of the selectivity spectrum, namely pan inhibition and highly selective HDAC6 inhibition.
Fig. 1 Pharmacophore model of selective HDACis: zinc-binding group (ZBG), hydrophobic alkyl or aryl linker and a sterically demanding, hydrophobic cap group.48,56 |
Except for class III HDACs, which are NAD+-dependent, all HDACs are Zn2+-dependent enzymes.20 Thus, the respective inhibitor usually possesses a zinc-binding group (ZBG), a hydrophobic linker and a cap group (Fig. 1, top left). Starting from the rather unselective FDA-approved hydroxamate-based drugs vorinostat, belinostat and panobinostat,20 within this pharmacophore model, HDAC6 isoform selectivity was shown to be connected to steric bulk and hydrophobicity of the cap group owing to the increased size of the outer channel opening of HDAC6 with 17 Å compared to 12 Å in HDAC1.17,21–23 High selectivity was achieved with cap groups bearing extended aromatic residues as in tubacin and tubastatin A as well as bulky cap groups as realized in ricolinostat-type HDACis.14,22,23 In order to increase selectivity, numerous hydrophobic (hetero)aromatic and, to some extent aliphatic cap groups were described in the literature.20,21 These carbon-based systems are for the most part limited to a two-dimensional shape as larger carbocyclic skeletons like adamantyl-, norbornenyl-, bicyclo-[2.2.2]octanyl- or fullerenyl-residues pose problems due to limited modification sites and possible rearrangement reactions or their sheer size.25 Therefore, when it comes to hydrophobic, sterically demanding scaffolds, an emerging building block in drug design is the carborane moiety.26 Replacement of (hetero)aryl moieties by these extremely versatile boron clusters has been widely applied and often resulted in dramatically increased potency, target selectivity and lower toxicity.25,27 Furthermore, carboranes play a key role in the development of boron neutron capture therapy (BNCT) agents and as carriers in nanomedicinal applications.29,30 Carboranes are a family of polyhedral clusters comprised of CH and BH vertices with the icosahedral closo-dicarbadodecaborane (C2B10H12) being the most prominent and well-examined member.27,28 For symmetrical reasons, the two carbon atoms can be arranged within the icosahedron separated by one (ortho-/1,2-), two (meta-/1,7-) or three (para-carborane/closo-1,12-dicarbadodecaborane) bonds (Fig. 1, bottom). The size of the cluster is in between a rotating phenyl ring and adamantane. The cluster fragments CH and BH are interconnected by a strong network of multi-electron-multi-center bonds and thus represent a completely delocalized three-dimensional σ-aromatic compound (deemed “superaromatic”) of very high stability.31 Bearing acidic (carbon-centered) and hydridic (boron-centered) hydrogen atoms, carboranes exhibit a peculiar interaction profile where besides the classic hydrophobic interactions, also hydrogen bonds C–H⋯X (X = O, N, S, F, π-system) and the so-called dihydrogen bonds B–H⋯H–X (with χX > χH > χB) are formed.32,33 As opposed to carbon cage structures, for carboranes manifold achiral, chiral and chiral-at-cage modifications at carbon- as well as at boron-vertices are described in the literature, thus paving the way also for the evaluation and design of target-specific asymmetric drugs.34,35
While boron-containing drugs like bortezomib, crisaborole and tavaborole are established therapeutics, carborane drugs are still scarce in clinical trials.30,36–38 Within this branch, a limited number of boron-containing HDACis have been investigated.39,40 However, only one of them was carborane-based and to the best of our knowledge none of them were tested for their HDAC6 selectivity. In the present publication we therefore seek to establish a systematic approach to carborane-based HDACis to provide a solid foundation for further investigations. By synthesizing a small compound library, we want to enable the direct head-to-head comparison of different carborane cap groups employing a set of known and successfully employed linker groups in combination with the well-investigated aryl or alkyl hydroxamate zinc-binding group. The hit compounds will be further compared to their respective phenyl analogues.
Following the synthetic protocol of Nekvinda and co-workers,43closo-1,2-dicarbadodecaborane (ortho-1) was monolithiated and treated with oxirane to afford hydroxymethyl carborane 2 which was subsequently oxidized with CrO3 to form the desired carboxymethyl carborane 3. As benzoic acid analogues of ortho-carborane prove problematic in amide couplings, we chose the meta isomer, closo-1,7-dicarbadodecaborane (meta-1), for our second cap group scaffold.44 The respective carboxylic acid 4 was obtained following a slightly modified literature procedure of Kasar and co-workers.46 Monolithiation of meta-1 and subsequent treatment with gaseous CO2 afforded 4 in quantitative yield (Scheme 1).
With the carborane cap group precursors 3 and 4 in hand, an elaborate synthetic approach for the linker and hydroxamate moieties was developed. To overcome the problems of the aforementioned possible deboronation that might occur under strongly basic conditions associated with the typical synthesis of hydroxamic acids through hydroxylaminolysis,40 we envisioned that a solid-phase synthesis as recently described by our group would be the suitable tool.47 Beside the prevention of harsh basic conditions, this methodology could offer a fast and efficient carborane-based HDACi library expansion through a parallel synthesis approach. Consequently, we treated in a first step the commercially available 2-chlorotritylchloride resin 5 with N-hydroxyphthalimide to form the modified resin 6. The phthaloyl residue in turn is cleaved using a 5% methanolic hydrazine solution to obtain a free hydroxylamine which can be coupled in a HATU-mediated amide coupling with an N-Fmoc protected linker to form the modified resins 7a–d (Scheme 2).
To cover a broad spectrum of possible HDACis, we chose linker moieties that were already employed in unselective and selective HDACis described in the literature in resemblance to ricolinostat, vorinostat and tubacin, two different chain lengths, namely C5 and C6 alkyl linkers (7a, 7b), as well as a tubastatin-type benzyl (7c) and a panobinostat-type vinylbenzyl linker (7d) were selected (Scheme 2).
After Fmoc deprotection of the modified resins 7a–d, amide coupling with the respective carboranyl carboxylic acid cap groups 3 and 4 and subsequent cleavage with TFA afforded the desired carborane-capped hydroxamates 8a–d and 9a–d in very good yields (Scheme 3).
Scheme 3 Fmoc deprotection of the modified resins 7a–d, amide coupling with the carboranyl cap groups and cleavage off the resin. Reagents: (a) 20% piperidine in DMF, rt; (b) 3, HATU, DIPEA, DMF, 18 h, rt, (c) 4, COMU, DIPEA, DMF, 18 h, rt; (d) 5% TFA, DCM, 1 h, rt, 74–99%. For definition of X see Scheme 2. |
While coupling of ortho-carborane derivative 3 was achieved flawlessly under standard conditions using HATU, DIPEA and DMF, the meta-carboranyl carboxylic acid 4 could only be coupled in moderate yields. However, following the procedure described by Scholz et al., by the replacement of HATU with the more active COMU, excellent yields were also achieved with the sterically demanding and less reactive 4.44b
With the eight hydroxamates 8a–d and 9a–d in hand, we set out to explore their inhibitory activity against HDAC1 and HDAC6 in a fluorogenic assay using ZMAL (Z-Lys(Ac)-AMC) as substrate and vorinostat (suberoylanilide hydroxamic acid, SAHA) as a control (Table 1) following our previously published protocol.45 To our delight, all new compounds showed IC50 values for HDAC6 in the double- or even single-digit nanomolar range comparable to those of SAHA (IC50: 0.031 μM). Two of the compounds (9a and 9d) showed improved HDAC6 inhibition compared to the control. Compound 9d was identified as the most potent HDAC6 inhibitor (IC50: 0.006 μM).
Structure | Target inhibition IC50a [μM] | MTT IC50b [μM] | Structure | Target inhibition IC50a [μM] | MTT IC50b [μM] | ||||
---|---|---|---|---|---|---|---|---|---|
HDAC1 | HDAC6 | SF6/1 | A2780 | HDAC1 | HDAC6 | SF6/1 | A2780 | ||
a Assays carried out with n ≥ 2 (each in duplicate wells), values are shown as means ± SD. b Assays carried out with n ≥ 3 (each in duplicate wells), values are shown as means ± SD. | |||||||||
0.451 ± 0.023 | 0.037 ± 0.006 | 12 | 5.78 ± 1.68 | 0.440 ± 0.166 | 0.024 ± 0.005 | 18 | 3.98 ± 0.95 | ||
0.101 ± 0.003 | 0.045 ± 0.011 | 2 | 4.64 ± 0.95 | 1.328 ± 0.039 | 0.078 ± 0.0003 | 17 | 5.28 ± 1.25 | ||
0.015 ± 0.002 | 0.031 ± 0.004 | 0.5 | 1.66 ± 0.81 | 0.081 ± 0.014 | 0.040 ± 0.008 | 2 | 3.07 ± 1.41 | ||
2.888 ± 0.049 | 0.059 ± 0.007 | 49 | 4.85 ± 1.64 | 1.715 ± 0.002 | 0.006 ± 0.0001 | 286 | 4.60 ± 1.67 | ||
0.176 ± 0.007 | 0.095 ± 0.008 | 2 | 1.39 ± 0.30 | 0.725 ± 0.041 | 0.025 ± 0.002 | 29 | 5.04 ± 1.21 | ||
0.094 ± 0.014 | 0.031 ± 0.009 | 3 | 1.03 ± 0.27 |
Due to a relatively even distribution of IC50 values no clear structure–activity relationship could be assigned for HDAC6 inhibition. For the activity against HDAC1, however, a clear selectivity trend was observed (Fig. 2).
Three of the compounds bearing an alkyl chain linker (8a, 9a and 9b) showed moderate (12- to 18-fold) preference for HDAC6 with slightly better inhibitory activity for the longer chains. The results of the short chain derivative 8b (SF6/1: 2) whose net chain length and carbonyl group positions match the parent aryl analogue pan inhibitor SAHA (SF6/1: 3) showed no clear preference.
Also, for the vinylbenzyl hydroxamates 8c and 9c, inhibition of both HDACs was very similar. Notably, inverse preference was observed between the two carborane cap groups shifting from either two-fold HDAC1 preference for 8c to two-fold HDAC6 preference for 9c. In contrast to this unselective inhibitory activity, the tubastatin A analogues with a benzyl linker (8d and 9d) showed excellent selectivity for HDAC6 whilst also maintaining very high inhibitory activity. Thus, we identified the meta-carboranyl hydroxamate 9d as the hit compound with an IC50 value of 0.006 μM and a more than 280-fold selectivity for HDAC6 (Table 1 and Fig. 2).
To investigate the influence of the carborane moiety, we synthesized aryl analogues for the best pan inhibitory compound 8c and the best selective HDAC6 inhibitor 9d. The synthesis was carried out employing our solid-phase protocol to afford vinylbenzyl derivative 10 from phenylacetic acid and benzyl derivative 11 using benzoic acid. Both phenyl analogues 10 and 11 showed good inhibitory activity against HDAC6 with IC50 values of 0.095 μM and 0.025 μM, respectively (Table 1). The HDAC6 selectivity (SF6/1: 2 and 29 resp.) was distributed the same way as observed for parent compounds 8c and 9d, yet the carborane compounds clearly outperformed their corresponding phenyl analogue in terms of HDAC6 inhibition and, in the case of 9d, selectivity.
To evaluate the anticancer properties of the carborane-capped HDACis, compounds 8a–d and 9a–d, as well as their phenyl analogues 10 and 11 were tested in MTT assays in the ovarian cancer cell line A2780, again using the FDA-approved HDACi vorinostat (SAHA) as positive control.45 The results are summarized in Table 1. All compounds demonstrated antiproliferative effects in the single-digit micromolar concentration range. The most potent carborane-capped HDACi was compound 8c with an IC50 value of 1.66 μM.
Encouraged by these promising results, we performed additional HDAC isoform profilings examining the HDAC6 selectivity of the two most potent compounds 8c and 9d over all class I isoforms with SAHA as control and in the case of HDAC8, also panobinostat (Table 2).
Compound | Target inhibition IC50 [μM]a | ||||
---|---|---|---|---|---|
HDAC1b | HDAC2b | HDAC3b | HDAC6b | HDAC8c | |
a Assays carried out with n ≥ 2 (each in duplicate wells), values are shown as means ± SD. b Z-Lys(Ac)-AMC was used as substrate. c Boc-Lys-(Tfa)-AMC was used as substrate. n.d.: not determined. | |||||
8c | 0.015 ± 0.002 | 0.020 ± 0.004 | 0.021 ± 0.001 | 0.031 ± 0.004 | 1.661 ± 0.126 |
9d | 1.715 ± 0.002 | 1.857 ± 0.167 | 1.346 ± 0.003 | 0.006 ± 0.0001 | 5.780 ± 0.192 |
SAHA | 0.094 ± 0.014 | 0.163 ± 0.015 | 0.113 ± 0.004 | 0.031 ± 0.009 | 3.410 ± 0.494 |
Panobinostat | n.d. | n.d. | n.d. | n.d. | 0.166 ± 0.015 |
Considering HDAC2 and HDAC3, again, 8c showed non-selective behavior while compound 9d exhibited substantial selectivity for HDAC6 over HDAC2 and HDAC3 (309- and 224-fold, respectively). The IC50 values for HDAC8 were greater than 1 μM and indicated neglectable inhibition for both compounds.
As already observed for tubastatin A and panobinostat, the potent HDAC6 inhibition may be the result of π–π stacking interactions between the linker aryl groups and the adjacent phenylalanine residues inside the hydrophobic channel.48,54 To further undermine this assumption, compounds 8c and 9d were docked into the human HDAC6 catalytic domain 2 (hHDAC6-CD2, Fig. 3A and B). Hydroxamate inhibitors are known to coordinate to the HDAC Zn2+ ion through mono- or bidentate binding modes.48,54,55 A monodentate Zn2+-binding was predicted for 8c and the amido NH group forms a hydrogen bond with S568. The phenyl ring of 8c forms a π–π stacking interaction with F680. On the other hand, a bidentate Zn2+-binding was predicted for 9d. The phenyl ring of 9d forms a π–π stacking interaction with F620 and the amido CO group can potentially form a hydrogen bond with S568.
Fig. 3 8c (A) and 9d (B) docked into the human HDAC6 catalytic domain 2 (PDB ID: 5EDU).48 HDAC6 is shown in gray. Residue numbering according to PDB ID: 5EDU. 8c (C) docked into human HDAC1 (PDB ID: 5ICN).49 HDAC1 is shown in grey. Residue numbering according to PDB ID: 5ICN. Superimposed HDAC6 backbone trace is shown in light blue for comparison (PDB ID: 5EDU). B: orange, C: gray/green, N: blue, O: red, Zn: purple. Hydrogen atoms are omitted. Ligand–receptor interaction and Zn2+ coordination is shown as dashed lines. |
The carborane clusters of both compounds form dihydrogen bonds with H500 and P501, but the ortho-carborane of 8c additionally interacts with F620. The dihydrogen bonds between the ligands and key hHDAC6-CD2 residues H500 and P501 certainly contribute to the observed potent inhibition of HDAC6. The sterically demanding hydrophobic carborane cap groups point outwards the hydrophobic channel and align to HDACs surface binding domain which leaves sufficient space for the carborane to act as isoform selection parameter without decreasing target inhibition.
As an attempt to understand the different selectivity trends, 8c and 9d were docked into human HDAC1.49 74% of the docked poses of 8c contained a monodentate Zn2+-binding. Fig. 3C shows a monodentate binding mode where the ortho-carborane interacts with P29 and the amido NH group forms a hydrogen bond with D99. In the case of 9d, only 39% of the docked poses chelated the zinc ion in the active site. The remaining poses of 9d did not enter the HDAC1 active site (see ESI, Fig. S1†). As a comparison, 100% of the docked HDAC6 poses for both 8c and 9d were chelating the zinc ion. It is reasonable to assume that compound 8c exhibiting a longer linker is capable of chelating the zinc ion of both HDAC1 and HDAC6, thus inhibiting both targets, whereas the shorter benzyl linker realized in 9d in combination with a sterically demanding carborane cap is unable to efficiently engage the zinc ion in HDAC1. With the combination of our linker library and meta- and ortho-carboranes as bioisosteric phenyl mimetics, thus, we could not only significantly improve the potency, but, in the case of 9d also the selectivity of the HDACis compared to their parent phenyl analogues.
In the next step, we investigated the selectivity profile of 8c and 9d in a cellular environment using the corresponding phenyl analogues 10 and 11 as well as vorinostat (SAHA) as control compounds. The acute myeloid leukemia (AML) cell line HL60 was treated for 24 h with 1 μM of each compound and the cell lysates were subsequently immunoblotted with antibodies against acetyl-histone H3 (a marker of HDAC1-3 inhibition), acetyl-α-tubulin (a marker of HDAC6 inhibition), and acetyl-SMC3 (a marker of HDAC8 inhibition). The results are presented in Fig. 4A. None of the compounds induced acetylation of SMC3. The unselective inhibitor 8c caused hyperacetylation of histone H3 indicating inhibition of HDAC1, 2 and/or 3. As expected, the highly selective HDAC6i 9d increased only the protein levels of acetyl-α-tubulin. Consequently, the biochemical and western blot data confirm that 9d is a selective HDAC6i both in a cell-free and cellular setting.
To elucidate the antiproliferative effects of compounds 8c and 9d against a selection of cell lines from different leukemia entities, we screened both compounds in a CellTiter Glo assay against the HPBALL (T-cell acute lymphoblastic leukemia, T-ALL), K562 (chronic myeloid leukemia, CML) and HL60 (AML) cell lines. The corresponding phenyl analogues 10 and 11 were used as controls (Fig. 4B).
Compound 9d displayed moderate antiproliferative activity with IC50 values ranging from 2.43 to 12.94 μM, while 8c demonstrated submicromolar activity against all three cell lines with the highest activity observed in HL60 cells (IC50: 0.12 μM). Of note, 8c outperformed its corresponding phenyl analogue 10 in all three cell lines and showed 7–8-fold lower IC50 values.
The relatively low cytotoxicity of the selective HDAC6i 9d in comparison to the unselective HDACi 8c is in good agreement with recent results from literature, which indicate that class I HDAC inhibition is important for single-agent activity.51 HDAC6 together with the motor protein dynein is essential to recruit ubiquitinated proteins to aggresomes. The dual inhibition of the aggresome (via HDAC6 inhibition) and proteasome pathways can consequently lead to accumulation of misfolded, cytotoxic proteins resulting in induction of apoptosis.52 Hence, HDAC6i and proteasome inhibitor exhibit synergistic anticancer properties and this synergism led to the approval of the combination therapy of panobinostat, the proteasome inhibitor bortezomib, and dexamethasone to treat multiple myeloma. This knowledge prompted us to investigate the synergistic effects of 8c, 9d and vorinostat as a control in combination with bortezomib, using 8 × 8 dose–response matrices (Fig. 5). Among the tested combinations, 8c displayed the highest synergy (ZIP synergy score: 10.52) with bortezomib, followed by vorinostat (ZIP synergy score: 7.18) and 9d (5.95). The comparative high synergistic activity of 8c with bortezomib is potentially due to its potent class I HDAC inhibitory activity, as targeting especially class I HDACs are critical targets to overcome bortezomib resistance.53
Fig. 5 (A)–(C) Illustrative synergy map of 8c, 9d and vorinostat (SAHA) after 72 h co-treatment of the acute myeloid leukemia (AML) cell line HL60 with bortezomib (proteasome inhibitor) at depicted concentrations. The mean synergy score calculations were based on the ZIP model and visualizations were performed using Synergy Finder webtool.50 |
To conclude, due to the tunable selectivity profile, this series of carborane-capped HDACis represents an encouraging starting point to develop unselective HDACis with improved single-agent activity as well as highly selective HDAC6 inhibitors for potential combination therapies.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1sc02268g |
This journal is © The Royal Society of Chemistry 2021 |