Skip to main content

Flubendazole inhibits PD-1 and suppresses melanoma growth in immunocompetent mice

Abstract

Background

Immune checkpoint inhibitor therapy has revolutionized the clinical management of a diverse range of cancer types, including advanced cutaneous melanoma. While immunotherapy targeting the PD-1/PD-L1 system has become standard of care, overall response rates remain unsatisfactory for most patients and there are no approved small molecule inhibitors of the PD-1/PD-L1 system. Flubendazole (FLU) is an anthelmintic that has been used to treat worm infections in humans and animals for decades.

Methods

Here we tested the anti-cancer activity of systemically delivered FLU with suppression of PD-1 in immunocompetent mice.

Results

In C57BL/6J mice bearing subcutaneous B16F10 melanoma, FLU reduced both tumor growth and PD-1 protein levels without affecting levels of PD-L1. FLU’s suppression of PD-1 was accompanied by increased CD3+ T cell infiltration. Western blotting with extracts from human Jurkat T cells showed that FLU inhibited PD-1 protein expression, findings confirmed by flow cytometry. To gain mechanistic insights on FLU’s ability to suppress PD-1 protein levels, we performed bulk RNA sequencing on extracts of Jurkat T cells exposed to the benzimidazole for 4 h. From a pool of 14,475 genes there were 1218 differentially-expressed genes; 687 with increased expression and 531 with decreased expression. Among the genes induced by FLU was the AP-1 family member, JUN and surprisingly, pdcd1. KEGG pathway analysis showed FLU up-regulated genes over-represented in multiple pathways (p < 0.01), the top hit being amoebiasis. FLU also affected the expression of genes in cancer-associated pathways, both through down-regulation and up-regulation. Gene set enrichment analysis revealed a large number of immunological signature gene sets correlated with FLU treatment, including gene sets associated with T cell differentiation, proliferation and function. The AP-1 inhibitor T5224 rescued PD-1 protein expression from inhibition by FLU.

Conclusion

This study is the first to show that FLU can inhibit melanoma growth with PD-1 suppression in immunocompetent mice.

Introduction

Immune checkpoint therapy has radically changed the management of a diverse range of cancer types, including advanced cutaneous melanoma [1]. For example, antibodies to programmed cell death protein 1 (PD-1) (such as nivolumab/Opdivo, pembrolizumab/Keytruda and cemiplimab/Libtayo), PD-L1 (atezolizumab/Tecentriq, avelumab/Bavencio and durvalumab/Imfinzi) and cytotoxic T-lymphocyte–associated antigen 4 (CTLA-4) (ipilimumab/Yervoy) have become standard of care in tumor immunotherapy. However, around 1 in 2 patients do not respond to PD-1 antibody-based therapies or the durability of this response is not sustained and immune-related adverse events are common [2,3,4,5]. There are no approved small molecule inhibitors of PD-1, PD-L1 or CTLA-4 that could be used as alternatives to antibodies. Small molecules offer potential advantages over antibodies such as greater stability, medicinal chemistry, low cost, penetrant effects across whole cell populations, intracellular signalling pathway targeting and potential oral administration meaning avoidance of in-clinic administration and reduced health costs [6].

Drug repurposing, the process by which older medicines can be redeveloped for new therapeutic indications offers a range of advantages over conventional drug discovery programs [7]. Repurposing can develop de-risked medicines with known pharmacology, tolerance and toxicity at potential lower cost. A recent study found that the median capitalized R&D cost of bringing a new drug to market was approximately $985 million and > 8 years per drug [8]. A well-known example of a repurposed drug is sildenafil, which was originally developed in the context of angina and is widely used for erectile dysfunction [9]. A more recent example is remdesivir which was developed for Ebola virus and is used for the management of SARS-CoV-2 infection [10].

PD-1 (CD279), encoded by the pdcd1 gene and first discovered in 1992 [11], is an inhibitory receptor expressed by T and B cells, natural killer cells and certain myeloid cell populations [12]. We and others have shown that PD-1 is also expressed by melanoma cells [13, 14]. PD-1, a type I transmembrane glycoprotein bound by glycoproteins PD-L1 (B7-H1, CD274) and PD-L2 (B7-DC, CD273) [15, 16] can promote tumor growth by preventing activation of cytotoxic T-cells [13, 14]. Given current limitations of antibody-based therapies, small molecule alternatives may provide an alternative means to disrupt the PD-1/PD-L1 axis. For example, Wang et al. recently reported their discovery of compound 17 which inhibits the PD-1/PD-L1 interaction and results in degradation of PD-L1 [17]. Studies in our laboratory show that flubendazole (FLU), a benzimidazole used as an anthelmintic for over 40 years [18], can inhibit human melanoma growth and metastasis in immunodeficient mice [14]. Others have also shown that FLU can cause mitotic catastrophe [19], disrupt cell cycle progression [20] and promote ferroptosis [21]. However, it is unclear whether FLU can affect PD-1 expression in T cells or influence the growth of melanoma in immunocompetent mice. In this paper we show that FLU can inhibit melanoma growth and PD-1 expression in immunocompetent mice and that FLU’s inhibition of PD-1 involves the transcription factor AP-1.

Materials and methods

Mouse tumor experiments

Six to eight week old female C57BL/6 J mice were sourced from Animal Resources Centre (Perth). Animals were inoculated subcutaneously with B16F10 cells (1 × 105 cells/animal in 100 µl DMEM containing 10% FBS and 50% Matrigel) into the shaven right flank on Day 0. These cells were obtained from ATCC. Mice were randomized on Day 5 into groups with mean tumor volume 40–45 mm3 per group. Treatments were started on Day 7, vehicle or FLU (Advanced Molecular Technologies, Scoresby) was administered intraperitoneally (20 ml/kg or 200 mg/kg, i.p.) once daily on a 5 days-on/2 days-off schedule. FLU was suspended at 10 mg/ml in vehicle (saline (0.9% NaCl) with 0.5% Tween 80 and 0.01% DMSO) and sonicated. Animals were weighed daily after Day 5. The general health condition and attitude of each animal was monitored daily throughout the study. Tumors were measured by length, height and width in millimeters daily. Tumor volumes were calculated using formula V = L × H × W × π/6. If a second tumor occurred in an animal, both tumor volumes were measured and their volumes were combined. After euthanasia by isoflurane overdose and cervical dislocation, tumors of each animal were fixed in 10% neutral buffered formalin. The protocol was approved by the UNSW Animal Care and Ethics Committee.

