Journal of Oncology

Journal of Oncology / 2021 / Article

Research Article | Open Access

Volume 2021 |Article ID 8862821 | https://doi.org/10.1155/2021/8862821

Feng Xiang, Linman Li, Jieling Lin, Shasha Li, Guiyuan Peng, "Network Pharmacology and Bioinformatics Methods Reveal the Mechanism of Zao-Jiao-Ci in the Treatment of LSCC", Journal of Oncology, vol. 2021, Article ID 8862821, 15 pages, 2021. https://doi.org/10.1155/2021/8862821

Network Pharmacology and Bioinformatics Methods Reveal the Mechanism of Zao-Jiao-Ci in the Treatment of LSCC

Academic Editor: Dan Zhao
Received10 Sep 2020
Accepted08 Jun 2021
Published28 Jun 2021

Abstract

Objective. Zao-Jiao-Ci (ZJC), a traditional Chinese medicine, is considered as a promising candidate to treat laryngeal squamous cell carcinoma (LSCC). However, the underlying molecular mechanism remains unclear. Methods. Gene expression profiles of GSE36668 were available from the GEO database, and differentially expressed genes (DEGs) of LSCC were obtained by R package; subsequently, enrichment analysis on KEGG and GO of DEGs was performed. The active ingredients of ZJC were screened from the TCMSP database, and the matched candidate targets were obtained by PharmMapper. Furthermore, we constructed protein-protein interaction (PPI) networks of DEGs and candidate targets, respectively, and we screened the core network from the merged network through combining the two PPI networks using Cytoscape 3.7.2. The key targets derived from the core network were analyzed to find out the associated KEGG signal enrichment pathway. By the GEPIA online website, Kaplan–Meier analysis was used to complete the overall survival and disease-free survival of the selected genes in the core module. Results. We identified 96 candidate targets of ZJC and 86 DEGs of LSCC, the latter including 50 upregulated genes and 36 downregulated genes. DEGs were obviously enriched in the following biological functions: extracellular structure organization, the extracellular matrix organization, and endodermal cell differentiation. The 60 key targets from the core network were enriched in the signal pathways including transcriptional misregulation cancer, cell cycle, and so on. We found that LSCC patients with high expression of HIST1H3J, HIST1H3F, and ITGA4 had worse overall survival, while higher expression of NTRK1, COPS5, HIST1H3A, and HIST1H3G had significantly worse disease-free survival. Conclusion. It suggested that the interaction between ZJC and LSCC was related to the signal pathways of transcriptional misregulation cancer and cell cycle, revealing that it may be the mechanism of ZJC in the treatment of LSCC.

1. Introduction

Laryngeal squamous cell carcinoma (LSCC) is the most common malignancy of the larynx, and its clinical manifestations are hoarseness, stridor, dyspnea, and even dysphagia [1, 2]. Disappointingly, despite various technologies such as surgery, laser therapy, and chemoradiation have advanced recently, and the survival rate has not improved because of a high rate of recurrence and metastasis [3, 4]. Therefore, in order to improve survival rates of the patients, there is an urgent need for effective treatment.

An increasing number of studies confirmed that traditional Chinese medicine (TCM) including multiple ingredients and targets play a critical role in the treatment of cancer. Zao-Jiao-Ci (ZJC), also known as Gleditsia sinensis, is a traditional Chinese medicine with a variety of bioactivities, especially antitumor activity, which has been widely used in clinic [5]. It was investigated that the ethanol extract of Gleditsia sinensis (EEGS) could suppress the growth of human colon cancer HCT116 cells in vitro and in vivo [6]. The extract of Gleditsia sinensis fruit performed inhibitory effects on esophageal squamous cell carcinoma (ESCC) cells, breast cancer MCF-7 cells, hepatoblastoma HepG2 cells, and so on [7, 8]. However, there is no study to investigate the anticancer effect of Gleditsia sinensis on LSCC, and the mechanism remains unclear.

Network pharmacology has exhibited specific utility in analyzing multicomponent and multitarget, consistent with the therapy hypothesis of complex diseases. By constructing a multilevel, multifaceted network model comprised of components, targets, pathways, and diseases, we can investigate TCM in the treatment of disease involved in the regulation of a variety of signaling pathways, key targets taxa, and biological process analysis, aiming to reveal the mechanism from the molecular level [9].

In this study, we used network pharmacology to investigate whether ZJC exerts anticancer effects on LSCC based on the GEO microarray dataset. And through the pathway enrichment analysis of the interaction targets between differentially expressed gene (DEGs) of LSCC and key node targets of ZJC, we further predicted the therapeutic mechanism of ZJC on LSCC. To our knowledge, this study is the first to explore the efficacy and mechanism of ZJC on LSCC, providing theoretical support and directions for further basic research.

