DOI QR코드

DOI QR Code

Discovering the Anti-cancer Effects of Ligusticum Chuanxiong through Network-based Pharmacology Analysis and Molecular Docking: An Inquiry into Natural Products

네트워크 기반 약리학 분석 및 분자 도킹을 통한 천궁의 항암 효과 예측: 천연물에 대한 탐구

  • Do Kyung Han (Department of Pharmacology, School of Korean Medicine, Pusan National University) ;
  • Jee Won Shon (Department of Pharmacology, School of Korean Medicine, Pusan National University) ;
  • Eui Suk Sung (Department of Otolaryngology-Head and Neck Surgery, Research Institute for Convergence of Biomedical Science and Technology, Pusan National University, Yangsan Hospital) ;
  • Youn Sook Kim (Research Institute for Longevity and Well-Being, Pusan National University) ;
  • Won G. An (Department of Pharmacology, School of Korean Medicine, Pusan National University)
  • 한도경 (부산대학교 한의학전문대학원) ;
  • 손지원 (부산대학교 한의학전문대학원) ;
  • 성의숙 (양산부산대학교병원 이비인후과) ;
  • 김윤숙 (부산대학교 장수웰빙연구소) ;
  • 안원근 (부산대학교 한의학전문대학원)
  • Received : 2023.10.30
  • Accepted : 2023.11.20
  • Published : 2023.11.30

Abstract

In some cases of head and neck cancers (HNC), surgical interventions may result in the loss of organs and/or changes to their functions, thereby significantly affecting the patient's quality of life. As a result, the surgical treatment of HNC patients is often limited to specific cases, and alternative treatment modalities, such as chemotherapy, are considered. However, serious adverse effects caused by chemotherapy, such as severe nausea and vomiting, necessitate the need for the development of adjunctive methods to minimize patient suffering. Chuanxiong, Ligusticum chuanxiong (L. chuanxiong), is a natural herb used in Eastern medicine to treat cerebrovascular disorders and headaches. This study aimed to predict the effect and potential of L. chuanxiong as an auxiliary anticancer drug through network-based pharmacology and molecular docking analysis. The study results showed that 40 out of 41 genes of L. chuanxiong shared common targets of HNC and their proteins could be used to target HNC cells to prevent cancer progression. The results of the functional enrichment analysis confirmed that L. chuanxiong is associated with the neuroactive-ligand metabolism and neurotransmitter pathways, indicating its potential medicinal value as an adjuvant in HNC treatment. Lastly, our findings demonstrated that the active ingredient of L. chuanxiong, (Z)-Ligustilide, has the ATP binding site of heat shock protein 90, a protein known to promote the activation of cancer cells. These results suggest that L. chuanxiong is a promising candidate for developing auxiliary anticancer drugs, and further research could potentially lead to the discovery of newer and safer anti-cancer agents.

두경부암(HNC)의 경우, 외과적 개입은 환자의 삶의 질에 심각한 영향을 미칠 수 있으며, 화학요법을 병행하게 된다. 그러나 화학요법에는 현저한 부작용이 있으므로 환자의 고통을 최소화하기 위한 보조 방법의 개발이 필요하다. 천궁(Ligusticum chuanxiong)은 동양 의학에서 뇌혈관 장애 및 두통에 사용되는 천연 허브이다. 본 연구에서는 네트워크 기반 약리학 및 분자 도킹 분석을 통해 천궁의 근본적인 항암기전을 예측하였다. 본 연구에서 HNC와 관련된 천궁의 공통 유전자를 밝혀내어 신경 활성 리간드의 대사 및 신경 전달 물질 경로와의 연관성을 확인했다. 본 연구는 천궁의 성분 중 하나인 (Z)-ligustilide 가 암세포 활성화에 관련된 heat shock protein 90의 ATP 결합 부위를 공유함을 입증했다. 이 결과는 천궁이 보조 항암제 개발을 위한 유망한 후보임을 시사하며, 향 후 더욱 새롭고 안전한 항암제의 연구개발에 과학적 근거를 제시하는 새로운 발견이다.

Keywords

Introduction

Head and neck cancer (HNC) is the sixth most frequent cancer in the world, according to Globocan 2020, an online database of the International Agency for Research on Cancer (IARC) that provides global cancer statistics [9]. The human papillomavirus (HPV), immunodeficiency, and exposure to carcinogens, such as those involved in heavy smoking and alcohol consumption are the well-known factors that raise the risk of HNC [29]. Among the cases of HNC, 90% are head and neck squamous cell carcinomas (HNSCC), which has the eighth highest mortality rate among other types of cancer.