Cell culture, Western blotting and flow cytometry

Human Jurkat T cells grown in RPMI 1640 with, pH 7.4, containing 10% FBS and penicillin/streptomycin were seeded into 6 well plates and treated with FLU or vehicle for the indicated times. Where indicated, T5224 (cat. S8966, Selleck Chem) or vehicle was added 2 h prior to the addition of FLU. Total cell extracts were prepared in RIPA buffer and lysates were resolved by SDS-PAGE and Western blot analysis with primary antibodies target PD-1 (1:1000, cat. ab214421, Abcam; 1:500, cat. 86163, CST, which was used for the T5224 studies), JUN (1:500, cat. ab32137; Abcam), FOS (1:500, cat. no. 2250; CST); or β-actin (1:3000, cat. A5316, Sigma) and followed by chemiluminescence detection using the Western Lightning Chemiluminescence system (ThermoFisher, cat. A38556) and visualized with an ImageQuant LAS 4000 Biomolecular Imager.

Alternatively, cells were washed and resuspended (5 × 106 cells/ml) in Stain Buffer. One hundred µl was added to 12 × 75 mm tubes and 5 µl of BV421 conjugated mouse anti-human CD279 (PD-1) (1:20, BD, cat. 562,516) or BV421 conjugated mouse IgG1 (1:20, BD, cat. 562,438) was added and then incubated for 45 min at 22 °C, protected from light. Cells were washed twice with 2 ml of Stain Buffer, centrifuged and pellets were resuspended in 0.5 ml Stain Buffer. Stained cell suspensions were analysed by flow cytometry using a BD LSR Fortessa X20.

Immunohistochemical processing and staining

Formalin-fixed, paraffin embedded sections were prepared from tumors. Heat-induced epitope retrieval was applied to all deparaffinized Sects. (4 μm Superfrost slides) with citrate buffer, pH 6.0 for 5 min at 110 °C. Sections were blocked with endogenous AP (levamisole) blocking agent (DAKO, S2003) for 10 min and then with 2% skim milk for 20 min. Rabbit monoclonal anti-PD-1 (cat. ab214421) and rabbit monoclonal anti-CD3 antibody (SP7) (cat. ab16669) were obtained from Abcam. Rabbit polyclonal anti-PD-L1 (cat. PA5-20,343) antibodies were obtained from ThermoFisher. Slides were incubated with primary antibody for 1 h at room temperature and then for 10 min with MACH3 Rabbit AP-Polymer Detection solution (probe incubation) (Biocare Medical, M3R533 G, H, L). After rinsing with buffer, the slides were incubated with MACH3 Rabbit AP-Polymer Detection solution (polymer incubation) (Biocare Medical, M3R533 G, H, L) for a further 10 min. Slides were incubated with red chromogen (Warp RedTM Chromogen Kit) for 5 min and counterstained in haematoxylin and Scott blue. Slides were dried with filter paper and dehydrated in xylene then coverslipped.

Immunostained slides were scanned using an Aperio ScanScope XT slide scanner (Leica Biosystems, Mt Waverley, Vic, Australia) and images were captured using ImageScope software (Leica Biosystems). Tissue area (µm2) and integrated optical density (IOD) of positive staining (red chromogen) using Image-Pro Plus software (Cybernetics, Bethesda, MD, USA). Positive staining cells numbers were manually counted for CD3 using Image-Pro Plus software (Cybernetics, Bethesda, MD, USA).

RNA-seq and bioinformatics analysis

Jurkat T cells grown to confluence in 100 mm plates in complete medium were incubated with 10 µM FLU or vehicle for 4 h. Total RNA was extracted using the RNeasy Mini Kit (Qiagen, cat. 74,004) with minor modification. Briefly, the cells were washed twice with cold 1 × PBS and TRIzol reagent (ThermoFisher, cat. 15,596,026) was used to lyse the cells. Chloroform was added to the mixture prior to microfuge centrifugation at 13,000 rpm for 15 min at 4 °C. The upper aqueous layer (containing total RNA) was transferred to microtubes, isopropanol was added and loaded into RNeasy columns. Columns were washed with buffers RPE and RW1. Total RNA was eluted using ribonuclease-free water. Six samples (3 biological replicates of each condition) were submitted to the UNSW Ramaciotti Centre for Genomics for TruSeq Stranded mRNA-seq preparation and sequencing by NextSeq 6000 to produce 75 bp single end reads.

