Literature DB >> 33643826

Development of hedgehog pathway inhibitors by epigenetically targeting GLI through BET bromodomain for the treatment of medulloblastoma.

Xiaohua Liu1,2,3, Yu Zhang4, Yalei Li5, Juan Wang4, Huaqian Ding1,3, Wenjing Huang4, Chunyong Ding1,2, Hongchun Liu5, Wenfu Tan4, Ao Zhang1,2,6,3,7.   

Abstract

Medulloblastoma (MB) is a common yet highly heterogeneous childhood malignant brain tumor, however, clinically effective molecular targeted therapy is lacking. Modulation of hedgehog (HH) signaling by epigenetically targeting the transcriptional factors GLI through bromodomain-containing protein 4 (BRD4) has recently spurred new interest as potential treatment of HH-driven MB. Through screening of current clinical BRD4 inhibitors for their inhibitory potency against glioma-associated oncogene homolog (GLI) protein, the BRD4 inhibitor 2 was selected as the lead for further structural optimization, which led to the identification of compounds 25 and 35 as the high potency HH inhibitors. Mechanism profiling showed that both compounds suppressed HH signaling by interacting with the transcriptional factor GLI, and were equally potent against the clinical resistant mutants and the wild type of smoothened (SMO) receptor with IC50 values around 1 nmol/L. In the resistant MB allograft mice, compound 25 was well tolerated and markedly suppressed tumor growth at both 5 mg/kg (TGI = 83.3%) and 10 mg/kg (TGI = 87.6%) doses. Although further modification is needed to improve the pharmacokinetic (PK) parameters, compound 25 represents an efficacious lead compound of GLI inhibitors, possessing optimal safety and tolerance to fight against HH-driven MB.
© 2021 Chinese Pharmaceutical Association and Institute of Materia Medica, Chinese Academy of Medical Sciences. Production and hosting by Elsevier B.V.

Entities:  

Keywords:  BCC, basal cell carcinoma; BET, bromo and extra C-terminal bromodomain proteins; BRD4; BRD4, bromodomain-containing protein 4; Drug resistance; GLI; HH, hedgehog; HTRF, homogeneous time-resolved fluorescence; Hedgehog signaling pathway; MB, medulloblastoma; Medulloblastoma; PK, pharmacokinetic; PTCH, patched; SAR, structure−activity relationship; SHH, Sonic hedgehog; SMO, smoothened; TGI, tumor growth inhibition; WNT, wingless; hERG, human ether-a-go-go-related gene; i.v., intravenous injection; p.o., per os

Year:  2020        PMID: 33643826      PMCID: PMC7893122          DOI: 10.1016/j.apsb.2020.07.007

Source DB:  PubMed          Journal:  Acta Pharm Sin B        ISSN: 2211-3835            Impact factor:   11.413


Introduction

Medulloblastoma (MB) is one of the most common childhood malignant brain tumors,, accounting for 6.6% of all pediatric central nervous system tumors (or 63.3% of all embryonal tumors) aging less than 19 years. It is a small round blue cell tumor of the cerebellum with high heterogeneity, and biological characterization and molecular classification of MB remain as an unsettled challenge. Historically, MB is divided into four distinct molecular subgroups, namely: wingless (WNT), Sonic hedgehog (SHH), group 3 and group 44. Among which, the SHH contains 30% of all MB diagnosis and is the most studied MB subgroup. The SHH MB is characterized by the somatic mutations in the SHH pathway genes (Fig. 1), including the 12-pass transmembrane receptor patched (PTCH1), G protein-coupled receptor smoothened (SMO), and the downstream transcriptional factors GLI (GLI1, GLI2, and GLI3),. In addition to MB, aberrant activations of SHH have also been observed in many other cancers, such as basal cell carcinoma (BCC). Currently, three SHH signaling pathway inhibitors (vismodegib,, sonidegib,, and glasdegib, Fig. 1) have been approved as SMO antagonists for clinical use of locally advanced or metastatic BCC, or acute leukemia. Unfortunately, clinical trials of these compounds with MB patients have not yet met success due to their transient effects as well as SMO mutation-induced drug resistance during treatment14, 15, 16. Therefore, the standard treatment for MB after surgery remains to be craniospinal irradiation, chemotherapy or their combinations, which are associated with many adverse effects17, 18, 19. Although more precise molecular classification of MB and the genomic and epigenetic alterations have been recently reported, clinically effective treatment of MB remains to be an unmet medical need.
Figure 1

The hedgehog signaling pathway and representative inhibitors of SMO, GLI and BRD4. Ac, acetylation; BRD4, bromodomain-containing protein 4; GLI, glioma-associated oncogene homolog protein; GLIA, activated GLI; HH, hedgehog; PTCH, patched; SMO, smoothened; SUFU, suppressor of fused.

The hedgehog signaling pathway and representative inhibitors of SMO, GLI and BRD4. Ac, acetylation; BRD4, bromodomain-containing protein 4; GLI, glioma-associated oncogene homolog protein; GLIA, activated GLI; HH, hedgehog; PTCH, patched; SMO, smoothened; SUFU, suppressor of fused. Since the transcriptional factor GLI is the final effector of the canonical HH signaling, acting downstream of SMO, inhibition of GLI function may overcome clinically observed drug resistance of SMO antagonist treatment20, 21, 22, 23. Meanwhile, Gli1 is the first identified HH pathway gene amplified in many HH-dependent cancers, and expression of Gli1 is a reliable marker of HH pathway activity,. Therefore, development of GLI-targeting inhibitors has become an attractive strategy to treat MB and many other cancers. Unfortunately, currently available GLI inhibitors such as GANT61, and As2O3 (Fig. 1) suffer from either toxicity or poor druglikeness or both. Recently, modulation of HH signaling by epigenetically targeting GLI through bromodomain-containing protein 4 (BRD4) has spurred new interest in the development of effective treatment for HH-dependent tumors, especially for MB,. BRD4, together with BRD2, BRD3 and the testis-specific BRDT constitutes the important family of bromo and extra C-terminal (BET) bromodomain proteins, all containing two bromodomains BD1 and BD2 at the N terminus,. As a serine kinase of RNA Pol II and an atypical histone acetyltransferase, BRD4 recruits transcriptional regulatory complexes to acetylate chromatin via recognition of acetylated lysines. Global mRNA expression profiling of clinical primary MB and normal cerebellar samples indicate that BRD2, BRD3, and BRD4 are all up-expressed in primary MB, especially SHH MB subgroups, whereas only marginal expression of BRD4 expressed in normal cerebellar tissues. These results indicate that BRD4 together with BRD2 and BRD3 is involved in the pathogenesis of MB. To decode the underlying mechanism of BRD4 involved in MB, Tang et al. and Long et al. recently disclosed that the bromodomains of BRD4 interact with GLI by occupying the cancer-specific proximal promoters of Gli1 and Gli2 genes. Notably, knockdown of Brd4 did not abrogate GLI activity completely, and knockdown of either Brd2 or Brd3 also resulted in substantial GLI reduction, suggesting pan-BET inhibitors might be more efficacious than selective ones to modulate GLI. Collectively, these pioneering studies provide the basis that inhibition of GLI can be achieved by suppression of BRD4, thus providing an alternative strategy for development of GLI inhibitors to treat HH-dependent MB. This strategy is especially appealing since development of highly potent GLI inhibitors is challenging due to the lack of targetable binding domains of GLI. Therefore, in this report we first conducted a drug repurposing campaign by screening current clinically investigational BRD4 inhibitors bearing different chemotypes for their inhibitory potency against the HH pathway. The BRD4 inhibitor ABBV-07534, 35, 36, 37, 38 is then selected as the lead for further structural optimization leading to identification of compounds 25 and 35 as highly potent GLI inhibitors, especially compound 25 which was well tolerated and showed significant tumor growth inhibition in the resistant MB allograft mice. These results not only confirm the effectiveness of development of GLI inhibitors through optimization of BRD4 inhibitors, but also provide an alternative utility for BRD4 inhibitors in the treatment of many highly heterogeneous and untreatable pediatric brain tumors, such as MB,,.

Results and discussions

Chemistry

The synthesis of carbamates 5 and 6 is outlined in Scheme 1. Pd-catalyzed coupling of 3-amino-5-bromo-1-methylpyridin-2(1H)-one (46) with borate 47 provided compound 48 in 84% yield, which was then subjected to aromatic substitution to give 2,4-difluorophenyl ether 49. Substitution of 49 with ethyl carbonochloridate or 2-methoxyethyl carbonochloridate afforded corresponding carbamate 50a or 50b in 66%–71% yields. Subsequent reduction of 50a or 50b under iron powder followed by substitution with ethanesulfonyl chloride gave target compounds 5 and 6 in 54% and 70% yields, respectively.
Scheme 1

Synthesis of compounds 5 and 6. Reagents and conditions: (a) 47, Pd(PPh3)4, 1,4-dioxane/2 mol/L Na2CO3 aq. (5:1), microwave 120 °C, 1 h, 84%; (b) 2,4-difluorophenol, Cs2CO3, DMF, r.t., 1 h; (c) ethyl carbonochloridate or 2-methoxyethyl carbonochloridate, NaH, THF, 0 °C to rt, 1 h, 66%–71%; (d) Fe, NH4Cl, EtOH/H2O (5:2), 70 °C, 3 h; (e) ethanesulfonyl chloride, pyridine, microwave 130 °C, 15 min, 54%–70%.

