Miaolun Yang
Qian Yan
Yuehua Luo1†
Baoqian Ye
Xiongwen Wang- 1The First Clinical Medical School, Guangzhou University of Chinese Medicine, Guangzhou, China
- 2The First Affiliated Hospital, Guangzhou University of Chinese Medicine, Guangzhou, China
Background: Hepatocellular carcinoma (HCC) is a malignant tumor harmful to human health. Ganji Fang (GJF) has good clinical efficacy in the treatment of HCC, but its mechanism is still unclear.
Objective: The aim of this study was to investigate the mechanism of action of GJF in the treatment of HCC through network pharmacology, molecular docking and in vitro experiments.
Methods: A series of network pharmacology methods were used to identify the potential targets and key pathways of GJF in the treatment of HCC. Then, molecular docking technology was used to explore the binding ability of key active ingredients and targets in GJF. Multiple external databases were used to validate the key targets. In in vitro experiments, we performed MTT assays, wound-healing assays, cell cycle assays, apoptosis assays and RT‒qPCR to verify the inhibitory effect of GJF on the Human hepatoma G2 (HepG2) cells.
Result: A total of 162 bioactive components and 826 protein targets of GJF were screened, and 611 potential targets of HCC were identified. Finally, 63 possible targets of GJF acting on HCC were obtained. KEGG enrichment analyses showed that the top five pathways were the cell cycle, cellular senescence, p53 signaling pathway, PI3K/Akt signaling pathway, and progesterone-mediated oocyte maturation. Among them, we verified the PI3K/Akt signaling pathway. CCNE1, PKN1, CCND2, CDK4, EPHA2, FGFR3, CDK6, CDK2 and HSP90AAI were enriched in the PI3K/Akt pathway. The molecular docking results showed that the docking scores of eight active components of GJF with the two targets were all less than -5.0, indicating that they had certain binding activity. In vitro cell experiments showed that GJF could inhibit the proliferation and migration of HepG2 cells, block the cell cycle and induce apoptosis of HepG2 cells, which may be related to the PI3K/Akt signaling pathway. In summary, EPHA2 may be an important target of GJF in HCC, and pachymic acid may be an important critical active compound of GJF that exerts anticancer activity.
Conclusion: In general, we demonstrated, for the first time, that the molecular mechanism of GJF in HCC may involve induction of G0/G1 phase cycle arrest through inhibition of the PI3K/Akt signaling pathway and promote apoptosis of hepatoma cell lines. This study provides a scientific basis for the subsequent clinical application of GJF and the in-depth study of its mechanism.
1 Introduction
Primary liver cancer (PLC) includes hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma as well as other rare types (Bray et al., 2018). HCC, the most prevalent type of PLC, accounts for approximately 90% of cases (Gunasekaran et al., 2020). HCC ranks fifth in terms of global cases and second in terms of deaths for males (Bray et al., 2018). Studies have shown that the incidence of HCC has stabilized in men after decades of steep increase but continues to rise in women by >2% annually (Siegel et al., 2021). Surgical removal involving resection and transplantation remains the only curative treatment for most types of liver tumors (Petrowsky et al., 2020). In addition, chemotherapy, radiotherapies, radiofrequency ablation (RFA), percutaneous ethanol injection (PEI), transarterial chemoembolization (TACE), targeted therapy, and immunological therapy are also important treatment methods (Petrowsky et al., 2020). However, all of these treatments have side effects that affect patients’ lives. As the use of immune checkpoint inhibitors continues to expand, this treatment will become more crucial than early detection, and management of these side effects becomes paramount to maximize the duration of treatment while minimizing toxicity for patients (Hussaini et al., 2021).
Although the effect of traditional Chinese medicine (TCM) on tumor shrinkage is not as obvious as that of modern Western medicine, reducing toxicity and increasing the efficacy of TCM can help to reduce the side effects, improve the curative effect and improve the quality of life of patients (Z. Ma et al., 2019; Ye et al., 2015). Evidence-based functions of Chinese herbal medicines against cancer can be summarized as enhancement of natural killer (NK) cell activity or relative percentage; prevention of tumor growth and metastasis; relief of side effects or complications of therapeutic strategies (i.e., chemotherapy, radiotherapy, and resection) (H. Liu et al., 2021). Some studies have shown that adjuvant therapy with TCM may prolong the median overall survival time and reduce the mortality of HCC patients (X. Liu et al., 2019). GJF is an established prescription drug with beneficial gas spleen detoxification and anticancer functions. GJF is commonly used to treat HCC in clinical practice, and some patients with advanced liver cancer even choose to be treated with pure Chinese medicine. Previous clinical and animal studies have preliminarily explored whether GJF may have anti-inflammatory effects, protect the liver, inhibit tumor cell proliferation, and induce tumor cell apoptosis. Therefore, this study will further explore the clinical application value and possible mechanism of GJF along the lines of previous research. TCM compounds are difficult to study because they are composed of multiple herbs, contain complex components, and their mechanism of action is unclear.
Network pharmacology uses drug, compound, gene, and disease database information to construct drug-target, target-disease, and drug-disease interaction networks to reveal the complex mechanisms of TCM formulations that have multiple targets and multicomponent characteristics (S. Li & Zhang, 2013), which provides a new strategy for the study of complex traditional Chinese medicine systems. The comprehensive, systematic and holistic nature of network pharmacology are consistent with the multicompound, multitarget and multipathway characteristics of Chinese medicine (Luo et al., 2020). Many researchers use network pharmacology to predict the active phytochemicals and molecular mechanisms of traditional Chinese medicine in the treatment of HCC. Many researchers use network pharmacology to predict the active phytochemicals and molecular mechanisms of traditional Chinese medicine in the treatment of HCC (Khan & Lee, 2022). In this paper, the possible pathways and targets of GIF for HCC treatment were predicted mainly through network pharmacology methods and then verified by molecular docking, external databases, experimental research and other methods. The detailed flow chart of the study is shown in Figure 1.
FIGURE 1. Overall workflow of this study (GJF, Ganji Fang; TCMSP, Traditional Chinese Medicine System Pharmacology database; STP, SwissTargetPrediction database; HCC, Hepatocellular carcinoma).
2 Materials and methods
2.1 Screening active compounds GJF
We obtained the compounds of each drug of GJF from the Traditional Chinese Medicine System Pharmacology Database and Analysis Platform (TCMSP, http://lsp.nwu.edu.cn/tcmsp.php, accessed on 12 May 2020). In addition, the potential active compounds of each Chinese medicine were also supplemented through a literature search. Oral bioavailability (OB) and drug-likeness (DL) are some of the most important pharmacokinetic parameters in the process of absorption, distribution, metabolism, and excretion (ADME). We screened compounds on the condition that oral utilization of compounds OB ≥ 30% and DL ≥ 0.18.
2.2 Prediction of drug targets for GJF
The 3D structures of the active ingredients of each compound and canonical SMILES format were downloaded from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/, accessed on 20 May 2020). The potentially effective action targets of GJF were obtained using the SwissTargetPrediction database (STP, http://www.swisstargetprediction, accessed on 20 May 2020) (Daina et al., 2019). We chose the target probability to be greater than 0, and if there were too many predicted targets, we chose the top 100 targets ranked according to the target probability. When SwissTargetPrediction could not predict relevant targets, we used the TMSCP database to make a supplement prediction.
2.3 Extracting the targets of hepatocellular carcinoma
We consolidated HCC-related targets through the Databases Liverome (http://liverome.kobic.re.kr/index.php, accessed on 23 May 2020), OncoDB. HCC databases (http://www.oncodb.hcc.ibms.sinica.edu.org/, accessed on 23 May 2020), GeneCards databases (https://www.genecards.org/, accessed on 23 May 2020), and OMIM databases (http://www.ncbi.nlm.nih.gov/, accessed on 23 May 2020) (Su et al., 2007; Lee et al., 2011; Paolacci et al., 2019).
Additionally, the mRNA sequencing data of HCC were downloaded for the PI3K/Akt signaling pathway on The Cancer Genome Atlas (TCGA) dataset (http://cancergenome.nih.gov/, accessed on 23 May 2020), including 370 primary HCC tissues and 50 normal tissues and the corresponding clinical follow-up data. To improve the accuracy of the data, we preprocessed the dataset, including deleting the sites where 70% methylation levels were unavailable and samples expressing more than 30% missing values. Genes with zero RPKM expression in the sample were also excluded (J. Xu et al., 2015).
The data were standardized using the RMA algorithm in the “Limma” package in R software (Smyth, 2004; Law et al., 2014; Ritchie et al., 2015). A false detection rate (FDR) method was used to regulate the p-value. The screening criteria for differentially methylated genes (DMGs) were FDR <0.05 and β values >0.2, HCC in cancerous tissue and normal liver tissue differentially expressed genes (Differentially expressed genes, DEGs) in line with∣1og 2 (fold change)∣> 1 and FDR <0.05 standard. We intersected HCC-related targets in the public database with DMGS and DEGs in the TCGA database, regarded this intersection as possible key targets in HCC, and mapped these targets to related targets predicted by compounds in GJF as potential targets of GJF for HCC.
2.4 Construction of networks and analysis
The key targets of GJF for HCC were imported into the STRING database to construct a PPI network, and the confidence level was set as ≥0.70. Cytoscape (Version: 3.9.1) software was used for visualization. The key genes and their relationship network of TCM, drug active ingredient, and drug and disease crossing were imported into Cytoscape 3.9.1 software to construct a herb-compound-target network of GJF.
2.5 Gene ontology and pathway enrichment analysis for HCC-Related targets of GJF
Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed using the “cluster profile” package in R software for better biological interpretation (Yu et al., 2012). GO and KEGG enrichment analyses of the potential targets were performed using the database for annotation, visualization, and integrated discovery (DAVID, https://david.ncifcrf.gov/, accessed on 26 May 2021). GO enrichment analysis included biological processes (BP), molecular functions (MF), and cellular components (CC).
2.6 Molecular docking
The RCSB PDB database (https://www.rcsb. org/, accessed on 10 June 2021) to retrieve and download 3D structure files for key target proteins, while using the PubChem database (https://pubchem.ncbi.nlm.nih.gov/, accessed on 12 June 2021) to download 3D structure files for active compounds. For docking analysis, all protein and molecular files were converted into PDBQT format with all water molecules excluded, and polar hydrogen atoms were added. The grid box was centered to cover the domain of each protein and to accommodate free molecular movement. The grid box was set to 30 Å × 30 Å × 30 Å, and the grid point distance was 0.05 nm. Molecular docking studies were performed by Autodock Vina 1.2.2 (http://autodock.scripps.edu/, accessed on 12 June 2021).
2.7 Validation of external databases
For the above key targets, we performed a series of validations in the TCGA database. The cBioPortal (http://cbioportal.org, accessed on 16 June 2021) genomics database (J. Gao et al., 2013) demonstrated the association of genetic mutation status, mRNA, and DNA methylation. The GEPIA (Gene Expression Profiling Interactive Analysis, http://gepia.cancer-pku.cn/, accessed on 16 June 2021) database was used to examine the prognostic differences between high and low gene levels. The (T. Li et al., 2017) database was used to verify the difference in mRNA level expression.
2.8 In Vitro validation cell culture
The HCC cell line HepG2, which was purchased from Guangzhou Cellcook Biotech Co., Ltd., was routinely cultured in DMEM (Gibco, Grand Island, NY, United States) containing 10% FBS (Thermo Fisher Scientific, Rochester, NY, United States). Cells were cultured at 37 °C in a 5% CO2 incubator, subcultured at a ratio of 1:2, and used for subsequent experiments after reaching the logarithmic growth phase. In this study, we set the 72-hour IC50 to a medium dose, half the 72-hour IC50 to a low dose, and two times the 72-hour IC50 to a high dose.
2.9 The source of GJF and preparation of lyophilized powder
GJF was purchased from the Outpatient Pharmacy of the First Affiliated Hospital of Guangzhou University of Traditional Chinese Medicine. The dose of 1 fu of GJF was 241 g (12 g bupleurum, 12 g paeoniae, 30 g codonopsis, 15 g poria, 15 g atractylodes, 20 g Rhizoma lucidum, 30 g ligustrum scutellariae, 12 g Scutellaria scutellariae, 20 g Rhizoma barbiliae, 30 g Augustaloza, 15 g Zedoary Zedoaria, 30 g Pellucidaria). Deionized water was added at 10x (approximately 2400 ml), sampled was cooked twice and centrifuged at 1000 r for 5 min. The supernatant was collected and refrigerated at -80 °C for 6–8 h until it was completely frozen. Then, the samples were quickly put into a lyophilized powder machine. After the preparation of lyophilized powder, complete culture medium was used to prepare 50 mg/ml liver-product solution.
2.10 MTT assay
The logarithmic phase cells were collected, the cell density was adjusted to 5 ˣ 104/ml, and the cells were inoculated in three 96-well plates. Then, 100 µL of GJF solution with different concentration gradients was added to the cell culture box. One plate was removed at 24 h, 48 h, and 72 h, and a mixture of MTT solution and basic medium (ratio: 1:4) was added under closed light conditions, 100 µL for each well. After further incubation for 4 h, 150 µL DMSO solution was added to each well. The wavelength of the enzyme plate analyzer was set at 490 nm, and the OD values of 3 96-well plates were measured at different times.
2.11 Wound-healing assay
HepG2 cells in logarithmic growth were cultured in a 6-well plate (2 × 105 cells/well). After the cells reached a subconfluency state (approximately 80%), a 2-mm scratch was made along the midline of the plate. Cells were assigned to the blank serum group (blank serum) (cells cultured with normal culture medium supplemented with blank serum), GJF low-dose group (GJF-L), and GJF high-dose group (GJF-H), with two replicate wells set up for each group. When cells were cultured for 0, 12, and 24 h in drug-containing or drug-free media, they were photographed under an inverted microscope (Olympus, Tokyo, Japan), and the rate of cell migration was calculated.
2.12 Cell cycle detection
The cells were inoculated in a 6-well plate, divided into control group, GJF of low, medium, and high concentration group, allowed to adhere and grown to a certain density, and then a different concentration of GJF frozen-dried powder solution was added to each well. Twenty-four hours later the cells were added to 1 ml of ice bath precooled with 70% ethanol and fixed at 4 °C for 2 h or more. Then, the cells were precipitated and resuspended. Next, 0.5 ml of cyclosome iodide dyeing liquid was added to each group, incubated at 37 °C in a bath for 30 min in the dark, and stored at 4 °C or on an ice bath for 24 h to complete the flow test.
2.13 Apoptosis detection
The cells at several stages were inoculated in a 6-well plate, divided into control group, GJF of low, medium, and high concentration group, 24 h after taking 5–10 million suspended cells, the cells were centrifuged at 1000r for 5 min, after discarding, 195 µL Annexin V- FITC binding fluid was added to gently resuspend the cells, 5 µL Annexin V-FITC was added, 10u1 iodide dye was added, and cells were incubated at room temperature (20°C–25°C) for 10–20 min, then placed in an ice bath. The cells were tested on a flow cytometer as soon as possible within an hour.
2.14 RT‒qPCR microarray analysis
Total RNA was extracted from experimental cells using a universal RNA extraction kit. Then, EVO-M-MLVRT Master Mix was used for reverse transcription into cDNA. The reaction system was prepared according to the instructions of the SYBR premixed RT‒qPCR kit for RT‒qPCR amplification. The RT‒qPCR conditions were as follows: pre-denaturation at 95 °C for 30 s, denaturation at 95°C for 30 s, and annealing at 60 °C for 30 s, for a total of 40 cycles. The experiment was repeated 3 times. Gene expression was evaluated by RT-qPCR, using the 2 -∆∆CT method.
2.15 Statistical analysis
Measurement data are expressed as the mean taxi standard deviation (X ± S) and were statistically processed using SPSS (version: 26.0) software. When the experimental data met the normal distribution and the variance was consistent, one-way analysis of variance and two-way analysis of variance were used, and a non-parametric test was used if the data did not meet the normal distribution. p < 0.05 indicates that the difference is statistically significant. The above experiments were repeated three times, with the results of the cell cycle and apoptosis experiments presented in FlowJo (version 10.7.1) and the rest in GraphPad Prism (version 8.0).
3 Results
3.1 Active ingredients and prediction of a protein target of GJF
A total of 162 bioactive components and 826 protein targets of GJF were obtained from both TCMSP databases and a literature search (Supplementary Table S1). Radix Bupleuri was supplemented with Saponin A (Wen-Sheng, 2003) and Saponin D (Ren et al., 2019), Atractylodide I, Atractylodide II, and Atractylodide III were added to Atractylodes macrocephala Koidz (Zhu et al., 2018), raffinose was added to Codonopsitis Radix (S. Gao et al., 2019), Paeoniae Radix Alba was supplemented with oxypaeoniflorin, (M. Lu et al., 2019), Eclalbasaponin I and ursolic acid were added to Eclipse Herba (Q. M. Liu et al., 2012), oleanic acid was added to Fructus Ligustri Lucidi, Phytodolor and rosmarinic acid were added to Herba Sarcandrae, Scutellariae Barbatae Herba was supplemented with scutellarin, curcumin was added to Curcumae Rhizoma, and Akebiae Frucyus was added to Calceolarioside B.
3.2 Prediction results of hepatocellular carcinoma targets
We screened 5,714 genes in four public databases. The mRNA and methylation sites differentially expressed in the TCGA were further screened. According to the screening criteria for the differentially expressed genes, we obtained 7645 differentially expressed genes, and after the conversion of the differentially methylated sites, we identified 8938 differentially methylated genes, which we display in the form of heatmaps and volcano plots. (Figures 2A–C).
FIGURE 2. (A) Heatmap of HCC-associated differential methylation sites in the TCGA database (B) Heatmap of HCC-related differentially expressed genes in the TCGA database (C) Volcano plot of HCC-related differentially expressed genes in the TCGA (D) Venn diagram of the intersection of HCC-related targets (DEGs, differentially expressed genes; DMGs, differentially methylated genes).
3.3 GO and KEGG enrichment analyses
The 611 potential targets associated with HCC were mapped to the targets predicted by GJF (Figure 2D), and 63 intersections were obtained, which were the possible targets of GJF acting on HCC. The 63 genes were analyzed by GO and KEGG (Figure 3). According to the number of genes enriched in each item, the results showed that the top five biological processes involved in these genes were epithelial cell proliferation, regulation of epithelial cell proliferation, regulation of cell cycle phase transition, gland development, and embryonic organ development; in terms of cell composition, they were mainly involved in cell-substrate junction, cell-substrate adherens junction, focal adhesion, adherens junction, and extracellular matrix composition. In the molecular function category, they were mainly associated with cell adhesion molecule binding, proximal promoter sequence-specific DNA binding, chromatin binding, RNA polymerase II proximal promoter sequence-specific DNA binding and laminin binding. In addition, the top five pathways were the cell cycle (Gotoh et al., 2003), cellular senescence, p53, signaling pathway, PI3K/Akt signaling pathway, and progesterone-mediated oocyte maturation (Supplementary Table S2).
FIGURE 3. GO and KEGG analysis diagram of the intersection of targets of GJF in HCC (A) The analysis diagrams of the biological process (B) Diagrams of the analysis of cell components (C) The analysis diagrams of molecular functions (D) KEGG analysis diagram (GO, gene ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes).
3.4 PPI network and herb-compound-target network of GJF
We performed PPI network analysis of 63 intersection genes between HCC (Figure 4A). GJF and the top 10 genes were HSP90AA1, SRC, CDK1, CCNA2, CDK4, PLK1, CDK2, TOP2A, AURKB and CDK6. Based on the results of the KEGG analysis, we verified the role of the PI3K/Akt signaling pathway in GJF. KEGG analysis showed that CCNE1, PKN1, CCND2, CDK4, EPHA2, FGFR3, CDK6, CDK2, and HSP90AA1 were enriched in the PI3K/Akt pathway. We constructed a network diagram of “TCM-compound-target-pathway” (Figure 4B). The network consisted of 235 nodes (12 TCM name nodes, 162 active ingredient nodes, 63 gene target nodes) and 1225 edges.
FIGURE 4. (A) The PPI network of 63 nodes (genes) (The darker the color, the larger the circle, the more connected; The 10 targets in the innermost layer are the key targets with the highest connection degree) (B) Herb-compound-target network of GJF (green represents the name of traditional Chinese medicine, orange represents the active ingredient, blue represents the key gene target, and red represents Ganji Fang).
3.5 Molecular docking
We conducted a molecular docking analysis between HSP90AA1, anhydroicaritin, cubebin, (2R)-7-hydroxy-5-methoxy-2-phenylchroman-4-one, paeoniflorin (PF), albiflorin_qt, and bisdemethoxycurcumin to screen the compounds with the strongest binding ability to HSP90AA1. In addition, PF, albiflorin_qt, pachymic acid (PA), and (7,9 (11)-dehydropachymic acid were used to screen the compounds that acted most closely with EPHA2. The results showed that the docking scores of the eight active components of GJF with the two targets were all less than -4.0, indicating that they had certain binding activity. As shown in Table 1, the strongest binding energy between PF and HSP90AA1 is -7.3 and that between PA and EPHA2 was -8.8. Therefore, these key components may be potential active compounds of GJF in the treatment of HCC (Figure 5).
TABLE 1. Docking scores of key compounds and targets in hepatocellular formulas.
FIGURE 5. Molecular docking results. HSP90AA1 binds with MOL004373 (2D and 3D) (A, B), MOL013187 (2D and 3D) (C, D), MOL000228 (2D and 3D) (E, F), MOL001924 (2D and 3D) (G, H), MOL001928 (2D and 3D) (I, J), MOL000940 (2D and 3D) (K, L). EPHA2 binds with MOL001924 (2D and 3D) (M, N), MOL001928 (2D and 3D) (O, P), MOL000289 (2D and 3D) (Q, R), MOL000276 (2D and 3D) (S, T).
3.6 Validation of EPHA2 and HAP90AA1 in the TCGA database
The TIMER database showed that EPHA2 and HSP90AA1 were differentially expressed in multiple tumors, and EPHA2 was expressed at low levels in HCC. The mRNA expression levels of EPHA2 and HSP90AA1 were strongly correlated with the DNA methylation level (Figure 6), suggesting that the regulatory effect of GJF on EPHA2 and HSP90AA1 might be caused by the regulation of their methylation status by methyltransferase and then the regulation of their expression. However, HSP90AA1 was hypermethylated and highly expressed in HCC, and there was a negative correlation between the two, which may require further experimental confirmation. Our study showed that low expression of EPHA2 may have a better prognosis, but the difference was not statistically significant (Figure 7A, p = 0.35). In multiple TCGA datasets, the mutation rate of EPHA2 was 1.9% and that of HSP90AA1 was 1.4% (Figures 7C, D), suggesting that these two molecules have a certain application value, but further basic experiments may be needed to confirm this.
FIGURE 6. (A, B) Differential expression of EPHA2 and HSP90AA1 in various tumors. (C, D) Correlation analysis of mRNA expression levels and DNA methylation of EPHA2 and HSP90AA1 in HCC.
FIGURE 7. (A, B) Survival analysis of high and low expression populations of EPHA2 and HSP90AA1. (In the figure, the red curve refers to the population with high expression in the TCGA database, and the blue curve refers to the population with low expression in the TCGA database. The P value represents the statistical difference in the survival curve between the two groups). (C, D) Genetic mutation analysis of EPHA2 and HSP90AA1 in various TCGA datasets. (C) The percentage of total mutations of EPHA2 and HSP90AA1. (D) The mutations of EPHA2 and HSP90AA1 in each sample.
3.7 Effects on cell proliferation and migration
The MTT experiment results showed that a certain concentration of GJF freeze-dried powder solution inhibited the proliferation of HepG2 cell strains. With increasing drug concentration and time, the inhibition gradually increased in a time- and concentration-dependent manner (Figure 8A). The IC50 values of 24H, 48H, and 72H in the HepG2 cell line were 1.91 mg/ml, 1.56 mg/ml, and 1.12 mg/ml, respectively. The GJF freeze-dried powder solution was divided into a low concentration group (0.5 mg/ml), medium concentration group (1.0 mg/ml), and high concentration group (3.0 mg/ml).
FIGURE 8. GJF inhibited HepG2 cells growth in vitro. (A) The effect of GJF on the survival rate of HepG2 cells. (B) Effect of GJF on the cell mobility of HepG2 cells. (C) GJF induces hepatocellular arrest of HepG2 cells cycles. (D) GJF induces apoptosis in HepG2 cells. (E) Changes in the mRNA expression levels of related molecules in the PI3K/Akt pathway in HepG2 cell lines. Gene expression was evaluated by RT-qPCR, using the 2 −∆∆CT method. (CON, negative control treatment; GJF-L, treated with low concentrations (0.5 mg/ml) of GJF; GJF-M, treated with medium concentrations (1.0 mg/ml) of GJF; GJF-H, treated with medium concentrations (1.5 mg/ml) of GJF. Data are presented as the mean ± SD, and analysis of variance was used to determine the significance. *p < 0.05; **p < 0.01; ***p < 0.001 compared with the control group).
The wound healing assay showed that compared with the control group, the cell mobility of the high concentration group and the low concentration group showed a downward trend in the HepG2 cell line. After 24 h of treatment, the inhibitory effect of the high concentration group was significantly stronger than that of the control group (Figure 8B).
3.8 Effects of the cell cycle and apoptosis
The cell cycle can be divided into the GOG1 phase, G phase, and G2/M phase, and the ratio of the GOG1 phase is usually used to represent its inhibitory effect on proliferation. As shown in the figure below, GJF blocked HepG2 cell lines in the GOG1 phase in a concentration-dependent manner (Figure 8C). Compared with the control group, the difference between the low-dose, medium-dose, and high-dose groups and the control group was statistically significant.
As shown in Figure 8D, GJF blocked HepG2 cell lines in the G0/G1 phase in a concentration-dependent manner. In the HepG2 cell line, the apoptosis rate of the medium-dose group was higher than that of the high-dose group. Compared with the control group, the differences between the low-dose group, medium-dose group, and high-dose group and the control group were statistically significant.
3.9 Validation of related indicators in the PI3K/Akt signaling pathway
We verified the relevant molecules in the PI3K/Akt signaling pathway in the HepG2 cell line (Supplementary Table S3). According to the RT‒qPCR results (Figure 8E), we observed that GJF reduced the expression levels of PIK3R1 and Akt1 in the HepG2 cell line, and the difference was statistically significant. These results suggest that GJF may inhibit the expression level of the PI3K/Akt signaling axis in HCC. On the other hand, GJF decreased the levels of cycling-related proteins (CCND1, CDK2, CDK4, and CDK6) in HepG2 cell lines and increased the expression of the tumor suppressor gene CDKN2A. For apoptosis-related proteins, GJF increased the expression levels of Bax and Bim while decreased the expression level of Bel-2. In the HepG2 cell line, GJF decreased the expression levels of EPHA2.
4 Discussion
HCC, the most prevalent type of PLC, is a common malignant tumor threatening human health and life (Siegel et al., 2020). Clinically, the combination of traditional Chinese and Western medicine is very important in the treatment of HCC. Preliminary clinical and animal studies have explored the effects of GJF on anti-inflammatory, antitumor, liver protection, inhibition of tumor cell proliferation, and induction of tumor cell apoptosis. Therefore, this study further explored the clinical application value and possible mechanism of action of GJF along with previous research ideas.
To further explore the mechanism of GJF intervention in the HCC process, we first screened the potential active compounds of GJF acting on HCC and clarified the molecular mechanism of GJF intervention in HCC through the PI3K/Akt signaling axis using network pharmacology. A total of 63 key targets of GJF acting on HCC were obtained, and the results showed that the top five BPs involved in these genes were epithelial cell proliferation, regulation of epithelial cell proliferation, regulation of cell cycle phase transition, gland development, and embryonic organ development, respectively; In terms of CC, the targets were mainly involved in cell-substrate junction, cell-substrate adherens junction, focal adhesion, adherens junction, and extracellular matrix composition. In the category of MF, they were mainly associated with cell adhesion molecule binding, proximal promoter sequence-specific DNA binding, chromatin binding, RNA polymerase II proximal promoter sequence-specific DNA binding and laminin-binding are related. In addition, the top five pathways were the cell cycle, cellular senescence, p53 signaling pathway, PI3K/Akt signaling pathway, and progesterone-mediated oocyte maturation. Studies have shown that telomere shortening and cell cycle checkpoint inactivation are characteristic of the development of human liver cancer (Plentz et al., 2007). Cellular senescence is the process leading to end-of-life growth arrest with characteristic morphological features, mediated by telomere dependent, oncogene-induced, and ROS-induced pathways, but persistent DNA damage is the most common cause. Aging arrest is mediated by p16-and p21-dependent pathways that both lead to retinoblastoma protein activation. p53 acts as a relay between DNA damage detection and p21 activation, and aging arrest and cellular immortality are likely to significantly promote the development of HCC. Aging in oncogene-induced precancerous lesions and reversible immortality of cancer cells, including HCCs, offer new potential for tumor prevention and treatment (Ozturk et al., 2009). Yan Wang et al. showed that circulating neutrophils can predict low HCC survival and promote HCC progression through the p53 and STAT3 signaling pathways (Wang et al., 2020). Several studies have shown that progestin-mediated oocyte maturation is strongly associated with the progression of HCC, but further experimental confirmation is needed (B. Jin et al., 2015; Ma et al., 2020). Xiao-Lu Ma et al. showed that CD73 could promote the progression and metastasis of HCC by inducing RAP1-mediated P110β membrane localization to activate PI3K/Akt signaling, which was associated with poor prognosis of HCC(X. L. Ma et al., 2019). Similarly, the results showed that GJF may regulate the mRNA expression levels of relevant molecules in the PI3K/Akt signaling pathway and play a role in regulating the mRNA expression levels of proteins related to proliferation and apoptosis downstream of the pathway.
HSP90AA1 plays an important role in the proliferation, differentiation, survival, and angiogenesis of tumor cells (Jego et al., 2013). HSP90AA1 is highly expressed in a variety of malignancies, including breast, endometrial, ovarian, colon, lung, and prostate cancers (Kang et al., 2010; Y; Li et al., 2011; Soroka et al., 2012), which is consistent with our conclusion based on the TCGA database. Qiuran Xu et al. showed that Hsp90 can regulate the abundance of PKM2 through phosphorylation of Thr-328 to promote cell glycolysis and proliferation and inhibit apoptosis of HCC cells (Q. Xu et al., 2017). Xiao Xiang et al. showed that the VEGF/VEGFR2 pathway may be related to the recurrence of HCC in patients with high expression of HSP90AA1 (Xiang et al., 2018). EPHA2 overexpression is associated with tumor progression, metastasis and prognosis of HCC. Pu Yang et al. showed that EPHA2 is involved in key mediators of angiogenesis and invasion (Yang et al., 2009). Qiao Jin et al. showed that EPHA2 can negatively regulate the radiosensitivity of MHCC97H cells (Jin et al., 2015), which is consistent with our findings. We verified the application value of EPHA2 by analyzing data from the TCGA database, and the results showed that EPHA2 was differentially expressed in a variety of tumors, and its high expression was associated with poor prognosis in HCC, but the difference between the high and low expression of EPHA2 was not statistically significant. The RT‒qPCR results showed that GJF decreased the mRNA expression levels of EPHA2 in HepG2 cells, and the expression levels decreased gradually with increasing hepatocellular formula concentration. The mRNA expression level of HSP90AA1 in HepG2 cells decreased gradually with increasing GJF concentration but only in the high-dose group (p < 0.05). Our hypothesis was further verified by RT‒qPCR results. GJF may inhibit the mRNA expression levels of related molecules in the PI3K/Akt signaling pathway and regulate the mRNA expression levels of proliferation- and apoptosis-related proteins downstream of the pathway. Combined with molecular docking technology, EPHA2 may be the key target of GJF in HCC through the PI3K/Akt signaling pathway.
Molecular docking results showed that PF and PA might be the key compounds in hepatocellular formation and have anticancer activity. PF is an important compound that protects the liver, reduces cholestasis, alleviates liver fibrosis, and prevents HCC(X. Ma et al., 2020). Yang Zhou et al. showed that PF may block the expression of the Wnt/β-catenin pathway by downregulating the expression of 5-HT1D, thus inhibiting the occurrence of HCC(Zhou et al., 2020). In addition, PF is an effective antimetastatic and anti-invasive agent for inhibiting HCC invasion and metastasis (J. T. Lu et al., 2014). PA is a triterpenoid compound with anticancer activity that has anti-inflammatory, anticancer, anti-aging, and insulin-like properties and other pharmacological effects (Miao et al., 2019). It induces apoptosis in human prostate cancer (Gapter et al., 2005), lung cancer (C. Lu et al., 2018), and bladder cancer (Jeong et al., 2015) cells. In conclusion, our results suggest that PA and PF may be the key compounds of hepatocellular development that exert anticancer activity.
However, our study had certain limitations. Because TCM has the characteristics of multitarget and multipathway effects, the pathway we verified is likely to be only one way for GJF to exert its effects. In this study, we made predictions based on data analysis and a small number of in vitro experiments. GJF is a well-known prescription with beneficial properties. Therefore, further in vivo and in vitro experimental verification is still necessary.
5 Conclusion
In conclusion, this study revealed a potential pharmacological mechanism by which GJF can treat HCC at the systemic level, possibly involving the synergistic effects of multiple mechanisms, such as cell proliferation, apoptosis, cell migration, immune regulation, and inflammatory induction. Our results suggested that GJF may induce G0/G1 phase cycle arrest by inhibiting the PI3K/Akt signaling pathway and promoting the apoptosis of liver cancer cell lines. EPHA2 may be an important target of GJF in HCC, and pachymic acid may be an important active compound of GJF that exerts anticancer activity. This study can provide a scientific basis for the subsequent clinical application of GJF and the in-depth study of its mechanism. Given that the study is primarily based on data analysis, further biological experiments are crucial to verify the results.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
Part of the data in this paper are from the dissertation of co-first author QY. The data presented in the current study were previously deposited as a dissertation in the archives of Guangzhou University of Chinese Medicine, accession number 20171101040. The first author MY participated in the main writing, modification and experimental assistance of this study. XW and QY put forward the idea of the study and supervised the study. MY, QY, and YL carried out the data preparation network construction. BW and SC participated in the experiment. LL and BY analyzed the results. MY and QY wrote the manuscript. XW participated and directed the revision of the manuscript. All authors read and approved the final version of the manuscript.
Funding
This work was supported by the General Project of National Natural Science Foundation of China (Receive No. 8227152807, Grant No. 82274605).
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.
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.2023.1016967/full#supplementary-material
Supplementary Table 1 | Top 5 compounds in GJF (ranked by OB %).
Supplementary Table 2 | Top 10 KEGG Pathways in which GJF acted on the key targets of HCC (ranked by number of genes enriched).
Supplementary Table 3 | mRNA expression level of each molecule in HepG2 cell line.
References
Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 68 (6), 394–424. doi:10.3322/caac.21492
Daina, A., Michielin, O., and Zoete, V. (2019). SwissTargetPrediction: Updated data and new features for efficient prediction of protein targets of small molecules. Nucleic Acids Res. 47 (W1), W357–W364. doi:10.1093/nar/gkz382
Gao, J., Aksoy, B. A., Dogrusoz, U., Dresdner, G., Gross, B., Sumer, S. O., et al. (2013). Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci. Signal 6 (269), pl1. doi:10.1126/scisignal.2004088
Gao, S., Liu, J., Wang, M., Liu, Y., Meng, X., Zhang, T., et al. (2019). Exploring on the bioactive markers of Codonopsis Radix by correlation analysis between chemical constituents and pharmacological effects. J. Ethnopharmacol. 236, 31–41. doi:10.1016/j.jep.2019.02.032
Gapter, L., Wang, Z., Glinski, J., and Ng, K.-y. (2005). Induction of apoptosis in prostate cancer cells by pachymic acid from Poria cocos. Biochem. Biophysical Res. Commun. 332 (4), 1153–1161. doi:10.1016/j.bbrc.2005.05.044
Gotoh, J., Obata, M., Yoshie, M., Kasai, S., and Ogawa, K. (2003). Cyclin D1 over-expression correlates with beta-catenin activation, but not with H-ras mutations, and phosphorylation of Akt, GSK3 beta and ERK1/2 in mouse hepatic carcinogenesis. Carcinogenesis 24 (3), 435–442. doi:10.1093/carcin/24.3.435
Gunasekaran, G., Bekki, Y., Lourdusamy, V., and Schwartz, M. (2020). Surgical treatments of hepatobiliary cancers. Hepatology 73 (S1), 128–136. doi:10.1002/hep.31325
Hussaini, S., Chehade, R., Boldt, R. G., Raphael, J., Blanchette, P., Maleki Vareki, S., et al. (2021). Association between immune-related side effects and efficacy and benefit of immune checkpoint inhibitors - a systematic review and meta-analysis. Cancer Treat. Rev. 92, 102134. doi:10.1016/j.ctrv.2020.102134
Jego, G., Hazoumé, A., Seigneuric, R., and Garrido, C. (2013). Targeting heat shock proteins in cancer. Cancer Lett. 332 (2), 275–285. doi:10.1016/j.canlet.2010.10.014
Jeong, J.-W., Lee, W. S., Go, S.-i., Nagappan, A., Baek, J. Y., Lee, J.-D., et al. (2015). Pachymic acid induces apoptosis of EJ bladder cancer cells by DR5 up-regulation, ROS generation, modulation of bcl-2 and IAP family members. Phytotherapy Res. 29 (10), 1516–1524. doi:10.1002/ptr.5402
Jin, B., Wang, W., Du, G., Huang, G. Z., Han, L. T., Tang, Z. Y., et al. (2015). Identifying hub genes and dysregulated pathways in hepatocellular carcinoma. Eur. Rev. Med. Pharmacol. Sci. 19 (4), 592–601.
Jin, Q., Li, X., and Cao, P. (2015). EphA2 modulates radiosensitive of hepatocellular carcinoma cells via p38/mitogen-activated protein kinase-mediated signal pathways. Kaohsiung J. Med. Sci. 31 (10), 510–517. doi:10.1016/j.kjms.2015.09.001
Kang, G. H., Lee, E. J., Jang, K. T., Kim, K. M., Park, C. K., Lee, C. S., et al. (2010). Expression of HSP90 in gastrointestinal stromal tumours and mesenchymal tumours. Histopathology 56 (6), 694–701. doi:10.1111/j.1365-2559.2010.03550.x
Khan, S. A., and Lee, T. K. W. (2022). Network-pharmacology-based study on active phytochemicals and molecular mechanism of cnidium monnieri in treating hepatocellular carcinoma. Int. J. Mol. Sci. 23 (10), 5400. doi:10.3390/ijms23105400
Law, C. W., Chen, Y., Shi, W., and Smyth, G. K. (2014). voom: Precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biol. 15 (2), R29. doi:10.1186/gb-2014-15-2-r29
Lee, L., Wang, K., Li, G., Xie, Z., Wang, Y., Xu, J., et al. (2011). Liverome: A curated database of liver cancer-related gene signatures with self-contained context information. BMC Genomics 12, S3. doi:10.1186/1471-2164-12-S3-S3
Li, S., and Zhang, B. (2013). Traditional Chinese medicine network pharmacology: Theory, methodology and application. Chin. J. Nat. Med. 11 (2), 110–120. doi:10.1016/s1875-5364(13)60037-0
Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). Timer: A web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 77 (21), e108–e110. doi:10.1158/0008-5472.can-17-0307
Li, Y., Zhang, T., Schwartz, S. J., and Sun, D. (2011). Sulforaphane potentiates the efficacy of 17-allylamino 17-demethoxygeldanamycin against pancreatic cancer through enhanced abrogation of Hsp90 chaperone function. Nutr. cancer 63 (7), 1151–1159. doi:10.1080/01635581.2011.596645
Liu, H., Wang, Z. Y., Zhou, Y. C., Song, W., Ali, U., and Sze, D. M. (2021). Immunomodulation of Chinese herbal medicines on NK cell populations for cancer therapy: A systematic review. J. Ethnopharmacol. 268, 113561. doi:10.1016/j.jep.2020.113561
Liu, Q. M., Zhao, H. Y., Zhong, X. K., and Jiang, J. G. (2012). Eclipta prostrata L. Phytochemicals: Isolation, structure elucidation, and their antitumor activity. Food Chem. Toxicol. 50 (11), 4016–4022. doi:10.1016/j.fct.2012.08.007
Liu, X., Li, M., Wang, X., Dang, Z., Yu, L., Wang, X., et al. (2019). Effects of adjuvant traditional Chinese medicine therapy on long-term survival in patients with hepatocellular carcinoma. Phytomedicine 62, 152930. doi:10.1016/j.phymed.2019.152930
Lu, C., Cai, D., and Ma, J. (2018). Pachymic acid sensitizes gastric cancer cells to radiation therapy by upregulating Bax through hypoxia. Am. J. Chin. Med. 46 (4), 875–890. doi:10.1142/s0192415x18500465
Lu, J. T., He, W., Song, S. S., and Wei, W. (2014). Paeoniflorin inhibited the tumor invasion and metastasis in human hepatocellular carcinoma cells. Bratisl. Lek. Listy 115 (7), 427–433. doi:10.4149/bll_2014_084
Lu, M., Hu, Q., Zhang, Y., Zhai, Y., Zhou, Y., and Jiang, J. (2019). Comparative chemical profiling of three TCM drugs in the Paeoniaceae family by UPLC-MS/MS combined with chemometric methods. Biochem. Syst. Ecol. 83, 121–129. doi:10.1016/j.bse.2019.02.002
Luo, T. T., Lu, Y., Yan, S. K., Xiao, X., Rong, X. L., and Guo, J. (2020). Network pharmacology in research of Chinese medicine formula: Methodology, application and prospective. Chin. J. Integr. Med. 26 (1), 72–80. doi:10.1007/s11655-019-3064-0
Ma, C., Kang, W., Yu, L., Yang, Z., and Ding, T. (2020). AUNIP expression is correlated with immune infiltration and is a candidate diagnostic and prognostic biomarker for hepatocellular carcinoma and lung adenocarcinoma. Front. Oncol. 10, 590006. doi:10.3389/fonc.2020.590006
Ma, X. L., Shen, M. N., Hu, B., Wang, B. L., Yang, W. J., Lv, L. H., et al. (2019). CD73 promotes hepatocellular carcinoma progression and metastasis via activating PI3K/AKT signaling by inducing Rap1-mediated membrane localization of P110β and predicts poor prognosis. J. Hematol. Oncol. 12 (1), 37. doi:10.1186/s13045-019-0724-7
Ma, X., Zhang, W., Jiang, Y., Wen, J., Wei, S., and Zhao, Y. (2020). Paeoniflorin, a natural product with multiple targets in liver diseases-A mini review. Front. Pharmacol. 11, 531. doi:10.3389/fphar.2020.00531
Ma, Z., Fan, Y., Wu, Y., Kebebe, D., Zhang, B., Lu, P., et al. (2019). Traditional Chinese medicine-combination therapies utilizing nanotechnology-based targeted delivery systems: A new strategy for antitumor treatment. Int. J. Nanomedicine 14, 2029–2053. doi:10.2147/IJN.S197889
Miao, G., Han, J., Zhang, J., Wu, Y., and Tong, G. (2019). Targeting pyruvate kinase M2 and hexokinase II, pachymic acid impairs glucose metabolism and induces mitochondrial apoptosis. Biol. Pharm. Bull. 42 (1), 123–129. doi:10.1248/bpb.b18-00730
Ozturk, M., Arslan-Ergul, A., Bagislar, S., Senturk, S., and Yuzugullu, H. (2009). Senescence and immortality in hepatocellular carcinoma. Cancer Lett. 286 (1), 103–113. doi:10.1016/j.canlet.2008.10.048
Paolacci, S., Precone, V., Acquaviva, F., Chiurazzi, P., Fulcheri, E., Pinelli, M., et al. (2019) Genetics of lipedema: New perspectives on genetic research and molecular diagnoses. Eur. Rev. Med. Pharmacol. Sci. 23 (13), 5581–5594. doi:10.26355/eurrev_201907_18292
Petrowsky, H., Fritsch, R., Guckenberger, M., De Oliveira, M. L., Dutkowski, P., and Clavien, P. A. (2020). Modern therapeutic approaches for the treatment of malignant liver tumours. Nat. Rev. Gastroenterol. Hepatol. 17 (12), 755–772. doi:10.1038/s41575-020-0314-8
Plentz, R. R., Park, Y. N., Lechel, A., Kim, H., Nellessen, F., Langkopf, B. H., et al. (2007). Telomere shortening and inactivation of cell cycle checkpoints characterize human hepatocarcinogenesis. Hepatology 45 (4), 968–976. doi:10.1002/hep.21552
Ren, M., McGowan, E., Li, Y., Zhu, X., Lu, X., Zhu, Z., et al. (2019). Saikosaponin-d suppresses COX2 through p-STAT3/C/EBPβ signaling pathway in liver cancer: A novel mechanism of action. Front. Pharmacol. 10, 623. doi:10.3389/fphar.2019.00623
Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43 (7), e47. doi:10.1093/nar/gkv007
Siegel, R. L., Miller, K. D., Fuchs, H. E., and Jemal, A. (2021). Cancer statistics, 2021. CA Cancer J. Clin. 71 (1), 7–33. doi:10.3322/caac.21654
Siegel, R. L., Miller, K. D., and Jemal, A. (2020). Cancer statistics, 2020. CA Cancer J. Clin. 70 (1), 7–30. doi:10.3322/caac.21590
Smyth, G. K. (2004). Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat. Appl. Genet. Mol. Biol. 3, Article3. doi:10.2202/1544-6115.1027
Soroka, J., Wandinger, S. K., Mäusbacher, N., Schreiber, T., Richter, K., Daub, H., et al. (2012). Conformational switching of the molecular chaperone Hsp90 via regulated phosphorylation. Mol. Cell 45 (4), 517–528. doi:10.1016/j.molcel.2011.12.031
Su, W. H., Chao, C. C., Yeh, S. H., Chen, D. S., Chen, P. J., and Jou, Y. S. (2007). OncoDB.HCC: An integrated oncogenomic database of hepatocellular carcinoma revealed aberrant cancer target genes and loci. Nucleic Acids Res. 35, D727–D731. doi:10.1093/nar/gkl845
Wang, Y., Yao, R., Zhang, D., Chen, R., Ren, Z., and Zhang, L. (2020). Circulating neutrophils predict poor survival for HCC and promote HCC progression through p53 and STAT3 signaling pathway. J. Cancer 11 (13), 3736–3744. doi:10.7150/jca.42953
Wen-Sheng, W. (2003). ERK signaling pathway is involved in p15INK4b/p16INK4a expression and HepG2 growth inhibition triggered by TPA and Saikosaponin a. Oncogene 22 (7), 955–963. doi:10.1038/sj.onc.1206237
Xiang, X., You, X. M., and Li, L. Q. (2018). Expression of HSP90AA1/HSPA8 in hepatocellular carcinoma patients with depression. Onco Targets Ther. 11, 3013–3023. doi:10.2147/ott.s159432
Xu, J., Li, Y., Lu, J., Pan, T., Ding, N., Wang, Z., et al. (2015). The mRNA related ceRNA-ceRNA landscape and significance across 20 major cancer types. Nucleic Acids Res. 43 (17), 8169–8182. doi:10.1093/nar/gkv853
Xu, Q., Tu, J., Dou, C., Zhang, J., Yang, L., Liu, X., et al. (2017). HSP90 promotes cell glycolysis, proliferation and inhibits apoptosis by regulating PKM2 abundance via Thr-328 phosphorylation in hepatocellular carcinoma. Mol. Cancer 16 (1), 178. doi:10.1186/s12943-017-0748-y
Yang, P., Yuan, W., He, J., Wang, J., Yu, L., Jin, X., et al. (2009). Overexpression of EphA2, MMP-9, and MVD-CD34 in hepatocellular carcinoma: Implications for tumor progression and prognosis. Hepatol. Res. 39 (12), 1169–1177. doi:10.1111/j.1872-034X.2009.00563.x
Ye, L., Jia, Y., Ji, K. E., Sanders, A. J., Xue, K., Ji, J., et al. (2015). Traditional Chinese medicine in the prevention and treatment of cancer and cancer metastasis. Oncol. Lett. 10 (3), 1240–1250. doi:10.3892/ol.2015.3459
Yu, G., Wang, L. G., Han, Y., and He, Q. Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 16 (5), 284–287. doi:10.1089/omi.2011.0118
Zhou, Y., Liu, X., Gao, Y., Tan, R., Wu, Z., Zhong, Q., et al. (2020). Paeoniflorin affects hepatocellular carcinoma progression by inhibiting Wnt/β-catenin pathway through downregulation of 5-HT1D. Curr. Pharm. Biotechnol. 22, 1246–1253. doi:10.2174/1389201021666201009153808
Keywords: Ganji Fang, hepatocellular carcinoma, network pharmacology, HepG2 cells, PI3K/Akt signaling pathway, EPHA2
Citation: Yang M, Yan Q, Luo Y, Wang B, Deng S, Luo H, Ye B and Wang X (2023) Molecular mechanism of Ganji Fang in the treatment of hepatocellular carcinoma based on network pharmacology, molecular docking and experimental verification technology. Front. Pharmacol. 14:1016967. doi: 10.3389/fphar.2023.1016967
Received: 11 August 2022; Accepted: 05 January 2023;
Published: 19 January 2023.
Edited by:
Khuloud Bajbouj, University of Sharjah, United Arab EmiratesReviewed by:
Yuhe Lei, Shenzhen Hospital of Guangzhou University of Chinese Medicine, ChinaWilfred W. K. Lin, HerbMiners Informtics Limited, Hong Kong SAR, China
Shakeel Ahmad Khan, Daegu Gyeongbuk Institute of Science and Technology (DGIST), South Korea
Copyright © 2023 Yang, Yan, Luo, Wang, Deng, Luo, Ye and Wang. 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: Xiongwen Wang, YXdlbjY4MTAyOUAxNjMuY29t
†These authors have contributed equally to this work and share first authorship