Research Paper Volume 13, Issue 24 pp 26003—26021

Network pharmacology reveals the potential mechanism of Baiying Qinghou decoction in treating laryngeal squamous cell carcinoma

Kun Gao1,2, , Yanan Zhu3, , Hui Wang4, , Xianwei Gong5, , Zhiyong Yue1,2, , Aiai Lv3, , Xuanchen Zhou1,2, ,

  • 1 Department of Otorhinolaryngology Head and Neck Surgery, Shandong Provincial Hospital Affiliated to Shandong First Medical University, Jinan 250021, Shandong, China
  • 2 Department of Otorhinolaryngology Head and Neck Surgery, Shandong Provincial Hospital Affiliated to Shandong University, Jinan 250021, Shandong, China
  • 3 Department of Internal Medicine, Shandong Provincial Chest Hospital Affiliated to Shandong University, Jinan 250013, Shandong, China
  • 4 Department of Ultrasound, The Fifth People’s Hospital of Jinan, Jinan 250022, Shandong, China
  • 5 Department of Pharmacy, Shandong Provincial Hospital Affiliated to Shandong First Medical University, Jinan 250021, Shandong, China

Received: June 24, 2021       Accepted: November 22, 2021       Published: December 20, 2021      

https://doi.org/10.18632/aging.203786
How to Cite

Copyright: © 2021 Gao et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Context: Baiying Qinghou as a traditional Chinese medicine decoction shows anticancer property on laryngeal squamous cell carcinoma. However, little is known about the precise mechanism of Baiying Qinghou detection against laryngeal squamous cell carcinoma.

Objective: This study was aimed to explore potential mechanism of therapeutic actions of Baiying Qinghou decoction on laryngeal squamous cell carcinoma.

Materials and Methods: The active chemical components of Baiying Qinghou decoction were predicted, followed by integrated analysis of network pharmacology and molecular docking approach. The network pharmacology approach included target protein prediction, protein-protein interaction network construction and functional enrichment analysis.

Results: Sitosterol and quercetin were predicted to be the overlapped active ingredients among three Chinese herbs of Baiying Qinghou decoction. The target proteins were closely associated with response to chemical, response to drug related biological process and cancer related pathways such as PI3K-Akt signaling, HIF-1 signaling and Estrogen signaling pathway. The target proteins of TP53, EGFR, PTGS2, NOS3 and IL1B as the key nodes in PPI network were cross-validated, among which EGFR, IL1B, NOS3 and TP53 were significantly correlated with the prognosis of patients with laryngeal squamous cell carcinoma. Finally, the binding modes of EGFR, IL1B, NOS3 and TP53 with quercetin were visualized.

Discussion and Conclusion: Quercetin of Baiying Qinghou decoction showed therapeutic effect against laryngeal squamous cell carcinoma by regulating TP53, EGFR, NOS3 and IL1B involved with drug resistance and PI3K-AKT signaling pathway. TP53, EGFR, NOS3 and IL1B may be the candidate targets for the treatment of laryngeal squamous cell carcinoma.

Introduction

Head and neck cancer (HNC) is the seventh most common malignant tumor worldwide and occurs in nasopharynx, larynx, and thyroid [1, 2]. Head and neck squamous cell carcinoma (HNSCC) accounts for over 90% cases of head and neck cancer (HNC) and mainly derives from the mucosal surfaces of the upper aerodigestive tract [3]. Laryngeal squamous cell carcinoma (LSCC) represents the second most common histological subtypes of HNC with an increasing incidence rate [4]. Patients with LSCC are usually diagnosed at a late clinical stage and the survival rate is lower due to regional or distant metastases. Thus, LSCC has caused noticeable medical and economic burden worldwide.

Currently, several mainstream options against LSCC include surgery, radiotherapy, chemotherapy and chemo-radiotherapy [5]. However, the overall survival of patients with LSCC has not been remarkably improved because of chemotherapy or drug resistance and undesired effects. Traditional Chinese medicines (TCM) has been considered as an attractive alternative therapy for conquering cancers in China [6]. TCM focuses on restoring body balance and boosting immunity by the synergistic effects of various active ingredients [7]. Some Chinese herbs have been frequently utilized for the treatment of different diseases, such as cardiovascular diseases, diabetes and caners [810]. A recent research highlights that erchen plus huiyanzhuyu decoction can inhibit the growth of laryngeal carcinoma by modulating STAT3/cyclin D1 signaling pathway [11]. TCM exerts promising therapeutic effect on various malignant tumors.

Network pharmacology is proposed by Hopkins et al. and aims to explore the multilevel interactions of diseases, genes, and drugs as a whole [12]. This systems pharmacology is based on systems biology, computational biology and omics theory to evaluate the therapeutic effects of Chinese medicines on several diseases [13]. Gao et al. indicate that 8 herbs regulate multiple hepatocellular carcinoma-related genes and are strongly associated with prognosis by a network pharmacology approach [14]. Huang et al. suggest Huanglian Jiedu decoction has therapeutic roles in cancers such as hepatocellular carcinoma by building the herb-compound, compound-protein, protein-pathway, and gene-disease networks [15]. Notably, Bai Ying Qing Hou (BYQH) decoction is a Chinese medicinal formula and widely used for the treatment of LSCC in China, which is primarily composed of five key traditional Chinese herbs: Solanum lyratum thumb (30 g), Scutellaria barbata (24 g), Duchesnea indica (24 g), Solanum nigrum (30 g), and Actinidia chinensis planch (30 g). However, the pharmacological mechanism of BYQH decoction in the treatment of LSCC has not been clarified.

In this study, we attempted to explore the mechanism of BYQH decoction in the treatment of LSCC by network pharmacology analysis. First, the key active chemical constituents and their target proteins of BYQH decoction were screened. Then, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of gene targets were carried out. Besides, protein-protein interaction (PPI) and target proteins were cross-validated followed by survival analysis. Finally, molecular docking was used to identify the binding mode of key active ingredients and target proteins of BYQH decoction.

Materials and Methods

Screening active constituents of BYQH decoction

The traditional Chinese medicine systems pharmacology database (TCMSP) is an herbal repository and provides the system information about Chinese herbal medicines [16]. The absorption, distribution, metabolism, and excretion (ADME) system has been successfully used to evaluate the pharmacokinetics characteristics of chemical compounds [17, 18]. Herein, the chemical composition of five key components of BYQH decoction was firstly searched from TCMSP, mainly including the number of the pharmaceutical ingredient, molecule name, molecular weight, lipid-water partition coefficient, the number of hydrogen-bonding donor or acceptor, oral bioavailability (OB), blood brain barrier, Caco-2 permeability (Caco-2), drug-likeness (DL), and half-life (HL). OB ≥ 30% and Caco-2 ≥ −0.4 indicate good OB and permeability of molecules. The mean DL value of drugs collected in database is 0.18. Then, active chemical compounds were screened in ADME system according to the criteria of OB ≥ O30%, DL ≥ 0.18 and Caco-2 ≥ −0.4 according the previous description [19].