HNC located in the oral cavity or oropharynx is visually observable with the naked eye. Further, examination of the pharynx and larynx through nasofibroscopy allows real-time observation and ease in performing biopsies. Despite this, HNC is usually diagnosed only in loco-regionally advanced stages compared to other solid tumors, and more than 50% of the patients die within five years of being diagnosed, making it one of the leading causes of death due to cancer worldwide [15].

Modern conventional treatment options for HNC include surgical intervention, radiotherapy, and chemotherapy. However, in some cases of HNC, the use of surgical interventions may be limited in order to preserve the patient’s quality of life. Therefore, radiotherapy and chemotherapy are often the preferred alternative treatment modalities. According to the National Comprehensive Cancer Network (NCCN), combining cisplatin and radiation therapy is highly effective and recommended for the treatment of HNC patients. NCCN also suggests the co-administration of cetuximab, an effective radiation sensitizer with radiotherapy, and the weekly administration of carboplatin, an alkylating agent similar to cisplatin [4]. However, the treatments mentioned above induce severe side effects: patients commonly experience dry mouth, anxiety, vomiting, and chronic fatigue. Besides, unlike other cancer patients, HNC patients may also experience difficulty breathing, eating, and speaking, depending on the anatomical location of the primary cancer [6, 11, 16]. The poor clinical outcomes and the high toxicity of the chemotherapeutic agents necessitate the need to develop less toxic therapeutic agents for HNC patients. This comprehensive study explores the medicinal effects and biological mechanisms of the natural product Chuanxiong, Ligusticum chuanxiong (L. chuanxiong), as a safe and efficacious alternative to currently available anti-cancer therapies.

For thousands of years, L. chuanxiong has been utilized to manage head and neck pain in traditional medicine. It is also the most frequently mentioned herb in the traditional Korean medical book, Donguibogam, for treating diseases associated with the cranial region [7, 26]. In Eastern medicine, this herb is commonly prescribed for treating virus-related diseases and various inflammatory reactions in the head [34]. Recent studies have reported that L. chuanxiong exerts its anticancer effects by regulating the oxidative stress pathways and inducing apoptosis in cancer cells [14]. However, its efficacy and safety as an anticancer agent is unclear, and further research is required.

This study predicted and validated the anticancer effects of (Z)-Ligustilide, an active ingredient of L. chuanxiong through systematic pharmacological analyses and molecular docking. The findings of this study suggest that this ingredient could contribute to modern anticancer therapeutics activity through the development of an agent which could function alone or synergistically with other chemotherapeutic agents through an alternate mechanism of action offering safety and efficacy for the treatment of HNC.

Materials and Methods

Screening active compounds and fishing the target genes of L. chuanxiong

To identify the major bioactive compounds in L. chuanxiong we used the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP, https://tcmsp-e.com/tcmsp.php) [22]. TCMSP provides the pharmacokinetic properties of each component of the herbal medicines along with information on the correlation between compounds, targets, and diseases. In this study, molecular weight (MW), oral bioavailability (OB), and drug-likeness (DL) values were referenced to identify the active compounds of L. chuanxiong [25, 36]. Compounds that satisfied the conditions of 180≤ MW ≤500, OB ≥50%, and DL ≥0.07 were considered the active compounds of L. chuanxiong. The predicted target proteins of these selected active compounds were also obtained from the TCMSP database, and the corresponding genes translating the target proteins were selected from the UniProt Knowledgebase (https://www.uniprot.org/) [32]. The compound-target network was constructed using Cytoscape 3.9.1. (https://cytoscape.org/).

Common target genes of HNC and L. chuanxiong

