Download PDF
Original Article  |  Open Access  |  7 May 2024

Isocucurbitacin B inhibits glioma growth through PI3K/AKT pathways and increases glioma sensitivity to TMZ by inhibiting hsa-mir-1286a

Views: 277 |  Downloads: 50 |  Cited:   0
Cancer Drug Resist 2024;7:16.
10.20517/cdr.2024.01 |  © The Author(s) 2024.
Author Information
Article Notes
Cite This Article

Abstract

Aim: Glioma accounts for 81% of all cancers of the nervous system cancers and presents one of the most drug-resistant malignancies, resulting in a relatively high mortality rate. Despite extensive efforts, the complete treatment options for glioma remain elusive. The effect of isocucurbitacin B (isocuB), a natural compound extracted from melon pedicels, on glioma has not been investigated. This study aims to investigate the inhibitory effect of isocuB on glioma and elucidate its underlying mechanisms, with the objective of developing it as a potential therapeutic agent for glioma.

Methods: We used network pharmacology and bioinformatics analysis to predict potential targets and associated pathways of isocuB in glioma. Subsequently, the inhibitory effect of isocuB on glioma and its related mechanisms were assessed through Counting Kit-8 (CCK-8), wound healing, transwell, Western blot (WB), reverse transcription-quantitative polymerase chain reaction (RT-qPCR), and other in vitro experiments, alongside tumor formation experiments in nude mice.

Results: Based on this investigation, it suggested that isocuB might inhibit the growth of gliomas through the PI3K-AKT and MAPK pathways. Additionally, we proposed that isocuB may enhance glioma drug sensitivity to temozolomide (TMZ) via modulation of hsa-mir-1286a. The CCK-8 assay revealed that isocuB exhibited inhibitory effects on U251 and U87 proliferation and outperformed TMZ. Wound healing and transwell experiments showed that isocuB inhibited the invasion and migration of U251 cells by suppressing the activity of MMP-2/9, N-cadherin, and Vimentin. The TdT-mediated dUTP-biotin nick end labeling (TUNEL) and flow cytometry (FCM) assays revealed that isocuB induced cell apoptosis through inhibition of BCL-2. Subsequently, we conducted RT-qPCR and WB experiments, which revealed that PI3K/AKT and MAPK pathways might be involved in the mechanism of the inhibition isocuB on glioma. Additionally, isocuB promoted the sensitivity of glioma U251 to TMZ by inhibiting hsa-mir-1286a. Furthermore, we constructed TMZ-resistant U251 strains and demonstrated effective inhibition by isocuB against these resistant strains. Finally, we confirmed that isocuB can inhibit tumor growth in vivo through experiments on tumors in nude mice.

Conclusion: IsocuB may protect against glioma by acting on the PI3K/AKT and MAPK pathways and promote the sensitivity of glioma U251 to TMZ by inhibiting hsa-mir-1286a.

Keywords

Isocucurbitacin B, glioma, network pharmacology, hsa-mir-1286a, drug sensitivity

INTRODUCTION

The most challenging malignancies worldwide encompass stomach, lung, liver, colorectal, esophageal, breast cancer, and glioma. Gliomas account for 81% of central nervous system (CNS) malignancies and pose significant therapeutic challenges. They are classified into four grades (1-4) ranging from low to high grade[1,2]. According to the fifth edition of the WHO classification of tumors of the central nervous system, adult-type diffuse gliomas, pediatric-type diffuse low-grade gliomas, pediatric-type diffuse high-grade gliomas, and circumscribed astrocytic gliomas were featured[3]. There are numerous therapeutic methods for gliomas, which are mainly categorized into surgery, immunotherapy, targeted therapy, electric field therapy, and chemophototherapy[4,5]. At present, the predominant approach globally involves maximal surgical resection supplemented with temozolomide (TMZ) and radiotherapy. However, this regimen generally yields a dismal[6].

Modern pharmaceutical systems commonly use natural remedies to extract active ingredients with pharmacological properties. 60,000 years ago, plants were used as medicines, reflecting a wide range of pharmacological properties, including anticancer activity. Numerous natural compounds[7], such as quercetin[8], isoliquiritigenin[9], and lycopene[10], have apparent efficacy in inhibiting tumorigenesis, inhibiting tumor metastasis, and controlling tumor growth. Isocucurbitacin B (IsocuB) is isolated from the stalked tips of melons in the angiosperm family Cucurbitaceae. Natural cucurbitines exhibit various biological and pharmacological activities. Almost all the Cucurbitaceae families consist of approximately 229 species. Cucurbitacin is a class of highly oxidized tetracyclic triterpenoids with potent anticancer activity. Among its eight components of cucurbitacins with notable anticancer activity are cucurbitacin B, D, E, I, C, II, A, L-glucoside, S and R[11]. Its enantiomeric isomer, Cucurbitacin B[12], has long been known for its role in lung cancer[13,14], colorectal cancer[15], colon cancer[16], and glioma[17]. Melon pedicle was first documented in Shennong’s Herbal Classic, which describes its function in eliminating dampness and inducing vomiting. Furthermore, many ancient texts, such as “Shengji General Record”, “Qianjin Yaofang”, and “Ancient and Modern Medical Systems” also documented the use of melon pedicels for treating toothache, malaria, and hemorrhoids. Modern research has revealed that melon stalks have antitumor, hepatoprotective, and other beneficial effects. Clinical reports also indicate the utilization of melon stalks in the treatment of acute jaundice, infectious hepatitis, chronic rhinitis, chronic hepatitis, primary liver cancer, and other conditions. However, the anti-glioma effects of isocuB have not been proven.

While much is known about the molecular, structural, energetic, and chemical aspects of drug-target interactions, challenges remain in the selection and definition of targets, hindering the advancement of pharmacology and pharmacotherapy to an exact science[18]. The development of network pharmacology and bioinformatics will not only reduce the cost of drug development, but also shorten the time required for drug development. Furthermore, molecular docking techniques enable the prediction of the binding affinities and conformations of receptors and ligands[19,20]. This has been extremely helpful in our research. A deeper comprehension of the role of microRNAs (miRNAs) in development and disease, particularly in cancer, renders them attractive tools and targets for innovative therapeutic approaches. Functional studies have indicated the causal role of miRNA dysregulation in many cancers. MiRNAs show promise in preclinical development as tumor suppressors or oncogenes miRNA mimics, and molecules targeting miRNAs[21].

The aim of our study was to predict the target genes and pathways of isocuB against glioma using network pharmacology and data analysis. Afterwards, we validated these predictions through Counting Kit-8 (CCK-8), wound healing assay, transwell invasion, TdT-mediated dUTP-biotin nick end labeling (TUNEL) staining, reverse transcription-quantitative polymerase chain reaction (RT-qPCR), WB, and tumor growth experiments in nude mice. Then, we further investigated and predicted the relationship between miRNA and resistance to the drug TMZ. Our study demonstrated that isocuB exerts its inhibitory effects on glioma by regulating specific pathways and gene expression, suggesting a novel and effective treatment approach for glioma.

METHODS

Network pharmacology and databases analysis

Predicting the target genes of isocuB

