Ginger compounds and their interacting targets
After the searching, filtering, and removal of the duplicates, 6 compounds (MOL002464, MOL002501, MOL002514, MOL000358, MOL000359, and MOL002467) were obtained from TCMSP database (listed in Table 1). It was noted that the representative component 6-gingerol was recruited in this study, since the DL value is relatively close to 0.18. Meanwhile, 288 targets in total interacting with these active compounds were collected, among which 254 were obtained from the Swiss TargetPrediction and 34 from TCMSP database.

Disease-associated targets.
After removal of the duplicates, there remained 1356 colon cancer-associated genes, in which 1165 genes were from the GeneCards database, 156 from the OMIM, and 35 from the Drugbank. After merging the disease and active compounds related targets, 114 candidate targets were collected for further mechanisms study of ginger on treatment of colon cancer (Fig. 2A).
PPI network analysis.
The 114 candidate targets were connected to establish an initial PPI network that included 103 nodes and 510 edges, and 4 isolated targets genes were removed (Fig. 2B). In addition, the top two protein-protein interacting clusters were constructed (Fig. 2C and 2D). The PIK3R1, MAPK3, and TP53 were as the core targets for one cluster (19 nodes and 47 edges), while another cluster (17 nodes and 42 edges) as the STAT 3 and HSP90AA1. As shown in Fig. 2, the color and size of the nodes reflected the degree value, and larger and redder nodes meant a greater “degree” value, such as PIK3CA, SRC, PIK3R1, and TP53. Subsequently, the median values of three topological indexes (degree, betweenness, and closeness) were calculated, and a total of 32 key targets, with SRC, PIK3R1, TP53, STAT3, and HSP90AA1 as the top ones, were collected for pathway enrichment analysis (Table 2).
Table 2
The 32 key targets and the counts of their interacting targets.
No. | Uniprot ID | Gene symbol | Target name | Edge |
1 | P12931 | SRC | Proto-oncogene tyrosine-protein kinase | 37 |
2 | P27986 | PIK3R1 | Phosphatidylinositol 3-kinase regulatory subunit alpha | 37 |
3 | P04637 | TP53 | Cellular tumor antigen p53 | 35 |
4 | P40763 | STAT3 | Signal transducer and activator of transcription 3 | 31 |
5 | P07900 | HSP90AA1 | Heat shock protein HSP 90-alpha | 27 |
6 | P45983 | MAPK8 | Mitogen-activated protein kinase 8 | 23 |
7 | O60674 | JAK2 | Tyrosine-protein kinase JAK2 | 20 |
8 | P00533 | EGFR | Epidermal growth factor receptor | 19 |
9 | P42338 | PIK3CB | Phosphatidylinositol 4,5-bisphosphate 3-kinase catalytic subunit beta isoform | 19 |
10 | P23458 | JAK1 | Tyrosine-protein kinase | 18 |
11 | P03372 | ESR1 | Estrogen receptor | 18 |
12 | Q05513 | PRKCZ | Protein kinase C zeta type | 18 |
13 | Q05655 | PRKCD | Protein kinase C delta type | 18 |
14 | Q04206 | RELA | Transcription factor p65 | 17 |
15 | Q05397 | PTK2 | Focal adhesion kinase 1 | 15 |
16 | P42574 | CASP3 | Caspase-3 | 15 |
17 | P04626 | ERBB2 | Receptor tyrosine-protein kinase erbB-2 | 15 |
18 | P08254 | JUN | Mitogen-activated protein kinase 8 | 14 |
19 | P52333 | JAK3 | Tyrosine-protein kinase JAK3 | 14 |
20 | Q02750 | MAP2K1 | Dual specificity mitogen-activated protein kinase kinase 1 | 14 |
21 | P06493 | CDK1 | Cyclin-dependent kinase 1 | 14 |
22 | P10275 | AR | Androgen receptor | 13 |
23 | P18031 | PTPN1 | Tyrosine-protein phosphatase non-receptor type 1 | 12 |
24 | P20248 | CCNA2 | Cyclin-A2 | 12 |
25 | P42345 | MTOR | Serine/threonine-protein kinase mTOR | 10 |
26 | P06401 | PGR | Progesterone receptor | 9 |
27 | P14780 | MMP9 | Matrix metalloproteinase-9 | 9 |
28 | P07949 | RET | Proto-oncogene tyrosine-protein kinase receptor Ret | 8 |
29 | P53350 | PLK1 | Serine/threonine-protein kinase PLK1 | 8 |
30 | P08254 | MMP3 | Stromelysin-1 | 8 |
31 | P49841 | GSK3B | Glycogen synthase kinase-3 beta | 8 |
32 | O14757 | CHEK1 | Serine/threonine-protein kinase Chk1 | 8 |
GO analysis.
As shown in Fig. 3, eventually, 20 markedly enriched BPs terms were identified, with transmembrane receptor protein tyrosine kinase GO: 0007169, regulation of cellular protein localization GO: 1903827, and positive regulation of kinase activity GO: 0045860 as the top ones. In addition, the size of the dot in bubble chart indicates the number of target genes in the corresponding function pathway, and the enrichment expresses the ratio of the number of target genes belonging to the number of all the annotated genes located in the pathway. Compared to other pathway, these three BPs terms mentioned above also enriched more counts of targets, while growth hormone receptor signaling pathway displayed a higher enrichment. For CCs, the items with significant enrichment were in glutamatergic synapse GO: 0098978, membrane raft GO: 0045121, and microtubule organizing center GO: 0005815; and for MFs in phosphotransferase activity, alcohol group receptor GO: 0016773, kinase binding GO: 0019900, and protein tyrosine kinase activity GO: 0004713.
KEGG pathway enrichment analysis.
The KEGG enrichment will be able to display how ginger acts on the signaling pathway, thereby playing a therapeutic role in colon cancer. Here, based on the 32 core targets, 10 significant signaling pathways were obtained (Fig. 4A), with pathways in cancer, hepatitis B, and estrogen signaling pathway as the top ones. As displayed in Fig. 4B, pathways in cancer was identified as a set of significantly key pathway with most targets enrichment, and bladder cancer showed a higher level of enrichment.
Drug-compound-disease-target-pathway network.
The drug-compound-disease-target-pathway network was shown in Fig. 5, which included 50 nodes (6 compounds, 32 targets, and 10 pathways) and 161 edges. The red rectangle node is the Chinese herbal medicine ginger; darker blue nodes is colon cancer, V nodes is 32 core target genes, purple ellipse nodes are 10 significant signaling pathway, blue ellipse node represents 6 effective active ingredients, while lines represent the interactions between them. According to the network analysis, each active compound acts on at least one target genes, and 6-gingerol was regarded as the most effective compound that interacts with 20 target genes. Most of genes were regulated by at least 2 active components, and the most significant pathway in cancer associated with 19 core targets. This network analysis indicated the characteristics of multiple components and multiple targets of ginger in the treatment of colon cancer.
Molecular docking result and analysis.
To validate the 8 core targets and their interacting compounds, the molecular docking were performed. The binding affinities of the testing compounds were compared to the corresponding ligands. As the LibDockScores displayed in Table 3 and Fig. 6, all the compounds had strong interactions than the prototype ligands, or similar effects to the ligands, except for targets SRC and PIK3R1. Among these compounds, 1-monolinolein (MOL002464), belong to fatty acid compounds, always displayed the highest binding affinities with the most testing targets protein, especially for the target HSP90AA1. The molecular docking results suggested that all these compounds can be considered as the active components of ginger against colon cancer, and meanwhile, TP53, HSP90AA1, MAPK8, JAK2, CASP3, and ERBB2 were speculated to be the potential targets for reaching this effect. The representative molecular docking result were exhibited in Fig. 7.
Table 3
The LibDock Scores of 8 core targets and their interacting compounds.
Target | PDB ID | Compounds | LibDock Score |
SRC | 2BDF | Ligand 1 | 141.019 |
MOL002464 | 126.807 |
MOL002501 | 108.312 |
MOL002467 | 107.09 |
MOL000358 | 113.863 |
MOL002514 | 94.4539 |
MOL000359 | 113.863 |
PIK3R1 | 4ZOP | Ligand 2 | 124.549 |
MOL002464 | 122.846 |
MOL002501 | 116.663 |
MOL002467 | 118.603 |
MOL000358 | 122.444 |
MOL002514 | 96.8916 |
MOL000359 | 120.734 |
TP53 | 5O1F | Ligand 3 | 105.7 |
MOL002464 | 135.803 |
MOL002501 | 135.516 |
MOL002467 | 122.954 |
MOL000358 | 115.344 |
MOL002514 | 112.737 |
MOL000359 | 115.344 |
HSP90AA1 | 4BQG | Ligand 4 | 89.7644 |
MOL002464 | 130.891 |
MOL002501 | 142.366 |
MOL002467 | 119.565 |
MOL000358 | 116.638 |
MOL002514 | 110.012 |
MOL000359 | 116.638 |
MAPK8 | 4E73 | Ligand 5 | 91.9872 |
MOL002464 | 98.8619 |
MOL002501 | 99.3929 |
MOL002467 | 99.8414 |
MOL000358 | 99.7914 |
MOL002514 | 92.4867 |
MOL000359 | 99.7914 |
JAK2 | 3KCK | Ligand 6 | 106.149 |
MOL002464 | 116.428 |
MOL002501 | 112.673 |
MOL002467 | 103.286 |
MOL000358 | 101.945 |
MOL002514 | 93.3826 |
MOL000359 | 102.536 |
CASP3 | 1RE1 | Ligand 7 | 85.4841 |
MOL002464 | 111.344 |
MOL002501 | 106.137 |
MOL002467 | 104.814 |
MOL000358 | 87.2257 |
MOL002514 | 85.8912 |
MOL000359 | 85.7436 |
ERBB2 | 3PPO | Ligand 8 | 103.488 |
MOL002464 | 139.158 |
MOL002501 | 149.566 |
MOL002467 | 128.035 |
MOL000358 | 136.813 |
MOL002514 | 103.498 |
MOL000359 | 136.762 |