Synthesis of compounds 5 and 6. Reagents and conditions: (a) 47, Pd(PPh3)4, 1,4-dioxane/2 mol/L Na2CO3 aq. (5:1), microwave 120 °C, 1 h, 84%; (b) 2,4-difluorophenol, Cs2CO3, DMF, r.t., 1 h; (c) ethyl carbonochloridate or 2-methoxyethyl carbonochloridate, NaH, THF, 0 °C to rt, 1 h, 66%–71%; (d) Fe, NH4Cl, EtOH/H2O (5:2), 70 °C, 3 h; (e) ethanesulfonyl chloride, pyridine, microwave 130 °C, 15 min, 54%–70%. The preparation of 3-aminopyridinones 7–11 was described in Scheme 2. Compounds 52a−e were obtained from substitution reaction of aminopyridinone 46 or C–N coupling reaction of 51 with appropriate amines. Coupling of 52a−e with borate 53 under Pd(PPh3)4 catalysis yielded compounds 54a−e in 59%–87% yields, which were readily converted to sulfonamides 7–11 by reacting with ethanesulfonyl chloride in the presence of Et3N.
Scheme 2

Synthesis of compounds 7–11. Reagents and conditions: (a) RH, Et3N, DCM, r.t., 1 h; (b) RH, Pd2(dba)3, xantphos, Cs2CO3, toluene, 100 °C, 3 h, 19%–71%; (c) 10, Pd(PPh3)4, 1,4-dioxane/2 mol/L Na2CO3 aq. (5:1), 100 °C, 1 h, 59%–87%; (d) i) ethanesulfonyl chloride, Et3N, DCM, r.t., 1 h; ii) NaOH, 1,4-dioxane, 90 °C, 30 min, 33%–62%.

Synthesis of compounds 7–11. Reagents and conditions: (a) RH, Et3N, DCM, r.t., 1 h; (b) RH, Pd2(dba)3, xantphos, Cs2CO3, toluene, 100 °C, 3 h, 19%–71%; (c) 10, Pd(PPh3)4, 1,4-dioxane/2 mol/L Na2CO3 aq. (5:1), 100 °C, 1 h, 59%–87%; (d) i) ethanesulfonyl chloride, Et3N, DCM, r.t., 1 h; ii) NaOH, 1,4-dioxane, 90 °C, 30 min, 33%–62%. As shown in Scheme 3, substituted arylamine 55, which was prepared by following literature procedures, was reacted with sulfurisocyanatidic chloride followed by cyclization with 1,2-dibromoethane or 1,3-dibromopropane to deliver 56a and 56b in 90% and 88% yields, respectively. Subsequent reduction of the esteric moiety by NaBH4 followed by N-Ts deprotection with NaOH afforded compounds 12 and 13 in 72% and 78% yields, respectively.
Scheme 3

Synthesis of compounds 12 and 13. Reagents and conditions: (a) methyl 2-(chlorosulfonyl)acetate, Et3N, DCM, r.t., 1 h, 88%; (b) 1,2-dibromoethane or 1-bromo-3-chloropropane, K2CO3, DMF, 70 °C, 88%–90%; (c) NaBH4, CaCl2, dry THF, 0 °C to r.t., 3 h; (d) NaOH, 1,4-dioxane, 90 °C, 30 min, 72%–78% over two steps.

Synthesis of compounds 12 and 13. Reagents and conditions: (a) methyl 2-(chlorosulfonyl)acetate, Et3N, DCM, r.t., 1 h, 88%; (b) 1,2-dibromoethane or 1-bromo-3-chloropropane, K2CO3, DMF, 70 °C, 88%–90%; (c) NaBH4, CaCl2, dry THF, 0 °C to r.t., 3 h; (d) NaOH, 1,4-dioxane, 90 °C, 30 min, 72%–78% over two steps. The synthesis of sulfuric diamides 14–23 is illustrated in Scheme 4. Sulfuric diamide 57 was generated from substitution reaction of 55 with freshly made tert-butyl(chlorosulfonyl) carbamate in 74% yield. Cyclization of compound 57 with 1,2-dibromoethane or 1,3-dibromopropane delivered 1,2,5-thiadiazolidine 1,1-dioxides 58a and 58b in 54% and 75% yields, respectively. N-Boc deprotection of 58a or 58b with TFA and further substitution with different alkyl bromides afforded compounds 59a−d in 68%–98% yields. Subsequent N-Ts deprotection of compounds 59a−c in the presence of NaOH afforded corresponding compounds 14, 15 and 23 in 89%–93% overall yields. It should be noted that compounds 16–19 were prepared by N-Boc deprotection of 58b with TFA, subsequent N-Ts deprotection and further substitution with 1,1-difluoro-2-iodoethane or different alkyl bromides in 28%–55% overall yields. In addition, substitution of bromide 59d with different amines followed by N-Ts deprotection with NaOH delivered compounds 20–22 in 75%–78% yields.
Scheme 4

Synthesis of compounds 14–23. Reagents and conditions: (a) sulfurisocyanatidic chloride, t-BuOH, Et3N, DCM, 0 °C to r.t., 2.5 h, 74%; (b) 1,2-dibromoethane or 1,3-dibromopropane, K2CO3, MeCN, reflux, overnight, 54%–75%; (c) TFA, DCM, r.t., 2 h; (d) RBr or (for 16) RI, K2CO3, DMF, 60 °C, 1–5 h; (e) NaOH, 1,4-dioxane, 90 °C, 30 min; (f) from 59h, amines, K2CO3, MeCN, 60 °C, 75%–78%.

Synthesis of compounds 14–23. Reagents and conditions: (a) sulfurisocyanatidic chloride, t-BuOH, Et3N, DCM, 0 °C to r.t., 2.5 h, 74%; (b) 1,2-dibromoethane or 1,3-dibromopropane, K2CO3, MeCN, reflux, overnight, 54%–75%; (c) TFA, DCM, r.t., 2 h; (d) RBr or (for 16) RI, K2CO3, DMF, 60 °C, 1–5 h; (e) NaOH, 1,4-dioxane, 90 °C, 30 min; (f) from 59h, amines, K2CO3, MeCN, 60 °C, 75%–78%. As shown in Scheme 5, compounds 24 and 25 bearing the 3-fluoro-6-methyl-1,6-dihydro-7H-pyrrolo[2,3-c]pyridin-7-one scaffold were prepared. Starting from 4-fluoro-1H-pyrrole-2-carboxylic acid 60, C-3 fluorinated bicyclic pyrrolopyridone core 62 was prepared in 68% overall yield by first condensation of 60 with 2,2-dimethoxyethanamine followed by self-cyclization in the presence of p-toluene sulfonic acid. Protection of pyrrolidine 62 with TsCl in the presence of NaH followed by bromination with NBS generated the intermediate 63 in 84% yield over two steps. Coupling of aryl bromide 63 with borate 53 under Pd(PPh3)4 catalysis afforded compound 64 in 85% yield. Treatment of compound 64 with 3-chloropropane-1-sulfonyl chloride or 4-chlorobutane-1-sulfonyl chloride yielded bis-substituted 65a and 65b in 64% and 68% yields, respectively. Subsequent N-Ts deprotection, N-sulfonyl removal and intramolecular cyclization of compounds 65a and 65b in one pot with NaOH afforded target compounds 24 and 25 in 95% and 91% yields, respectively.
Scheme 5

Synthesis of compounds 24 and 25. Reagents and conditions: (a) 2,2-dimethoxyethanamine, DIPEA, propanephosphonic acid cyclic anhydride, THF, −5–60 °C, 20 h, 75%; (b) TsOH·H2O, THF, 60 °C, overnight, 90%; (c) TsCl, NaH, DMF, 0 °C, 1 h; (d) NBS, TsOH·H2O, THF, r.t., 1 h, 84% over 2 steps; (e) 53, Pd(PPh3)4, 1,4-dixone/2 mol/L Na2CO3 aq. (5:1), 100 °C, 1 h, 85%; (f) 3-chloropropane-1-sulfonyl chloride or 4-chlorobutane-1-sulfonyl chloride, Et3N, DCM, r.t., 1 h, 64%–68%; (g) NaOH, 1,4-dioxane, 90 °C, 0.5 h, 91%–95%.

Synthesis of compounds 24 and 25. Reagents and conditions: (a) 2,2-dimethoxyethanamine, DIPEA, propanephosphonic acid cyclic anhydride, THF, −5–60 °C, 20 h, 75%; (b) TsOH·H2O, THF, 60 °C, overnight, 90%; (c) TsCl, NaH, DMF, 0 °C, 1 h; (d) NBS, TsOH·H2O, THF, r.t., 1 h, 84% over 2 steps; (e) 53, Pd(PPh3)4, 1,4-dixone/2 mol/L Na2CO3 aq. (5:1), 100 °C, 1 h, 85%; (f) 3-chloropropane-1-sulfonyl chloride or 4-chlorobutane-1-sulfonyl chloride, Et3N, DCM, r.t., 1 h, 64%–68%; (g) NaOH, 1,4-dioxane, 90 °C, 0.5 h, 91%–95%. As shown in Scheme 6, carboxylic acids 26 and 27 and amides 28–32 were obtained using 6-methyl-1-tosyl-1,6-dihydro-7H-pyrrolo[2,3-c]pyridin-7-one 66 as the starting material. Substitution of the pyrrole C-2 of 66 with ethyl chloroformate delivered ethyl ester 67, which was then converted to compounds 68a and 68b through a similar reaction as that for 65a and 65b. Subsequent N-Ts deprotection, N-sulfonyl removal, intramolecular cyclization and ethyl ester hydrolysis were conducted in one pot in the presence of NaOH leading to acids 26 and 27 in 83% and 89% yields, respectively. Condensation of 26 or 27 with NH4Cl, CH3CH2NH2 or 3-(4-methylpiperazin-1-yl)propan-1-amine in the presence of HOBT, EDCI and Et3N produced compunds 28–32 in 65%–78% yields.
Scheme 6