In this study, we used the Comparative Toxicogenomics database (http://ctdbase.org/), the PharmMapper database (http://www.lilab-ecust.cn/pharmmapper/), and the SwissTargetPrediction database (https://www.sib.swiss/). We set our filter to “number of interactions > 1” to identify isocuB targets. In total, 300 potential targets predicted for isocuB were identified.

Search related therapeutic target genes of glioma

We simultaneously used “glioma” as a keyword in the Online Catalog of Human Genes and Disorders (OMIM) database (https://www.omim.org/), the CTD database, and the GeneCards database (https://www.genecards.org). We set the filter to “Relevance score ≥ 10” to search for therapeutic targets associated with glioma. After collating the results, we identified 22,320 relevant therapeutic targets by removing duplicate values.

Obtaining intersection target genes of isocuB and glioma

We utilized the Bioinformatics & Evolutionary Genomics website (http://bioinformatics.psb.ugent.be/webtools/Venn/) to identify common targets at the intersection of the aforementioned glioma therapeutic targets and isocuB target genes. Subsequently, we generated a Venn diagram.

Protein-protein interaction network construction

We imported the previously obtained intersection target genes of isocuB and glioma into the Search Tool for the Retrieval of Interacting Genes (STRING) database (https://string-db.org/). We specified the species as “Homo sapiens” and selected the parameter with the highest confidence level of 0.900 to retrieve a related file in the “.tsv” format. We then conducted a topological analysis using the Cytoscape software (https://cytoscape.org/) to identify the core genes with the strongest correlation.

Analysis of biological functions and pathways

To predict the correlation pathway and function of isocuB in glioma, we used the Database for Annotation, Visualization, and Integrated Discovery (DAVID) (https://david.ncifcrf.gov/) to conduct Gene Ontology (GO) process and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses (P < 0.05). Based on the above results, we utilized the R programming language to visualize the GO and KEGG pathway analyses. GO enrichment involves analyses of molecular function (MF), cell component (CC), and biological process (BP). KEGG is a bioinformatics resource used to identify a comprehensive list of genes that significantly impact metabolic pathways.

Molecular docking analysis

To improve the prediction of the relationship between isocuB and the top five related genes in the protein-protein interaction (PPI) analysis, we performed molecular association prediction. First, we obtained the structure of the gene from the Protein Data Bank (PDB) (http://www.rcsb.org) and acquired the 3D structure of isocuB from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/). Next, we removed the hydrogen and ligands and dehydrated the hub gene proteins using the PyMOL 1.7.x software (http://www.pymol.org/). The AutoDock Tools 1.5.6 software (https://autodock.scripps.edu/) was used to convert the file to the PDBQT format. The PyMOL software was used to visualize the docking results.

Research based on database data

The hub genes (RXRα, AKT1, ESR1, MAPK1, and HSP90AA1) were identified, and their mRNA expression levels in various tumor tissues were analyzed using the University of Alabama at Birmingham (UALCAN) platform (https://ualcan.path.uab.edu/analysis.html). We analyzed the mRNA expression levels in glioma and related factors using the Chinese Glioma Genome Atlas (CGGA) database (http://www.cgga.org.cn/). The top five gene mutations and copy number variants were analyzed using the cBioPortal for Cancer Genomics database (https://www.cbioportal.org/). The CGGA database was utilized to analyze the correlation between the mRNA expression levels of each gene and the clinical parameters and prognosis associated with glioma. The GEPIA database, which includes all glioma data and group cutoffs from The Cancer Genome Atlas (TCGA), was selected as the reference. The mRNA expression levels of each gene in the gliomas were confirmed using GEPIA. The mRNA expression levels of each gene in gliomas were verified, and their relationship with prognosis was analyzed.

Experimental verification

Cell culture and drug

Human glioma cells (U251 and U87, obtained from Shanghai Zhong Qiao Xin Zhou Biotechnology Co., Ltd., China) were cultured in complete DMEM (Shanghai Zhong Qiao Xin Zhou Biotechnology Co., Ltd., China). The environment for cultured cells needs to meet three conditions simultaneously: water saturation, a temperature of 37 °C, and 5% CO2 in the air. IsocuB, with a purity of ≥ 99.0%, was obtained from Chengdu Must Bio-technology Co., Ltd., China. TMZ was purchased from MCE in China with a purity of ≥ 99.0%. The drugs were dissolved in dimethyl sulfoxide (DMSO). The DMSO was diluted to less than 0.1% to ensure that it had no effect on the cells.

Reagents

CCK-8 (Dojindo, Japan), Matrigel Matrix (BD Biosciences, USA), paraformaldehyde (Solarbio, China), crystal violet (Sigma, USA), BeyoECL Plus (PIERCE, USA), TUNEL Assay Kit (Yesen, China), Annexin v-FITC/PI Apoptosis Kit (Boster, China), TRIzol Solution (Invitrogen, USA), Reverse Transcription System (Promega, USA), DEPC (Sigma, Germany), GoTaq® qPCR Master Mix (Promega, USA), Primers (Sangon Biotech, China), Bovine Serum Albumin (BSA) and BAC Protein Quantifier Kit (Beyotime Biotechnology, China), N-cadherin, Vimentin, BCl-2, Hsp90, p38-MAPK, p-MAPK, MMP-2, MMP-9, PDK1, anti-AKT, anti-p-AKT (CST, USA), p-STAT3 and STAT3 (Abclonal, China), β-actin (Servicebio, China), Bulge-Loop miRNA ORT-PCR Starter Kit, Human Bulge-Loop hsa-miR-1268a Primer Set, Bulge-Loop U6 gPCR Primer Set, Human microFF hsa-miR-1286 Inhibitor and microFF Inhibitor NC (RIBOBIO, China), micrON hsa-miR-1286 mimic (Genepharma, China) were used in this study.

Cell proliferation assay

We evaluated the effect of isocuB on glioma cell proliferation using a CCK-8 assay. We initially used 0.25% trypsin to digest U251 and U87 cells at the logarithmic growth stage and subsequently diluted them to a seeding density of 7 × 104 cells/mL in 96-well plates. IsocuB was then diluted to seven concentrations (0.001, 0.01, 0.1, 0.5, 1, 5, and 25 µmol/L). TMZ was then diluted to seven concentrations (10, 50, 100, 200, 500, 1,000, and 2,000 µmol/L). These were added 12 and 24 h later. Finally, the CCK-8 reagent was added at 12 and 24 h for absorbance detection at 450 nm using a microplate reader. Triplicate experiments were performed independently {Cell viability = [(ODExperiment group - ODBlank)/(ODControl group - ODBlank)] × 100%}.

Wound healing assay

We assessed the impact of isocuB on glioma cell migration by conducting a wound-healing assay. We diluted the cell density of U251 to 2.5 × 105 cells/mL and plated it into a 6-well plate. Five parallel lines (2 mL per well) were drawn beside each well with a marker. Then, 24 h later, the lines were drawn with a 200 µL pipette tip tilted at a 45° angle, washed three times with PBS, and DMEM dilution containing 1% FBS in three concentrations of isocuB (0, 0.01, and 0.1 µmol/L) was added. Finally, five visual fields were randomly selected at 12, 24, and 36 h to capture images and observe the migration of U251 cells. Triplicate experiments were performed independently.

Transwell invasion

To assess the invasive impact of isocuB on glioma cells, a transwell invasion assay was conducted. The diluted matrix adhesive was added to the upper chamber and placed in an incubator for 8 h to allow the matrix to form and absorb excess fluid. After 12 h of cell starvation, 5 × 104 cells were seeded into the upper chamber with 100 μL of 0.2% BSA DMEM containing isocuB (0, 0.01, and 0.1 µmol/L), and 500 μL of complete DMEM medium was added to the lower chamber. After 24 h, non-migrated cells on the filter side of the upper chamber were removed using a cotton swab, fixed with 4% paraformaldehyde for 1 h, rinsed three times with PBS, and stained with 1 mL of crystal violet for 1 h. The transwell membrane was then covered with a cover glass, and the migrated cells were counted under a microscope. Triplicate experiments were performed independently.

TUNEL staining

The U251 cells in the growth stage were seeded in 6-well plates at a density of 2.5 × 105 cells/mL, with 2 mL per well. The drug (0, 0.1, 0.5, and 1 µmol/L) was added after cell adhesion. After a 24-hour incubation, the cells were harvested using ethylenediaminetetraacetic acid-free pancreatic enzymes and washed twice with PBS. Each sample was supplemented with 50 µL of FITC-12-dUTP and incubated at 37 °C for 1 h, followed by testing using a machine. Triplicate experiments were performed independently.

Flow cytometry analysis of cell apoptosis

The U251 cells in the growth stage were seeded in 6-well plates at a density of 2.5 × 105 cells/mL, with 2 mL per well. The drug (0, 0.1, 0.5, and 1 µmol/L) was added after cell adhesion. After a 24-hour incubation, the cells were harvested using ethylenediaminetetraacetic acid-free pancreatic enzymes and washed twice with PBS. 500 µL of 1× Annexin V Binding Buffer was added to each sample. Subsequently, 5 µL of annexin V-FITC staining solution and 5 µL of PI staining solution were added to each sample, incubated at room temperature for 5-15 min, and then analyzed by machine. Triplicate experiments were performed independently.

RT-qPCR analysis

To investigate the mechanism of action of isocuB in glioma, we used RT-qPCR. First, RNA was extracted from U251 cells treated with isocuB (0, 0.1, 0.5, 1 µmol/L) using TRIzol. RNA concentration and purity were determined using a miRNA protein quantification assay. The assay was performed according to the instructions provided with the SYBR® Green Real-time PCR Master Mix kit. Cycle threshold (Ct) was determined using quantitative PCR. The primer sequences for MMP-2, MMP-9, PDK1, RXRα, PPARα, Bcl-2, and β-actin are provided in Table 1. Additionally, we analyzed the RNA content of hsa-mir-1286a using U6 as the internal reference. Hsa-miR-1286a inhibitor and microRNA inhibitor NC were then added to the 6-well plate at a concentration of 100 µmol/L. RNA was extracted 24 h later, and the level of hsa-miR-1286a was detected. Moreover, hsa-miR-1286a mimic was added to the 6-well plate at a concentration of 100 µmol/L. RNA was extracted 24 h later, and the level of hsa-miR-1286a was detected. Triplicate experiments were performed independently.

Table 1

Primers used for RT-qPCR analysis

GenePrimer sequence
Forward (5’-3’)Reverse (5’-3’)
MMP-2GATACCCCTTTGACGGTAAGACCTTCTCCCAAGGTCCATAGC
MMP-9GGGACGCAGACATCGTCATCTCGTCATCGTCGAAATGGGC
PDK1CGTACGGCAATGGCTTTATCAATCCCCTCCTGCAACTTCT
RXRαATGGACACCAAACATTTCCTGGGGAGCTGATGACCGAGAAAG
PPARαTTCGCAATCCATCGGCGAGCCACAGGATAAGTCACCGAGG
Bcl-2GGTGGGGTCATGTGTGTGGCGGTTCAGGTACTCAGTCATCC
β-actinAAAGCGGCTGTTAGTCACTGGCGAGTCATTGCATACTGTCCAT

WB analysis

To further investigate the mechanism of action of isocuB on glioma, the proteins in U251 cells treated with different concentrations of isocuB (0, 0.1, 0.5, 1 µmol/L) for 24 h were extracted using RIPA lysis. The extracted protein solution was then assayed using the BCA kit and denatured using SDS. Proteins were separated at a constant voltage of 110 V for 90 min and then transferred to a PVDF membrane (Mannheim, Germany). After blocking with 5% skim milk powder, the primary antibodies for MMP-2, MMP-9, p-PI3K, PI3K, p-AKT(S473), AKT, p-STAT3, STAT3, RXRα, PDK1, p-MAPK1/3, MAPK1/3, N-cadherin, Vimentin and Bcl-2 (1:1000) were added and incubated overnight at 4 °C. After washing with TBST, appropriate secondary antibodies were added, incubated at 37 °C for 1 h, and then washed with TBST and TBS. Luminescent agents were then added to facilitate detection in the chemiluminescence imaging system, and further analysis was performed using ImageJ software. Triplicate experiments were performed independently.

Establishment of U251/TMZ resistant strains

Cells in the logarithmic growth phase (80%-90%) were treated with drugs at a low initial concentration (1/10-1/5 of IC50 of the parent cell line is recommended), and cultured in an incubator at 37 °C and 5% CO2. When the cell density reached 50%, the culture medium was abandoned, PBS was cleaned twice, and the drug-free medium was replaced for further culture. When the cell growth density recovered to 80%-90%, the above drug treatment was repeated 6-8 times. The final drug concentration was stable and the drug-resistant cell line was obtained. The IC50 values of drug-resistant cell lines were determined, and the resistance index (RI) was calculated using the formula: RI = IC50 of drug-resistant cell lines/IC50 of parental cell lines. The RI grades range from 1 to 5, indicating low drug resistance, 5 to 15, indicating moderate resistance, and above 15, indicating high resistance. Triplicate experiments were performed independently.

Inhibitory effect of isocuB on U251/TMZ resistant strains

The primary physiological characteristics of tumor cells include proliferation, invasion, and migration. We evaluated the effect of isocuB on glioma cell proliferation using the CCK-8 assay. We initially used 0.25% trypsin to digest U251/TMZ-resistant strain cells at the logarithmic growth stage and subsequently diluted them to a seeding density of 7 × 104 cells/mL in 96-well plates. TMZ was then diluted to seven concentrations (10, 50, 100, 200, 500, 1,000, and 2,000 µmol/L). IsocuB was then diluted to seven concentrations (0.001, 0.01, 0.1, 0.5, 1, 5, and 25 µmol/L). These were added 24 h later. Finally, the CCK-8 reagent was added at 24 h for absorbance detection at 450 nm using a microplate reader. Triplicate experiments were performed independently {Cell viability = [(ODExperiment group - ODBlank)/(ODControl group - ODBlank)] × 100%}.

Transfection of hsa-mir-1286a inhibitor increased the sensitivity of glioma U251 to TMZ

We first used 0.25% trypsin to digest U251 cells at the logarithmic growth stage, and then diluted them to a seeding density of 7 × 104 cells/mL in 96-well plates. Human microRNA hsa-mir-1286a inhibitor was then added to the 6-well plate at a concentration of 100 µmol/L, and RNA was extracted 24 h later. TMZ was then diluted to seven concentrations (10, 50, 100, 200, 500, 1,000, and 2,000 µmol/L) and added 24 h later. Finally, the CCK-8 reagent was added after 24 h for absorbance measurement at 450 nm using a microplate reader. Triplicate experiments were performed independently {Cell viability = [(ODExperiment group - ODBlank)/(ODControl group - ODBlank)] × 100%}.

Transfection of hsa-miR-1286a mimic inhibited the sensitivity of glioma U251 to TMZ

We initially used 0.25% trypsin to digest U251 cells at the logarithmic growth stage, and subsequently diluted them to a seeding density of 7 × 104 cells/mL in 96-well plates. A hsa-miR-1286a mimic was then added to the 6-well plate at a concentration of 100 µmol/L, and RNA was extracted 24 h later. TMZ was then diluted to seven concentrations (10, 50, 100, 200, 500, 1,000, and 2,000 µmol/L) and added 24 h later. Finally, the CCK-8 reagent was added after 24 h for absorbance measurement at 450 nm using a microplate reader. Triplicate experiments were performed independently {Cell viability = [(ODExperiment group - ODBlank)/(ODControl group - ODBlank)] × 100%}.

Tumor formation experiment in nude mice

In our study, 16 six-week-old specific pathogen-free (SPF) BALB/c nude mice were initially transplanted with 5 × 106 U251 cells in a 100 µL PBS volume in the right upper axilla. The mice were then divided into two groups (IsocuB group at 2 mg/kg and DMSO group) on day 6 after transplantation. During the experiment, tumor size was measured daily in both groups of nude mice. On day 18, the masses were removed and weighed. We assessed and compared the supervisors. These studies were conducted in compliance with the ARRIVE guidelines[22,23] and the National Research Council Guide for the Care and Use of Laboratory Animals. All animal experiments followed the regulations of ethics committee of the Experimental Animal Administration of Sichuan University (NO. K2024006).

Statistical analysis

To ensure the accuracy of the experimental results, each set of experiments was repeated at least three times. The final data were presented as the mean ± standard deviation (SD), and the differences were statistically analyzed using GraphPad software (P < 0.05 was considered statistically significant).

RESULTS

Network pharmacology and databases analysis

In this study, the targets of isocuB [Figure 1A] were predicted using the CTD, Swiss Target Prediction, and PharmMapper databases. A total of 301 predicted targets were identified as potential targets of isocuB. In this study, we identified potential therapeutic target genes for glioma using the GeneCards, OMIM, and CTD databases. The results were compiled, yielding 22,320 relevant therapeutic targets. We imported the 280 common genes [Figure 1B, Supplementary Table 1] into the STRING database and obtained the PPI information for 280 nodes and 631 edges [Figure 1C]. The corresponding PPI information was imported into Cytoscape software. The analysis revealed that the top five genes affected by isocuB in glioma were RXRα, AKT1, ESR1, MAPK1, and HSP90AA1 [Figure 1D, Supplementary Table 2]. A total of 826 GO analyses were performed, and the first 10 enrichment results are visualized in Figure 1E based on the P value. The analyses were mainly divided into three categories: (1) 278 BP, including response to protein binding, peptidyl-tyrosine phosphorylation, and protein autophosphorylation; (2) 279 CC, including cytosol, extracellular exosome, and extracellular region; (3) 279 MF, including signaling receptor activator activity and identical protein binding activity. KEGG analysis revealed 152 major pathways, with the 20 most significantly enriched pathways shown in Figure 1F. These pathways included the PI3K/AKT signaling pathway, MAPK signaling pathway, and proteoglycans in cancer signaling pathways. The top five targets were found to be highly connected in the network according to the PPI analysis, located at the core of the network, and identified as the most important nodes. The binding affinities of isocuB to RXRα (3HOA), AKT1 (4EJN), ESR1 (6V8T), MAPK1 (6G54), and HSP90AA1 (4BQG) were -7.41, -6.67, -7.65, -7.95, and -7.68 kcal/mol [Figure 1G]. In addition, a coupling fraction of less than 0 kcal/mol indicates that the component can spontaneously bind to the target, less than -4.25 kcal/mol indicates a good affinity coupling, and less than -7 kcal/mol is considered as strong affinity coupling[24]. The results indicated a strong binding activity between the main components and the hub gene (affinity < -6.00 kcal/mol).

Isocucurbitacin B inhibits glioma growth through PI3K/AKT pathways and increases glioma sensitivity to TMZ by inhibiting hsa-mir-1286a

Figure 1. (A) IsocuB, drawn using Chemdraw; (B) Obtain 280 intersection target genes for isocuB and glioma; (C) STRING database of the 280 potential therapeutic targets; (D) PPI network of the 280 potential therapeutic targets; (E) The top 10 significant enrichment terms in BP, CC, and MF are shown in a GO enrichment bubble map containing 280 intersecting targets; (F) The top 20 most significantly enriched pathways Bar diagram of KEGG pathway enrichment analysis with 280 cross-target sites; (G) The result diagram of molecular docking. The yellow 3D structural formula represents isocuB; the blue bond represents the hydrogen bond at the binding site. IsocuB: Isocucurbitacin B; STRING: search tool for the retrieval of interacting genes; PPI: protein-protein interaction; BP: biological process; CC: cell component; MF: molecular function; GO: Gene Ontology; KEGG: Kyoto encyclopedia of genes and genomes.

Based on RNA-Seq data from the UALCAN database, a total of 7932 case samples were included. The mRNA expression of RXRα, AKT1, ESR1, MAPK1, and HSP90AA1 was analyzed in different tumor tissues. The mRNA expression levels of AKT1 and MAPK1 were relatively high in glioma, while those of RXRα, ESR1, and HSP90AA1 were relatively low [Figure 2A]. The mRNA expression of RXRα, AKT1, ESR1, MAPK1, and HSP90AA1 in tumor and normal tissues was analyzed using the GEPIA database. The mRNA expression of RXRα, ESR1, and HSP90AA1 did not show significant differences between tumor tissues and normal tissues [Figure 2B]. AKT1 and MAPK1 mRNA expression levels were significantly different in tumor and normal tissues in gliomas (P < 0.05) [Figure 2B]. According to the KEGG analysis, the top two pathways were PI3K-AKT and MAPK. Molecular docking showed a strong binding effect. The mRNA expression levels of AKT1 and MAPK1 were significantly different between normal and glioma tissues. Next, we analyzed the mRNA expression levels of MAPK1 and AKT1. Gene mutations and copy number variations of MAPK1 and AKT1 in gliomas were analyzed using the cBioPortal for Cancer Genomics database. Based on an analysis of the cBioPortal for Cancer Genomics database, seven studies included 3436 cases of glioma.

Isocucurbitacin B inhibits glioma growth through PI3K/AKT pathways and increases glioma sensitivity to TMZ by inhibiting hsa-mir-1286a

Figure 2. (A) mRNA expression levels in common tumor tissues; (B) The mRNA expression levels in common tumor tissues (*P < 0.05 vs. control group); (C) Mutation, deep deletion, and amplification of AKT1 in gliomas. Green, blue, and red indicate mutations, deep deletions, and amplifications, respectively. AKT1 methylation is negatively correlated with the expression of AKT1 in glioma tissues (****P < 0.0001 vs. control group). The AKT1 gene has a low methylation level in gliomas. (****P < 0.0001 vs. control group); (D) Mutations, deep deletions, and MAPK1 amplification in glioma cells. Green indicates mutations. Different types of MAPK1 mutations. Different types of MAPK1 mutation in glioma at different levels. AKT1: AKT serine/threonine kinase 1; MAPK1: mitogen-activated protein kinase 1.

AKT1 contains mutations, deep deletions, and amplifications in gliomas. AKT1 exhibited the highest mutation rate of 3.81%, amplification rate of 2.02%, and deep deletion rate of 0.34% [Figure 2C]. A total of 189 samples from two studies were included in the methylation analysis. The results showed a negative correlation between AKT1 methylation and AKT1 gene expression in glioma tissues (P < 0.001), indicating that AKT1 gene expression may be epigenetically regulated. Copy number amplification of AKT1 was detected in three cases, while mild AKT1 copy number loss was detected in 29 cases. According to the pathological type, all 189 cases of GBM (100%) had low AKT1 methylation. Mutations in MAPK1 were identified in gliomas. MAPK1 exhibited the highest mutation rate (1.9%) without any apparent significant deletions or amplifications [Figure 2D]. Different mutations and different types of gliomas associated with MAPK1 are shown in Figure 2D.

To investigate the relationship between AKT1 and MAPK1 mRNA expression and clinicopathologic parameters, we used the CGGA database, which comprises 325 cases from the glioma mRNA microarray database, for correlation analysis. The cutoff values for AKT1 and MAPK1 mRNA expression were 70.48 and 44.84, respectively. Values higher than the cutoff indicated high expression, and vice versa. The results indicated that AKT1 and MAPK1 were present in different tumor stages. There were statistically significant differences in the expression of glioma grades (WHO II, WHO III, and WHO IV) (P < 0.001) [Figure 3A and B], as well as statistically significant differences among different states of isocitrate dehydrogenase (IDH) (wild type, mutant type) (P < 0.001). There was a statistically significant difference among the various age groups (P < 0.05) [Figure 3A and B]. Expression of AKT1 and MAPK1 was correlated with WHO classification, IDH status, and age. The TCGA glioma data were retrieved from the GEPIA database. The group cutoff was selected as the median to divide the high expression group and the low expression group (cutoff-high = 50%; cutoff-low = 50%) for analyzing the mRNA expression of AKT1 and MAPK1 in glioma to verify their prognostic relationship. To validate the analysis results from the CGGA database and assess the prognostic significance of AKT1 and MAPK1 in glioma, survival analyses were performed using the GEPIA database. The results, using the Kaplan-Meier method, showed that there were 676 cases of TCGA glioma. The prognosis of the low AKT1 and MAPK1 expression groups was better than that of the high expression group, and the differences in overall survival (OS) between the two groups were statistically significant (P < 0.01, P < 0.0001). Disease-free survival of the low AKT1 and MAPK1 expression groups was better than that of the high expression group. The differences were significant (P < 0.05) [Figure 3C and D].

Isocucurbitacin B inhibits glioma growth through PI3K/AKT pathways and increases glioma sensitivity to TMZ by inhibiting hsa-mir-1286a

Figure 3. (A) Relationship between mRNA expression levels of AKT1 and clinicopathological features in patients with glioma. Significant differences in the expression of (WHO II, WHO III, and WHO IV) (****P < 0.0001 vs. control group). Significant differences among different states of IDH (wild type, mutant type) (****P < 0.0001 vs. control group). Significant differences among different age status (</≥ 42) (*P < 0.05 vs. control group); (B) Relationship between mRNA expression levels of MAPK1 and clinicopathological features in patients with glioma. Significant differences in the expression of (WHO II, WHO III, and WHO IV) (****P < 0.0001 vs. control group). Significant differences among different states of IDH (wild type, mutant type) (****P < 0.0001 vs. control group); Significant differences among different age status (</≥ 42) (****P < 0.0001 vs. control group); (C and D) The prognosis of the low AKT1 and MAPK1 expression groups was better than that of the high expression group (*P < 0.05, ****P < 0.0001 vs. control group). The disease-free survival of the low AKT1 and MAPK1 expression groups was better thanthat of the high expression group (*P < 0.05 vs. control group). AKT1: AKT serine/threonine kinase 1; IDH: isocitrate dehydrogenase; MAPK1: mitogen-activated protein kinase 1.

Experimental verification results

IsocuB inhibits proliferation, migration and invasion and increases apoptosis

Our experimental results demonstrated that isocuB had a concentration- and time-dependent effect on glioma cell proliferation in the CCK8 assay [Figure 4A]. The IC50 of isocuB-inhibited U251 cells is 0.79 µmol/L at 24 h and 10.54 µmol/L at 12 h. The IC50 of isocuB-inhibited U87 cells is 2.12 µmol/L at 24 h [Figure 4B]. Thus, the inhibitory effect of isocuB on U251 was much greater than that on U87; therefore, U251 was primarily used in our subsequent experiments. Our results demonstrate that isocuB inhibited glioma migration in a time- and dose-dependent manner, as illustrated in the wound healing assay [Figure 4C and D].

Isocucurbitacin B inhibits glioma growth through PI3K/AKT pathways and increases glioma sensitivity to TMZ by inhibiting hsa-mir-1286a

Figure 4. (A) The viability of U251 cells after 12 and 24 h of isocuB treatment (**P < 0.01, ***P < 0.001; ****P < 0.0001 vs. control group); (B) The viability of U87 cells after 24 h of isocuB treatment (***P < 0.001; ****P < 0.0001 vs. control group); (C and D) Cell mobility of U251 cells after12, 24 and 36 h of isocuB treatment (*P < 0.05, **P < 0.01, ***P < 0.001 vs. control group). isocuB: Isocucurbitacin B.

In addition, isocuB inhibited the invasion of glioma cells in the transwell invasion assay [Figure 5A and B]. The MMP family is associated with cancer cell survival, proliferation, apoptosis, invasion, and metastasis[25]. The mRNA and protein expression levels were significantly decreased after interference with isocuB in U251 cells [Figure 5C-E]. Epithelial-mesenchymal transition (EMT) is a process in which cells acquire invasive mesenchymal competence through a series of events that include loss of cellular connections, cytoskeletal reorganization, and remodeling of the extracellular matrix[26]. The results of the WB experiment results showed a significant decrease in the protein levels of N-cadherin and Vimentin. This indicated that isocuB could inhibit EMT and consequently inhibit glioma invasion [Figure 5C and D].

Isocucurbitacin B inhibits glioma growth through PI3K/AKT pathways and increases glioma sensitivity to TMZ by inhibiting hsa-mir-1286a

Figure 5. (A and B) The transwell invasion process is measured by the number of cells after 24 h of treatment (**P < 0.01, ***P < 0.001 vs. control group); (C and D) Protein expressions of MMP-2/9, N-cadherin, and Vimentin were analyzed using WB. Data are represented as mean ± SD (*P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001 vs. control group vs. control group); (E) Expression of MMP-2 and MMP-9 RNA was detected using RT-qPCR (*P < 0.05, **P < 0.01, ***P < 0.001 vs. control group). MMP-2/9: Matrix metalloproteinase 2/9; WB: Western blot; RT-qPCR: reverse transcription-quantitative polymerase chain reaction.

The TUNEL staining assay results indicated that isocuB promoted the apoptosis of U251 cells [Figure 6A and B]. Annexin V was used to detect phosphatidylserine (PS) exposure on the outer membrane of early apoptotic cells. However, annexin V cannot distinguish between poor cell death (middle and late apoptotic cells) and early apoptotic cells. PI can enter the necrotic cells (middle and late apoptotic cells) but is excluded from early apoptotic cells. This finding indicated that isocuB could induce the apoptosis of glioma U251 cells. The proportion of late apoptosis also increased with increasing drug concentration [Figure 6C and D]. WB analysis showed that the proportion of BCL-2 increased in U251 cells after administration [Figure 6E-G].

Isocucurbitacin B inhibits glioma growth through PI3K/AKT pathways and increases glioma sensitivity to TMZ by inhibiting hsa-mir-1286a

Figure 6. (A and B) Dye content of FITC-12-dUTP was detected by FCM (****P < 0.0001 vs. control group); (C and D) The apoptosis rate was increased by FCM. Q3-1 represents normal cells, Q3-4 represents viable apoptotic cells, Q3-2 represents late apoptotic cells, and Q3-1 represents mechanical destruction of cells (**P < 0.01, ***P < 0.001, ****P < 0.0001 vs. control group); (E and F) Protein expression of BCL-2 was analyzed using WB. Data are represented as mean ± SD (*P < 0.05, **P < 0.01 vs. control group); (G) Expression of BCL-2 RNA was detected using RT-qPCR (n = 3) (*P < 0.05, **P < 0.01, ***P < 0.001 vs. control group). FCM: Flow cytometry; BCL-2: B-cell lymphoma-2; WB: western blot; RT-qPCR: reverse transcription-quantitative polymerase chain reaction.

IsocuB inhibited the activation of the PI3K/AKT, MAPK, and STAT3 signaling pathways

We demonstrated that isocuB inhibited the glioma pathway based on WB and PCR results. IsocuB was found to inhibit glioma growth by suppressing PDK1, Bcl-2, PI3K-AKT, and MAPK signaling pathways, as well as the activation of MMP-2/9. Our results showed that the level of p-MAPK1/3 protein decreased significantly, while the level of total MAPK1/3 protein remained unchanged in the groups treated with different concentrations. In addition, the levels of total PI3K and t-AKT protein decreased significantly, while the levels of p-PI3K and p-AKT (S473) did not change significantly in the groups treated with different concentrations. The protein levels of PDK1 were significantly decreased [Figure 7A]. Furthermore, p-STAT3 protein levels were significantly reduced in the groups treated with different concentrations, but total STAT3 protein levels remained unchanged [Figure 7B]. IsocuB reduced the mRNA expression levels of PDK1, RXRα, PPARα, and Bcl-2 in U251 cells as determined by an RT-qPCR assay [Figure 8A]. These results suggested that isocuB played an anticancer role through the PI3K-AKT and MAPK pathways.

Isocucurbitacin B inhibits glioma growth through PI3K/AKT pathways and increases glioma sensitivity to TMZ by inhibiting hsa-mir-1286a

Figure 7. (A) Protein expressions of p-PI3K, PI3K, p-AKT(S473), t-AKT, RXRα, PDK1, p-MAPK1/3, MAPK1/3 and Bcl-2 were analyzed using WB. Data are represented as mean ± SD (*P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001 vs. control group); (B) Protein expressions of p-STAT3, STAT3 were analyzed using WB. Data are represented as mean ± SD (**P < 0.01 vs. control group). PI3K: Phosphoinositide 3-kinase; AKT: AKT serine/threonine kinase; PDK1: 3-phosphoinositide-dependent protein kinase 1; MAPK: mitogen-activated protein kinase; WB: western blot; STAT3: signal transducer and activator of transcription 3.

Isocucurbitacin B inhibits glioma growth through PI3K/AKT pathways and increases glioma sensitivity to TMZ by inhibiting hsa-mir-1286a

Figure 8. (A) Expression of pdk1, RXRα, PPAα, and Bcl-2 RNA was detected using RT-qPCR (*P < 0.05, **P < 0.01 vs. control group); (B) The viability of U251 cells after 12 and 24 h of TMZ treatment (*P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001 vs. control group); (C) The viability of U251/TMZ resistant strains after 24 h of TMZ treatment (*P < 0.05, ***P < 0.001, ****P < 0.0001 vs. control group); (D) The viability of U251/TMZ resistant strains after 24 h of isocuB treatment (*P < 0.05, ***P < 0.001, ****P < 0.0001 vs. control group). RT-qPCR: Reverse transcription-quantitative polymerase chain reaction; TMZ: temozolomide; isocuB: isocucurbitacin B.

IsocuB inhibited U251/TMZ resistant strains

We investigated the inhibitory effects of 12-hour and 24-hour TMZ on U251 glioma cells. The IC50 values for 12-hour and 24-hour TMZ were 1,055 and 819.4 µmol/L, respectively. These values were significantly higher than the IC50 of isocuB at 12 and 24 h. Therefore, we could conclude that isocuB was more effective than TMZ [Figure 8B]. We established a resistant strain of U251/TMZ with an IC50 of 2,413 µmol/L after 24 h of TMZ administration [Figure 8C]. The RI grades range from 1 to 5, indicating low drug resistance, 5 to 15, indicating moderate resistance, and above 15, indicating high resistance. Thus, the RI is 2.95 and we have established a low drug resistance. Then, we verified the inhibitory effect of isocuB on U251/TMZ-resistant strains with an IC50 of 1.01µmol/L [Figure 8D].

IsocuB increased the pharmacological sensitivity to TMZ with the regulation of hsa-mir-1286a expression levels

Based on the literature, we hypothesized that hsa-mir-1286a would be associated with TMZ resistance after predicting the related miRNAs of AKT1 and MAPK1. Using the CGGM database, we examined the prognostic significance of hsa-mir-1286a in grade 4 gliomas and found that reduced expression was linked to better survival for these patients [Figure 9A]. The level of hsa-mir-1286a decreased steadily as the drug concentration increased according to the RT-qPCR data [Figure 9B]. Following the addition of micrOFF inhibitor NC and hsa-miR-1286a inhibitor, the RNA levels of hsa-mir-1286a decreased [Figure 9C]. Furthermore, the IC50 of U251 cells was significantly decreased by TMZ with the addition of the microRNA hsa-mir-1286a inhibitor; the IC50 was determined to be 441.2 µmol/L after 24 h [Figure 9D]. The RNA levels of hsa-mir-1286a decreased after the addition of micrON hsa-miR-1286a mimic [Supplementary Figure 1]. Moreover, after the addition of the micrON hsa-miR-1286a mimic, the IC50 of U251 cells was significantly reduced by TMZ, with the IC50 measured at 1,151 µmol/L after 24 h [Supplementary Figure 2].

Isocucurbitacin B inhibits glioma growth through PI3K/AKT pathways and increases glioma sensitivity to TMZ by inhibiting hsa-mir-1286a

Figure 9. (A) Prognosis of low and high expression levels of hsa-mir-1286a in grade 4 gliomas; (B) Expression of hsa-mir-1286a was detected using RT-qPCR (*P < 0.05, **P < 0.01 vs. control group); (C) Expression of hsa-mir-1286a after using hsa-miR-1286a inhibitor was detected using RT-qPCR (**P < 0.01 vs. control group); (D) Cell viability of U251 cells after 24 h of TMZ treatment after microFF hsa-miR-1286 inhibitor added (****P < 0.0001 vs. control group). RT-qPCR: Reverse transcription-quantitative polymerase chain reaction; TMZ: temozolomide.

IsocuB inhibited tumor growth in nude mice

Previous studies of the toxicity of isocuB in C57BL/6 mice showed that the drug had no effect on the mice at doses below 2 mg/kg. Consequently, isocuB was administered to this experimental species at a dose of 2 mg/kg. After 18 days of intraperitoneal injection of 2 mg/kg/day of isocuB, the mass size was measured daily and the weight on the last day was recorded. “V (mm3) = L × W2/2” is the formula used to calculate the mass. Our results demonstrated that isocuB strongly suppressed cell proliferation. After administration of isocuB for 18 days (P < 0.0001), the weight of the mass decreased dramatically compared to the control group [Figure 10A-F]. Taken together, isocuB inhibited tumors in vivo.

Isocucurbitacin B inhibits glioma growth through PI3K/AKT pathways and increases glioma sensitivity to TMZ by inhibiting hsa-mir-1286a

Figure 10. IsocuB inhibits tumor growth in nude mice. (A) Bossing of the mass in nude mice: the control group is shown above, and the treatment group is shown below; (B) Bossing: control group at the top and treatment group at the bottom; (C) Body weight changes in nude mice; (D) Daily tumor volume changes after administration (**P < 0.01, ****P < 0.0001 vs. control group); (E) The bossing weight (****P < 0.0001 vs. control group); (F) Bossing weight/body weight (****P < 0.0001 vs. control group). IsocuB: Isocucurbitacin B.

DISCUSSION

Accounting for 81% of all CNS malignancies, glioma is one of the most challenging malignancies to treat effectively. Recent studies have underscored the utility of exploring active compounds in natural medicinal materials as a promising strategy for developing novel medications[7]. IsocuB belongs to the class of terpenoids, which have been found to interact with a variety of signaling pathways, such as the MAPK, nuclear factor κB (NF-κB), and glycolytic enzyme pathways[27]. Namely, it is interesting to verify whether IsocuB may have anti-glioma properties.

Network pharmacology integrates systems biology, pharmacology, and computational analysis techniques to investigate the complex relationships between compounds, diseases, and targets[28]. Our study revealed 280 overlapping target genes of isocuB and glioma. Through PPI analysis, the top five genes were RXRα, AKT1, ESR1, MAPK1, and HSP90AA1. The GO and KEGG results revealed that the primary pathways targeted by isocuB against glioma were the PI3K-AKT and MAPK signaling pathways. AKT1 and MAPK1 were identified as the pivotal gene targets within these pathways, respectively. Previous studies have shown that the development, progression, and deterioration of various tumors are associated with the activation of the PI3K/Akt signaling pathway[2,29]. The MAPK signaling pathway, comprising several crucial signaling components and phosphorylation events, plays a pivotal role in tumorigenesis, including breast cancer[30] and glioma[31]. Therefore, isocuB may exert its anti-glioma effect through the PI3K/AKT and MAPK signaling pathways. Molecular docking is a technique used in drug design to simulate receptor-drug interaction patterns. In recent years, the utilization of molecular docking to elucidate relevant mechanisms of action has become a trend in new drug development. Our results showed that isocuB exhibited a strong binding affinity for the top five genes, implying its potential in the growth of glioma through the action of these genes, which should be further analyzed and investigated. In this study, we analyzed the potential targets and related pathways of isocuB from a network pharmacology perspective.

Next, the genes and pathways through which isocuB acts were verified. A series of experiments confirmed that isocuB inhibited the proliferation, migration, and invasion of U251 and U87 cells, and induced apoptosis of U251 cells. Furthermore, the inhibitory effect of isocuB on U251 cells was significantly greater than that of TMZ. We hypothesized that it functions through the MMP family. The MMP family is associated with the survival, proliferation, apoptosis, invasion, and metastasis of cancer cells[25]. The WB and RT-qPCR results indicated a significant decrease in the mRNA and protein expression levels of MMP-2 and MMP-9. EMT is a process in which cells acquire invasive mesenchymal potential by undergoing a series of events such as loss of cell junctions, alteration of the cytoskeleton, and remodeling of the extracellular matrix[26]. The results showed a significant reduction in the protein levels of N-cadherin and Vimentin, indicating that isocuB could inhibit EMT and thus inhibit glioma invasion. Several studies have shown that the expression of MMP2/9 is regulated by the STAT3 signaling pathway in a variety of solid tumors, suggesting that increased activation of STAT is the cause of upregulation of MMP2/9 expression in cancer cells[32,33]. Our results further demonstrated that isocuB could impede glioma cell proliferation and migration by inhibiting STAT3 to inhibit MMP2/9. For specific pathways, WB and RT-qPCR experiments were performed for verification. Our study showed that isocuB could inhibit MAPK and STAT3 phosphorylation and PI3K-AKT non-phosphorylated expression, but not PI3K-AKT phosphorylated expression. Our results demonstrated that isocuB could inhibit U251 cells with the modulation of PI3K/AKT and MAPK signaling pathways. To further investigate the mechanism of action of isocuB, we verified the expression of pdk1 and Bcl-2 in vitro using RT-qPCR. The experiments in vitro showed that isocuB downregulated the protein and mRNA expressions by inhibiting PDK1, RXRα, PPARα, and Bcl-2. Different PPARs have significant effects on cancer[34]: PPARα is implicated in promoting proliferation, whereas PPARγ exhibits the opposite effect. Previous studies have linked PPARα to the promotion of the PI3K-AKT pathway[35]. Furthermore, RXRα[36], PDK1[37,38], and Bcl-2[39,40] genes are involved in the PI3K/AKT and MAPK signaling pathways. These results further demonstrated that isocuB acts with the regulation of the PI3K/AKT and MAPK signaling pathways.

Hsa-mir-1268a is upregulated in both cells and exosomes of highly metastatic melanoma CSCs. In addition, the OL-ScS-derived exosome hsa-mir-1268a, after being taken up by OL cells, promoted the transfer and colonization ability of OL cells in vitro and in vivo. Therefore, inhibiting hsa-mir-1286a can inhibit tumor metastasis and colonization[41]. Then, through database search and literature review, we found that hsa-mir-1286a has great research potential in glioma resistance[42], and hsa-mir-1268a was inhibited, and its inhibitors enhanced adriamycin-induced liver cancer cell death[43]. Our RT-qPCR experiment observed a gradual reduction in hsa-mir-1286a RNA content with increasing drug concentration. After inhibiting the expression of hsa-mir-1286a, the drug sensitivity of TMZ in U251 can be effectively improved. Then, we further established the U251/TMZ resistant strains, and the CCK-8 test found that isocuB could inhibit the growth of the resistant strains. In vivo experiments showed that isocuB inhibited tumor growth in nude mice. Our results collectively suggested that isocuB inhibited glioma cell proliferation, migration, and invasion, and induced apoptosis via the PI3K/AKT, MAPK, and STAT3 pathways. Our experiments showed that isocuB could reduce the level of hsa-mir-1286a in glioma U251 cells and increase the drug sensitivity to TMZ.

Although the above results exerted the inhibitory effect of isocuB on glioma, further study could be conducted to explore more detail about the anti-glioma effect and mechanism of isocuB. In the first place, the modulation of isocuB on hub genes in glioma is worthwhile to investigate further. We found the potential core targets of isocuB through network pharmacology and molecular docking, and verified the regulation of isocuB on hub genes via western blot and real-time PCR. However, the direct binding of isocuB to the core target has not been confirmed yet. Cellular thermal shift assay and MicroScale thermophoresis could be utilized for the confirmation of whether isocuB binds directly to the core target. Secondly, exploring the mechanism of isocuB in reversing drug resistance in glioma is necessary. We established the U251/TMZ-resistant strain, and found hsa-mir-1286a might play an important role in sensitizing drug-resistant cells to chemotherapy. Anyway, the combination effect of isocuB and TMZ, and the regulatory effect of hsa-mir-1286a on TMZ-resistant xenograft, as well as PI3K-AKT and MAPK pathways, are still unclear. In this situation, more experiments (e.g., calculating combined index) could be conducted with the combination of isocuB and TMZ in U251/TMZ-resistant strain. Adopting more resistant strains is beneficial for the study. For the animal study, the effect of isocuB on U251/TMZ-resistant model could be conducted in established TMZ-resistant xenograft. Transfection of hsa-mir-1286a mimic or inhibitor and addition of activator of PI3K-AKT and MAPK pathways benefit the investigation of the underlying mechanism of isocuB in TMZ-resistance. Thirdly, more animal models could be used in the study to give a more comprehensive evaluation of the effect of isocuB on glioma. We mainly utilized nude mice to establish the xenograft with a single dose administration. Actually, administration with different dosages in groups, at least dividing into high-dose group, medium-dose group, and low-dose group, could assess the dose-dependent effect of isocuB in vivo. Immunohistochemistry could be utilized to detect the protein expression in vivo and in situ hybridization could exhibit RNA expression, especially for hsa-mir-1286a, in vivo. Establishing a glioma-bearing mouse model via orthotopic transplantation of GL261 cells in C57BL/6 mice is also interesting for further study. Finally, according to our previous experiments, we found that isocuB was slightly toxic to mice. Ensuring the safety of innovative drugs is necessary for preclinical research. Thus, it is notable to conduct in vitro and in vivo experiments for evaluation and investigation of the toxicology of isocuB, and search for novel approaches for detoxification of isocuB, such as structural modification and addition of detoxicant.

Altogether, the underlying mechanism of anti-glioma activity of isocuB and detoxification of isocuB deserve further investigation, providing novel insights into the drug candidate in glioma treatment.

Conclusions

The potential top five genes of the anti-glioma effect of isocuB were identified by network pharmacology as RXRα, AKT1, ESR1, MAPK1, and HSP90AA1. Molecular docking data further indicated that isocuB could bind to the corresponding protein of these genes with a good affinity. Moreover, AKT1 and MAPK1 mRNA were expressed at a high level in glioma, according to the clinical database analysis results. Furthermore, GO and KEGG analyses revealed that the PI3K/Akt and MAPK pathways were the top two connected pathways, as AKT1 and MAPK1 are the key genes of the PI3K/AKT and MAPK pathways. IsocuB exhibited superior efficacy over TMZ in inhibiting glioma proliferation, as evidenced by the CCK8 assay. IsocuB suppressed the invasion and migration of U251 cells by blocking the activity of MMP2 and MMP9 in wound healing and transwell assays. The TUNEL and FCM assays demonstrated the inhibitory effect of isocuB on BCL-2, thereby inducing cell death. Subsequently, we conducted RT-qPCR and WB assays to determine that isocuB could improve the sensitivity of U251 glioma to TMZ by inhibiting hsa-mir-1286a and reduce glioma growth with the modulation of PI3K/AKT and MAPK pathways. Ultimately, our findings illustrate that isocuB exerts an inhibitory effect on glioma growth and EMT process. The mechanism is associated with the regulation of PI3K/AKT and MAPK pathway. Additionally, it is interesting to find that isocuB enhances the drug sensitivity to TMZ via the decrease of hsa-mir-1286a RNA. The comprehensive evidence from experimental demonstration, database analysis, and network pharmacology supports the therapeutic potential of isocuB for glioma in the future.

DECLARATIONS

Authors’ contributions

Completed the network pharmacology and molecular docking experiments and prepared the original draft: Han M

Completed the GO and KEGG enrichment analysis and revised the manuscript: An J

Conducted the experimental verification: Li S, Fan H, Wang L

Contributed to the animal experiment: Du Q

Assisted in conceptualization: Du J

Supervised and revised the manuscript: Yang Y, Song Y

Administrated funding and revised the manuscript: Peng F

This manuscript was conceptualized by all the authors.

All authors have read and agreed to the published version of the manuscript.

Availability of data and materials

The following supporting information can be downloaded.

Financial support and sponsorship

The study was supported by National Natural Science Foundation of China (No. 82003879), Project of Science and Technology Department of Sichuan Province (No. 2023NSFSC1928; 2023NSFSC1992), Young Elite Scientists Sponsorship Program by China Association for Science and Technology (No. CACM-2020-QNRC1-01), and the Fundamental Research Funds for the central Universities.

Conflicts of interest

All authors declared that there are no conflicts of interest.

Ethical approval and consent to participate

All animal experiments followed the regulations of ethics committee of the Experimental Animal Administration of Sichuan University (NO. K2024006).

Consent for publication

Not applicable.

Copyright

© The Author(s) 2024.

REFERENCES

1. Berger TR, Wen PY, Lang-Orsini M, Chukwueke UN. World Health Organization 2021 classification of central nervous system tumors and implications for therapy for adult-type gliomas: a review. JAMA Oncol 2022;8:1493-501.

2. Wang L, Li S, Fan H, et al. Bifidobacterium lactis combined with Lactobacillus plantarum inhibit glioma growth in mice through modulating PI3K/AKT pathway and gut microbiota. Front Microbiol 2022;13:986837.

3. Louis DN, Perry A, Wesseling P, et al. The 2021 WHO classification of tumors of the central nervous system: a summary. Neuro Oncol 2021;23:1231-51.

4. Cai X, Wang M, Mu P, et al. Sequence-defined nanotubes assembled from IR780-conjugated peptoids for chemophototherapy of malignant glioma. Research 2021;2021:9861384.

5. Xu S, Tang L, Li X, Fan F, Liu Z. Immunotherapy for glioma: current management and future application. Cancer Lett 2020;476:1-12.

6. Stupp R, Mason WP, van den Bent MJ, et al. European Organisation for Research and Treatment of Cancer Brain Tumor and Radiotherapy Groups; National Cancer Institute of Canada Clinical Trials Group. Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. Oncol Times 2005;27:15-6.

7. Han M, Li S, Fan H, An J, Peng C, Peng F. Regulated cell death in glioma: promising targets for natural small-molecule compounds. Front Oncol 2024;14:1273841.

8. Tamtaji OR, Razavi ZS, Razzaghi N, Aschner M, Barati E, Mirzaei H. Quercetin and glioma: which signaling pathways are involved? Curr Mol Pharmacol 2022;15:962-8.

9. Liu A, Jiang B, Song C, et al. Isoliquiritigenin inhibits circ0030018 to suppress glioma tumorigenesis via the miR-1236/HER2 signaling pathway. MedComm 2023;4:e282.

10. Buyel JF. Plants as sources of natural and recombinant anti-cancer agents. Biotechnol Adv 2018;36:506-20.

11. Li Y, Li Y, Yao Y, et al. Potential of cucurbitacin as an anticancer drug. Biomed Pharmacother 2023;168:115707.

12. Dai S, Wang C, Zhao X, et al. Cucurbitacin B: a review of its pharmacology, toxicity, and pharmacokinetics. Pharmacol Res 2023;187:106587.

13. Yuan R, Zhao W, Wang QQ, et al. Cucurbitacin B inhibits non-small cell lung cancer in vivo and in vitro by triggering TLR4/NLRP3/GSDMD-dependent pyroptosis. Pharmacol Res 2021;170:105748.

14. Liu JH, Li C, Cao L, Zhang CH, Zhang ZH. Cucurbitacin B regulates lung cancer cell proliferation and apoptosis via inhibiting the IL-6/STAT3 pathway through the lncRNA XIST/miR-let-7c axis. Pharm Biol 2022;60:154-62.

15. Zhang H, Zhao B, Wei H, Zeng H, Sheng D, Zhang Y. Cucurbitacin B controls M2 macrophage polarization to suppresses metastasis via targeting JAK-2/STAT3 signalling pathway in colorectal cancer. J Ethnopharmacol 2022;287:114915.

16. Dandawate P, Subramaniam D, Panovich P, et al. Cucurbitacin B and I inhibits colon cancer growth by targeting the Notch signaling pathway. Sci Rep 2020;10:1290.

17. Yin D, Wakimoto N, Xing H, et al. Cucurbitacin B markedly inhibits growth and rapidly affects the cytoskeleton in glioblastoma multiforme. Int J Cancer 2008;123:1364-75.

18. Nogales C, Mamdouh ZM, List M, Kiel C, Casas AI, Schmidt HHHW. Network pharmacology: curing causal mechanisms instead of treating symptoms. Trends Pharmacol Sci 2022;43:136-50.

19. Pinzi L, Rastelli G. Molecular docking: shifting paradigms in drug discovery. Int J Mol Sci 2019;20:4331.

20. Ren L, Zheng X, Liu J, et al. Network pharmacology study of traditional Chinese medicines for stroke treatment and effective constituents screening. J Ethnopharmacol 2019;242:112044.

21. He B, Zhao Z, Cai Q, et al. miRNA-based biomarkers, therapies, and resistance in cancer. Int J Biol Sci 2020;16:2628-47.

22. WMA. WMA Declaration of Helsinki - Ethical principles for medical research involving human subjects. Available from: https://www.wma.net/policies-post/wma-declaration-of-helsinki-ethical-principles-for-medical-research-involving-human-subjects/. [Last accessed on 26 Apr 2024].

23. Percie du Sert N, Hurst V, Ahluwalia A, et al. The ARRIVE guidelines 2.0: updated guidelines for reporting animal research. PLoS Biol 2020;18:e3000410.

24. Zhou Y, Zhou B, Pache L, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun 2019;10:1523.

25. Jabłońska-Trypuć A, Matejczyk M, Rosochacki S. Matrix metalloproteinases (MMPs), the main extracellular matrix (ECM) enzymes in collagen degradation, as a target for anticancer drugs. J Enzyme Inhib Med Chem 2016;31:177-83.

26. Polyak K, Weinberg RA. Transitions between epithelial and mesenchymal states: acquisition of malignant and stem cell traits. Nat Rev Cancer 2009;9:265-73.

27. Zheng Z, Zhang L, Hou X. Potential roles and molecular mechanisms of phytochemicals against cancer. Food Funct 2022;13:9208-25.

28. Yuan Z, Pan Y, Leng T, et al. Progress and prospects of research ideas and methods in the network pharmacology of traditional Chinese medicine. J Pharm Pharm Sci 2022;25:218-26.

29. Zhou J, Xu N, Liu B, et al. lncRNA XLOC013218 promotes cell proliferation and TMZ resistance by targeting the PIK3R2-mediated PI3K/AKT pathway in glioma. Cancer Sci 2022;113:2681-92.

30. An J, Fan H, Han M, Peng C, Xie J, Peng F. Exploring the mechanisms of neurotoxicity caused by fuzi using network pharmacology and molecular docking. Front Pharmacol 2022;13:961012.

31. Reinhardt A, Stichel D, Schrimpf D, et al. Anaplastic astrocytoma with piloid features, a novel molecular class of IDH wildtype glioma with recurrent MAPK pathway, CDKN2A/B and ATRX alterations. Acta Neuropathol 2018;136:273-91.

32. Yu H, Lee H, Herrmann A, Buettner R, Jove R. Revisiting STAT3 signalling in cancer: new and unexpected biological functions. Nat Rev Cancer 2014;14:736-46.

33. Piperi C, Papavassiliou KA, Papavassiliou AG. Pivotal role of STAT3 in shaping glioblastoma immune microenvironment. Cells 2019;8:1398.

34. Venkata NG, Robinson JA, Cabot PJ, Davis B, Monteith GR, Roberts-Thomson SJ. Mono(2-ethylhexyl)phthalate and mono-n-butyl phthalate activation of peroxisome proliferator activated-receptors alpha and gamma in breast. Toxicol Lett 2006;163:224-34.

35. Leng J, Li H, Niu Y, et al. Low-dose mono(2-ethylhexyl) phthalate promotes ovarian cancer development through PPARα-dependent PI3K/Akt/NF-κB pathway. Sci Total Environ 2021;790:147990.

36. Zhang R, Li H, Zhang S, et al. RXRα provokes tumor suppression through p53/p21/p16 and PI3K-AKT signaling pathways during stem cell differentiation and in cancer cells. Cell Death Dis 2018;9:532.

37. Gagliardi PA, Puliafito A, Primo L. PDK1: at the crossroad of cancer signaling pathways. Semin Cancer Biol 2018;48:27-35.

38. Vara JA, Casado E, de Castro J, Cejas P, Belda-Iniesta C, González-Barón M. PI3K/Akt signalling pathway and cancer. Cancer Treat Rev 2004;30:193-204.

39. Tang J, Yao C, Liu Y, et al. Arsenic trioxide induces expression of BCL-2 expression via NF-κB and p38 MAPK signaling pathways in BEAS-2B cells during apoptosis. Ecotoxicol Environ Saf 2021;222:112531.

40. Hou Y, Wang K, Wan W, Cheng Y, Pu X, Ye X. Resveratrol provides neuroprotection by regulating the JAK2/STAT3/PI3K/AKT/mTOR pathway after stroke in rats. Genes Dis 2018;5:245-55.

41. Li X, Liu D, Chen H, et al. Melanoma stem cells promote metastasis via exosomal miR-1268a inactivation of autophagy. Biol Res 2022;55:29.

42. Li Y, Liu Y, Ren J, et al. miR-1268a regulates ABCC1 expression to mediate temozolomide resistance in glioblastoma. J Neurooncol 2018;138:499-508.

43. Lu YL, Yao JG, Huang XY, et al. Prognostic significance of miR-1268a expression and its beneficial effects for post-operative adjuvant transarterial chemoembolization in hepatocellular carcinoma. Sci Rep 2016;6:36104.

Cite This Article

Export citation file: BibTeX | RIS

OAE Style

Han M, An J, Li S, Fan H, Wang L, Du Q, Du J, Yang Y, Song Y, Peng F. Isocucurbitacin B inhibits glioma growth through PI3K/AKT pathways and increases glioma sensitivity to TMZ by inhibiting hsa-mir-1286a. Cancer Drug Resist 2024;7:16. http://dx.doi.org/10.20517/cdr.2024.01

AMA Style

Han M, An J, Li S, Fan H, Wang L, Du Q, Du J, Yang Y, Song Y, Peng F. Isocucurbitacin B inhibits glioma growth through PI3K/AKT pathways and increases glioma sensitivity to TMZ by inhibiting hsa-mir-1286a. Cancer Drug Resistance. 2024; 7: 16. http://dx.doi.org/10.20517/cdr.2024.01

Chicago/Turabian Style

Han, Mingyu, Junsha An, Sui Li, Huali Fan, Li Wang, Qing Du, Junrong Du, Yuxin Yang, Yuqin Song, Fu Peng. 2024. "Isocucurbitacin B inhibits glioma growth through PI3K/AKT pathways and increases glioma sensitivity to TMZ by inhibiting hsa-mir-1286a" Cancer Drug Resistance. 7: 16. http://dx.doi.org/10.20517/cdr.2024.01

ACS Style

Han, M.; An J.; Li S.; Fan H.; Wang L.; Du Q.; Du J.; Yang Y.; Song Y.; Peng F. Isocucurbitacin B inhibits glioma growth through PI3K/AKT pathways and increases glioma sensitivity to TMZ by inhibiting hsa-mir-1286a. Cancer Drug Resist. 2024, 7, 16. http://dx.doi.org/10.20517/cdr.2024.01

About This Article

Special Issue

This article belongs to the Special Issue Role of Non-Coding RNA in Cancer Drug Resistance
© The Author(s) 2024. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, sharing, adaptation, distribution and reproduction in any medium or format, for any purpose, even commercially, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Data & Comments

Data

Views
277
Downloads
50
Citations
0
Comments
0
6

Comments

Comments must be written in English. Spam, offensive content, impersonation, and private information will not be permitted. If any comment is reported and identified as inappropriate content by OAE staff, the comment will be removed without notice. If you have any queries or need any help, please contact us at support@oaepublish.com.

0
Download PDF
Cite This Article 0 clicks
Like This Article 6 likes
Share This Article
Scan the QR code for reading!
See Updates
Contents
Figures
Related
Cancer Drug Resistance
ISSN 2578-532X (Online)

Portico

All published articles will preserved here permanently:

https://www.portico.org/publishers/oae/

Portico

All published articles will preserved here permanently:

https://www.portico.org/publishers/oae/