Screening protein targets of active chemical compounds

TCMSP and DrugBank (https://www.drugbank.ca/) databases can predict the relationships between drugs and corresponding targets. Moreover, DrugBank database provides the detailed information on experimental and investigational drugs [20]. The active chemical compounds were subjected to TCMSP and DrugBank databases to search their potential targets. Then, the target proteins were mapped to corresponding gene symbols using string database. For those proteins that did not map to any gene symbols, they were manually retrieved by Universal Protein (UniProt) database to map to corresponding gene symbols. These potential gene targets were used for following analysis.

Functional enrichment analysis of target genes

To investigate the biological function of the target genes, the gene ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were carried out. There are three categories for GO terms, including biological process (BP), cellular component (CC) and molecular function (MF). Gene symbols were input into string database. The parameters were set as the default values. The significantly enriched GO and pathway terms with p < 0.05 were collected. Then, the top 10 GO terms in each category and top 10 significant pathways were visualized in bubble chart by using R language.

Target protein-protein interaction (PPI) analysis

The underlying relationships among target proteins were analyzed based on the string database according to the default parameters [21]. The Cytoscape software was used to construct and visualize the PPI network. Following this, the topological characteristics of the PPI network were also evaluated. The node degree (connectivity) was calculated.

Cross-validation of target proteins by the database retrieval

DisGeNET is a versatile platform that comprehensively collects human gene-disease associations (over 380,000 associations between more than 16,000 genes and 13,000 diseases) for the validation of computationally predicted genes for human diseases [22, 23]. In this study, the genes related to LSCC were firstly searched from this database using “carcinoma of larynx” as the keyword. Phenopedia database also provides systematical genetic association studies and summarizes the information about the association between genes studied and a particular disease. Therefore, the known genes associated with LSCC were obtained by retrieving Phenopedia database with the keyword of “laryngeal neoplasms”. Subsequently, the overlapping genes among target genes of BYQH decoction, LSCC-related genes in DisGeNET or Phenopedia database were extracted by jvenn tool [24, 25]. These genes may be associated with the therapeutic effects of BYQH decoction against LSCC.

Survival analysis of candidate genes

Firstly, expression profiles of HNSCC samples deposited in the Cancer Genome Atlas (TCGA) dataset were downloaded from UCSC Xena (https://xenabrowser.net/datapages/). Based on the clinical information, the HNSCC samples were selected according to following criteria: 1) samples were primary tumor tissues and the primary site was at the larynx; 2) the clinical information of samples with regard to TNM stage, clinical stage, grade and follow-up were complete; 3) the samples without survival time were deleted. Finally, a total of 105 laryngeal cancer samples were retained for the next prognostic analysis.

To investigate whether key candidate genes are associated with prognosis of patients with laryngeal cancer, the survival analyses of these gene targets were carried out. In brief, the raw data of 105 laryngeal cancer samples were firstly extracted and then standardized by log2 (count+1). After that, the optimal cutoff was determined using survminer package in R language based on the expressions of candidate genes. Subsequently, the samples were divided into high- and low- expression group according to the optimal cutoff. Following this, the survival analyses of gene targets were respectively preformed with survival package in R language. P < 0.05 was considered significant.

Predicting the binding mode between active compounds and key protein targets of BYQH decoction

The Protein Data Bank (PDB) offers the experimental data for determined three-dimensional (3D) biological macromolecules structure [26]. Herein, we would further explore the underlying molecular mechanism how active compounds interacted with target proteins of BYQH decoction. The PDB ID for the key protein was firstly acquired from PDB. Subsequently, the binding mode between key protein and corresponding active chemical ingredient was predicted by the LeDock tool. Finally, we used the Pymol software to visualize the predicted binding modes of key active ingredients to their corresponding targets.

Ethics approval and consent to participate

This study was approved by Ethics Committee of Provincial Hospital Affiliated to Shandong First Medical University and Shandong Provincial Hospital affiliated to Shandong University.

Availability of data and materials

The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.

Highlights

The therapeutic mechanism of BYQH decoction on LSCC is analyzed by network pharmacology. Sitosterol and quercetin were key active ingredients BYQH decoction. TP53, EGFR, NOS3 and IL1B may be therapeutic targets of LSCC.

Results

Active chemical compounds of BYQH decoction

The active chemical compounds of BYQH were retrieved from TCMSP database by using five keywords of “Solanum lyratum thumb” or “Scutellaria barbata” or “Duchesnea indica” or “Solanum nigrum” or “Actinidia chinensis Planch”. Totally, 174 active ingredients were obtained: 15 in Solanum lyratum thumb, 39 in Solanum nigrum, 94 in Scutellaria barbata and 26 in Actinidia chinensis Planch, respectively. Notably, the active components of Duchesnea indica were not included in TCMSP database. Subsequently, 41 key active compounds (7 in Solanum nigrum, 28 in Scutellaria barbata and 6 in Actinidia chinensis Planch) were further screened according to the ADME parameters of OB ≥ 30%, DL ≥ 0.18 and Caco-2 ≥ −0.4 (Table 1; Figure 1). Interestingly, sitosterol and quercetin were identified in three Chinese herbs (Solanum nigrum, Scutellaria barbata and Actinidia chinensis Planch; Table 1). In addition, beta-sitosterol was identified in Scutellaria barbata and Actinidia chinensis Planch, and cholesterol was identified in Solanum nigrum and Scutellaria barbata (Table 1).

Table 1. The list of active components of Baiying Qinghou decoction.

CompoundsMol IDMolecule NameMWAlogPHdonHaccOB (%)Caco-2BBBDLHL
S.b.D.DMOL001040(2R)-5,7-dihydroxy-2-(4-hydroxyphenyl)chroman-4-one272.272.33542.360.38−0.480.2116.83
S.b.D.DMOL0122455,7,4′-trihydroxy-6-methoxyflavanone302.32.283636.630.43−0.320.2716.12
S.b.D.DMOL0122465,7,4′-trihydroxy-8-methoxyflavanone302.32.283674.240.37−0.430.2616.85
S.b.D.DMOL0122485-hydroxy-7,8-dimethoxy-2-(4-methoxyphenyl)chromone328.342.821665.820.850.070.3316.41
S.b.D.DMOL0122507-hydroxy-5,8-dimethoxy-2-phenyl-chromone298.312.841543.720.960.220.2516.77
S.b.D.DMOL012251Chrysin-5-methylether268.282.851437.270.910.160.217.24
S.b.D.DMOL0122529,19-cyclolanost-24-en-3-ol426.87.551138.691.451.160.785.41
S.b.D.DMOL012254campesterol400.767.631137.581.340.980.714.63
S.b.D.DMOL000953CLR386.737.381137.871.431.130.684.52
S.b.D.DMOL000358beta-sitosterol414.798.081136.911.320.990.755.36
S.b.D.DMOL012266rivularin344.342.552737.940.65−0.130.3716.25
S.b.D.DMOL001973Sitosterol acetate456.838.460240.391.391.110.856.34
S.b.D.DMOL012269Stigmasta-5,22-dien-3-ol-acetate454.818.020246.441.411.060.866.77
S.b.D.DMOL012270Stigmastan-3,5,22-triene394.758.430045.031.91.810.716.21
S.b.D.DMOL000449Stigmasterol412.777.641143.831.4410.765.57
S.b.D.DMOL000173wogonin284.282.592530.680.790.040.2317.75
S.b.D.DMOL001735Dinatin300.282.323630.970.48−0.490.2716.44
S.b.D.DMOL00175524-Ethylcholest-4-en-3-one412.778.180136.081.461.220.765.49
S.b.D.DMOL002714baicalein270.252.333533.520.63−0.050.2116.25
S.b.D.DMOL0027196-Hydroxynaringenin288.272.034633.230.27−0.270.2415.67
S.b.D.DMOL002915Salvigenin328.342.821649.070.86−0.030.3315.87
S.b.D.DMOL000351Rhamnazin330.312.013747.140.53−0.320.3413.54
S.b.D.DMOL000359sitosterol414.798.081136.911.320.870.755.37
S.b.D.DMOL005190eriodictyol288.272.034671.790.17−0.540.2415.81
S.b.D.DMOL005869daucostero_qt414.798.081136.911.320.870.755.08
S.b.D.DMOL000006luteolin286.252.074636.160.19−0.840.2515.94
S.b.D.DMOL008206Moslosooflavone298.312.841544.091.010.540.2517.02
S.b.D.DMOL000098quercetin302.251.55746.430.05−0.770.2814.4
S.n.LMOL00205840957-99-1388.452.122757.20.49−0.290.622.04
S.n.LMOL002773beta-carotene536.96120037.182.251.520.584.36
S.n.LMOL000359sitosterol414.798.081136.911.320.870.755.37
S.n.LMOL000546diosgenin414.694.631380.880.820.270.814.14
S.n.LMOL007356solanocapsine430.753.494452.940.39−0.220.677.86
S.n.LMOL000953CLR386.737.381137.871.431.130.684.52
S.n.LMOL000098quercetin302.251.55746.430.05−0.770.2814.4
A.c.PMOL000358beta-sitosterol414.798.081136.911.320.990.755.36
A.c.PMOL000359sitosterol414.798.081136.911.320.870.755.37
A.c.PMOL000471aloe-emodin270.251.673583.38−0.12−1.070.2431.49
A.c.PMOL000492(+)-catechin290.291.925654.83−0.03−0.730.240.61
A.c.PMOL000073ent-Epicatechin290.291.925648.960.02−0.640.240.63
A.c.PMOL000098quercetin302.251.55746.430.05−0.770.2814.4
Abbreviations: S.b.D.D: Scutellaria barbata D. Don; S.n.L: Solanum nigrum L; A.c.P: Actinidia chinensis Planch; Mol ID: molecule ID; MW: molecular weight; AlogP: lipid-water partition coefficient; Hdon/Hacc: the number of hydrogen-bonding donor or acceptor; OB: oral bioavailability; BBB: blood brain barrier; Caco-2: Caco-2 permeability; DL: drug-likeness; HL: half-life.
Venn diagrams for active ingredients of Chinese herbal medicine composition in Baiying Qinghou decoction. Abbreviations: S.b.D.D: Scutellaria barbata D. Don; S.n.L: Solanum nigrum L; A.c.P: Actinidia chinensis Planch.

Figure 1. Venn diagrams for active ingredients of Chinese herbal medicine composition in Baiying Qinghou decoction. Abbreviations: S.b.D.D: Scutellaria barbata D. Don; S.n.L: Solanum nigrum L; A.c.P: Actinidia chinensis Planch.

Functional enrichment analysis of targets of the active Ingredients of BYQH decoction

Totally, 137 gene targets of the active ingredients of BYQH decoction were predicted according to the methods mentioned above (Supplementary Table 1). Afterwards, the functional enrichment analyses of these genes were conducted to investigate their potential biological roles in LSCC treatment. The results showed that they were enriched in 1444 GO-BP terms, 114 GO-CC terms and 203 GO-MF terms. Figure 2 showed top 10 significantly enriched GO terms in MF, CC, and BP. At the cellular level, the gene targets were dramatically enriched in cellular response to chemical stimulus (involved with TP53, EGFR, NOS3 and IL1B et al.), response to drug and response to chemical (involved with TP53, EGFR, NOS3 and IL1B et al.). The membrane raft, plasma membrane region and plasma membrane raft were the three GO-CC terms enriched by gene targets. For GO-MF category, these genes were significantly involved in protein binding, enzyme biding and molecular transducer activity. In addition, the KEGG enrichment analysis revealed that the target genes were enriched in 172 pathways. Moreover, we found they were primarily correlated with several cancer-related pathways, such as PI3K-Akt signaling pathway (involved with BCL2, EGFR, NOS3 and TP53 et al.), estrogen signaling pathway (involved with BCL2, EGFR, MMP2, NOS3 et al.) and HIF-1 signaling pathway (involved with EGF, EGFR, IL6 and NOS3 et al.). These findings suggested that BYQH decoction may have essential effect on the treatment of LSCC.

Functional analyses of target genes of Baiying Qinghou decoction. The top 10 significantly enriched GO-BP/MF/CC terms and KEGG pathways were displayed. Abbreviations: GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; BP: biological process; CC: cellular component; MF: molecular function.

Figure 2. Functional analyses of target genes of Baiying Qinghou decoction. The top 10 significantly enriched GO-BP/MF/CC terms and KEGG pathways were displayed. Abbreviations: GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; BP: biological process; CC: cellular component; MF: molecular function.

Target PPI network analysis

The potential interactions among target proteins of BYQH decoction was constructed and visualized by STRING database. As shown in Figure 3, there were 128 protein nodes and 1282 edges in PPI network. The protein nodes with higher degree may have closer biological connection with other nodes. Herein, the top 15 nodes were regarded as hub genes and listed in Table 2, including IL6 (interleukin-6 degree = 67), VEGFA (vascular endothelial growth factor A; degree = 63), TP53 (degree = 61), JUN (degree = 58), EGF (degree = 57), MAPK1 (degree = 54), EGFR (degree = 53), PTGS2 (degree = 52), ESR1 (degree = 51), CAT (degree = 49), NOS3 (degree = 48), IL1B (degree = 47), HSP90AA1 (degree = 44), CCL2 (degree = 44), and AR (degree = 43).

The protein-protein interaction (PPI) network of key target genes of Baiying Qinghou decoction. Dots represent target proteins/genes and the lines represents interactions among proteins/targets. The larger dot size shows the higher degree.

Figure 3. The protein-protein interaction (PPI) network of key target genes of Baiying Qinghou decoction. Dots represent target proteins/genes and the lines represents interactions among proteins/targets. The larger dot size shows the higher degree.

Table 2. The top 15 hub proteins in protein-protein interaction network.

NameDegreeBetweenness centralityCloseness centrality
IL6670.0800.658
VEGFA630.0500.638
TP53610.0310.608
JUN580.0430.626
EGF570.0580.623
MAPK1540.0220.593
EGFR530.0350.608
PTGS2520.0300.599
ESR1510.0310.585
CAT490.0580.591
NOS3480.0630.599
IL1B470.0210.593
HSP90AA1440.0230.561
CCL2440.0130.559
AR430.0290.555

Cross-validation of target proteins

To further narrow the range of potential target genes of BYQH decoction, the candidate target genes were cross-validated. The genes related to LSCC were respectively acquired from Phenopedia and DisGeNET database, which provided information on gene-disease interactions. Then, the intersecting target genes were obtained between predicted genes and gene targets of BYQH decoction. Finally, a total of 10 intersection-associated genes were extracted, including PTGS2, NOS3, BCL2, ADH1C, TP53, MPO, EGFR, GSTP1, IL1B and GSTM1 (Figure 4). Among them, five genes (TP53, EGFR, PTGS2, NOS3 and IL1B) served as the hub genes in PPI network. Therefore, these genes were considered to be key targets of BYQH decoction against LSCC.

Venn diagrams for cross-validation gene targets by Phenopedia and DisGeNET databases. The intersecting genes were extracted between gene targets of Baiying Qinghou decoction and genes related to laryngeal cancer in Phenopedia and DisGeNET databases.

Figure 4. Venn diagrams for cross-validation gene targets by Phenopedia and DisGeNET databases. The intersecting genes were extracted between gene targets of Baiying Qinghou decoction and genes related to laryngeal cancer in Phenopedia and DisGeNET databases.

Survival analysis

To explore the underlying impacts of 10 key gene targets on the clinical prognosis of LSCC, we conducted the survival analyses using TCGA data of 105 HNSCC patients. The expression patterns of these genes were firstly determined among tissue samples of 105 HNSCC patients. Then, we calculated the optimal cutoff value of each gene to stratify all patients into high- and low-expression groups. Our results showed that the optimal cutoff value was 2.58 for ADH1C, 7.03 for BCL2, 15.13 for EGFR, 13.45 for GSTM1, 16.06 for GSTP1, 10.53 for IL1B, 5.29 for MPO, 8.5 for NOS3, 9.36 for PTGS2, and 10.95 for TP53 (Figure 5 and Supplementary Figure 1). Survival analysis revealed that seven genes (ADH1C, EGFR, GSTM1, GSTP1, IL1B, NOS3 and TP53) were significantly associated with the prognosis of HNSCC patients (P < 0.05; Figure 5). Moreover, HNSCC patients with a high expression level of ADH1C had a better prognosis (P = 0.012; Figure 5). However, the high expression levels of the remaining six genes were closely related to the unfavorable prognosis of HNSCC patients (EGFR: P = 0.014; GSTM1: P = 0.004; GSTP1: P = 0.034; IL1B: P = 0.03; NOS3: P = 0.039 and TP53: P = 0.02; Figure 5). Notably, TP53, EGFR, NOS3 and IL1B also acted as hub genes in PPI network.

Survival analysis of key target genes of Baiying Qinghou decoction using a TCGA dataset. The optimal cutoff value of each gene was calculated to stratify all patients into high- and low-expression groups (the left of the A–G). Seven genes (ADH1C, EGFR, GSTM1, GSTP1, IL1B, NOS3 and TP53) were significantly associated with the prognosis of HNSCC patients (the right of the A–G). The high expression level of ADH1C (A) had a better prognosis of HNSCC patients (P = 0.012). However, the high expression levels of EGFR (B), GSTM1 (C), GSTP1 (D), IL1B (E), NOS3 (F) and TP53 (G) showed the poor prognosis of HNSCC patients. Abbreviations: HNSCC: head and neck squamous cell carcinoma; TCGA: The Cancer Genome Atlas.

Figure 5. Survival analysis of key target genes of Baiying Qinghou decoction using a TCGA dataset. The optimal cutoff value of each gene was calculated to stratify all patients into high- and low-expression groups (the left of the AG). Seven genes (ADH1C, EGFR, GSTM1, GSTP1, IL1B, NOS3 and TP53) were significantly associated with the prognosis of HNSCC patients (the right of the AG). The high expression level of ADH1C (A) had a better prognosis of HNSCC patients (P = 0.012). However, the high expression levels of EGFR (B), GSTM1 (C), GSTP1 (D), IL1B (E), NOS3 (F) and TP53 (G) showed the poor prognosis of HNSCC patients. Abbreviations: HNSCC: head and neck squamous cell carcinoma; TCGA: The Cancer Genome Atlas.

Prediction of the binding mode between key bioactive compounds and proteins

The binding characteristics of key active ingredients of BYQH decoction and four key protein targets (TP53, EGFR, NOS3 and IL1B) were investigated. Herein, we focused on analyzing the binding modes between two overlapped active compounds (sitosterol and quercetin) among three Chinese herbs and four target proteins. Our analysis indicated that the corresponding protein targets of sitosterol were PGR, NCOA2, NR3C2, which were not targets studied. For quercetin, 75 target proteins were obtained, including four key protein targets. Subsequently, the chemical structure formula of quercetin was acquired from TSMCP database, and four target proteins were also searched from PDB database, including TP53 (PDB ID: 3ZME), EGFR (PDB ID: 1XKK), NOS3 (PDB ID: 6PP1) and IL1B (PDB ID: 5R86). Finally, the molecular docking was carried out using LeDock tool. Notably, a closer binding between proteins and small biological molecules indicated more energy released and a lower ΔG value. Figure 6 depicted the optimal binding modes of quercetin and its four target proteins.

The binding modes. (A) The binding modes of quercetin and EGFR. The optimal binding modes of quercetin and EGFR (ΔG = −6.87 kcal/mol). Left: two-dimensional image; Right: three-dimensional image. (B) The binding modes of quercetin and IL1B. The optimal binding modes of quercetin and IL1B (ΔG = −4.84 kcal/mol). Left: two-dimensional image; Right: three-dimensional image. (C) The binding modes of quercetin and NOS3. The optimal binding modes of quercetin and NOS3 (ΔG = −7.12 kcal/mol). Left: two-dimensional image; Right: three-dimensional image. (D) The binding modes of quercetin and TP53. The optimal binding modes of quercetin and TP53 (ΔG = −6.15 kcal/mol). Left: two-dimensional image; Right: three-dimensional image.

Figure 6. The binding modes. (A) The binding modes of quercetin and EGFR. The optimal binding modes of quercetin and EGFR (ΔG = −6.87 kcal/mol). Left: two-dimensional image; Right: three-dimensional image. (B) The binding modes of quercetin and IL1B. The optimal binding modes of quercetin and IL1B (ΔG = −4.84 kcal/mol). Left: two-dimensional image; Right: three-dimensional image. (C) The binding modes of quercetin and NOS3. The optimal binding modes of quercetin and NOS3 (ΔG = −7.12 kcal/mol). Left: two-dimensional image; Right: three-dimensional image. (D) The binding modes of quercetin and TP53. The optimal binding modes of quercetin and TP53 (ΔG = −6.15 kcal/mol). Left: two-dimensional image; Right: three-dimensional image.

Discussion

TCM is an accepted medical practice and has been used for the treatment of complicated diseases such as cancers over past few decades in China [27, 28]. BYQH decoction, a TCM prescription, has been clinically proven to be effective as adjuvant treatment against LSCC. In this study, we performed a network pharmacological analysis to explore pharmacodynamic effects and therapeutic mechanism of BYQH decoction for LSCC. We identified 41 key active ingredients of BYQH decoction. Of these, sitosterol and quercetin were intersected among three Chinese herbs of BYQH decoction. Moreover, these crucial components were corresponded to 137 target proteins. Top 15 genes (IL6, VEGFA, TP53, JUN, EGF, MAPK1, EGFR, PTGS2, ESR1, CAT, NOS3, IL1B, HSP90AA1, CCL2 and AR) were considered as hub nodes in PPI network. Among them, five proteins (TP53, EGFR, PTGS2, NOS3 and IL1B) were cross-validated. Additionally, survival analysis showed the expressions of ADH1C, EGFR, GSTM1, GSTP1, IL1B, NOS3 and TP53 were markedly correlated with the prognosis of HNSCC patients. Finally, the binding modes of quercetin and four target proteins (TP53, EGFR, NOS3 and IL1B) were constructed.

Quercetin (3,3′,4′,5,7-pentahydroxyflavone), a common flavonol, is widely distributed in plant species, such as vegetables and grains. Our results showed that quercetin was a key active chemical compound of BYQH decoction. Extensive evidence suggests that quercetin exerts diverse biological functions, such as antioxidant, anti-inflammatory, anti-bacterial, antiviral and anticarcinogenic properties [2931]. Numerous studies have evaluated the anticarcinogenic effects of quercetin and indicated that quercetin was involved in regulating several cancer-related pathways, such PI3K/Akt/ signaling pathways and MAPK/ERK1/2 pathways [32, 33]. Sharma et al. report that quercetin induced human laryngeal HeP2 cells death and synergistically enhanced the antiproliferative ability of cisplatin in these cells [34]. Numerous researchers demonstrate that quercetin (50 μM) can significantly increase photodynamic therapy-induced cytotoxicity via reducing the cell viability of human larynx carcinoma cells (HEp-2) [35]. In addition, the apoptotic events are increased stronger in HEp-2 cells with the combination of quercetin than the single drug treatment [36]. These evidences reveals that quercetin may be a promising therapeutic molecule for LSCC treatment.

TP53, EGFR, NOS3 and IL1B were key hub genes in PPI network of target proteins of BYQH decoction and cross-validated in two databases (Phenopedia and DisGeNET). More notably, these four proteins were also target proteins of quercetin. p53 protein encoded by TP53 gene functions as a tumor-suppressive factor by regulating multiple cellular processes such as cell-cycle progression and apoptosis [37]. Overwhelming evidence suggests several mutations in TP53 can influence the progression of HNSCC and clinical treatment response [38, 39]. Clemente-Soto AF et al. point out that quercetin can induce G2 phase cell cycle arrest and apoptosis in human cervical cancer cells, accompanied by upregulating p53 level [40]. Our molecular docking analysis suggested that quercetin can strongly bind with TP53. EGFR is frequently reported to be associated with HNSCC [41]. Chan et al. indicate that quercetin may suppress cell migration and invasion of HNSCC cells overexpressing EGFR by down-regulating the expression of MMP-2 and MMP-9 [42]. Similarly, Chung et al. indicate that there is a high frequency of EGFR copy number in HNSCC, which acts as a poor predictor for the prognosis of HNSCC patients [43]. IL1B is an inflammatory cytokine gene and acts as a therapeutic target for the treatment of head and neck cancer [44, 45]. NOS3 is located on chromosome 7 (7q36) and regulates at transcriptional and post-transcriptional levels [46]. Previous studies identified numerous polymorphic sites of NOS3 such as single nucleotide polymorphism and insertion or deletion [47, 48]. Moreover, the NOS3 polymorphisms play crucial roles in the molecular mechanism of cancers and clinical survivals of patients undergoing cancers, including laryngeal cancer [4951]. Guo et al. recently report that quercetin has a binding interaction with NOS3, which was consistent with our finding [52]. Besides, we also found the elevated level of NOS3 represented an unfavorable prognosis.

Besides, in this study, GO enrichment analysis showed that the four genes of interest (TP53, EGFR, NOS3 and IL1B) were significantly enriched in cellular response to chemical stimulus, response to chemical, response to drug related biological processes. TP53 is a short-lived protein that plays a central role in mediating cellular response to stressful and genotoxic stimuli, such as anticancer drugs exposure [53]. EGFR has been proposed to be the prognostic marker for HNSCC that closely related with radiation sensitivity, tumor size and recurrence [54]. Additionally, our survival data showed that the high expression of TP53, EGFR, NOS3 and IL1B in HNSCC patients exhibited a poor prognosis, which was consistent with the previous reports mentioned above. Thus, we speculated that the quercetin might reduce drug resistance of HNSCC by targeting TP53, EGFR, NOS3 and IL1B.

Furthermore, PI3K-AKT signaling is considered to be a classical pathway related with various cancers [55, 56]. The descending activation of PI3K-AKT signaling is closely associated with radiosensitivity enhancement of LSCC patients [57]. The reduced protein expressions of PI3K and AKT are accompanied by suppressed LSCC tumor growth [58]. In this study, PI3K-AKT signaling was identified to a significant pathway involved with TP53, EGFR, NOS3. Our findings supported the significant role of PI3K-AKT signaling pathway in LSCC and proposed the target role of the key proteins identified in this study.

Although our study identified several key components and target proteins of BYQH decoction, relevant experimental assays are needed to validate our findings. Although were considered as the prognostic indicators for patients with LSCC, The expression patterns of the four candidate prognostic indicators in LSCC have not been investigated in clinical samples. Moreover, the relationships of these prognostic genes and some clinical parameters are also not analyzed. Additionally, many chemical compounds are possibly not identified because the limited data is available in public databases, which may lead to some false positive results.

Conclusions

Network pharmacology analysis was carried out to elucidate the therapeutic role of BYQH decoction on LSCC. Quercetin was predicted to be the active compounds of BYQH decoction by targeting TP53, EGFR, NOS3 and IL1B involved in drug resistance and PI3K-AKT signaling pathway. Moreover, the four genes of TP53, EGFR, NOS3 and IL1B exerted promising prognostic potential for LSCC. Our results offered an important clue for uncovering the underlying mechanism of BYQH decoction in the treatment of LSCC. However, more experimental evidences about the effects of quercetin on LSCC in vivo and in vitro by targeting TP53, EGFR, NOS3 and IL1B are warranted in the near future.

Abbreviations

BYQH: Baiying Qinghou; LSCC: laryngeal squamous cell carcinoma; PPI: protein-protein interaction; HNC: Head and neck cancer; HNSCC: Head and neck squamous cell carcinoma; TCM: Traditional Chinese medicines; GO: Gene Ontology; OB: oral bioavailability, blood brain barrier; Caco-2: Caco-2 permeability; DL: drug-likeness; HL: half-life; BP: biological process; CC: cellular component; MF: molecular function.

Author Contributions

Aiai Lv, Xuanchen Zhou and Kun Gao carried out the Conception and design of the research, Kun Gao and Xianwei Gong participated in the Acquisition of data. Kun Gao and Yanan Zhu carried out the Analysis and interpretation of data. Kun Gao and Hui Wang participated in the design of the study and performed the statistical analysis. Zhiyong Yue participated in Obtaining funding. Aiai Lv, Xuanchen Zhou, Kun Gao and Aiai Lv conceived of the study, and participated in its design and coordination and helped to draft the manuscript and revision of manuscript for important intellectual content. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflicts of interest related to this study.

Funding

This work was supported by 2018 the Key R&D Program of Shandong Province (Program No. 2018GSF118192).

References

  • 1. Chow LQM. Head and Neck Cancer. N Engl J Med. 2020; 382:60–72. https://doi.org/10.1056/NEJMra1715715 [PubMed]
  • 2. Marur S, Forastiere AA. Head and Neck Squamous Cell Carcinoma: Update on Epidemiology, Diagnosis, and Treatment. Mayo Clin Proc. 2016; 91:386–96. https://doi.org/10.1016/j.mayocp.2015.12.017 [PubMed]
  • 3. Pedregal-Mallo D, Hermida-Prado F, Granda-Díaz R, Montoro-Jiménez I, Allonca E, Pozo-Agundo E, Álvarez-Fernández M, Álvarez-Marcos C, García-Pedrero JM, Rodrigo JP. Prognostic Significance of the Pluripotency Factors NANOG, SOX2, and OCT4 in Head and Neck Squamous Cell Carcinomas. Cancers (Basel). 2020; 12:1794. https://doi.org/10.3390/cancers12071794 [PubMed]
  • 4. Cossu AM, Mosca L, Zappavigna S, Misso G, Bocchetti M, De Micco F, Quagliuolo L, Porcelli M, Caraglia M, Boccellino M. Long Non-coding RNAs as Important Biomarkers in Laryngeal Cancer and Other Head and Neck Tumours. Int J Mol Sci. 2019; 20:3444. https://doi.org/10.3390/ijms20143444 [PubMed]
  • 5. Li Y, Wang K, Yin S, Zheng H, Min D. Xanthohumol inhibits proliferation of laryngeal squamous cell carcinoma. Oncol Lett. 2016; 12:5289–94. https://doi.org/10.3892/ol.2016.5313 [PubMed]
  • 6. So TH, Chan SK, Lee VH, Chen BZ, Kong FM, Lao LX. Chinese Medicine in Cancer Treatment - How is it Practised in the East and the West? Clin Oncol (R Coll Radiol). 2019; 31:578–88. https://doi.org/10.1016/j.clon.2019.05.016 [PubMed]
  • 7. Li J, Wu Y, Wang D, Zou L, Fu C, Zhang J, Leung GP. Oridonin synergistically enhances the anti-tumor efficacy of doxorubicin against aggressive breast cancer via pro-apoptotic and anti-angiogenic effects. Pharmacol Res. 2019; 146:104313. https://doi.org/10.1016/j.phrs.2019.104313 [PubMed]
  • 8. Guo R, Luo X, Liu J, Liu L, Wang X, Lu H. Omics strategies decipher therapeutic discoveries of traditional Chinese medicine against different diseases at multiple layers molecular-level. Pharmacol Res. 2020; 152:104627. https://doi.org/10.1016/j.phrs.2020.104627 [PubMed]
  • 9. Gerber B, Scholz C, Reimer T, Briese V, Janni W. Complementary and alternative therapeutic approaches in patients with early breast cancer: a systematic review. Breast Cancer Res Treat. 2006; 95:199–209. https://doi.org/10.1007/s10549-005-9005-y [PubMed]
  • 10. Zhong Y, Lee K, Deng Y, Ma Y, Chen Y, Li X, Wei C, Yang S, Wang T, Wong NJ, Muwonge AN, Azeloglu EU, Zhang W, et al. Arctigenin attenuates diabetic kidney disease through the activation of PP2A in podocytes. Nat Commun. 2019; 10:4523. https://doi.org/10.1038/s41467-019-12433-w [PubMed]
  • 11. Tan X, Luo Q, Zhou S, Huang W, Feng X, Chen W, Yang C, Li Y. Erchen Plus Huiyanzhuyu Decoction Inhibits the Growth of Laryngeal Carcinoma in a Mouse Model of Phlegm-Coagulation-Blood-Stasis Syndrome via the STAT3/Cyclin D1 Pathway. Evid Based Complement Alternat Med. 2020; 2020:2803496. https://doi.org/10.1155/2020/2803496 [PubMed]
  • 12. Hopkins AL. Network pharmacology. Nat Biotechnol. 2007; 25:1110–1. https://doi.org/10.1038/nbt1007-1110 [PubMed]
  • 13. Fang J, Liu C, Wang Q, Lin P, Cheng F. In silico polypharmacology of natural products. Brief Bioinform. 2018; 19:1153–71. https://doi.org/10.1093/bib/bbx045 [PubMed]
  • 14. Gao L, Wang XD, Niu YY, Duan DD, Yang X, Hao J, Zhu CH, Chen D, Wang KX, Qin XM, Wu XZ. Molecular targets of Chinese herbs: a clinical study of hepatoma based on network pharmacology. Sci Rep. 2016; 6:24944. https://doi.org/10.1038/srep24944 [PubMed]
  • 15. Huang J, Guo W, Cheung F, Tan HY, Wang N, Feng Y. Integrating Network Pharmacology and Experimental Models to Investigate the Efficacy of Coptidis and Scutellaria Containing Huanglian Jiedu Decoction on Hepatocellular Carcinoma. Am J Chin Med. 2020; 48:161–82. https://doi.org/10.1142/S0192415X20500093 [PubMed]
  • 16. Ru J, Li P, Wang J, Zhou W, Li B, Huang C, Li P, Guo Z, Tao W, Yang Y, Xu X, Li Y, Wang Y, Yang L. TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J Cheminform. 2014; 6:13. https://doi.org/10.1186/1758-2946-6-13 [PubMed]
  • 17. Dong J, Wang NN, Yao ZJ, Zhang L, Cheng Y, Ouyang D, Lu AP, Cao DS. ADMETlab: a platform for systematic ADMET evaluation based on a comprehensively collected ADMET database. J Cheminform. 2018; 10:29. https://doi.org/10.1186/s13321-018-0283-x [PubMed]
  • 18. Ertl P, Rohde B, Selzer P. Fast calculation of molecular polar surface area as a sum of fragment-based contributions and its application to the prediction of drug transport properties. J Med Chem. 2000; 43:3714–7. https://doi.org/10.1021/jm000942e [PubMed]
  • 19. Huang J, Cheung F, Tan HY, Hong M, Wang N, Yang J, Feng Y, Zheng Q. Identification of the active compounds and significant pathways of yinchenhao decoction based on network pharmacology. Mol Med Rep. 2017; 16:4583–92. https://doi.org/10.3892/mmr.2017.7149 [PubMed]
  • 20. Wishart DS, Feunang YD, Guo AC, Lo EJ, Marcu A, Grant JR, Sajed T, Johnson D, Li C, Sayeeda Z, Assempour N, Iynkkaran I, Liu Y, et al. DrugBank 5.0: a major update to the DrugBank database for 2018. Nucleic Acids Res. 2018; 46:D1074–82. https://doi.org/10.1093/nar/gkx1037 [PubMed]
  • 21. von Mering C, Jensen LJ, Snel B, Hooper SD, Krupp M, Foglierini M, Jouffre N, Huynen MA, Bork P. STRING: known and predicted protein-protein associations, integrated and transferred across organisms. Nucleic Acids Res. 2005; 33:D433–7. https://doi.org/10.1093/nar/gki005 [PubMed]
  • 22. Piñero J, Bravo À, Queralt-Rosinach N, Gutiérrez-Sacristán A, Deu-Pons J, Centeno E, García-García J, Sanz F, Furlong LI. DisGeNET: a comprehensive platform integrating information on human disease-associated genes and variants. Nucleic Acids Res. 2017; 45:D833–9. https://doi.org/10.1093/nar/gkw943 [PubMed]
  • 23. Piñero J, Queralt-Rosinach N, Bravo À, Deu-Pons J, Bauer-Mehren A, Baron M, Sanz F, Furlong LI. DisGeNET: a discovery platform for the dynamical exploration of human diseases and their genes. Database (Oxford). 2015; 2015:bav028. https://doi.org/10.1093/database/bav028 [PubMed]
  • 24. Yu W, Clyne M, Khoury MJ, Gwinn M. Phenopedia and Genopedia: disease-centered and gene-centered views of the evolving knowledge of human genetic associations. Bioinformatics. 2010; 26:145–6. https://doi.org/10.1093/bioinformatics/btp618 [PubMed]
  • 25. Bardou P, Mariette J, Escudié F, Djemiel C, Klopp C. jvenn: an interactive Venn diagram viewer. BMC Bioinformatics. 2014; 15:293. https://doi.org/10.1186/1471-2105-15-293 [PubMed]
  • 26. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE. The Protein Data Bank. Nucleic Acids Res. 2000; 28:235–42. https://doi.org/10.1093/nar/28.1.235 [PubMed]
  • 27. Chan KK, Yao TJ, Jones B, Zhao JF, Ma FK, Leung CY, Lau SK, Yip MW, Ngan HY. The use of Chinese herbal medicine to improve quality of life in women undergoing chemotherapy for ovarian cancer: a double-blind placebo-controlled randomized trial with immunological monitoring. Ann Oncol. 2011; 22:2241–9. https://doi.org/10.1093/annonc/mdq749 [PubMed]
  • 28. Parekh HS, Liu G, Wei MQ. A new dawn for the use of traditional Chinese medicine in cancer therapy. Mol Cancer. 2009; 8:21. https://doi.org/10.1186/1476-4598-8-21 [PubMed]
  • 29. Colunga Biancatelli RML, Berrill M, Catravas JD, Marik PE. Quercetin and Vitamin C: An Experimental, Synergistic Therapy for the Prevention and Treatment of SARS-CoV-2 Related Disease (COVID-19). Front Immunol. 2020; 11:1451. https://doi.org/10.3389/fimmu.2020.01451 [PubMed]
  • 30. Russo M, Spagnuolo C, Tedesco I, Bilotto S, Russo GL. The flavonoid quercetin in disease prevention and therapy: facts and fancies. Biochem Pharmacol. 2012; 83:6–15. https://doi.org/10.1016/j.bcp.2011.08.010 [PubMed]
  • 31. Massi A, Bortolini O, Ragno D, Bernardi T, Sacchetti G, Tacchini M, De Risi C. Research Progress in the Modification of Quercetin Leading to Anticancer Agents. Molecules. 2017; 22:1270. https://doi.org/10.3390/molecules22081270 [PubMed]
  • 32. Reyes-Farias M, Carrasco-Pozo C. The Anti-Cancer Effect of Quercetin: Molecular Implications in Cancer Metabolism. Int J Mol Sci. 2019; 20:3177. https://doi.org/10.3390/ijms20133177 [PubMed]
  • 33. Tang SM, Deng XT, Zhou J, Li QP, Ge XX, Miao L. Pharmacological basis and new insights of quercetin action in respect to its anti-cancer effects. Biomed Pharmacother. 2020; 121:109604. https://doi.org/10.1016/j.biopha.2019.109604 [PubMed]
  • 34. Sharma H, Sen S, Singh N. Molecular pathways in the chemosensitization of cisplatin by quercetin in human head and neck cancer. Cancer Biol Ther. 2005; 4:949–55. https://doi.org/10.4161/cbt.4.9.1908 [PubMed]
  • 35. de Paula Rodrigues R, Tini IR, Soares CP, da Silva NS. Effect of photodynamic therapy supplemented with quercetin in HEp-2 cells. Cell Biol Int. 2014; 38:716–22. https://doi.org/10.1002/cbin.10251 [PubMed]
  • 36. Bądziul D, Jakubowicz-Gil J, Paduch R, Głowniak K, Gawron A. Combined treatment with quercetin and imperatorin as a potent strategy for killing HeLa and Hep-2 cells. Mol Cell Biochem. 2014; 392:213–27. https://doi.org/10.1007/s11010-014-2032-4 [PubMed]
  • 37. Lane D, Levine A. p53 Research: the past thirty years and the next thirty years. Cold Spring Harb Perspect Biol. 2010; 2:a000893. https://doi.org/10.1101/cshperspect.a000893 [PubMed]
  • 38. Zhou G, Liu Z, Myers JN. TP53 Mutations in Head and Neck Squamous Cell Carcinoma and Their Impact on Disease Progression and Treatment Response. J Cell Biochem. 2016; 117:2682–92. https://doi.org/10.1002/jcb.25592 [PubMed]
  • 39. Poeta ML, Manola J, Goldenberg D, Forastiere A, Califano JA, Ridge JA, Goodwin J, Kenady D, Saunders J, Westra W, Sidransky D, Koch WM. The Ligamp TP53 Assay for Detection of Minimal Residual Disease in Head and Neck Squamous Cell Carcinoma Surgical Margins. Clin Cancer Res. 2009; 15:7658–65. https://doi.org/10.1158/1078-0432.CCR-09-1433 [PubMed]
  • 40. Clemente-Soto AF, Salas-Vidal E, Milan-Pacheco C, Sánchez-Carranza JN, Peralta-Zaragoza O, González-Maya L. Quercetin induces G2 phase arrest and apoptosis with the activation of p53 in an E6 expression-independent manner in HPV-positive human cervical cancer-derived cells. Mol Med Rep. 2019; 19:2097–106. https://doi.org/10.3892/mmr.2019.9850 [PubMed]
  • 41. Kozaki K, Imoto I, Pimkhaokham A, Hasegawa S, Tsuda H, Omura K, Inazawa J. PIK3CA mutation is an oncogenic aberration at advanced stages of oral squamous cell carcinoma. Cancer Sci. 2006; 97:1351–8. https://doi.org/10.1111/j.1349-7006.2006.00343.x [PubMed]
  • 42. Chan CY, Lien CH, Lee MF, Huang CY. Quercetin suppresses cellular migration and invasion in human head and neck squamous cell carcinoma (HNSCC). Biomedicine (Taipei). 2016; 6:15. https://doi.org/10.7603/s40681-016-0015-3 [PubMed]
  • 43. Chung CH, Ely K, McGavran L, Varella-Garcia M, Parker J, Parker N, Jarrett C, Carter J, Murphy BA, Netterville J, Burkey BB, Sinard R, Cmelak A, et al. Increased epidermal growth factor receptor gene copy number is associated with poor prognosis in head and neck squamous cell carcinomas. J Clin Oncol. 2006; 24:4170–6. https://doi.org/10.1200/JCO.2006.07.2587 [PubMed]
  • 44. Huang YH, Chang CY, Kuo YZ, Fang WY, Kao HY, Tsai ST, Wu LW. Cancer-associated fibroblast-derived interleukin-1β activates protumor C-C motif chemokine ligand 22 signaling in head and neck cancer. Cancer Sci. 2019; 110:2783–93. https://doi.org/10.1111/cas.14135 [PubMed]
  • 45. Chen F, Zheng A, Li F, Wen S, Chen S, Tao Z. Screening and identification of potential target genes in head and neck cancer using bioinformatics analysis. Oncol Lett. 2019; 18:2955–66. https://doi.org/10.3892/ol.2019.10616 [PubMed]
  • 46. Oliveira-Paula GH, Lacchini R, Tanus-Santos JE. Endothelial nitric oxide synthase: From biochemistry and gene structure to clinical implications of NOS3 polymorphisms. Gene. 2016; 575:584–99. https://doi.org/10.1016/j.gene.2015.09.061 [PubMed]
  • 47. Cooke GE, Doshi A, Binkley PF. Endothelial nitric oxide synthase gene: prospects for treatment of heart disease. Pharmacogenomics. 2007; 8:1723–34. https://doi.org/10.2217/14622416.8.12.1723 [PubMed]
  • 48. Choi JY, Barlow WE, Albain KS, Hong CC, Blanco JG, Livingston RB, Davis W, Rae JM, Yeh IT, Hutchins LF, Ravdin PM, Martino S, Lyss AP, et al. Nitric oxide synthase variants and disease-free survival among treated and untreated breast cancer patients in a Southwest Oncology Group clinical trial. Clin Cancer Res. 2009; 15:5258–66. https://doi.org/10.1158/1078-0432.CCR-09-0685 [PubMed]
  • 49. Marisi G, Petracci E, Raimondi F, Faloppi L, Foschi FG, Lauletta G, Iavarone M, Canale M, Valgiusti M, Neri LM, Ulivi P, Orsi G, Rovesti G, et al. ANGPT2 and NOS3 Polymorphisms and Clinical Outcome in Advanced Hepatocellular Carcinoma Patients Receiving Sorafenib. Cancers (Basel). 2019; 11:1023. https://doi.org/10.3390/cancers11071023 [PubMed]
  • 50. Sonveaux P, Brouet A, Havaux X, Grégoire V, Dessy C, Balligand JL, Feron O. Irradiation-induced angiogenesis through the up-regulation of the nitric oxide pathway: implications for tumor radiotherapy. Cancer Res. 2003; 63:1012–9. [PubMed]
  • 51. Yanar K, Çakatay U, Aydın S, Verim A, Atukeren P, Özkan NE, Karatoprak K, Cebe T, Turan S, Ozkök E, Korkmaz G, Cacına C, Küçükhüseyin O, Yaylım İ. Relation between Endothelial Nitric Oxide Synthase Genotypes and Oxidative Stress Markers in Larynx Cancer. Oxid Med Cell Longev. 2016; 2016:4985063. https://doi.org/10.1155/2016/4985063 [PubMed]
  • 52. Guo MF, Dai YJ, Gao JR, Chen PJ. Uncovering the Mechanism of Astragalus membranaceus in the Treatment of Diabetic Nephropathy Based on Network Pharmacology. J Diabetes Res. 2020; 2020:5947304. https://doi.org/10.1155/2020/5947304 [PubMed]
  • 53. Sohn TA, Bansal R, Su GH, Murphy KM, Kern SE. High-throughput measurement of the Tp53 response to anticancer drugs and random compounds using a stably integrated Tp53-responsive luciferase reporter. Carcinogenesis. 2002; 23:949–57. https://doi.org/10.1093/carcin/23.6.949 [PubMed]
  • 54. Psyrri A, Yu Z, Weinberger PM, Sasaki C, Haffty B, Camp R, Rimm D, Burtness BA. Quantitative determination of nuclear and cytoplasmic epidermal growth factor receptor expression in oropharyngeal squamous cell cancer by using automated quantitative analysis. Clin Cancer Res. 2005; 11:5856–62. https://doi.org/10.1158/1078-0432.CCR-05-0420 [PubMed]
  • 55. Martini M, De Santis MC, Braccini L, Gulluni F, Hirsch E. PI3K/AKT signaling pathway and cancer: an updated review. Ann Med. 2014; 46:372–83. https://doi.org/10.3109/07853890.2014.912836 [PubMed]
  • 56. Xu W, Yang Z, Lu N. A new role for the PI3K/Akt signaling pathway in the epithelial-mesenchymal transition. Cell Adh Migr. 2015; 9:317–24. https://doi.org/10.1080/19336918.2015.1016686 [PubMed]
  • 57. Tang T, Xiao ZY, Shan G, Lei HB. Descending-SHIP2-mediated radiosensitivity enhancement through PI3K/Akt signaling pathway in laryngeal squamous cell carcinoma. Biomed Pharmacother. 2019; 118:109392. https://doi.org/10.1016/j.biopha.2019.109392 [PubMed]
  • 58. Ye D, Zhou C, Deng H, Lin L, Zhou S. MicroRNA-145 inhibits growth of laryngeal squamous cell carcinoma by targeting the PI3K/Akt signaling pathway. Cancer Manag Res. 2019; 11:3801–12. https://doi.org/10.2147/CMAR.S199291 [PubMed]