Synthesis of compounds 26–32. Reagents and conditions: (a) n-BuLi, ethyl carbonochloridate, dry THF, N2, −10 °C, 1 h, 60%; (b) NBS, TsOH·H2O, THF, r.t., 1 h; (c) 53, Pd(PPh3)4, 1,4-dixone/2 mol/L Na2CO3 aq. (5:1), 100 °C, 1 h, 76%; (d) 3-chloropropane-1-sulfonyl chloride or 4-chlorobutane-1-sulfonyl chloride, Et3N, DCM, r.t., 1 h, 84%–88%; (e) NaOH, 1,4-dioxane, 90 °C, 0.5 h, 89%–92%; (f) HOBT, EDCI, Et3N, NH4Cl or CH3CH2NH2 or 3-(4-methylpiperazin-1-yl)propan-1-amine, DCM, r.t., 3 h, 65%–78%.

Synthesis of compounds 26–32. Reagents and conditions: (a) n-BuLi, ethyl carbonochloridate, dry THF, N2, −10 °C, 1 h, 60%; (b) NBS, TsOH·H2O, THF, r.t., 1 h; (c) 53, Pd(PPh3)4, 1,4-dixone/2 mol/L Na2CO3 aq. (5:1), 100 °C, 1 h, 76%; (d) 3-chloropropane-1-sulfonyl chloride or 4-chlorobutane-1-sulfonyl chloride, Et3N, DCM, r.t., 1 h, 84%–88%; (e) NaOH, 1,4-dioxane, 90 °C, 0.5 h, 89%–92%; (f) HOBT, EDCI, Et3N, NH4Cl or CH3CH2NH2 or 3-(4-methylpiperazin-1-yl)propan-1-amine, DCM, r.t., 3 h, 65%–78%. Compounds 33–39 bearing benzofuran or cycloalkyl amine moieties were synthesized as described in Scheme 7. Aromatic substitution of compound 69 with benzofuran-5-ol or different amines furnished compounds 70a−e in 33%–98% yields, which were then reduced by hydrogen atmosphere under Pd/C or by iron powder in the presence of NH4Cl to afford corresponding anilines 71a−e. It is worth noting that N-Ts protecting group was partially taken off during substitution reaction with trans-4-methylcyclohexyl amine. In this case, additional N-Ts protection under NaH was necessary to afford compound 70b. Treatment of compounds 71a−e with 3-chloropropane-1-sulfonyl chloride or 4-chlorobutane-1-sulfonyl chloride provided di-substituted intermediate 72a and mono-substituted 72b−g in 34%–78% yields. It should be pointed out that for compound 71e, the two hydroxyl groups need to be protected with TMSCl before treatment with 4-chlorobutane-1-sulfonyl chloride. Finally, N-deprotection, N-sulfonyl removal and intramolecular cyclization of compounds 72a−g were conducted in one pot in the presence of NaOH to provide target compounds 33–39 in 85%–93% overall yields.
Scheme 7

Synthesis of compounds 33–39. Reagents and conditions: (a) RH, DIPEA, DMSO, 80 °C, 6 h, 33%–51%; or (for 70a) RH, K2CO3, DMF, r.t., 30 min, 98%; or (for 70b): i) RH, DMSO, 80 °C, 6 h; ii) TsCl, NaH, DMF, 0 °C, 1 h, 72% over two steps. (b) Pd/C, H2, 45 °C, CHCl3/MeOH (1:1), 8 h; or (for 71a) Fe, NH4Cl, EtOH/H2O (5:2), 70 °C, 1 h. (c) 3-Chloropropane-1-sulfonyl chloride or 4-chlorobutane-1-sulfonyl chloride, Et3N, DCM, r.t., 1 h, 77%–88%; or (for 72g): i) TMSCl, Et3N, DMF, 0 °C, 0.5 h; ii) 4-chlorobutane-1-sulfonyl chloride, Et3N, DCM, r.t., 1 h; (d) NaOH, 1,4-dioxane, 90 °C, 0.5 h, 85%–93%.

Synthesis of compounds 33–39. Reagents and conditions: (a) RH, DIPEA, DMSO, 80 °C, 6 h, 33%–51%; or (for 70a) RH, K2CO3, DMF, r.t., 30 min, 98%; or (for 70b): i) RH, DMSO, 80 °C, 6 h; ii) TsCl, NaH, DMF, 0 °C, 1 h, 72% over two steps. (b) Pd/C, H2, 45 °C, CHCl3/MeOH (1:1), 8 h; or (for 71a) Fe, NH4Cl, EtOH/H2O (5:2), 70 °C, 1 h. (c) 3-Chloropropane-1-sulfonyl chloride or 4-chlorobutane-1-sulfonyl chloride, Et3N, DCM, r.t., 1 h, 77%–88%; or (for 72g): i) TMSCl, Et3N, DMF, 0 °C, 0.5 h; ii) 4-chlorobutane-1-sulfonyl chloride, Et3N, DCM, r.t., 1 h; (d) NaOH, 1,4-dioxane, 90 °C, 0.5 h, 85%–93%. The synthesis of C-2 or C-3 fluorinated pyrropyridinones 40–45 bearing a 4-methylcyclohexan-1-amino motif was described in Scheme 8. Fluorination of pyrrole C-2 of 66 with NFSI under n-BuLi delivered compound 73 in 45% yield, which was brominated to yield intermediate 74 through a similar reaction as described in Scheme 5. Borate 77 was obtained in 80% overall yield with 1-fluoro-4-nitrobenzene (75) as the starting material by substitution with trans-4-methylcyclohexyl amine followed by coupling with bis(pinacolato)diboron. The coupling of borate 77 with aryl bromide 63 or 74 delivered corresponding 78a and 78b in 57% and 68% yields, respectively, which were then converted to target compounds 40–45 in 65%–85% overall yields through a similar reaction procedure as described in Scheme 7.
Scheme 8

Synthesis of compounds 40–45. Reagents and conditions: (a) n-BuLi, NFSI, dry THF, N2, −20 °C, 2.5 h, 45%; (b) NBS, TsOH·H2O, THF, r.t., 1 h; (c) trans-4-methylcyclohexyl amine, DMSO, 100 °C, 2 h; (d) bis(pinacolato)diboron, Pd(dppf)Cl2, AcOK, 1,4-dioxane, N2, 100, 5 h, 80% over two steps; (e) 63 or 73, Pd2(dba)3, meCgPPh, K3PO4, dioxane/H2O (4:1), 80 °C, 5 h, 57%–68%; (f) Pd/C, H2, 45 °C, CHCl3/MeOH (1:1), 8 h; (g) ethanesulfonyl chloride or 3-chloropropane-1-sulfonyl chloride or 4-chlorobutane-1-sulfonyl chloride, Et3N, DCM, r.t., 1 h, 63%–82%; (h) NaOH, 1,4-dioxane, 90 °C, 0.5 h, 65%–85%.

Synthesis of compounds 40–45. Reagents and conditions: (a) n-BuLi, NFSI, dry THF, N2, −20 °C, 2.5 h, 45%; (b) NBS, TsOH·H2O, THF, r.t., 1 h; (c) trans-4-methylcyclohexyl amine, DMSO, 100 °C, 2 h; (d) bis(pinacolato)diboron, Pd(dppf)Cl2, AcOK, 1,4-dioxane, N2, 100, 5 h, 80% over two steps; (e) 63 or 73, Pd2(dba)3, meCgPPh, K3PO4, dioxane/H2O (4:1), 80 °C, 5 h, 57%–68%; (f) Pd/C, H2, 45 °C, CHCl3/MeOH (1:1), 8 h; (g) ethanesulfonyl chloride or 3-chloropropane-1-sulfonyl chloride or 4-chlorobutane-1-sulfonyl chloride, Et3N, DCM, r.t., 1 h, 63%–82%; (h) NaOH, 1,4-dioxane, 90 °C, 0.5 h, 65%–85%.

Screening of clinical BRD4 inhibitors against HH signaling pathway

To gain insights on the activity relationship between BRD4 and GLI, we selected seven clinically studied BRD4 inhibitors and tested their inhibitory activity in the HH signaling pathway in our dual luciferase reporter assays using light II cells. These cells are NIH-3T3 cells stably transfected with a GLI-responsive firefly luciferase reporter and renilla-luciferase expression vector,. As shown in Table 1, the seven clinically studied BRD4 inhibitors were found to show high but quite different potency against the HH signaling pathway with IC50 values ranging between 0.15 and 40.6 nmol/L. Among them, the phase II BRD4 inhibitor 1 (BMS-986158) is the most potent, whereas INCB-057643 is the least potent HH pathway inhibitor, with IC50 values of 0.15 and 40.6 nmol/L, respectively. The overall trend of activity is in agreement with the BRD4 activity where the phase II compound 1 remains the most potent BRD4 inhibitor, and compound 4 (AZD5153) is much less potent, with IC50 values of 0.16 and 18.2 nmol/L, respectively. This result indicates that development of high potency GLI-targeting HH pathway inhibitors might be realized by optimization of BRD4 activity through structural modification, since direct structural information on the transcriptional factor GLI is limited. Considering the potential toxicity liability and pharmacokinetic (PK) disadvantage of the tricyclic compounds (e.g., 1 and 3, Table 1),, we decided to choose compound 2 (ABBV-075, Table 1) as our lead HH inhibitor for structural optimization. This compound shows high inhibitory activity (0.41 nmol/L) against HH pathway, but is nearly 10-fold less potent against BRD4 (3.4 nmol/L). Our objective is to identify a high potency and well tolerated HH pathway inhibitor for treatment of MB, especially those resistant to current SMO-targeting HH inhibitors.
Table 1

