ORIGINAL RESEARCH article
Computational Drug Repositioning and Experimental Validation of Ivermectin in Treatment of Gastric Cancer
- 1Department of Clinical and Molecular Medicine, Norwegian University of Science and Technology (NTNU), Trondheim, Norway
- 2The Central Norway Regional Health Authority (RHA), Stjørdal, Norway
- 3Surgical Clinic, St. Olavs Hospital, Trondheim University Hospital, Trondheim, Norway
- 4Division of Digestive and Liver Diseases, Columbia University College of Physicians and Surgeons, New York, NY, United States
Objective: The aim of the present study was repositioning of ivermectin in treatment of gastric cancer (GC) by computational prediction based on gene expression profiles of human and mouse model of GC and validations with in silico, in vitro and in vivo approaches.
Methods: Computational drug repositioning was performed using connectivity map (cMap) and data/pathway mining with the Ingenuity Knowledge Base. Tissue samples of GC were collected from 16 patients and 57 mice for gene expression profiling. Additional seven independent datasets of gene expression of human GC from the TCGA database were used for validation. In silico testing was performed by constructing interaction networks of ivermectin and the downstream effects in targeted signaling pathways. In vitro testing was carried out in human GC cell lines (MKN74 and KATO-III). In vivo testing was performed in a transgenic mouse model of GC (INS-GAS mice).
Results: GC gene expression “signature” and data/pathway mining but not cMAP revealed nine molecular targets of ivermectin in both human and mouse GC associated with WNT/β-catenin signaling as well as cell proliferation pathways. In silico inhibition of the targets of ivermectin and concomitant activation of ivermectin led to the inhibition of WNT/β-catenin signaling pathway in “dose-depended” manner. In vitro, ivermectin inhibited cell proliferation in time- and concentration-depended manners, and cells were arrested in the G1 phase at IC50 and shifted to S phase arrest at >IC50. In vivo, ivermectin reduced the tumor size which was associated with inactivation of WNT/β-catenin signaling and cell proliferation pathways and activation of cell death signaling pathways.
Conclusion: Ivermectin could be recognized as a repositioning candidate in treatment of gastric cancer.
Drug repositioning (also called drug repurposing) is a strategy for identifying new uses for approved or investigational drugs that are outside the scope of the original medical indications. Repositioned drugs may reveal new targets and pathways that can be further exploited (Ashburn and Thor, 2004; Pushpakom et al., 2019). Advantages of drug repositioning are related to the drugs that have known mechanisms of action, pharmacological properties, such as pharmacokinetics, pharmacodynamics, posology (the appropriate doses of drugs) and toxicity (Verbaanderd et al., 2017). Of note is that both the pre-clinical and clinical safety data are available (Verbaanderd et al., 2017; Antoszczak et al., 2020). Thus, compared to traditional methods of drug development, drug repositioning requires drastically shortened development time and reduced costs while providing similar therapeutic benefits. Approaches of drug repositioning include computational methods, such as connectivity map (cMap), data mining, pathway mining, ontology modeling, in silico and biological experimental validations (e.g., in vitro and in vivo). The computational drug repositioning can be conducted as repurposing with a defined purpose, repurposing with a strategy, and repurposing with confidence by utilizing reference datasets which are disease-based, drug-based, or knowledge-based (Liu et al., 2013; Subramanian et al., 2017; Xue et al., 2018).
Ivermectin was identified in late 1960s, first approved as veterinary medicine and then human medicine in 1980s for the control of parasitic infection. The discovery and development of ivermectin by William C. Campbell and Satoshi Ōmura were recognized by Nobel Prize in Physiology or Medicine in 2015 (Callaway and Cyranoski, 2015). Ivermectin acts on γ-aminobutyric acid (GABA)-gated chloride channels in the interneuronic synapses of a parasite, whereas in humans, the nerves that are sensitive to GABA, are protected by the blood/brain barrier (Sutherland and Campbell, 1990; Davis et al., 1999; Ikeda, 2003; Chen and Kubo, 2018). However, ivermectin has been repositioned as a broad-spectrum antiviral and antimicrobial agent (Andersen et al., 2020). Interestingly, it is also known that ivermectin can be widely distributed in humans because of its high lipophilicity and thus might exhibit anti-tumor activity in colorectal cancer, breast cancer, glioma, head and neck cancer, leukemia, melanoma, pancreatic cancer, and prostate cancer (Chiou et al., 1987; Juarez et al., 2018).
Gastric cancer (GC) is the fifth most common malignant disease worldwide with the third highest incidence and mortality rate among all cancers (Rawla and Barsouk, 2019). The 5 years overall survival rate is 10–30% except for Japan (50–70%) (Parkin et al., 2002; Matsuda and Saika, 2013). Gastrectomy combined with platinum-based chemotherapy is the most beneficial approach in patient care, and novel targeted therapy, including PD-1 inhibitor in first and second-line setting for advanced GC, are under development (Sitarz et al., 2018; Selim et al., 2019; GBD 2017 Stomach Cancer Collaborators, 2020). However, new drugs and drug repositioning are needed particularly in consideration of the global burden of this deadly disease.
Previously, we have showed repositioning of botulinum toxin type A (also known as botox), everolimus (RAD001) and devimistat (CPI-613) in treatment of GC (Zhao et al., 2014; Rabben et al., 2016; Rabben et al., 2021). The aim of the present study was to reposition ivermectin in treatment of GC. To this end, we have developed and/or utilized the approaches from computational drug repositioning to in silico, in vitro and in vivo validations (Figure 1).
FIGURE 1. Flow chart of study design. Computational drug repositioning was carried out by using gene expression signatures representing gastric cancer of both patients and mouse model and connectivity map (cMap) and data/pathway mining with the Ingenuity Knowledge Base. Validation included in silico, in vitro and in vivo methods of ivermectin treatment. The rationale of using human samples of GC (without ivermectin treatment) was i) to perform computational prediction and data mining and ii) to make a comparison with the animal model in order to demonstrate that the animal study could be relevant in the design of clinical trial of ivermectin in the future.
Materials and Methods
Patients and Animals
Surgical biopsies were collected from 16 patients who underwent total/subtotal or distal gastrectomy because of GC since 2012 at St. Olav’s Hospital, Trondheim, Norway (Table 1). Four biopsies per patients were taken from tumor and normal tissue and used for clinical pathological evolution and gene expression profiling. The study was approved by the Regional Committees for Medical and Health Research Ethics Central Norway (REK 2012-1029). Furthermore, seven independent datasets of human GC from the TCGA database were used (Table 2).
TABLE 2. WNT/β-catenin signaling in human gastric cancer (one in the present study and 7 STAD datasets deposited in the TCGA database).
The mouse model of GC, i.e., the transgenic INS-GAS mice which spontaneously develop gastric cancer, was used (Zhao et al., 2014). Stomachs were collected from 26 mice, i.e., six females and four males INS-GAS mice at age of 15 months and eight females and eight males wild-type (WT) mice at age of 12 months, for gene expression profiling. In addition, 31 INS-GAS mice, i.e., 12 females and nine males at age of 10 months, and 10 WT mice, i.e., 6 females and four males at age of 10 months, were used for in vivo testing.
Animals were housed as four to five mice per cage on wood chip bedding with a 12 h light/dark cycle in a specific pathogen free environment with room temperature of 22°C and 40–60% relative humidity. Animals were inspected daily by investigators and authorized veterinarian using a scoring sheet. Animals should be euthanized at score of 10 if they are emaciated, underconditioned in five consecutive days, or show poor clinical signs (e.g., body weight, appearance, and behavior) before end of the study. This was done according to the Directive 2010/63/EU in which human primary endpoint are defined as “the earliest indicator in an animal experiment of (potential) pain and/or distress that, within the context of moral justification and scientific endpoints to be met, can be used to avoid or limit pain and/or distress by taking actions such as humane killing or terminating or alleviating the pain and distress (Hendriksen and Morton, 1999)” (https://www.humane-endpoints.info/en/council-directive-2010-63-eu). The study was approved by The Norwegian Food Safety Authority (Mattilsynet).
Total RNA was extracted from the surgical biopsies of patients and harvested stomachs of mice. RNA quality and quantity were obtained using NanoDrop One (Thermo Scientific, Norway) and Agilent Bioanalyser. For human samples, RNA microarray of GC samples, including 24 tumors of intestinal, diffuse and mixed types from seven patients and 37 normal tissue from six patients, was performed using Illumina platform as described earlier (Zhao et al., 2014). Illumina microarray data was analyzed using Lumi on the log2 scale and analyzed using the empirical Bayesian method implemented in Limma. The data is accessible via Mendeley Data repository with DOI link at http://dx.doi.org/10.17632/hzmfshy7hp.1. Illumina identifiers (ILMN) were uploaded to Ingenuity Pathway Analysis (IPA, QIAGEN, Hilden, Germany) together with log2-fold change, p-values and q-values (false discovery rates). For mouse samples, RNA sequencing was performed using Illumina HiSeqNS500 instrument (NextSeq 500) at 75 bp with paired end (PE) reads using NS500H flow cells with 25M reads/sample. Paired end forward read length (R1): 81, reverse read length (R2): 81. Downstream processing and analysis of the data was performed in the Bioconductor environment in R. For humans, a total of 47,323 transcripts was assigned to analysis in which 37,489 transcripts were mapped and 9,834 transcripts unmapped by Ingenuity Pathway Analysis (IPA) (QIAGEN, Hilden, Germany). For mice, a total of 54,460 transcripts was loaded in which 54,162 were mapped/298 transcripts were unmapped in IPA. For mouse GC after ivermectin treatment, 54,416 transcripts were loaded in which all were mapped in IPA. Filtering of datasets included species (mouse or human) and p-value cut-off (p < 0.05). Gene expression was analyzed using a t-test between tumor and normal tissue in patients, between INS-GAS and WT mice and between INS-GAS mice with and without ivermectin. Genes with a p-value of less than 0.05 were considered to be differentially expressed. Transcriptomics datasets were analyzed using IPA. Molecular networks and canonical pathways were algorithmically constructed based on known connectivity and relationships among genes/proteins/metabolites using Ingenuity Knowledge Base. Local and regulatory z-scores for canonical pathways and diseases and biofunctions that overlapped with the experimental data of the present study were calculated using the formula described previously (Sitarz et al., 2018). IPA has sophisticated algorithms to calculate predicted functional activation/inhibition of canonical pathways, diseases and functions, transcription regulators and regulators based on their downstream molecule expressions (QIAGEN Inc., https://www.qiagenbioinformatics.com/products/ingenuitypathway-analysis). Fischer’s exact test was used to calculate a p-value determining the probability that the association between the genes in the datasets from human GC and mouse GC and the canonical pathway or disease/function by chance alone.
Connectivity Map and Data/Pathway Mining
The concept of a Connectivity Map (cMap) was recently developed, whereby genes, drugs, and disease states are connected by virtue of common gene expression signatures (Qu and Rajpal, 2012; Subramanian et al., 2017; Musa et al., 2018). To identify candidate drugs, the gene expression signature of GC was generated based on the gene expression profile of human GC. A positive cMap score indicates there is a positive similarity between a given perturbagen’s signature, i.e., genes that are increased by treatment (in reference datasets) are also upregulated in the human GC dataset, while a negative score indicates that the two signatures are opposing. cMap was performed using the gene expression signature of human GC (n = 7 GC vs. n = 6 normal tissue).
Data mining was performed using the gene expression profile data of 61 samples from 16 patients, 26 samples from 26 mice, and 324 samples from seven independent datasets from the TCGA database. Furthermore, knowledge-based pathway mining was used based on previous studies that showed WNT/β-catenin signaling pathway as one of the important pathways in gastric tumorigenesis (Zhao et al., 2014; Rabben et al., 2021). Custom-made molecular networks were generated using the Ingenuity Knowledge Base. Networks were then algorithmically generated based on their interrelationships. WNT/β-catenin signaling pathway was constructed based on the transcriptomic data of INS-GAS mice and were then entered into the “Pathway” module of the IPA to obtain the nodes in every corresponding signaling pathway. Nodes from pathways were added as entries into the “My list”-function and all entries in the list were added to the “My pathway” in IPA. My pathway was used to produce a network of nodes/genes from the WNT/β-catenin signaling pathway that matched with our experimental data from INS-GAS vs. WT mice. The build-tool was used to connect nodes using edges, i.e., relationships including both direct and indirect interactions like chemical-protein interactions, ubiquitination, molecular cleavage, translocation, localization, phosphorylation, expression, protein-protein interactions, activation, regulation of binding, inhibition, membership, reaction, protein-DNA interactions, transcription and modification. The Canonical Pathway overlay-tool was used to arrange the entries into clusters based on pathway. Local z-scores were calculated in IPA based on the dataset’s correlation with the activated state. Negative z-scores indicate a decrease in activity, positive z-scores indicate an increase in activity. Canonical pathways were identified using statistical cut-offs at p < 0.05.
In Silico Testing
The expression data from mouse GC was compared to all genes in the pathway. The molecule activity predictor (MAP)-function was used to predict activation/inhibition between the nodes in the network. The in silico tool integrated with the MAP-function was employed to predict effects on the network after gene inhibition and/or stimulation in the ivermectin cluster. Connections between genes were then algorithmically generated based on their interrelationships including both direct and indirect interactions like chemical-protein interactions, ubiquitination, molecular cleavage, translocation, localization, phosphorylation, expression, protein-protein interactions, activation, regulation of binding, inhibition, membership, reaction, protein-DNA interactions, transcription and modification. Network clusters of WNT/β-catenin pathway was constructed based on the transcriptomic data of INS-GAS mice (i.e., limited to and built on genes from the dataset). The build-tool was used to connect nodes using edges, i.e., relationships. Categorical values were set to each gene/node using a semi-quantitative method to quantify the color-change resulting from in silico inhibition. Local z-scores were calculated in IPA based on the dataset’s correlation with the activated state. Negative z-scores indicate a decrease in activity, positive z-scores indicate an increase in activity. Canonical pathways were identified using statistical cut-offs at p < 0.05.
In Vitro Testing
GC cell lines included human gastric cancer cells MKN74 (intestinal type) and KATO-III (diffuse type) (for detailed information on molecular characteristic, see Yokozaki, 2000). It should be noticed that both cell lines overexpress β-catenin (Asciutti et al., 2011). Cells were maintained in RPMI-1640 medium (Sigma Aldrich, Oslo, Norway) supplemented with fetal bovine serum (10%, FBS), Sodium pyruvate and penicillin streptomycin solution (1%) (Sigma Aldrich, Oslo, Norway) in a humidified incubator holding 5% CO2 and 37°C. For proliferation assay, MKN74 and KATO-III were seeded in 96-well plates (2,500 cells/well and 3,000 cells/well, respectively) and incubated overnight. Ivermectin (MW: 875.09 g/mol) was dissolved in DMSO (100%) to 50 mM stock solution. Cells were treated with ivermectin (0–50 µM) or vehicle control (0.45% v/v DMSO) for 24, 48, and 72 h. Proliferation was measured using a commercial CCK-8 Kit (Sigma Aldrich, Oslo, Norway) with absorbance read at 450 nm. For cell cycle analysis, KATO-III cells were seeded as 3.0 × 105 cells/well in 6-well plates and incubated for 72 h with medium change after 48 h. Ivermectin was added to the wells as final concentrations of 12, 15 or 18 µM for 24 h. Cells were harvested by trypsin, washed twice in room tempered PBS, resuspended in ice cold ethanol (70%) and kept at −20°C for minimum 15 min. Cells were washed twice in cold PBS and centrifuged (1,500 rpm, 5 min, 4°C), and resuspended in freshly prepared PI staining solution (0.25% Triton- X-100, 50 μg/ml propidium iodide (PI) and 200 μg/ml RNAase) for minimum 30 min. Cell cycle distribution was analyzed using FACS. Single cells were gated to exclude doublets and clustered cells. 2.0 × 104 cells were counted per sample, and percentage cell distribution was derived from obtained histograms in the FACSDiva software. Results are presented as means of n = 3 replicates/treatments. Data was analyzed using Microsoft Excel 2010.
In Vivo Testing
Thirty-one INS-GAS mice were randomly divided into two groups: ivermectin treatment (12 females and nine males at age of 10 months) and controls (no treatment, six females and four males at age of 10 months). Ivermectin was reconstituted from lyophilized powder in DMSO to 50 mM solution and then diluted in saline before use. The treatment regimen was designed to let the mice tolerate the procedure easily, i.e., intraperitoneal injection at a dose of 10 mg/kg in a volume of about 0.5 ml/mouse with 27G needle once per day for 5 days, followed by no treatment for 5 days and then injection once per day for 10 days. This regimen was repeated 10 days later. The total duration of treatment was 2 months (2 × 30 days). Vehicle treatment was not performed because neither vehicle per se nor procedure would lead to any significant stress response. The mice were euthanized under isoflurane inhalation anesthesia (2–3%), and stomachs were collected as described previously (Zhao et al., 2014). Tumor volume density (% of glandular area of the stomach occupied by tumor) was measured using a point count method (Zhao et al., 2014). The tissue samples were collected for transcriptomics as aforementioned.
Values are expressed as means ± SEM or SD (stated in individual figure legend). For comparison of two independent groups, student independent t-test was used. For comparison of multiple groups, one-way ANOVA with Tukey’s or Dunnett’s post-hoc tests were used (stated in individual figure legend). SPSS version 26.0 for Windows (SPSS Inc., Chicago, IL, United States) was used and a p-value < 0.05 was considered to be statistically significant. Other methods are stated in corresponding figure legends.
Computational Drug Repositioning Suggests the WNT/β-Catenin Signaling as Potential Target of Ivermectin
A cMap was created according to the gene expression signature (Figures 2A,B). A total of 2,428 drugs were categorized into 47 groups of inhibitors for, e.g., DNA synthesis, murine double minute (MDM) and lactate dehydrogenase, including ivermectin. Additionally, drugs we have demonstrated previously, including gemcitabine, paclitaxel, everolimus, and scopolamine, were also found (Zhao et al., 2014; Rabben et al., 2021).
FIGURE 2. Gene expression signature and connectivity map (cMAP). (A) Heatmap of human GC gene expression signature that constitutes an activation of cancer disease based on differential expression of >22,000 genes. Size of square is proportional to the number of genes contained in the specific function and color represent activity state (z-score; orange: activated, blue: decreased). (B) Connectivity map (cMap) showing associations between a large-scale compendium of functional perturbations in cancer cell lines coupled to the human GC gene expression signature based on the L1000 assay (Subramanian et al., 2017). Note: Ivermectin and other known drugs are visualized.
Data/pathway mining revealed activation of the WNT/β-catenin signaling in human as well as mouse GC (Table 2). In addition to our own human GC and mouse GC data, seven independent datasets of human stomach adenocarcinoma (STAD) were found to have activation of the WNT/β-catenin (Figure 3A and Table 2). Using the knowledge-based repositioning strategy, nine targets were identified in connection with WNT/β-catenin signaling pathway and proliferation in both human and mouse GC (Figure 3B), i.e., ATP-dependent translocase (Abcb1b), retinol-binding proteins (Rbp), ATP binding cassette subfamily B member 1 (ABCB1), ATP Binding cassette subfamily G member 2 (ABCG2), cytochrome P450 family 3 subfamily A member 4 (CYP3A4), P-glycoprotein (also known as multi-drug resistant protein, MDRP), ATP binding cassette subfamily B member 4 (ABCB4), cytokine, and P2X purinoceptor 7 (P2RX7). Each gene/protein connected to a subset of algorithmically chosen genes based on the Ingenuity Knowledge Base. These genes were collectively activating both WNT/β-catenin signaling and proliferation, resulting in locally activated z-scores (shown in orange) (Table 3).
FIGURE 3. Data/pathway mining of WNT/β-catenin signaling pathway. (A) Activation of WNT/β-catenin signaling pathway in eight datasets of human GC (including one used in the present study as indicated in yellow) and one dataset of mouse GC (in blue). (B) Hierarchical network representation showing ivermectin and drug targets with downstream signaling pathways of WNT/β-catenin and proliferation. The schema was created in IPA using the grow-tool and the Ingenuity Knowledge Base. The molecular entities (genes and proteins) as well as molecular functions and interactive networks were connected based on interrelationships identified by the Ingenuity Knowledge Base. Expression levels from mouse GC vs. WT. p < 0.05. See also Table 3 In silico testing shows ivermectin inhibits WNT/β-catenin signaling.
TABLE 3. Gene expression in networks comprised of ivermectin interactions and WNT/β-catenin signaling pathways in mouse gastric cancer (as presented in Figure 4A).
An in silico interaction network of ivermectin with WNT/β-catenin signaling pathway was constructed (Figure 4A). Inhibition of the targets of ivermectin led to the inhibition of downstream nodes in a “dose-dependent manner” (Figure 4B). It should be noticed that the inhibition of single molecules was not enough to have inhibitory or stimulatory effects on the signaling pathway.
FIGURE 4. In silico interaction network of ivermectin with WNT/β-catenin signaling pathway. (A) An interaction network of ivermectin with WNT/β-catenin signaling pathway was made using the Ingenuity Knowledge Base. (B) Two-/three-/four-/five combinations of inhibition/stimulation of nodes in the ivermectin cluster (A, upper cluster) had “dose-dependent” effect of the sum network scores of WNT/β-catenin (including HIPPO) pathway (A, lower cluster). Categorical values were set to each gene/node using a semi-quantitative method to quantify the color-change resulting from in silico inhibition. Dark blue colored nodes were represented by − 2, light blue as − 1, white as 0, light orange as +1 and dark orange as +2. Values are represented of means of n = 3–8 experiments per node/gene. Multiple comparison using one-way ANOVA with Tukey post hoc test was used. Created in IPA (QIAGEN).
Ivermectin Inhibits Cell Proliferation and Induces Cell Cycle Arrest
Testing of ivermectin in MKN74 and KATO-III cells showed that there were time- and concentration-depended inhibitions of proliferation by the drug with similar IC50 values for the periods of 24, 48 and 72 h (Figures 5A,B). Furthermore, ivermectin induced cell cycle arrest in a concentration-depended manner (Figure 5C). It should be noticed that ivermectin at IC50 did not affect the cells in S-phase but increased percentage of cells in G1 while reducing percentage of cells in G2/M phases. By contrast, higher concentration of ivermectin increased percentage of cells in S phase while reducing the percentage of cells in G1 and G2/M phases, suggesting that ivermectin arrested cells at the G1 phase at IC50 and higher dose of the drug shifted cells to S phase.
FIGURE 5. Time- and concentration-dependent effects of ivermectin on proliferation and cell cycle. (A,B) Concentration-response curves upon 24–72 h ivermectin-treatment in MKN74 and KATO-III cells. Proliferation was assessed using CCK-8 Kit at 450 nm. IC50 values were calculated from sigmoidal regression curve fitting using variable slope on normalized response from log10-transformed x-values (GraphPad Prism v.6). (C) Cell cycle analysis of KATO-III cells. Means of n = 3 replicates/treatment with SD. One-way ANOVA with Dunnett’s post hoc test (2-sided) compared to respective control groups.
Ivermectin Reduces Tumor Size Which Was Associated With Inactivation of WNT/β-Catenin Signaling, Down Regulation of Cell Proliferation and Upregulation of Cell Death Signaling Networks
A treatment regimen using ivermectin at 10 mg/kg for 2 months was established based on the in silico and pilot experiments. Mice tolerated the treatment well, although some mice had weight loss during treatment (<15%, p > 0.05, two-tailed). The mice had no serious side effects of ivermectin and no mice that were treated with ivermectin were killed according to the human primary endpoints which include stressful behavior, abdominal pain and impaired physical activity. The tumor size was reduced by ivermectin treatment (Figure 6A). Comparison analysis between mouse GC with and without ivermectin treatment revealed 4,112 differentially expressed genes (Figure 6B). The genes involved in WNT/β-catenin signaling pathway were particularly inhibited by ivermectin treatment, as shown by a change in z-scores from 1.151 (mouse GC without treatment) down to −1.789 (mouse GC after ivermectin treatment) (Figure 6C and Table 2) and log2 fold-changes (Figure 6D vs. 6E).
FIGURE 6. Tumor size and gene expression profiles in response to ivermectin. (A) Tumor size (% of glandular area of stomach occupied by tumor) in age-matched controls (AMC, n = 14) and ivermectin-treated mice (n = 17) (Ivermectin). Independent t-test (2-sided) between group means (normality assumption met). Error bars represent SEM. (B) Global gene expression profile of mouse GC with and without ivermectin treatment (created in RStudio using heatmap.2 function). Only differentially expressed genes with p < 0.05 are included (4,112 genes). (C) WNT/β-catenin pathway was activated (z-score = 1.151) in mouse GC without treatment but inhibited in mouse GC with ivermectin treatment (z-score = −1.789). (D,E) WNT/β-catenin gene expressions in mouse GC mice without treatment (D) and with ivermectin treatment (E). Note: same orders of individual genes in (D) and (E).
Expression analysis in IPA revealed that cell proliferation was activated in mouse GC without treatment and inactivated in mouse GC with treatment. On the other hand, cell death including apoptosis was inactivated in mouse GC without treatment but activated in mouse GC with treatment (Figures 7A–D).
FIGURE 7. Representative networks of cell proliferation and cell death in mice with or without ivermectin treatment. (A) Mouse GC cell proliferation network was created using BioProfiler in IPA. Only genes differentially expressed (p < 0.05, Log2FC ± 1.0) between mouse GC vs. WT were included. Overexpression of genes in the network results in activation of cell proliferation (orange; z-score > 0). (B) Genes in the cell proliferation network was downregulated in GC mice with ivermectin treatment (blue; z-score < 0). (C) Mouse GC cell death network including apoptotic markers was created using BioProfiler in IPA. Only genes differentially expressed (p < 0.05, Log2FC ± 1.0) between mouse GC vs. WT were included. Overexpression of genes in the network results in inhibition of cell death (blue; z-score < 0). (D) Cell death network was aberrantly expressed in GC mice with ivermectin treatment (orange; z-score > 0).
The next generation connectivity map (cMap) has been recently developed and should be acknowledged that the cMap methods and data are available without restriction to the research community (Subramanian et al., 2017). As pointed out in the original paper, a future comprehensive cMap might expand in multiple dimensions, e.g., new cell types, patient-derived induced pluripotent stem cells and genome-edited isogenic cell lines (Subramanian et al., 2017). Using this method, we found that the scores of the known drugs in treatment of GC (including ivermectin) were too low to indicate strong associations between these drugs and human GC gene expression signature, which was most likely due to the fact that the reference profile catalogue of cMap has been built to date on 12,328 genes of various cancer cell lines (including AGS which is a moderately differentiated human gastric adenocarcinoma hyperdiploid cell line) but not tumor tissues (https://clue.io/connectopedia/l1000_gene_space and https://clue.io/connectopedia/core_cmap_cell_panel).
In addition to the hypothesis generation approach by cMap, we further utilized data mining and pathway mining of knowledge-based datasets to identify the potential drugs in connection with a broad concept ranging from molecular entities (such as genes and proteins) to biological phenomena (such as molecular functions, pathways and phenotypes). Based on a better understanding of GC biology and signaling pathways, in the present study we focused on the WNT/β-catenin pathway by utilizing the algorithms of IPA which is built on a comprehensive, manually curated content of the QIAGEN Knowledge Base (over 57,000 publicly available datasets and continuously updated).
The results of the present study showed that the potential molecular targets and the mechanisms of action of ivermectin in GC differed from those in parasites in which ivermectin causes an influx of Cl-ions through the cell membrane of invertebrates by activation of specific ivermectin-sensitive ion channels (Laing et al., 2017; Chen and Kubo, 2018). In the present study, we identified ivermectin in connection with cell proliferation, particularly towards the genes (e.g., members of the adenosine triphosphate (ATP)-binding cassette (ABC) transporters). ABC are a superfamily of membrane proteins which play significant roles in transporting various exogenous and endogenous substances across membranes against concentration gradients through ATP hydrolysis, and many of these transporters are known as multidrug resistance proteins (MRPs) (Mao et al., 2019). As showed in the present study, ivermectin also acted on the ABC and the signaling pathways, leading to inhibition of cell proliferation by deactivating LXR/RXR signaling (Beltowski, 2011; Saito-Hakoda et al., 2015).
It has been shown that activation of the WNT/β-catenin signaling pathway plays a pivotal role in many types of cancer (Clevers, 2006; Zhan et al., 2017). Previously, we and other research groups have demonstrated that the tumorigenesis of gastric cancer involves the WNT/β-catenin signaling pathway and the inhibition of the signaling pathways by means of denervation can suppress the tumorigenesis (Zhao et al., 2014; Chiurillo, 2015; Koushyar et al., 2020; Rabben et al., 2021). In the present study, we applied in silico modelling to show that ivermectin could inhibit the WNT/β-catenin signaling pathway including HIPPO signaling pathway, which is known to interact each other (Hayakawa et al., 2017; Li et al., 2019). We then employed in vitro and in vivo approaches to show that ivermectin could inhibit cell proliferation and reduce tumor size, which was associated with the inhibition of the WNT/β-catenin signaling pathway. Thus, we may suggest that ivermectin could target the WNT/β-catenin singling pathway, leading to a reduced tumorigenesis. This was also in line with possible antitumor actions of ivermectin in other types of cancer cells, such as breast, colon, lung, prostate and bladder (Melotti et al., 2014; Diao et al., 2019; Nappi et al., 2020).
Control of cell proliferation generally occurs during the G1 phase and multiple signals, ranging from growth factors to DNA damage to developmental cues, influence the decision to enter S phase, when DNA is replicated (Duronio and Xiong, 2013). The results of the present study showed that ivermectin altered cell cycle in a concentration-dependent manner, which is consistent with a previous report showing accumulation of cells in the G1/S phases (Zhang et al., 2019). In the present study, IC50-dose of ivermectin caused cell cycle arrest at G1 phase, whereas at higher doses, it caused S phase arrest. It has been suggested that WNT/β-catenin activation triggered cells in S phase, and HIPPO signaling might involve in G1 phase (Benham-Pyle et al., 2016; Kim et al., 2019). The evidence of possible link between the cell cycle arrest and inhibition of WNT/β-catenin and/or HIPPO singling pathways is needed to be further investigated, particularly in the context of ivermectin for GC.
There were several limitations of the present study. The cell proliferation and apoptosis in the in vitro experiment were not evaluated further by flow cytometry nor specific assays, e.g., annexin V staining or caspase activity. However, the gene expression profiling confirmed the association between the activities of networks of cell proliferation and cell death in mice, namely increased in cell proliferation and decrease in cell death in GC mice without treatment, and reversed activities in GC mice treated with ivermectin. It should be noticed that the decrease in tumor size 2 months after ivermectin treatment was modest. As a matter of fact, in a separate experiment, we found that chemotherapy with 5-FU and oxaliplatin at the maximal dosage given to GC mice at the same age as ones in this study was without inhibition on the tumor size during 2 months of treatment (as same as in this study) (data not shown). However, the impacts of ivermectin treatment after a longer period of treatment alone and/or in combination with chemotherapy on resistance, migration and invasion could be worthwhile for future investigation. The results of the present study showed evidence of possible involvement of WNT/β-catenin signaling pathway in connection with the anti-cancer effect of ivermectin. For instance, prediction of ivermectin was successfully made by the WNT/β-catenin signaling pathway mining but not cMap. Validation of ivermectin was significant in silico model of the WNT/β-catenin signaling pathway. Up-regulation of WNT/β-catenin signaling pathway took place in patients, human cell lines and mouse model of GC. Ivermectin treatment induced downregulation of the WNT/β-catenin signaling pathway in the mouse GC. However, additional evidence is needed to demonstrate that the effect of ivermectin is dependent on WNT/β-catenin signaling pathway. For instance, it would be worthwhile to further investigate how modulation of the WNT/β‐ catenin signaling pathway with specific inhibitors and activators will affect the response to ivermectin treatment in vitro and in vivo.
The results of the present study show that ivermectin is a promising drug candidate for treatment of GC. The results may indicate an alternative mechanism of action of ivermectin, i.e., inhibition of the WNT/β-catenin signaling pathway in mammals rather than it acts on glutamate-gated chloride channels, which are common in nematodes, insects and ticks, thereby paralysing pharyngeal and somatic muscles. As ivermectin is exceptionally safe for mammals because of the blood/brain barrier, further pre-clinical and clinical studies of repositioning ivermectin for GC are warranted.
Data Availability Statement
The original contributions presented in the study are publicly available. This data can be found here: The mouse RNA seq/microarray data have been deposited in the NCBI Bioproject database under the accession number PRJNA690520 which can be accessed using the following link: http://www.ncbi.nlm.nih.gov/bioproject/690520. The human microarray data is available online via Mendeley Data repository with DOI link at http://dx.doi.org/10.17632/hzmfshy7hp.1. The RNAseq data in mouse GC after ivermectin treatment is available from the authors upon request.
The studies involving human participants were reviewed and approved by Regional Committees for Medical and Health Research Ethics Central Norway (REK 2012-1029). The patients/participants provided their written informed consent to participate in this study. The animal study was reviewed and approved by Mattilsynet. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
H-LR: In vitro and in vivo experiments, data analysis, manuscript writing; GTA: Patient samples, in vivo experiment and manuscript writing; AI and DK: cMap and manuscript writing; MKO: In vivo experiments and manuscript writing; JEG: Patient samples and manuscript writing; TCW: In vivo model, manuscript writing; DC: Study idea and design, manuscript writing; C-MZ: Study idea and design, in vivo experiment and patient samples; manuscript writing.
We thank the research grants and PhD fellowships supported by the Liaison Committee between the Central Norway Regional Health Authority (Helse-Midt Norge RHF) and Norwegian University of Science and Technology (NTNU) (grant numbers 46056636/46056928/90061700/90061701), Joint Program of the Medical Faculty of NTNU and St. Olavs University Hospital, the Cancer Foundation of St. Olavs Hospital (Kreftfondet ved St. Olavs hospital), and the technical support by Genomics Core Facility (GCF) which is funded by the Faculty of Medicine and Health Sciences at NTNU and RHF.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Andersen, P. I., Ianevski, A., Lysvand, H., Vitkauskiene, A., Oksenych, V., Bjørås, M., et al. (2020). Discovery and development of safe-in-man broad-spectrum antiviral agents. Int. J. Infect. Dis. 93, 268–276. doi:10.1016/j.ijid.2020.02.018
Asciutti, S., Akiri, G., Grumolato, L., Vijayakumar, S., and Aaronson, S. A. (2011). Diverse mechanisms of Wnt activation and effects of pathway inhibition on proliferation of human gastric carcinoma cells. Oncogene 30, 956–966. doi:10.1038/onc.2010.475
Benham-Pyle, B. W., Sim, J. Y., Hart, K. C., Pruitt, B. L., and Nelson, W. J. (2016). Increasing β-catenin/Wnt3A activity levels drive mechanical strain-induced cell cycle progression through mitosis. Elife 5, e19799. doi:10.7554/eLife.19799
Chen, I. S., and Kubo, Y. (2018). Ivermectin and its target molecules: shared and unique modulation mechanisms of ion channels and receptors by ivermectin. J. Physiol. 596, 1833–1845. doi:10.1113/JP275236
Chiou, R., Stubbs, R. J., and Bayne, W. F. (1987). Determination of ivermectin in human plasma and milk by high-performance liquid chromatography with fluorescence detection. J. Chromatogr. 416, 196–202. doi:10.1016/0378-4347(87)80502-9
GBD 2017 Stomach Cancer Collaborators (2020). The global, regional, and national burden of stomach cancer in 195 countries, 1990-2017: a systematic analysis for the global burden of disease study 2017. Lancet Gastroenterol. Hepatol. 5, 42–54. doi:10.1016/S2468-1253(19)30328-0
Diao, H., Cheng, N., Zhao, Y., Xu, H., Dong, H., Thamm, D. H., et al. (2019). Ivermectin inhibits canine mammary tumor growth by regulating cell cycle progression and WNT signaling. BMC Vet. Res. 15, 276. doi:10.1186/s12917-019-2026-2
Hayakawa, Y., Sakitani, K., Konishi, M., Asfaha, S., Niikura, R., Tomita, H., et al. (2017). Nerve growth factor promotes gastric tumorigenesis through aberrant cholinergic signaling. Cancer Cell 31, 21–34. doi:10.1016/j.ccell.2016.11.005
Hendriksen, C. F. M., and Morton, D. B.(Editors) (1999). “Humane endpoints in animal experiments for biomedical research,” in Proceedings of the International Conference, Zeist, The Netherlands, November 22–25 (1998), (London: Royal Society of Medicine Press Limited, Laboratory Animals Supplement).
Hollingshead, M. G., Stockwin, L. H., Alcoser, S. Y., Newton, D. L., Orsburn, B. C., Bonomi, C. A., et al. (2014). Gene expression profiling of 49 human tumor xenografts from in vitro culture through multiple in vivo passages--strategies for data mining in support of therapeutic studies. BMC Genomics 15, 393. doi:10.1186/1471-2164-15-393
Ikeda, T. (2003). Pharmacological effects of ivermectin, an antiparasitic agent for intestinal strongyloidiasis: its mode of action and clinical efficacy. Nippon Yakurigaku Zasshi 122, 527–(538.) [in Japanese, with English summary]. doi:10.1254/fpj.122.527
Kim, W., Cho, Y. S., Wang, X., Park, O., Ma, X., Kim, H., et al. (2019). Hippo signaling is intrinsically regulated during cell cycle progression by APC/CCdh1. Proc. Natl. Acad. Sci. U.S.A. 116, 9423–9432. doi:10.1073/pnas.1821370116
Liu, Z., Fang, H., Reagan, K., Xu, X., Mendrick, D. L., Slikker, W., et al. (2013). In silico drug repositioning: what we need to know. Drug Discov. Today 18, 110–115. doi:10.1016/j.drudis.2012.08.005
Mao, X., He, Z., Zhou, F., Huang, Y., and Zhu, G. (2019). Prognostic significance and molecular mechanisms of adenosine triphosphate-binding cassette subfamily C members in gastric cancer. Medicine 98, e18347. doi:10.1097/md.0000000000018347
Melotti, A., Mas, C., Kuciak, M., Lorente-Trigos, A., Borges, I., and Ruiz i Altaba, A. (2014). The river blindness drug Ivermectin and related macrocyclic lactones inhibit WNT-TCF pathway responses in human cancer. EMBO Mol. Med. 6, 1263–1278. doi:10.15252/emmm.201404084
Musa, A., Ghoraie, L. S., Zhang, S. D., Glazko, G., Yli-Harja, O., Dehmer, M., et al. (2018). A review of connectivity map and computational approaches in pharmacogenomics. Brief. Bioinform. 19, 506–523. doi:10.1093/bib/bbw112
Nappi, L., Aguda, A. H., Nakouzi, N. A., Lelj-Garolla, B., Beraldi, E., Lallous, N., et al. (2020). Ivermectin inhibits HSP27 and potentiates efficacy of oncogene targeting in tumor models. J. Clin. Invest. 130, 699–714. doi:10.1172/JCI130819
Pushpakom, S., Iorio, F., Eyers, P. A., Escott, K. J., Hopper, S., Wells, A., et al. (2019). Drug repurposing: progress, challenges and recommendations. Nat. Rev. Drug Discov. 18, 41–58. doi:10.1038/nrd.2018.168
Rabben, H. L., Andersen, G. T., Olsen, M. K., Øverby, A., Ianevski, A., Kainov, A., et al. (2021). Neural signaling modulates metabolism of gastric cancer. iScience 24, 102091. doi:10.1016/j.isci.2021.102091
Saito-Hakoda, A., Uruno, A., Yokoyama, A., Shimizu, K., Parvin, R., Kudo, M., et al. (2015). Effects of RXR agonists on cell proliferation/apoptosis and ACTH secretion/pomc expression. PLoS One 10, e0141960. doi:10.1371/journal.pone.0141960
Sitarz, R., Skierucha, M., Mielko, J., Offerhaus, J., Maciejewski, R., and Polkowski, W. (2018). Gastric cancer: epidemiology, prevention, classification, and treatment. Cancer Manag. Res. 10, 239–248. doi:10.2147/cmar.s149619
Subramanian, A., Narayan, R., Corsello, S. M., Peck, D. D., Natoli, T. E., Lu, X., et al. (2017). A next generation connectivity map: L1000 platform and the first 1,000,000 profiles. Cell 171, 1437–1452.e17. doi:10.1016/j.cell.2017.10.049
Verbaanderd, C., Maes, H., Schaaf, M. B., Sukhatme, V. P., Pantziarka, P., Sukhatme, V., et al. (2017). Repurposing drugs in oncology (ReDO)-chloroquine and hydroxychloroquine as anti-cancer agents. Ecancermedicalscience 11, 781. doi:10.3332/ecancer.2017.781
Yang, Z., Li, C., Yan, C., Li, J., Yan, M., Liu, B., et al. (2019). KIF14 promotes tumor progression and metastasis and is an independent predictor of poor prognosis in human gastric cancer. Biochim. Biophys. Acta Mol. Basis Dis. 1865, 181–192. doi:10.1016/j.bbadis.2018.10.039
Zhang, P., Zhang, Y., Liu, K., Liu, B., Xu, W., Gao, J., et al. (2019). Ivermectin induces cell cycle arrest and apoptosis of HeLa cells via mitochondrial pathway. Cell Prolif. 52, e12543. doi:10.1111/cpr.12543
Keywords: gastric cancer, ivermectin, drug repositioning, ingenuity pathway analysis, in silico
Citation: Rabben H-L, Andersen GT, Ianevski A, Olsen MK, Kainov D, Grønbech JE, Wang TC, Chen D and Zhao C-M (2021) Computational Drug Repositioning and Experimental Validation of Ivermectin in Treatment of Gastric Cancer. Front. Pharmacol. 12:625991. doi: 10.3389/fphar.2021.625991
Received: 04 November 2020; Accepted: 10 February 2021;
Published: 31 March 2021.
Edited by:Klara Gyires, Semmelweis University, Hungary
Reviewed by:Cedric Coulouarn, INSERM U1242 Laboratoire COSS, France
Pedro Miguel Rodrigues, Biodonostia Health Research Institute (IIS Biodonostia), Spain
Copyright © 2021 Rabben, Andersen, Ianevski, Olsen, Kainov, Grønbech, Wang, Chen and Zhao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chun-Mei Zhao, email@example.com