Multitarget and Multipathway Regulation of Zhenqi Fuzheng Granule against Non-Small Cell Lung Cancer Based On Network Pharmacology and Molecular Docking

Zhou, Yanqing; Wu, Chenxi; Qian, Xian; Zhou, Jin; Li, Chengjian; Jiao, Yang; Li, Yue Yue; Zhao, Liang

doi:https://doi.org/10.1155/2022/5967078

Evidence-Based Complementary and Alternative Medicine

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Abbreviations Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 5967078 | https://doi.org/10.1155/2022/5967078

Multitarget and Multipathway Regulation of Zhenqi Fuzheng Granule against Non-Small Cell Lung Cancer Based On Network Pharmacology and Molecular Docking

Yanqing Zhou,¹Chenxi Wu,²Xian Qian,¹Jin Zhou,¹Chengjian Li,¹Yang Jiao,²Yue Yue Li,³and Liang Zhao¹

Academic Editor: Bhushan Patwardhan

Received25 Nov 2021

Revised31 Mar 2022

Accepted17 Oct 2022

Published17 Nov 2022

Abstract

Background and Objective. The morbidity and mortality rates of non-small cell lung cancer (NSCLC) remain high. Zhenqi Fuzheng (ZQFZ) granule, which consists of Astragali Radix and Ligustri Lucidi Fructus, is commonly used to improve the immunity of cancer patients. However, the mechanism of ZQFZ granule against NSCLC is still unclear. In this study, the network pharmacology and molecular docking approaches were used to investigate the potential mechanism of ZQFZ granule on NSCLC. Methods. The ingredients in the ZQFZ granule were considered in one study based on UPLC, and the potential targets were predicted in the SwissTargetPrediction database. NSCLC targets were gathered from GeneCards, OMIM, and TTD databases. The ingredient-target-NSCLC network was drawn by Cytoscape. The protein–protein interaction was obtained from the STRING database, and the gene function and biological pathways were analyzed by Metascape. AutoDock Vina was used to verify the molecular docking between the key compounds and core targets, and PyMol visualized the results. Results. 244 targets were related to 13 candidate compounds and 1904 targets were related to NSCLC, of which a total of 106 anti-NSCLC targets were predicted. The compound-target-NSCLC network indicated that sinapinic acid, ferulic acid, asiatic acid, pratensein, and glycitein might be the key components for treating NSCLC. The 41 vital targets (out of 106 targets) above the median calculated by PPI degree were selected for bioinformatics analysis. The top 10 targets out of 41 ranked by MCC were IL-6, SRC, CTNNB1, STAT3, CASP3, TNF, EGFR, MAPK8, HSP90AA1, and PTGS2. ZQFZ granule treatment for NSCLC involved many pathways through KEGG analyses, which included pathways in cancer (hsa05200), proteoglycans in cancer (hsa05205), endocrine resistance (hsa01522), microRNAs in cancer (hsa05206), PI3K-Akt signaling pathway (hsa04151), and IL-17 signaling pathway (hsa04657). Molecular docking studies revealed that sinapinic acid, ferulic acid, asiatic acid, pratensein, and glycitein had good infinity with most core targets. Conclusions. This study indicated that ZQFZ granule with multicompounds could treat NSCLC through multitargets and multipathways.

1. Introduction

GLOBOCAN 2020 showed that there were an estimated 19 million new cancer cases and 10 million deaths worldwide in 2020. The cancer morbidity and mortality rates remained high, with breast cancer ranking first among women and lung cancer being the leading cause of cancer death [1]. Surgery, radiotherapy, chemotherapy, targeted therapy, immunotherapy, and adjuvant therapy with Chinese medicine are the common means of cancer treatment. It is well known that Chinese medicine compounding has a long history with multicomponent and multitarget characteristics. In complementary treatment, Chinese medicine compounding has a great effect on increasing chemotherapy drugs' chemosensitivity, reducing adverse drug reactions and toxicity, relieving patients’ pain, and improving life quality particularly [2].

As a traditional Chinese medicine (TCM), ZQFZ granule has the ability to improve immunity and protect bone marrow and adrenal cortex functions [3]. Combined with surgery, chemotherapy, and radiotherapy, ZQFZ granule could effectively reduce adverse effects caused by chemotherapy and promote recovery of immune function, subsequently improving the patient’s quality of life and ultimately decreasing the recurrence rate. A prospective, open-label, randomized controlled trial confirmed that TCM treatment could promote NSCLC patients’ quality of life, relieve symptoms, and reduce adverse events. The TCM was composed of the immuno-boosting drug called Kangai injection, herbal decoction, and ZQFZ capsule [4].

