Utilizing Bioinformatics Technology to Explore the Potential Mechanism of Danggui Buxue Decoction Against NSCLC

doi:10.21203/rs.3.rs-970748/v1

Download PDF

Research

Utilizing Bioinformatics Technology to Explore the Potential Mechanism of Danggui Buxue Decoction Against NSCLC

https://doi.org/10.21203/rs.3.rs-970748/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 25 Feb, 2022

Read the published version in Disease Markers →

Version 1

posted

You are reading this latest preprint version

Background: While lung cancer poses a serious threat to human health, non–small cell lung cancer (NSCLC) is the most common type of lung cancer. Danggui Buxue Decoction (DBD) is a classic traditional Chinese medicine (TCM) commonly prescribed for anti-tumor use in China. However, the potential mechanism of DBD against NSCLC has not yet been expounded. Hence, this study clarified its potential mechanism through bioinformatics technology.

Methods: First, the active ingredients of DBD and the relevant targets were searched and screened by using the TCMSP database. Then the targets of NSCLC were obtained through searching GeneCards and OMIM database. Venny diagram was drawn to find the core targets, and an ingredient-disease-target network diagram was made by Cytoscape software, followed by GO functional enrichment analysis and KEGG pathway enrichment analysis via R software, and molecular docking via the Autodock Vina software , respectively. Finally, a survival analysis was performed to evaluate the prognosis of NSCLC patients through key genes.

Results: A total of 22 active ingredients of DBD were obtained. The targets found were 182 relevant to DBD and 5,773 to NSCLC, respectively, among which 140 common targets were obtained from Venny diagram. According to the enrichment degree value, five main active ingredients and ten core targets were identified. The enrichment analysis further revealed that DBD could achieve the purpose of treating NSCLC through the AGE-RAGE signaling pathway in diabetic complications. The findings from molecular docking demonstrated that quercetin and hederagenin have the best working mechanisms with PDE3A and PTGS1, and those from the survival analysis indicated that high PDE3A gene expression has a relatively poor prognosis for NSCLC patients (p < 0.05), PDE3A may be mainly distributed in the LU65 cell line originated from Asian population.

Conclusions: The present study demonstrated that DBD could treat NSCLC through the synergistic correlation between multiple ingredients, multiple targets, and multiple pathways; thus, effectively improving NSCLC prognosis.

Clinical Pharmacology

Danggui Buxue Decoction

NSCLC

Bioinformatics Technology

According to the estimates from World Health Organization's "Global Cancer Report 2020", breast cancer (2.26 million cases) has now replaced lung cancer (2.21 million cases) as the world’s most commonly diagnosed cancer [1]. The report also predicts that by 2040, due to the incidence rate being 20%, the number of new cancer cases will exceed 27 million worldwide, with one in every five people suffering from cancer. In 2020, an estimated 4.57 million new cancer cases were reported in China, of which lung cancer is by far the leading cause of the highest number of reported cases at 23.7% and is still the most frequently occurring cancer threatening the health of the Chinese people while breast cancer ranks fourth in the order of occurrence [1–3]. Among all the lung cancer subtypes, NSCLC is the most common form, accounting for 80% of all lung cancer diagnoses with a 5-year survival rate of 18% [4, 5]. As the NSCLC mortality rate is reaching towering proportions every year [6], the conventional treatment regime is still based on chemotherapy or radiotherapy, which predicts a poorer prognosis, reduces the quality of life, as well as increases the economic burden on patients and their families through expensive treatment costs [7]. Although the emergence of targeted therapy is now being considered a boon for cancer patients, the resultant increase in drug resistance and unexplained incidence of allergic reactions make it inappropriate for some patients [8, 9]. In the present scenario, several other studies involved in research on traditional Chinese medicine (TCM) have also postulated that TCM has its unique advantages while treating tumors with significant curative effects, fewer side effects, low drug resistance, and can also manifest excellent outcomes in combination with other conventional drugs [10–13].

As described in Li Gao's "Differentiation of Endogenous and Exogenous Disorders", DBD is a fixed-ratio (5:1) co-formulation consisting of Chinese medicines Danggui and Huangqi, which is also considered as a classic formula for invigorating qi and producing blood; thus, increasing hematopoietic function [14]. Duanggui is the dried root of a perennial leguminous plant Astragalus membranaceus (Fisch.) Bge.var.Mongholicus (Bge.)Hsiao(AMB) [15] whereas Huangqi is the dry root of Angelica sinensis (Oliv.) Diels (family Umbellaceae) [16]. A plethora of pharmacological studies revealed that DBD regulates immune responses [17], promotes hematopoiesis [17, 18], protects the cardiovascular system [19], shields the liver [20], as well as demonstrates anti-tumor properties [21]. DBD has a good clinical effect on NSCLC and can effectively improve the quality of life of patients with advanced NSCLC; when used alone or in combination with other conventional drugs, it can also reduce the side effects of radiotherapy and chemotherapy as well as improving their sensitivity [22–25]. Although it was observed in clinical practice that DBD is highly effective in the NSCLC treatment, the complex composition of DBD and the unclear mechanism of action makes it quite difficult to conduct further research. Henceforth, TCM being a traditional medicine can also treat diseases through multi-component, multi-target, and multi-pathway synergistic effects to achieve the purpose of curative healing [26]. However, the interconnected complexity of traditional Chinese medicine components and their indistinct mechanism of action duly limit further development and usage of such traditional medicines [27].

As a novel research approach, network pharmacology is based on high-throughput omics data analysis, computer virtual computing, and network database retrieval, which acts as an effective means for discovering active drug substances and revealing their pharmacological mechanisms of action from a crucial perspective of molecular biological networks [28]. The multi-component, multi-target, and multi-pathway characteristics of Chinese herbal compound prescriptions are the prime focus domain of network pharmacology research [29]. Although molecular docking is mainly used to study the interaction between the molecules and predict salient binding modes and relationships, it is also a mathematical simulation of drug-receptor interactions to predict the most active drug component [30].