Inhibition of clinical BRD4 inhibitors against BRD4 BD1 and the HH pathwaya.

Compd.Clinical stageStructureIC50 (nmol/L)
BRD4 (BD1)GLI-luc
BMS-986158 (1)I/IIImage 2<0.160.15 ± 0.05
ABBV-075 (2)IImage 33.4 ± 0.20.41 ± 0.34
PLX51107IImage 40.49 ± 0.145.5 ± 1.3
OTX-015IIImage 57.7 ± 0.210.64 ± 0.02
JQ1 (3)DiscoveryImage 69.8 ± 2.22.6 ± 1.7
INCB-057643I/IIImage 711.6 ± 1.840.6 ± 4.9
AZD5153 (4)IImage 818.2 ± 2.732.1 ± 7.71

IC50 values are shown as the mean ± SD (nmol/L) from three separate experiments.

Inhibition of clinical BRD4 inhibitors against BRD4 BD1 and the HH pathwaya. IC50 values are shown as the mean ± SD (nmol/L) from three separate experiments.

Structure-based drug design

4-Aryl-1,6-dihydro-7H-pyrrolo[2,3-c]pyridin-7-one 2 (ABBV-075, Fig. 2), developed by AbbVie, is a highly potent non-selective BRD4 inhibitor undergoing phase I clinical trials for the treatment of patients with acute myeloid leukemia, multiple myeloma and advanced solid tumors (NCT02391480). According to the reported co-crystal structure of its close analogue with BRD4 (BD2), the pyrrolopyridone core is located in the KAc binding pocket and provides critical interactions with Asn433 residue. The meta-difluorophenoxyl motif fills the WPF hydrophobic shelf (a hydrophobic region constructed by residues W81, P82 and F83), and the ethyl sulfonamide chain occupies the ZA channel of the BRD4 (BD2) protein. In addition, the pyrrole fragment and the ethyl sulfonamide chain are both exposed to the solvent interaction region, which may tolerate various substituents. Inspired by this analysis and in view of the lack of targetable binding domains of the transcriptional factors GLI, we have designed several series of analogues by modifying these interaction sites through three approaches, including: (1) opening the pyrropyridinone bicyclic core to afford flexible 3-aminopyridiones; (2) cyclizing the ethanesulfonamide motif; (3) replacing the hydrophobic difluorophenyl with a water soluble moiety (Fig. 1).
Figure 2

Structural analysis and design of new HH pathway inhibitors.

Structural analysis and design of new HH pathway inhibitors.

Structure−activity relationship (SAR) study

To investigate how the structural alterations impact activities of GLI, we first cleaved the bicyclic pyrrolopyridone core of compound 2 in 3-aminopyridinones 5–11. As shown in Table 2, all these compounds are much less potent than compound 2. Carbamate 5 has an IC50 value of 27.6 nmol/L against BRD4, which is 8-fold less potent than that of 2. Both carbamate 6 with a longer chain and the reverse amide 7 show more reduction of potency. 3-Ethanesulfonamido pyridinone 8 and 3-((1H-imidazol-2-yl)amino) pyridinone 9 display modest potency against BRD4 with IC50 value of 42.8 and 47.6 nmol/L, respectively. N-Heterocycle-substituted 3-aminopyridinones 10 and 11 show much reduced potency with IC50 values greater than 100 nmol/L. These results indicate that although the substituted 3-amino analogues may maintain the H-bonding with Asn433 residue in the KAc binding pocket of BRD4 (BD1), the electronic or steric property potentially affects the interaction with BRD4 as well. Disappointingly, all the compounds show markedly reduced potency against the HH pathway with IC50 values all greater than 100 nmol/L in the GLI-luciferase assay.
Table 2

Inhibition of compounds against BRD4 (BD1) and the HH pathwaya.

Compd.R1IC50 (nmol/L)
BRD4 (BD1)GLI-luc reporter
5Image 1027.6 ± 5.7157 ± 1.3
6Image 11>100>100
7Image 12>100>100
8Image 1342.8 ± 9.6>100
9Image 1447.6 ± 2.9>100
10Image 15>100>100
11Image 16>100>100
ABBV-075 (2)3.4 ± 0.20.41 ± 0.34

IC50 values are shown as the mean ± SD (nmol/L) from three separate experiments. –Not applicable.

Inhibition of compounds against BRD4 (BD1) and the HH pathwaya. IC50 values are shown as the mean ± SD (nmol/L) from three separate experiments. –Not applicable. Next, considering that there is a large space in the ZA channel interaction region in the binding mode of 2 with BRD4 (PDB code: 5uvx), we evaluated compounds derived from cyclization of the ethanesulfonamide. As shown in Table 3, the five- and six-membered cyclic sulfonamides 12 and 13 display much higher potency against BRD4 with IC50 values of 0.25 and 0.60 nmol/L, respectively. However, their potency against the HH pathway is significantly reduced with IC50 values of 3.29 and 20.5 nmol/L, respectively. In order to increase aqueous solubility (Table 7), we incorporated an additional N-atom moiety with various substituents leading to compounds 14–23. Most of these N-substituted cyclic aminosulfonamides except 16 and 22 were found much more potent than 2, by showing low- or sub-nanomolar potency against BRD4. However, all compounds are significantly less potent than 2 against the HH pathway with IC50 values between 17.7 and 100 nmol/L. Compound 22 bearing a substituted pyrrolidinylethyl moiety has the least potency in both assays. The five-membered analogues are slightly more potent than the six-membered congeners, especially against BRD4. These results indicate that the ZA channel of BRD4 BD1 can tolerate various substitutions, whereas larger substituents are not beneficial for interaction with the HH pathway.
Table 3

Inhibition of compounds against BRD4 (BD1) and the HH pathwaya.

Compd.XnR2IC50 (nmol/L)
BRD4 (BD1)GLI-luc reporter
12C1CH2OH0.25 ± 0.053.29 ± 0.14
13C2CH2OH0.60 ± 0.0720.5 ± 7.74
14N1Image 180.20 ± 0.0326.1 ± 5.17
15N1Image 191.19 ± 0.1071.9 ± 5.43
16N2Image 209.39 ± 0.8028.9 ± 11.9
17N2Image 210.65 ± 0.3633.9 ± 11.1
18N2Image 221.52 ± 0.7256.1 ± 8.25
19N2Image 231.47 ± 0.0317.7 ± 11.0
20N2Image 241.01 ± 0.3442.8 ± 30.7
21N2Image 250.99 ± 0.1939.4 ± 4.76
22N2Image 2671.9 ± 5.43>100
23N2Image 270.58 ± 0.1045.8 ± 19.2
23.4 ± 0.20.41 ± 0.34

IC50 values are shown as the mean ± SD (nmol/L) from three separate experiments. –Not applicable.

Table 7

hERG inhibition and predicted physico-chemical parameters of representative compoundsa.

Compd.GLI-luc reporter IC50 (nmol/L)BRD4 IC50 (nmol/L)hERG IC50 (μmol/L)tPSAb (Å)LogSb (μg/mL)LogPappb (cm/s)VDb (L/kg)
430.062 ± 0.03914.28 ± 1.2787.217.88−5.228−0.24
300.05 ± 0.010.60 ± 0.416.411321.35−5.135−0.795
340.12 ± 0.10.52 ± 0.03>4087.217.60−5.194−0.034
350.33 ± 0.130.99 ± 0.01>4087.217.94−5.235−0.078
250.48 ± 0.081.51 ± 0.4327.384.49.14−4.79−0.551
240.93 ± 0.160.50 ± 0.06>4084.418.14−4.773−0.518
310.81 ± 0.240.13 ± 0.00511.711330.16−5.164−0.799
331.03 ± 0.282.85 ± 0.5097.59.35−5.312−0.345
281.5 ± 0.390.09 ± 0.04>4012751.80−5.189−0.754
20.41 ± 0.343.4 ± 0.235.593.219.60−5.01−0.547

IC50 values are shown as the mean ± SD (nmol/L) from three separate experiments.

Calculated through ADMET lab software viahttp://admet.scbdd.com; tPSA, topological polar surface area; logS, solubility; Papp, Caco-2 permeability; VD, volume distribution. –Not applicable.