ZQFZ granule consists of Astragali Radix (AsR, dried root of Astragalus membranaceus, known as HuangQi in Chinese) and Ligustri Lucidi Fructus (LLF, dried ripe fruit of Ligustrum lucidum, known as NvZhenZi in Chinese). Triterpenoids and phenylethanoid glycosides are the two major types of constituents present in AsR [5]. As one of the most commonly used herbal remedies for cancer in China, AsR can activate immune regulation, inhibit cell proliferation, and attenuate adverse effects caused by cytotoxic therapy [6]. Similarly, the chemical constituents of LLF mainly include saponins, flavonoids, and nitrogen-containing chemicals [7]. Studies have shown that LLF can enhance the sensitivity of adriamycin-induced apoptosis, improve immunity, and promote apoptosis, so as to inhibit tumor growth [8]. The classical formulation of ZQFZ is a (2 : 1, w/w) mixture of AsR and LLF. The mixture was decocted twice with 10–30 times water, and the filtrate was concentrated into a thick paste, mixed evenly with sucrose, and dried into granules [9].

However, due to the complexity of the components of Chinese herbal compounds, the mechanisms and targets of reaction are difficult to elucidate. The mechanisms are being studied by examining the impacts of monomeric components, which are not exhaustive. Based on the whole perspective, network pharmacology combines network biology and multipharmacology to compensate for the shortage of multicomponent system research of TCM. In this study, we focused primarily on the active components from a recent study based on ultra-performance liquid chromatography (UPLC) [10] and explored the molecular mechanism of ZQFZ granule for the treatment of NSCLC. The flowchart is shown in Figure 1.

2. Materials and Methods

2.1. Collection of the Main Ingredients from ZQFZ Granule and ADME Evaluation

