Collection of candidate compounds of IJTand its corresponding targets
A total of 69 compounds in IJT were acquired from TCMSP, TCM-MESH, and ETCM databases after removing duplicates. 20 compounds in IJT accorded with the criteria of OB and DL at >30% and >0.18, respectively. In addition, 4 compounds were disregard owing to having no corresponding targets. Finally, 16 candidate compounds in IJT hit 280 targets altogether via the target prediction systems of TCMSP and SwissTargetPrediction databases . Thereinto, 16 compounds were constituted by 8 flavonoids, 1 triterpenoids, 4 sesquiterpenoids, 2 phytosterols, and 1 organic acid (Table 1).The detailed information of putative targets related to the candidate compounds of IJT is listed in Table S1.
Table 1. Detailed information about candidate compounds of IJT
Collection of NSCLC-related targets of IJT
The quantity of NSCLC-related therapeutic targets from the TTD and DisGeNET databases was 63 and 65, respectively. After the duplicates were removed, 110 NSCLC-related targets were ultimately obtained. The details are depicted in Table S2. A total of 24 NSCLC-related targets of IJT were obtained by overlapping dramatically of compound target and NSCLC-related targets (Table S3 and Figure 2), which might be the potential therapeutic targets of IJT against NSCLC.
C–T network
Extensive biological and pharmacological activities were exerted by TCM via multi-compounds and multi-targets. Hence, we built the C–T network to explore the sophisticated interactions based on the compounds and their related targets on NSCLC at a systems level. As displayed in Figure 3, the C–T network consisted of 40 nodes (16 compounds and 24 NSCLC-related targets of IJT ) and 60 edges. Meanwhile, the mean degree value of candidate compounds was 3.75 (Table S4). As previously reported [36], the node is regarded as a major hub if its degree is twice higher than the average degree of all other nodes. Quercetin (degree = 20) and luteolin (degree = 13) possessed degrees greater than 7.5, which became crucial active compounds for IJT. In addition, 10 compounds were corelated with equal to or more than two targets. Similarly, 50% targets hit at least 2 compounds, which demonstrated that IJT exhibited synergistic or additive effects on NSCLC.
T–T network for IJT against NSCLC
To clarify the potential mechanisms of IJT against NSCLC, we built the T–T network for IJT against NSCLC. First, 677 protein-protein interactions of NSCLC-related targets of IJT were obtained from HINT and STRING databases. Among them, 527 were identified from HINT, 168 were obtained from STRING, and 18 were overlapped. Detailed information is listed in Table S5. Furthermore, the linkages of the targets of NSCLC-related IJT, known therapeutic target, and interactional human proteins were visualized by network. A total of 476 nodes and 677 edges were formed in the T–T network (Figure 4). The hub targets measured the significance in the network screened by degree. Nodes with degrees that were twice higher than the average degree (3) of all other nodes were regarded as candidate targets. Consequently, 23 nodes that satisfied the criteria of degrees greater than 6 were encoded as candidate hubs. The detailed information is provided in Table S6.
In order to accurately screen the target of ILT in treatment of NSCLC, the hub of PPI network was constructed with 23 nodes and 167 edges (Figure 5). The major hubs were filtrated with the three topological parameters containing “degree,” “betweenness,” and “closeness.” To ensure that all of the targets of NSCLC-related IJT were contained, the major hubs should satisfy the standard, which was “degree,” “betweenness,” and “closeness” more than 15, 0.02, and 0.8, respectively (Table S7). We ultimately determined nine major hubs for the subsequent GO and KEGG pathway enrichment analysis, including cellular tumor antigen p53 (TP53), epidermal growth factor receptor (EGFR), cyclin-dependent kinase inhibitor 2A (CDKN2A), G1/S-specific cyclin-D1 (CCND1), VEGFA, MYC, MDM2, RAC-alpha serine/threonine-protein kinase (AKT1), and receptor tyrosine-protein kinase erbB-2 (ERBB2).
GOBP enrichment and KEGG pathway analysis
The GOBP and KEGG pathway enrichment analysis were used to analyze the nine major targets of IJT against NSCLC. A total of 21 GOBP enrichment terms with statistical significance at P < 0.05 were obtained, incorporating the positive regulation of protein phosphorylation, cellular response to epidermal growth factor stimulus, response to UV-A, positive regulation of epithelial cell proliferation, cell proliferation, and so on. The top five significantly enriched terms of BP (Figure 6) might be significant biological process of IJT in the treatment of NSCLC.
According to the results of KEGG pathway enrichment analysis, 9 major proteins were mapped onto 32 significant pathways (P < 0.05, Table S8). The top 10 KEGG pathways are depicted in Figure 7, which shows that IJT combined significantly with multiple cancer-related pathways in different types of cancers, incorporating bladder cancer, pathways in cancer, endometrial cancer, microRNAs in cancer, NSCLC, glioma, melanoma, chronic myeloid leukemia, prostate cancer, and central carbon metabolism in cancer. Therefore, the potential mechanisms of IJT in the treatment of NSCLC might be mainly attributed to its synergistical modulation on the relevant pathways of cancers. In addition, NSCLC enriched TP53, CCND1, ERBB2, CDKN2A, EGFR, and AKT1 and had minimum p value. Hence, NSCLC was considered the major pathway. TP53, CCND1, ERBB2, CDKN2A, EGFR, and AKT1 were ultamately reconsidered as significant targets.
Molecular docking validation
To verify the credibility of the predicted results, we implemented the docking process to verify the binding affinity of the targets with the compounds of IJT. The two compounds (quercetin and luteolin) and major NSCLC-related targets of IJT (TP53, CCND1, ERBB2, CDKN2A, EGFR, and AKT1) were validated by molecular docking. The results demonstrated that 10 of 12 pairs of C–T had binding efficiencies. The detailed information regarding the results of molecular docking is described in Table S9.