Owing to the fact that there is no network pharmacology study in the literature to date on NSCLC treatment with DBD, this study explored the functional molecular mechanism of DBD in NSCLC treatment through network pharmacology and molecular docking as well as analyzed the survival of important targets through the Cancer Genome Atlas (TCGA) database to duly assess the target’s impact on NSCLC prognosis. In addition, this study also conducted an in-depth analysis of the expression of important targets in NSCLC tumor cell lines through the Cancer Cell Line Encyclopedia (CCLE) database. Our research plan is shown in Figure 1.

Collection of active ingredients and chemical structure of DBD

Traditional Chinese medicine systems pharmacology (TCMSP) database (https://tcmsp-e.com/) was used to collect the active ingredients of DBD. For identifying the active ingredients in DSD, “Danggui” was heavily explored in the database, followed by scheduling screening conditions through Oral Bioavailability (OB) and Drug-Likeness (DL) parameters to acquire the needed active ingredients and downloading them. Since OB represents the percentage of traditional Chinese medicine’s absorption in the human circulatory system, while DL is the similarity between the ingredients and known drugs [31], the following values ( OB > 30% and DL > 0.18) were decided as the foremost screening conditions of all active ingredients. “Huangqi” was also profoundly searched through the TCMSP database and was filtered and downloaded according to these values (OB > 30% and DL > 0.18) set as the screening conditions.

The TCMSP database contains the chemical structure of all the active ingredients. The corresponding chemical structure of obtained active ingredients can also be downloaded through the TCMSP database. Additionally, active ingredients’ Tripos molecule structure format (MOL2) files were also obtained for further utilization in the molecular docking procedure. Due to the non-availability of information about the main active ingredients, MOL2 files of all active ingredients were downloaded for further screening processes. However, in the case of active ingredients lacking MOL2 files, ChemBioDraw Ultra and ChemBio3D Ultra software were used to draw them and later saved them in MOL2 format for further usage.

The potential targets database of NSCLC

We used the GeneCards database (https://www.genecards.org/) and Online Mendelian Inheritance in Man (OMIM database) (https://omim.org/) to search the NSCLC related targets. For instance, after inserting “Non-Small Cell Lung Cancer” as a keyword in two databases, the relevant targets downloaded in databases were merged and then deduplicated to acquire the NSCLC relevant targets. After the utilization of R software to draw a Venn diagram of DBD and NSCLC targets, finally, our core targets were achieved.

Protein-Protein Interaction construction

The construction of the Protein-Protein Interaction (PPI) network was completed on the string (https://string-db.org/) website. Firstly, the "Multiple Proteins" option was selected, followed by mentioning the core targets in the List of Names and inserting "homo sapiens" in the Organism inquiry. Lastly, after inserting the default values for other parameters, the PPI network map was obtained for our perusal.

Drug-ingredient-target-disease network construction

A Cytoscape software was employed to build a drug-ingredient-target-disease visualization network diagram. After inserting the obtained drugs, ingredients, core targets, diseases, and their mutual correspondence mechanisms into the Cytoscape software, the topological characteristics of network nodes were calculated through the Cytoscape plug-in Network Analyzer that included various degrees of correlation, namely, Degree, Betweenness and Closeness, among which the Degree value was an important evaluation index [32]. Subsequently, ten main targets were also acquired through the Degree correlation factor of the Cytoscape software.

GO functional analysis and KEGG enrichment pathway analysis

In order to further understand the function of the core targets and the key DBD pathway against NSCLC, a Gene Ontology (GO) functional enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment pathway analysis were conducted. Go functional enrichment analysis was performed by using Cytoscape software and R software while the results consisted of an amalgamation of biological process, cellular component, and molecular function and were presented in a bubble chart and a histogram with a p-value where p < 0.05 indicated statistical significance. KEGG enrichment pathway analysis was also performed through Cytoscape software and R software, while their results included a bubble chart, a histogram, and some signal path diagrams with a p-value; p < 0.05 indicated statistical significance.

Molecular docking

AutoDock Vina software was used to verify the molecular docking mechanism between the top ten key targets in the PPI network and the five main active DBD ingredients. The three-dimensional structure of the ten target proteins from the RCSB PDB protein structure database (https://www.rcsb.org/) and the chemical structure of Tripos molecule structure format (MOL2) file of active ingredients from TCMSP platform was obtained. Additionally, AutoDock tools were also utilized to process the above-mentioned protein receptors and ligands, whereas its plug-in Autogrid was used to detect the docking active sites, and a molecular docking process was conducted to obtain the affinity. In our study, the binding energy ≤ –5.0 kJ·mol⁻¹ was used as the key molecule possessing a better binding affinity with the key targets [33].

Survival analysis

In order to better understand the influence of core targets on NSCLC prognosis, a precise analysis of the influence of the ten main targets on the prognosis of all squamous cell carcinomas and adenocarcinomas of NSCLC from the TCGA database (https://www.cancer.gov/about-nci/organization/ccg/research/structural-genomics/tcga) was carried out. The target gene expression was ranked from high to low and used high (50%) and low (50%) cut-off values as expression thresholds to divide the high expression group and the low expression group. The resultant worst prognosis of the high expression group relatively reflected that this specific gene might promote tumor development, and vice versa was considered as a protective factor. Data analysis was conducted on R software, and log-rank test and Kaplan-Meier method were utilized to analyze 1014 samples to obtain survival maps; p < 0.05 was considered as a statistically significant difference. Moreover, R software was also used to draw the Sankey diagram in which the width of the extended branch in the figure corresponded to the data flow size, which could be further used to depict the relationship between a certain gene's clinical characteristics, such as smoking, radiation, and racial predilection in NSCLC tumor samples and patient survival, as well as the distribution trend of high and low expressions.

Expression of core target genes in NSCLC tumor cell lines

Since NSCLC has many cell lines, the core genes are distributed differently in various cell lines. For understanding the expression of core genes in the same tumor, we used the Cancer Cell Line Encyclopedia (CCLE) database (https://sites.broadinstitute. org/ccle) to obtain the cell line gene expression matrix of NSCLC tumors and analyzed the expression of core target genes in different NSCLC tumor cell lines by R software, to select suitable cell lines for further verification.

The active ingredients and chemical structure of DBD

According to the screening conditions of active ingredients, after obtaining 22 active DBD ingredients from the TCMSP database 182 related key targets were acquired. Since the chemical structures of all active ingredients were collected through the TCMSP database, their MOL2 files were also downloaded by the TCMSP database, which is shown as detailed information in Table 1.

The NSCLC related targets

The NSCLC related targets were obtained through GeneCards and OMIM database. In the GeneCards database, 5404 related targets of NSCLC were received, while 468 related targets of NSCLC were gained from the OMIM database. After the merger of two data sets, deduplication was performed, and 5773 related targets of NSCLC were obtained. Additionally, the intersection of the DBD active ingredient targets Venn diagram and the obtained targets of NSCLC through R software yielded 140 core targets. The results are displayed in Figure 2.

Protein-Protein Interaction construction

Several interaction data sets were obtained by incorporating 140 core targets into the string platform and selecting “homo sapiens” to generate the PPI network map, and subsequently, associated protein interaction relationships were acquired. The result is shown in Figure 3. In the PPI network map, there were inherent 140 nodes and 2351 edges in which the nodes represented proteins, edges represented their innate relationships, while the colors from yellow to red represented small to large values.

Drug-ingredient-target-disease network construction

The drug-ingredient-target-disease network diagram was constructed by Cytoscape software, and the results are displayed in Figure 4. The active ingredients and maximum key targets had numerous interactions through computational inference; thus, indicating that DBD can treat NSCLC through a precise analysis of multi-component, multi-target, and multi-pathway synergistic interactions, which is quite evident from Figure 5. Using Network Analyzer of Cytoscape software, a degree analysis was conducted, as well as five main DBD ingredients and ten key targets were screened out. The five main active ingredients were quercetin, kaempferol, formononetin, isorhamnetin, and hederagenin, while the ten core targets were HSP90AA1, NCOA2, PPARG, PRKACA, NOS2, PDE3A, PTGS1, PTGS2, ADRB2, and ESR1. Simultaneously, the Cytoscape software was again utilized to investigate the potential collaborative connections between the ten core targets and thereafter demonstrated the concept that the redder the color, the closer was its connection with other key targets. The result is shown in Figure 5.

GO and KEGG enrichment analysis

Cytoscape and R software were utilized to perform GO functional enrichment analysis on DBD treatment of NSCLC to explore the possible mechanisms of 140 candidate targets for the treatment of NSCLC, and a total of 131 related biological processes, molecular functions, and cellular components were acquired with a significance value of p< 0.05. The bubble chart and histogram of the top 20 items are displayed in Figure 6. It was revealed that the co-action targets were mainly enriched in DNA-binding transcription activator activity, RNA polymerase II-specific activity, ubiquitin-like protein ligase binding, endopeptidase activity, protein serine/threonine kinase activity, cytokine receptor binding, ubiquitin-protein ligase binding, etc.

Meanwhile, a similar method was used to perform KEGG enrichment analysis on the DBD treatment regime for NSCLC to explore the possible mechanisms of 140 candidate targets for NSCLC treatment and obtained a total of 174 items with p < 0.05. The bubble chart and histogram of the top 20 items are shown in Figure 7. Our study results revealed that the co-action targets were mainly enriched in several diseases: Lipid induced atherosclerosis, Hepatitis B, shear stress-induced atherosclerosis, etc., and also participated in different signaling pathways, included AGE-RAGE signaling pathway in diabetic complications, IL-17 signaling pathway, TNF signaling pathway, etc. Figure 8 shows the relevant key targets in both the DBD and AGE-RAGE signaling pathways in diabetic complications.

Molecular Docking

The top five target proteins, as well as the top ten active DBD ingredients in the drug-ingredient-target-disease network, were verified by molecular docking to explore the intertwined interactions between the receptors and ligands. Moreover, if the binding energy was less than zero, it indicated that the ligand molecules could bind to the receptor protein spontaneously while the affinity less than–5.0 kJ·mol-1revealed that the binding affinity was good. Therefore, it can be easily discerned that the smaller the binding affinity, the greater the molecular docking process. Our study results exhibited that all active ingredients bound well to the target proteins, and every affinity was less than –5.0 kJ·mol-1 as depicted in Table 2. In particular, quercetin and hederagenin had a better binding affinity to other target proteins. It was also expressed that quercetin possessed an excellent binding affinity to PTGS1 (Affinity = –9.7 kJ·mol-1) due to the formation of hydrogen bond interactions with the PTGS1 active site, which was the main force promoting its binding capability. Our results also disclosed that these five main active ingredients had better binding affinities towards these protein targets for proper intermolecular interactions.

Survival analysis

All squamous cell carcinomas and adenocarcinomas of NSCLC cases were categorized into high and low expression groups based on the expression levels of ten core genes, and the correlation between the individual ten core gene’s expressions and the prognosis of NSClC patients was investigated by exploring the TCGA database. Our study results revealed that the highly expressed genes were associated with poor prognosis, PDE3A was associated with overall survival in all NSCLC patients (p = 0.0384), while the overall survival analysis of other nine core genes with high and low expression did not have any statistical significance (p < 0.05) which is shown in Figure 9. As (p = 0.0384) was utilized to draw the Sankey diagram (Figure 10), the resultant finding unveiled the interwoven relationship between smoking, radiation, racial predilection, the role of PDE3A, and patient survival while showing the potential connection between the width of the extended branch and patient survival.

Expression of core target genes in NSCLC tumor cell lines

The survival analysis detected that PDE3A had a close connection with patient survival. PDE3A expression in different NSCLC tumor cell lines was analyzed by R software by exploring the CCLE database. Our study results depicted that PDE3A expression was the highest in LU65 and NCI-H810, albeit noting that the abscissa represented the PDE3A expression level while the ordinate represented different NSCLC tumor cell lines, as shown in Figure 11.

While lung cancer is still the leading cause of cancer-related deaths globally, Non–small cell lung cancer (NSCLC) is the most common type of lung cancer, accounting for a large proportion of all lung cancer cases [34, 35]. Therefore, choosing the appropriate medicine integrated with standard oncologic care is particularly important for NSCLC patients. Traditional Chinese medicine is an ancient popular medicine that originated and developed in China. Since it has several unique attributes like the presence of multiple active substances, multiple key targets as well as low toxicity [36], it also possesses superior anti-tumor properties that can be applied in treating different tumors [37]. Our study determined the key active ingredients and possible detailed molecular mechanisms of DBD in the NSCLC treatment through bioinformatics technology to further improve the desired treatment effects as well as to prolong the survival rate of NSCLC patients. In this study, the active ingredients and related targets of DBD were screened from the TCMSP database and further analyzed by GeneCards and OMIM databases to screen out the potential targets for NSCLC. Additionally, the intersection of both DBD ingredients-related targets and the NSCLC related targets was employed to obtain the core targets. Subsequently, a PPI network and a drug-ingredient-target-disease network were constructed through core targets while using R software later to perform GO function and KEGG pathway enrichment analysis. The top five active DBD ingredients involved in degree analysis were quercetin, kaempferol, formononetin, isorhamnetin, and hederagenin, while the top ten main core targets included HSP90AA1, NCOA2, PPARG, PRKACA, NOS2, PDE3A, PTGS1, PTGS2, ADRB2, and ESR1. Subsequently, the analysis of five main active ingredients with ten core targets by molecular docking process discovered that the combination of quercetin and PTGS1 had the best binding affinity. Additionally, the analysis of ten main target genes survival and the NSCLC prognosis through the TCGA database proved that due to the major distribution of PDE3A in the LU65 and NCI-H810 NSCLC tumor cell lines, it has a statistical significance for the prognostic survival time of NSCLC patients,

According to the drug-ingredient-target-disease network, it was observed that DBD has five main active ingredients: quercetin, kaempferol, formononetin, isorhamnetin, and hederagenin. Since quercetin, a bioflavonoid, can induce Hsp70 inhibition involved in growth inhibition of lung cancer cells, it has a great potential as a chemosensitizer in lung cancer treatment as well as the incorporation of dietary quercetin can also be a promising option for cancer prevention [38, 39]. Kuo et al., in their study, verified that kaempferol could be used as a radiosensitizer for NSCLC in vitro and in vivo while significantly improving the lethality of tumor cells [40]. Formononetin is a novel herbal isoflavonoid, which when isolated from herbal medicine might act as a potential chemopreventive drug for lung cancer therapy through induction of cell cycle arrest and apoptosis in NSCLC cells [41]. Isorhamnetin, a traditional Chinese medicine used to treat angina pectoris and acute myocardial infarction, displays a series of anti-tumor activities [42]. Past literary insights on A549 lung cancer cells discovered that isorhamnetin at a concentration of 20 µg/mL can induce apoptosis in A549 cells, upregulates the expression of apoptotic genes Bax, Caspase-3, and p53, as well as down-regulates the expression of Bcl-2, cyclin D1, and PCNA proteins. Isorhamnetin’s mechanism of action may involve apoptosis initiation induced by down-regulation of oncogenes and up-regulation of apoptotic genes, thus proving that it can significantly inhibit the growth of A549 cells by inducing cell apoptosis [42]. Hederagenin, an oleanolic acid derivative isolated from ivy leaves by displaying potential anti-tumor activity, might become a promising therapeutic candidate for human colon cancer [43]. A study by Wang et al. reckoned that hederagenin can also induce ROS accumulation and enhance cisplatin and paclitaxel cytotoxicity in lung cancer cells by blocking autophagic flux [44]. Therefore, hederagenin also has a potential synergistic effect in the treatment of lung cancer.

After exploring 5,773 candidate NSCLC targets from GeneCards and OMIM databases, 140 common targets were obtained between the NSCLC and DBD, which were subsequently considered as potential targets for the NSCLC treatment. Through the visual analysis of the drug-ingredient-target-disease network and degree ranking, ten core genes were identified. The degree of association between these genes is shown in Figure 6. These ten core genes included HSP90AA1, NCOA2, PPARG, PRKACA, NOS2, PDE3A, PTGS1, PTGS2, ADRB2, and ESR1and may play an important role in tumor cells, particularly in the process of proliferation, migration, and apoptosis. A previous study disclosed that upregulated KCNQ1OT1 levels in NSCLC tissues and cell lines affirmed that higher KCNQ1OT1 levels were related to the poor progression-free survival of NSCLC patients [45]. Additionally, it was found that the HSP90AA1 expression was reduced after downregulating KCNQ1OT1 levels, which proved that KCNQ1OT1 positively regulated the expression of HSP90AA1 by the formation of the miR-27b-3p sponge. These data revealed the role of KCNQ1OT1 as an oncogene through the modulation of the miR-27b-3p/HSP90AA1 axis during the progression of NSCLC. Although it was suggested that HSP90AA1 could be a potential target for NSCLC treatment, another protein-coding gene, NCOA2, might be relevant for gastric cancer, liver cancer, or prostate cancer [46–48]. Although many previous studies have hypothesized that NCOA2 might become a potential therapeutic target for NSCLC, further verification is still needed to validate it [49]. PPARG expression may also act as a potential therapeutic agent for NSCLC, especially for lung squamous cell carcinoma (LSCC). Since the activation of PPARG expression can inhibit LSCC development and progress by regulating the upstream regulator and downstream marker genes, which are involved in tumor cell proliferation and protein polyubiquitination/ubiquitination [50]. Another study suggested that the loss of NOS2 reduces the growth of lung tumors and the inflammation caused by oncogenic KRAS while stating that KRAS and NOS2 jointly promote the occurrence and inflammation of lung tumors [51]. Moreover, inhibition of NOS2 may have therapeutic value for lung cancer with oncogenic KRAS mutations. It was also evident that highly methylated DNA, down-regulated PDE3A in chemoresistant NSCLC cells by forcing PDE3A expression to make A549/Cis cells sensitive to cisplatin. This result indicated that high PDE3A expression might promote the NSCLC treatment by increased efficacy of combination therapies [52]. A network pharmacology study by Wang et al. on NSCLC did not find any effect of PDE3A on NSCLC [53]. Moreover, in our study, the survival analysis of 10 main targets suggested that only PDE3A had statistical significance for the prognosis of NSCLC patients. The survival rate of PDE3A in the low expression group was higher than that in the high expression group, further confirming that reducing the expression of PDE3A may improve the quality of life of NSCLC patients. Meanwhile, through the Sankey diagram, it was evident that PDE3A was mainly distributed in the LU65 tumor cell line of NSCLC. Interestingly, the LU65 cell line originated in the Asian populations, especially East Asia, accounting for over 79.97% of inhabitants [54]. Therefore, PDE3A might be a more suitable potential therapeutic target for Asian NSCLC patients. Wang et al. analyzed 12 human plasma samples by using RNA-Seq and bioinformatics techniques and found seven key targets related to lung tumorigenesis: COX1, COX2, COX3, ND1, ND2, ND4L, and ATP6 [55]. A previous study demonstrated that the SCC pathogenesis caused by COPD is regulated by HSP90AA1, ADRB2, TBL1XR1, and HSPB1. Therefore, these genes can be used as potential therapeutic targets for the treatment of COPD-related SCC patients [56]. A 2008 study observed that the use of real-time PCR to assess the methylation of the ESR1 promoter in the blood proved very useful for the diagnosis of lung cancer, as these methylated genes might become crucial biomarkers for the early detection of lung cancer. The results also indicated that a comparative evaluation of methylation ratios before and after surgery might be a powerful tool for predicting the prognosis of lung cancer patients [57]. Simultaneously, it was suggested that ESR1 mRNA overexpression is innately associated with NSCLC prognosis [58]. Therefore, ESR1 could also become an important key target for treating NSCLC, which is consistent with the results of Wang et al. [53]. It was also estimated that as PRKACA was a tumor target, it might be helpful as therapy, but due to lack of evidence, its involvement in the development of lung cancer is debatable.

GO annotation is an important means to examine the function of gene products [59]. Through GO functional enrichment analysis, this study inferred all the important biological processes, their cellular compositions, and intermolecular functions involved in the core targets. The details are displayed in Figure 7. In KEGG enrichment analysis, the results were mainly related to the AGE-RAGE signaling pathway in diabetic complications. Accumulation of AGEs and upregulated expression of RAGE is associated with various pathological conditions, including diabetes, cardiovascular diseases, neurodegenerative disorders, and cancer. The role of AGE-RAGE signaling has been previously demonstrated in the progression of various types of cancer and other pathological disorders [60]. Therefore, the regulation of the AGE-RAGE signaling pathway is closely related to tumor development. In the molecular docking process, analysis of five main active ingredients with ten key core targets while calculating the affinity revealed that the affinity of quercetin and PTGS1 was the greatest (Affinity = –9.7 kJ·mol⁻¹), while the binding energy of hederagenin and PDE3A was also good (Affinity = –8.9 kJ·mol⁻¹). These results observed that, although PTGS1 and PDE3A might be the key targets in DBD for treating NSCLC, quercetin, and hederagenin may also be effective ingredients for the potential treatment of NSCLC.

However, it should be noted that our study was based on public databases, which had finite information and need to be continuously improved. Additionally, our study also overlooked the influence of concentration, picking time, processing method, and medication time of DBD on NSCLC. But overall, our study is worthy of further exploration and verification to extract the intricate molecular mechanisms governing therapeutic targets for NSCLC treatment.

Through several bioinformatics methods, our study results predicted the key active ingredients and related targets of DBD as well as explored the potential targets of NSCLC to devise effective drug therapies for the same. Through the enrichment analysis, it was revealed that the potential working mechanism of DBD in NSCLC treatment mainly focused on three crucial interactions: DNA-binding transcription activator activity, RNA polymerase II-specific activity, and chemical carcinogenesis-receptor activation activity. Additionally, few important observations were also discovered. Namely, the AGE-RAGE signaling pathway in diabetic complications, IL-17 signaling pathway, and TNF signaling pathway are the three main DBD signaling pathways, PDE3A and PTGS1 are the potential key genes, while quercetin and hederagenin are the two main active ingredients crucial for NSCLC treatment. Our study results also depicted that DBD can achieve the goal of treating NSCLC while improving the prognosis of NSCLC patients through a synergistic interplay between multiple ingredients, multiple targets, and multiple signaling pathways. Nevertheless, further experimental validation is also required to ascertain effective therapeutic approaches for NSCLC patients. At a later stage, our team might conduct a following pharmacological study on important signaling pathways and targets of DBD against NSCLC, and try to provide a functional basis and rational references toward a new future for drug development and targeted therapy.

CCLE: Cancer Cell Line Encyclopedia; DBD: Danggui Buxue Decoction; DL: Druglikeness; GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; NSCLC: Non-Small Cell Lung Cancer; OB: Oral Bioavailability; OMIM: Online Mendelian Inheritance in Man; PPI: Protein-Protein Interaction; TCM: Traditional Chinese Medicine; TCMSP: Traditional Chinese Medicine Systems Pharmacology; TCGA: The Cancer Genome Atlas

Acknowledgements

Not applicable.

Authors’ contributions

BY and GYL designed the study; HX supervised the study; BY, GYL and YLL performed the research, analyzed the data and wrote the manuscript; HX revised the manuscript; ZYL, GYL and YLL prepared figures and tables. All authors reviewed the manuscript. And all authors read and approved the final version of the manuscript.

Funding

This study was supported by Hypertension Disease (Taige) Special Scientific Research Project of Sichuan Medical Association (2019TG28)

Availability of data and materials

All data can be fully available without any restrictions.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

All authors declare that they have no conflicts of interest.

Latest global cancer data: Cancer burden rises to 19.3 million new cases and 10.0 million cancer deaths in 2020. Dec 15, 2020. from https://www.iarc.who.int/news-events/latest-global-cancer-data-cancer-burden-rises-to-19-3-million-new-cases-and-10-0-million-cancer-deaths-in-2020/
Cao W, Chen HD, Yu YW, Li N, Chen WQ. Changing profiles of cancer burden worldwide and in China: a secondary analysis of the global cancer statistics 2020. Chin Med J (Engl). 2021; 134(7): 783–791.
Zeng H, Chen W, Zheng R, Zhang S, Ji JS, Zou X, Xia C, Sun K, Yang Z, Li H, Wang N, Han R, Liu S, Li H, Mu H, He Y, Xu Y, Fu Z, Zhou Y, Jiang J, Yang Y, Chen J, Wei K, Fan D, Wang J, Fu F, Zhao D, Song G, Chen J, Jiang C, Zhou X, Gu X, Jin F, Li Q, Li Y, Wu T, Yan C, Dong J, Hua Z, Baade P, Bray F, Jemal A, Yu XQ, He J. Changing cancer survival in China during 2003–15: a pooled analysis of 17 population-based cancer registries. Lancet Glob Health. 2018; 6(5): e555-e567.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020; 70(1): 7–30.
Calabrese L, Mariette X. The evolving role of the rheumatologist in the management of immune-related adverse events (irAEs) caused by cancer immunotherapy. Ann Rheum Dis. 2018; 77(2): 162–164.
Ye Z, Huang Y, Ke J, Zhu X, Leng S, Luo H. Breakthrough in targeted therapy for non-small cell lung cancer. Biomed Pharmacother. 2021; 133: 111079.
Zhou C, Wang J, Bu H, Wang B, Han B, Lu Y, Wang Z, Zhu B, Wang Z, Song Q, Ren S, Lin D, He Y, Hu X, Zhao H, Qin S. Chinese Experts Consensus on Immune Checkpoint Inhibitors for Non-small Cell Lung Cancer (2019 version). Zhongguo Fei Ai Za Zhi. 2020; 23(2): 65–76.
Lee YT, Tan YJ, Oon CE. Molecular targeted therapy: Treating cancer with specificity. Eur J Pharmacol. 2018; 834: 188–196.
Ransohoff JD, Kwong BY. Cutaneous Adverse Events of Targeted Therapies for Hematolymphoid Malignancies. Clin Lymphoma Myeloma Leuk. 2017; (12): 834–851.
Liu Y, Yang S, Wang K, Lu J, Bao X, Wang R, Qiu Y, Wang T, Yu H. Cellular senescence and cancer: Focusing on traditional Chinese medicine and natural products. Cell Prolif. 2020; 53(10): e12894.
Yan Z, Lai Z, Lin J. Anticancer Properties of Traditional Chinese Medicine. Comb Chem High Throughput Screen. 2017; 20(5): 423–429.
Wang S, Long S, Deng Z, Wu W. Positive Role of Chinese Herbal Medicine in Cancer Immune Regulation. Am J Chin Med. 2020; 48(7): 1577–1592.
Su XL, Wang JW, Che H, Wang CF, Jiang H, Lei X, Zhao W, Kuang HX, Wang QH. Clinical application and mechanism of traditional Chinese medicine in treatment of lung cancer. Chin Med J (Engl). 2020; 133(24): 2987–2997.
Wang P, Liang YZ. Chemical composition and inhibitory effect on hepatic fibrosis of Danggui Buxue Decoction. Fitoterapia. 2010; 81(7): 793–8.
Lan S, Duan J, Zeng N, Yu B, Yang X, Ning H, Huang Y, Rao Y. Network pharmacology-based screening of the active ingredients and mechanisms of Huangqi against aging. Medicine (Baltimore). 2021; 100(17): e25660.
Fu J, Wang Z, Huang L, Zheng S, Wang D, Chen S, Zhang H, Yang S. Review of the botanical characteristics, phytochemistry, and pharmacology of Astragalus membranaceus (Huangqi). Phytother Res. 2014; 28(9): 1275-83.
Shi XQ, Yue SJ, Tang YP, Chen YY, Zhou GS, Zhang J, Zhu ZH, Liu P, Duan JA. A network pharmacology approach to investigate the blood enriching mechanism of Danggui buxue Decoction. J Ethnopharmacol. 2019; 235: 227–242.
Li F, Tang R, Chen LB, Zhang KS, Huang XP, Deng CQ. Effects of Astragalus Combined with Angelica on Bone Marrow Hematopoiesis Suppression Induced by Cyclophosphamide in Mice. Biol Pharm Bull. 2017; 40(5): 598–609.
Li YD, Ma YH, Zhao JX, Zhao XK. Protection of ultra-filtration extract from Danggui Buxue Decoction on oxidative damage in cardiomyocytes of neonatal rats and its mechanism. Chin J Integr Med. 2011; 17(11): 854–9.
Chen Y, Tao YY, Li FH. Effects of Danggui Buxue Decoction on liver fibrosis and hepatic lipid peroxidation in rats. Zhongguo Zhong Xi Yi Jie He Za Zhi. 2008; 28(1): 39-42.
Feng SH, Zhao B, Zhan X, Motanyane R, Wang SM, Li A. Danggui Buxue Decoction in the Treatment of Metastatic Colon Cancer: Network Pharmacology Analysis and Experimental Validation. Drug Des Devel Ther. 2021; 15: 705–720.
Du QC, Yang KZ, Sun XF. Efficacy of auxiliary therapy with Danggui Buxue Decoction No.1 in treating patients of non-small cell lung cancer at peri-operational stage. Chin J Integr Med. 2009; 15(3): 184–8.
Sun X, Xu X, Chen Y, Guan R, Cheng T, Wang Y, Jin R, Song M, Hang T. Danggui Buxue Decoction Sensitizes the Response of Non-Small-Cell Lung Cancer to Gemcitabine via Regulating Deoxycytidine Kinase and P-glycoprotein. Molecules. 2019; 24(10): 2011.
Liu Y, Chang M, Hu Z, Xu X, Wu W, Ning M, Hang T, Song M. Danggui Buxue Decoction enhances the anticancer activity of gemcitabine and alleviates gemcitabine-induced myelosuppression. J Ethnopharmacol. 2021; 273: 113965.
Shi HY, Zhang LY, Tian YZ, Shao LL. Effect of Jiawei Danggui Buxue Decoction on immune function and quality of life in patients with advanced non-small cell lung cancer chemotherapy. Journal of New Chinese Medicine. 2017, 49(03): 106–108.
Ma YM, Zhang XZ, Su ZZ, Li N, Cao L, Ding G, Wang ZZ, Xiao W. Insight into the molecular mechanism of a herbal injection by integrating network pharmacology and in vitro. J Ethnopharmacol. 2015; 173: 91–9
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(1): 161–182.
Luo TT, Lu Y, Yan SK, Xiao X, Rong XL, Guo J. Network Pharmacology in Research of Chinese Medicine Formula: Methodology, Application and Prospective. Chin J Integr Med. 2020; 26(1): 72–80.
Li S, Zhang B, Jiang D, Wei Y, Zhang N. Herb network construction and co-module analysis for uncovering the combination rule of traditional Chinese herbal formulae. BMC Bioinformatics. 2010; 11 Suppl 11(Suppl 11): S6.
Kaur T, Madgulkar A, Bhalekar M, Asgaonkar K. Molecular Docking in Formulation and Development. Curr Drug Discov Technol. 2019; 16(1): 30–39.
Wan Y, Xu L, Liu Z, Yang M, Jiang X, Zhang Q, Huang J. Utilising network pharmacology to explore the underlying mechanism of Wumei Pill in treating pancreatic neoplasms. BMC Complement Altern Med. 2019; 19(1): 158.
Ge Q, Chen L, Tang M, Zhang S, Liu L, Gao L, Ma S, Kong M, Yao Q, Feng F, Chen K. Analysis of mulberry leaf components in the treatment of diabetes using network pharmacology. Eur J Pharmacol. 2018; 833: 50–62.
Ruan X, Du P, Zhao K, Huang J, Xia H, Dai D, Huang S, Cui X, Liu L, Zhang J. Mechanism of Dayuanyin in the treatment of coronavirus disease 2019 based on network pharmacology and molecular docking. Chin Med. 2020, 15: 62.
Wang M, Herbst RS, Boshoff C. Toward personalized treatment approaches for non-small-cell lung cancer. Nat Med. 2021; 27(8): 1345–1356.
Judd J, Borghaei H. Combining Immunotherapy and Chemotherapy for Non-Small Cell Lung Cancer. Thorac Surg Clin. 2020; 30(2): 199–206.
Wang K, Chen Q, Shao Y, Yin S, Liu C, Liu Y, Wang R, Wang T, Qiu Y, Yu H. Anticancer activities of TCM and their active components against tumor metastasis. Biomed Pharmacother. 2021; 133: 111044.
Fan DC, Qi JY, Zhang MZ. Insights of Chinese medicine on ventricular remodeling: Multiple-targets, individualized-treatment. Chin J Integr Med. 2017; 23(9): 643-647.
Lee SH, Lee EJ, Min KH, Hur GY, Lee SH, Lee SY, Kim JH, Shin C, Shim JJ, In KH, Kang KH, Lee SY. Quercetin Enhances Chemosensitivity to Gemcitabine in Lung Cancer Cells by Inhibiting Heat Shock Protein 70 Expression. Clin Lung Cancer. 2015; 16(6): e235–43.
Murakami A, Ashida H, Terao J. Multitargeted cancer prevention by quercetin. Cancer Lett. 2008; 269(2): 315–25.
Kuo WT, Tsai YC, Wu HC, Ho YJ, Chen YS, Yao CH, Yao CH. Radiosensitization of non-small cell lung cancer by kaempferol. Oncol Rep. 2015; 34(5): 2351–6.
Yang Y, Zhao Y, Ai X, Cheng B, Lu S. Formononetin suppresses the proliferation of human non-small cell lung cancer through induction of cell cycle arrest and apoptosis. Int J Clin Exp Pathol. 2014; 7(12): 8453–61.
Li Q, Ren FQ, Yang CL, Zhou LM, Liu YY, Xiao J, Zhu L, Wang ZG. Anti-proliferation effects of isorhamnetin on lung cancer cells in vitro and in vivo. Asian Pac J Cancer Prev. 2015; 16(7): 3035–42.
Liu BX, Zhou JY, Li Y, Zou X, Wu J, Gu JF, Yuan JR, Zhao BJ, Feng L, Jia XB, Wang RP. Hederagenin from the leaves of ivy (Hedera helix L.) induces apoptosis in human LoVo colon cells through the mitochondrial pathway. BMC Complement Altern Med. 2014; 14: 412.
Wang K, Liu X, Liu Q, Ho IH, Wei X, Yin T, Zhan Y, Zhang W, Zhang W, Chen B, Gu J, Tan Y, Zhang L, Chan MT, Wu WK, Du B, Xiao J. Hederagenin potentiated cisplatin-and paclitaxel-mediated cytotoxicity by impairing autophagy in lung cancer cells. Cell Death Dis. 2020; 11(8): 611.
Dong Z, Yang P, Qiu X, Liang S, Guan B, Yang H, Li F, Sun L, Liu H, Zou G, Zhao K. KCNQ1OT1 facilitates progression of non-small-cell lung carcinoma via modulating miRNA-27b-3p/HSP90AA1 axis. J Cell Physiol. 2019; 234(7): 11304–11314.
Lin Z, Yang F, Lu D, Sun W, Zhu G, Lan B. Knockdown of NCOA2 Inhibits the Growth and Progression of Gastric Cancer by Affecting the Wnt Signaling Pathway-Related Protein Expression. Technol Cancer Res Treat. 2020; 19: 1533033820928072.
O'Donnell KA, Keng VW, York B, Reineke EL, Seo D, Fan D, Silverstein KA, Schrum CT, Xie WR, Mularoni L, Wheelan SJ, Torbenson MS, O'Malley BW, Largaespada DA, Boeke JD. A Sleeping Beauty mutagenesis screen reveals a tumor suppressor role for Ncoa2/Src-2 in liver cancer. Proc Natl Acad Sci U S A. 2012; 109(21): E1377–86.
Silva MP, Barros-Silva JD, Vieira J, Lisboa S, Torres L, Correia C, Vieira-Coimbra M, Martins AT, Jerónimo C, Henrique R, Paulo P, Teixeira MR. NCOA2 is a candidate target gene of 8q gain associated with clinically aggressive prostate cancer. Genes Chromosomes Cancer. 2016; 55(4): 365–74.
Liu J, He Z, Li S, Huang W, Ren Z. Network Pharmacology-Based Analysis of the Effects of Corydalis decumbens (Thunb.) Pers. in Non-Small Cell Lung Cancer. Evid Based Complement Alternat Med. 2021; 2021: 4341517.
Shi S, Yu G, Huang B, Mi Y, Kang Y, Simon JP. PPARG Could Work as a Valid Therapeutic Strategy for the Treatment of Lung Squamous Cell Carcinoma. PPAR Res. 2020; 2020: 2510951.
Okayama H, Saito M, Oue N, Weiss JM, Stauffer J, Takenoshita S, Wiltrout RH, Hussain SP, Harris CC. NOS2 enhances KRAS-induced lung carcinogenesis, inflammation and microRNA-21 expression. Int J Cancer. 2013; 132(1): 9–18.
Tian FM, Zhong CY, Wang XN, Meng Y. PDE3A is hypermethylated in cisplatin resistant non-small cell lung cancer cells and is a modulator of chemotherapy response. Eur Rev Med Pharmacol Sci. 2017 Jun;21(11):2635–2641.
Wang Y, Zhang Y, Wang Y, Shu X, Lu C, Shao S, Liu X, Yang C, Luo J, Du Q. Using Network Pharmacology and Molecular Docking to Explore the Mechanism of Shan Ci Gu (Cremastra appendiculata) Against Non-Small Cell Lung Cancer. Front Chem. 2021; 9: 682862.
Dutil J, Chen Z, Monteiro AN, Teer JK, Eschrich SA. An Interactive Resource to Probe Genetic Diversity and Estimated Ancestry in Cancer Cell Lines. Cancer Res. 2019; 79(7): 1263–1273.
Wang L, Wang J, Jia E, Liu Z, Ge Q, Zhao X. Plasma RNA sequencing of extracellular RNAs reveals potential biomarkers for non-small cell lung cancer. Clin Biochem. 2020; 83: 65–73.
Wang L, Zhao H, Zhang L, Luo H, Chen Q, Zuo X. HSP90AA1, ADRB2, TBL1XR1 and HSPB1 are chronic obstructive pulmonary disease-related genes that facilitate squamous cell lung cancer progression. Oncol Lett. 2020; 19(3): 2115–2122.
Suga Y, Miyajima K, Oikawa T, Maeda J, Usuda J, Kajiwara N, Ohira T, Uchida O, Tsuboi M, Hirano T, Kato H, Ikeda N. Quantitative p16 and ESR1 methylation in the peripheral blood of patients with non-small cell lung cancer. Oncol Rep. 2008; 20(5): 1137–42.
Teng Y, Ding Y, Zhang M, Chen X, Wang X, Yu H, Liu C, Lv H, Zhang R. Genome-wide haplotype association study identifies risk genes for non-small cell lung cancer. J Theor Biol. 2018; 456: 84–90.
Leale G, Baya AE, Milone DH, Granitto PM, Stegmayer G. Inferring Unknown Biological Function by Integration of GO Annotations and Gene Expression Data. IEEE/ACM Trans Comput Biol Bioinform. 2018; 15(1): 168–180.
Waghela BN, Vaidya FU, Ranjan K, Chhipa AS, Tiwari BS, Pathak C. AGE-RAGE synergy influences programmed cell death signaling to promote cancer. Mol Cell Biochem. 2021; 476(2): 585–598.

Due to technical limitations, table 1 and 2 docx are only available as a download in the Supplemental Files section.

Table1.docx
Table 1 Information of active ingredients and related targets from DBD
Table2.docx
Table 2 Molecular docking results of 5 main active ingredients and 10 core targets

Download PDF

Journal Publication

published 25 Feb, 2022

Read the published version in Disease Markers →

Version 1

posted

You are reading this latest preprint version

Utilizing Bioinformatics Technology to Explore the Potential Mechanism of Danggui Buxue Decoction Against NSCLC

Status:

Journal Publication

Version 1

Abstract

Figures

Background

Methods

Results

Discussion

Conclusion

Abbreviations

Declarations

References

Tables

Supplementary Files

Status:

Journal Publication

Version 1