Inhibition of compounds against BRD4 (BD1) and the HH pathwaya. IC50 values are shown as the mean ± SD (nmol/L) from three separate experiments. –Not applicable. Inhibition of compounds against BRD4 (BD1) and the HH pathwaya. IC50 values are shown as the mean ± SD (nmol/L) from three separate experiments. –Not applicable. Inhibition of compounds against BRD4 (BD1) and the HH pathwaya. IC50 values are shown as the mean ± SD (nmol/L) from three separate experiments. –Not applicable. Inhibition of compounds against BRD4 (BD1) and the HH pathwaya. IC50 values are shown as the mean ± SD (nmol/L) from three separate experiments. –Not applicable. hERG inhibition and predicted physico-chemical parameters of representative compoundsa. IC50 values are shown as the mean ± SD (nmol/L) from three separate experiments. Calculated through ADMET lab software viahttp://admet.scbdd.com; tPSA, topological polar surface area; logS, solubility; Papp, Caco-2 permeability; VD, volume distribution. –Not applicable. As shown in Table 4, our next effort is to incorporate fluoro or other electron-withdrawing substituent on the C2- or C3-position of the pyrrolopyridone skeleton to the cyclic analogues of ethanesulfonamide motif. Compounds 24 and 25, both with a fluoro-substituent at the C3-position of the pyrrole component, display high potency against BRD4 and the HH pathway with IC50 values ranging between 0.48 and 1.5 nmol/L. However, compounds 26 and 27 bearing a carboxylic group at the C2-position of the pyrrole component show >4-fold decrease of potency against BRD4, together with markedly reduced potency (>100 nmol/L) against the HH pathway. Pleasantly, high potency was observed for compounds 28–32 containing an amide moiety at the C2-position of the pyrrole component. The unsubstituted carboxamides 28 and 29 show an identical IC50 value of 0.09 nmol/L against BRD4, whereas their potency against the HH pathway is slightly lower by showing IC50 values of 1.5 and 2.46 nmol/L, respectively. N-Ethyl substituted amide analogues 30 and 31 also display high potency, especially compound 30 having IC50 values of 0.6 and 0.05 nmol/L, respectively against BRD4 and the HH pathway, which are much more potent than compound 2. However, amide 32 with a longer N-substituent shows much reduced potency (>100 nmol/L) against the HH pathway, despite of the high potency for BRD4.
Table 4

Inhibition of compounds against BRD4 (BD1) and the HH pathwaya.

Compd.nR3R4IC50 (nmol/L)
BRD4 (BD1)GLI-luc reporter
241HF0.50 ± 0.060.93 ± 0.16
252HF1.51 ± 0.430.48 ± 0.08
261COOHH21.0 ± 11.0>100
272COOHH14.3 ± 4.52>100
281CONH2H0.09 ± 0.041.5 ± 0.39
292CONH2H0.09 ± 0.082.46 ± 0.3
301CONHCH2CH3H0.60 ± 0.410.05 ± 0.01
312CONHCH2CH3H0.13 ± 0.0050.81 ± 0.24
321Image 29H0.75 ± 0.07>100
23.4 ± 0.20.41 ± 0.34

IC50 values are shown as the mean ± SD (nmol/L) from three separate experiments. –Not applicable.

Meanwhile, structural modification of the WPF interaction region of BRD4 (BD1) was conducted based on the ethanesulfonamide motif cyclized analogues. As shown in Table 5, compound 33 containing a benzofuran moiety to replace the original 2,4-difluorophenyl group maintains good potency against BRD4 (2.85 nmol/L), whereas the potency for the HH pathway was reduced (1.03 nmol/L). Compounds 34–39 all contain a water soluble substituent to replace the original phenylether motif. Compared to 2, these compounds generally retain good potency against BRD4, but the activity against the HH pathway is distinctly different. Notably, compounds 34 and 35 both bearing a 4-methylcyclohexyl amino substituent show extremely high potency in both assays, although the five-membered analogue 35 is slightly less potent than the six-membered analogue 34. Indeed, compound 34 is the most potent analogue among this series and shows IC50 values of 0.52 and 0.12 nmol/L, respectively against BRD4 and the HH pathway, and is 6.5- and 3.4-fold more potent than compound 2. Incorporation of an additional C4-hydroxyl group on the cyclohexyl (compounds 36 and 37) or reversing the chirality of the resulting quaternary carbon (compound 38) led to no improvement with IC50 values ranging between 1.2 and 3.4 nmol/L for both BRD4 and the HH pathway. In addition, compound 39 bearing a (3-hydroxy-4-hydroxymethyl) cyclopentylamino substituent lost the activity against the HH pathway significantly, despite remaining high potency for BRD4.
Table 5

Inhibition of compounds against BRD4 (BD1) and the HH pathwaya.

Compd.nR5IC50 (nmol/L)
BRD4 (BD1)GLI-luc reporter
331Image 312.85 ± 0.501.03 ± 0.28
341Image 320.52 ± 0.030.12 ± 0.1
352Image 330.99 ± 0.010.33 ± 0.13
361Image 341.24 ± 0.091.58 ± 0.66
372Image 352.60 ± 0561.47 ± 0.66
382Image 363.22 ± 0.111.34 ± 0.48
392Image 373.42 ± 0.25>100
23.4 ± 0.20.41 ± 0.34

IC50 values are shown as the mean ± SD (nmol/L) from three separate experiments. –Not applicable.

Based on the optimization results, we collected all the advantageous elements and incorporated them into a small series of multisubstituted pyrropyridinone analogues bearing a 4-methylcyclohexan-1-amino motif in the WPF binding site. As shown in Table 6, compared to the high potency of compound 34, all these compounds show reduced potency against BRD4 with IC50 values ranging between 9.51 and 20.54 nmol/L. The potency for the HH pathway is also decreased but distinctly different. Similarly modest potency was observed for compounds 40 and 41 in both assays indicating that there is no significant difference for the fluoro substituent on either the C2- or C3-position of the pyrrole component. Dramatically, compounds 42–45 with the ethanesulfonamido component cyclized to five- or six-membered sulfonamides show higher potency for the HH pathway, though their potency against BRD4 is modest, suggesting that the interaction mode for the HH signaling is different from that for BRD4. Among these, compound 43 is the most potent HH pathway inhibitor with IC50 values of 0.062 nmol/L, and shows 230-fold selectivity against BRD4.
Table 6

Inhibition of compounds against BRD4 (BD1) and the HH pathwaya.

Compd.R2R3R4IC50 (nmol/L)
BRD4 (BD1)GLI-luc reporter
40Image 39HF14.25 ± 2.762.06 ± 1.16
41Image 40FH17.57 ± 0.555.88 ± 2.72
42Image 41HF9.51 ± 1.130.33 ± 0.056
43Image 42FH14.28 ± 1.270.062 ± 0.039
44Image 43HF19.10 ± 3.550.58 ± 0.35
45Image 44FH20.52 ± 1.181.75 ± 0.63
23.4 ± 0.20.41 ± 0.34

IC50 values are shown as the mean ± SD (nmol/L) from three separate experiments. –Not applicable.

hERG inhibition and predicted physico-chemical properties of potent compounds

From the SAR, we concluded that the inhibitory activity against BRD4 does not completely parallel with the potency of the HH pathway. Overall, BRD4 is well tolerant to various substitutions, whereas the activity for the HH pathway is much different and highly dependent on different substitution patterns. As listed in Table 7, Table 9 new compounds turn out as the high potency HH inhibitors with IC50 values less than 1.5 nmol/L. The most potent HH inhibitor 43 has an IC50 value of 0.062 nmol/L and is 230-fold more potent than against BRD4 (14.28 nmol/L), whereas compound 28, though highly potent against BRD4 (0.09 nmol/L), is 16-fold less potent for the HH pathway (1.5 nmol/L).
Table 9

Antiproliferative inhibition of various cancer cells.

GI50a (μmol/L)H3132N87HT1080Colo680NHSC-2KYSE450TOV-21GCAL27SKM-1H1975DOHH-2MOLM-14MV4-11TMD8REC-1CHL
250.140.131.02.30.771.60.0230.110.0190.0340.00660.00810.00350.00930.015>10
350.330.460.884.00.321.40.0080.120.00630.00540.00210.0020.00130.00220.00221.1

Antiproliferative effects were initially tested in 10 μmol/L, and then diluted in eight diluted concentrations. GI50 values were obtained from two separate experiments.

To investigate the safety and druglikeness of these new HH inhibitors, we tested their inhibition on hERG channel and predicted their physico-chemical properties. As shown in Table 7, most compounds except 30 and 31 show negligible inhibition on the hERG with IC50 values greater than 20 μmol/L, indicating their low liability of cardiac toxicity. Compounds 30 and 31 have IC50 values of 6.4 and 11.7 μmol/L, indicating their slightly safety concerns on hERG, likely due to their more basic nature since both compounds contain a carboxamide substituent on the pyrrolopyridone component. Meanwhile, most compounds have similar predicated physico-chemical properties, including acceptable tPSA (<140 Å), modest aqueous solubility (logS), and low Caco-2 permeability (Table 7).

Pharmacokinetic parameters of potent compounds

To select optimal compound for further in vitro and in vivo antitumor study, the 7 highly potent HH inhibitors with IC50 values less than 1.0 nmol/L were selected for further evaluation for their pharmacokinetic (PK) parameters in Sprague–Dawley rats dosed intravenously (i.v., 1 mg/kg) and orally (p.o., 3 mg/kg). As shown in Table 8, these compounds generally have short half-life (t1/2) and modest plasma exposure (AUC). Compounds 24, 25 and 35 have relatively acceptable oral bioavailability (20%–30%) and modest clearance. Unfortunately, the most potent HH inhibitors 43 (0.062 nmol/L) and 30 (0.05 nmol/L) display poor PK parameters, especially with scarce oral bioavailability. Therefore, after consideration of the HH inhibitory potency and PK profile, we chose compounds 25 and 35 for further investigation.
Table 8

PK parameters of potent compounds after p.o. (3 mg/kg) and i.v. (1 mg/kg) administration.

