3.1 Ingredients and Targets Analysis
In our study, the TCMSP database was used to predicted the ingredients of Piwei-Peiyuan Decoction initially. We obtained 55 ingredients for Bai Zhu, 87 ingredients for Huang Qi, 104 ingredients for Xiang Fu, 220 ingredients for Gui Zhi, 85 ingredients for Bai Shao and 46 ingredients for Liu Jinu. Then we deleted the active ingredients without matched targets and selected the active ingredients from 6 herbs according to oral bioavailability (OB) and drug likeness (DL) values (OB≥30% and DL≥0.18). Eventually, 67 active ingredients such as quercetin, isorhamnetin, luteolin and kaempferol were discovered (Table.1), and they all have hepatoprotective, anti-inflammatory, antibacterial, analgesic and antitumor effects 26-29.
Table 1. Active ingredients of Piwei-Peiyuan Decoction
Medicine
|
MOL ID
|
MOL NAME
|
OB(%)
|
DL
|
Bai Zhu
|
MOL000033
|
(3S,8S,9S,10R,13R,14S,17R)-10,13-dimethyl-17-[(2R,5S)-5-propan-2-yloctan-2-yl]-2,3,4,7,8,9,11,12,14,15,16,17-dodecahydro-1H-cyclopenta[a]phenanthren-3-ol
|
36.23
|
0.78
|
Bai Zhu
|
MOL000020
|
12-senecioyl-2E,8E,10E-atractylentriol
|
62.4
|
0.22
|
Bai Zhu
|
MOL000021
|
14-acetyl-12-senecioyl-2E,8E,10E-atractylentriol
|
60.31
|
0.31
|
Bai Zhu
|
MOL000022
|
14-acetyl-12-senecioyl-2E,8Z,10E-atractylentriol
|
63.37
|
0.30
|
Bai Zhu
|
MOL000049
|
3β-acetoxyatractylone
|
54.07
|
0.22
|
Bai Zhu
|
MOL000072
|
8β-ethoxy atractylenolide Ⅲ
|
35.95
|
0.21
|
Bai Zhu
|
MOL000028
|
α-Amyrin
|
39.51
|
0.76
|
Huang Qi
Huang Qi
Huang Qi
Huang Qi
Huang Qi
Huang Qi
Huang Qi
Huang Qi
Huang Qi
Huang Qi
Huang Qi
Huang Qi
Huang Qi
Huang Qi
Huang Qi
Huang Qi
Huang Qi
Huang Qi
Huang Qi
Huang Qi
|
MOL000438
|
(3R)-3-(2-hydroxy-3,4-dimethoxyphenyl)chroman-7-ol
|
67.67
|
0.26
|
MOL000033
|
(3S,8S,9S,10R,13R,14S,17R)-10,13-dimethyl-17-[(2R,5S)-5-propan-2-yloctan-2-yl]-2,3,4,7,8,9,11,12,14,15,16,17-dodecahydro-1H-cyclopenta[a]phenanthren-3-ol
|
36.23
|
0.78
|
MOL000380
|
(6aR,11aR)-9,10-dimethoxy-6a,11a-dihydro-6H-benzofurano[3,2-c]chromen-3-ol
|
64.26
|
0.42
|
MOL000442
|
1,7-Dihydroxy-3,9-dimethoxy pterocarpene
|
39.05
|
0.48
|
MOL000371
|
3,9-di-O-methylnissolin
|
53.74
|
0.48
|
MOL000374
|
5'-hydroxyiso-muronulatol-2',5'-di-O-glucoside
|
41.72
|
0.69
|
MOL000378
|
7-O-methylisomucronulatol
|
74.69
|
0.30
|
MOL000379
|
9,10-dimethoxypterocarpan-3-O-β-D-glucoside
|
36.74
|
0.92
|
MOL000387
|
Bifendate
|
31.10
|
0.67
|
MOL000417
|
Calycosin
|
47.75
|
0.24
|
MOL000433
|
FA
|
68.96
|
0.71
|
MOL000392
|
Formononetin
|
69.67
|
0.21
|
MOL000296
|
Hederagenin
|
36.91
|
0.75
|
MOL000398
|
Isoflavanone
|
109.99
|
0.30
|
MOL000439
|
Isomucronulatol-7,2'-di-O-glucosiole
|
49.28
|
0.62
|
MOL000354
|
Isorhamnetin
|
49.60
|
0.31
|
MOL000239
|
Jaranol
|
50.83
|
0.29
|
MOL000422
|
Kaempferol
|
41.88
|
0.24
|
MOL000211
|
Mairin
|
55.38
|
0.78
|
MOL000098
|
Quercetin
|
46.43
|
0.28
|
Bai Shao
Bai Shao
Bai Shao
Bai Shao
Bai Shao
Bai Shao
Bai Shao
Bai Shao
Bai Shao
Bai Shao
Bai Shao
Bai Shao
Bai Shao
|
MOL001930
|
Benzoyl paeoniflorin
|
31.27
|
0.75
|
MOL000359
|
Sitosterol
|
36.91
|
0.75
|
MOL000358
|
Beta-sitosterol
|
36.91
|
0.75
|
MOL000422
|
Kaempferol
|
41.88
|
0.24
|
MOL001919
|
(3S,5R,8R,9R,10S,14S)-3,17-dihydroxy-4,4,8,10,14-pentamethyl-2,3,5,6,7,9-hexahydro-1H-cyclopenta[a]phenanthrene-15,16-dione
|
43.56
|
0.53
|
MOL001921
|
Lactiflorin
|
49.12
|
0.80
|
MOL001924
|
Paeoniflorin
|
53.87
|
0.79
|
MOL000492
|
(+)-Catechin
|
54.83
|
0.24
|
MOL000211
|
Mairin
|
55.38
|
0.78
|
MOL001910
|
11alpha,12alpha-epoxy-3beta-23-dihydroxy-30-norolean-20-en-28,12beta-olide
|
64.77
|
0.38
|
MOL001928
|
albiflorin_qt
|
66.64
|
0.33
|
MOL001925
|
Paeoniflorin_qt
|
68.18
|
0.40
|
MOL001918
|
Paeoniflorgenone
|
87.59
|
0.37
|
Xiang Fu
Xiang Fu
Xiang Fu
Xiang Fu
Xiang Fu
Xiang Fu
Xiang Fu
Xiang Fu
Xiang Fu
Xiang Fu
Xiang Fu
Xiang Fu
Xiang Fu
Xiang Fu
Xiang Fu
Xiang Fu
Xiang Fu
Xiang Fu
|
MOL003044
|
Chryseriol
|
35.85
|
0.27
|
MOL000354
|
Isorhamnetin
|
49.60
|
0.31
|
MOL003542
|
8-Isopentenyl-kaempferol
|
38.04
|
0.39
|
MOL000358
|
Beta-sitosterol
|
36.91
|
0.75
|
MOL000359
|
Sitosterol
|
36.91
|
0.75
|
MOL004027
|
1,4-Epoxy-16-hydroxyheneicos-1,3,12,14,18-Pentaene
|
45.10
|
0.24
|
MOL004053
|
Isodalbergin
|
35.45
|
0.20
|
MOL004058
|
Khell
|
33.19
|
0.19
|
MOL004059
|
Khellol glucoside
|
74.96
|
0.72
|
MOL010489
|
Resivit
|
30.84
|
0.27
|
MOL004068
|
Rosenonolactone
|
79.84
|
0.37
|
MOL004071
|
Hyndarin
|
73.94
|
0.64
|
MOL004074
|
Stigmasterol glucoside_qt
|
43.83
|
0.76
|
MOL004077
|
Sugeonyl acetate
|
45.08
|
0.20
|
MOL000422
|
Kaempferol
|
41.88
|
0.24
|
MOL000449
|
Stigmasterol
|
43.83
|
0.76
|
MOL000006
|
Luteolin
|
36.16
|
0.25
|
MOL000098
|
Quercetin
|
46.43
|
0.28
|
Gui Zhi
Gui Zhi
Gui Zhi
Gui Zhi
Gui Zhi
Gui Zhi
|
MOL000073
|
ent-Epicatechin
|
48.96
|
0.24
|
MOL000492
|
(+)-Catechin
|
54.83
|
0.24
|
MOL001736
|
(-)-Taxifolin
|
60.51
|
0.27
|
MOL004576
|
Taxifolin
|
57.84
|
0.27
|
MOL000359
|
Sitosterol
|
36.91
|
0.75
|
MOL011169
|
Peroxyergosterol
|
44.39
|
0.82
|
Liu Ji’nu
Liu Ji’nu
Liu Ji’nu
Liu Ji’nu
Liu Ji’nu
|
MOL001733
|
EUPATORIN
|
30.23
|
0.37
|
MOL000358
|
Beta-sitosterol
|
36.91
|
0.75
|
MOL000006
|
Luteolin
|
36.16
|
0.25
|
MOL008127
|
Ermanin
|
58.95
|
0.30
|
MOL008135
|
3,4-di-O-caffeoylquinic acid
|
49.62
|
0.69
|
3.2 Targets of Gastric Carcinoma and Piwei-Peiyuan Decoction
9877 targets of GC were obtained by Genecards database. We selected targets with hit scores greater than median 1.55. Finally, after integrating the target information in the database, 4938 targets were saved.
3.3 Construction of Interaction Network Maps
Based on the targets of Piwei-Peiyuan Decoction and GC, 181 common targets shown in the venn diagram were recognized as targets of Piwei-Peiyuan Decoction in treatment of GC (Figure 1). A PPI network was constructed based on STRING database (Figure 2), and a Targets-Pathways Interaction Network was built based on Cytoscape 3.7.1(Figure 3). The obtained targets, such as MAPK1, TP53, IL-1, Fox and Bax, are all involved in anti-inflammatory, cell proliferation-promoting, angiogenic, and anti-tumor processes30-33.
3.4 GO and KEGG Pathway Enrichment Analysis
GO enrichment analysis consists three parts: Biological Process, Molecular Function and Cellular Component. We sent 181 targets information to the STRING database for GO analysis. Finally, 2311 biological process terms, 139 cellular component terms and 259 molecular function terms were found (Figure 4). Based on the correct P value, we found the top 5 terms in cellular component were cell (174), intracellular (166), organelle (151), intracellular organelle (150), membrane-bounded organelle (149). And top 5 terms in biological processes were cellular process (178), biological regulation (174), response to stimulus (173), regulation of biological process (168), regulation of cellular process (165). Top 5 terms in molecular funtions were binding (173), protein binding (145), ion binding (111), catalytic activity (96), organic cyclic compound binding (92).
To further uncover the potential pharmacological mechanisms of Piwei-Peiyuan Decoction against GC, pathway analysis was conducted to explore the potential pathways affected by Peiwei-Peiyuan Decoction. 190 pathways were found by KEGG database, we selected top 20 pathways and draw a bubble diagram as follows (Figure 5). Combining the pathogenesis of GC, the pathways which have no association with GC were removed. Finally, 6 remarkable terms were found to be the related pathways in treatment of GC. Results demonstrated that “Pathway in cancer”, “PI3K-Akt signaling pathway”, “MAPK signaling pathway”, “Ras signaling pathway”, “IL-17 signaling pathway”, “HIF-1alpha signaling pathway” and “TNF signal pathway” were obviously enriched. According to the KEGG analysis, We found “Pathway in cancer” is the signal pathway that contains the most targets (72 targets), the network of “Pathway in cancer” is shown below (Figure 6). These signaling pathways are closely related to cell differentiation, proliferation, apoptosis and angiogenesis, most of which play a key role in the development and progression of cancer. The molecular functions and biological processes were closely related to the occurrence and development of GC, which indicated that Piwei-Peiyuan Decoction can treat GC through multiple targets and pathways. Our research shows that Piwei-Peiyuan Decoction acts as a treatment for GC mainly through the coordinated regulation of cancer-related signal pathways. For example, PI3K-Akt signaling pathway is a key pathway that regulates cell proliferation, differentiation and metastasis during the development of cancer34, and activated AKT can lead to apoptosis by participating in the regulation of cellular protein expression, then ultimately inhibiting cell proliferation35. It has been shown that the abnormal expression of P13K-Akt signaling pathway in GC tissues is closely related to the development and prognosis of tumor, P13K-AKT signal pathway plays an important role in regulating the proliferation, invasion and metastasis of tumor cells. Meanwhile, the high expression of P13K-Akt signaling pathway in GC is related to the degree of tumor differentiation, its expression level was positively correlated with the malignancy of GC36-38. MAPK signal pathway is expressed in almost all eukaryotes. It participates in cellular activities such as gene expression, cell proliferation and apoptosis, playing a key role in cellular activities39. ERK signaling pathway is a significant part of MAPK signal pathway, the gene mutations or abnormal activation in ERK signal pathway can lead to development of cancer40. Ras is a binding kinase in the upstream of ERK/MAPK signal pathway, its point mutations can lead to dysregulation of ERK/MAPK signal pathway and abnormal cellular activity, which in turn lead to migration and invasion of cancer cells and ultimately induces the development of tumor41 42. IL-17 is an inflammatory factor with great anti-inflammatory effects. It participates in the processes of cell proliferation and differentiation, immune regulation, and tumor growth. A study43 showed that the median serum IL-17 level in GC patients was significantly higher than that in controls (9.04 VS 8.07 pg/ml, p=0.01), which shows that serum IL-17 level can be used as a new potential indicator for GC diagnosis. Xu J et al44 found that IL-17-mediated downregulation of LCN2 expression inhibited the proliferation, migration, and invasion of GC cells by regulating SLPI. In addition to IL-17A, a member of the IL-17 family, several experiments have confirmed that IL-17B/IL-17 receptor B (IL-17RB) can inhibit growth and progression of tumor. A recent study45 by Bastid J et al showed that IL-17B promotes tumor progression by promoting the secretion of chemokines and cytokines and thereby dramatically altering the tumor microenvironment. Furthermore, n GC patients, higher serum IL-17B levels have been shown to be strongly associated with poor prognostic outcomes. Therefore, it is clear that IL-17 levels are closely associated with the development of GC. However, few pharmacological findings have been reported for the treatment of GC with Piwei-Peiyuan Decoction. Therefore, more experiments should be conducted to validate our findings in the future.
3.5 Results of Molecular Docking Experiment
Molecular docking was used for verification of interaction between ingredient and its target gene. The results of molecular docking scores of the ingredients and target genes we selected in this study were shown in Table.2, and the optimal schematic of molecular docking was shown in Figure.7.
From the results, the lowest binding free energy of caspase-3 and kaempferol is -6.1 kcal/mol. There is one hydrogen bond between amino acid THR62 and kaempferol. The lowest binding free energy of caspase-3 and luteolin is -6.1 kcal/mol. There are hydrophobic interactions between amino acid THR62, SER251 and luteolin. In addition, the lowest binding free energy of MAPK1 and luteolin is -7.9 kcal/mol. There are hydrophobic interactions between amino acid LYS54, GLN105, MET108, ASP167 and luteolin. The lowest binding free energy of caspase-3 and quercetin is -6.0 kcal/mol. There are two hydrophobic bond forces between amino acid THR62 and quercetin. Moreover, The lowest binding free energy of MAPK1 and quercetin is -8.1 kcal/mol. There are two hydrophobic bond forces between amino acid ASP106, MET108 and quercetin. All compounds showed a compact binding pattern to the protein active pocket, and these interactions enabled the proteins to form stable complexes with all compounds. In this study, molecular docking experiment gives explanation for the way of protein-compound interactions and lays the theoretical foundation for further studies of Piwei-Periyuan Decoction in treatment of GC.
Table.2 Scores of Molecular Docking
Receptor_Name
|
Ligand_Name
|
Scores (kcal/mol)
|
CASP3
|
kaempferol.pdbqt
|
-6.1
|
CASP3
|
luteolin.pdbqt
|
-6.1
|
CASP3
|
quercetin.pdbqt
|
-6.0
|
MAPK1
|
luteolin.pdbqt
|
-7.9
|
MAPK1
|
quercetin.pdbqt
|
-8.1
|
Note: A higher absolute value of scores means that the binding of receptor protein and active ingredient is more stable.