Mechanism underlying the effect of Pulsatilla decoction in hepatocellular carcinoma treatment: a network pharmacology and in vitro analysis
BMC Complementary Medicine and Therapies volume 23, Article number: 405 (2023)
Currently, hepatocellular carcinoma (HCC) is associated with a poor prognosis. Moreover, there exist limited strategies for treating HCC. Pulsatilla decoction (PD), a traditional Chinese medicine formula, has been used to treat inflammatory bowel disease and several cancer types. Accordingly, we explored the mechanism of PD in HCC treatment via network pharmacology and in vitro experiments.
Online databases were searched for gene data, active components, and potential target genes associated with HCC development. Subsequently, bioinformatics analysis was performed using protein–protein interaction and Network Construction and Kyoto Encyclopedia of Genes and Genomes (KEGG) to screen for potential anticancer components and therapeutic targets of PD. Finally, the effect of PD on HCC was further verified by in vitro experiments.
Network pharmacological analysis revealed that 65 compounds and 180 possible target genes were associated with the effect of PD on HCC. These included PI3K, AKT, NF-κB, FOS, and NFKBIA. KEGG analysis demonstrated that PD exerted its effect on HCC mainly via the PI3K-AKT, IL-17, and TNF signaling pathways. Cell viability and cell cycle experiments revealed that PD could significantly inhibit cancer cell proliferation and kill HCC cells by inducing apoptosis. Furthermore, western blotting confirmed that apoptosis was mediated primarily via the PI3K-AKT, IL-17, and TNF signaling pathways.
To the best of our knowledge, this is the first study to elucidate the molecular mechanism and potential targets of PD in the treatment of HCC using network pharmacology.
Hepatocellular carcinoma (HCC) is a highly prevalent cancer globally and is associated with poor prognosis, owing to which it has resulted in a considerable economic burden [1, 2]. Currently, HCC is the second leading cause of death among malignant tumors, with the 5-year survival rate of HCC being extremely low [3, 4]. The primary treatment strategy for early HCC is surgical treatment, and there exist no effective strategies for treating advanced HCC [5, 6]. In recent years, with the advent of targeted therapy and immunotherapy, the prognosis of HCC has improved significantly. For example, immunotherapy has shown encouraging findings in the treatment of HCC . However, immune escape and the emergence of multidrug resistance have limited the application of immunotherapy in HCC. Thus, further research is required to develop new therapeutic agents and targets.
Traditional Chinese medicine (TCM) compounds contain a variety of active antitumor components, owing to which these compounds have been shown to improve the prognosis of patients with several cancer types and have lower side effects than conventional treatments . TCM compounds can inhibit the proliferation, invasion, metastasis, and angiogenesis of tumor cells. In addition, they can enhance the anti-tumor ability of the body and enhance the efficacy of other conventional drugs administered in combination . Pulsatilla decoction (PD) is a TCM compound consisting of Radix Pulsatillae, Cortex Phellodendri, Rhizoma Coptidis, and Cortex Fraxini and is currently used to treat colitis and intestinal tumors [10,11,12]. A few studies have also assessed the application of PD in the treatment of malignant tumors. One such study demonstrated that PD can enhance the antitumor effect of colon cancer by enhancing tumor sensitivity to 5-FU . Another study by Wu et al. revealed that PD can inhibit the proliferation of pancreatic cancer cells as well as induce apoptosis, as assessed using the cell counting kit (CCK-8) and via animal experiments . In addition, PD may exert anticancer effects in colon cancer via multiple components, targets, and pathways, as determined using network pharmacology . Rhizoma Coptidis, one of the main components of PD, has been proven to induce inhibitory and lethal effects on HCC cells [15,16,17]. However, the antitumor mechanism of PD, especially HCC, is still unknown.
In recent years, bioinformatic studies of large online databases provide an ideal platform for researchers to study the active components and therapeutic targets of TCM compounds. This will facilitate the further discovery of new active antitumor components and new therapeutic targets . Network pharmacology is another research technique that is conducive to the development of new therapeutic targets and strategies. It functions by elucidating the relationship between drug, gene, and tumor as well as identifying the relationship between active drug ingredients and target genes . More importantly, network pharmacology can be used to enhance the efficacy of antitumor drugs and reduce the side effects by analyzing the occurrence and development of diseases as well as identifying disease regulatory pathways .
This study aimed to assess the efficacy of PD in the treatment of HCC as well as identify the pathways involved in the antitumor effect of PD, using network pharmacology.
Materials and methods
Composition of PD
PD is composed of Radix Pulsatillae, Cortex Phellodendri, Rhizoma Coptidis, and Cortex Fraxini. The active compounds of PD, with an oral bioavailability (OB) of ≥ 30% and drug-likeness (DL) of ≥ 0.18, were screened from the Chinese Medicine Systematic Pharmacology Database and Analysis Platform (TCMSP: https://old.tcmsp-e.com/tcmsp.php) and Chinese Medicine Molecular Mechanism Biological Information Analysis Tool (http://bionet.ncpsb.org.cn/batman-tcm/).
Predicted targets of PD and HCC-related genes
The putative targets of the four components of PD were collected from TCMSP and DRUGBANK (https://go.drugbank.com/), after that input all putative targets into the uniport database (https://www.uniprot.org/) to transfer corresponding gene symbols of the target proteins. The hepatocellular carcinoma as key term search for two GeneCards (https://www.genecards.org/, ver, 5.11) and OMIM (http://www.omim.org/; updated: July 15, 2022) databases to acquire HCC-related targets. Based on results above two steps, a Venn diagram was drawn to exhibit the intersection targets of PD and HCC. The intersection targets of PD and HCC as hub gene for subsequent study.
Protein–protein interaction (PPI) analysis
In order to reveal relationship between targets and identify core target, the PPI analysis were performed. The overlapping targets of compounds in PD and HCC, as hub genes, was entered into the online STRING database (https://cn.string-db.org/) for PPI analysis. The species was limited to Homo sapiens, and the interaction score was retained as > 0.990, hidden disconnect targets. After that, the TSV files, which containing results of PPI analysis were analyzed using the STRING tool, was download. Subsequently, the TSV files was put into the Cytoscape 3.7.1 to construct and visualize the PPI network. The CytoNCA plugin in the Cytoscape software were used to calculate degree, that a key topological parameter characterizes the most significant nodes in the network. Higher degree quantitative values of targets, more crucial in PD for treating HCC.
Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses
After transferring official gene symbols of hub genes to related Entrez IDs, Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were conducted based on R 3.6.3 and plugins in the R-package such as clusterProfiler, ggplot2, enrichplot, pathview, colorspace, DOSE, org.Hs.eg.db, etc.). Only GO functional and KEGG pathway term with p value < 0.05 were identified statistically significant and stayed for subsequent study.
The network construction is aim to clearly elucidated relationship between nodes. In this study, three network constructions were established and visualized through Cytoscape 3.7.1. (1) bioactive component-disease target network (C-D-T) describe relationship among ingredients, common targets, and disease. The node shape and color represent ingredients, targets, and disease, respectively, while edges stand for interactions between nodes. The network elucidates the interaction between the active components of PD and the targets anti-HCC. (2) PPI network and (3) Bioactive ingredient-pathway-target network. The bioactive ingredient, top 20 signal pathways with the smallest p value of KEGG enrichment analysis, and the corresponding targets of each pathway were submitted into Ctyoscape3.7.1, and drawn bioactive ingredient-pathway-target network diagram of PD against HCC.
Molecular docking is widely used for network pharmacology as it can precisely predict the conformation of small-molecule ligands within the target binding site and assess the binding affinity. For molecular docking, candidate target proteins and compounds were selected based on the results of the PPI network, C-D-T network, and KEGG pathway analyses. Molecular docking was used for validating the compound–target interaction using the Molecular Operating Environment (MOE) (v2015.10) software, as follows. Firstly, the 3D structure of the candidate protein was downloaded from Protein Data Bank (PDB) (http://www.rcsb.org/) with the species limited to “Homo sapiens.” The structure of the compound was obtained from PubChem (https://pubchem.ncbi.nlm.nih.gov/). Subsequently, the 3D structure of the candidate protein was imputed into the MOE to construct a mating pocket after deleting water molecules, preparing the protein structure, and minimizing energy. Lastly, the mating pocket was docked with the candidate compound. Figure 1 shows the experimental methods used in this research.
Preparation of PD extract
Radix Pulsatillae (Bai Tou Weng, Pulsatilla chinensis (Bge.) Regel), Cortex Phellodendri (Huang Bai, Phellodendron amurense Rupr.), Rhizoma Coptidis (Huang Lian, Coptis chinensis Franch.), and Cortex Fraxini (Qin Pi, Fraxinus chinensis Roxb.) were obtained from the Second Xiangya Hospital of Central South University and were used to prepare PD aqueous extract, as per the Chinese Pharmacopoeia (2020 edition) . Briefly, 20 g of Radix Pulsatillae, 16 g of Cortex Phellodendri, 8 g of Rhizoma Coptidis, and 16 g of Cortex Fraxini were soaked with twice the volume of pure water for 30 min and then boiled for 1.5 h. After cooling to room temperature, the mixture was centrifuged at 10,000 rpm for 30 min to collect the supernatant, and the process was repeated twice. The collected supernatant was vacuum dried to obtain a powder. Then, the powder was solubilized with dimethyl sulfoxide (DMSO; 200 mg/ml), filtered, and stored at − 20℃ until further use.
HCC-LM3 and MHCC-97 H cell lines were purchased from the Institute of Biochemistry and Cell Biology, Chinese Academy of Sciences. HCC cell lines were cultured in Dulbecco’s Modified Eagle Medium supplemented with 10% fetal bovine serum and cultured in Petri dishes at 37 °C with 5% CO2 and 95% air.
Cell viability assay
After cell counting, HCC-LM3 and MHCC-97 H cells (logarithmic growth stage) were seeded in 96-well plates at 10,000 cells/well. After incubation for 24 h, the cells were treated with different concentrations of PD for 24, 48, and 72 h. Then, all cells were treated with 0.5% DMSO for 24, 48, and 72 h. Then, the fresh medium was replaced and 10 µl of CCK-8 (Genview, USA) reagent was added to each well. After 2 h incubation, cell viability in each well was determined using a spectrophotometer.
Colony formation assay
After cell counting, HCC-LM3 and MHCC-97 H cells with a cell confluence of 70–80% were seeded in 6-well plates with 3000 cells/well and incubated for 24 h at 37 °C. Based on the experimental results of the cell viability assay. Subsequently, 100 and 200 µg/ml of PD were added to the 6-well plates and incubated for 2 weeks at 37 °C with 5% CO2 and 95% air. Then, the medium was removed, and 4% paraformaldehyde was added and fixed for 30 min. The samples were then stained with crystal violet for 2 h, cleaned, and then photographed for analysis.
Flow cytometry for cell cycle analysis
After cell counting, HCC-LM3 and MHCC-97 H cells in the logarithmic stage of growth (cell confluence: 70–80%) were seeded in 6-well plates at 3 × 105 cells/well and incubated at 37℃ for 24 h. Based on the experimental results of the cell viability assay. Subsequently, 100, 200, and 400 µg/ml of PD were added to the 6-well plates. After 48 h, the cells were harvested and fixed with 70% ethanol overnight. The fixed cells were then stained with propidium iodide (PI, 50 µg/ml; Genview, USA) for 45 min and were analyzed using a Canto II flow cytometer (BD Bioscience, USA).
Flow cytometry analysis of cell apoptosis
HCC-LM3 and MHCC-97 H cells in the logarithmic growth stage were seeded in 6-well plates at 3 × 105 cells/well and incubated at 37℃ for 24 h. Based on the experimental results of the cell viability assay. Subsequently, 200, 400, and 600 µg/ml of PD were added to the 6-well plates for 48 h. Then, cells were digested with trypsin without EDTA. The harvested cells were then resuspended in 200 µl 1× binding buffer and incubated for 15 min with 5 µl Annexin V-FITC and 5 µl PI (Genview, USA). The cells were then analyzed using a Canto II flow cytometer (BD Bioscience, USA). The obtained data were analyzed using the FlowJo ver.7.6 Software (De Novo Software, USA).
After cell counting, HCC-LM3 and MHCC-97 H cells in the logarithmic growth stage (cell confluence: 70–80%) were seeded in 6-well plates at 3 × 105 cells/well and incubated for 37℃ for 24 h. PD (100, 200, and 400 µg /ml) was then added to the plates and treated for 48 h. The protein was denatured, electrophoresed, and transferred to the PVDF membrane. The PVDF membrane was then blocked with 5% BSA for 2 h at room temperature. And the blots were cut prior to hybridisation with antibodies during blotting. Subsequently, antibodies were added and incubated overnight at 4 °C. The primary antibodies used were those against PI3K (1:1000, ABclonal, Cat No: A0265), AKT (1:1000, ABclonal, Cat No: A17909), NF-κB (1:1000, ABclonal, Cat No: A6667), FOS (1:1000, Proteintech, Cat No: 66590-1-Ig), Caspase 3 (1:1000, Proteintech, Cat No: 19677-1-AP), NFKBIA (1:1000, ABclonal, Cat No: A11397), ERK1/2 (1:1000, Proteintech, Cat No: 16443-1-AP), MMP3 (1:1000, Proteintech, Cat No: 17873-1-AP), and GAPDH (1:10000, Proteintech, Cat No: 60004-1-Ig). The membrane was then incubated with the secondary antibody (1:5000, ABclonal) for 1 h. The enhanced chemiluminescence detection kit (Genview, USA) was used to detect the fluorescence of the bands, which were quantitatively analyzed using the ImageJ software (Version 11).
This project was approved by the Medical Ethics Committee of the Second Xiangya Hospital of Central South University.
The student’s t-test and one-way ANOVA analysis were used for statistically analyzing all data, and the GraphPad Prism 6 software (La Jolla, CA, USA) was used for data processing. All data were expressed as means ± standard deviations, and results with P < 0.05 were considered statistically significant.
Compounds and putative targets of PD
PD consists of Radix Pulsatillae, Cortex Phellodendri, Rhizoma Coptidis, and Cortex Fraxini. A total of 266 compounds were authenticated in PD, including 57 in Radix Pulsatillae, 140 in Cortex Phellodendri, 48 in Rhizoma Coptidis, and 21 in Cortex Fraxini. Based on the OB of ≥ 30% and DL of ≥ 0.18, 65 bioactive compounds (11 in Radix Pulsatillae, 37 in Cortex Phellodendri, 14 in Rhizoma Coptidis, and 3 in Cortex Fraxini) were examined further (Table 1).
Prediction targets and C-D-T network analysis
In total, 180 target genes from the TCSP and DRUGBANK databases were screened. Furthermore, using “HCC” as a term for research, 6297 related genes were retrieved from GeneCards with a relevance score of ≥ 1 and 58 related genes were obtained from OMIM. Eventually, 6310 HCC-related genes were obtained from both GeneCards and OMIM (Fig. 2A). Comparing the 180 putative target genes of the compounds with the 6310 HCC-related genes, 152 common genes were identified as hub genes for subsequent study (Fig. 2B).
As shown in Fig. 2C, the C-D-T network comprised 192 nodes and 585 edges. The green circle represents genes, the pink triangle represents compounds, the blue diamond represents the four herbs of PD, and the red hexagon represents HCC. The top three compounds with the highest degree were quercetin (MOL000098 127), stigmasterol (MOL000449, 22), and isorhamnetin (MOL000354, 21), which may play key roles in mediating the anticancer effect in HCC.
PPI network analysis
To elucidate the interaction between HCC and PD-related targets, we analyzed the interaction between the two by establishing a PPI network. In addition, we analyzed 152 PD therapeutic targets using the STRING database, and the obtained data were analyzed by Cytoscape 3.7.1. As shown in Fig. 3A, the constructed PPI network contains 70 nodes and 109 edges. Table 2 lists the top 10 intensity data of the red nodes, which indicates that PD plays an important role in the treatment of HCC.
GO enrichment analysis
The GO and KEGG analyses were used to reveal functions and pathways related to the 152 anti-HCC PD-related target genes. In total, 2508 GO terms were found to be associated with the 152 anti-HCC PD target genes. Among these, 2281 terms were related to biological process, 67 to cell composition, and 160 terms to molecular function. After screening the top ten most obvious differences in GO analysis, the targets of PD were found mainly located in the membrane raft and membrane microdomain and were characterized by nuclear receptor activity. In addition, these targets of PD were involved in response to lipopolysaccharides, molecules of bacterial origin, a steroid hormone, the metabolic process of reactive oxygen species, and oxygen levels for the treatment of HCC (Fig. 3B).
C-P-T and KEGG enrichment analyses
To clarify the pathways associated with the anti-HCC effect of PD-related target genes, KEGG pathway analysis was performed. The results demonstrated that 145 of the 152 target genes were enriched in 168 pathways (P < 0.05). This indicated that the anti-HCC mechanism of PD could involve multiple target genes and pathways. The top 20 KEGG pathways with lower p values were visualized using a bar plot (Fig. 3C). Among the top 20 pathways, the TNF signaling pathway with a lower p-value was found to have a high association with carcinoma and was selected for subsequent validation (Supplementary Fig. 1). Besides, to elucidate the interrelationship among compounds, pathways, and targets of the top 20 pathways, the C-P-T network was constructed (Fig. 4), which included 154 nodes and 685 edges. The red circle indicates genes, the pink triangle indicates compounds, the blue diamond indicates the pathways, the green hexagon indicates the four herbs of PD, and the yellow octagon indicates HCC.
Molecular docking analysis
Combining the results of PPI, C-D-T network analysis, and the TNF signaling pathway, the targets JUN (PDB code:2P33) and AKT1 (PDB code:1h10) with a higher degree were selected for docking with quercetin. The results of quercetin docking on JUN and AKT1 are shown in Fig. 5A, B. Quercetin is linked via four hydrogen bonds and one H-π bond with Ser72, Gly147, Lys93, and Val78 in JUN, (Fig. 5C), whereas it formed two hydrogen bonds with Lys14 and Asn53 in AKT1 (Fig. 5D).
In summary, the results of molecular docking illustrated that compounds and protein targets interacted well via different bonds. The matching degree between compounds and target proteins assessed using binding energy revealed that the lower the binding energy, the greater the stability (Table 3). This suggested that compounds could bind well with the active site of protein targets.
PD inhibits HCC cell proliferation
HCC-LM3 and MHCC-97 H cells were used to analyze the inhibitory effect of PD on HCC proliferation. Different concentrations of PD (100–1000 µg/ml) were tested for different incubation periods (24, 48, and 72 h). The DMSO toxicity assay confirmed that at the maximum concentration of 0.5%, DMSO did not exert a significant cytotoxic effect on HCC-LM3 and MHCC-97 H cells (Supplementary Fig. 2A). The survival rate of HCC cells decreased with the dose and exposure time (Fig. 6A). The half inhibitory concentration of PD at 24, 48, and 72 h is shown in Table 4. Cell viability analysis showed that PD had a significant killing effect on HCC cells, and the effect was stronger after 48 and 72 h of incubation than that after 24 h, whereas there was no significant difference between the effect at 48 and 72 h. For the colony formation assay, HCC-LM3 and MHCC-97 H cells were treated with concentrations of 0,100, and 200 µg/ml (Fig. 6B). The results revealed that PD significantly reduced colony formation in both cell types and could inhibit the proliferation of HCC cells. In addition, cell cycle analysis revealed a delayed G1/S transition and cycle arrest in HCC-LM3 and MHCC-97 H cells following treatment with PD (Fig. 6C). Therefore, PD can effectively inhibit the occurrence and development of liver cancer.
PD promotes HCC cell apoptosis
Cell viability experiments showed that PD could induce the death of HCC cells. Therefore, the role of PD on HCC-LM3 and MHCC-97 H cell apoptosis was assessed using flow cytometry. Increasing the concentration of PD increased the apoptosis rate of HCC cells, reaching the maximum at 600 µg/ml. Moreover, apoptosis mainly comprised early apoptosis and late apoptosis (Fig. 7, Supplementary Fig. 2B). These findings demonstrate that PD can induce apoptosis in HCC cells; however, the specific mechanism requires further elucidation.
PD inhibits PI3K/AKT, TNF, and IL-17 signaling pathways in HCC cells
We further analyzed the potential molecular mechanisms underlying the effect of PD and signaling pathways of PD in HCC treatment using western blotting: PI3K/AKT, TNF, and IL-17 pathway. As shown in Fig. 8A, PI3K and AKT proteins in the PI3K/AKT pathways were significantly downregulated following PD treatment. Moreover, NF-κB and FOS proteins in the TNF pathway were down-regulated, whereas cleaved caspase 3 was upregulated following PD treatment (Fig. 8B). With regard to the IL-17 pathway, the levels of NFKBIA, ERK1/2, and MMP3 proteins were significantly downregulated following PD treatment (Fig. 8C). The results revealed that PI3K/AKT, TNF, and IL-17 signaling pathways play a key role in tumor proliferation and that PD may affect the activation of these pathways to inhibit the proliferation of HCC cells. However, further studies are needed to validate these findings.
The global prevalence of HCC is rising; however, its prognosis remains poor [4, 22]. As liver cancer is not easily detected during the early stages, patients often present with advanced disease. Moreover, there are limited effective treatments for advanced HCC . Recently, TCM is garnering increasing attention for its obvious pharmacological effects and low side effects. Specifically, TCM compounds may have potential applications in the treatment of cancer because they contain several effective active ingredients and have excellent anticancer effects. Therefore, it is essential to analyze the role of TCM in cancer. Several tools and methods have been developed to assess the role of treatments in diseases, particularly cancer. Network pharmacology is a discipline that studies the relationship between drugs and disease targets and is being widely used to analyze and predict new drug targets .
PD is a TCM compound that has been proven to have effects in ulcerative colitis [11, 24, 25], Crohn’s disease , and the immunogenic death of colorectal cancer cells in combination with 5-FU . PD has been used in Chinese medicine for thousands of years to treat diseases caused by bacteria. The primary active ingredient of PD, Radix Pulsatillae, mainly contains  saponins [27,28,29]. The saponin pulchinenoside C is the main effective ingredient of Pulsatilla. Moreover, pulchinenoside D is well known for its remarkable anticancer activity against H22 cells , HeLa cells , and colon cancer . Reportedly, Pulsatilla saponins can promote tumor cell death. Xue et al.  reported that pulchinenoside C can induce tumor cell apoptosis via the Bcl-2-caspase-3 pathway. Furthermore, pulchinenoside D can inhibit the PI3K–Akt signaling pathway to inhibit tumor growth and induce apoptosis . Cortex Phellodendri, a herbal medicine, has been used in China for more than 1,000 years to treat dysentery, jaundice, urinary tract infections, and rheumatoid arthritis . However, it is receiving increasing attention for the prevention of prostate cancer and osteoarthritis. Limonin, extracted from Cortex Phellodendri, has been proven to have significant anticancer efficacy , particularly by inducing apoptosis in human colon cancer cells [37, 38] and liver cancer cells . In addition, limonin can inhibit glycolysis in liver cancer cells, thus inducing apoptosis . Rhizoma Coptidis contains coptisine as its main active ingredient . Coptisine has anticancer, anti-inflammatory, coronary artery disease (CAD) improvement, and antibacterial potential . In addition, it has been found to downregulate the expression of PI3K and AKT and induce the apoptosis of colon cancer cells . Similarly, the anticancer effects of coptisine have been observed in liver, breast, pancreatic, and lung cancers [43,44,45]. Cortex Fraxini is a commonly used TCM compound that has been found to have disease-resistant anti-microbial, anti-tumor, and vascular protection functions . These results indicate that the main active ingredients of PD have anti-tumor effects. However, the antitumor efficacy of PD and the mechanism underlying this effect have not been studied yet.
Although researchers have assessed the antitumor efficacy of the individual components of PD, their efficacy has been limited. The present study is novel as it systematically analyzed the anticancer components of PD through network pharmacology to further study the mechanism and potential targets of PD for the treatment of HCC. Network pharmacology is a new research technology integrating systems biology, bioinformatics, and network science. It can analyze the molecular relationship between drugs and treatment objects from the perspective of the system and the whole biological network to reveal the systematic pharmacological mechanism of drugs . Network pharmacology involves drug targets, disease targets, and their pharmacological relationships and can be used to systematically study the effects of drugs on complex diseases. It can also be used to evaluate the drug compound–target–disease network, which is a novel way to study the mechanism of action of TCM compounds.
Tumor occurrence and development are often characterized by the induction of signaling pathways such as proliferation and invasion. As a result, these signaling pathways are promising targets for cancer treatment . The present study identified 266 compounds, of which 65 active components were identified through network pharmacology analysis. These included 11 compounds in Radix Pulsatillae, 37 in Cortex Phellodendri, 14 in Rhizoma Coptidis, and 3 in Cortex Fraxini. Of the 6310 HCC-related genes identified from GeneCards and OMIM, 152 common genes were selected for further analysis. Using PPI network analysis, the potential target genes of PD for the treatment of liver cancer were screened. Furthermore, the active ingredients and molecular mechanisms of PD were analyzed in the treatment of HCC.
The analysis demonstrated that the PI3K/AKT, TNF, and IL-17 pathways may be the key signaling pathways via which PD induces its effect in HCC. AKT signaling has been shown to promote tumor proliferation and metastasis. In addition, the PI3K/AKT signaling pathway can promote the growth and metastasis of tumor cells [48,49,50]. It additionally regulates cell differentiation, proliferation, migration, and apoptosis by inhibiting or activating a series of downstream substrates through phosphorylation [51, 52]. Alternatively, the TNF signaling pathway was also activated, as revealed by network pharmacology. TNF has been reported to trigger the activation of several pathways, including the NF-κB and MAPK signaling. Western blotting revealed that the expression of NF-κB and FOS was downregulated whereas that of Caspase 3 was upregulated. Reportedly, the upregulation of the NF-κB signaling pathway is related to tumor proliferation and metastasis . NF-κB has anti-apoptotic and pro-angiogenic functions in many cancers . In addition, FOS is a proto-oncogene that promotes tumor proliferation and metastasis . Therefore, NF-κB may be a key molecule that mediates the anticancer effect of PD in HCC, by regulating the proliferation and apoptosis of HCC. The study further revealed that PD treatment downregulated the expression of NFKBIA, ERK1/2, and MMP3. As a member of the MAPK signaling pathway, ERK1/2 can promote the proliferation and differentiation of tumor cells, leading to poor prognosis. Specifically, MMP3 is upregulated in malignant tumors and promotes tumor invasion and metastasis . Therefore, PD can downregulate the expression of NFKBIA, ERK1/2, and MMP3, thus inhibiting the proliferation of HCC and promoting apoptosis. Moreover, IL-17 and TNF signaling pathways are primarily involved in inflammation, tumor cell proliferation, and immune regulation. Dong et al. have demonstrated that downregulating IL17 mRNA levels inhibited the proliferation of tumor cells and reduced the expression of inflammatory cytokines . Moreover, Xun et al. have found that the IL-17 and TNF signaling pathways regulate the expression of PD-L1 in tumor cells, thus regulating the efficacy of drug therapy for tumors [58, 59]. Thus, PD may exert its anti-HCC effects via several pathways, including tumor proliferation, apoptotic, and inflammatory pathways. In addition, many other signaling pathways have been discovered and their role in PD-mediated treatment of HCC needs to be further explored. In conclusion, network pharmacology can effectively predict the antitumor effect of PD.
Based on network pharmacological analysis, this study confirmed the antitumor efficacy of PD components, which was further validated by in vitro assays. In addition, the study revealed the pathways involved in mediating the anti-HCC effects of PD. Further in vitro experiments confirmed that PI3K/AKT, TNF, and IL-17 pathways mediated the apoptotic effect of PD in HCC cells. These results clarify the mechanism of PD in the treatment of HCC and provide a new therapeutic potential for the treatment of malignant tumors with TCM.
The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.
McGlynn KA, Petrick JL, El-Serag HB. Epidemiol Hepatocellular Carcinoma Hepatol. 2021;73(1):4–13.
Kong FH, Ye QF, Miao XY, Liu X, Huang SQ, Xiong L, et al. Current status of sorafenib nanoparticle delivery systems in the treatment of hepatocellular carcinoma. Theranostics. 2021;11(11):5464–90.
Jemal A, Ward EM, Johnson CJ, Cronin KA, Ma J, Ryerson B et al. Annual Report to the Nation on the Status of Cancer, 1975–2014, Featuring Survival. J Natl Cancer Inst. 2017;109(9).
Chen W, Zhang Z, Fang X, Xiong L, Wen Y, Zhou J, et al. Prognostic value of the ALBI grade among patients with single hepatocellular carcinoma without macrovascular invasion. Med (Baltim). 2021;100(24):e26265.
Ahn JC, Teng PC, Chen PJ, Posadas E, Tseng HR, Lu SC, et al. Detection of circulating Tumor cells and their implications as a biomarker for diagnosis, prognostication, and therapeutic monitoring in Hepatocellular Carcinoma. Hepatology. 2021;73(1):422–36.
Kong FH, Miao XY, Zou H, Xiong L, Wen Y, Chen B, et al. End-stage Liver Disease score and future liver remnant volume predict post-hepatectomy Liver Failure in hepatocellular carcinoma. World J Clin Cases. 2019;7(22):3734–41.
Cao J, Kong FH, Liu X, Wang XB. Immunotherapy with dendritic cells and cytokine-induced killer cells for hepatocellular carcinoma: a meta-analysis. World J Gastroenterol. 2019;25(27):3649–63.
Hu B, Wang SS, Du Q. Traditional Chinese medicine for prevention and treatment of hepatocarcinoma: from bench to bedside. World J Hepatol. 2015;7(9):1209–32.
Hu B, An HM, Wang SS, Chen JJ, Xu L. Preventive and therapeutic effects of Chinese Herbal compounds against Hepatocellular Carcinoma. Molecules. 2016;21(2):142.
Jie Y, Yang X, Chen W. Pulsatilla Decoction combined with 5-Fluorouracil triggers immunogenic cell death in the Colorectal Cancer cells. Cancer Biother Radiopharm. 2021.
Wang X, Xu L, Wang T, Xu J, Fan F, Zhang Y et al. Pulsatilla decoction alleviates Colitis by enhancing autophagy and regulating PI3K–Akt–mTORC1 signaling pathway. Mol Med Rep. 2022;25(3).
Wang XM, Tian G, Duan QJ, Wu DQ, Shao J, Wang TM, et al. Therapeutic potential of n-butanol extract of Pulsatilla decoction in a murine model of ulcerative Colitis induced by DSS combined with Candida albicans colonization. Zhongguo Zhong Yao Za Zhi. 2018;43(14):2979–84.
Wu R, Xi Z, Liu M, Ren H, Dai R, Jiang X, et al. Pulsatilla Decoction and its bioactive component β-peltatin induce G2/M cell cycle arrest and apoptosis in Pancreatic cancer. Chin Med. 2023;18(1):61.
Liu H, Hu Y, Qi B, Yan C, Wang L, Zhang Y, et al. Network pharmacology and molecular docking to elucidate the mechanism of pulsatilla decoction in the treatment of colon Cancer. Front Pharmacol. 2022;13:940508.
Kim SY, Park C, Kim MY, Ji SY, Hwangbo H, Lee H et al. ROS-Mediated Anti-tumor Effect of Coptidis Rhizoma against Human Hepatocellular Carcinoma Hep3B Cells and xenografts. Int J Mol Sci. 2021;22(9).
Zhu M, Wang N, Tsao SW, Yuen MF, Feng Y, Wan TS, et al. Up-regulation of microRNAs, miR21 and miR23a in human Liver cancer cells treated with Coptidis rhizoma aqueous extract. Exp Ther Med. 2011;2(1):27–32.
Wang N, Feng Y, Lau EP, Tsang C, Ching Y, Man K, et al. F-actin reorganization and inactivation of rho signaling pathway involved in the inhibitory effect of Coptidis Rhizoma on hepatoma cell migration. Integr Cancer Ther. 2010;9(4):354–64.
Chen Y, Kern TS, Kiser PD, Palczewski K. Eyes on systems pharmacology. Pharmacol Res. 2016;114:39–41.
Jing C, Sun Z, Xie X, Zhang X, Wu S, Guo K, et al. Network pharmacology-based identification of the key mechanism of Qinghuo Rougan Formula acting on uveitis. Biomed Pharmacother. 2019;120:109381.
Huang S, Zhang Z, Li W, Kong F, Yi P, Huang J, et al. Network Pharmacology-Based Prediction and Verification of the active ingredients and potential targets of Zuojinwan for treating Colorectal Cancer. Drug Des Devel Ther. 2020;14:2725–40.
Xu X, Xu H, Shang Y, Zhu R, Hong X, Song Z, et al. Development of the general chapters of the Chinese pharmacopoeia 2020 edition: a review. J Pharm Anal. 2021;11(4):398–404.
Llovet JM, De Baere T, Kulik L, Haber PK, Greten TF, Meyer T, et al. Locoregional therapies in the era of molecular and immune treatments for hepatocellular carcinoma. Nat Rev Gastroenterol Hepatol. 2021;18(5):293–313.
Grimstad T, Kvivik I, Kvaløy JT, Aabakken L, Omdal R. Heat-shock protein 90α in plasma reflects severity of fatigue in patients with Crohn’s Disease. Innate Immun. 2020;26(2):146–51.
Hua YL, Ma Q, Zhang XS, Jia YQ, Peng XT, Yao WL, et al. Pulsatilla Decoction can treat the dampness-heat Diarrhea Rat Model by regulating glycerinphospholipid metabolism based Lipidomics Approach. Front Pharmacol. 2020;11:197.
Hua YL, Ma Q, Li W, Zhang XS, Cheng XH, Jia YQ, et al. Metabolomics analysis of Pulsatilla decoction on treatment of wetness-heat-induced diarrhea in rats based on UPLC-Q/TOF-MS/MS. Biomed Chromatogr. 2019;33(11):e4629.
Liu J, Zhang L, Wang Z, Chen S, Feng S, He Y, et al. Network Pharmacology-based strategy to identify the pharmacological mechanisms of Pulsatilla Decoction against Crohn’s Disease. Front Pharmacol. 2022;13:844685.
Zhong J, Tan L, Chen M, He C. Pharmacological activities and molecular mechanisms of Pulsatilla saponins. Chin Med. 2022;17(1):59.
Jin MM, Zhang WD, Xu YM, Du YF, Sun Q, Guo W, et al. Simultaneous determination of 12 active components in the roots of Pulsatilla chinensis using tissue-smashing extraction with liquid chromatography and mass spectrometry. J Sep Sci. 2017;40(6):1283–92.
Ling Y, Lin Z, Zha W, Lian T, You S. Rapid Detection and Characterisation of Triterpene Saponins from the Root of Pulsatilla chinensis (Bunge) Regel by HPLC-ESI-QTOF-MS/MS. Phytochem Anal. 2016;27(3–4):174–83.
Fang Y, Hu D, Li H, Hu J, Liu Y, Li Z, et al. Synthesis, Biological evaluation, and Mode of Action of Pulsatilla Saponin D Derivatives as promising Anticancer agents. Front Pharmacol. 2019;10:1208.
Zhang Y, Bao J, Wang K, Jia X, Zhang C, Huang B, et al. Pulsatilla Saponin D inhibits autophagic flux and synergistically enhances the Anticancer activity of chemotherapeutic agents against HeLa cells. Am J Chin Med. 2015;43(8):1657–70.
Son MK, Jung KH, Hong SW, Lee HS, Zheng HM, Choi MJ, et al. SB365, Pulsatilla saponin D suppresses the proliferation of human colon Cancer cells and induces apoptosis by modulating the AKT/mTOR signalling pathway. Food Chem. 2013;136(1):26–33.
Xue S, Zhou Y, Zhang J, Xiang Z, Liu Y, Miao T, et al. Anemoside B4 exerts anti-cancer effect by inducing apoptosis and autophagy through inhibiton of PI3K/Akt/mTOR pathway in hepatocellular carcinoma. Am J Transl Res. 2019;11(4):2580–9.
Hong SW, Jung KH, Lee HS, Choi MJ, Son MK, Zheng HM, et al. SB365 inhibits angiogenesis and induces apoptosis of hepatocellular carcinoma through modulation of PI3K/Akt/mTOR signaling pathway. Cancer Sci. 2012;103(11):1929–37.
Sun H, Wang H, Zhang A, Yan G, Han Y, Li Y, et al. Chemical discrimination of Cortex Phellodendri amurensis and Cortex Phellodendri Chinensis by Multivariate Analysis Approach. Pharmacogn Mag. 2016;12(45):41–9.
Fan S, Zhang C, Luo T, Wang J, Tang Y, Chen Z et al. Limonin: a review of its Pharmacology, Toxicity, and Pharmacokinetics. Molecules. 2019;24(20).
Qian P, Jin HW, Yang XW. New limonoids from Coptidis Rhizoma-Euodiae Fructus couple. J Asian Nat Prod Res. 2014;16(4):333–44.
Chidambara Murthy KN, Jayaprakasha GK, Kumar V, Rathore KS, Patil BS. Citrus limonin and its glucoside inhibit colon adenocarcinoma cell proliferation through apoptosis. J Agric Food Chem. 2011;59(6):2314–23.
Rahman A, Siddiqui SA, Jakhar R, Kang SC. Growth inhibition of various human cancer cell lines by imperatorin and limonin from poncirus trifoliata rafin. Seeds Anticancer Agents Med Chem. 2015;15(2):236–41.
Yao J, Liu J, Zhao W. By blocking hexokinase-2 phosphorylation, limonin suppresses Tumor glycolysis and induces cell apoptosis in hepatocellular carcinoma. Onco Targets Ther. 2018;11:3793–803.
Wu J, Luo Y, Deng D, Su S, Li S, Xiang L, et al. Coptisine from Coptis chinensis exerts diverse beneficial properties: a concise review. J Cell Mol Med. 2019;23(12):7946–60.
Cao Q, Hong S, Li Y, Chen H, Shen Y, Shao K, et al. Coptisine suppresses Tumor growth and progression by down-regulating MFG-E8 in Colorectal cancer. RSC Adv. 2018;8(54):30937–45.
Huang Z, Dong F, Li S, Chu M, Zhou H, Lu Z, et al. Berberine-induced inhibition of adipocyte enhancer-binding protein 1 attenuates oxidized low-density lipoprotein accumulation and foam cell formation in phorbol 12-myristate 13-acetate-induced macrophages. Eur J Pharmacol. 2012;690(1–3):164–9.
Zhou L, Yang F, Li G, Huang J, Liu Y, Zhang Q, et al. Coptisine induces apoptosis in human hepatoma cells through activating 67-kDa laminin Receptor/cGMP signaling. Front Pharmacol. 2018;9:517.
Han B, Jiang P, Xu H, Liu W, Zhang J, Wu S, et al. 8-Cetylcoptisine, a new coptisine derivative, induces mitochondria-dependent apoptosis and G0/G1 cell cycle arrest in human A549 cells. Chem Biol Interact. 2019;299:27–36.
Fang LH, Lv Y, Du GH. [Progress in study of pharmacological effect of cortex Fraxini]. Zhongguo Zhong Yao Za Zhi. 2008;33(23):2732–6.
Moeini A, Cornella H, Villanueva A. Emerging Signaling pathways in Hepatocellular Carcinoma. Liver Cancer. 2012;1(2):83–93.
ShuoWang, Song Z, Gong X, Ou C, Zhang W, Wang J et al. Chloroform extract from Sophora Tonkinensis Gagnep. Inhibit proliferation, migration, invasion and promote apoptosis of nasopharyngeal carcinoma cells by silencing the PI3K/AKT/mTOR signaling pathway. J Ethnopharmacol. 2021:113879-.
Yu X, Li S, Pang M, Du Y, Xu T, Bai T et al. TSPAN7 exerts Anti-tumor effects in Bladder Cancer through the PTEN/PI3K/AKT pathway. Front Oncol. 2021;10.
Zhang M, Lin H, Ge X, Xu Y. Overproduced CPSF4 promotes Cell Proliferation and Invasion via PI3K-AKT signaling pathway in oral squamous cell carcinoma. Journal of oral and maxillofacial Surgery. official journal of the American Association of Oral and Maxillofacial Surgeons; 2021.
Danielsen SA, Eide PW, Nesbakken A, Guren T, Leithe E, Lothe RA. Portrait of the PI3K/AKT pathway in Colorectal cancer. Biochim Biophys Acta. 2015;1855(1):104–21.
Bishnupuri KS, Alvarado DM, Khouri AN, Shabsovich M, Chen B, Dieckgraefe BK, et al. IDO1 and kynurenine pathway metabolites activate PI3K-Akt signaling in the neoplastic Colon epithelium to Promote Cancer Cell Proliferation and inhibit apoptosis. Cancer Res. 2019;79(6):1138–50.
Mitrakas L, Gravas S, Papandreou C, Koukoulis G, Karasavvidou F, Dimakopoulos G, et al. Primary high-Grade non-muscle-invasive Bladder Cancer: high NFκB expression in Tumor specimens distinguishes patients who are at risk for Disease Progression. Pathol Oncol Res. 2019;25(1):225–31.
Mortezaee K, Najafi M, Farhood B, Ahmadi A, Shabeeb D, Musa AE. NF-κB targeting for overcoming Tumor resistance and normal tissues toxicity. J Cell Physiol. 2019;234(10):17187–204.
Kong F, Ye S, Zhong Z, Zhou X, Zhou W, Liu Z, et al. Single-cell transcriptome analysis of chronic antibody-mediated rejection after renal transplantation. Front Immunol. 2021;12:767618.
Liang M, Wang J, Wu C, Wu M, Hu J, Dai J, et al. Targeting matrix metalloproteinase MMP3 greatly enhances oncolytic virus mediated Tumor therapy. Transl Oncol. 2021;14(12):101221.
Pan D, Huang B, Gan Y, Gao C, Liu Y, Tang Z. Phycocyanin ameliorates Colitis-Associated Colorectal Cancer by regulating the gut microbiota and the IL-17 signaling pathway. Mar Drugs. 2022;20(4).
Wang X, Yang L, Huang F, Zhang Q, Liu S, Ma L, et al. Inflammatory cytokines IL-17 and TNF-α up-regulate PD-L1 expression in human prostate and colon Cancer cells. Immunol Lett. 2017;184:7–14.
Yaghoubi N, Soltani A, Ghazvini K, Hassanian SM, Hashemy SI. PD-1/ PD-L1 blockade as a novel treatment for Colorectal cancer. Biomed Pharmacother. 2019;110:312–8.
Hunan Innovation Platform and Talent Plan-2020 Huxiang Young Talents: 2020RC3062.
Ethics approval and consent to participate
This project was approved by the Medical Ethics Committee of the Second Xiangya Hospital of Central South University. All experiments and strategies follow ethical principles.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Liu, K., Cao, Z., Huang, S. et al. Mechanism underlying the effect of Pulsatilla decoction in hepatocellular carcinoma treatment: a network pharmacology and in vitro analysis. BMC Complement Med Ther 23, 405 (2023). https://doi.org/10.1186/s12906-023-04244-w