Compd.Routet1/2(h)Tmax(h)Cmax(ng/mL)AUClast(h·ng/mL)AUCINF_obs(h·ng/mL)CL_obs(mL/min/kg)MRTINF_obs(h)Vss_obs(mL/kg)F
24p.o.0.950.504744734731.2630.9
i.v.0.5650950933.20.761513
25p.o.0.960.255034064060.9522.0
i.v.0.5461461427.20.771258
30p.o.1.350.2515.29.919.910.891.2
i.v.0.6827827861.50.7942891
31p.o.1.550.3315.518.322.12.382.9
i.v.0.6921321777.60.573512
34p.o.2.061.0029.197.31073.4912.7
i.v.1.0925526176.41.668781
35p.o.1.000.51101901901.5421.0
i.v.1.0630030059.21.204211
43p.o.2.020.423468619243.184.71
i.v.1.05609961352.781.46239
2p.o.3.342.75401203020525.4557.0
i.v.2.741186119517.63.112772

aValues are the average of three runs. Vehicle: p.o., DMSO/0.5% HPMC (5/95, v/v); i.v., EtOH/PEG300/NaCl (10/40/50, v/v/v). CL, clearance; Vss, volume of distribution; t1/2, half-life; AUC, area under the plasma concentration time curve; F, oral bioavailability. –Not applicable.

PK parameters of potent compounds after p.o. (3 mg/kg) and i.v. (1 mg/kg) administration. aValues are the average of three runs. Vehicle: p.o., DMSO/0.5% HPMC (5/95, v/v); i.v., EtOH/PEG300/NaCl (10/40/50, v/v/v). CL, clearance; Vss, volume of distribution; t1/2, half-life; AUC, area under the plasma concentration time curve; F, oral bioavailability. –Not applicable.

Antiproliferative effects of compounds 25 and 35 against various cancer cells

Since the HH signaling pathway is a major regulator governing cell proliferation and differentiation, the selected HH inhibitors 25 and 35 were tested for their antiproliferative effects against a wide spectrum of cancer cell lines, including nine solid and six hematologic tumor cell lines, together with the CHL normal cell line. As shown in Table 9, no or low inhibition on the normal CHL cell lines was observed, especially for compound 25 with an GI50 value greater than 10 μmol/L, indicating the high selectivity against cancer cells. Notably, compared to the solid tumor cells, both compounds show markedly higher potency in the hematologic tumor cell lines, especially for SKM-1, DOHH-2, MOLM-14, MV4-11, MOLM-13 and TMD8 cell lines with GI50 value less than 10 nmol/L. Antiproliferative inhibition of various cancer cells. Antiproliferative effects were initially tested in 10 μmol/L, and then diluted in eight diluted concentrations. GI50 values were obtained from two separate experiments.

Inhibition of HH signaling pathway at the level of GLI

To investigate the exact mechanism of compounds 25 and 35 within the HH signaling pathway, we investigated the suppressing effects of both compounds in the HH pathway by overexpressing GLI1 or GLI2. As shown in Fig. 3, compound 25 significantly inhibited GLI-luciferase activity provoked by forced expression of GLI1 or GLI2, with IC50 values of 0.97 and 2.2 nmol/L, respectively (Fig. 3A and B). Meanwhile, compound 25 dramatically reduced the abundance of endogenous GLI2 proteins in light2 cells at concentration of both 10 and 100 nmol/L, and the effect on GLI2 proteins of compound 25 (100 nmol/L) is superior to JQ1 (1 μmol/L), whereas the SMO inhibitor GDC-0449 (vismodegib) has no effect on the expression of endogenous GLI2 (Fig. 3C). In addition, compound 25 failed to inhibit the expression of exogenous GLI2-MYC (Fig. 3D).
Figure 3

Compounds 25 and 35 inhibit HH activity at the level of GLI1/2. (A), (B) and (E), (F): GLI-luciferase activity in light2 cells treated with 25 or 35 after transfection with GLI1 or GLI2 construct; (C) and (G): Western blots evaluating endogenous GLI2 expression in light2 cells after treatment with 25, 35, JQ1 (1 μmol/L), and GDC-0449 (100 nmol/L); (D) and (H): GLI2-MYC protein expressions in light2 cells transfected with GLI2-MYC; (I): Western blots analysis of exogenous GLI1-FLAG.

Compounds 25 and 35 inhibit HH activity at the level of GLI1/2. (A), (B) and (E), (F): GLI-luciferase activity in light2 cells treated with 25 or 35 after transfection with GLI1 or GLI2 construct; (C) and (G): Western blots evaluating endogenous GLI2 expression in light2 cells after treatment with 25, 35, JQ1 (1 μmol/L), and GDC-0449 (100 nmol/L); (D) and (H): GLI2-MYC protein expressions in light2 cells transfected with GLI2-MYC; (I): Western blots analysis of exogenous GLI1-FLAG. Similarly, compound 35 also inhibited GLI-luciferase activity provoked by forced expression of GLI1 or GLI2, with IC50 values of 0.30 and 0.43 nmol/L, respectively (Fig. 3E and F). Meanwhile, compound 35 down-regulated endogenous GLI2 proteins in light2 cells (Fig. 3G), whereas no effect was observed on the expression of exogenous GLI1-FLAG and GLI2-MYC (Fig. 3H and I). Given the fact that significant suppression on the HH activation provoked by forced expression of GLI1 or GLI2 was observed by compounds 25 and 35, but not by SMO inhibitor, as well as the robust inhibition of 25 and 35 on the activity of BRD proteins, we conclude that compounds 25 and 35 epigenetically inhibit HH signaling pathway at the level of GLI, not at the upstream targets of the HH signaling pathway, e.g., SMO.

Compounds 25 and 35 can circumvent drug resistance caused by SMO mutations

Having characterized that compounds 25 and 35 function as HH inhibitors by acting at the level of GLI, we continued to evaluate its in vitro ability of combating the resistance to SMO inhibitors caused by SMO mutation. As expected, both compounds markedly blocked GLI-luciferase activity in light2 cells induced by ectopic expression of wild type SMO (SMOWT, Fig. 4A and E), as well as SMOD473H (Fig. 4B and F) and SMOW535L (Fig. 4C and G), two most predominant mutants of SMO responsible for the resistance to SMO inhibitors, with respective IC50 values of 1.16, 1.44 and 1.13 nmol/L for 25, and 0.39, 0.53 and 0.63 nmol/L for 35 (Fig. 4D and H). These findings indicate that compounds 25 and 35 are promising HH inhibitors with capacity of overcoming drug resistance caused by SMO mutations.
Figure 4

Compounds 25 and 35 possess in vitro ability of circumventing the resistance to SMO inhibitors caused by SMO mutations. GLI-luciferase analysis of the inhibitory effects of compounds in light2 cells followed by transfection with SMOWT (A) and (E), SMOD473H (B) and (F), or SMOW535L (C) and (G), respectively. (D) and (H) IC50 values of compounds in light2 cells with forced expression of SMOWT, SMOD473H and SMOW535L.

Compounds 25 and 35 possess in vitro ability of circumventing the resistance to SMO inhibitors caused by SMO mutations. GLI-luciferase analysis of the inhibitory effects of compounds in light2 cells followed by transfection with SMOWT (A) and (E), SMOD473H (B) and (F), or SMOW535L (C) and (G), respectively. (D) and (H) IC50 values of compounds in light2 cells with forced expression of SMOWT, SMOD473H and SMOW535L.

In vivo antitumor effects in SMOA1 MB allograft model

To translate the in vitro ability of the GLI inhibitors 25 and 35 in combating the resistance caused by SMO inhibition into in vivo antitumor efficacy, we established spontaneous SMOA1 MB model, of which tumors contain mutation at the amino acid residue 539 of SMO (SMOW539L), and are resistant to current SMO inhibitor vismodegib (GDC-0449) in ND2:SMOA1 mice. As shown in Fig. 5A, treatment with compound 25 at 10 mg/kg by oral administration (p.o., bid) realized 88.2% tumor growth inhibition (TGI), whereas the SMO inhibitor GDC-0449 failed to affect the tumor growth at 20 mg/kg (p.o., bid), a dosage resulting in complete tumor growth inhibition in SMO wild MB mice. However, despite its higher potency of 35 in vitro than compound 25, treatment of this compound at the same dosage was not tolerated leading to significant casualty (data are not shown). However, a lower dose of 35 at 1 mg/kg was well tolerant with TGI of 68.7%. No obvious body weight loss was observed in both the 10 mg/kg group of 25 and the 1 mg/kg group of 35 during treatment (Fig. 5B). These results suggest that both compounds 25 and 35 are capable of combating SMO inhibitor-resistant MB tumors.
Figure 5

Antitumor activity of compounds 25 and 35 in SMOA1 MB allograft model. Nude mice with subcutaneously transplanted SMOA1 MB allografts were treated with vehicle control (0.5% CMC-Na, p.o., bid), compounds 25 and 35 (p.o., bid) or GDC-0449 (p.o., bid). (A) Tumor volume and (B) body-weight change over time in treated SMOA1 MB allograft. Data represent group means ± SEM (control, n = 4; GDC-0449, n = 4; 35, n = 4; 25, n = 3).

