Open Access Article
Marion
Poirier‡
a,
Jonai
Pujol-Giménez‡
bcd,
Cristina
Manatschal
e,
Sven
Bühlmann
a,
Ahmed
Embaby
a,
Sacha
Javor
a,
Matthias A.
Hediger
*bcd and
Jean-Louis
Reymond
*a
aDepartment of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, 3012 Bern, Switzerland. E-mail: jean-louis.reymond@dcb.unibe.ch
bInstitute of Biochemistry and Molecular Medicine, University of Bern, Bühlstrasse 28, 3012 Bern, Switzerland
cMembrane Transport Discovery Lab, Department of Nephrology and Hypertension, Inselspital, University of Bern Kinderklinik, Freiburgstrasse 15, 3010 Bern, Switzerland. E-mail: matthias.hediger@ibmm.unibe.ch
dDepartment of Biomedical Research, University of Bern, Murtenstrasse 35, 3008 Bern, Switzerland
eDepartment of Biochemistry, University of Zürich, Winterthurerstrasse 190, Zürich, Switzerland
First published on 2nd June 2020
Solute carrier proteins (SLCs) control fluxes of ions and molecules across biological membranes and represent an emerging class of drug targets. SLC11A2 (hDMT1) mediates intestinal iron uptake and its inhibition might be used to treat iron overload diseases such as hereditary hemochromatosis. Here we report a micromolar (IC50 = 1.1 μM) pyrazolyl-pyrimidone inhibitor of radiolabeled iron uptake in hDMT1 overexpressing HEK293 cells acting by a non-competitive mechanism, which however does not affect the electrophysiological properties of the transporter. Isothermal titration calorimetry, competition with calcein, induced precipitation of radioactive iron and cross inhibition of the unrelated iron transporter SLC39A8 (hZIP8) indicate that inhibition is mediated by metal chelation. Mapping the chemical space of thousands of pyrazolo-pyrimidones and similar 2,2′-diazabiaryls in ChEMBL suggests that their reported activities might partly reflect metal chelation. Such metal chelating groups are not listed in pan-assay interference compounds (PAINS) but should be checked when addressing SLCs.
Two families of small molecule hDMT1 inhibitors have been reported in the literature as the results of high-throughput screening campaigns, namely bis-cationic isothioureas such as dibenzofurans 1 and mesitylene 2,14 as well as pyrazolyl-pyridine 3 (Fig. 1).15 In our own investigations on hDMT1, we used the inhibitors mentioned above as seeds for a ligand-based virtual screening campaign guided by 3D-shape and pharmacophore similarity16 and discovered bis-isothiourea 4 and pyrazolyl-pyrimidone 5 as two additional hDMT1 inhibitors.17 Kinetic studies, an X-ray structure of a related brominated bis-isothiourea inhibitor in complex with a bacterial analog of the transporter and mutational studies recently showed that bis-isothiourea-based compounds act as competitive inhibitors of hDMT1.18 On the other hand, pyrazolyl-pyrimidone 5, whose small size and better drug-like properties made it an attractive candidate compound, acted as a non-competitive inhibitor. Here we set out to investigate the pyrazolyl series (3 and 5) closer and understand its mechanism of action.
![]() | ||
| Fig. 1 Structure and reported activity of DMT1 inhibitors. Data for inhibition of Fe2+ uptake into transfected HEK cells measured by calcein fluorescence assay (1–3) from ref. 14 and 15 or radioactive 55Fe2+ uptake assay (4 and 5) from ref. 17. | ||
![]() | ||
| Fig. 2 Chemical structures and X-ray crystallography of analogs of pyrazole 5. Carbon is in gray, nitrogen in blue, oxygen in red and chlorine in green. | ||
We measured hDMT1 activity using a radioactive iron uptake assay in HEK293 cells stably overexpressing hDMT1 as described previously.17 Dibenzofuran inhibitor 1 showed only micromolar potency in this assay, which is 5-fold weaker than the originally reported value based on a calcein assay, while the mesitylene inhibitor 2 showed a comparable value to the originally reported inhibition (Table 1). For the present study, we selected dibenzofuran inhibitor 1 as positive control. However, we were unable to determine an IC50 value for the known pyrazolyl-pyridine 3 due to a lack of inhibition at high compound concentration.
| Compound | IC50a (μM) | Ligand efficiency |
|---|---|---|
| a IC50 values were calculated using the radiolabeled iron uptake assay in HEK293 cells (source: American Type Culture Collection, catalog no. CRL-1573) stably overexpressing hDMT1. Absorbed radioactive iron (1 μM) was measured after 15 min of incubation in the presence of the indicated compound at different concentrations at extracellular pH 5.5. b Partial inhibition was detected at 10 μM, however apparent uptake increased at higher concentrations. Compounds 8–12, 14–15, 20–25 showed less than 20% inhibition at 10 μM in this assay. | ||
| 1 | 1.76 ± 0.06 | 0.35 |
| 2 | 0.35 ± 0.04 | 0.25 |
| 3 | ≫10 μMb | |
| 5 | 13.1 ± 0.7 | 0.46 |
| 6 | 10.0 ± 0.4 | 0.44 |
| 7 | 4.2 ± 0.4 | 0.40 |
| 13 | 1.1 ± 0.01 | 0.45 |
| 16 | 3.8 ± 0.2 | 0.45 |
| 17 | 2.26 ± 0.3 | 0.42 |
| 18 | 0.94 ± 0.2 | 0.42 |
| 19 | 12.5 ± 0.3 | 0.43 |
Activity screening of the twenty synthesized analogs of 5 for inhibition of hDMT1 revealed that seven pyrazoles (6, 7, 13, 16–19) showed similar or improved inhibition compared to 5. On the other hand, pyrazoles 8–12, 14–15 and 20 as well as pyrazolones 21–25 did not show significant inhibition of iron uptake in our assay. Across all compounds tested, pyrazole 13 stood out as one of the most potent compounds with IC50 = 1.1 ± 0.01 μM, representing a 10-fold improvement over our initial inhibitor 5 while preserving the original ligand efficiency. Note that the bis-cyclohexane analog 18 was slightly more potent, however at the cost of higher hydrophobicity.
On the other hand, the uptake of radioactive iron in Xenopus oocytes overexpressing hDMT1 was totally inhibited by the bis-thiourea inhibitor 1, while pyrazolyl-pyrimidone 13 only inhibited 50% of the uptake (Fig. 3a). Electrophysiological recordings in this system furthermore showed that only inhibitor 1 blocked hDMT1-induced current with a potency comparable to the iron uptake experiment (IC50 = 0.14 μM), but that pyrazolyl-pyrimidone 13 had no measurable effect on the recorded currents (Fig. 3b/c).
Taken together, these experiments showed that 13 was indeed a more potent version of the non-competitive inhibitor 5. However, in contrast to the well-behaved bis-isothiourea 1, the inhibitory effect of 13 was only detected in iron uptake experiments and was not observed in electrophysiological recordings. A similar effect had been reported for the literature inhibitor 3.15
A first indication that pyrazolyl-pyrimidone 13 might interact with divalent metal ions was obtained by isothermal titration calorimetry (ITC) with Cd2+, which we used as a redox stable analog of Fe2+ less prone to hydroxide formation at pH = 7.4 and which is also transported by DMT1. At a pH of 7.4, the ITC data indicated an apparent dissociation constant of KD = 6.4 μM with n = 0.54, consistent with the formation of a 2:1 ligand–metal complex (Fig. 4a). At pH 5.5 under which the hDMT1 transport experiment is carried out, the ITC data could not be interpreted with certainty. Upon titration of Cd2+ into a solution of 13, the exothermic signal seems to saturate after the first injections, but later results in a noisy signal, indicating precipitation (Fig. 4b). The less potent pyrazolyl-pyrimidone DMT1 inhibitor 5 also interacted with Cd2+ at pH 7.4, although binding was much weaker (KD ∼ 250 μM, n = 0.74).
A calcein competition assay to detect Fe2+ chelation at pH 5.5 also indicated a significant level of iron chelation by 13, although the effect was weaker than with the positive control bipyridine or the literature inhibitor 3, and pyrazolyl-pyrimidone 5 did not show any iron chelation in this assay (Fig. 4c). Tracking the radioactivity of 55Fe2+ in solutions after incubation with the various inhibitors indicated that pyrazolyl-pyrimidone 13 as well as the literature inhibitor 3 induced significant precipitation of Fe2+ at pH 7.4, while bis-isothiourea 1, bipyridine or pyrazolyl-pyrimidone 5 had no effect. On the other hand, at pH 5.5, we only observed precipitation with the literature inhibitor 3 corresponding to the pH of the assay with hDMT1 (Fig. 4d and e). Such precipitation might indicate chelation of iron to form an insoluble complex. Taken together, these data pointed to potential iron chelating abilities of our inhibitors, however the evidence was inconclusive and did not match well with the observed inhibitory potencies.
Indeed, activity screening showed that bipyridine, a well-known iron chelator, non-specifically inhibited both transporters. Similarly, the literature pyrazolyl-pyridine 3, our optimized pyrazolyl-pyrimidone 13, and to a lesser extent pyrazolyl-pyrimidone 5, all inhibited both transporters, although the inhibition by 3 was not consistently observed at higher concentrations, probably due to adsorption of the complex to the cell material and assay plate, resulting in apparent uptake. On the other hand, bis-isothiourea 1, which inhibits hDMT1 by competing in the iron binding site, only inhibited hDMT1 and had no measurable activity on hZIP8 (Fig. 5a). A precise determination of the inhibitory effects with pyrazolyl-pyrimidone 13 and with the known iron chelator bipyridine showed that these compounds had comparable IC50 values on both transporters (13: IC50 ∼ 1.1 μM, bipyridine: IC50 ∼ 6.5 μM, Fig. 5b/c).
However, while all pyrazolyl-pyrimidones and pyrazolyl-pyridines tested here should have comparable metal chelating abilities, only a few of them actually showed activity in the iron uptake assays with hDMT1 and hZIP8. We therefore postulate that inhibition only occurs either if metal complexation is tight, as for the reference chelator bipyridine, or if metal complexation results in the formation of an insoluble complex, as observed using radioactive 55Fe2+ with 3 and 13. Note that, while iron uptake into hDMT1 expressing cells was partially reduced with the literature inhibitor 3, we did not observe inhibition at higher concentrations (Table 1), probably because complexed iron was in this case absorbed on the cell material and assay plate, and resulted in apparent iron uptake. This is in line with the Fe2+ precipitation experiments at pH 5.5, where we only observe high precipitation of iron in the presence of compound 3. The lack of inhibition by most of the compounds tested might therefore either reflect partial iron chelation or formation of an iron complex that partially absorbs on the cells, resulting in apparent iron uptake and no net inhibition.
585 molecules) contain a potentially metal chelating 2,2′-diazabiaryl group, including 1717 pyrazolyl-pyridines related to 3 and 169 pyrazolyl-pyrimidones related to 5 and 13 and annotated with bioactivities (Table 2). Among the 9.2 million screening compounds available in stock from commercial providers as listed in the ZINC database,24 we similarly found 0.87% (80
153 molecules) 2,2′-diazabiaryls, however with relatively fewer pyrazolyl-pyridines (4673 molecules) but more abundant pyrazolyl-pyrimidones (1960 molecules), which suggests that such 2,2′-diazabiraryl compounds are relatively easy to synthesize and drug-like. The abundance of 2,2′-diazabiaryls was significantly lower in theoretically enumerated chemical space databases such as FDB17,25 GDBMedChem26 and GDBChEMBL,27 reflecting the fact that these databases contain a much higher fraction of non-aromatic molecules.
| Database | ChEMBL24 | ZINCinStock | FDB17 | GDBMedChem | GDBChEMBL |
|---|---|---|---|---|---|
a Compounds were extracted using the corresponding SMARTS patterns. A TMAP of the 15 585 2,2′-diazabiaryls from ChEMBL24 is available at http://tm.gdb.tools/dmt1_inhibitors/diazabiaryls/. See also Fig. 6c–f and Table S1† for SMARTS patterns.
|
|||||
| Number of molecules | 1 820 035 |
9 238 092 |
10 101 204 |
9 994 112 |
9 978 095 |
| 2,2′-Diazabiaryls | 15 585 |
80 153 |
13 880 |
6241 | 2629 |
| % of database | 0.86% | 0.87% | 0.14% | 0.06% | 0.03% |
| Pyrazolyl-pyridines | 1717 | 4673 | 110 | 72 | 65 |
| % of 2,2′-diazabiaryls | 11% | 5.8% | 0.8% | 1.2% | 2.5% |
| Pyrazolyl-pyrimidones | 169 | 1960 | 0 | 26 | 0 |
| % of 2,2′-diazabiaryls | 1.1% | 2.4% | 0.0% | 0.4% | 0.0% |
We recently developed interactive chemical space visualization tools that provide useful insights into polypharmacology and structure–activity relationships.28,29 To obtain a closer insight into the selected 2,2′-diazabiaryls, we created a tree-map (TMAP)30 representing pyrazolyl-pyridines and pyrazolyl-pyrimidones from ChEMBL (Fig. 6c). The TMAP can be interactively analyzed in-browser by loading its layout using the program Faerun.31 In this map each compound is represented by a point that can be color-coded according to selected properties, and the molecular structure is shown upon mouse-over by the program Smilesdrawer.32 The compounds are connected in a tree topology according to their structural similarities as measured by the MHFP6 molecular fingerprint.33 The TMAP illustrates the diversity of pyrazolyl-pyrimidones and pyrazolyl-pyridines in ChEMBL. The pyrazolyl-pyrimidones investigated here belong to the more drug-like molecules in the series, as can be appreciated by color-coding for drug-likeness using the quantitative estimate of drug-likeness (QED).34
The activities reported for these molecules in ChEMBL cover a diversity of target types, spanning from transporters and ion channels to enzymes and uncharacterized targets. We found 45 2,2′-diazabiaryl including 18 pyrazolyl-pyrimidones closely related to our compounds and reported with transporter activity in ChEMBL (Table S2†). The ChEMBL record indicates, among various targets, inhibition of glucose and hexose transporters identified in antiparasitic screening campaigns against Leishmania mexicana, tuberculosis, and Plasmodium falciparum.35 However, the screens performed were not transport assays but simple cytotoxicity tests. Considering that metal chelation is known to be a mechanism of toxicity against these parasites,36,37 one can speculate that the reported transporter inhibitory activity of these pyrazolyl-pyrimidones might in fact reflect toxicity as a result of metal chelation. Further activities reported for pyrazolyl-pyrimidones and pyrazolyl-pyridines might similarly be the consequence of metal chelation rather than the attributed activity, including a recently reported activity of 3 and 5 against cancer cells.38 Note that potentially metal chelating 2,2′-diazabiaryls such as pyrazolyl-pyrimidones are not listed in PAINS (pan-assay interference compounds).39–41 Indeed a systematic checking of our compounds with the PAINS filter from RDkit does not flag any molecule in the series.
Our combined experiments suggest that this non-selective inhibition reflects chelation of Fe2+ in the assay medium by these molecules to form an insoluble precipitate. Although metal chelation is often exploited as a mechanism of drug action,10–12 metal chelation is in principle an assay interference. This effect might occur across thousands of bioactive molecules reported in ChEMBL that share the same potentially chelating substructure and explain some of their reported bioactivities. Pyrazolyl-pyrimidones and related potentially metal chelating 2,2′-diazabiaryls are not considered in the list of pan-assay interference compounds (PAINS), however these should be checked whenever screening SLC targets.
Footnotes |
| † Electronic supplementary information (ESI) available: Methods for bioassays, details of kinetic studies, SMARTS pattern analysis, chemical synthesis procedures and compound characterization, X-ray deposition information, HPLC purity table, NMR spectra of final compounds, and SMILES list of final compounds with activities and of compounds extracted from ChEMBL. CCDC 1976759, 1976832 and 1976862 For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/d0md00085j |
| ‡ These authors contributed equally to the work. |
| This journal is © The Royal Society of Chemistry 2020 |