The ingredients in ZQFZ granule were considered from one study based on UPLC [10]. The main active compounds were selected according to Lipinski’s rule of five (RO5) [11]. The criteria of Lipinski’s RO5 were as follows: a molecule weight (MW) should be less than 500, the number of hydrogen bond donors (Hdon) should be less than 5, the number of hydrogen bond acceptors (Hacc) should be less than 10, lipid-water partition coefficient (LogP) no more than 5, and the number of rotatable bonds (Rbon) no more than 10. The chemical structures and canonical SMILES were searched in PubChem (https://pubchem.ncbi.nlm.nih.gov/) and the ADME of compounds were uploaded into SwissADME (https://www.swissadme.ch/) to evaluate the druggability.

2.2. Screening of Targets of the Vital Ingredients from ZQFZ Granule in NSCLC

The SwissTargetPrediction (https://www.swisstargetprediction.ch/) was used to predict the potential targets of vital ingredients from the ZQFZ granule. The targets with a probability greater than zero were included. NSCLC targets were gathered from GeneCards (https://www.genecards.org/), OMIM (https://omim.org/), and TTD (https://db.idrblab.net/ttd/) by the keywords “non-small cell lung cancer” and “non-small cell lung carcinoma.” The common genes, probably the targets for ZQFZ granule to treat NSCLC, were assessed by Venny 2.1.0 (https://bioinfogp.cnb.csic.es/tools/venny/). An ingredient-disease-target network was visualized by Cytoscape 3.7.2 software.

2.3. Protein–Protein Interaction (PPI) Analysis and Vital Target Screening

The PPI of common genes was analyzed in STRING (https://string-db.org/). The organism was set to Homo sapiens and the confidence of 0.4 was selected as significant. The active interaction sources for PPI analysis and median calculation included known, predicted, and other interactions. The PPI network was analyzed and visualized using the Cytoscape software. The node size represents the degree value and the edge represents the connection between proteins. In this study, the vital genes, which were above the median calculated by degree, were then used for biological function analysis. Moreover, the top 10 proteins ranked by maximal clique centrality (MCC) using the cytoHubba plug-in were selected as core targets.

2.4. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathways Enrichment Analyses

GO analysis of the vital genes, including the biological process (BP), cellular component (CC), and molecular function (MF), was enriched and analyzed by the Metascape database (https://metascape.org/gp). The screening thresholds were a minimum of 3 genes, and an enrichment factor greater than 1.5. The same method and settings were used for KEGG pathway enrichment. Finally, the results were plotted as bubble charts using the online bioinformatics tool (https://bioinformatics.com.cn/). Furthermore, the Cytoscape plug-in ClueGo was used for analyzing the relevant KEGG pathway.

2.5. Molecular Docking

To verify the binding of the vital ingredients of ZQFZ granule with predicted core genes, compounds in mol2 format were downloaded from TCMSP (https://tcmsp-e.com/) and the core protein structures were acquired from RCSB PDB (https://www.rcsb.org/). The compounds and proteins were imported to AutoDockTools 1.5.6 software for dehydration, hydrogenation, and other pretreatments. Then, molecular docking for analyzing the binding activity was performed using AutoDock Vina 1.1.2 [12] software, and some great results were visualized using PyMol 4.6.0 software.

3. Results

3.1. The Chemical Structure and ADME Properties of the Vital Ingredients from ZQFZ Granule

A recent study has identified 95 chemical components from ZQFZ granule including 15 batches from 3 producers based on UPLC [10]. The chemical structures of the compounds were searched in the PubChem database, and ADME properties were evaluated in SwissADME. A total of 15 compounds met Lipinski’s RO5, indicating that these components had good drug-like properties. Moreover, potential target genes were predicted based on the chemical structure via the SwissTargetPrediction database, and two of them had no targets. Therefore, we selected 13 compounds (Figure 2 and Table 1) as candidate compounds.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

(m)

3.2. Acquisition of Targets of the Candidate Ingredients from ZQFZ Granule in NSCLC

Based on the SwissTargetPrediction database, a total of 244 targets for the 13 candidate compounds were acquired. In addition, results from GeneCards (https://www.genecards.org/), Online Mendelian Inheritance in Man (OMIM, https://omim.org/), and Therapeutic Target Database (TTD, https://db.idrblab.net/ttd/) for non-small cell lung cancer identified 1904 targets relevant to NSCLC. Furthermore, the Venn diagram reflected the common 106 target genes of the candidate compounds and NSCLC for further analysis (Figure 3(a)). The 106 common targets of ZQFZ granule against NSCLC are listed in Supplementary Table S1. Moreover, these 106 target proteins were grouped into seven different classes according to cellular function via PANTHER (https://www.pantherdb.org/), with protein-modifying enzyme (PC00260, 31.7%) being the most enriched class (Figure 3(b)). These protein-modifying enzymes include protease (PC00190), which covers ACE, CASP3, ECE1, F2, MME, MMP1, MMP12, MMP13, MMP2, MMP9, PLAT, PLAU, and PREP, protein phosphatase (PC00195) that covers ACP1, CDC25A, CDC25C, DUSP1, PTPN1, PTPN11, PTPN6, and PTPRC, and nonreceptor serine/threonine protein kinase (PC00167) including CDK4, ILK, MAPK14, MAPK8, and RAF1 (Figure 3(c)). The results above implied that the candidate ingredients from the ZQFZ granule exerted anti-NSCLC effects through a variety of targets and biological functions.

(a)

(b)

(c)

3.3. Construction of a Candidate Ingredients from the ZQFZ Granule Target NSCLC Network

To study the mechanism of ZQFZ granule in the treatment of NSCLC, 106 common targets and 13 candidate compounds were used for constructing the potential ZQFZ granule target NSCLC network. As shown in Figure 4, all compounds were associated with multiple targets, generating a total of 344 component-target connections between the 13 compounds and 106 targets. The network revealed sinapinic acid (degree = 40) had the most targets, followed by ferulic acid (degree = 33), asiatic acid (degree = 30), pratensein (degree = 24), and glycitein (degree = 21), indicating that these compounds from ZQFZ granule were highly likely to become key components for treating NSCLC.

3.4. Protein–Protein Interaction (PPI) Analysis

The 106 common targets of the predicted ZQFZ granule and NSCLC were explored in the STRING database and protein interaction relationships were gained (Figure 5(a)). To further study the relationship between these targets, 41 vital genes above the median calculated by degree were selected and the protein interaction relationships were analyzed via Cytoscape software. Table 2 shows the details of the 41 vital targets. The PPI network generated 41 nodes, which represent proteins, and 552 edges which represent the connections between proteins. The larger node means the greater degree (Figure 5(b)). Furthermore, the top 10 targets ranked by MCC using the cytoHubba plug-in were IL-6, SRC, CTNNB1, STAT3, CASP3, TNF, EGFR, MAPK8, HSP90AA1, and PTGS2 (Figure 5(c)) and selected as core targets for interaction with NSCLC.

(a)

(b)

(c)

Figure 5

Protein–protein interaction (PPI) network: (a) PPI network of the 106 common targets between the 13 candidate compounds from ZQFZ granule and NSCLC. The PPI was analyzed in the STRING database and visualized by Cytoscape software. Nodes represent target proteins and edges represent interactions among targets. Yellow nodes represent the 41 vital targets above the median calculated by degree. (b) PPI network of the 41 vital targets. The larger the node is, the greater the degree is. (c) The top 10 targets in the PPI network ranked by maximal clique centrality (MCC) using the cytoHubba plug-in.

3.5. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathways Enrichment Analyses

The 41 vital targets between ZQFZ granule and NSCLC were further analyzed by GO enrichment analysis and KEGG pathway analysis by Metascape. The GO results showed 421 entries for BP, 49 entries for CC, and 82 entries for MF. The top 5 BPs, CCs, and MFs are shown in bubble charts (Figure 6(a)). The top 5 biological processes were related to oxidative stress (GO: 0006979), protein phosphorylation (GO: 0001934), and transmembrane receptor protein-tyrosine kinase signaling pathway (GO: 0007169). Cellular components analysis revealed including membrane raft (GO: 0045121), receptor complex (GO: 0043235), and side of the membrane (GO: 0098552). For molecular functions, the targets were enriched in protein kinase binding (GO: 0019901), transmembrane receptor protein-tyrosine kinase activity (GO: 0004714), and protein domain-specific binding (GO: 0019904). The complex list of the top 5 BP, CC, and MF pathway terms enriched in ZQFZ granule against NSCLC is shown in Supplementary Table S2.

(a)

(b)

In addition, the KEGG pathway revealed that 206 pathways were involved, and the bubble chart (Figure 6(b)) reflects the top 20 terms. Table 3 displays detailed information about the top 20 pathways related to pathways in cancer (hsa05200), proteoglycans in cancer (hsa05205), endocrine resistance (hsa01522), microRNAs in cancer (hsa05206), PI3K-Akt signaling pathway (hsa04151), and IL-17 signaling pathway (hsa04657). In particular, 8 out of 10 core targets (IL-6, CTNNB1, STAT3, CASP3, EGFR, MAPK8, HSP90AA1, and PTGS2) involved in the first-term pathways in cancer (hsa05200), and 7 core targets (IL-6, SRC, CTNNB1, STAT3, CASP3, TNF, and EGFR) involved in the second-term proteoglycans in cancer (hsa05205). A compound-target-pathway network (Figure 7) shows the connection between the top 20 pathways, targets, and compounds. The predictive diagram of the mechanism of the proteoglycans pathway is displayed in Figure 8.

3.6. Molecular Docking Simulation

According to the potential compounds from the ZQFZ granule target NSCLC network (Figure 4), sinapinic acid, ferulic acid, asiatic acid, pratensein, and glycitein had a high number of targets against NSCLC and were the top 5 compounds ranked by degree. The PDBIDs and resolutions of the core target proteins are listed in Supplementary Table S3. Thus, molecular docking predicted the binding of the five compounds to core targets related to NSCLC. The affinity is usually regarded as the binding effect of the compound to the target; the lower the absolute value of affinity is, the more stable the binding is. As shown in Figure 9, sinapinic acid showed a strong binding to TNF, with an affinity value of −7.1 kcal/mol; ferulic acid showed a strong binding to TNF and HSP90AA1, with affinity values of −7.2 and −7.1 kcal/mol; asiatic acid showed the strong binding to IL-6, SRC, CASP3, EGFR, and HSP90AA1 with affinity values of −7.6, −8.7, −8.7, −7.5, and −7.3 kcal/mol; pratensein showed the strong binding to SRC, CASP3, TNF, EGFR, MAPK8, and HSP90AA1 with affinity values of −8.1, −7.8, −8.9, −8.2, −9.0, and −8.3 kcal/mol; and glycitein showed a strong binding to SRC, CASP3, TNF, EGFR, MAPK8, HSP90AA1, and PTGS2, with affinity values of −8.1, −7.4, −8.9, −8.2, −8.3, −8.3, and −7.1 kcal/mol. In general, an affinity value of less than −4.25 kcal/mol means some binding capacity, the affinity value of less than −5.0 kcal/mol means good binding capacity, and the affinity value of less than −7.0 kcal/mol means strong binding capacity. Table 4 shows some representative docking results with the core targets.

As shown in Table 4, the structure of asiatic acid could form one hydrogen bond, two hydrogen bonds, and two hydrogen bonds with Tyr103, Tyr32, and Asp56, respectively, in IL-6. Respectively, asiatic acid could interact with Leu273 and Asp404 via one hydrogen bond in SRC. The structure of pratensein could interact with Asp207 and Lys281 through one hydrogen bond in CTNNB1. The structure of asiatic acid could form one hydrogen bond and two hydrogen bonds with Glu638, and Gln644, respectively, in STAT3. The structure of asiatic acid could form one hydrogen bond with Arg207, Phe250, and Phe256, respectively, in CASP3. Sinapinic acid could form one hydrogen bond with Tyr195 and Tyr227 in TNF. Glycitein could interact with Lys745 and Gly796 through one hydrogen bond in EGFR. The structure of pratensein could interact with Ala36, Gln37, and Asn114 through one hydrogen bond, respectively, and Met111 via two hydrogen bonds in MAPK8. The structure of ferulic acid was linked to Asn51, Tyr139, and Trp162 through one hydrogen bond, respectively, in HSP90AA1. One hydrogen bond could combine pratensein with Tyr139 in HSP90AA1. Glycitein could interact with His90, Tyr385, and Arg513 via one hydrogen bond, respectively, and Ser530 via two hydrogen bonds in PTGS2. Furthermore, some known inhibitors for the core targets were used to validate the molecular docking. It is shown in Table 5. The result indicated the affinities of the inhibitors were comparable to the respective compounds of ZQFZ granule. The structures of the known inhibitors for the core targets are listed in Supplementary Table S4.

4. Discussion

Lung cancer remains the leading cause of cancer death, accounting for 18% of all cancer deaths [1]. To explore the active ingredients and potential mechanisms of ZQFZ granule for the treatment of NSCLC, this study comprehensively investigated the therapeutic effect of ZQFZ granule in NSCLC.

Generally, compounds for network pharmacology studies are obtained from public databases, while the effect of the content may be ignored, resulting in inaccurate predictions. Especially, active compounds were derived from the study that analyzed 15 batches from 3 producers, which significantly increased the credibility of data [10]. Based on Lipinski’s RO5 and SwissTargetPrediction databases, 13 compounds were selected as candidate compounds and 106 potential targets were screened against NSCLC. Moreover, a potential compound from the ZQFZ granule target NSCLC network implied that sinapinic acid, ferulic acid, asiatic acid, pratensein, and glycitein were probably the critical components for NSCLC treatment. Sinapinic acid, one of the most common hydroxycinnamic acids, has an anticancer effect on prostate cancer cells. Treatment with sinapinic acid (1 mM for 72 h) could eliminate half the number of PC-3 and LNCaP cells [13]. Ferulic acid, as a 4-hydroxy-3-methoxycinnamic acid, can inhibit phorbol-12-myristate-13-acetate (PMA)-stimulated invasion of A549 cells at a concentration of ≥100 μM [14]. Asiatic acid, as a natural triterpene, has been proven to enhance the sensitivity of multidrug-resistant A549 cells to cisplatin by downregulating P-glycoprotein (MDR1) and its targets [15]. Glycitein had significant cytotoxic effects on gastric cancer cells by inducing cell apoptosis and G0/G1 phase cell cycle arrest via the ROS-related MAPK/STAT3/NF-κB signaling pathway [16]. Therefore, these compounds play an important role in the treatment of tumors. Notably, there are few studies on the treatment of NSCLC with sinapinic acid, pratensein, and glycitein.

The 106 targets in common between ZQFZ granule and NSCLC were considered potential targets. The top 10 core targets shown in Figure 5(c) included IL-6, SRC, CTNNB1, STAT3, CASP3, TNF, EGFR, MAPK8, HSP90AA1, and PTGS2. These target genes are associated with the regulation of tumorigenesis. IL-6, STAT3, and TNF are involved in immunoactivation and play critical roles in the development of NSCLC. CTNNB1 and EGFR are highly expressed in carcinomas such as lung and breast cancers. SRC is related to endocrine regulation and CASP3 is commonly regarded as the predominant terminal shear enzyme in cell apoptosis. Furthermore, the relationship between the potential compounds and the core targets could be validated by the literature. For instance, sinapinic acid significantly repressed TNF levels and activated CASP3 in acute doxorubicin-induced cardiotoxicity rats [17]. Ferulic acid inhibited UVB-induced TNF and IL-6 protein expression in mice skin tissue, which might be related to the activation of MAPK and NF-kappa B signaling pathways [18]. Asiatic acid induced HT-29 cell apoptosis via CASP3 activation and inhibited the growth and metastasis of breast cancer in mice by downregulating SRC protein expression [19, 20].

The core targets of ZQFZ granule for the treatment of NSCLC were applied to obtain an enrichment map of GO and KEGG pathway analyses via the Metascape database. The results showed BP was involved in oxidative stress, protein phosphorylation, protein-tyrosine kinase signaling pathway, and inflammatory response. CC was primarily concerned with various cell bodies such as membrane, complexes, and cytoplasm. Additionally, MF was mainly associated with protein kinase binding. Furthermore, the main signaling pathways were endocrine resistance, PI3K/Akt, IL-17, and NF-kappa B signaling pathways, while the diseases involved were cancers, atherosclerosis, and measles.

IL-6 is a pleiotropic cytokine with various biological functions in immunity, tissue regeneration, and metabolism. It is highly expressed in lung cancer and negatively correlated with survival [21]. Blocking IL-6 expression can inhibit lung cancer promotion, cell proliferation, angiogenesis markers, and tumor cell-intrinsic STAT3 activation [22]. STAT3 is mainly expressed in naive CD4+ T cells and activated p-STAT3 can modulate the transcriptional activity of target genes associated with tumor cell migration and invasion [23]. IL-6/STAT3 signaling is activated in lung tumorigenesis and metastasis. IL-6 can activate STAT3 signaling and escape host immunity by upregulating and recruiting granulocyte-likemyeloid-derived suppressor cells and type II macrophage polarization [24]. EGFR, a member of the receptor tyrosine kinase ERBB family, is an important oncogene in NSCLC progression. EGFR-TKIs are used for treating NSCLC by modulating the immune microenvironment. However, EGFR-TKI resistance is widely present in NSCLC and promotes immune escape via increased PD-L1 expression and, subsequently, activates the PI3K/AKT/mTOR pathway aberrantly [25, 26]. The PI3K/Akt signaling pathway, a classical mediator that regulates cell growth and inflammation, was found to be activated in this study. There are 17 targets identified in this pathway, and IL-6, TNF, EGFR, MAPK8, and HSP90AA1 are core targets. Tumor necrosis factor (TNF), a cytokine, can combine with extracellular death receptors to activate the apoptosis process and regulate the inflammatory microenvironment via the PI3K/Akt pathway in NSCLC [27, 28]. Upregulation of HSP90AA1 is related to poor overall survival in cancer patients [29]. This could be due to the fact that HSP90AA1 overexpression reduces immune surveillance and resists foreign substances [30].

This is a novel way to study the mechanism of ZQFZ granule in the treatment of NSCLC via a network pharmacologic method. The selected compounds of ZQFZ granule were derived from a recent study based on UPLC. The results may be useful to better understand the multipathway and multitarget regulation of ZQFZ granule against NSCLC.

In total, 13 compounds in the ZQFZ granule and 106 common targets against NSCLC were screened by network pharmacology analysis. GO and KEGG enrichment analyses implied the 13 compounds in the ZQFZ granule could treat NSCLC via multiple NSCLC pathological processes. However, the theoretical observations still need to be validated by clinical studies in the future. This study provided a theoretical basis for the clinical application of ZQFZ granule in NSCLC.

5. Conclusion

This study revealed the pharmacological mechanism of ZQFZ granule in NSCLC based on network pharmacology and molecular docking. GO and KEGG enrichment analyses implied that the 13 compounds in the ZQFZ granule played an important role in the treatment of NSCLC. This study provided a theoretical basis for further research on the clinical application of ZQFZ granule in NSCLC.

Abbreviations

NSCLC:	Non-small cell lung cancer
ZQFZ:	Zhenqi Fuzheng
TCM:	Traditional Chinese medicine
AsR:	Astragali Radix
LLF:	Ligustri Lucidi Fructus
UPLC:	Ultra-performance liquid chromatography
RO5:	Rule of five
Hdon:	Hydrogen bond donors
Hacc:	Hydrogen bond acceptors
LogP:	Lipid-water partition coefficient
Rbon:	Rotatable bonds
PPI:	Protein–protein interaction
MCC:	Maximal clique centrality
GO:	Gene ontology
KEGG:	Kyoto Encyclopedia of Genes and Genomes
BP:	Biological process
CC:	Cellular component
MF:	Molecular function
LUAD:	Lung adenocarcinoma
LUSD:	Lung squamous cell carcinoma
OMIM:	Online Mendelian Inheritance in Man
TTD:	Therapeutic target database.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Yanqing Zhou and Chenxi Wu contributed equally to this work.

Acknowledgments

This study was financially supported by the Research Fund of Medical Specialties Construction Project in the Field of Healthcare in Baoshan District, Shanghai (BSZK-2019-C03) and the Research Fund of Key Specialities Construction Project in Shanghai Baoshan Luodian Hospital, Shanghai (TSZK-2020-03).

Supplementary Materials

Table S1. The 106 common targets of ZQFZ granule against NSCLC. Table S2. Top 5 BP, CC, and MF pathway terms enriched in ZQFZ granule against NSCLC. Table S3. Molecular docking of core targets with top 5 compounds ranked by degree. Table S4. Structures of the known inhibitors for the core targets. (Supplementary Materials)

References

H. Sung, J. Ferlay, R. L. Siegel et al., “Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: A Cancer Journal for Clinicians, vol. 71, no. 3, pp. 209–249, 2021.
View at: Google Scholar
X. L. Su, J. W. Wang, H. Che et al., “Clinical application and mechanism of traditional Chinese medicine in treatment of lung cancer,” Chinese Medical Journal, vol. 133, no. 24, pp. 2987–2997, 2020.
View at: Google Scholar
M. Miao, L. Cao, M. Bai, and G. Chen, “Effect of Zhen Qi Fu Zheng granules on the bone marrow depression model induced by Zidorf,” Saudi Journal of Biological Sciences, vol. 25, no. 2, pp. 220–225, 2018.
View at: Publisher Site | Google Scholar
Z. Xiao, Z. Chen, R. Han et al., “Comprehensive TCM treatments combined with chemotherapy for advanced non-small cell lung cancer: a randomized, controlled trial,” Medicine (Baltimore), vol. 100, no. 18, Article ID e25690, 2021.
View at: Publisher Site | Google Scholar
C. H. Zhang, X. Yang, J. R. Wei et al., “Ethnopharmacology, phytochemistry, pharmacology, toxicology and clinical applications of radix astragali,” Chinese Journal of Integrative Medicine, vol. 27, no. 3, pp. 229–240, 2021.
View at: Publisher Site | Google Scholar
Y. Jung, U. Jerng, and S. Lee, “A systematic review of anticancer effects of radix astragali,” Chinese Journal of Integrative Medicine, vol. 22, no. 3, pp. 225–236, 2016.
View at: Publisher Site | Google Scholar
B. Chen, L. Wang, L. Li et al., “Fructus ligustri lucidi in osteoporosis: a review of its pharmacology, phytochemistry, pharmacokinetics and safety,” Molecules, vol. 22, no. 9, p. 1469, 2017.
View at: Publisher Site | Google Scholar
R. Fabiani, P. Rosignoli, A. De Bartolomeo et al., “Oxidative DNA damage is prevented by extracts of olive oil, hydroxytyrosol, and other olive phenolic compounds in human blood mononuclear cells and HL60 cells,” Journal of Nutrition, vol. 138, no. 8, pp. 1411–1416, 2008.
View at: Publisher Site | Google Scholar
M. Lu, The preparation method of Zhenqi Fuzheng granule, 2016.
Z. Xue, L. Xu, Z. Shang, X. Shi, M. Ye, and X. Qiao, “Discovery of minor quality evaluation marker compounds for Chinese patent medicine products using a two-leveled metabolomics strategy,” Journal of Chromatography A, vol. 1652, Article ID 462354, 2021.
View at: Publisher Site | Google Scholar
C. A. Lipinski, “Lead- and drug-like compounds: the rule-of-five revolution,” Drug Discovery Today: Technologies, vol. 1, no. 4, pp. 337–341, 2004.
View at: Publisher Site | Google Scholar
O. Trott and A. J. Olson, “AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading,” Journal of Computational Chemistry, vol. 31, no. 2, pp. 455–461, 2010.
View at: Publisher Site | Google Scholar
C. Eroglu, E. Avci, H. Vural, and E. Kurar, “Anticancer mechanism of sinapic acid in PC-3 and LNCaP human prostate cancer cell lines,” Gene, vol. 671, pp. 127–134, 2018.
View at: Publisher Site | Google Scholar
C. M. Tsai, G. C. Yen, F. M. Sun, S. F. Yang, and C. J. Weng, “Assessment of the anti-invasion potential and mechanism of select cinnamic acid derivatives on human lung adenocarcinoma cells,” Molecular Pharmaceutics, vol. 10, no. 5, pp. 1890–1900, 2013.
View at: Publisher Site | Google Scholar
Q. Cheng, M. Liao, H. Hu, H. Li, and L. Wu, “Asiatic acid (AA) sensitizes multidrug-resistant human lung adenocarcinoma A549/DDP cells to cisplatin (DDP) via downregulation of P-glycoprotein (MDR1) and its targets,” Cellular Physiology and Biochemistry, vol. 47, no. 1, pp. 279–292, 2018.
View at: Publisher Site | Google Scholar
Y. Q. Zang, Y. Y. Feng, Y. H. Luo et al., “Glycitein induces reactive oxygen species-dependent apoptosis and G0/G1 cell cycle arrest through the MAPK/STAT3/NF-κB pathway in human gastric cancer cells,” Drug Development Research, vol. 80, no. 5, pp. 573–584, 2019.
View at: Publisher Site | Google Scholar
Y. A. Bin Jardan, M. A. Ansari, M. Raish et al., “Sinapic acid ameliorates oxidative stress, inflammation, and apoptosis in acute doxorubicin-induced cardiotoxicity via the NF-κB-mediated pathway,” BioMed Research International, vol. 2020, Article ID 3921796, 10 pages, 2020.
View at: Publisher Site | Google Scholar
K. Ambothi, N. R. Prasad, and A. Balupillai, “Ferulic acid inhibits UVB-radiation induced photocarcinogenesis through modulating inflammatory and apoptotic signaling in Swiss albino mice,” Food and Chemical Toxicology, vol. 82, pp. 72–78, 2015.
View at: Publisher Site | Google Scholar
P. Bunpo, K. Kataoka, H. Arimochi et al., “Inhibitory effects of asiatic acid and CPT-11 on growth of HT-29 cells,” Journal of Medical Investigation, vol. 52, no. 1-2, pp. 65–73, 2005.
View at: Publisher Site | Google Scholar
M. Tian, K. Chen, J. Huang et al., “Asiatic acid inhibits angiogenesis and vascular permeability through the VEGF/VEGFR2 signaling pathway to inhibit the growth and metastasis of breast cancer in mice,” Phytotherapy Research, vol. 35, no. 11, pp. 6389–6400, 2021.
View at: Publisher Site | Google Scholar
L. Barrera, E. Montes-Servin, A. Barrera et al., “Cytokine profile determined by data-mining analysis set into clusters of non-small-cell lung cancer patients according to prognosis,” Annals of Oncology, vol. 26, no. 2, pp. 428–435, 2015.
View at: Publisher Site | Google Scholar
M. S. Caetano, H. Zhang, A. M. Cumpian et al., “IL6 blockade reprograms the lung tumor microenvironment to limit the development and progression of K-ras-Mutant lung cancer,” Cancer Research, vol. 76, no. 11, pp. 3189–3199, 2016.
View at: Publisher Site | Google Scholar
J. E. Yeh, S. Kreimer, S. R. Walker et al., “Granulin, a novel STAT3-interacting protein, enhances STAT3 transcriptional function and correlates with poorer prognosis in breast cancer,” Genes Cancer, vol. 6, no. 3-4, pp. 153–168, 2015.
View at: Publisher Site | Google Scholar
B. Jing, T. Wang, B. Sun et al., “IL6/STAT3 signaling orchestrates premetastatic niche formation and immunosuppressive traits in lung,” Cancer Research, vol. 80, no. 4, pp. 784–797, 2020.
View at: Publisher Site | Google Scholar
S. Peng, R. Wang, X. Zhang et al., “EGFR-TKI resistance promotes immune escape in lung cancer via increased PD-L1 expression,” Molecular Cancer, vol. 18, no. 1, p. 165, 2019.
View at: Publisher Site | Google Scholar
C. Fumarola, M. A. Bonelli, P. G. Petronini, and R. R. Alfieri, “Targeting PI3K/AKT/mTOR pathway in non small cell lung cancer,” Biochemical Pharmacology, vol. 90, no. 3, pp. 197–207, 2014.
View at: Publisher Site | Google Scholar
R. Zhang, Y. Dong, M. Sun et al., “Tumor-associated inflammatory microenvironment in non-small cell lung cancer: correlation with FGFR1 and TLR4 expression via PI3K/Akt pathway,” Journal of Cancer, vol. 10, no. 4, pp. 1004–1012, 2019.
View at: Publisher Site | Google Scholar
F. Balkwill, “TNF-alpha in promotion and progression of cancer,” Cancer and Metastasis Reviews, vol. 25, no. 3, pp. 409–416, 2006.
View at: Publisher Site | Google Scholar
G. Fan, Y. Tu, N. Wu, and H. Xiao, “The expression profiles and prognostic values of HSPs family members in head and neck cancer,” Cancer Cell International, vol. 20, no. 1, p. 220, 2020.
View at: Publisher Site | Google Scholar
X. Xiang, X. M. You, and L. Q. Li, “Expression of HSP90AA1/HSPA8 in hepatocellular carcinoma patients with depression,” OncoTargets and Therapy, vol. 11, pp. 3013–3023, 2018.
View at: Publisher Site | Google Scholar

Copyright

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

PDF Download Citation

Download other formats

Order printed copies

Views

328

Downloads

493

Citations