RNA-seq reads were assessed for quality using FastQC (v0.11.8) www.bioinformatics.babraham.ac.uk/projects/fastqc/). Subread was used for aligning reads to the Ensembl human genome (GRCh38) and then the featureCounts function of Subread was used to quantify reads (http://subread.sourceforge.net/). Reads assigned to gene features ranged from 46.6 M to 53.8 M reads per sample. The R packages edgeR and limma (voom) [22] were used to identify differentially-expressed genes (DEG) comparing the FLU treated Jurkat T cells with vehicle-treated Jurkat T cells. Lowly expressed genes were filtered out leaving those genes with at least 1 count per million (CPM) in at least 3 samples (n = 14,475) for further analysis. The functions lmFit, eBayes and treat were used to identify genes which differed in the FLU treated cells with a least a log2 (1.2) fold change. Significantly DEG were taken as having an adjusted p < 0.05. R plotting functions including plotMDS (limma) and pheatmap R package were used to prepare plots. The kegga function (limma) was used to perform KEGG pathway analysis.

Gene set enrichment analysis (GSEA)

Gene counts for 27,729 genes with a symbol were normalized using the TMM function of edgeR [23] and used as input into the desktop version of GSEA (v4.3.2, Broad Institute). GSEA then generated a ranked list and analyzed the list to ascertain the degree to which a gene set was overrepresented at the top or bottom of the list, thus generating a normalized enrichment score (NES) for each gene set. Gene set analysis was run using the C7 immunologic signature gene set collection and C2 curated gene sets from MSigDb. Gene set enrichment was ranked by NES and gene sets with a false discovery rate (FDR) < 0.25 were considered significant.

Statistics

Analysis was performed using Graphpad PRISM v9. If distribution was not normal, Mann–Whitney or Kruskal–Wallis test was performed. Normally distributed data was analyzed by t test or one-way ANOVA. Plotted data represent mean ± SEM. Differences were considered significant when p ≤ 0.05. Where indicated, *p ≤ 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001.

Results

FLU suppresses melanoma growth in immunocompetent mice

We recently determined that FLU can inhibit melanoma growth as xenografts in immunocompromised mice [14]. To explore whether FLU can influence melanoma growth in immunocompetent animals, we used a tumor isograft model in which FLU was administered systemically (i.p.) to C57BL/6J mice bearing subcutaneous B16F10 tumors and commenced treatment after the tumors were established. The C57BL/6-B16F10 system is a PD-1 antibody resistant mouse model of melanoma growth [13]. Ten days after melanoma inoculation, tumors grew rapidly in the vehicle group (Fig. 1A, B). On the other hand, FLU delayed melanoma growth by several days and caused near cessation of growth after 15 days, with tumors in the vehicle group measuring ~ 1000 mm3 while those treated with FLU measuring ~ 200 mm3 (Fig. 1B). FLU reduced tumor growth rate between Day 15 to 17 (Fig. 1C). FLU enabled survival in 100% of melanoma bearing mice after 17 days, whereas survival of tumor-bearing mice treated with vehicle after this time was 0% (Fig. 1D).

Fig. 1
figure 1

FLU suppresses melanoma growth in immunocompetent mice. A Schedule of mouse treatments. B C57BL/6J mice bearing s.c. B16F10 tumors were administered with FLU or vehicle (200 mg/kg or 20 ml/kg, respectively, i.p.) once daily on a 5 days-on/2 days-off schedule. Treatment commenced on Day 7. Tissues were collected for further analysis on Day 17. Data represent mean ± SEM. n = 6 mice/group. Statistical significance was assessed by t test. ***p < 0.001. C Tumor growth rates were determined by calculating the natural log (ln) of individual tumor volumes on Days 15–17 and calculating slopes in MS Excel. Data represent mean ± SEM. n = 5–6 mice/group. Statistical significance was assessed by Mann Whitney test. D Kaplan–Meier survival analysis was performed with a tumor size limit of 500 mm3 as previously described [51]. n = 6 mice/group. Statistical significance was assessed by log-rank (Mantel-Cox) test

FLU inhibits PD-1 levels without affecting levels of PD-L1 and stimulates CD3 + cell accumulation

Immunohistochemical staining of the tumors revealed that systemic FLU treatment suppressed PD-1 expression, both peritumorally (Fig. 2A) and intratumorally (Fig. 2B). FLU also caused the dramatic accumulation of CD3+ cells as compared with animals treated with vehicle (Fig. 2C). In contrast, FLU treatment did not affect levels of PD-L1 (Fig. 2D).

Fig. 2
figure 2figure 2

FLU inhibits PD-1 levels without affecting PD-L1 and stimulates CD3+ cell accumulation. Immunohistochemical analysis was performed on Day17 B16F10 tumors from vehicle or FLU-treated mice. Quantification and representative photomicrographs of immunostaining for A peritumoral, B intratumoral PD-1+, C peritumoral skin CD3+ and D intratumoral PD-L1+. IOD and tissue area were assessed using Image-Pro Plus® and mean of IOD/µm2 was determined. Representative immunohistochemical staining is shown (right). Data was analyzed by t-test or Mann Whitney test (n = 5–6 mice/group) (left). Error bars represent SEM. S denotes skin, PT denotes peritumoral area, T denotes tumor, IT denotes intratumoral. Arrows indicate positive staining. Dashed lines separate PT and T

FLU inhibits PD-1 protein expression in Jurkat T cells

Western blotting with extracts of Jurkat cells, a model human T cell line, revealed that FLU reduced PD-1 protein expression within 24 h (Fig. 3A). These findings were supported by separate experiments using flow cytometry (Fig. 3B).

Fig. 3
figure 3

FLU inhibits PD-1 expression in Jurkat T cells. A Western blotting was performed with extracts of Jurkat T cells treated with FLU (10 µM) or vehicle for 24 h. Membranes were treated with antibodies to PD-1 or β-actin then with secondary antibodies for detection. B Alternatively, cells treated with FLU (10 µM) or vehicle for 24 h were tested for PD-1 immunopositivity by flow cytometry. Data is representative of 3 biologically-independent experiments. Error bars represent SEM. Statistical significance was assessed by t test

RNA-seq reveals that FLU modulates gene expression

To gain insights into FLU’s mechanism of action, we exposed Jurkat T cells to FLU or vehicle for 4 h prior to harvest and bulk next generation RNA sequencing (RNA-seq). We chose this early time point to help identify regulatory factors that regulate PD-1 expression controlled by FLU. Multidimensional scaling (MDS) (Fig. 4A) showed clear separation between treatment groups (3 biological replicates per group). RNA-seq revealed that from a pool of 14,475 genes (12,259 protein coding, 2216 non-protein coding) (Additional file 1: Table S1) there were 1218 DEG (Additional file 2: Table S2); 687 genes increased expression and 531 reduced expression by FLU with an adjusted p value of < 0.05 (Additional file 2: Table S2). Surprisingly, RNA-seq revealed that at the 4 h timepoint, FLU increased pdcd-1 (PD-1) mRNA levels by 2.3 fold (Additional file 2: Table S2). Among genes with the largest fold increase in expression with FLU compared to vehicle were galectin 12 (LGALS12, 51.1 fold), serine protease 56 (PRSS56, 25.6 fold), endothelin converting enzyme like 1 (ECEL1, 20.1 fold), formin 1 (FMN1, 18.5 fold) and ecto-nucleotide pyrophosphatase/phosphodiesterase (ENPP1, 14.3 fold) (Additional file 2: Table S2, Fig. 4B, C). On the other hand, genes with the largest fold reduction in expression with FLU were interleukin 26 (IL26, 0.22 fold), c–c chemokine receptor type 2 (CCR2, 0.23 fold), zinc finger and BTB domain containing 16 (ZBTB16, 0.28 fold) and tripartite motif containing 67 (TRIM67, 0.31 fold) (Additional file 2: Table S2, Fig. 4B, C).

Fig. 4
figure 4

FLU modulates gene expression in Jurkat T cells. RNA-seq was performed with total RNA prepared from Jurkat T cells treated with 10 µM FLU for 4 h. A Multidimensional scaling (MDS) analysis. Con1, Con2, Con3 and Treat1, Treat2, Treat3 refer to vehicle and FLU treatment, respectively, where each is a biologically-independent replicate. B Heat map showing the top 100 DEG in our analysis. Heatmap was constructed based on z-scores of normalized expression values for the DEG and display expression levels between the biologically-independent samples. Arrow indicates JUN. C Volcano plot showing representative genes that down-regulated (left) and up-regulated (right) by FLU

Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was performed to determine whether DEGs following FLU treatment were over-represented in specific pathways. The up-regulated genes were over-represented in 25 pathways (p < 0.01), the top hit was amoebiasis, not surprising for an anthelmintic. FLU also up-regulated genes in pathways associated with hematopoeisis and PI3K-Akt signaling (Additional file 3: Table S3A and Fig. 5). Pathways associated with cancer, MAPK signaling and cytokine-cytokine receptor interactions were among those which had a significant over-representation in both up-regulated and down-regulated genes (Additional file 3: Tables S3A, B and Fig. 5).

Fig. 5
figure 5

Heat map representation of genes significantly up- or down-regulated by FLU in Jurkat T cells. These genes are associated with the A MAPK pathway, B PI3K-Akt signaling pathway, C pathways in cancer, D JAK-STAT signaling pathway, and E cytokine-cytokine receptor interactions. RNA-seq was performed with total RNA prepared from Jurkat T cells treated with 10 µM FLU for 4 h. Heatmaps were constructed based on z-scores of normalized expression values for the DEG and display expression levels between the biologically-independent samples. Each row is scaled using the scale function of the heat map tool which takes the z-score across each row

Given FLU’s effects on PD-1 expression in Jurkat T cells and its ability to induce CD3+ T cell accumulation in tumors whose growth it inhibits, we performed GSEA on gene sets in the C7 collection representing cell states and perturbations within the immune system. This identified a large number of immunologic signature gene sets that correlated with FLU treatment (605 sets FDR < 0.25) (Additional file 4: Table S4), including genes associated with T cell function, differentiation and proliferation (Fig. 6A–C). GSEA performed with the C2 collection identified 210 gene sets correlated with FLU treatment, including those associated with cell metabolism and cell signaling (Additional file 5: Table S5 & Fig. 6D) providing further biological insight into the potential actions of FLU.

Fig. 6
figure 6

Representative GSEA of gene sets enriched by FLU in Jurkat T cells. Representative C7 gene sets enriched by FLU were those associated with T cell A function, B differentiation and C proliferation. Representative C2 gene sets enriched by FLU were those associated with cell metabolism and cell signaling (D)

FLU stimulates JUN and FOS expression

FLU caused after 4 h a 2.2 fold increase in mRNA levels of the transcription factor, JUN (Additional file 2: Table S2, Fig. 5A, C), a member of the basic leucine zipper (bZIP) transcription factor family. Interestingly, recent studies show that overexpression of JUN in T cells can reinvigorate T cells and improve anti-tumor potency [24]. JUND and JUNB expression was modulated 1.61 and 0.67 fold, respectively, by FLU, while another immediate-early gene EGR1 [25,26,27,28] was modulated 0.45 fold (Additional file 2: Table S2). FLU induction of JUN at the level of protein was confirmed by Western blotting which showed that JUN levels were induced within 2 h and peaked at 6 h (Fig. 7A). This, to the best of our knowledge, is the first demonstration of FLU’s induction of JUN in any cell type. FLU also caused a transient induction of FOS protein, peaking at 2 h, with substantially earlier induction kinetics compared with JUN (Fig. 7A), which may explain why RNA-seq detected increased levels of JUN but not FOS (Additional file 2: Table S2, Fig. 5A, C). Others have also noted more acute and less sustained induction of FOS protein than JUN [29, 30]. These experiments also showed that FLU inhibited PD-1 protein expression within 2–4 h (Fig. 7B).

Fig. 7
figure 7

FLU activates JUN and FOS expression in Jurkat T cells and its suppression of PD-1 expression is rescued by T5224. Western blotting was performed with extracts of Jurkat T cells treated with FLU (10 µM) or vehicle for various times. Membranes were treated with antibodies to A JUN, FOS or β-actin, and B PD-1 or β-actin then with secondary antibodies for detection. C Western blotting was performed with antibodies to PD-1 (pretreated with vehicle or 10 µM T5224 for 2 h) then treated with FLU (10 µM) for a further 24 h. Densitometry was performed using NIH Image J. Error bars represent SEM. Data is representative of 3 biologically-independent experiments

Finally we determined whether FLU’s suppression of PD-1 was reliant upon AP-1. We exposed Jurkat T cells to FLU after preincubation with T5224 (3-(5-(4-(cyclopentyloxy)-2-hydroxybenzoyl)-2-((3-hydroxybenzo[d]isoxazol-6-yl)methoxy)phenyl)propanoic acid), a small molecule benzophenone-based selective inhibitor of AP-1 (FOS/JUN) binding to DNA [31, 32]. FLU suppressed PD-1 protein expression which was prevented by T5224 (Fig. 7C) suggesting that PD-1 inhibition by FLU relies, at least in part, upon AP-1.

Discussion

Immune checkpoint inhibitor therapy has fundamentally changed the clinical management of advanced melanoma and a range of other cancers. Monoclonal antibody blockade of the PD-1/PD-L1 system seeks to re-orientate immunoediting toward immunosurveillance and immune-mediated tumor recognition and lysis. Adjuvant use of immune checkpoint inhibitors in melanoma is now directed by clinical practice guidelines [33]. Recent research has demonstrated survival advantage from neoadjuvant-adjuvant immune checkpoint inhibitor therapy in melanoma, suggesting benefit in supporting the immune system to recognise the tumor prior to excision [34]. Notwithstanding this, ablative surgical intervention remains the primary treatment modality for local or locoregional melanoma given that objective response rates and response durability remain unsatisfactory for most patients. Toxicity remains a concern and dosing require hospital attendance. Repurposing existing therapeutic agents may improve response, durability and offer additional “concurrent” therapies; whilst mitigating adverse events and providing for an “at home” per oral regime. Given the potential for these repurposed agents to target signaling upstream of PD-1/PD-L1, a non-surgical approach to melanoma may be realised.

This study provides the first demonstration that the anti-parasitic drug FLU inhibits melanoma growth in immunocompetent mice with PD-1 inhibition through AP-1. RNA-seq analysis revealed that hundreds of genes were differentially expressed in a human T cell line exposed to FLU. One of these genes was the prototypic AP-1 family member, JUN. T5224, a small molecule that blocks AP-1 binding to DNA [31, 32], prevented FLU suppression of PD-1 expression. That FLU inhibits both tumor growth and PD-1 protein expression in immunocompetent mice suggests the potential use of this anthelmintic as a repurposed, small molecule inhibitor for melanoma and possibly other cancer types.

Previous studies have determined that FLU can inhibit breast cancer, gastrointestinal cancer, melanoma, neuroblastoma, oral carcinoma, leukemia and myeloma in mice [18]. However, to the best of our knowledge, the present study is the first to show that FLU can inhibit melanoma growth with suppression of PD-1 in immunocompetent mice and builds upon our earlier demonstration in immunodeficient mice [14]. Our previous study [14] showed that FLU can abrogate growth and metastasis of human melanoma grown as xenografts in SCID mice. FLU inhibited levels of phospho-(Tyr705) STAT3, PD-1 and myeloid-derived suppressor cell (MDSC) accumulation within the tumors [14]. On the other hand, the present study investigated the effect of FLU on murine melanoma growth in immunocompetent mice, we found FLU inhibited tumor growth and suppressed PD-1 levels in the tumor microenvironment. While we did not investigate the involvement of MDSC or macrophages, we found that FLU increased immunostaining for the pan T cell marker CD3+ and inhibited PD-1 expression in Jurkat cells, a model human T cell line. Moreover, unlike our previous work [14], this study incorporates the results of differential gene expression and bioinformatics analysis which provide new insights on transcriptional changes effected by FLU.

Suppression of peritumoral and intratumoral PD-1 by FLU was accompanied by increased CD3+ T cell accumulation in peritumoral skin. Park et al.also observed changes in CD3+ T cell populations in peritumoral skin [35], while Haywood et al. and Halse et al. reported changes in peritumoral PD-1 expression [36, 37]. FLU inhibited tumor growth and PD-1 levels in this study without affecting levels of PD-L1. KEGG and GSEA provided novel insight into the potential functions of FLU. For example, FLU up-regulated genes associated with PI3K-Akt signaling, a pathway linked with the development of T cell effector function [38] and reduced regulatory T cell phenotype and metabolic state [39]. This is in line with GSEA that identified a range of immunologic signature gene sets correlated with FLU treatment such as those underpinning T cell function, differentiation and proliferation. Future studies should determine the direct effects of FLU on T cell state and function.

PD-1 expression is regulated at the level of transcription by factors such as signal transducer and activator of transcription (STAT) and members of the AP-1 family [40,41,42]. AP-1 is induced by TCR signalling, cooperates with NFAT in T cells and drives IL-2 expression and effector function. The present study demonstrates that FLU increases JUN expression and that AP-1 is needed for FLU’s suppression of PD-1. T5224 has been used previously to demonstrate an important role for AP-1 in prolyl-isomerase expression in Jurkat T cells [43] and in pancreatic beta-like cells where it rescued mutant CDKAL1-associated beta cell defects [44]. Effects on PD-1 levels using T5224 and FLU are supported by recent findings by Carnevale et al. showing that increased AP-1 activity prevents T cell exhaustion and/or T cell anergy in Jurkat T cells [45]. Our findings showing differences in pdcd-1 mRNA and PD-1 protein expression levels at the 4 h timepoint suggest complex regulation. This could arise from altered pdcd-1 transcription or mRNA stability as a compensatory response to reduced levels of PD-1 protein in cells exposed to FLU. For example, PD-1 can undergo dynamic ubiquitination and proteasomal degradation in T cells [46]. Conversely, inducible FOS protein but not mRNA levels at this timepoint could reflect shorter lived FOS mRNA in cells exposed to FLU.

Repurposing FLU as a compound that can stimulate anti-cancer immunity is attractive since the benzimidazole has established safety profiles in humans and animals. FLU has been used for the treatment of gastrointestinal nematode infections in humans and animals since the 1980s. The dose of FLU for use as an anthelmintic in humans (oral 100 mg dose twice a day for ascariasis, trichuriasis, ancylostomiasis, mixed infections) is substantially below the dose used in this cancer study (200 mg/kg), which is that we used previously in immunodeficient mice [14]. Spagnuolo et al. administered FLU at 200 mg/kg daily i.p. for 2 weeks and found that it did not induce neuropathy [47]. FLU is a compound with an extremely high LD50 (> 5000 mg/kg) in mice, rats and guinea pigs [48, 49] and is metabolized. In human liver, FLU is reduced in cytosol by carbonyl reduction and NADPH is the preferred coenzyme [50]. This study provides the first evidence that FLU can inhibit tumor growth with reduced PD-1 expression in immunocompetent mice. Our study suggests the potential use of FLU as a repurposed small molecule inhibitor in melanoma and potentially other cancers involving PD-1.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article.

References

  1. Pauken KE, Torchia JA, Chaudhri A, Sharpe AH, Freeman GJ. Emerging concepts in PD-1 checkpoint biology. Semin Immunol. 2021;52:101480.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Larkin J, Chiarion-Sileni V, Gonzalez R, Grob JJ, Cowey CL, Lao CD, Schadendorf D, Dummer R, Smylie M, Rutkowski P, et al. Combined nivolumab and ipilimumab or monotherapy in untreated melanoma. N Engl J Med. 2015;373:23–34.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Robert C, Schachter J, Long GV, Arance A, Grob JJ, Mortier L, Daud A, Carlino MS, McNeil C, Lotem M, et al. Pembrolizumab versus ipilimumab in advanced melanoma. N Engl J Med. 2015;372:2521–32.

    Article  CAS  PubMed  Google Scholar 

  4. Wolchok JD, Chiarion-Sileni V, Gonzalez R, Rutkowski P, Grob JJ, Cowey CL, Lao CD, Wagstaff J, Schadendorf D, Ferrucci PF, et al. Overall survival with combined nivolumab and ipilimumab in advanced melanoma. N Engl J Med. 2017;377:1345–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Callahan MK, Kluger H, Postow MA, Segal NH, Lesokhin A, Atkins MB, Kirkwood JM, Krishnan S, Bhore R, Horak C, et al. Nivolumab plus ipilimumab in patients with advanced melanoma: updated survival, response, and safety data in a phase I dose-escalation study. J Clin Oncol. 2018;36:391–8.

    Article  CAS  PubMed  Google Scholar 

  6. Zhong L, Li Y, Xiong L, Wang W, Wu M, Yuan T, Yang W, Tian C, Miao Z, Wang T, Yang S. Small molecules in targeted cancer therapy: advances, challenges, and future perspectives. Signal Transduct Target Ther. 2021;6:201.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Li YY, Jones SJ. Drug repositioning for personalized medicine. Genome Med. 2012;4:27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Wouters OJ, McKee M, Luyten J. Estimated research and development investment needed to bring a new medicine to market, 2009–2018. JAMA. 2020;323:844–53.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Ghofrani HA, Osterloh IH, Grimminger F. Sildenafil: from angina to erectile dysfunction to pulmonary hypertension and beyond. Nat Rev Drug Discov. 2006;5:689–702.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Chakraborty C, Sharma AR, Bhattacharya M, Agoramoorthy G, Lee SS. The drug repurposing for COVID-19 clinical trials provide very effective therapeutic combinations: lessons learned from major clinical studies. Front Pharmacol. 2021;12:704205.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Ishida Y, Agata Y, Shibahara K, Honjo T. Induced expression of PD-1, a novel member of the immunoglobulin gene superfamily, upon programmed cell death. EMBO J. 1992;11:3887–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Sharpe AH, Pauken KE. The diverse functions of the PD1 inhibitory pathway. Nat Rev Immunol. 2018;18:153–67.

    Article  CAS  PubMed  Google Scholar 

  13. Kleffel S, Posch C, Barthel SR, Mueller H, Schlapbach C, Guenova E, Elco CP, Lee N, Juneja VR, Zhan Q, et al. Melanoma cell-intrinsic PD-1 receptor functions promote tumor growth. Cell. 2015;162:1242–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Li Y, Acharya G, Elahy M, Xin H, Khachigian LM. The anthelmintic flubendazole blocks human melanoma growth and metastasis and suppresses programmed cell death protein-1 and myeloid-derived suppressor cell accumulation. Cancer Lett. 2019;459:268–76.

    Article  CAS  PubMed  Google Scholar 

  15. Freeman GJ, Long AJ, Iwai Y, Bourque K, Chernova T, Nishimura H, Fitz LJ, Malenkovich N, Okazaki T, Byrne MC, et al. Engagement of the PD-1 immunoinhibitory receptor by a novel B7 family member leads to negative regulation of lymphocyte activation. J Exp Med. 2000;192:1027–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Latchman Y, Wood CR, Chernova T, Chaudhary D, Borde M, Chernova I, Iwai Y, Long AJ, Brown JA, Nunes R, et al. PD-L2 is a second ligand for PD-1 and inhibits T cell activation. Nat Immunol. 2001;2:261–8.

    Article  CAS  PubMed  Google Scholar 

  17. Wang T, Cai S, Cheng Y, Zhang W, Wang M, Sun H, Guo B, Li Z, Xiao Y, Jiang S. Discovery of small-molecule inhibitors of the PD-1/PD-L1 axis that promote PD-L1 internalization and degradation. J Med Chem. 2022;65:3879–93.

    Article  CAS  PubMed  Google Scholar 

  18. Khachigian LM. Emerging insights on functions of the anthelmintic flubendazole as a repurposed anticancer agent. Cancer Lett. 2021;522:57–62.

    Article  CAS  PubMed  Google Scholar 

  19. Rudolf K, Rudolf E. An analysis of mitotic catastrophe induced cell responses in melanoma cells exposed to flubendazole. Toxicol In Vitro. 2020;68:104930.

    Article  CAS  PubMed  Google Scholar 

  20. Hou ZJ, Luo X, Zhang W, Peng F, Cui B, Wu SJ, Zheng FM, Xu J, Xu LZ, Long ZJ, et al. Flubendazole, FDA-approved anthelmintic, targets breast cancer stem-like cells. Oncotarget. 2015;6:6326–40.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Zhou X, Zou L, Chen W, Yang T, Luo J, Wu K, Shu F, Tan X, Yang Y, Cen S, et al. Flubendazole, FDA-approved anthelmintic, elicits valid antitumor effects by targeting P53 and promoting ferroptosis in castration-resistant prostate cancer. Pharmacol Res. 2021;164:105305.

    Article  CAS  PubMed  Google Scholar 

  22. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Robinson MD, Oshlack A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010;11:R25.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Lynn RC, Weber EW, Sotillo E, Gennert D, Xu P, Good Z, Anbunathan H, Lattin J, Jones R, Tieu V, et al. c-Jun overexpression in CAR T cells induces exhaustion resistance. Nature. 2019;576:293–300.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Khachigian LM. Early growth response-1, an integrative sensor in cardiovascular and inflammatory disease. J Am Heart Assoc. 2021;10:e023539.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Li D, Ilnytskyy Y, Kovalchuk A, Khachigian LM, Bronson RT, Wang B, Kovalchuk O. Crucial role for early growth response-1 in the transcriptional regulation of miR-20b in breast cancer. Oncotarget. 2013;4:1373–87.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Day FL, Rafty LA, Chesterman CN, Khachigian LM. Angiotensin II (ATII)-inducible platelet-derived growth factor A-chain gene expression is p42/44 extracellular signal-regulated kinase-1/2 and Egr-1-dependent and mediated via the ATII type 1 but not type 2 receptor: induction by ATII antagonized by nitric oxide. J Biol Chem. 1999;274:23726–33.

    Article  CAS  PubMed  Google Scholar 

  28. Santiago FS, Atkins DG, Khachigian LM. Vascular smooth muscle cell proliferation and regrowth after mechanical injury in vitro are Egr-1/NGFI-A-dependent. Am J Pathol. 1999;155:897–905.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Kick G, Messer G, Plewig G, Kind P, Goetz AE. Strong and prolonged induction of c-Jun and c-Fos proto-oncogenes by photodynamic therapy. Br J Cancer. 1996;74:30–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Candeliere GA, Prud’homme J, St-Arnaud R. Differential stimulation of Fos and Jun family members by calcitriol in osteoblastic cells. Mol Endocrinol. 1991;5:1780–8.

    Article  CAS  PubMed  Google Scholar 

  31. Aikawa Y, Morimoto K, Yamamoto T, Chaki H, Hashiramoto A, Narita H, Hirono S, Shiozawa S. Treatment of arthritis with a selective inhibitor of c-Fos/activator protein-1. Nat Biotechnol. 2008;26:817–23.

    Article  CAS  PubMed  Google Scholar 

  32. Kumar A, Bora U. In silico inhibition studies of Jun-Fos-DNA complex formation by curcumin derivatives. Int J Med Chem. 2012;2012:316972.

    PubMed  PubMed Central  Google Scholar 

  33. Long GV, Luke JJ, Khattak MA, de la Cruz Merino L, Del Vecchio M, Rutkowski P, Spagnolo F, Mackiewicz J, Chiarion-Sileni V, Kirkwood JM, et al. Pembrolizumab versus placebo as adjuvant therapy in resected stage IIB or IIC melanoma (KEYNOTE-716): distant metastasis-free survival results of a multicentre, double-blind, randomised, phase 3 trial. Lancet Oncol. 2022;23:1378–88.

    Article  CAS  PubMed  Google Scholar 

  34. Patel SP, Othus M, Chen Y, Wright GP Jr, Yost KJ, Hyngstrom JR, Hu-Lieskovan S, Lao CD, Fecher LA, Truong TG, et al. Neoadjuvant-adjuvant or adjuvant-only pembrolizumab in advanced melanoma. N Engl J Med. 2023;388:813–23.

    Article  CAS  PubMed  Google Scholar 

  35. Park SL, Buzzai A, Rautela J, Hor JL, Hochheiser K, Effern M, McBain N, Wagner T, Edwards J, McConville R, et al. Tissue-resident memory CD8(+) T cells promote melanoma-immune equilibrium in skin. Nature. 2019;565:366–71.

    Article  CAS  PubMed  Google Scholar 

  36. Haywood S, Garioch J, Ramaiya A, Moncrieff M. Quantitative and spatial analysis of CD8+/PD-1 tumor-infiltrating lymphocytes as a predictive biomarker for clinical response of melanoma in-transit metastases to topical immunotherapy. Ann Surg Oncol. 2021;28:1029–38.

    Article  PubMed  Google Scholar 

  37. Halse H, Colebatch AJ, Petrone P, Henderson MA, Mills JK, Snow H, Westwood JA, Sandhu S, Raleigh JM, Behren A, et al. Multiplex immunohistochemistry accurately defines the immune context of metastatic melanoma. Sci Rep. 2018;8:11158.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Kim EH, Suresh M. Role of PI3K/Akt signaling in memory CD8 T cell differentiation. Front Immunol. 2013;4:20.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Han JM, Patterson SJ, Levings MK. The role of the PI3K signaling pathway in CD4(+) T cell differentiation and function. Front Immunol. 2012;3:245.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Bally AP, Austin JW, Boss JM. Genetic and epigenetic regulation of PD-1 expression. J Immunol. 2016;196:2431–7.

    Article  CAS  PubMed  Google Scholar 

  41. Xiao G, Deng A, Liu H, Ge G, Liu X. Activator protein 1 suppresses antitumor T-cell function via the induction of programmed death 1. Proc Natl Acad Sci USA. 2012;109:15419–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Lu P, Youngblood BA, Austin JW, Mohammed AU, Butler R, Ahmed R, Boss JM. Blimp-1 represses CD8 T cell expression of PD-1 using a feed-forward transcriptional circuit during acute viral infection. J Exp Med. 2014;211:515–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Kumar S, Singh SK, Viswakarma N, Sondarva G, Nair RS, Sethupathi P, Sinha SC, Emmadi R, Hoskins K, Danciu O, et al. Mixed lineage kinase 3 inhibition induces T cell activation and cytotoxicity. Proc Natl Acad Sci USA. 2020;117:7961–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Zeng H, Guo M, Zhou T, Tan L, Chong CN, Zhang T, Dong X, Xiang JZ, Yu AS, Yue L, et al. An isogenic human ESC platform for functional evaluation of genome-wide-association-study-identified diabetes genes and drug discovery. Cell Stem Cell. 2016;19:326–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Carnevale J, Shifrut E, Kale N, Nyberg WA, Blaeschke F, Chen YY, Li Z, Bapat SP, Diolaiti ME, O’Leary P, et al. RASA2 ablation in T cells boosts antigen sensitivity and long-term function. Nature. 2022;609:174–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Hu X, Wang J, Chu M, Liu Y, Wang ZW, Zhu X. Emerging role of ubiquitination in the regulation of PD-1/PD-L1 in cancer immunotherapy. Mol Ther. 2021;29:908–19.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Spagnuolo PA, Hu J, Hurren R, Wang X, Gronda M, Sukhai MA, Di Meo A, Boss J, Ashali I, Beheshti Zavareh R, et al. The antihelmintic flubendazole inhibits microtubule function through a mechanism distinct from Vinca alkaloids and displays preclinical activity in leukemia and myeloma. Blood. 2010;115:4824–33.

    Article  CAS  PubMed  Google Scholar 

  48. Fuchs R. Flubendazole, vol 31. Geneva, Switzerland: World Health Organization (WHO); 1993.

    Google Scholar 

  49. FLUBENDAZOLE_EMEA/CVMP/33128/2006-FINAL. http://www.ema.europa.eu/docs/en_GB/document_library/Maximum_Residue_Limits_-_Report/2009/11/WC500014292.pdf.

  50. Kubicek V, Skalova L, Skarka A, Kralova V, Holubova J, Stepankova J, Subrt Z, Szotakova B. Carbonyl reduction of flubendazole in the human liver: strict stereospecificity, sex difference, low risk of drug interactions. Front Pharmacol. 2019;10:600.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Haynes J, McKee TD, Haller A, Yadong Wang Y, Leung C, Gendoo DMA, Lima-Fernandes E, Kreso A, Wolman R, Szentgyorgyi E, et al. Administration of hypoxia-activated prodrug evofosfamide after conventional adjuvant therapy enhances therapeutic outcome and targets cancer-initiating cells in preclinical models of colorectal cancer. Clin Cancer Res. 2018;24:2116–27.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank A/Prof Anthony M. Joshua (Garvan Institute of Medical Research, Sydney) for critical review of the manuscript and Prof Christopher R. Parish (John Curtin School of Medical Research, Australian National University, Canberra) for helpful comments. Data presented in this work was acquired by personnel and/or instruments at the Mark Wainwright Analytical Centre (MWAC) of UNSW Sydney, which is in part funded by the Research Infrastructure Scheme at UNSW.

Funding

This work was supported by grants from the National Health and Medical Research Council of Australia (LMK) and from the NSW State Government RAAP scheme (SC and MRW).

Author information

Authors and Affiliations

Authors

Contributions

LMK wrote the first draft. All authors provided intellectual input and/or data and contributed to the writing of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Levon M. Khachigian.

Ethics declarations

Ethics approval and consent to participate

The protocol was approved by the UNSW Animal Care and Ethics Committee.

Consent for publication

N/A.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Table S1. Pool of 14475 genes in bioinformatics analysis of RNA-seq data showing 12259 protein coding and 2216 non-protein coding genes. RNA-seq was performed with total RNA prepared from Jurkat T cells treated with 10 µM FLU for 4 h

Additional file 2

: Table S2. Complete list of DEG following RNA-seq and bioinformatics analysis. RNA-seq was performed with total RNA prepared from Jurkat T cells treated with 10 µM FLU for 4 h.

Additional file 3

: Table S3. KEGG pathway analysis of RNA-seq data. RNA-seq was performed with total RNA prepared from Jurkat T cells treated with 10 µM FLU for 4 h. Over-representation of DEG following FLU treatment in specific pathways that were (A) up-regulated and (B) down-regulated

Additional file 4

: Table S4. C7 gene sets enriched by FLU. RNA-seq was performed with total RNA prepared from Jurkat T cells treated with 10 µM FLU for 4 h. There were 3738 of 5082 C7 gene sets upregulated in the FLU group, 605 sets significant at FDR < 0.25, and 1344 of the 5082 gene sets correlated with the vehicle, 381 were significant at FDR < 0.25.

Additional file 5

: Table S5. C2 gene sets enriched by FLU. RNA-seq was performed with total RNA prepared from Jurkat T cells treated with 10 µM FLU for 4 h. There were 3058 of 4691 C2 gene sets upregulated in the FLU group, 210 sets significant at FDR < 0.25, and 1633 of the 4691 gene sets correlated with the vehicle group, 345 sets significant at FDR < 0.25.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, Y., Wu, B., Hossain, M.J. et al. Flubendazole inhibits PD-1 and suppresses melanoma growth in immunocompetent mice. J Transl Med 21, 467 (2023). https://doi.org/10.1186/s12967-023-04289-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-023-04289-y

Keywords