2. Methods

2.1. Active Ingredients Screening and Targets Prediction for ZJC

Through the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP), all components of ZJC could be found by searching the term “Zao-Jiao-Ci.” We set oral bioavailability (OB) > 30% and drug-likeness (DL) > 0.18 as screening conditions supported by the published literatures to obtain the final active ingredients [10, 11]. PharmMapper server is the first webserver for potential drug targets identification through large-scale reverse pharmacophore mapping strategy [12]. The MOL structure of active ingredients provided by TCMSP was input into PharmMapper server (http://lilab.ecust.edu.cn/pharmmapper/) to get the targets of the pharmacophore model. The first 15 targets sorted by the fit score were seemed as candidate targets of ZJC.

2.2. Active Ingredient-Target PPI Network Construction

To explore the association between ingredients and targets, we established an interaction network. Cytoscape 3.7.2, one of the most favorite open-source software tools, provides visually biomedical interaction networks composed of protein, gene, and other types of interactions. It was used to develop an active ingredient-target PPI network to visualize the relationship between the active ingredients and their targets of ZJC.

2.3. GEO Data Collection and DEGs Identification

The original data series GSE84957 was downloaded from the Gene Expression Omnibus (GEO) microarray dataset, which contained gene expression profiles of 18 tissue samples (9 LSCC tumor tissues and 9 normal tissues). The R language was used to process the original data sets, and the RMA algorithm of Affy software package was used to perform background correction and quartile standardization of the expression matrix. The gene ID, the gene probe name of the expression matrix, was replaced by the gene symbol provided by the GPL17843 Agilent-042818 Human lncRNA Microarray 8_24_v2 platform, and the average value of multiple probes for the same gene was used for analysis. Limma package was used to identify the significant differentially expressed genes (DEGs) according to , |log2 (FC)| > 3.

The screened DEGs were mapped into a volcano map using the R language heatmap package for intuitive vision; finally, the clusterProfiler package was used to carry out GO enrichment analysis and KEGG pathway enrichment analysis for DEGs.

2.4. PPI Network Construction

BisoGenet plugin, comprising of six available PPI databases (the Biological General Repository for Interaction Datasets (BioGRID), Biomolecular Interaction Network Database (BIND), Molecular Interaction Database (MINT), Human Protein Reference Database (HPRD), and Database of Interacting Proteins (DIP)), was used to build the PPI network for DEGs and candidate target genes, respectively [13]. Then, the merged network was conducted for the two PPI networks. We filtered the output nodes with degrees of freedom greater than 2 times the median of all nodes according to the indicators of degree and betweenness centrality. Then, a core PPI network was constructed using CytoNCA, a Cytoscape plugin. The ClueGO plugin was used for the KEGG signaling pathway enrichment analysis. was taken as the inclusion standard for pathway items. The results of enrichment analysis were presented in the form of the pie chart and nodes.

2.5. Cluster of the Core PPI Network

The MCODE plugin in Cytoscape software was used to screen the highly clustered important modules in the core PPI network. We set the parameters as degree cutoff = 2 and κ-core = 2 and conducted KEGG signaling pathway enrichment analysis for the most significantly clustered modules.

2.6. Gene Expression Data of the Core Cluster for LSCC

The correlation between survival rates of LSCC patients (disease-free survival rate and overall survival rate) and the gene expression levels (NTRK1, COPS5, HIST1H3A, HIST1H3G, HIST1H3J, HIST1H3J, HIST1H3F, and ITGA4) were calculated using GEPIA online database (http://GEPIA.cancer pku.cn/) [14].

3. Results

3.1. Active Ingredients and Targets of ZJC

Searching for all the reported components in the TCMSP database, 30 active ingredients of ZJC were collected. Sequentially, only 11 active ingredients were retained which conformed to OB > 30% and DL > 0.18, such as fisetin, fustin, flavanone, and kaempferol (Table 1). Then, 96 candidate targets of the above 11 active components were obtained after the duplicate targets were excluded.


No.ComponentOB (%)DL

1Fisetin52.600.24
2Fustin50.910.24
3(-)-Taxifolin60.510.27
4Flavanone41.350.24
5Beta-sitosterol36.910.75
6Sitosterol36.910.75
7Kaempferol41.880.24
8Stigmasterol43.830.76
9Stigmast-4-ene-3,6-dione39.120.79
10Ent-epicatechin48.960.24
11Quercetin46.430.28

3.2. Active Ingredients-Targets PPI Network Construction

A PPI network of the active components and relevant targets, containing 107 nodes and 165 edges, was constructed by the network graphing tool Cytoscape 3.7.2. The 11 active ingredients could be connected with multiple targets, respectively, and each target also could be connected with multiple active ingredients, which directly demonstrated the relationship between active ingredients and targets of ZJC (Figure 1).

3.3. LSCC Differentially Expressed Genes (DEGs)

By analyzing the gene chip of GSE84957, a total of 81 genes with significant different expression of the LSCC tissues compared with adjacent nonneoplastic tissues were obtained, among which 50 genes were upregulated and 31 genes were downregulated in tumor tissues (Table 2; Figure 2).


GeneSamplelogFCAveExprt valueAdj. valueB

CST14.2555180227.41775987814.097761571.49E − 113.71E − 0715.51198443
XLOC_0044263.0442282113.37072678311.592583324.32E − 103.61E − 0612.75595229
MMP113.9033253336.98687714411.585730974.36E − 103.61E − 0612.74750381
GPRIN13.8577450445.65357341111.295657436.67E − 104.15E − 0612.38486185
COL7A13.1917376449.25824273310.780420021.45E − 097.22E − 0611.71597451
FAM3D−5.7925450677.2950313−10.097924.24E − 091.51E − 0510.77877059
LRP123.5808215564.0053048679.3161606821.54E − 084.26E − 059.629111691
CTHRC14.7408064337.1280871729.1964674281.89E − 084.30E − 059.445581444
LOC1005060274.1432875784.5686711898.9993474882.66E − 084.97E − 059.138857621
TJP3−3.3646346337.834813517−8.9040232113.14E − 084.97E − 058.9885142
IGFBP33.60413816710.552446928.8937423133.20E − 084.97E − 058.972220302
ARSI3.1526468673.9869231898.2984788129.26E − 080.0001213138.002153147
TMEM1583.7181060334.6970891178.2228941871.06E − 070.0001323537.875178695
PLAUR3.0274720337.8246327618.0514980151.46E − 070.0001651327.584034636
MSR13.9185479224.4979709067.7056722.80E − 070.0002486826.98290287
TGFBI3.4155966339.2350727947.5342420433.89E − 070.0003119416.678084284
SPP14.8530972227.3187210117.2962833876.17E − 070.0004150166.247440183
CYP4B1−3.9384062895.174530178−7.1454131338.31E − 070.0004927165.96986977
SH3BGRL2−4.0715280896.438326422−7.129994368.57E − 070.0004955825.941304558
LOC1006531493.2681193565.4438078117.0780460149.50E − 070.0005250055.844794163
MMP14.8624599226.7830569947.0022326911.10E − 060.0005906015.70320204
TM4SF194.0527956785.816474656.874358071.43E − 060.0006965695.462380919
COL5A23.5043230787.5844092836.7590328971.80E − 060.0008157425.243053839
HOXD113.5509908565.0875631726.546721592.78E − 060.0011173494.834020589
CXCL114.7623212114.6527693946.5197582652.94E − 060.001151714.781590785
LOC1006528323.2538606566.4650064396.5025262463.05E − 060.001151714.748026804
CXCL104.5191928445.6402736446.4987693643.07E − 060.001151714.740703422
XLOC_006053−3.6358996339.156396294−6.4757564063.22E − 060.001151714.69579809
KRT174.40471211111.515995616.4642350693.30E − 060.001151714.673287016
PTHLH4.0988200337.9496798946.3412124654.26E − 060.0013151844.431701824
CXCL12−3.4765157678.330393117−6.240929915.26E − 060.0014225964.233143024
XLOC_l2_0060214.37245925611.502697826.2299638095.38E − 060.0014247814.211342414
COL8A13.2111216676.0609220786.1774309726.01E − 060.0015264384.106669064
CLCA4−5.9193416568.217143061−6.1260823616.70E − 060.0016612544.003977039
HMGA23.7063291223.7034292616.1170339716.83E − 060.001666153.985842669
WISP13.0219815565.5954507896.0392662448.06E − 060.0018387783.829512486
NRG13.4305063116.6157730676.0134031378.51E − 060.0019077963.777336092
ANKRD20A9P−3.0271043116.964110233−5.9912595778.92E − 060.0019473363.732590438
GCNT3−3.4833602225.808183467−5.9680767979.38E − 060.0019654853.685672938
MYOC−3.2909366443.378123278−5.9642616629.45E − 060.0019654853.677944815
ODZ23.0308502117.3428003395.9556987689.63E − 060.0019654853.660592194
DNAPTP34.1625125898.3550088175.9480575039.79E − 060.0019792243.645098841
FUT3−3.4528091337.136305−5.8322216461.25E − 050.002213453.409270286
CRNN−6.5247332789.72241585−5.8219412971.28E − 050.0022270343.388254529
CFD−3.5200374569.369048272−5.8005101251.34E − 050.0022868883.344398953
COL1A23.09644633311.786843065.8004025321.34E − 050.0022868883.34417863
CXCL93.0673486226.6156631785.7597094151.47E − 050.0023695493.260741185
PDPN3.17027573310.04381885.7092849881.64E − 050.002529113.157054467
INHBA3.0998870785.8097533175.6560436611.84E − 050.0026404633.04722486
MMP124.1625211673.7365494175.6388645841.91E − 050.0026875473.011710913
SCARA5−3.8610034677.106386856−5.6348062531.92E − 050.0026875473.003315853
COL4A13.29399372211.038862925.6234307391.97E − 050.0027092322.979773593
COL5A13.2051796228.6743017445.6044076532.05E − 050.002792952.940368581
ANKRD20A5P−3.1188266115.357065783−5.5742108892.19E − 050.0028731172.877727019
TNXB−3.3901475448.386793894−5.464372672.79E − 050.0033455262.648949986
CCDC25−3.3405721337.218539933−5.4610950792.81E − 050.0033455262.6421014
FN13.58053544411.632547565.4279580083.02E − 050.0035259592.572791248
MAL−6.3603311338.9132381−5.3804985083.36E − 050.0037267742.473306046
FBN23.7020987674.6470011285.3175403913.86E − 050.0039636982.340945176
KRT4−5.78358505612.44990214−5.275014984.24E − 050.0043045412.251297748
MSC3.0214818115.7789504945.2730635354.25E − 050.0043045412.247179302
MYZAP−3.9666216.975986556−5.2189243654.80E − 050.004537582.132761995
CHI3L13.4477558568.4549251835.2127064094.86E − 050.0045659362.119601631
XLOC_008370−4.0868010565.071876417−5.1957743855.05E − 050.0046186852.083744902
AMY1C−3.0661629224.171807861−5.1811246825.22E − 050.0047026732.052698071
FAM107 A−3.3669290446.613174211−5.1531774925.55E − 050.0048294061.993410827
CA94.0047194894.8324982675.1057409886.17E − 050.0051995051.89260426
SFI1−3.1726544897.898643289−5.0389752457.16E − 050.0056576811.75036129
KRT24−4.3074977566.663276378−4.9034972099.71E − 050.0068028321.460515802
FAM3B−3.0028641787.472917967−4.8510356310.0001092360.0073850341.347873623
PSCA−4.4241612676.177034822−4.8366972730.0001128240.0075253051.317050516
TREM13.1321565675.3291403064.8024127530.0001218940.0077559971.243287728
MMP74.00776958.4338795834.7974473820.0001232670.0078085671.232597743
XLOC_l2_007931−3.33594024410.48608082−4.7891591730.0001255950.0078751821.214750103
KRT6B4.0766163229.7039405064.7731840490.0001302090.0079789951.180336064
SERPINE13.50385767.44660574.7529196260.000136310.0081520531.136656679
ABCA8−3.29184085.542601178−4.7498403020.0001372620.0081697171.130016856
SNX31−3.2943623333.714450967−4.7270654350.0001445170.0083636251.080888859
KRT163.95689196710.740062524.711194290.0001498010.0085675711.046633089
RSAD23.4109529115.31107054.6673431360.0001654420.0090861520.951904596
CEACAM5−3.9827577119.060092478−4.6492017320.0001723870.0093438390.912681384

3.4. GO Enrichment Analysis and KEGG Pathway Analysis for DEGs

GO enrichment analysis was used to explore the molecular mechanism of DEGs. The results were given as follows: (i) in the BP category, DEGs were mostly enriched in the extracellular structure organization, the extracellular matrix organization, endodermal cell differentiation, endoderm formation, and endoderm development; (ii) in the category of CC, DEGs were mainly enriched in the extracellular matrix, collagen-containing extracellular matrix, endoplasmic reticulum lumen, collagen trimer, and extracellular matrix component; (iii) in the MF category, extracellular matrix structural constituent, cytokine activity, and receptor ligand activity were selected for main MF. The results of KEGG pathway analysis showed that ECM-receptor interaction, protein digestion and absorption, focal adhesion, Staphylococcus aureus infection, and viral protein interaction with cytokine and cytokine receptor were the major pathways involved in DEGs (Figure 3).

3.5. PPI Network Construction and Key Targets Screening

A PPI network based on the targets of ZJC active ingredients was constructed. It showed that ZJC had direct or indirect correlation with the 1572 targets, and there were 29,098 interconnections between these targets. At the same time, the PPI network was mapped for DEGs, and 2,262 targets were directly or indirectly related to LSCC, with 50,181 interconnections between these targets. Then, the intersections of the two PPI networks were used to construct a merged network with 510 nodes and 10950 edges (Figures 4(a)4(c)). Furthermore, we analyzed the topological properties of the nodes in the merged network of the protein interactions to find the key nodes. Finally, 60 key nodes were identified through the network topology analysis (Figure 4(d) and Table 3).


Name (the key targets)DegreeBetweennessBetweenness centralityClosenessCloseness centralityTopological coefficient

YWHAZ1553001.3775710.011839190.5853658540.585365850.14373835
YWHAG87816.39124450.003220330.5373134330.537313430.16424877
MCM51011080.1330080.004260680.5425188370.542518840.14457089
STAU11141407.4858810.005551950.5502183410.550218340.13468492
ITGA42087021.251650.027695930.6168910650.616891060.10712371
APP1255666.7773490.022353090.5575221240.557522120.10929128
CUL32117900.5389770.031164360.6268656720.626865670.1083315
CUL21212193.4908540.008652410.5569060770.556906080.1311238
CUL4A841104.3888690.004356360.5190525230.519052520.14961657
CUL11582737.2700540.01079740.5799769850.579976990.12774443
COPS51734596.0813680.018129640.5929411760.592941180.11812865
YWHAQ1162340.926770.009233990.55323820.55323820.13955671
EEF1A11112303.8651010.00908780.5538461540.553846150.15726594
OBSL11512614.0517490.010311350.5720771850.572077190.11367624
TARDBP99902.42125420.003559680.5401929260.540192930.14812918
HSPA81021110.4840510.00438040.5490196080.549019610.1617004
HSPA51242786.828780.010992890.561872910.561872910.15308876
HSP90AA11153058.2417450.01206350.5581395350.558139530.13921995
EIF4A394923.52523690.003642930.5350318470.535031850.14405585
HSP90AB11112024.4665080.007985680.5550660790.555066080.14751182
MYC881227.8986150.004843550.5367412140.536741210.11599255
CCDC81452628.8101610.010369570.5688487580.568848760.11392955
HDAC1841756.6798560.006929380.5255474450.525547450.11154264
RPA11061443.9100920.005695630.5448648650.544864860.13668578
RPA2981508.9057010.005952010.5413533830.541353380.13008096
HNRNPA11301803.5675030.007114330.5656565660.565656570.14709052
EGFR1434724.0666420.018634490.5727272730.572727270.1103327
HIST1H3F86437.52732050.001725860.5294117650.529411760.14705882
UBC1333538.6552040.013958530.5669291340.566929130.12972445
CUL71744304.2315670.016978410.5943396230.594339620.10695547
HNRNPK90914.37291760.003606820.5356004250.535600430.17748918
TUBB851017.8854620.004015140.5378868730.537886870.16372149
HNRNPU1062391.6987180.009434260.5478260870.547826090.15722622
CDK21725391.3868680.021266790.5929411760.592941180.11382953
EWSR1851514.3168620.005973350.5322069690.532206970.15183716
CDC5L1132162.9549890.008531960.5508196720.550819670.13024984
NEDD886533.45167890.002104250.5283018870.528301890.14769424
HIST1H3A86437.52732050.001725860.5294117650.529411760.14705882
HIST1H3D86437.52732050.001725860.5294117650.529411760.14705882
TRAF61002757.0459310.010875410.5384615380.538461540.09652163
HIST1H3C86437.52732050.001725860.5294117650.529411760.14705882
NPM11563506.4629250.013831550.584009270.584009270.14115178
HIST1H3E86437.52732050.001725860.5294117650.529411760.14705882
HIST1H3I86437.52732050.001725860.5294117650.529411760.14705882
HIST1H3G86437.52732050.001725860.5294117650.529411760.14705882
FUS1091065.0238990.004201080.5514223190.551422320.157963
HIST1H3J86437.52732050.001725860.5294117650.529411760.14705882
HIST1H3H86437.52732050.001725860.5294117650.529411760.14705882
HIST1H3B86437.52732050.001725860.5294117650.529411760.14705882
NTRK127016394.028250.064667660.6737967910.673796790.09989618
FLNA891820.0584480.007179380.5373134330.537313430.15216314
FN130727515.326650.108536580.7088607590.708860760.09048518
SIRT71282444.8010550.009643730.55323820.55323820.11177503
GRB21364533.8476830.017884150.5701357470.570135750.12538522
FBXO61292240.283710.008836990.5581395350.558139530.11844315
VCP1031651.5245160.006514580.5448648650.544864860.14998166
CAND11522355.4295270.00929120.5786452350.578645240.12845528
XPO11062782.0813880.010974160.5431034480.543103450.12465874
MCM22076633.5227070.02616650.6214549940.621454990.11887293
YWHAE931174.4769770.004632830.5373134330.537313430.16260551

3.6. KEGG Pathway Analysis and Main Module of the Core PPI Network

The KEGG signaling pathways analysis suggested that 60 key targets were mainly enriched in cell cycle, central carbon metabolism in cancer, and DNA replication, indicating the mechanisms of ZJC in the treatment of LSCC. The other signaling pathways included prostate cancer, protein processing in endoplasmic reticulum, spliceosome, transcriptional misregulation in cancer, and ubiquitin-mediated proteolysis (Figure 5). Through the MCODE plugin, two main modules of the core PPI network were obtained, one of which was functionally enriched in alcoholism, transcriptional misregulation in cancer, and systemic lupus erythematosus (Figure 6).

3.7. LSCC Survival Analysis

To demonstrate the relationship between key genes and LSCC, we analyzed the genes in core module through the GEPIA online database and Kaplan–Meier curve. We found that LSCC patients with high expression of HIST1H3J, HIST1H3F, and ITGA4 had worse overall survival, while LSCC patients with high expression of NTRK1, COPS5, HIST1H3A, and HIST1H3G had significantly worse disease-free survival (Figure 7).

4. Discussion

Based on the network pharmacology analysis of drug and disease target, collateral relationship can effectively reveal the mechanism of ZJC in the treatment of LSCC. Here, we found 96 candidate targets of ZJC and 81 DEGs of LSCC. Then, we constructed the PPI network for them separately. The huge genes involved in the interacted PPI network were analyzed to derive the possible mechanisms of anti-LSCC of ZJC, including transcriptional misregulation cancer, alcoholism, and cell cycle.

In our study, we identified 11 active ingredients of ZJC, which synergistically regulated 96 candidate targets. A large number of published literatures showed that the 11 active ingredients had anticancer activities, respectively. As reported, fisetin could inhibit the proliferation and migration of human laryngeal cancer via ERK1/2 and AKT/NF-KB/mTOR signaling pathways and induce apoptosis in human lung cancer through the MAPK signaling pathway [15, 16]. It was also revealed that kaempferol and quercetin were potential to inhibit cell migration and invasion in human head and neck squamous cell carcinoma [17, 18]. Li et al. emphasized that taxifolin may arrest aggressive breast cancer by promoting the MET progress through decreasing the expression of β-catenin [19]. Additionally, the inhibitory potency of flavanone on human breast cancer and gastric cancer has been reported previously [20, 21]. To our knowledge, no previous studies have explored the synergistic effect of the 11 active ingredients deriving from ZJC in suppressing LSCC development.

To investigate the possible mechanism of anti-LSCC of ZJC at a system level, we applied GlueGO to complete KEGG enrichment signaling pathway analysis, through analyzing the huge targets of the core PPI network in tight corresponding to LSCC and ZJC. We identified 11 items, in particular, transcriptional misregulation in cancer, alcoholism, cell cycle, and central carbon metabolism in cancer (all ). It is apparent that both signal pathways of transcriptional misregulation in cancer and central carbon metabolism in cancer were closely associated with cancer [22, 23]. Sequentially, transcriptional misregulation in cancer was the most significant pathway following ZJC acting on LSCC (). As reported, cancer is more likely to occur in the mucous membrane in direct contact with alcohol; therefore, an intermediate increase in the risk of laryngeal cancer was found among alcoholics [24]. Aberrant cell cycle results in uncontrolled proliferation of cells, which is the common nature of cancer [25]. Zhou et al. demonstrated that Erchen decoction plus Huiyanzhuyu decoction was promising medicine in treatment of LSCC through inhibiting the cell cycle and inducing apoptosis of LSCC cells [26]. Protein processing in endoplasmic reticulum (ER) is crucial for the pathogenesis of cancer, with severe ER stress closely related to the development and invasion of cancer [27, 28]. These findings were consistent with the network pharmacology analysis.

5. Conclusion

Our study revealed that the anti-LSCC mechanism of ZJC was closely connected to transcriptional misregulation cancer, alcoholism, and cell cycle signaling pathway, which provided an important basis for further discussion of the follow-up experiment al design, making the experimental research more reasonable and more instructive.

Data Availability

The original data series GSE84957 used to support the findings of this study is downloaded from the Gene Expression Omnibus (GEO) microarray dataset.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Feng Xiang conceived the study and drafted the manuscript. Guiyuan Peng analyzed and interpreted the data and reviewed the manuscript. Shasha Li instructed the research and reviewed the manuscript. Linman Li and Jieling Lin contributed to conception and design of the study. All authors revised, read, and approved the manuscript.

References

  1. H. R. J. Pukander and J. Pukander, “Symptoms of laryngeal carcinoma and their prognostic significance,” Acta Oncologica, vol. 39, no. 2, pp. 213–216, 2009. View at: Publisher Site | Google Scholar
  2. L. D. R. Thompson and D. R. Lester, “Laryngeal dysplasia, squamous cell carcinoma, and variants,” Surgical Pathology Clinics, vol. 10, no. 1, pp. 15–33, 2017. View at: Publisher Site | Google Scholar
  3. K. D. Hunter, E. K. Parkinson, and P. R. Harrison, “Profiling early head and neck cancer,” Nature Reviews Cancer, vol. 5, no. 2, pp. 127–135, 2005. View at: Publisher Site | Google Scholar
  4. J. Zhu, S. Fedewa, and A. Y. Chen, “The impact of comorbidity on treatment (chemoradiation and laryngectomy) of advanced, nondistant metastatic laryngeal cancer,” Archives of Otolaryngology-Head & Neck Surgery, vol. 138, no. 12, pp. 1120–1128, 2012. View at: Publisher Site | Google Scholar
  5. J.-P. Zhang, X.-H. Tian, Y.-X. Yang et al., “Gleditsia species: an ethnomedical, phytochemical and pharmacological review,” Journal of Ethnopharmacology, vol. 178, no. 3, pp. 155–171, 2016. View at: Publisher Site | Google Scholar
  6. S. J. Moon, Y. H. Cho, H. Kim et al., “Inhibitory effects of the ethanol extract of Gleditsia sinensis thorns on human colon cancer HCT116 cells in vitro and in vivo,” Oncology Reports, vol. 22, no. 06, 2009. View at: Publisher Site | Google Scholar
  7. L. M. C. Chow, J. C. O. Tang, I. T. N. Teo et al., “Antiproliferative activity of the extract of Gleditsia sinensis fruit on human solid tumour cell lines,” Chemotherapy, vol. 48, no. 6, pp. 303–308, 2002. View at: Publisher Site | Google Scholar
  8. K. C. Tang, K. Y. Lam, S. Law et al., “The inhibitory effect of Gleditsia sinensis on cyclooxygenase-2 expression in human esophageal squamous cell carcinoma,” International Journal of Molecular Medicine, vol. 23, no. 1, pp. 121–129, 1998. View at: Publisher Site | Google Scholar
  9. X. Liang, H. Li, and S. Li, “A novel network pharmacology approach to analyse traditional herbal formulae: the Liu-Wei-Di-Huang pill as a case study,” Molecular BioSystems, vol. 10, no. 5, pp. 1014–1022, 2014. View at: Publisher Site | Google Scholar
  10. N. Wang, Y. Zheng, J. Gu et al., “Network-pharmacology-based validation of TAMS/CXCL-1 as key mediator of XIAOPI formula preventing breast cancer development and metastasis,” Scientific Reports, vol. 7, no. 1, p. 14513, 2017. View at: Publisher Site | Google Scholar
  11. A. Y. Lee, W. Park, T.-W. Kang, M. H. Cha, and J. M. Chun, “Network pharmacology-based prediction of active compounds and molecular targets in Yijin-tang acting on hyperlipidaemia and atherosclerosis,” Journal of Ethnopharmacology, vol. 221, no. 15, pp. 151–159, 2018. View at: Publisher Site | Google Scholar
  12. X. Liu, S. Ouyang, B. Yu et al., “PharmMapper server: a web server for potential drug target identification using pharmacophore mapping approach,” Nucleic Acids Research, vol. 38, no. 2, pp. W609–W614, 2010. View at: Publisher Site | Google Scholar
  13. A. Martin, M. E. Ochagavia, L. C. Rabasa, J. Miranda, J. Fernandez-de-Cossio, and R. Bringas, “BisoGenet: a new tool for gene network building, visualization and analysis,” Bmc Bioinformatics, vol. 11, no. 1, pp. 91–99, 2010. View at: Publisher Site | Google Scholar
  14. Z. Tang, C. Li, B. Kang, G. Gao, C. Li, and Z. Zhang, “GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses,” Nucleic Acids Research, vol. 45, no. W1, pp. W98–W102, 2017. View at: Publisher Site | Google Scholar
  15. X.-J. Zhang and S.-S. Jia, “Fisetin inhibits laryngeal carcinoma through regulation of AKT/NF-κB/mTOR and ERK1/2 signaling pathways,” Biomedicine & Pharmacotherapy, vol. 83, pp. 1164–1174, 2016. View at: Publisher Site | Google Scholar
  16. A. Klimaszewska-Wisniewska, M. Halas-Wisniewska, T. Tadrowski, M. Gagat, D. Grzanka, and A. Grzanka, “Paclitaxel and the dietary flavonoid fisetin: a synergistic combination that induces mitotic catastrophe and autophagic cell death in A549 non-small cell lung cancer cells,” Cancer Cell International, vol. 16, no. 1, 2016. View at: Publisher Site | Google Scholar
  17. C.-Y. Chan, C.-H. Lien, M.-F. Lee, and C.-Y. Huang, “Quercetin suppresses cellular migration and invasion in human head and neck squamous cell carcinoma (HNSCC),” Biomedicine, vol. 6, no. 3, p. 15, 2016. View at: Publisher Site | Google Scholar
  18. H. I. Swanson, E.-Y. Choi, W. B. Helton, C. G. Gairola, and J. Valentino, “Impact of apigenin and kaempferol on human head and neck squamous cell carcinoma,” Oral Surgery, Oral Medicine, Oral Pathology and Oral Radiology, vol. 117, no. 2, pp. 214–220, 2014. View at: Publisher Site | Google Scholar
  19. J. Li, L. Hu, T. Zhou et al., “Taxifolin inhibits breast cancer cells proliferation, migration and invasion by promoting mesenchymal to epithelial transition via β-catenin signaling,” Life Sciences, vol. 232, p. 116617, 2019. View at: Publisher Site | Google Scholar
  20. C. H. Park, E. R. Hahm, J. H. Lee, K. C. Jung, and C. H. Yang, “Inhibition of β-catenin-mediated transactivation by flavanone in AGS gastric cancer cells,” Biochemical and Biophysical Research Communications, vol. 331, no. 4, pp. 1222–1228, 2005. View at: Publisher Site | Google Scholar
  21. E. Kim, J. Lee, and G. Kim, “Anti-carcinogenic effect of a new analogue 4’-chloroflavanone from flavanone in human breast cancer cells,” International Journal of Molecular Medicine, vol. 25, no. 2, pp. 293–298, 2009. View at: Publisher Site | Google Scholar
  22. T. I. Lee and R. A. Young, “Transcriptional regulation and its misregulation in disease,” Cell, vol. 152, no. 6, pp. 1237–1251, 2013. View at: Publisher Site | Google Scholar
  23. T. L. Wong, N. Che, and S. Ma, “Reprogramming of central carbon metabolism in cancer stem cells,” Biochimica Et Biophysica Acta Molecular Basis of Disease, vol. 1863, no. 7, pp. 1728–1738, 2017. View at: Publisher Site | Google Scholar
  24. A. J. McMichael, “Increases in laryngeal cancer in Britain and Australia in relation to alcohol and tobacco consumption trends,” Lancet, vol. 311, no. 8076, pp. 1244–1247, 1978. View at: Publisher Site | Google Scholar
  25. T. Otto and S. Piotr, “Cell cycle proteins as promising targets in cancer therapy,” Nature Reviews Cancer, vol. 17, no. 2, pp. 93–115, 2017. View at: Publisher Site | Google Scholar
  26. S. Q. Zhou, Q. L. Luo, X. Tan et al., “Erchen decoction plus huiyanzhuyu decoction inhibits the cell cycle, migration and invasion and induces the apoptosis of laryngeal squamous cell carcinoma cells,” Journal of Ethnopharmacology, vol. 256, 2020. View at: Publisher Site | Google Scholar
  27. H. W. Moon, H. G. Han, and Y. J. Jeon, “Protein quality control in the endoplasmic reticulum and cancer,” International Journal of Molecular Sciences, vol. 19, no. 10, p. 3020, 2018. View at: Publisher Site | Google Scholar
  28. E. Mohamed, Y. Cao, and P. C. Rodriguez, “Endoplasmic reticulum stress regulates tumor growth and anti-tumor immunity: a promising opportunity for cancer immunotherapy,” Cancer Immunology Immunotherapy Cii, vol. 66, pp. 1069–1078, 2017. View at: Publisher Site | Google Scholar

Copyright © 2021 Feng Xiang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views201
Downloads378
Citations

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.