Abstract
Background:
Glioma is a highly aggressive brain tumor with a poor prognosis. Recent studies have demonstrated the anti-tumor potential of Liujunzi decoction (LJZD), but its specific effects and mechanisms in glioma remain unclear. This study aimed to elucidate the anti-glioma effects of LJZD and explore its underlying therapeutic mechanisms.
Methods:
The active ingredients of LJZD were retrieved from the TCMSP database and screened using ADME analysis. Potential targets of both the ingredients and glioma were retrieved from various databases. The overlapping targets were identified as LJZD’s therapeutic targets for glioma. To elucidate the biological functions and potential association among the overlapping targets, protein-protein interaction (PPI) network analysis, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were conducted. Molecular docking was employed to evaluate the binding affinities between core targets and LJZD’s main active ingredients. The inhibitory effect of LJZD on U251 glioma cells was validated by in vitro experiments.
Results:
In total, we identified 76 active ingredients and 752 potential targets of LJZD, while 1,456 glioma-related targets were collected from databases, yielding 174 overlapping targets through intersection. PPI network analysis revealed 40 hub targets, with TP53, AKT1, HIF1A, TNF and IL6 ranking as the top five core targets based on degree. GO and KEGG enrichment analyses demonstrated that LJZD’s anti-glioma effects are mediated through genes related to apoptosis and key pathways including platinum drug resistance, EGFR tyrosine kinase inhibitor resistance, and PD-LI expression and PD-1 checkpoint pathway in cancer. Molecular docking confirmed good to strong binding affinities between LJZD’s main active ingredients (glabridin, pinocembrin and chrysophanol) and all top 10 core targets. In vitro experiments showed that glabridin inhibited U251 cell viability and promoted apoptosis. Glabridin significantly downregulated phosphorylation levels of AKT, STAT3 and c-Jun, while reducing expression of c-Myc and HIF-1α in U251 cells.
Conclusion:
LJZD may exert anti-glioma effects by inducing tumor cell apoptosis and overcoming drug resistance, potentially serving as an adjuvant to enhance conventional therapies. Additionally, its active ingredients (e.g., glabridin, pinocembrin and chrysophanol) provide novel leads for developing small-molecule therapeutics.
1 Introduction
Glioma is a common malignant intracranial tumor, accounting for 80% of all malignant intracranial tumors (Schaff and Mellinghoff, 2023). Current treatment options primarily include surgery, radiotherapy, and chemotherapy, with the addition of tumor-treating fields showing some survival benefit (Stupp et al., 2017). Even with combined surgery and postoperative radiochemotherapy, the median survival is only around 14 months (Stupp et al., 2005; Stupp et al., 2009). Although immunotherapy and targeted therapy have recently emerged as potential treatments for glioma, the overall prognosis for patients remains poor, with many tumors eventually recurring. Therefore, the development of more effective glioma treatment strategies remains a key area of research.
Traditional Chinese medicine (TCM) and formulas are a treasure of Chinese civilization, with a history of thousands of years of application in China. Liujunzi decoction (LJZD), derived from the “Yixue Zhengzhuan”, is composed of six Chinese herbal ingredients, which are Panax Ginseng C. A. Mey., Atractylodes Macrocephala Koidz., Poria Cocos (Schw.) Wolf., licorice, Arum Ternatum Thunb., and Citrus Reticulata. It is based on the four ingredients of Sijunzi decoction with the addition of Arum Ternatum Thunb. and Citrus Reticulata. Therefore, in addition to its original functions of invigorating Qi and strengthening the spleen, LJZD also has the effects of drying dampness, resolving phlegm, regulating Qi and harmonizing the stomach. Historically, LJZD has been used to treat syndromes of spleen-stomach qi deficiency accompanied by phlegm-dampness or qi stagnation. Recent studies have revealed its potential antitumor effects. For instance, according to , LJZD has the ability to enhance the expression of miR-122-3p in hepatocellular carcinoma cells. This action leads to the downregulation of UBE2I and the inhibition of the NF-κB/PD-L1 pathway, effectively suppressing SUMOylation in cancer cells and safeguarding PD-1+ T cells co-cultured with tumor cells from undergoing apoptosis. Additionally, it has been demonstrated that LJZD can suppress the growth of esophageal cancer cells while boosting antitumor immune responses through elevated IFN-γ levels and enhanced infiltration of CD8+ tumor-infiltrating lymphocytes (). Furthermore, LJZD has been reported to improve anorexia associated with tumors or chemotherapy (Wu et al., 2022; Sun et al., 2020), suggesting its potential as an adjuvant in cancer therapy. However, whether LJZD has similar antitumor effects on glioma remains unclear.
To further elucidate the potential role of LJZD in glioma treatment, we employed a network pharmacology approach to analyze the main active ingredients of LJZD and their targets. By intersecting these targets with glioma-related targets, we constructed a target network elucidating LJZD’s anti-glioma mechanism. Core targets identified from this network were subjected to molecular docking with LJZD’s active ingredients to evaluate binding stability based on binding energy. Finally, we evaluated the effects of LJZD’s key active ingredients on glioma cell biological behaviors (cell viability, clonogenicity and apoptosis) as well as the expression patterns of core targets and their downstream molecules in vitro.
2 Materials and methods
2.1 Screening active ingredients and targets of LJZD
Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP, version 2.3) (Ru et al., 2014) serves as a specialized systems pharmacology platform for Chinese herbal medicines, illuminating the complex relationships between drugs, their targets, and associated diseases, and it has become a cornerstone in TCM research. The classic LJZD formula consists of six herbs: Panax Ginseng C. A. Mey., Atractylodes Macrocephala Koidz., Poria Cocos (Schw.) Wolf., licorice, Arum Ternatum Thunb., and Citrus Reticulata. We searched for the active ingredients of these six herbs from TCMSP. The SMILES of the active ingredients were input into the SwissADME website () for ADME pharmacokinetic analysis. Initially, we selected ingredients with “High” oral bioavailability and “Yes” blood-brain barrier permeability. Subsequently, we implemented a modified Lipinski-based screening strategy: while retaining the classic thresholds for molecular weight (≤500 Da), lipophilicity (XLOGP ≤5), hydrogen bond donors (≤5), and acceptors (≤10), we introduced an additional molecular weight cutoff (>200 Da) to exclude excessively small molecules with potentially insufficient target engagement. Additionally, ingredients that met two or more of the Extended Drug-Likeness Rules (Ghose, Veber, Egan, and Muegge) were defined as the active ingredients of LJZD. The targets of the active ingredients of LJZD were predicted using SwissTargetPrediction (), Herb (version 2.0) (), and BATMAN-TCM (version 2.0) () platforms. We incorporated all targets predicted by Herb, retained targets with probability >0.6 from SwissTargetPrediction, and selected both known targets and predicted targets with score >80 from BATMAN. The targets obtained from the three websites were combined and deduplicated to obtain the potential targets for LJZD.
2.2 Glioma-related targets screening
We searched for glioma-related targets in public databases using “glioma” as the keyword. The databases and screening criteria are as follows: In the GeneCards database (Stelzer et al., 2016), we included targets with a relevance score >5. In the Online Mendelian Inheritance in Man (OMIM) database (), all targets were included. In the DisGeNET database (Pinero et al., 2020), we included targets with scoreGDA >0.6. In the Open Targets Platform (), we included targets with globalscore >0.6. In the Comparative Toxicogenomics Database (), we included targets with inference score >40. All targets were combined and deduplicated to form the set of glioma-related targets.
2.3 Construction of protein-protein interaction (PPI) network and bioinformatics analysis
To identify the potential targets of LJZD against glioma, we intersected the targets of LJZD with the glioma-related targets and constructed a Venn diagram (). The intersecting targets were input into the STRING database (Version 12.0) (www.string-db.org) for PPI analysis, with the species set as “Homo sapiens” and a medium confidence score of 0.4. The PPI network was visualized using Cytoscape 3.10.3. We analyzed the network using Centiscape 2.2 to obtain the betweenness, closeness, and degree values for each node, which reflect the importance of the nodes in the network. We calculated the average betweenness, closeness, and degree values for all intersecting targets and selected targets with values above the average as hub targets. The interaction network of hub targets was constructed using Cytoscape 3.10.3.
To further understand the functions of the targets related to LJZD’s treatment of glioma and the affected signaling pathways, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses on the intersecting targets using the DAVID bioinformatics analysis platform (Sherman et al., 2022; ). The GO analysis covered three categories: biological process (BP), cellular component (CC), and molecular function (MF). For the KEGG pathway enrichment analysis, we used an FDR <0.05 as the cut-off value. The results of the GO and KEGG pathway enrichment analyses were visualized using the bioinformatics platform (www.bioinformatics.com.cn).
2.4 Molecular docking
To further elucidate the interactions between the active ingredients of LJZD and core targets, we constructed a network of intersection targets and active ingredients using Cytoscape 3.10.3. We selected the top 10 active ingredients with the highest degree as the main active ingredients for glioma treatment. Similarly, we identified the top 10 targets with the highest degree from the PPI network analysis as core targets. Subsequently, we performed molecular docking between the top 10 active ingredients and core targets to assess their binding affinities. The 3D structure files of the active ingredients were obtained from PubChem (). The 3D crystal structures of the top 10 core target proteins were retrieved from the RCSB Protein Data Bank (rcsb.org) (). Molecular docking was conducted using AutoDock Vina (Version 1.2.0) (). We recorded the binding energies of the active ingredients and target proteins and visualized these interactions using a binding energy heatmap. Lower binding energies indicate higher binding affinities. Generally, binding energies below −5.0 kcal/mol are considered good binding, while those below −7.0 kcal/mol are considered strong binding. The PDB files of the optimal binding conformations of the active ingredients and core targets were submitted to the Protein-Ligand Interaction Profiler website for intermolecular interaction analysis (). The results were visualized using PyMOL software (Version 3.11).
2.5 Experimental validation of anti-glioma effect of LJZD
2.5.1 Chemicals and reagents
The active ingredients of LJZD, including pinocembrin (HY-N0575), glabridin (HY-N0393), and chrysophanol (HY-13595), were purchased from MedChemExpress (Monmouth Junction, NJ, USA). CCK-8 kit was obtained from Sangon Biotech (Shanghai, China; Cat. No. E606335-0500), while the apoptosis detection kit (Cat. No. A211-01) was sourced from Vazyme (Jiangsu, China). The antibodies used in the Western blot were purchased from HUABIO (Hangzhou, China), with the following details: AKT (ET1609-51, 1:1,000), phosphorylated AKT (p-AKT, ET1607-73, 1:1,000), phosphorylated STAT3 (p-STAT3, ET1603-40, 1:1,000), phosphorylated c-Jun (p-c-Jun, HA722475, 1:1,000), c-Myc (HA722895, 1:1,000), HIF-1α (HA721997, 1:1,000), Cleaved caspase 3 (HA722367, 1:1,000). Antibody for GAPDH was purchased from Cell Signaling Technology (#5174, CST, 1:1,000).
2.5.2 Cell lines and culture
The human glioma U251 cell line was obtained from Shanghai Fuheng Biology (Shanghai, China). The cells were cultured in high-glucose Dulbecco’s Modified Eagle Medium supplemented with 10% fetal bovine serum and 1% penicillin/streptomycin at 37 °C with 5% CO2. Cell dissociation was performed using 0.25% trypsin (). Log-phase cells were used for subsequent experiments.
2.5.3 Cell viability assay
Logarithmically growing cells were digested, washed and then resuspended in culture medium. After cell counting, the cells were seeded into 96-well plates at a density of 20,000 cells per well for the CCK8 assay. To determine the half-maximal inhibitory concentrations (IC50) of the selected ingredients from LJZD against glioma, U251 cells were treated with gradient concentrations. After 24 h of incubation, 10 μL of CCK8 reagent was added to each well, and the plates were further incubated at 37 °C in the dark for 40 min. Absorbance was measured at 450 nm. The IC50 values of the selected ingredients were calculated using GraphPad Prism software (). Each concentration was tested in quintuplicate.
To evaluate the inhibitory effects of glabridin on the viability of U251 cells, low (40 μM), medium (80 μM), and high (120 μM) concentrations of glabridin were added to the wells based on its IC50 value, with an equivalent volume of DMSO serving as the control. After incubation for 24, 48, and 72 h, the culture medium was replaced, and 10 μL of CCK8 reagent was added to each well. Absorbance was measured at 450 nm following further incubation.
2.5.4 Colony formation assay
U251 cells were seeded into 6-well plates at a density of 1,000 cells per well and cultured for 2 weeks. During this period, cells were treated with low (40 μM), medium (80 μM), and high (120 μM) concentrations of glabridin, with cells treated with DMSO serving as the control group. At the end of the culture period, cell colonies were fixed with paraformaldehyde and stained with crystal violet. The stained colonies were photographed and counted under a microscope to assess the impact of glabridin on cell colony-forming ability.
2.5.5 Apoptosis assay
To assess apoptosis, U251 cells were plated in 6-well plates. The cells were then exposed to varying concentrations of glabridin—low (40 μM), medium (80 μM), and high (120 μM)—with DMSO used as a control. Following 24 h of treatment, the cells were gathered, centrifuged to discard the supernatant, and resuspended in binding buffer. Annexin V-PE-CY7 and propidium iodide were introduced and mixed gently. After a 15-min incubation at room temperature in the dark, the cells were subjected to flow cytometry analysis. Each experiment was conducted in triplicate.
2.5.6 Reverse transcription polymerase chain reaction (RT-PCR)
U251 cells were exposed to 80 μM glabridin. Following a 24-h incubation period, total RNA was isolated using the TRIzol reagent kit (Invitrogen; Carlsbad, CA, USA). The RNA concentration was quantified via spectrophotometry. Subsequently, the RNA was reverse transcribed into cDNA. For the amplification process, specific primers for the target genes (Supplementary Table S1) were incorporated into the cDNA samples. The amplification protocol involved an initial denaturation at 95 °C for 2 min, followed by 40 cycles consisting of 15 s at 95 °C and 30 s at 60 °C ((Ma et al., 2025)). GAPDH was used as an internal control (Tang et al., 2024). The mRNA expression levels were normalized to GAPDH using the 2−ΔΔCt method.
2.5.7 Western blotting
U251 cells were treated with either glabridin (80 μM) or DMSO for 24 h. After washing the cells with PBS, lysis buffer was added, and the cells were incubated before centrifugation to obtain total protein. Protein concentration was measured using a BCA protein assay kit. The protein samples were mixed with 5× SDS-PAGE loading buffer and boiled at 100 °C for 5 min to denature the proteins. The proteins were then separated on a 10% SDS-PAGE gel and transferred to a PVDF membrane. The membrane was blocked with 5% non-fat milk and incubated overnight at 4 °C with primary antibodies, using GAPDH as an internal control. After washing to remove excess primary antibodies, the membrane was incubated with secondary antibodies at 37 °C. Chemiluminescence was detected using an ECL kit. The gray values of the Western blot bands were analyzed using ImageJ software.
2.6 Statistical analysis
All data are presented as the mean ± standard deviation. Statistical analyses were conducted using GraphPad Prism 10 (GraphPad Software, Boston, MA, USA). Comparisons between two groups were performed using the t-test, while comparisons among multiple groups were analyzed using one-way analysis of variance (ANOVA), followed by Sidak’s test for post-hoc pairwise comparisons. A P-value of less than 0.05 was considered statistically significant.
3 Results
3.1 Targets associated with LJZD’s anti-glioma activity
Using the TCMSP database, we identified 190 active ingredients in Panax Ginseng C. A. Mey., 55 in Atractylodes Macrocephala Koidz., 34 in Poria Cocos (Schw.) Wolf., 63 in Citrus Reticulata., 280 in licorice, and 116 in Arum Ternatum Thunb. We then analyzed the pharmacokinetic ADME properties of each active ingredient using SwissADME and selected 90 ingredients that met our criteria. Subsequently, we predicted the targets of these active ingredients using SwissTargetPrediction, Herb, and BATMAN, excluding 14 ingredients with no predicted targets. In total, 76 active ingredients and 752 targets were identified for LJZD (Figure 1; Supplementary Table S2).
FIGURE 1
To identify glioma-related targets, we searched public databases including GeneCards, OMIM, DisGeNET, Open Targets, and CTD. Applying the screening criteria outlined in Section 2.2, we retrieved 1,022 targets from GeneCards, 196 from OMIM, 138 from DisGeNET, 41 from Open Targets, and 431 from CTD. After merging, we identified 1,456 glioma-related targets. Finally, we intersected the LJZD targets with the glioma-related targets, resulting in 174 intersection targets (Figure 2). To elucidate the relationships among these intersection targets, we constructed a PPI network using the STRING database (Figure 3A). This network comprises 172 nodes and 4932 edges, with mean betweenness, closeness, and degree values of 120.105, 0.0035, and 57.349, respectively. Among the 174 targets, 40 targets had betweenness, closeness, and degree values above the mean (Figure 3B) and were designated as hub targets. TP53 exhibited the highest betweenness, closeness, and degree, followed by AKT1, HIF1A, and TNF. Additionally, JUN, IL6, STAT3, and MYC ranked in the top 10, highlighting their significance in the hub network. The PPI network of hub targets is shown in Figure 3C.
FIGURE 2
FIGURE 3
3.2 GO and KEGG analysis of intersection targets
To further reveal the intracellular biological functions of intersection targets and their involvement in signaling pathways during LJZD treatment of glioma, we performed GO and KEGG analyses on the 174 intersection targets using the DAVID platform. The GO analysis yielded a total of 2556 significant terms (FDR <0.05), comprising 2164 BP terms, 243 MF terms, and 149 CC terms. Using the bioinformatics platform, we visualized the top 10 terms in each category (Figures 4A,B; Supplementary Table S3). As shown in Figure 4, LJZD’s effects on glioma primarily involve several intracellular biological processes, including enzyme binding, identical protein binding, transcription factor binding, ubiquitin protein ligase binding, protein domain-specific binding, and protein kinase binding. These findings highlight the potential of LJZD to regulate protein function and signal transduction mechanisms. The main intracellular ingredients implicated include the cytoplasm, cytosol, protein-containing complex, organelle lumen, and membrane-enclosed lumen. Molecular functions primarily involve cellular responses to stimuli (e.g., cellular response to chemical stimulus, response to chemical, response to oxygen-containing compound) and regulation of cell death (e.g., regulation of programmed cell death and regulation of apoptotic process).
FIGURE 4
KEGG pathway enrichment analysis identified a total of 170 pathways. As shown in Figure 5A, the top 50 pathways with the highest fold enrichment scores (Supplementary Table S4) include apoptosis and the p53 signaling pathway in cellular processes, suggesting that LJZD may induce apoptosis in glioma cells. Pathways in environmental information processing involve VEGF, HIF-1, ErbB, TNF, and FoxO signaling pathways. Additionally, several pathways are associated with human diseases, including various common cancers such as non-small cell lung cancer, colorectal cancer, and acute myeloid leukemia. To more clearly illustrate the specific molecular pathways affected by LJZD in glioma, we removed pathways related to other tumors or specific diseases and presented the top 20 pathways in a bubble chart (Figure 5B). Among the top 20 enriched pathways, those related to current cancer drug treatments, such as platinum drug resistance, EGFR tyrosine kinase inhibitor resistance, PD-L1 expression and PD-1 checkpoint pathway in cancer, and antifolate resistance, have high enrichment scores, indicating the potential of LJZD as a therapeutic agent for glioma. Based on these findings, we constructed a LJZD ingredient-target-pathway network (Figure 5C).
FIGURE 5
3.3 Molecular docking of core targets and major active ingredients of LJZD
To identify the major active ingredients of LJZD that contribute to its anti-glioma effects, we constructed an intersection target-ingredient-formula network (Figure 6A). Network analysis revealed that the top five active ingredients (ranked by their degree) are isoliquiritigenin (degree = 44), pinocembrin (degree = 39), lauric acid (degree = 36), formononetin (degree = 32), and shogaol (degree = 28). In the interaction network, a higher degree indicates more connections with other nodes. In the intersection target-ingredient-formula network, the degree of active ingredients indicates their influence on the intersection targets. Based on this, we selected the top 10 active ingredients with the highest degree as the major active ingredients of LJZD for glioma treatment. Similarly, we identified the top 10 targets with the highest degree from the PPI network of intersection targets as core targets. To further elucidate the interactions between these major active ingredients and core targets (Tables 1, 2), we performed molecular docking analysis to assess their binding affinities and binding modes.
FIGURE 6
TABLE 1
| Rank | Ingredient code | Name | PubChem CID | SMILES | Degree |
|---|---|---|---|---|---|
| 1 | GC8 | Isoliquiritigenin | 638278 | C1 = CC(=CC = C1/C=C/C (=O)C2 = C(C=C(C=C2)O)O)O | 44 |
| 2 | GC7 | Pinocembrin | 68071 | C1 [C@H](OC2 = CC(=CC(=C2C1 = O)O)O)C3 = CC = CC = C3 | 39 |
| 3 | D | Lauric acid | 3893 | CCCCCCCCCCCC(=O)O | 36 |
| 4 | GC10 | Formononetin | 5280378 | COC1 = CC = C(C=C1)C2 = COC3 = C(C2 = O)C=CC(=C3)O | 32 |
| 5 | BX3 | Shogaol | 5281794 | CCCCC/C=C/C (=O)CCC1 = CC(=C(C=C1)O)OC | 28 |
| 6 | ChP2 | Tangeretin | 68077 | COC1 = CC = C(C=C1)C2 = CC(=O)C3 = C(O2)C (=C(C(=C3OC)OC)OC)OC | 27 |
| 7 | GC27 | Glabridin | 124052 | CC1(C=CC2 = C(O1)C=CC3 = C2OC [C@H](C3)C4 = C(C=C(C=C4)O)O)C | 20 |
| 8 | BX2 | Chrysophanol | 10208 | CC1 = CC2 = C(C(=C1)O)C (=O)C3 = C(C2 = O)C=CC = C3O | 19 |
| 9 | B | Dibutyl Phthalate | 3026 | CCCCOC(=O)C1 = CC = CC = C1C(=O)OCCCC | 19 |
| 10 | GC9 | 7-Methoxy-2-methyl isoflavone | 354368 | CC1 = C(C(=O)C2 = C(O1)C=C(C=C2)OC)C3 = CC = CC = C3 | 16 |
The top 10 main active ingredients of LJZD.
TABLE 2
| Rank | Core protein | PDB ID | Degree |
|---|---|---|---|
| 1 | TP53 | 5O1C | 151 |
| 2 | AKT1 | 3O96 | 144 |
| 3 | HIF1A | 4H6J | 136 |
| 4 | TNF | 6OOY | 134 |
| 5 | JUN | 5FV8 | 131 |
| 6 | IL6 | 1P9M | 130 |
| 7 | CASP3 | 1NM | 128 |
| 8 | BCL2 | 8HOG | 128 |
| 9 | STAT3 | 6NJS | 125 |
| 10 | MYC | 1NKP | 124 |
The top 10 core targets between LJZD and glioma.
As shown in Figure 6B, except for lauric acid and dibutyl phthalate, the remaining eight active ingredients exhibited strong binding affinities (binding energy < −5 kcal/mol) with all 10 core proteins. For core targets such as TP53 (degree = 151), AKT1 (degree = 144), and TNF (degree = 134), most active ingredients demonstrated very strong binding (binding energy < −7 kcal/mol), highlighting their significant roles in LJZD’s anti-glioma effects. Among all 10 active ingredients, glabridin showed the highest binding affinity with eight out of 10 core targets (TP53, AKT1, HIF1A, JUN, CASP3, BCL2, STAT3, MYC). Chrysophanol exhibited the lowest binding energy with TNF (BE = −10.5 kcal/mol), while pinocembrin had the strongest binding affinity with IL6 (BE = −7.4 kcal/mol). We visualized the binding conformations and interactions of each core target-active ingredient pair with the lowest binding energy (Figure 7).
FIGURE 7
As illustrated in Figure 7, the hydroxyl group at the 4′-position of glabridin participates in the formation of hydrogen bonds within the binding conformations of multiple core proteins. This hydroxyl group acts as a hydrogen bond donor, forming hydrogen bonds with GLN320 in HIF-1α, GLN99 in BCL2, and ASP228 in TP53. It also acts as a hydrogen bond acceptor, forming hydrogen bonds with GLU612 and SER613 in STAT3, VAL147 in TP53, LYS268 and SER205 in AKT1, and ARG739 in MYC. In addition to these interactions, glabridin forms parallel π-stacking interactions with TYR340 in HIF-1α, TYR202 in BCL2, and TRP80 in AKT1. It also engages in π-cation interactions with LYS539 and ARG739 in MYC. In the binding conformation of pinocembrin with IL6, pinocembrin forms four hydrogen bonds with LYS46, ASP160, and THR43 of IL6. Additionally, it forms hydrophobic interactions with PHE105 and LEU39 of IL6, which together constitute the structural basis for the stable binding of pinocembrin to IL6 (binding energy = −7.4 kcal/mol). In the binding conformation of TNF with chrysophanol, chrysophanol is located within the channel formed by the TNF trimer. It forms hydrogen bonds with TYR119 in chain B, TYR151 in chain A, and TYR59 in chain A of TNF. It also forms hydrophobic interactions with multiple tyrosine residues in chains C and A of TNF, ensuring a stable binding interaction.
3.4 Active ingredients of LJZD inhibited U251 glioma cell viability, colony formation and trigger apoptosis
As demonstrated by molecular docking analysis, glabridin, chrysophanol, and pinocembrin - exhibited strong binding affinity with the top 10 core targets (Supplementary Figure S1). Consequently, we selected these three compounds for further validation of their anti-glioma effects. Figure 8A shows that all three components significantly inhibited the viability of U251 glioma cells, with IC50 values of 80.35 μM for glabridin, 686.2 μM for chrysophanol, and 292.5 μM for pinocembrin. Given that glabridin not only demonstrated the lowest IC50 value but also showed the highest affinity with 8/10 of LJZD’s top anti-glioma core targets, we selected glabridin for subsequent experiments. Our results revealed that glabridin inhibited U251 cell viability in both concentration- and time-dependent manners (Figure 8B). Moreover, glabridin treatment markedly suppressed the colony-forming ability of U251 cells, with the inhibitory effect strengthening as concentrations increased (Figure 8C). Flow cytometry analysis further confirmed that glabridin effectively promoted apoptosis in U251 cells (Figure 8D). It is noteworthy that among LJZD’s core anti-glioma targets, BCL2 and CASP3 are key regulators of apoptosis. Elevated BCL2 expression enhances cellular resistance to apoptosis, while CASP3-encoded caspase-3 serves as the key executioner protease in the apoptotic cascade. Our findings demonstrated that treatment with 80 μM glabridin significantly downregulated BCL2 mRNA levels (Figure 8E) and markedly increased cleaved caspase-3 protein expression (Figure 8F) in U251 cells. These results collectively demonstrate that glabridin exerts anti-glioma effects by inducing apoptosis in glioma cells.
FIGURE 8
3.5 Active ingredients of LJZD modulates core targets related to drug resistance in glioma
To elucidate the molecular mechanisms by which glabridin modulates LJZD’s core anti-glioma targets, we systematically analyzed the expression changes of key signaling molecules in glabridin-treated U251 cells, including total AKT, p-AKT (Ser473), p-STAT3 (Tyr705), p-c-Jun (Ser73), as well as transcriptional regulators c-Myc and HIF-1α. Western blot analysis revealed (Figure 9A) that glabridin treatment significantly inhibited the phosphorylation activation of AKT, STAT3, and c-Jun, while reducing the protein expression levels of c-Myc and HIF-1α in U251 cells. KEGG pathway enrichment analysis demonstrated that LJZD’s core targets were significantly enriched in various signaling pathways related to drug resistance. Moreover, among these core targets, Akt, STAT3, and c-Jun have all been shown to transcriptionally regulate the expression of drug efflux pump ABCB1. Our results showed that glabridin treatment markedly downregulated ABCB1 mRNA expression in U251 cells (Figure 9B). Furthermore, based on the established roles of c-Myc and HIF-1α in regulating key glycolytic enzymes and promoting tumor drug resistance (), we assessed the mRNA expression levels of HK2 and LDHA. The data demonstrated that glabridin treatment significantly suppressed the transcriptional activity of these two critical glycolytic enzymes (Figure 9C). These findings collectively demonstrate that glabridin, as one of the principal active ingredients of LJZD, can influence glioma drug resistance mechanisms through a multi-target regulatory network.
FIGURE 9
4 Discussion
Although there have been no previous reports on the application of LJZD in glioma treatment, it has demonstrated promising therapeutic effects in esophageal cancer () and liver cancer (). In this study, we identified 76 active ingredients from LJZD and obtained 174 intersection targets between LJZD and glioma. These intersection targets were enriched in 2164 BP terms, 243 MF terms, 149 CC terms, and 170 pathways. These results strongly suggest that LJZD may exert its anti-glioma effects through multiple active ingredients acting on multiple targets and pathways, which aligns with the multi-ingredient, multi-target, and multi-pathway characteristics of TCM(()). The malignant proliferation, invasion, and treatment resistance of glioma are driven by a complex interplay of signaling pathways, including the p53, receptor tyrosine kinase, and Wnt signaling pathways, as well as cellular processes such as apoptosis, autophagy, senescence, and necrosis (Sottoriva et al., 2013; Patel et al., 2014). Given the intricate molecular network underlying these processes, multi-target and multi-pathway regulation is likely more effective than single-target drug therapy in controlling tumors, which is one of the advantages of TCM formulas. These formulas, composed of various herbs, can simultaneously modulate the same target or different nodes of the same pathway, thereby enhancing therapeutic efficacy through synergistic effects. As shown in Figure 5A, 71 out of 76 active ingredients in LJZD collectively influence 174 common targets, with targets such as PTGS2, ESR1, and NOS2 being regulated by over 30 active ingredients, highlighting the networked regulatory effects of LJZD on glioma-related targets and pathways.
PPI network analysis identified 40 key targets of LJZD in glioma, with TP53, AKT1, and HIF1A occupying central positions in the network topology. These targets exhibit both high degree and high betweenness, indicating that they serve as hubs for multiple pathways and critical control points for information flow within the network. They are also key molecules in the development, progression, and treatment resistance of glioma. It has been well established that the p53 pathway and the PI3K/Akt/mTOR pathway are central to glioblastoma tumorigenesis, pathophysiology, and resistance to treatment. Studies have shown that dysregulation of the p53-ARF-MDM2 pathway and PI3K/Akt pathway were observed in most glioma cases (Obrador et al., 2024; ). While disrupting p53-MDM2 interaction () or inhibiting PI3K/Akt pathway () effectively suppressed glioma progression in vitro and in vivo. HIF-1α, encoding by HIF1A, has been elucidated that its expression level is positively correlated with the pathological grade and microvascular density of glioma and negatively correlated with overall patient survival (). Targeting HIF-1α is also one of the promising strategies for glioma treatment (Peng et al., 2021). Overall, given the pivotal roles of these targets in glioma, LJZD holds promise as a potential therapeutic strategy.
Our KEGG enrichment results revealed that apoptosis was the pathway with the highest fold enrichment score, while the GO analysis showed both “regulation of apoptotic process” and “regulation of programmed cell death” were among the top 10 BP terms, suggesting that inducing apoptosis is likely one of the mechanisms by which LJZD exerts its anti-glioma effects. In vitro experiments further demonstrated that all the three active ingredients (glabridin, chrysophanol, and pinocembrin) from LJZD significantly inhibited the viability of U251 glioma cells. While chrysophanol has previously been reported to suppress proliferation and induce apoptosis in U251 and SHG-44 cells (), the effects of glabridin on glioma were unclear until the present study demonstrated that glabridin triggers apoptosis in U251 cells in a concentration-dependent manner, substantiating the pro-apoptotic activity of LJZD. Additionally, key targets involved in the core mechanisms of LJZD, such as TP53, BCL2, and CASP3, are closely related to apoptosis. Subsequent molecular docking further illuminated potential multi-target pro-apoptotic mechanisms. Small molecules that occupy the TP53-Y220C pocket stabilize the mutant protein and restore transcriptional activity (), and six of the top ten LJZD constituents displayed binding energies below −7 kcal/mol to this pocket, with glabridin achieving the lowest at −8.2 kcal/mol, suggesting its potential to rescue mutant p53 function. The P2 hydrophobic pocket of BCL2 is the canonical inhibitory binding site (Liu et al., 2024), and both glabridin and chrysophanol docked here with favourable energies of −7.8 kcal/mol, implying neutralization of BCL2’s anti-apoptotic activity. Further validation revealed that glabridin exposure elevated cleaved-CASP3 and concurrently reduced BCL2 mRNA levels, consistent with previous reports of glabridin-mediated BCL2 repression in gastric (Liu et al., 2024) and breast cancers (), indicating a possible dual inhibition of BCL2 and concomitant apoptosis induction. Together, these findings position LJZD as a multi-ingredient, multi-target regimen that induces glioma apoptosis through stabilisation of mutant TP53, inhibition of BCL2 and activation of CASP3. Altogether, the convergence of network pharmacology, molecular docking and experimental validation frames LJZD as a multi-target pro-apoptotic regimen and underscores its translational promise for glioma therapy.
Among the top 10 KEGG pathways, several are associated with tumor pharmacotherapy, including platinum drug resistance and EGFR tyrosine kinase inhibitor resistance. Previous studies identified ABCB1-mediated efflux as a common mechanism underlying resistance to both platinum drugs and EGFR-TKIs (; Zhang et al., 2025; ). ABCB1 transcription is governed by PI3K/AKT, STAT3, HIF-1α and JNK pahways (Salaroglio et al., 2022; Schulz et al., 2023). We found that glabridin simultaneously reduced the protein levels of p-AKT, p-STAT3, HIF-1α and p-c-Jun in U251 cells, leading to a decrease in ABCB1 mRNA (Figure 9). These data suggest that glabridin may downregulates ABCB1 through convergent inhibition of multiple drug resistance-associated signaling axes. In parallel, glabridin suppressed the glycolytic enzymes HK2 and LDHA—two established HIF-1α targets that drive the Warburg effect. Consistent with reports that inhibition of Warburg effect sensitizes tumors to platinum agents and EGFR-TKIs (Tyagi et al., 2021; ), our findings indicate that glabridin may augment drug efficacy via the HIF-1α–HK2/LDHA axis.
Temozolomide (TMZ) remains the first-line chemotherapeutic for glioblastoma, yet intrinsic or acquired resistance limits its efficacy. Although the present study did not directly examine TMZ, the glabridin-perturbed network we delineated aligns with a growing body of evidence linking each of these targets to TMZ sensitization. Studies show that blockade of LDHA ((Yue et al., 2024)) or HK2((Yang et al., 2024)) restores TMZ sensitivity in glioma cells, mirroring the reduced LDHA and HK2 expression we observed after glabridin treatment. Similarly, silencing HIF-1α increases glioma TMZ chemosensitivity (Wang et al., 2025), and STAT3 inhibition downregulates MGMT to overcome TMZ resistance (), both of which are consistent with glabridin-mediated suppression of HIF-1α and p-STAT3. Furthermore, AKT/mTOR pathway inhibition enhances TMZ efficacy (Yin et al., 2024), paralleling the decreased p-AKT levels we report. Additionally, ABCB1 and related efflux transporters are established drivers of TMZ resistance, and their down-modulation restores TMZ sensitivity (Phon et al., 2025); reduction of ABCB1 expression by glabridin in U251 cells therefore positions it as a candidate for reversing this resistance mechanism. Collectively, these convergent literature findings suggest that the multi-pathway suppression effected by glabridin could meaningfully enhance TMZ responsiveness in glioma.
Finally, additional constituents of LJZD, such as isoliquiritigenin (Wang et al., 2019), shagaol (), and tangeretin (Xie et al., 2022), have independently been shown to restore TMZ, cisplatin or EGFR-TKIs sensitivity. Collectively, these observations position LJZD as a promising multi-ingredient regimen for overcoming chemo- and targeted-therapy resistance in glioma.
5 Conclusion
In conclusion, integrating network pharmacology with molecular docking and cellular assays, we delineated the multi-component, multi-target and multi-pathway mechanisms through which LJZD may counteract glioma. The study provided a cell-level findings for its use as an adjuvant sensitizer in combination with chemotherapy or targeted therapy, which should be validated fin animal models and clinical settings. Furthermore, in-depth exploration and necessary structural optimization of the core active ingredients screened in this study (such as pinocembrin, glabridin, or chrysophanol) may lead to the development of novel small-molecule drugs, thereby advancing clinical treatment progress for glioma.
6 Limitation
First, we applied Lipinski’s Rule of Five in combination with oral bioavailability and blood-brain barrier permeability to screen active ingredients in LJZD, aiming to prioritize candidate molecules with potential for oral drug development. However, this approach may exclude certain pharmacologically active components, such as polysaccharides and saponins, due to their molecular weight or polarity exceeding the predefined criteria. Although some active ingredients demonstrate suboptimal oral absorption or blood-brain barrier penetration, their therapeutic efficacy could potentially be enhanced through functional group modifications or alternative administration routes (e.g., intravenous injection or localized delivery). Therefore, the therapeutic potential of other active ingredients excluded by our screening criteria should not be disregarded. Second, to elucidate LJZD’s effects on core targets, we selected glabridin for experimental validation based on molecular docking results, which yielded promising results. It should be noted that glabridin was identified from 76 active ingredients in LJZD, and other ingredients may also possess considerable anti-glioma activity warranting further investigation. Meanwhile, further identification of LJZD active ingredients by mass-spectrometric analysis may yield additional novel findings. Third, the mechanistic validation reported here was conducted primarily in U251 cells; future work should extend these observations to additional glioma cell lines, primary cultures, and in vivo models. Moreover, the synergistic interactions among multiple ingredients represent a crucial aspect of TCM’s multi-ingredient, multi-target mechanism. Building upon single-component studies, exploring the network-based regulatory effects of multiple ingredients on key targets and pathways will provide a more comprehensive understanding of LJZD’s anti-glioma mechanisms and advance glioma treatment development.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
Ethical approval was not required for the studies on humans in accordance with the local legislation and institutional requirements because only commercially available established cell lines were used.
Author contributions
JH: Methodology, Writing – original draft, Investigation, Formal Analysis, Software. SL: Writing – original draft, Visualization, Software, Funding acquisition, Methodology. LL: Formal Analysis, Writing – original draft, Funding acquisition. LX: Conceptualization, Writing – review and editing, Resources, Project administration. DK: Supervision, Writing – review and editing, Conceptualization, Project administration, Resources.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Natural Science Foundation of Fujian Province 2021J011393 and 2023J011744.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2025.1682550/full#supplementary-material
References
1
AdasmeM. F.LinnemannK. L.BolzS. N.KaiserF.SalentinS.HauptV. J.et al (2021). PLIP 2021: expanding the scope of the protein-ligand interaction profiler to DNA and RNA. Nucleic Acids Res.49 (1), 530–534. 10.1093/nar/gkab294
2
AmbergerJ. S.BocchiniC. A.SchiettecatteF.ScottA. F.HamoshA. (2015). OMIM.org: online Mendelian Inheritance in Man (OMIM®), an online catalog of human genes and genetic disorders. Nucleic Acids Res.43, D789–D798. 10.1093/nar/gku1205
3
BardouP.MarietteJ.EscudieF.DjemielC.KloppC. (2014). Jvenn: an interactive Venn diagram viewer. BMC Bioinforma.15 (1), 293. 10.1186/1471-2105-15-293
4
BaudM. G. J.BauerM. R.VerduciL.DinglerF. A.PatelK. J.Horil RoyD.et al (2018). Aminobenzothiazole derivatives stabilize the thermolabile P53 cancer mutant Y220c and show anticancer activity in P53-Y220c cell lines. Eur. J. Med. Chem.152, 101–114. 10.1016/j.ejmech.2018.04.035
5
BermanH. M.WestbrookJ.FengZ.GillilandG.BhatT. N.WeissigH.et al (2000). The protein Data Bank. Nucleic Acids Res.28 (1), 235–242. 10.1093/nar/28.1.235
6
BunielloA.SuvegesD.Cruz-CastilloC.LlinaresM. B.CornuH.LopezI.et al (2025). Open targets platform: facilitating therapeutic hypotheses building in drug discovery. Nucleic Acids Res.53 (1), D1467–D1475. 10.1093/nar/gkae1128
7
Cancer Genome Atlas ResearchN. (2008). Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature455 (7216), 1061–1068. 10.1038/nature07385
8
ChaoJ.DaiY.VerpoorteR.LamW.ChengY. C.PaoL. H.et al (2017). Major achievements of evidence-based traditional Chinese medicine in treating major diseases. Biochem. Pharmacol.139, 94–104. 10.1016/j.bcp.2017.06.123
9
CuiX.YaoA.LvJ.ZhangC.ChenY.MeiT.et al (2025). Siramesine induced cell death of glioblastoma through inactivating the stat3-Mgmt signaling pathway. J. Transl. Med.23 (1), 780. 10.1186/s12967-025-06693-y
10
DainaA.MichielinO.ZoeteV. (2017). Swissadme: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry Friendliness of small molecules. Sci. Rep.7, 42717. 10.1038/srep42717
11
DainaA.MichielinO.ZoeteV. (2019). Swisstargetprediction: updated data and new features for efficient prediction of protein targets of small molecules. Nucleic Acids Res.47 (W1), W357-W364–W64. 10.1093/nar/gkz382
12
DavisA. P.WiegersT. C.SciakyD.BarkalowF.StrongM.WyattB.et al (2025). Comparative Toxicogenomics database's 20th anniversary: update 2025. Nucleic Acids Res.53 (1), D1328–D1334. 10.1093/nar/gkae883
13
DuchatelR. J.JacksonE. R.ParackalS. G.KiltschewskijD.FindlayI. J.MannanA.et al (2024). Pi3k/Mtor is a therapeutically targetable genetic dependency in diffuse intrinsic pontine glioma. J. Clin. Invest134 (6), e170329. 10.1172/JCI170329
14
DyrstadS. E.LotsbergM. L.TanT. Z.PettersenI. K. N.HjellbrekkeS.TusubiraD.et al (2021). Blocking aerobic glycolysis by targeting pyruvate dehydrogenase kinase in combination with egfr tki and ionizing radiation increases therapeutic effect in non-small cell lung cancer cells. Cancers (Basel)13 (5), 941. 10.3390/cancers13050941
15
EberhardtJ.Santos-MartinsD.TillackA. F.ForliS. (2021). Autodock Vina 1.2.0: new docking methods, expanded force field, and Python bindings. J. Chem. Inf. Model61 (8), 3891–3898. 10.1021/acs.jcim.1c00203
16
FengJ.LiJ.WuL.YuQ.JiJ.WuJ.et al (2020). Emerging roles and the regulation of aerobic glycolysis in hepatocellular carcinoma. J. Exp. Clin. Cancer Res.39 (1), 126. 10.1186/s13046-020-01629-4
17
GaoK.LiuL.LeiS.LiZ.HuoP.WangZ.et al (2025). Herb 2.0: an updated database integrating clinical and experimental evidence for traditional Chinese medicine. Nucleic Acids Res.53 (D1), D1404–D1414. 10.1093/nar/gkae1037
18
GeW.YuanG.WangD.DongL. (2024). Exploring the therapeutic mechanisms and prognostic targets of biochanin a in glioblastoma via integrated computational analysis and in vitro experiments. Sci. Rep.14 (1), 3783. 10.1038/s41598-024-53442-0
19
German-CortesJ.HerreroR.TorroglosaN.PumarolaA.Fischer-AlbiolN.Campos-MorenoS.et al (2025). Preclinical evaluation of several polymeric micelles identifies Soluplus®-docetaxel as the most effective candidate in multiple glioblastoma models. J. Control Release381, 113616. 10.1016/j.jconrel.2025.113616
20
GuJ.RauniyarS.WangY.ZhanW.YeC.JiS.et al (2021). Chrysophanol induced glioma cells apoptosis via activation of mitochondrial apoptosis pathway. Bioengineered12 (1), 6855–6868. 10.1080/21655979.2021.1972079
21
GuoZ.WangY.QinW.HengY.ChenX.LiuN.et al (2024). Mir-122-3p targets Ube2i to regulate the immunosuppression of liver cancer and the intervention of Liujunzi formula. J. Ethnopharmacol.329, 118081. 10.1016/j.jep.2024.118081
22
HanY.FanX.FanL.WuY.ZhouZ.WangG.et al (2023). Liujunzi decoction exerts potent antitumor activity in oesophageal squamous cell carcinoma by inhibiting mir-34a/stat3/il-6r feedback loop, and modifies antitumor immunity. Phytomedicine111, 154672. 10.1016/j.phymed.2023.154672
23
HaoX.BahiaR. K.CsehO.BozekD. A.BlakeS.RinnenthalJ.et al (2023). Bi-907828, a novel potent Mdm2 inhibitor, inhibits glioblastoma brain tumor stem cells in vitro and prolongs survival in orthotopic xenograft mouse models. Neuro Oncol.25 (5), 913–926. 10.1093/neuonc/noac271
24
HsuC. M.SuH. C.YangM. Y.TsaiY. T.TsaiM. S.YangY. H.et al (2023). 6-Shogaol is a potential treatment for head and neck squamous cell carcinoma. Int. J. Med. Sci.20 (2), 238–246. 10.7150/ijms.80542
25
HuB.ZouT.QinW.ShenX.SuY.LiJ.et al (2022). Inhibition of egfr overcomes acquired lenvatinib resistance driven by stat3-abcb1 signaling in hepatocellular carcinoma. Cancer Res.82 (20), 3845–3857. 10.1158/0008-5472.CAN-21-4140
26
Huang daW.ShermanB. T.LempickiR. A. (2009). Systematic and integrative analysis of large gene lists using david bioinformatics Resources. Nat. Protoc.4 (1), 44–57. 10.1038/nprot.2008.211
27
KimS.ChenJ.ChengT.GindulyteA.HeJ.HeS.et al (2023). Pubchem 2023 update. Nucleic Acids Res.51 (1), D1373–D1380. 10.1093/nar/gkac956
28
KodanA.FutamataR.KimuraY.KiokaN.NakatsuT.KatoH.et al (2021). Abcb1/Mdr1/P-Gp employs an atp-dependent twist-and-Squeeze mechanism to Export hydrophobic drugs. FEBS Lett.595 (6), 707–716. 10.1002/1873-3468.14018
29
KongX.LiuC.ZhangZ.ChengM.MeiZ.LiX.et al (2024). Batman-tcm 2.0: an enhanced integrative database for known and predicted interactions between Traditional Chinese medicine ingredients and target proteins. Nucleic Acids Res.52 (D1), D1110–D1120. 10.1093/nar/gkad926
30
LiuQ.CaoP. (2015). Clinical and prognostic significance of HIF-1α in glioma patients: a meta-analysis. Int. J. Clin. Exp. Med.8 (12), 22073–22083.
31
LiuJ.LiS.WangQ.FengY.XingH.YangX.et al (2024). Sonrotoclax overcomes Bcl2 G101v mutation-induced venetoclax resistance in preclinical models of hematologic malignancy. Blood143 (18), 1825–1836. 10.1182/blood.2023019706
32
MaC.WangL.ZhangR.LiT.LiP.DingY.et al (2025). Linc00601 promotes the progression of glioma via the P-Stat3 signaling pathway. Cancer Cell Int.25 (1), 220. 10.1186/s12935-025-03735-9
33
ObradorE.Moreno-MurcianoP.Oriol-CaballoM.Lopez-BlanchR.PinedaB.Gutierrez-ArroyoJ. L.et al (2024). Glioblastoma therapy: past, present and future. Int. J. Mol. Sci.25 (5), 2529. 10.3390/ijms25052529
34
PatelA. P.TiroshI.TrombettaJ. J.ShalekA. K.GillespieS. M.WakimotoH.et al (2014). Single-cell rna-seq highlights intratumoral heterogeneity in primary glioblastoma. Science344 (6190), 1396–1401. 10.1126/science.1254257
35
PengG.WangY.GeP.BaileyC.ZhangP.ZhangD.et al (2021). The HIF1α-PDGFD-PDGFRα axis controls glioblastoma growth at normoxia/mild-hypoxia and confers sensitivity to targeted therapy by echinomycin. J. Exp. Clin. Cancer Res.40 (1), 278. 10.1186/s13046-021-02082-7
36
PhonB. W. S.ChelliahS. S.OsmanD. E.BhuvanendranS.RadhakrishnanA. K.KamarudinM. N. A. (2025). Revisiting abc transporters and their clinical significance in glioblastoma. Pharm. (Basel)18 (1), 102. 10.3390/ph18010102
37
PineroJ.Ramirez-AnguitaJ. M.Sauch-PitarchJ.RonzanoF.CentenoE.SanzF.et al (2020). The disgenet knowledge platform for disease genomics: 2019 update. Nucleic Acids Res.48 (1), D845-D855–D55. 10.1093/nar/gkz1021
38
RuJ.LiP.WangJ.ZhouW.LiB.HuangC.et al (2014). Tcmsp: a database of systems pharmacology for drug discovery from herbal medicines. J. Cheminform6, 13. 10.1186/1758-2946-6-13
39
SalaroglioI. C.BelisarioD. C.AkmanM.La VecchiaS.GodelM.AnobileD. P.et al (2022). Mitochondrial ros drive resistance to chemotherapy and immune-killing in hypoxic non-small cell lung cancer. J. Exp. Clin. Cancer Res.41 (1), 243. 10.1186/s13046-022-02447-6
40
SchaffL. R.MellinghoffI. K. (2023). Glioblastoma and other primary brain malignancies in adults: a review. JAMA329 (7), 574–587. 10.1001/jama.2023.0023
41
SchulzJ. A.HartzA. M. S.BauerB. (2023). Abcb1 and Abcg2 regulation at the blood-brain barrier: potential new targets to improve brain drug delivery. Pharmacol. Rev.75 (5), 815–853. 10.1124/pharmrev.120.000025
42
ShermanB. T.HaoM.QiuJ.JiaoX.BaselerM. W.LaneH. C.et al (2022). David: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. 50 (1), W216–W221. 10.1093/nar/gkac194
43
SottorivaA.SpiteriI.PiccirilloS. G.TouloumisA.CollinsV. P.MarioniJ. C.et al (2013). Intratumor heterogeneity in human glioblastoma reflects cancer evolutionary dynamics. Proc. Natl. Acad. Sci. U. S. A.110 (10), 4009–4014. 10.1073/pnas.1219747110
44
StelzerG.RosenN.PlaschkesI.ZimmermanS.TwikM.FishilevichS.et al (2016). The genecards suite: from gene data mining to disease Genome sequence analyses. Curr. Protoc. Bioinforma.54 (1), 30 1–1 3. 10.1002/cpbi.5
45
StuppR.MasonW. P.van den BentM. J.WellerM.FisherB.TaphoornM. J.et al (2005). Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. N. Engl. J. Med.352 (10), 987–996. 10.1056/NEJMoa043330
46
StuppR.HegiM. E.MasonW. P.van den BentM. J.TaphoornM. J.JanzerR. C.et al (2009). Effects of radiotherapy with concomitant and adjuvant temozolomide versus radiotherapy alone on survival in glioblastoma in a randomised phase iii study: 5-year analysis of the eortc-ncic trial. Lancet Oncol.10 (5), 459–466. 10.1016/S1470-2045(09)70025-7
47
StuppR.TaillibertS.KannerA.ReadW.SteinbergD.LhermitteB.et al (2017). Effect of tumor-treating fields plus maintenance temozolomide vs maintenance temozolomide alone on survival in patients with glioblastoma: a randomized clinical trial. JAMA318 (23), 2306–2316. 10.1001/jama.2017.18718
48
SunL. L.LaiH. Z.ChenZ. Z.ZhuX. S.LinL. Z. (2020). Modified Liujunzi decoction alleviates chemotherapy-induced anorexia in advanced non-small cell lung cancer: a propensity score matched case-control study. Chin. J. Integr. Med.26 (4), 256–262. 10.1007/s11655-020-3185-5
49
TangH.YangX.LiG.PengK.SunY.JiangL.et al (2024). Development and experimental validation of dephosphorylation-related biomarkers to assess prognosis and immunotherapeutic response in gliomas. Front. Immunol.15, 1488894. 10.3389/fimmu.2024.1488894
50
TyagiK.MandalS.RoyA. (2021). Recent advancements in therapeutic targeting of the Warburg effect in refractory ovarian cancer: a promise towards disease remission. Biochim. Biophys. Acta Rev. Cancer1876 (1), 188563. 10.1016/j.bbcan.2021.188563
51
WangC.ChenY.WangY.LiuX.LiuY.LiY.et al (2019). Inhibition of cox-2, mpges-1 and Cyp4a by isoliquiritigenin blocks the angiogenic Akt signaling in glioma through cerna effect of mir-194-5p and lncrna Neat1. J. Exp. Clin. Cancer Res.38 (1), 371. 10.1186/s13046-019-1361-2
52
WangP.LiaoB.GongS.GuoH.ZhaoL.LiuJ.et al (2025). Temozolomide promotes glioblastoma stemness expression through senescence-associated reprogramming via HIF1α/HIF2α regulation. Cell Death Dis.16 (1), 317. 10.1038/s41419-025-07617-w
53
WuX.DaiY.NieK. (2022). Research progress of Liujunzi decoction in the treatment of tumor-associated anorexia. Drug Des. Devel Ther.16, 1731–1741. 10.2147/DDDT.S365292
54
XieY.FengS. L.HeF.YanP. Y.YaoX. J.FanX. X.et al (2022). Down-regulating nrf2 by tangeretin reverses multiple drug resistance to both chemotherapy and egfr tyrosine kinase inhibitors in lung cancer. Pharmacol. Res.186, 106514. 10.1016/j.phrs.2022.106514
55
YangL.LiS.YuL.LengJ.LiN. (2024). Targeting glycolysis: exploring a new frontier in glioblastoma therapy. Front. Immunol.15, 1522392. 10.3389/fimmu.2024.1522392
56
YinH. T.HuiL.YangJ. H.LiQ.LiM.ZhaoQ. Q.et al (2024). Daurisoline suppress glioma progression by inhibiting autophagy through pi3k/akt/mtor pathway and increases tmz sensitivity. Biochem. Pharmacol.223, 116113. 10.1016/j.bcp.2024.116113
57
YueQ.WangZ.ShenY.LanY.ZhongX.LuoX.et al (2024). Histone H3k9 lactylation confers temozolomide resistance in glioblastoma via luc7l2-mediated Mlh1 intron retention. Adv. Sci. (Weinh)11 (19), e2309290. 10.1002/advs.202309290
58
ZhangY.QiuJ. G.WangW.SunF. L.WangX.LiuW. J.et al (2025). Suppression of cyld by Her3 confers ovarian cancer platinum resistance via inhibiting apoptosis and by inducing drug efflux. Exp. Hematol. Oncol.14 (1), 21. 10.1186/s40164-025-00620-z
Summary
Keywords
glioma, Liujunzi decoction, molecular docking, network pharmacology, drug resistance
Citation
Huang J, Lin S, Lin L, Xu L and Kang D (2025) Multi-pathway mechanisms of Liujunzi decoction in promoting glioma apoptosis and reversing drug resistance: network pharmacology and experimental validation. Front. Pharmacol. 16:1682550. doi: 10.3389/fphar.2025.1682550
Received
09 August 2025
Accepted
10 September 2025
Published
24 September 2025
Volume
16 - 2025
Edited by
Zhenyu Gong, Sun Yat-sen University, China
Updates
Copyright
© 2025 Huang, Lin, Lin, Xu and Kang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Luning Xu, xlning@fjmu.edu.cn; Dezhi Kang, kdz986688@163.com
† These authors have contributed equally to this work and share first authorship
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.