The target genes for the treatment of HNC were selected using the DisGeNET (https://www.disgenet.org/) [20, 21], GeneCards database (https://www.genecards.org/) [23, 28] and Comparative Toxicogenomics Database (CTD, https://ctdb ase.org/) [5]. A search was conducted on each database using keywords such as 'Head and neck cancer', 'Head and neck neoplasms' and 'Head and neck carcinoma'. Using the FunRich 3.1.4 software tool (http://www.funrich.org) [10, 18, 19], a Venn diagram was drawn to identify genes that were duplicated twice or more among all the genes found in these three databases. From this selection, the genes that intersected with the genes encoding the targets of active compounds of L. chuanxiong were chosen as the common target genes.

Protein-protein interaction network construction

To investigate the Protein-Protein Interaction (PPI) between the Common target Genes of L. chuanxiong and HNC (CGLHs), we utilized the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database (https://stringdb.org/). We set the species as Homo sapiens and by inputting CGLHs into the database under these criteria, we sought to investigate their PPIs.

Functional enrichment analysis

The functional analysis of the intersecting genes between HNC and L. chuanxiong was conducted using Gene ontology (GO) and the Kyoto encyclopedia of genes and genomes (KEGG) pathways. We utilized the database for annotation, visualization, and integrated discovery (DAVID) bioinformatics resources for mapping and visualized the GO analysis using the Science and research online plot (SRplot, https:// www.bioinformatics.com.cn/en). The genes were analyzed at the biological process (BP), cellular component (CC), and molecular function (MF) levels to verify their functions.

Molecular docking

Docking analysis of molecular-level interactions between ligands and proteins was performed using AutoDock Vina (https://github.com/ccsb-scripps/AutoDock-Vina) [8, 31]. The protein data bank (PDB) structures of (Z)-Ligustilide (Pubchem CID 5319022) and adenosine triphosphate (ATP, Pubchem CID 5957) were downloaded from the Pubchem compound database (https://pubchem.ncbi.nlm.nih.gov). The PDB ID of human heat shock protein 90 (HSP90), 1YES (resolution 2.20 Å) [27], was obtained from the Research Collaboratory for Structural Bioinformatics (RCSB) Protein Data Bank website (https://www.rcsb.org/). Converting the PDB files to Protein Data Bank, Partial Charge (Q), & Atom Type (T) (PDBQT) files and creating grid boxes for docking was done using the AutoDockTool-1.5.7. The visualization of the results of the protein-ligand binding was performed using the Biovia Discovery Studio Visualizer v17.2, 2021 (BIOVIA, Dssault Systèmes, Waltham, USA).

Results

The active ingredients of L. chuanxiong

A total of 22 active compounds of L. chuanxiong were identified by setting the criteria of 180≤ MW ≤500, OB ≥ 50%, and DL ≥0.07 in the TCMSP database (Table 1).

Table 1. Active compounds of L. chuanxiong with 180≤ MW ≤500, OB ≥50%, and DL ≥0.07 in the TCMSP database

SMGHBM_2023_v33n11_876_t0001.png 이미지

Identification of disease-drug common target genes

A total of 41 target proteins of the active compounds of L. chuanxiong were identified and converted into gene IDs using the UniProt database (Table 2). The network of the active compounds and target genes are visualized in Fig. 1. The average number of neighbors within this network was 6.476 for each compound. Three compounds had the highest degree score and are assumed to play a central role: Methyl 2-pentanoylbenzoate (22 targets), 4-Hydroxy-3-butylphthalide (21 targets), and (Z)-Ligustilide (18 targets). The centrality rating of targets, by degree score, are ranked as follows: Gamma-aminobutyric acid receptor subunit alpha-1 (GABRA1), Prostaglandin G/H synthase 2 (PTGS2), Muscarinic acetylcholine receptor M1, 2, 3 (CHRM1, CHRM2, CHRM3).

Table 2. 41 Targets of Ligusticum chuanxiong

SMGHBM_2023_v33n11_876_t0002.png 이미지

SMGHBM_2023_v33n11_876_f0001.png 이미지

Fig. 1. Network of active compounds of L. chuanxiong and their targets. The network shows the multiple targeting characteristics of L. chuanxiong. Blue squares represent active compounds, and red circles represent the target genes. The size of the shape of nodes indicates the degree score.

To find the CGLHs, we searched the DisGeNET, Gene Cards, and CTD databases to collect genes associated with HNC. We obtained 786 genes from DisGeNET, 1,418 genes from GeneCards, and 35,682 genes from CTD. From these, we selected genes that overlapped across multiple databases and were also related to L. chuanxiong (Fig. 2). Out of the 41 L. chuanxiong-related genes, 40 genes (except for Dopamine transporter1 [DAT1]) overlapped with the target genes of HNC.

SMGHBM_2023_v33n11_876_f0002.png 이미지

Fig. 2. Identification of the common target genes of L. chuanxiong and head and neck cancer. Venn Intersection Analysis by FunRich3.14 showed there were 40 target genes common to L. chuanxiong and head and neck cancer.

PPI network construction

We confirmed that there were 40 CGLHs and put these 40 targets into the STRING database to proceed with the PPI network analysis. In the PPI network, 40 nodes and 166 edges were observed. The thickness of the network edges represents the confidence score from data sources. The active interaction source categories settings for the confidence score mainly focused on the experimental/biochemical data, association in curated databases, and textmining evidence through mentioning in relevant publications (Fig. 3A). As a result of calculating the importance of genes through the bottleNeck of cytohubba plugin in cytoscpe, PTGS2, GABRA1, Estrogen receptor (ESR1), Dopamine D1 receptor (DRD1), GABRA2, PTGS1, Sodium-dependent serotonin transporter (SLC6A4), Choline O-Acetyltransferase (CHAT), Beta-1 adrenergic receptor (ADRB2), and HSP90AB1 were in the top 10. Intermediation represents the number of the shortest path passing through a node, and nodes with such high mediation nodes are called BottleNecks and are one of the conditions for deriving a hub gene along with the degree (Fig. 3B) [37].

SMGHBM_2023_v33n11_876_f0003.png 이미지

Fig. 3. Protein-protein interaction network construction for common genes between L. chuanxiong and head and neck cancer. (A) There were 40 nodes and 166 edges found. The PPI enrichment p-value was less than 1.0e-16. The average node degree was 8.9. The proteins with notably high concentration scores are highlighted in red. (B) Key genes generated from cytoHubba plugin in Cytoscape ranked by BottleNeck.

Enrichment analysis for HNC - L. chuanxiong common targets

Using the 40 CGLHs, GO and KEGG pathway analyses were conducted. According to the GO analysis, 1,409 results were obtained in the perspective of BP, 142 results in the perspective of CC, and 200 results in the perspective of MF. The most observed terms in all three perspectives were dopamine, gamma-aminobutyric acid (GABA) and synaptic transmission. This indicated that the active compounds in L. chuanxiong were associated with the treatment of diseases associated with the head and neck region (Fig. 4A).

SMGHBM_2023_v33n11_876_f0004.png 이미지

Fig. 4. Enrichment analysis of 40 common targets. (A) Gene Ontology (GO) enrichment analysis at the level of Biological Process (BP), Cellular Component (CC), and Molecular Function (MF). (B) Kyoto Encyclopaedia of Genes and Genomes (KEGG) enrichment pathways.

In the KEGG pathway analysis, pathways related to neuro-active ligand-receptor interaction were found to be significant. This indicates that the CGLHs are closely related to the signalling pathway of diseases associated with the human brain (Fig. 4B).

Molecular docking analysis for (Z)-Ligustilide and HSP 90

To validate the molecular binding between one of the active compounds of L. chuanxiong, (Z)-Ligustilide and the molecular chaperone HSP90, we employed the PDB code 1YES as the receptor. 1YES is HSP90 N-terminal domain and reported as human HSP90 geldanamycin binding domain (HSP90-GBD), open conformation (residues 9-232). Geldanamycin mimics ATP structure binding in the HSP90 N-terminal nucleotide-binding pocket, and the pocket is made up by the residues such as Leucine (Leu)48, Asparagine (Asn)51, Aspartic acid (Asp)54, Alanine (Ala)55, Lysine (Lys)58, Isoleucine (Ile)91, Asp93, Ile96, Glycine (Gly)97, Methionine (Met)98, Asn106, Leu107, Lys112, Gly135, Phenylalanine (Phe)138, Valine (Val)150, Threonine (Thr)184, and Val186 [27]. The human HSP90 N-terminal ATP-binding sites are reported in the crystal structure of HSP90 N-terminal-ATP complex (PDB code 3T0Z), and the interacting residues are Asn51, Asp54, Asp93, Glu47, Gly135, Gly137, Gly97, Leu48, Met98, Phe138, Ser52, Thr184, and Val136 [13]. Molecular docking results showed a calculated ligand-receptor binding affinity of -6.1 kcal/mol, with the formation of conventional hydrogen bonds, hydrophobic bonds and van der Waals contacts between (Z)-Ligustilide and amino acid residues of 1YES (Fig. 5A). (Z)-Ligustilide interacted with the following amino acid residues: Asp93, Ala55, Asn51, Met98, Gly97, Ile96, Phe138, Leu107, Thr184, Ser52, Asn106, Val150, and Val186. O2 of (Z)-Ligustilide made a conventional hydrogen bond to Gly97. (Z)-Ligustilide made hydrophobic bindings to Asn51, Ala55, Met98, and Val 186, and made van der Waals contacts with Phe138, ASP93, Ser52, Thr184, Ile96, Asn106, Leu107, and Val150 (Fig. 5B). To further compare the (Z)-Ligustilide binding sites and ATP binding sites for Hsp90, we performed molecular docking between ATP and 1YES (Fig. 5C), and the interacting modes are visualized in Fig. 5D.

SMGHBM_2023_v33n11_876_f0005.png 이미지

Fig. 5. Molecular docking analysis between HSP90 (PDB code 1YES) and (Z)-Ligustilide and ATP. (A) Molecular docking of (Z)-Ligustilide to 1YES model. (B) Two-dimensional binding diagram of (Z)-Ligustilide and amino acid residues of 1YES. (C) Molecular docking of ATP to 1YES. (D) Two-dimensional visualization of ATP-1YET docking result.

Using the alignments of HSP-GBD sequence, the binding sites of (Z)-Ligustilide and ATP were compared. Thirteen amino acid residues of (Z)-Ligustilide-HSP90 binding site were in the nucleotide-binding pocket of the HSP90 N-terminal domain. Eleven residues were the same as the ATP-binding site of the ATP-1YES binding result. Eight residues were the same as the ATP-binding site of 3T0Z (Fig. 6).

SMGHBM_2023_v33n11_876_f0006.png 이미지

Fig. 6. (Z)-Ligustilide and HSP90 binding site analysis using alignment of HSP90 N-terminal domain sequences (1YES). The amino acid sequences of 1YES were downloaded from UnitProt database. The residues that make up the nucleotide binding pocket in 1YES [27], the ATP-binding sites of 3T0Z [13], the ATP-binding sites of 1YES, and the (Z)-Ligustilide binding sites of 1YES are highlighted (cyan blue).

Lastly, the results of the compound-target network analysis showed that (Z)-Ligustilide has 18 targets (Fig. 1). According to our molecular docking results between (Z)-Ligustilide and HSP90, HSP90 should be added as a target of (Z)-Ligustilide.

Discussion

The 22 active compounds of L. chuanxiong were identified from the network pharmacology analyses (Table 1). The relevant genes of the compounds were compared with the HNC genes to identify the common genes. It was found that 40 out of the 41 genes relating to L. chuanxiong overlapped with the HNC–associated genes (Fig. 2). The key genes generated from PPIs of CGLHs were PTGS2, GABRA1, ESR1, DRD1, GABRA2, PTGS1, SLC6A4, CHAT, ADRB2, and HSP90AB1 (Fig. 3B). We conducted GO and KEGG enrichment analyses on the CGLHs, further exploring the underlying biological mechanisms of L. chuanxiong. We approached the studies with three main areas of focus: BP, CC and MF. The GO results showed that the target genes mainly enriched the biological functions that regulate the postsynaptic membrane potential, synaptic transmission GABAnergic, and the regulation of the membrane potential (Fig. 4A). The results support the correlation between L. chuanxiong and neurotransmitter transmission – especially the inhibitory signalling pathways. Further research could discover valuable medicinal insights into L. chuanxiong to provide a supporting mechanism to reduce the unwanted side effects of modern conventional anti-cancer therapies.

Notably, this study further elucidated the anti-cancer mechanism of L. chuanxiong in treating HNC through (Z)-Ligustilide, one of its main active compounds. We focused on its correlation with HSP90 which is an ATP-dependent protein that indirectly affects the processing of cell cycle regulation, proliferation, migration, and apoptosis by folding various regulatory proteins in response to external stress on the cell [30]. Among the isoforms of HSP90, HSP90AB1 was highly expressed in HNSC, regulating HNSC cell proliferation, migration, and glycolysis through the AKT pathway [38]. HSP90 acts as a chaperone, assisting in proper protein folding and maintaining the structure of proteins [24]. It controls the cell cycle and aids in DNA repair in response to damage to the DNA [2]. In addition, HSP90 can also stabilize overexpressed oncogenic proteins that contribute to cancers [17]. In some cases, it also acts as a factor that promotes cancer development. These characteristics have made HSP90 a potential target for cancer therapy [33, 35].

Disrupting the structure of HSP90 complexes, which directly aid in promoting oncoproteins, should facilitate the degradation of target proteins. This is also an important factor that validates the function of anti-cancer drugs [1, 3, 12]. Studies have reported that when Baicalein, an HSP90 ATPase inhibitor, binds to the ATP binding site of HSP90, the HSP90/ PTGS2 complex dissociates, leading to the inactivation of HSP90 [39]. Here, when HSP90 activates as it binds with PTGS2, to deactivate HSP90, a drug should work as an ATPase domain inhibitor of HSP90. In this study, (Z)-Ligustilide bound to the HSP90 chaperone, interfering with the function of the signalling molecules. We performed 3D molecular docking using AutoDock Vina, generated a 2D visualization model, and verified the association by examining (Z)-Ligustilide and HSP90 (Fig. 5A, 5B). The ATP binding sites of HSP90 and (Z)-Ligustilide were found to be identical based on sequences comparison with an existing ATP binding model (3T0Z) and molecular docking results of the ATP binding site of 1YES (Fig. 6). This suggests that the underlying anti-cancer mechanisms of L. chuanxiong can be further elucidated and future research on this natural product could lead to advancement in anti-cancer therapeutics. However, these results are limited to laboratory-based systematic pharmacology. Further practical validation is needed to precisely analyse the underlying molecular mechanisms for use in nature-based adjuvant therapies or develop a new method in conventional therapeutics with fewer side effects.

We provided insights into the medicinal value of L. chuanxiong. Based on system pharmacology, we have suggested that the therapeutic role of the herb in the treatment of cerebrovascular diseases is an important factor that could lead to advancements in the treatment of other diseases including modern anticancer therapeutics by reducing the serious side effects associated with chemotherapeutic agents. We also elucidated the anti-cancer mechanisms of (Z)-Ligustilide, one of the main active compounds of L. chuanxiong and we confirmed that (Z)-Ligustilide binds to the ATP binding site of HSP90. Therefore, we suggest that further research may unravel the role of L. chuanxiong as a potential agent in adjuvant therapy for the treatment of cancer.

Acknowledgements

This research was supported by [Basic Science Research Program] through the [National Research Foundation of Korea (NRF)] funded by the [Ministry of Education] under grant [NRF-2021R1I1A3A04037158] and [NRF-2021R1I1A 1A01058697].

The Conflict of Interest Statement

The authors declare that they have no conflicts of interest with the contents of this article.

References

  1. An, W.G., Schulte, T.W. and Neckers, L.M. 2000. The heat shock protein 90 antagonist geldanamycin alters chaperone association with p210bcr-abl and v-src proteins before their degradation by the proteasome. Cell Growth Differ. 11, 355-360.
  2. Arlander, S.J., Felts, S.J., Wagner, J.M., Stensgard, B., Toft, D.O. and Karnitz, L.M. 2006. Chaperoning checkpoint kinase 1 (Chk1), an Hsp90 client, with purified chaperones. J. Biol. Chem. 281, 2989-2998. https://doi.org/10.1074/jbc.M508687200
  3. Canella, A., Welker, A.M., Yoo, J.Y., Xu, J., Abas, F.S., Kesanakurti, D., Nagarajan, P., Beattie, C.E., Sulman, E.P. and Liu, J. 2017. Efficacy of onalespib, a long-acting second-generation HSP90 inhibitor, as a single agent and in combination with temozolomide against malignant gliomas. Clin. Cancer Res. 23, 6215-6226. https://doi.org/10.1158/1078-0432.CCR-16-3151
  4. Caudell, J.J., Gillison, M.L., Maghami, E., Spencer, S., Pfister, D.G., Adkins, D., Birkeland, A.C., Brizel, D.M., Busse, P.M. and Cmelak, A.J., et al. 2022. NCCN Guidelines® Insights: Head and neck cancers, version 1.2022. J. Natl. Compr. Canc. Netw. 20, 224-234. https://doi.org/10.6004/jnccn.2022.0016
  5. Davis, A.P., Wiegers, T.C., Johnson, R.J., Sciaky, D., Wiegers, J. and Mattingly, C.J. 2023. Comparative toxicogenomics database (CTD): update 2023. Nucleic Acids Res. 51, D1257-D1262. https://doi.org/10.1093/nar/gkac833
  6. Deganello, A., Battat, N., Muratori, E., Cristofaro, G., Buongiorno, A., Mannelli, G., Picconi, M., Giachetti, R., Borsotti, G. and Gallo, O. 2016. Acupuncture in shoulder pain and functional impairment after neck dissection: A prospective randomized pilot study. Laryngoscope. 126, 1790-1795. https://doi.org/10.1002/lary.25921
  7. Dong yi bao jian. 1613. Seoul: Royal Hospital.
  8. Eberhardt, J., Santos-Martins, D., Tillack, A.F. and Forli, S. 2021. Autodock vina 1.2. 0: New docking methods, expanded force field, and python bindings. J. Chem. Inf. Model. 61, 3891-3898. https://doi.org/10.1021/acs.jcim.1c00203
  9. Fang, S., Dong, L., Liu, L., Guo, J., Zhao, L., Zhang, J., Bu, D., Liu, X., Huo, P. and Cao, W., et al. 2021. HERB: a high-throughput experiment- and reference-guided database of traditional Chinese medicine. Nucleic Acids Res. 49, D1197-D1206. https://doi.org/10.1093/nar/gkaa1063
  10. Fonseka, P., Pathan, M., Chitti, S.V., Kang, T. and Mathivanan, S. 2021. FunRich enables enrichment analysis of OMICs datasets. J. Mol. Biol. 433, 166747.
  11. Goodman, J.F. and Wang, M.B. 2022. Complementary and integrative medicine in head and neck cancer. Otolaryngol Clin North Am. 55, 993-1006. https://doi.org/10.1016/j.otc.2022.06.007
  12. Kryeziu, K., Bruun, J., Guren, T.K., Sveen, A. and Lothe, R.A. 2019. Combination therapies with HSP90 inhibitors against colorectal cancer. Biochim. Biophysica. Acta. Rev. Cancer. 1871, 240-247. https://doi.org/10.1016/j.bbcan.2019.01.002
  13. Li, J., Sun, L., Xu, C., Yu, F., Zhou, H., Zhao, Y., Zhang, J., Cai, J., Mao, C. and Tang, L. 2012. Structure insights into mechanisms of ATP hydrolysis and the activation of human heat-shock protein 90. Acta. Biochim. Biophys. Sin. 44, 300-306. https://doi.org/10.1093/abbs/gms001
  14. Li, W., Tang, Y., Chen, Y. and Duan, J.A. 2012. Advances in the chemical analysis and biological activities of chuanxiong. Molecules. 17, 10614-10651. https://doi.org/10.3390/molecules170910614
  15. Li, Y., Li, S., Meng, X., Gan, R.Y., Zhang, J.J. and Li, H.B. 2017. Dietary natural products for prevention and treatment of breast cancer. Nutrients. 9, 728.
  16. Margalit, D.N., Salz, T., Venchiarutti, R., Milley, K., McNamara, M., Chima, S., Wong, J., Druce, P. and Nekhlyudov, L. 2022. Interventions for head and neck cancer survivors: Systematic review. Head Neck. 44, 2579-2599. https://doi.org/10.1002/hed.27142
  17. Neckers, L. and Ivy, S.P. 2003. Heat shock protein 90. Curr. Opin. Oncol. 15, 419-424. https://doi.org/10.1097/00001622-200311000-00003
  18. Pathan, M., Keerthikumar, S., Ang, C.S., Gangoda, L., Quek, C.Y., Williamson, N.A., Mouradov, D., Sieber, O.M., Simpson, R.J. and Salim, A. 2015. FunRich: An open access standalone functional enrichment and interaction network analysis tool. Proteomics. 15, 2597-2601. https://doi.org/10.1002/pmic.201400515
  19. Pathan, M., Keerthikumar, S., Chisanga, D., Alessandro, R., Ang, C.S., Askenase, P., Batagov, A.O., Benito-Martin, A., Camussi, G. and Clayton, A. 2017. A novel community driven software for functional enrichment analysis of extracellular vesicles data. J. Extracell. Vesicles. 6, 1321455.
  20. Pinero, J., Bravo, A., Queralt, R.N., Gutierrez, S.A, Deu, P. J., Centeno, E., Garcia, G. J., Sanz, F. and Furlong, L.I. 2017. DisGeNET: a comprehensive platform integrating information on human disease-associated genes and variants. Nucleic Acids Res. 45, D833-D839. https://doi.org/10.1093/nar/gkw943
  21. Pinero, J., Ramirez, A. J.M., Sauch-Pitarch, J., Ronzano, F., Centeno, E., Sanz, F. and Furlong, L.I. 2020. The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 48, D845-d855.
  22. Ru, J., Li, P., Wang, J., Zhou, W., Li, B., Huang, C., Li, P., Guo, Z., Tao, W. and Yang, Y., et al. 2014. TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J. Cheminform. 6, 13.
  23. Safran, M., Rosen, N., Twik, M., BarShir, R., Stein, T.I., Dahary, D., Fishilevich, S. and Lancet, D. 2021. Practical guide to life science databases, pp. 27-56, The genecards suite.
  24. Schopf, F.H., Biebl, M.M. and Buchner, J. 2017. The HSP90 chaperone machinery. Nat. Rev. Mol. Cell Biol. 18, 345-360. https://doi.org/10.1038/nrm.2017.20
  25. Shen, M., Tian, S., Li, Y., Li, Q., Xu, X., Wang, J. and Hou, T. 2012. Drug-likeness analysis of traditional Chinese medicines: 1. property distributions of drug-like compounds, non-drug-like compounds and natural compounds from traditional Chinese medicines. J. Cheminform. 4, 31.
  26. Song, B.K., Won, J.H. and Kim, S. 2016. Historical medical value of Donguibogam. J. Pharmacopuncture. 19, 16.
  27. Stebbins, C.E., Russo, A.A., Schneider, C., Rosen, N., Hartl, F.U. and Pavletich, N.P. 1997. Crystal structure of an Hsp90-geldanamycin complex: targeting of a protein chaperone by an antitumor agent. Cell 89, 239-250. https://doi.org/10.1016/S0092-8674(00)80203-2
  28. Stelzer, G., Rosen, N., Plaschkes, I., Zimmerman, S., Twik, M., Fishilevich, S., Stein, T.I., Nudel, R., Lieder, I., Mazor, Y., et al. 2016. The genecards suite: from gene data mining to disease genome sequence analyses. Curr. Protoc. Bioinformatics 54, 1.30.31-31.30.33.
  29. Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A. and Bray, F. 2021. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA. Cancer. J. Clin. 71, 209-249. https://doi.org/10.3322/caac.21660
  30. Trepel, J., Mollapour, M., Giaccone, G. and Neckers, L. 2010. Targeting the dynamic HSP90 complex in cancer. Nat. Rev. Cancer. 10, 537-549. https://doi.org/10.1038/nrc2887
  31. Trott, O. and Olson, A.J. 2010. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 31, 455-461. https://doi.org/10.1002/jcc.21334
  32. UniProt C. 2019. UniProt: a worldwide hub of protein knowledge. Nucleic. Acids. Res. 47, D506-D515. https://doi.org/10.1093/nar/gky1049
  33. Wang, H., Deng, G., Ai, M., Xu, Z., Mou, T., Yu, J., Liu, H., Wang, S. and Li, G. 2019. Hsp90ab1 stabilizes LRP5 to promote epithelial-mesenchymal transition via activating of AKT and Wnt/β-catenin signaling pathways in gastric cancer progression. Oncogene. 38, 1489-1507. https://doi.org/10.1038/s41388-018-0532-5
  34. Wang, H., Zhang, H., Zhang, X., Yin, Y., Ding, G., Tang, X., Hou, P., Sun, S. and Wang, W. 2023. Identification of coniferyl ferulate as the bioactive compound behind the xanthine oxidase inhibitory activity of Chuanxiong Rhizome. J. Funct. Foods. 100, 105378.
  35. Xu, Q., Tu, J., Dou, C., Zhang, J., Yang, L., Liu, X., Lei, K., Liu, Z., Wang, Y. and Li, L. 2017. HSP90 promotes cell glycolysis, proliferation and inhibits apoptosis by regulating PKM2 abundance via Thr-328 phosphorylation in hepatocellular carcinoma. Mol. Cancer 16, 1-16.
  36. Xu, X., Zhang, W., Huang, C., Li, Y., Yu, H., Wang, Y., Duan, J. and Ling, Y. 2012. A novel chemometric method for the prediction of human oral bioavailability. Int. J. Mol. Sci. 13, 6964-6982. https://doi.org/10.3390/ijms13066964
  37. Yin, L., Wang, Y., Lin, Y., Yu, G. and Xia, Q. 2019. Explorative analysis of the gene expression profile during liver regeneration of mouse: a microarray-based study. Artif. Cells Nanomed. Biotechnol. 47, 1113-1121. https://doi.org/10.1080/21691401.2019.1593851
  38. Zhang, H., Yin, X., Zhang, X., Zhou, M., Xu, W., Wei, Z., Song, C., Han, S. and Han, W. 2022. HSP90AB1 promotes the proliferation, migration, and glycolysis of head and neck squamous cell carcinoma. Technol. Cancer Res. Treat. 21, 15330338221118202.
  39. Zhang, M., Cui, J., Shen, F., Ye, L., Cheng, C., Li, Y., Zhang, Q., Niu, L., Hou, Y. and Bai, G. 2022. A novel mode of action for COX-2 inhibition: Targeting ATPase domain of HSP90 induces ubiquitin degradation of new client protein COX-2. Clin. Transl. Med. 12.