Network pharmacology-based strategy to investigate the effect and mechanism of α-solanine against glioma
BMC Complementary Medicine and Therapies volume 23, Article number: 371 (2023)
An anti-tumour activity has been demonstrated for α-solanine, a bioactive compound extracted from the traditional Chinese herb Solanum nigrum L. However, its efficacy in the treatment of gliomas and the underlying mechanisms remain unclear. The aim of this study was to investigate the inhibitory effects of α-solanine on glioma and elucidate its mechanisms and targets using network pharmacology, molecular docking, and molecular biology experiments.
Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP) was utilized to predict the potential targets of α-solanine. GeneCards was used to gather glioma-related targets, and the STRING online database was used to analyze protein–protein interaction (PPI) networks for the shared targets. Hub genes were identified from the resulting PPI network and further investigated using Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Additionally, prognostic and gene set enrichment analyses (GSEA) were carried out to identify potential therapeutic targets and their underlying mechanisms of action in relation to the prognosis of gliomas. In vitro experiments were conducted to verify the findings from the network pharmacology analysis.
A total of 289 α-solanine targets and 1149 glioma-related targets were screened, of which 78 were common targets. 11 hub genes were obtained, including SRC, HRAS, HSP90AA1, IGF1, MAPK1, MAPK14, KDR, STAT1, JAK2, MAP2K1, and IGF1R. The GO and KEGG pathway analyses unveiled that α-solanine was strongly associated with several signaling pathways, including positive regulation of MAP kinase activity and PI3K-Akt. Moreover, α-solanine (10 µM and 15 µM) inhibited the proliferation and migration but promoted the apoptosis of glioma cells. Finally, STAT1 was identified as a potential mediator of the effect of α-solanine on glioma prognosis.
α-Solanine can inhibit the proliferation and migration of gliomas by regulating multiple targets and signalling pathways. These findings lay the foundation for the creation of innovative clinical anti-glioma agents.
Glioma, a primary tumor most frequently in the central nervous system (CNS), has become a significant global public health concern owing to its high mortality rate . Although numerous contemporary treatments are available for glioma, it continues to be one of the most severe afflictions of the human nervous system, marked by a poor prognosis, a propensity for recurrence, and a significant occurrence of disability . If left untreated, patients usually have a survival period of only a few months. As a result, there is a pressing requirement to explore innovative and effective medicinal treatments to combat or manage gliomas.
α-Solanine is a glycoalkaloid that occurs naturally and is synthesized by different plants, including Solanum nigrum L. (Longkui), a kind of traditional Chinese medicine (TCM), as well as several edible plants such as potatoes, cherries, and tomatoes . Numerous studies have indicated that α-solanine offers a multitude of health benefits, including anti-inflammatory properties and the capacity to boost the immune system [4, 5]. It has also been reported to exhibit anti-cancer activity against different types of tumors, such as those present in the breast, liver, pancreas, colon, and other regions . Although the potential of α-solanine for treating glioma has been investigated, its precise mechanism of action remains unclear.
Network pharmacology is a rapidly evolving field of research that aims to unravel the intricacies of disease and drug mechanisms within complex biological networks . This approach has revolutionized investigations by providing a systematic framework to address scientific challenges at multiple levels. In recent years, there has been a growing scholarly emphasis on network pharmacology research with the objective of utilizing Traditional Chinese Medicine (TCM) for the management of tumors. This surge in interest arises from the methodology's capacity to discern regulators possessing multi-targeting capabilities and synergistic effects, thereby elucidating the fundamental mechanisms of TCM in a manner that aligns with the holistic principles of the traditional medical system [8,9,10]. In particular, the accuracy and reliability of this approach deserve recognition . In the present study, we employed a network pharmacology method to determine potential targets and underlying mechanisms by which α-solanine inhibits glioma. The results were subsequently validated through molecular docking and biological experiments. A diagram of the presented workflow is shown in Fig. 1.
Materials and methods
Acquisition target genes of α-solanine and glioma
To identify genes targeted by α-solanine in glioma, we conducted a keyword search for "α-solanine" using the TCMSP (http://tcmspw.com/). To obtain genes related to glioma, we used the GeneCards database (https://www.genecards.org) and searched for the keyword "glioma".
Identifying the common targets between drugs and disease
To identify the common targets between α-solanine and glioma, we applied the intersection method and represented the common targets using a Venn diagram. Furthermore, we constructed a network using Cytoscape v3.6.1 to visualize the connections between the component, common targets and diseases.
Construction of PPI network
The intersection targets were integrated into the STRING database (https://string-db.org/) to obtain information on their protein–protein interaction (PPI) network . We established the screening parameters for the organism as "Homo sapiens" with a minimum required interaction score of "highest confidence (0.9)". Subsequently, the resulting PPI information was visualized using Cytoscape, and a PPI network was constructed.
Cytoscape can calculate various parameters of individual nodes in network diagrams, including degree centrality (DC), betweenness centrality (BC), closeness centrality (CC), and average shortest path length (ASPL), which facilitate in-depth analysis of the characteristics of nodes in an interaction network . In this study, the identification of hub targets for the treatment of glioma with α-solanine was accomplished by selecting nodes that exhibited DC, BC, and CC values exceeding the corresponding median values in the PPI network.
GO and KEGG enrichment analysis
Gene Ontology (GO) enrichment analysis is a common method for identifying biological processes, molecular functions, and cellular components. Meanwhile, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis can reveal significant signalling pathways associated with biological processes. In this study, we conducted GO and KEGG enrichment analyses of the hub targets identified using the DAVID database (https://DAVID.ncifcrf.gov) . The results were further visualized using Bioinformatics software (https://www.bioinformatics.com.cn).
Molecular docking analysis
The application of molecular docking has become essential in computational drug design . In the study, we employed molecular docking to investigate the mode and strength for the interaction between α-solanine and the hub targets. To perform a docking analysis, we adopted CB-Dock (http://cao.labshare.cn/cb-dock/), a recently developed blind docking technology that can autonomously identify the sites for binding to a specific protein and conduct molecular docking using AutoDock . CB-Dock has been proven to surpass other cutting-edge blind docking techniques in accurately forecasting binding locations and binding structures . The crystal structures of the hub protein targets were downloaded from the Protein Data Bank. (http://www.rcsb.org). For the ligand, the 3D structure of α-solanine was retrieved from the PubChem Compound database (https://www.ncbi.nlm.nih.gov). The protein (receptor) and α-solanine (ligand) files were then submitted to the CB-Dock website for molecular docking analysis.
Prognostic and GSEA analysis of hub targets
To reveal the prognostic value of α-solanine in individuals with glioma, we conducted survival analysis on its hub targets. The hub targets were entered into the GEPIA database (http://gepia.cancer-pku.cn) to screen for potential therapeutic targets relevant to the prognosis of glioma patients .
To investigate the mechanisms associated with the differential expression of therapeutic targets in glioma and gain insight into how α-solanine affects the glioma prognosis, we used GSEA software (version 3.0) to examine the biological functions related to the differences in the expression of therapeutic targets between normal and glioma tissues. The reference gene sets were based on the c2 sets of the Molecular Signatures Database (KEGG gene sets,c2.all.v6.2.symbols.gmt), and the permutation number was set at 1,000. The gene sets enriched through GSEA with a nominal P-value below 0.05 and a false discovery rate (FDR) below 0.25 were deemed to be of statistical significance .
The human glioma-derived U87MG cell line, a typical grade IV glioma cell line, was purchased from Wuhan Pu-nuo-sai Life Technology Co. Ltd. in Wuhan, China. This particular cell line is renowned for its vigorous proliferation and tumorigenic properties . The cells were cultivated in RPMI-1640 medium supplemented with 10% fetal bovine serum, 100 µg/mL ampicillin, and 100 µg/mL streptomycin.
Cell morphology visualization
Stock solutions of α-solanine (1 mmol/l) were prepared in dimethyl sulfoxide (DMSO) and diluted with double-distilled water containing 10% DMSO to obtain various concentrations. To investigate the effects of α-solanine on U87MG cells, three different concentrations of the compound (5 µM, 10 µM and 15 µM) were applied to the cells for a period of 24 h. The resulting morphological changes in the cells were then analyzed using a light microscope (Leica, Germany).
Cell viability assay
To evaluate the impact of α-solanine on cell viability, U87MG cells were cultured in 96-well plates at a density of 4 × 103 cells per well and subjected to treatment with α-solanine at varying concentrations of 5 µM, 10 µM and 15 µM for 24 h. The Cell Counting Kit-8 (CCK-8) assay was performed to measure cell proliferation, with absorbance measurements recorded at a wavelength of 450 nm. GraphPad Prism 7.0 software was then used for data analysis and plotting.
U87MG cells were cultured in 6-well plates at a density of 3 × 105 cells per well, allowing them to form a monolayer. To create a wound, three lines were drawn vertically in each well using a sterile pipette aspirator. The cells were then treated with a drug-free medium (serum-free) or a drug-containing medium (serum-free with 10 µM or 15 µM α-solanine) for 24 h. Following treatment, the cells were imaged using a light microscope (Leica, Germany), and the distance between the markers on the cells was measured at both 0 and 24 h. The relative migration distance was calculated using the formula: Relative migration distance (mm) = distance within the scratch at 0 h—distance within the scratch at 24 h. Each experiment was performed in triplicate.
Flow cytometry analysis for apoptosis
U87MG cells were subjected to treatment with α-solanine at concentrations of 10 and 15μΜ for 48 h after being cultured in 6-well plates at 37 °C. To determine the extent of apoptosis, the Annexin V-FITC early apoptosis kit (Shanghai Yeasen Biotechnology Co., Ltd., Shanghai, China) and a flow cytometer (Guava Technologies, Hayward, CA, USA) were used. After washing the treated cells with cold phosphate-buffered saline, they were resuspended in 400 μL of 1 × binding buffer, and 5 μL annexin V-FITC and 5 μL PI staining solution were added. After incubation in the dark, flow cytometry analysis was performed within one hour.
Quantitative Polymerase Chain Reaction (qPCR) assay
The RNA Easy Fast Cell Kit (Tiangen, China) was used to extract total RNA from U87MG cells according to the manufacturer's instructions. The quality of the RNA was evaluated using a SpectraMax Quick Drop reader (Molecular Devices, USA). Reverse transcription was performed using the iScript cDNA synthesis kit (Bio-Rad) on 1 μg of RNA. RT-qPCR was conducted using SYBR Green Real-Time PCR Master Mix (Toyobo, Japan) to measure the mRNA expression levels of four genes: STAT1, P53, BAX, and BCL-2. The β-actin gene served as the internal control, and the expression of each gene was calculated using the 2−ΔΔCt method. The primer sequences used for qPCR analysis are listed in Additional file 1.
Continuous data that followed a normal distribution were represented as the mean value ± the standard deviation (SD). GraphPad Prism 7.0 software was used for conducting statistical analysis, and the two-tailed unpaired Student's t-test was used to determine differences between the groups. Statistical significance was indicated by an asterisk (*: P < 0.05) or two asterisks (**: P < 0.01).
Identifying gene targets for α-Solanine in Glioma
Based on the TCMSP database, we identified 290 genes that are targeted by α-solanine. To ensure consistency in gene names and species information, we standardized the data using the UniProt database, leading to a final set of 289 genes. Simultaneously, 1149 disease-related genes were obtained from the GeneCards database (Additional file 2). The two datasets were merged for co-analysis, resulting in the identification of 78 overlapping genes, representing 5.7% of the total genes analyzed, as shown in Fig. 2a.
Construction of Compound-Target-Disease (CTD) network
The pharmacological potential of α-solanine for glioma treatment was visualized using the Cytoscape software to construct a network of α-solanine-common targets-glioma (Fig. 2b). The network consists of 78 nodes representing common target genes (depicted in green) and two supplemental nodes (depicted in yellow) corresponding to α-solanine and glioma. The network diagram showed 155 edges representing the connections between the nodes. The findings indicate that α-solanine exerts a inhibitory effect on glioma cell proliferation by acting on multiple targets.
Construction and analysis of PPI network
PPI analysis was performed for the 78 genes commonly targeted by α-solanine based on data in the STRING database (Fig. 3). The resulting TSV data were imported into Cytoscape, which generated a network diagram consisting of 78 nodes and 885 edges, as presented in Fig. 4a. The size and colour intensity are directly proportional to the degree value and target probability. Finally, 11 hub targets were identified based on their DC, BC, and CC values (Fig. 4b and Table 1). A list of attribute values for the remaining 67 common targets is included in Additional file 3.
GO and KEGG pathway enrichment analysis
A total of 1680 GO terms were enriched, with 1525 belonging to BP, 37 to CC, and 118 to MF. The enriched BP terms were predominantly related to the regulation of phosphatidylinositol 3-kinase signalling, response to reactive oxygen species, and positive regulation of MAP kinase activity. The CC terms were mainly related to membrane raft, membrane microdomain, and membrane region. Moreover, the MF category had enriched terms in protein tyrosine kinase activity, MAP kinase activity, integrin binding, and protein phosphatase binding. The top 10 GO terms ranked by their adjusted p-values are displayed in Fig. 4c. The terms with lower adjusted p-values are coloured red, and this indicate higher enrichment.
KEGG pathway analysis revealed 119 pathways, comprising 11 disease signalling pathways and 39 molecular signalling pathways among the top 50 enriched KEGG pathways. Figure 4d demonstrates the top 10 disease signalling pathways with high counts, including the PI3K-Akt signalling pathway, mTOR signalling pathway, apoptosis, microRNAs in cancer, glioma, PD-L1 expression, and PD-1 checkpoint pathway in cancer. Our findings suggest that α-solanine may exert its therapeutic effects on glioma through various potential biological mechanisms.
Verification of compound‑target interaction
We selected 11 hub target proteins for molecular docking to verify the interaction between α-solanine and its potential targets. The structure of α-solanine was then uploaded to CB-Dock to analyse its potential for docking with SRC, HRAS, HSP90AA1, IGF1, MAPK1, MAPK14, KDR, STAT1, JAK2, MAP2K1, and IGF1R. The docking scores for each target protein were documented in Table 2. According to the Vina scoring system, a score of -5.0 kcal/mol or lower is indicative of favorable binding interactions between ligands and receptors . Therefore, the molecular docking analysis revealed a significant affinity between α-solanine and the 11 hub target proteins. These findings not only validate the accuracy of the selected hub targets but also suggest the potential effectiveness of α-solanine in glioma treatment. The docking sketches of target proteins with α-solanine are presented in Fig. 5 and Additional file 4.
To display the active site crystal structure of the protein, we utilized a color-coding scheme where carbon atoms are represented by white, oxygen atoms by red, nitrogen atoms by blue, and sulfur atoms by yellow. Meanwhile, the crystal pose of the ligand is depicted using a different color scheme, where hydrogen atoms are white, carbon atoms are grey, and oxygen atoms are red.
Survival analysis and GSEA analysis of the key targets
An investigation was conducted utilizing data from the GEPIA database to examine the expression and prognostic relevance of 11 hub targets in glioma. It was observed that there was a significant increase in STAT1 mRNA transcription in LGG and GBM tissues compared to normal tissues, according to our findings. Moreover, STAT1 was the only hub target significantly correlated with the overall survival (OS) and disease-free survival (DFS) of glioma patients (Fig. 6a-c).
To better understand how α-solanine exerts its therapeutic effects through the hub target STAT1 and improves the prognosis of glioma patients, we conducted a GSEA analysis to identify potential biological signalling pathways affected by the overexpression of STAT1 in glioma. As depicted in Fig. 6d, our results revealed an inverse correlation between STAT1 overexpression in glioma cells and the p53 pathway, which has been identified as a beneficial factor for enhancing the overall outlook of patients with glioma .
α-Solanine suppressed the proliferation of glioma cells
In order to investigate the effects of α-solanine on glioma cell proliferation, U87MG cells were treated with various concentrations of α-solanine. The morphology of cells was examined using an optical microscope. It was observed that α-solanine had no significant effect on the viability of U87MG cells when administered at a concentration of 5 µM (Fig. 7a–b). Consequently, only the concentrations of 10 µM and 15 µM were employed in all subsequent experimental procedures. At a concentration of 10 µM, we observed a reduction in the number of glioma cells, an increase in cell spacing, and a change in cell shape towards an oval and contracted morphology. When the α-solanine concentration was further increased to 15 µM, we observed pronounced cell wall wrinkling and nuclear consolidation, accompanied by a marked decrease in the glioma cell count and the appearance of turbidity in the cell culture medium (Fig. 7c-d).
Further statistical analysis showed that concentrations of 10 and 15 μM significantly inhibited the proliferation of glioma cells (Fig. 7e). Nonlinear regression analysis of the data obtained from the CCK-8 experiments was then performed (GraphPad Prism 7.0) to obtain the IC50 valuesof α-solanine against glioma cells (Fig. 7f). Overall, these results suggested that α-solanine suppresses the proliferation of glioma cells in a dose-dependent manner.
α-Solanine inhibited migration of glioma cells
In vitro cell scratching experiments were conducted to investigate the effect of α-solanine on the migratory ability of U87MG cells. The results, including the initial scratch displacement for each cohort, as well as the inter-group scratch distance after 24 h with or without α-solanine treatment, are shown in Fig. 8. We calculated the relative migration distances of glioma cells in each group and compared the findings between the α-solanine treatment and the control group. The results of the statistical analysis clearly revealed that α-solanine concentrations of 10 µM and 15 µM effectively decreased the relative migration distances of glioma cells.
α-Solanine promoted apoptosis of glioma cells
The effect of α-solanine on apoptosis in U87MG cells was investigated using flow cytometry. Our results showed that treatment with α-solanine at concentrations of 10 µM and 15 µM significantly increased early apoptosis (Fig. 9), demonstrating the efficacy of α-solanine in inducing apoptosis of U87MG cells.
Effect of α-solanine on mRNA expression levels of STAT1, P53, BAX, and BCL-2 in glioma cells
Real-time PCR analysis was performed on U87MG cells to validate the GSEA analysis results. Treatment with α-solanine significantly decreased the transcription of STAT1mRNA, but markedly enhanced P53mRNA expression levels (Fig. 10). In addition, we observed a considerable increase in the ratio of Bax mRNA to Bcl-2 mRNA after α-solanine treatment.
Glioma is a highly aggressive neoplasm that originates from glial cells, with a poor prognosis, and it commonly affects the brain and spinal cord. Standard treatment modalities for glioma include surgical resection, chemotherapy, and radiation therapy . Despite these options, a significant proportion of patients experience recurrence within 6–10 months of surgical treatment, resulting in an overall 5-year survival rate of less than 3% . Recent evidence indicates that even with standardized treatment, including surgery, chemotherapy, and radiation therapy, the 5-year survival rate for glioma patients remains below 5% , which has led to the notion that gliomas are incurable .
Traditional Chinese Medicine (TCM) has a long history of application in treating various diseases and is now considered a complementary or alternative approach to conventional medicine . In contemporary times, it has witnessed significant advancements, leading to an increase in the use of herbal medications to treat tumors . Solanum ni grum Linn. (Longkui) is a widely used traditional Chinese medicinal herb with anticancer properties. Its main steroidal glycoalkaloid, α-solanine, is highly effective against various types of cancer in all three stages: carcinogenesis, growth, and metastasis [4, 30, 31]. These studies have offered valuable indications for detecting the correlation between α-solanine and glioma. In this study, we applied network pharmacology to identify the possible targets of α-solanine and elucidate its underlying mechanisms against glioma.
The present study employed a systematic approach to identify potential target genes for α-solanine treatment of glioma based on relevant databases and Venn diagrams. Subsequently, 78 candidate genes were identified, and their interactions were visualized using protein–protein interaction (PPI) analysis. Through further network analysis using Cytoscape,11 hub target genes, including SRC, HRAS, HSP90AA1,IGF1, MAPK1, MAPK14, KDR, STAT1, JAK2, MAP2K1, and IGF1R, were selected for subsequent investigation. Molecular docking analysis was performed to determine whether α-solanine binds with the hub targets. The results revealed that α-solanine exhibited high binding affinity with the 11 core target genes, suggesting its potential as a therapeutic agent for glioma.
SRCs are critical members of the SRC family of protein tyrosine kinases and play a crucial role in various cellular processes, including morphology, motility, proliferation, and survival . Inhibiting SRC and upregulating PTEN significantly reduces the migration and invasion of glioma stem cells , leading to a growing interest in developing novel SRC inhibitors for glioma treatment . Additionally, the expression levels of HRAS, MAPK1, IGF1, and IGF1R in glioma cells have been closely linked to the growth, invasion, and metastasis of glioma [35,36,37,38]. Previous studies have identified MAPK14 (p38) and KDR as potential tumour suppressors in glioma development and potential antigens for vaccine advancement [39, 40]. HSP90AA1 is consistent with MAPK14 in playing a vital role in regulating apoptosis . Therefore, targeting these hub targets may have important implications for the treatment of glioma.
The GO analysis revealed that α-solanine participates in several critical biological processes in glioma, including the regulation of reactive oxygen species (ROS) response, mitogen-activated protein (MAP) kinase activity, and phosphatidylinositol 3-kinase (PI3K) signalling, among others. ROS plays a crucial role in inducing apoptosis of tumour cells, and various medicinal plants have been found to increase ROS levels to induce cancer cell apoptosis [3, 42, 43]. Moreover, both MAP kinase activity and PI3K signalling are associated with the mitogenic process , and the cessation of mitosis can trigger apoptosis in tumorigenic cells .
The KEGG pathway analysis suggested that α-solanine inhibits the proliferation of glioma through several pathways, such as the Ras signalling pathway, PD-L1 expression and the PD-1 checkpoint pathway, the PI3K-Akt signalling pathway, MAPK signalling pathways, and VEGF signalling pathways. The Ras signalling pathway is an essential component of cellular signalling pathways and plays a crucial role in oncogenic signalling . The PD-1 and PD-L1 pathways serve as a crucial checkpoint in the immune system, regulating its activity and preventing over activation. Recent studies have shown that blocking these pathways could suppress the progression of several cancers . The PI3K-Akt pathway has been the subject of extensive research in tumorigenesis and has been implicated in various aspects of glioma development, including formation, migration, invasion, and apoptosis . Studies have also demonstrated that Girdin, an actin-binding protein, regulates the migration and invasion of glioma cells through the PI3K-Akt signalling pathway . Additionally, the signalling pathways of MAPK and VEGF have crucial functions in the processes of glioma cell migration, proliferation, and angiogenesis [50, 51].
STAT1, as a transcription factor, plays a significant role in cytokine-mediated signalling pathways that have a range of biological effects such as antiviral activities, suppression of tumour cell growth, and promotion of apoptosis . Although STAT1 is commonly regarded as a tumour suppressor, several studies have highlighted its oncogenic role in the pathogenesis of tumours [53, 54]. A recent study reported that STAT1 has a protumorigenic role in gliomas . Additionally, high expression of STAT family members predicts poor glioma prognosis . In agreement with these studies, our findings showed that STAT1 was upregulated in glioma specimens compared with matched normal tissue specimens, and its high expression was associated with poor OS and DFS in glioma patients, supporting its potential as a reliable prognostic indicator and a promising therapeutic target for glioma.
Despite the well-established knowledge of solanine's inhibitory effects on multiple signaling pathways such as NF-κB, ERK1/2, AKT and STAT1, its precise effect on gliomas through STAT1 remains uncertain . According to the GSEA results, STAT1 was found to negatively regulate the p53 pathway in glioma cells. Two in vitro experiments conducted on human fibrosarcoma cells and lymphocytes confirmed this finding that STAT1 in tumor cells has a regulatory effect on the p53 pathway, demonstrating that the overexpression of STAT1 suppressed MDM2 and upregulated P53 expression [58, 59]. We hypothesized that α-solanine could decrease STAT1 expression in glioma cells, which, in turn, could activate the p53 signalling pathway. Subsequently, our qRT-PCR results validated this hypothesis, which showed a significant decrease in STAT1 expression and a substantial increase in P53 expression with an increase in α-solanine concentration. Furthermore, considering that the Bax/Bcl-2 pathway is a pro-apoptotic pathway positively regulated by P53 [60, 61], the significant increase in the Bax/Bcl-2 ratio following α-solanine intervention indicates that p53 and its downstream signaling pathways were activated, inducing U87MG cell apoptosis. Therefore, we conclude that α-solanine enhances the activity of the P53 signalling pathway by down-regulating STAT1 expression, initiating apoptosis of glioma cells, and ultimately contributing to a favorable prognosis for glioma patients.
In recent years, network pharmacology, which combines principles from systems and computational biology, has emerged as a valuable research tool for investigating the network regulatory effects of TCM or formulas in the context of various diseases . Molecular docking is another powerful tool that can validate the accuracy of network pharmacology results and predict biological experiments by simulating the interaction between small molecules and biological macromolecules [62, 63]. Therefore, the amalgamation of network pharmacology and molecular docking trials can amplify the efficacy of TCM investigation. In this study, we also applied experimental methods to explore the possible mechanisms of α-solanine against glioma at the molecular and cellular levels, complementing the results of network pharmacology. However, our research has some limitations. First, some potential α-solanine targets may have been overlooked due to the inherent challenges in obtaining relevant information from the Network Pharmacology database. Second, while GO and KEGG enrichment analysis predicted multiple signaling pathways targeted by α-solanine in glioma cells, the detailed mechanisms require further investigation. Third, although our findings provide a foundation for future research into the potential treatment of gliomas with α-solanine, further investigation are necessary to translate these findings into clinical applications.
α-Solanine is a promising and safe TCM for the treatment of glioma. Our network pharmacological analysis revealed that α-solanine exerts an anti-glioma effect by modulating the expression or activity of multiple targets, pathways and biological processes crucial for cellular metabolism and apoptosis. Specifically, α-solanine improves the glioma prognosis by down-regulating the expression of STAT1. Further exploration of the potential benefits of α-solanine glioma treatment is warranted to validate our findings.
Dolecek TA, Propp JM, Stroup NE, Kruchko C. Cbtrus statistical report: primary brain and central nervous system tumors diagnosed in the united states in 2005–2009. Neuro Oncol. 2012;14 Suppl 5(Suppl 5):1–49.
Davis ME. Epidemiology and overview of gliomas. Semin Oncol Nurs. 2018;34(5):420–9.
Meng XQ, Zhang W, Zhang F, Yin SY, Xie HY, Zhou L, et al. Solanine-induced reactive oxygen species inhibit the growth of human hepatocellular carcinoma hepg2 cells. Oncol Lett. 2016;11(3):2145–51.
Lv C, Kong H, Dong G, Liu L, Tong K, Sun H, et al. Antitumor efficacy of alpha-solanine against pancreatic cancer in vitro and in vivo. PLoS One. 2014;9(2):e87868.
Gao J, Ying Y, Wang J, Cui Y. Solanine inhibits immune escape mediated by hepatoma treg cells via the tgfbeta/smad signaling pathway. Biomed Res Int. 2020;2020:9749631.
Zhong WF, Liu SP, Pan B, Tang ZF, Zhong JG, Zhou FJ. Solanine inhibits prostate cancer du145 xenograft growth in nude mice by inducing cell cycle arrest in g1/s phase. Nan Fang Yi Ke Da Xue Xue Bao. 2016;36(5):665–70.
Fan L, Guo Y, Wu Q, Hu T, Chen X, Guo J, et al. Mechanism of xiezhuo huayu yiqi tongluo formula in the treatment of uric acid nephropathy based on network pharmacology, molecular docking, and in vivo experiments. Evid Based Complement Alternat Med. 2023;2023:6931644.
Choudhary N, Singh V. Insights about multi-targeting and synergistic neuromodulators in ayurvedic herbs against epilepsy: integrated computational studies on drug-target and protein-protein interaction networks. Sci Rep. 2019;9(1):10565.
Choudhary N, Choudhary S, Kumar A, Singh V. Deciphering the multi-scale mechanisms of tephrosia purpurea against polycystic ovarian syndrome (pcos) and its major psychiatric comorbidities: studies from network pharmacological perspective. Gene. 2021;773:145385.
Choudhary N, Singh V. Multi-scale mechanism of antiviral drug-alike phytoligands from ayurveda in managing covid-19 and associated metabolic comorbidities: insights from network pharmacology. Mol Divers. 2022;26(5):2575–94.
Hopkins AL. Network pharmacology: the next paradigm in drug discovery. Nat Chem Biol. 2008;4(11):682–90.
Qi YS, Xiao MY, Xie P, Xie JB, Guo M, Li FF, et al. Comprehensive serum metabolomics and network analysis to reveal the mechanism of gypenosides in treating lung cancer and enhancing the pharmacological effects of cisplatin. Front Pharmacol. 2022;13:1070948.
Du J, Zhang L. Pathway deviation-based biomarker and multi-effect target identification in asbestos-related squamous cell carcinoma of the lung. Int J Mol Med. 2017;39(3):579–86.
Li S, Zhang B. Traditional chinese medicine network pharmacology: theory, methodology and application. Chin J Nat Med. 2013;11(2):110–20.
Liu C, Li H, Wang K, Zhuang J, Chu F, Gao C, et al. Identifying the antiproliferative effect of astragalus polysaccharides on breast cancer: coupling network pharmacology with targetable screening from the cancer genome atlas. Front Oncol. 2019;9:368.
Liu Y, Yang X, Gan J, Chen S, Xiao ZX, Cao Y. CB-dock2: improved protein-ligand blind docking by integrating cavity detection, docking and homologous template fitting. Nucleic Acids Res. 2022;50(W1):W159–64.
Liu Y, Grimm M, Dai WT, Hou MC, Xiao ZX, Cao Y. Cb-dock: a web server for cavity detection-guided protein-ligand blind docking. Acta Pharmacol Sin. 2020;41(1):138–44.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. Gepia: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–102.
Lin Y, Qu L, Wu J, Pu M, Huang Y, Cao Y. Identification of adipogenesis subgroups and immune infiltration characteristics in diabetic peripheral neuropathy. J Immunol Res. 2023;2023:3673094.
Qu M, Yu J, Liu H, Ren Y, Ma C, Bu X, et al. The candidate tumor suppressor gene slc8a2 inhibits invasion, angiogenesis and growth of glioblastoma. Mol Cells. 2017;40(10):761–72.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.
Wang D, Zhang Y, Wang X, Zhang L, Xu S. Construction and validation of an aging-related gene signature predicting the prognosis of pancreatic cancer. Front Genet. 2023;14:1022265.
Pollack IF, Finkelstein SD, Woods J, Burnham J, Holmes EJ, Hamilton RL, et al. Expression of p53 and prognosis in children with malignant gliomas. N Engl J Med. 2002;346(6):420–7.
Li W, Zhang R, Yang J, Wang R. Efficacy and prognosis of surgery combined with (125)i seed implantation in treatment of recurrent glioma. Oncol Lett. 2017;14(6):7201–6.
Alapati K, Kesanakurti D, Rao JS, Dasari VR. Upar and cathepsin b-mediated compartmentalization of jnk regulates the migration of glioma-initiating cells. Stem Cell Res. 2014;12(3):716–29.
Huang L, Li X, Ye H, Liu Y, Liang X, Yang C, et al. Long non-coding rna nck1-as1 promotes the tumorigenesis of glioma through sponging microrna-138-2-3p and activating the trim24/wnt/beta-catenin axis. J Exp Clin Cancer Res. 2020;39(1):63.
Wang X, Zhao Y, Dong S, Lee RJ, Yang D, Zhang H, et al. Correction: wang et al. Cell-penetrating peptide and transferrin co-modified liposomes for targeted therapy of glioma. Molecules. 2019, 24, 3540. Molecules. 2022;27(17):5729.
Zhou W, Wang J, Wu Z, Huang C, Lu A, Wang Y. Systems pharmacology exploration of botanic drug pairs reveals the mechanism for treating different diseases. Sci Rep. 2016;6:36985.
Meng C, Teng Y, Jiang X. Raddeanin a induces apoptosis and cycle arrest in human hct116 cells through pi3k/akt pathway regulation in vitro and in vivo. Evid Based Complement Alternat Med. 2019;2019:7457105.
Shen KH, Liao AC, Hung JH, Lee WJ, Hu KC, Lin PT, et al. Alpha-solanine inhibits invasion of human prostate cancer cell by suppressing epithelial-mesenchymal transition and mmps expression. Molecules. 2014;19(8):11896–914.
Yan X, Li M, Chen L, Peng X, Que ZJ, An HM, et al. Alpha-solanine inhibits growth and metastatic potential of human colorectal cancer cells. Oncol Rep. 2020;43(5):1387–96.
Yang J, Wang L, Wu MX. 830 nm photobiomodulation therapy promotes engraftment of human umbilical cord blood-derived hematopoietic stem cells. Sci Rep. 2020;10(1):19671.
Jaraiz-Rodriguez M, Tabernero MD, Gonzalez-Tablas M, Otero A, Orfao A, Medina JM, et al. A short region of connexin43 reduces human glioma stem cell migration, invasion, and survival through src, pten, and fak. Stem Cell Reports. 2017;9(2):451–63.
Yang L, Li H, Yang M, Zhang W, Li M, Xu Y, et al. Exploration in the mechanism of kaempferol for the treatment of gastric cancer based on network pharmacology. Biomed Res Int. 2020;2020:5891016.
Xie J, Sun J, Feng J, Yang F, Wang J, Wen T, et al. Kernel differential subgraph analysis to reveal the key period affecting glioblastoma. Biomolecules. 2020;10(2):318.
Tang D, Gao W, Yang J, Liu J, Zhao J, Ge J, et al. Mir-181d promotes cell proliferation via the igf1/pi3k/akt axis in glioma. Mol Med Rep. 2020;22(5):3804–12.
Yu M, Yu S, Gong W, Chen D, Guan J, Liu Y. Knockdown of linc01023 restrains glioma proliferation, migration and invasion by regulating igf-1r/akt pathway. J Cancer. 2019;10(13):2961–8.
Li Q, Feng Y, Chao X, Shi S, Liang M, Qiao Y, et al. Hotair contributes to cell proliferation and metastasis of cervical cancer via targetting mir-23b/mapk1 axis. Biosci Rep. 2018;38(1):BSR20171563.
Tang T, Wang LX, Yang ML, Zhang RM. Lncrna tptep1 inhibits stemness and radioresistance of glioma through mir-106a-5p-mediated p38 mapk signaling. Mol Med Rep. 2020;22(6):4857–67.
Zhou Q, Yan X, Zhu H, Xin Z, Zhao J, Shen W, et al. Identification of three tumor antigens and immune subtypes for mrna vaccine development in diffuse glioma. Theranostics. 2021;11(20):9775–90.
Sheng J, Zhao J, Xu Q, Wang L, Zhang W, Zhang Y. Bioinformatics analysis of srsf1-controlled gene networks in colorectal cancer. Oncol Lett. 2017;14(5):5393–9.
Li T, Sun W, Dong X, Yu W, Cai J, Yuan Q, et al. Total ginsenosides of chinese ginseng induces cell cycle arrest and apoptosis in colorectal carcinoma ht-29 cells. Oncol Lett. 2018;16(4):4640–8.
Eom KS, Kim HJ, So HS, Park R, Kim TY. Berberine-induced apoptosis in human glioblastoma t98g cells is mediated by endoplasmic reticulum stress accompanying reactive oxygen species and mitochondrial dysfunction. Biol Pharm Bull. 2010;33(10):1644–9.
Jiao P, Feng B, Xu H. Mapping mkp-3/foxo1 interaction and evaluating the effect on gluconeogenesis. PLoS One. 2012;7(7):e41168.
Zhao J, Liu N, Hao X, Que L, Liu J, Tang X. Association between integrin-linked kinase and hyperthermia in oral squamous cell carcinoma. Oncol Lett. 2017;14(6):7705–14.
Dunnett-Kane V, Burkitt-Wright E, Blackhall FH, Malliri A, Evans DG, Lindsay CR. Germline and sporadic cancers driven by the ras pathway: parallels and contrasts. Ann Oncol. 2020;31(7):873–83.
Wang Z, Zhang C, Liu X, Wang Z, Sun L, Li G, et al. Molecular and clinical characterization of pd-l1 expression at transcriptional level via 976 samples of brain glioma. Oncoimmunology. 2016;5(11):e1196310.
Tang Z, He Z. Tigar promotes growth, survival and metastasis through oxidation resistance and akt activation in glioblastoma. Oncol Lett. 2019;18(3):2509–17.
Ni W, Fang Y, Tong L, Tong Z, Yi F, Qiu J, et al. Girdin regulates the migration and invasion of glioma cells via the pi3k-akt signaling pathway. Mol Med Rep. 2015;12(4):5086–92.
Zhang Y, Xu J, Zhu X. A 63 signature genes prediction system is effective for glioblastoma prognosis. Int J Mol Med. 2018;41(4):2070–8.
Nalluri JJ, Barh D, Azevedo V, Ghosh P. Mirsig: a consensus-based network inference methodology to identify pan-cancer mirna-mirna interaction signatures. Sci Rep. 2017;7:39684.
Verhoeven Y, Tilborghs S, Jacobs J, De Waele J, Quatannens D, Deben C, et al. The potential and controversy of targeting stat family members in cancer. Semin Cancer Biol. 2020;60:41–56.
Wu S, Wu Y, Lu Y, Yue Y, Cui C, Yu M, et al. Stat1 expression and hpv16 viral load predict cervical lesion progression. Oncol Lett. 2020;20(4):28.
Li X, Wang Z, Zhang S, Yao Q, Chen W, Liu F. Ruxolitinib induces apoptosis of human colorectal cancer cells by downregulating the jak1/2-stat1-mcl-1 axis. Oncol Lett. 2021;21(5):352.
Guo G, Wang S, Hao Y, Ren Y, Wu Y, Zhang J, et al. Stat1 is a modulator of the expression of frequently rearranged in advanced t-cell lymphomas 1 expression in u251 cells. Oncol Lett. 2020;20(1):248–56.
Ji W, Liu Y, Xu B, Mei J, Cheng C, Xiao Y, et al. Bioinformatics analysis of expression profiles and prognostic values of the signal transducer and activator of transcription family genes in glioma. Front Genet. 2021;12:625234.
Zhao L, Wang L, Di SN, Xu Q, Ren QC, Chen SZ, et al. Steroidal alkaloid solanine a from solanum nigrum linn. Exhibits anti-inflammatory activity in lipopolysaccharide/interferon gamma-activated murine macrophages and animal models of inflammation. Biomed Pharmacother. 2018;105:606–15.
Townsend PA, Scarabelli TM, Davidson SM, Knight RA, Latchman DS, Stephanou A. Stat-1 interacts with p53 to enhance dna damage-induced apoptosis. J Biol Chem. 2004;279(7):5811–20.
Baran-Marszak F, Feuillard J, Najjar I, Le Clorennec C, Bechet JM, Dusanter-Fourt I, et al. Differential roles of stat1alpha and stat1beta in fludarabine-induced cell cycle arrest and apoptosis in human b cells. Blood. 2004;104(8):2475–83.
Liu J, Wang J, Ning Y, Chen F. The inhibition of mir-101a-3p alleviates h/r injury in h9c2 cells by regulating the jak2/stat3 pathway. Mol Med Rep. 2020;21(1):89–96.
Shi Y. Meg3 regulates apoptosis of adipose-derived stem cells. Mol Med Rep. 2020;21(6):2435–42.
Li Y, Zhang C, Ma X, Yang L, Ren H. Identification of the potential mechanism of radix pueraria in colon cancer based on network pharmacology. Sci Rep. 2022;12(1):3765.
Sharma R, Jadhav M, Choudhary N, Kumar A, Rauf A, Gundamaraju R, et al. Deciphering the impact and mechanism of trikatu, a spices-based formulation on alcoholic liver disease employing network pharmacology analysis and in vivo validation. Front Nutr. 2022;9:1063118.
We thank Bullet Edits (http://www.bulletedits.cn/) for editing this manuscript.
Ethics approval and consent to participate
We confirm that all experiments followed the Helsinki Declaration. Since the in vitro experiments were conducted on cell lines, the ETHICS Committee did not have to approve this study.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
PCR primers and their sequences used in this study.
Targets of disease and drug action and their intersection.
The attribute values of 78 common targets.
The key amino acid residues of STAT1interact with α-solanine. Blue dashed lines represent hydrogen bonds, gray dashed lines represent hydrophobic interactions, and yellow dashed lines represent salt bridges.
About this article
Cite this article
Wang, C., Liu, X. & Guo, S. Network pharmacology-based strategy to investigate the effect and mechanism of α-solanine against glioma. BMC Complement Med Ther 23, 371 (2023). https://doi.org/10.1186/s12906-023-04215-1