Antitumor activity of compounds 25 and 35 in SMOA1 MB allograft model. Nude mice with subcutaneously transplanted SMOA1 MB allografts were treated with vehicle control (0.5% CMC-Na, p.o., bid), compounds 25 and 35 (p.o., bid) or GDC-0449 (p.o., bid). (A) Tumor volume and (B) body-weight change over time in treated SMOA1 MB allograft. Data represent group means ± SEM (control, n = 4; GDC-0449, n = 4; 35, n = 4; 25, n = 3). Since compound 25 functioned as an efficacious and more tolerant GLI inhibitor than compound 35 to overcome the resistant MB tumors caused by SMOW539L mutation, we continued to evaluate this compound at lower dosage. As shown in Fig. 6A, compound 25, even at a lower dosage of 5 mg/kg (p.o., bid) still shows 83.3% tumor growth inhibition, similar to the effect of 10 mg/kg dose (TGI = 87.6%). Meanwhile, the new HH inhibitor 25 was well tolerant in both doses without significant overall toxicity (Fig. 6B).
Figure 6

Antitumor activity of compound 25 in SMOA1 MB allograft model. Nude mice with subcutaneously transplanted SMOA1 MB allografts were treated with vehicle control (0.5% CMC-Na, p.o., bid), or compound 25 (p.o., bid). (A) Tumor volume and (B) body-weight change over time in treated SMOA1 MB allograft. Data represent group means ± SEM (control, n = 4; 25, 5 mg/kg, n = 4; 25, 10 mg/kg, n = 5).

Antitumor activity of compound 25 in SMOA1 MB allograft model. Nude mice with subcutaneously transplanted SMOA1 MB allografts were treated with vehicle control (0.5% CMC-Na, p.o., bid), or compound 25 (p.o., bid). (A) Tumor volume and (B) body-weight change over time in treated SMOA1 MB allograft. Data represent group means ± SEM (control, n = 4; 25, 5 mg/kg, n = 4; 25, 10 mg/kg, n = 5).

Conclusions

In summary, by repurposing a series of the clinically investigational BRD4 inhibitors as the HH signaling pathway inhibitors, the AbbVie's phase I clinical pan-BET inhibitor 2 (ABBV-075) was selected for a systemic medicinal chemistry optimization. The established SAR indicates that the inhibitory activity against BRD4 does not parallel with the potency for the HH pathway. Notably, BRD4 is well tolerant to various substitution patterns, whereas the activity for the HH pathway is much different and highly dependent on the steric and electrostatic natures of the substituents. A number of compounds were identified showing IC50 values less than 1.0 nmol/L against the HH signaling pathway, and further investigation of their inhibition on hERG and PK parameters elected compounds 25 and 35 as the high potency HH inhibitors. Similar to compound 2, these new compounds are non-selective BET inhibitors (Supporting Information Fig. S2 and Table S2) that might be more effective to modulate GLI activity. Further mechanism profiling showed that the two new compounds suppressed HH signaling by interaction with the transcriptional factor GLI, and are equally potent against both the two clinically identified resistant SMO mutants (SMOD473H and SMOW535L) and SMO wild type, with IC50 values around 1 nmol/L or less. Although the PK parameters are not optimal, both compounds showed significant antitumor efficacy in vivo in the resistant MB allograft mice. Compared to the dose-limited toxicity of compound 35, compound 25 was more tolerant and significantly suppressed tumor growth at both 5 mg/kg (TGI = 83.3%) and 10 mg/kg (TGI = 87.6%) doses. Therefore, the new GLI inhibitor 25 represents not only a high potent HH inhibitor both in vitro and in vivo, but also a relatively safe and tolerant lead compound capable of suppressing HH-driven MB, that is resistant to current SMO antagonists.

Experimental

All reactions were performed in glassware containing a teflon-coated stir bar. All commercially purchased reagents and solvents were chemical pure and used without futher purification. 1H NMR and 13C NMR spectra were recorded on a Varian Mercury 300, 400 or 500 NMR spectrometer (Bruker Biospin AG, Romanshorn, Switzerland) and referenced to deuterium dimethyl sulfoxide (DMSO-d6), deuterium chloroform (CDCl3), deuterium methanol (CD3OD) or deuterium dichloromethane (CD2Cl2). Chemical shifts (δ) were reported in ppm downfield from an internal TMS standard. High-resolution mass spectrometry (HRMS) analysis was recorded at anionizing voltage of 70 eV on a Finnigan/MAT95 spectrometer (Agilent, Palo Alto, CA, USA). Flash column chromatography on silica gel (200–300 mesh, Yucheng chemical, Shanghai, China) was used for the routine purification of reaction products. All reactions were monitored by TLC on silica gel plates (50 mm × 15 mm) and spots were visualized under UV light (Changbo, Shanghai, China). Melting points were determined using a SGW X-4 hot stage microscope (INESA, Shanghai, China) and are uncorrected. HPLC analysis was conducted for all bioassayed compounds on an Agilent Technologies 1260 series LC system (ZORBAX-C18 (150 mm × 4.6 mm, Daicel, Shanghai, China), 5 μmol/L, MeOH/H2O or MeCN/H2O, r.t.) with two ultraviolet wavelengths (UV 254 and 214 nm). The purities of these compounds were above 95%.

Synthetic procedures

Synthesis of target compounds and key intermediates was shown in Supporting Information. Compounds 46, 51, 60a, 66 and 75 were commericially purchased, and compounds 47, 53, 55, 60 and 69 were prepared according to corresponding literature procedures.

Biological studies

Homogeneous time-resolved fluorescence

The binding of compounds to BRD4 was assessed using an homogeneous time-resolved fluorescence assay. Recombinant BRD4 proteins (Active Motif, Cat#31380 and 31446, Carlsbad, CA, USA), compounds and peptide were diluted to desired working concentrations with EPIgenerous binding buffer (Cisbio, Cat#62DLBDDF, Codolet, France). 4 μL of 5 × BRD4 protein at 165 ng/mL, 4 μL of 5 × biotin labeled acetylated H4K5/K8/K12/K16 at 200 ng/mL and 2 μL of serial diluted test compounds were added to the wells of a 384 ProxiPlate (PerkinElmer, Cat#6008289, Waltham, MA, USA). After 30 min incubation at 37 °C, a 10 μL mixture of IgG antibody labeled with Eu and streptavidin labeled with XL665 were prepared and added to each well for 3 h to reach equilibrium at room temperature. The signal of fluorescence was detected on Envison (PerkinElmer) and analyzed in GraphPad Prism (San Diego, CA, USA).

Dual luciferase reporter assay

Light2 cells (ATCC, Manassas, VA, USA) were treated with compounds as indicated for 36 h. Relative-luciferase assays were performed using dual-luciferase reporter assay kit (Promega, Madison, WI, USA) according to manufacturer's instructions. Relative luciferase values were normalized to TK-renilla value.

Cell growth inhibition assays

SKM-1 cells (JCRB, Osaka, Japan) were seeded into a 96-well plate at a suitable density in a volume of 100 μL medium. After incubation overnight, compounds dissolved in DMSO stock solutions were thawed at room temperature and diluted to the desired concentrations with saline. The compounds were added and cultured for 72 h, the IC50 (or GI50) value was measured with the Cell Counting Kit-8 (CCK-8) assay (Dojindo, Kumamoto, Japan).

GLI protein expression assays

Light2 cells were transfected with MYC-GLI2-FL or FLAG-GLI1 (Addgene, Cambridge, MA, USA) using Lipo2000 (Invitrogen, Carlsbad, CA, USA) as manufacturer's protocol. Cells were treated with indicated compounds 24 h after transfection.

Western blot analysis

Cells were lysed with LDS Sample Buffer (NP0008, Life, Waltham, MA, USA). The protein samples were analyzed with SDS-PAGE gel (Cell Signaling Technology, Boston, MA, USA) and immunoblotted with the following primary antibodies: anti-HA, anti-FLAG, anti-GLI2 (1:1000, Cell Signaling Technology) or GAPDH (1:5000, Santa Cruz Biotechnology, Santa Cruz, CA, USA).

Medulloblastoma allograft model

SMOA1 medulloblastoma model was derived from ND2:SMOA1 transgenic mice (Shanghai SLAC Laboratory Animal Co., Shanghai, China) as previously described. Briefly, primary SMOA1 allografts were subcutaneously grafted into nude mice. Mice allografted with SMOA1 MB were administered compounds GDC-0449, 25 and 35 or vehicle twice a day by oral administration (p.o., n ≥ 3) once the tumor size reached 100–150 mm3. Tumor volumes were measured with calipers every 3 days and calculated as 0.5 × length × width. All animal procedures were approved by Animal Care and Use Committee of Fudan University, Shanghai, China.

Acknowledgments

This work was supported by grants from (Grants Nos. 81773565, 81703327, 81430080, 81573452, and 81773767). Supporting grants from the of the Frontier Science (Grant No. 160621, China) of the , the Strategic Priority Research Program of the (Grant No. XDA12020374, China) are also appreciated. We also thank the support from the State Key Laboratory of Esophageal Cancer Prevention and Treatment, , (Grant No. K2020-0012, China). In addition, we thank Prof. Zhaobing Gao and Dr. Xueqin Chen of Shanghai Institute of Materia Medica (SIMM, Shanghai, China), and Prof. Qingsong Liu and Dr. Aoli Wang of High Magnetic Field Laboratory, Chinese Academy of Sciences (Hefei, China) for the test of hERG inhibition and cellular antiproliferative activity, respectively.

Author contributions

Ao Zhang, Wenfu Tan and Hongchun Liu conceived and designed the study. Xiaohua Liu performed the synthesis of the target compounds and analyzed the data. Yu Zhang, Yalei Li and Juan Wang performed the biological experiments. Huaqian Ding and Wenjing Huang assisted with carrying out experiments. Ao Zhang and Xiaohua Liu wrote the manuscript. Chunyong Ding, Wenfu Tan, Hongchun Liu and Juan Wang revised the manuscript. All the authors approved the final version of the manuscript.

Conflicts of interest

The authors have no conflicts of interest to declare.
  45 in total

Review 1.  Discovery and preclinical development of vismodegib.

Authors:  Stephen E Gould; Jennifer A Low; James C Marsters; Kirk Robarge; Lee L Rubin; Frederic J de Sauvage; Daniel P Sutherlin; Harvey Wong; Robert L Yauch
Journal:  Expert Opin Drug Discov       Date:  2014-05-23       Impact factor: 6.098

Review 2.  Medulloblastoma.

Authors:  Paul A Northcott; Giles W Robinson; Christian P Kratz; Donald J Mabbott; Scott L Pomeroy; Steven C Clifford; Stefan Rutkowski; David W Ellison; David Malkin; Michael D Taylor; Amar Gajjar; Stefan M Pfister
Journal:  Nat Rev Dis Primers       Date:  2019-02-14       Impact factor: 52.329

3.  Optimization of a Series of Bivalent Triazolopyridazine Based Bromodomain and Extraterminal Inhibitors: The Discovery of (3R)-4-[2-[4-[1-(3-Methoxy-[1,2,4]triazolo[4,3-b]pyridazin-6-yl)-4-piperidyl]phenoxy]ethyl]-1,3-dimethyl-piperazin-2-one (AZD5153).

Authors:  Robert H Bradbury; Rowena Callis; Gregory R Carr; Huawei Chen; Edwin Clark; Lyman Feron; Steve Glossop; Mark A Graham; Maureen Hattersley; Chris Jones; Scott G Lamont; Gilles Ouvry; Anil Patel; Joe Patel; Alfred A Rabow; Craig A Roberts; Stephen Stokes; Natalie Stratton; Graeme E Walker; Lara Ward; David Whalley; David Whittaker; Gail Wrigley; Michael J Waring
Journal:  J Med Chem       Date:  2016-08-24       Impact factor: 7.446

4.  Smoothened variants explain the majority of drug resistance in basal cell carcinoma.

Authors:  Scott X Atwood; Kavita Y Sarin; Ramon J Whitson; Jiang R Li; Geurim Kim; Melika Rezaee; Mina S Ally; Jinah Kim; Catherine Yao; Anne Lynn S Chang; Anthony E Oro; Jean Y Tang
Journal:  Cancer Cell       Date:  2015-03-09       Impact factor: 31.743

5.  Discovery of N-(4-(2,4-Difluorophenoxy)-3-(6-methyl-7-oxo-6,7-dihydro-1H-pyrrolo[2,3-c]pyridin-4-yl)phenyl)ethanesulfonamide (ABBV-075/Mivebresib), a Potent and Orally Available Bromodomain and Extraterminal Domain (BET) Family Bromodomain Inhibitor.

Authors:  Keith F McDaniel; Le Wang; Todd Soltwedel; Steven D Fidanze; Lisa A Hasvold; Dachun Liu; Robert A Mantei; John K Pratt; George S Sheppard; Mai H Bui; Emily J Faivre; Xiaoli Huang; Leiming Li; Xiaoyu Lin; Rongqi Wang; Scott E Warder; Denise Wilcox; Daniel H Albert; Terrance J Magoc; Ganesh Rajaraman; Chang H Park; Charles W Hutchins; Jianwei J Shen; Rohinton P Edalji; Chaohong C Sun; Ruth Martin; Wenqing Gao; Shekman Wong; Guowei Fang; Steven W Elmore; Yu Shen; Warren M Kati
Journal:  J Med Chem       Date:  2017-10-12       Impact factor: 7.446

6.  Exploitation of Castration-Resistant Prostate Cancer Transcription Factor Dependencies by the Novel BET Inhibitor ABBV-075.

Authors:  Emily J Faivre; Denise Wilcox; Xiaoyu Lin; Paul Hessler; Maricel Torrent; Wei He; Tamar Uziel; Daniel H Albert; Keith McDaniel; Warren Kati; Yu Shen
Journal:  Mol Cancer Res       Date:  2016-10-05       Impact factor: 5.852

7.  Discovery of Novel Macrocyclic Hedgehog Pathway Inhibitors Acting by Suppressing the Gli-Mediated Transcription.

Authors:  Gang Liu; Wenjing Huang; Juan Wang; Xiaohua Liu; Jun Yang; Yu Zhang; Yong Geng; Wenfu Tan; Ao Zhang
Journal:  J Med Chem       Date:  2017-09-20       Impact factor: 7.446

Review 8.  Targeting GLI factors to inhibit the Hedgehog pathway.

Authors:  Paola Infante; Romina Alfonsi; Bruno Botta; Mattia Mori; Lucia Di Marcotullio
Journal:  Trends Pharmacol Sci       Date:  2015-06-10       Impact factor: 14.819

9.  Intertumoral Heterogeneity within Medulloblastoma Subgroups.

Authors:  Florence M G Cavalli; Marc Remke; Ladislav Rampasek; John Peacock; David J H Shih; Betty Luu; Livia Garzia; Jonathon Torchia; Carolina Nor; A Sorana Morrissy; Sameer Agnihotri; Yuan Yao Thompson; Claudia M Kuzan-Fischer; Hamza Farooq; Keren Isaev; Craig Daniels; Byung-Kyu Cho; Seung-Ki Kim; Kyu-Chang Wang; Ji Yeoun Lee; Wieslawa A Grajkowska; Marta Perek-Polnik; Alexandre Vasiljevic; Cecile Faure-Conter; Anne Jouvet; Caterina Giannini; Amulya A Nageswara Rao; Kay Ka Wai Li; Ho-Keung Ng; Charles G Eberhart; Ian F Pollack; Ronald L Hamilton; G Yancey Gillespie; James M Olson; Sarah Leary; William A Weiss; Boleslaw Lach; Lola B Chambless; Reid C Thompson; Michael K Cooper; Rajeev Vibhakar; Peter Hauser; Marie-Lise C van Veelen; Johan M Kros; Pim J French; Young Shin Ra; Toshihiro Kumabe; Enrique López-Aguilar; Karel Zitterbart; Jaroslav Sterba; Gaetano Finocchiaro; Maura Massimino; Erwin G Van Meir; Satoru Osuka; Tomoko Shofuda; Almos Klekner; Massimo Zollo; Jeffrey R Leonard; Joshua B Rubin; Nada Jabado; Steffen Albrecht; Jaume Mora; Timothy E Van Meter; Shin Jung; Andrew S Moore; Andrew R Hallahan; Jennifer A Chan; Daniela P C Tirapelli; Carlos G Carlotti; Maryam Fouladi; José Pimentel; Claudia C Faria; Ali G Saad; Luca Massimi; Linda M Liau; Helen Wheeler; Hideo Nakamura; Samer K Elbabaa; Mario Perezpeña-Diazconti; Fernando Chico Ponce de León; Shenandoah Robinson; Michal Zapotocky; Alvaro Lassaletta; Annie Huang; Cynthia E Hawkins; Uri Tabori; Eric Bouffet; Ute Bartels; Peter B Dirks; James T Rutka; Gary D Bader; Jüri Reimand; Anna Goldenberg; Vijay Ramaswamy; Michael D Taylor
Journal:  Cancer Cell       Date:  2017-06-12       Impact factor: 31.743

10.  An Investigator-Initiated Open-Label Trial of Sonidegib in Advanced Basal Cell Carcinoma Patients Resistant to Vismodegib.

Authors:  Christina Danial; Kavita Y Sarin; Anthony E Oro; Anne Lynn S Chang
Journal:  Clin Cancer Res       Date:  2015-11-06       Impact factor: 12.531

View more
  6 in total

1.  Challenges and Recent Advances of Novel Chemical Inhibitors in Medulloblastoma Therapy.

Authors:  Anand Maurya; Upendra Kumar Patel; Jitendra Kumar Yadav; Virender Pratap Singh; Alka Agarwal
Journal:  Methods Mol Biol       Date:  2022

Review 2.  Small molecules in targeted cancer therapy: advances, challenges, and future perspectives.

Authors:  Lei Zhong; Yueshan Li; Liang Xiong; Wenjing Wang; Ming Wu; Ting Yuan; Wei Yang; Chenyu Tian; Zhuang Miao; Tianqi Wang; Shengyong Yang
Journal:  Signal Transduct Target Ther       Date:  2021-05-31

Review 3.  Hedgehog Pathway Inhibitors against Tumor Microenvironment.

Authors:  Silpa Gampala; Jer-Yen Yang
Journal:  Cells       Date:  2021-11-12       Impact factor: 6.600

Review 4.  Transcription Factors with Targeting Potential in Gliomas.

Authors:  Angeliki-Ioanna Giannopoulou; Dimitrios S Kanakoglou; Christina Piperi
Journal:  Int J Mol Sci       Date:  2022-03-28       Impact factor: 5.923

Review 5.  Molecular Mechanisms Involving the Sonic Hedgehog Pathway in Lung Cancer Therapy: Recent Advances.

Authors:  Chao Ma; Kang Hu; Irfan Ullah; Qing-Kang Zheng; Nan Zhang; Zhi-Gang Sun
Journal:  Front Oncol       Date:  2022-04-01       Impact factor: 5.738

Review 6.  Trafficking and retention of protein antigens across systems and immune cell types.

Authors:  Thu A Doan; Tadg Forward; Beth A Jirón Tamburini
Journal:  Cell Mol Life Sci       Date:  2022-05-03       Impact factor: 9.207

  6 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.