Network pharmacology-based identification of key pharmacological pathways of Yin-Huang-Qing-Fei capsule acting on chronic bronchitis.

For decades in China, the Yin-Huang-Qing-Fei capsule (YHQFC) has been widely used in the treatment of chronic bronchitis, with good curative effects. Owing to the complexity of traditional Chinese herbal formulas, the pharmacological mechanism of YHQFC remains unclear. To address this problem, a network pharmacology-based strategy was proposed in this study. At first, the putative target profile of YHQFC was predicted using MedChem Studio, based on structural and functional similarities of all available YHQFC components to the known drugs obtained from the DrugBank database. Then, an interaction network was constructed using links between putative YHQFC targets and known therapeutic targets of chronic bronchitis. Following the calculation of four topological features (degree, betweenness, closeness, and coreness) of each node in the network, 475 major putative targets of YHQFC and their topological importance were identified. In addition, a pathway enrichment analysis based on the Kyoto Encyclopedia of Genes and Genomes pathway database indicated that the major putative targets of YHQFC are significantly associated with various pathways involved in anti-inflammation processes, immune responses, and pathological changes caused by asthma. More interestingly, eight major putative targets of YHQFC (interleukin [IL]-3, IL-4, IL-5, IL-10, IL-13, FCER1G, CCL11, and EPX) were demonstrated to be associated with the inflammatory process that occurs during the progression of asthma. Finally, a molecular docking simulation was performed and the results exhibited that 17 pairs of chemical components and candidate YHQFC targets involved in asthma pathway had strong binding efficiencies. In conclusion, this network pharmacology-based investigation revealed that YHQFC may attenuate the inflammatory reaction of chronic bronchitis by regulating its candidate targets, which may be implicated in the major pathological processes of the asthma pathway.

Introduction

Chronic bronchitis, a very common disease associated with chronic inflammation, is caused by a variety of etiological factors, such as air pollution, smoking, and pathogen infection.1,2 It has been predicted to be the fourth leading cause of death by 2030.3,4 There are two stages of episodes for chronic bronchitis: recurrent reaction and acute exacerbation.5 The main symptoms of recurrent reaction include cough, asthma, and expectoration, while acute exacerbation is diagnosed when patients present with increased sputum volume, cough, and worsening dyspnea.5 Although current therapies, including bronchodilators, anti-tussive agents, steroids, and leukotriene (LT) receptor antagonists, can alleviate some major symptoms of chronic bronchitis,6 the treatments based on these drugs often lead to a series of adverse reactions.7 Therefore, it is of great clinical significance to develop more efficient therapeutic strategies to both improve the curative therapy of patients with chronic bronchitis and reduce the incidence of side effects of current treatments.2

Traditional Chinese medicine (TCM), as a major component of complementary and alternative medicine remedies, has attracted increasing attention worldwide because of its satisfactory clinical efficacy.8 It has also played an important role in preventing people from contracting chronic bronchitis for centuries.9 In recent years, the Yin–Huang–Qing–Fei capsule (YHQFC), based on the famous classical formula for the Ma-Xing-Shi-Gan Decoction from the TCM book Treatise on Febrile Diseases (Shang Han Lun), has become a typical formula used to treat chronic bronchitis and has produced satisfactory therapeutic effects clinically. One study reported that the effective rate of treatment for patients with acute exacerbation of chronic obstructive pulmonary disease in the group receiving YHQFC was 97.5%, compared to 85.0% in the control group.10 YHQFC contains 14 Chinese herbs including Ephedra sinica Stapf (Mi Ma Huang [MMH]), Lepidium apetalum Willd. (Bei Ting Li Zi [BTLZ]), and Semen Armeniacae Amarum (Ku Xing Ren [KXR]). Zhou et al11 identified 54 compounds in YHQFC using the high performance liquid chromatography coupled with quadrupole time-of-flight mass spectrometry method and preliminarily explored the bioactive components of this formula. Additionally, a serum containing YHQFC showed a strong anti-respiratory syncytial virus effect in vitro and was found to inhibit influenza virus replication in cells to a certain degree.12,13 Because YHQFC contains a great number of chemical compounds and regulates diverse targets, precisely determining the pharmacological mechanisms of its therapeutic actions and deciphering the relationships between the herbs and diseases are challenges.

Network pharmacology is emerging as a promising strategy, one that is closely related to the application of multiple omics- and systems biology-based technologies.14 It is a valuable tool for achieving a holistic view of and comprehensive and systematic insight into the mechanisms of multi-ingredient medicine.15 Various molecular networks of complex ingredients and multilevel target-based protein and gene interactions have been constructed for predicting their functions and promoting discovery of active compounds.16 Because TCM herbal formulas are considered to be multi-component and multi-target therapeutics, the methodologies of network pharmacology are suitable for pursuing a priori knowledge about the combinations of rules embedded in these formulas.17 Thus, the application of network pharmacology to TCM provides new opportunities to understand the interactions between active compounds and relevant targets, which in turn highlight the mechanisms of action.8,18 Using this context, the aim of this study was to develop a comprehensive network pharmacology-based approach to investigate the pharmacological mechanisms of YHQFC acting on chronic bronchitis. Figure 1 depicts a flowchart of the experimental procedures of our study.

Figure 1

Whole framework of this study based on network pharmacology for deciphering pharmacological mechanisms of YHQFC acting on chronic bronchitis.

Abbreviations: YHQFC, Yin–Huang–Qing–Fei capsule; OMIM, Online Mendelian Inheritance in Man; PPI, protein–protein interaction.

Materials and methods

Ethics

All experimental procedures were carried out according to the Guideline for Experimentation of the China Academy of Chinese Medical Sciences, and the protocol was reviewed and approved by the Ethics Committee of the institution.

Data preparation

Construction of the chemical information database of YHQFC

The chemical components of each herb contained in YHQFC – Ephedra sinica Stapf (MMH), Lepidium apetalum Willd. (BTLZ), Semen Armeniacae Amarum (KXR), Fritillaria thunbergii Miq. (Zhe Bei Mu [ZBM]), Eriobotrya japonica (Thunb.) Lindl. (Pi Pa Ye [PPY]), Gypsum Fibrosum (Sheng Shi Gao [SSG]), Acorus tatarinowii Schott (Shi Chang Pu [SCP]), Isatis indigotica Fort. (Da Qing Ye [DQY]), Artemisia rupestris L. (Xin Jiang Yi Zhi Hao [XJYZH]), Dioscorea nipponica Makino (Chuan Shan Long [CSL]), Schisandra chinensis (Turcz.) Baill. (Wu Wei Zi [WWZ]), Ginkgo biloba L (Yin Xing Ye [YXY]), Citrus aurantium L. (Zhi Shi [ZS]), and Glycyrrhiza uralensis Fisch. (Gan Cao [GC]) – were collected from the Chinese Academy of Sciences chemistry database (http://www.organchem.csdb.cn/scdb/main/slogin.asp, updated on July 20, 2015). This database displays chemical information, including chemical and crystal structures, spectra, reactions, syntheses, and thermo-physical properties.19 The available structural information for YHQFC included 29 compounds in MMH, 12 compounds in BTLZ, 7 compounds in KXR, 34 compounds in ZBM, 24 compounds in PPY, 1 compound in SSG (as CaSO4⋅2H2O), 18 compounds in SCP, 14 compounds in DQY, 35 compounds in XJYZH, 7 compounds in CSL, 89 compounds in WWZ, 55 compounds in YXY, 78 compounds in ZS, and 190 compounds in GC. Detailed information on the constituent compounds of each herb contained in YHQFC is provided in Table S1. The molecular files of all the constituent compounds were downloaded from the ChemSpider database (http://www.chemspider.com/, updated on November 29, 2015) and were saved in .mol format.

Known therapeutic targets of drugs in the treatment of chronic bronchitis

The known therapeutic targets of drugs used in the treatment of chronic bronchitis were acquired from two sources: DrugBank20 (http://www.drugbank.ca/, version 4.3) and the Online Mendelian Inheritance in Man (OMIM) database21 (http://www.omim.org/, updated on April 21, 2016). Using drug–target interactions mentioned in DrugBank, only drugs approved by the Food and Drug Administration (FDA) were selected for the treatment of chronic bronchitis and whose targets are human genes/proteins. In total, we chose 238 known therapeutic targets for treating chronic bronchitis were chosen. In addition, the main symptoms of chronic bronchitis, such as “cough,” “asthma,” and “bronchitis,” were used as key words when searching the OMIM database, and the resulting 286 known therapeutic targets were enrolled in this study. Detailed information about these known therapeutic targets is provided in Table S2. After removing redundant entries, 381 known therapeutic targets for the treatment of chronic bronchitis were used for the data analysis.

Prediction of putative targets of YHQFC

A drug similarity search tool MedChem Studio (version 3.0; Simulations Plus, Inc., Lancaster, CA, USA, 2012),22 effective for quick identification of all molecules with structure similarities, was used to identify known drugs that are structurally similar to the chemical components of YHQFC. Then, the therapeutic targets of these known drugs, obtained from the DrugBank database, were considered as putative targets of YHQFC. To improve the reliability of the target prediction results, only drugs with a high structural similarity score (>0.80, moderately similar to very similar) were selected.

Protein–protein interaction data

Protein–protein interaction (PPI) data were derived from eight existing PPI databases, including HAPPI,23 Reactome,24 OPHID,25 InAct,26 HPRD,27 MINT,28 DIP,29 and PDZBase.30 Detailed information about these PPI databases is provided in Table S3. The differing ID types of the proteins were converted to UniProt IDs.

Network construction and analysis

Putative YHQFC target-known therapeutic targets of the chronic bronchitis network and putative YHQFC target networks were constructed to understand the associations between the herbs and chemical compounds containing YHQFC, putative targets of YHQFC, and known therapeutic targets of chronic bronchitis. Network visualization was performed using Navigator software (version 2.2.1, Toronto, ON, USA).

Putative YHQFC target-known therapeutic targets of the chronic bronchitis network

Putative YHQFC target-known therapeutic targets of the chronic bronchitis network were constructed using the links between putative targets of YHQFC and known therapeutic targets of chronic bronchitis. As previously reported,31 if the degree of a node is more than twofold the median degree of all the nodes in a network, the node may function as a hub. Subsequently, the hub network was constructed based on the direct interactions among hubs. Moreover, four topological properties, ”degree”, “betweenness”, “closeness”, and “coreness”, were calculated to screen the putative targets for topological importance. The definitions of the four topological features listed earlier are provided in Supplementary File S1. Only the topological features higher than the corresponding median values were identified as the major hubs.

Pathway enrichment analysis

To clarify the pathways that are involved in putative YHQFC targets, a pathway enrichment analysis was performed using the Database Visualization and Integrated Discovery software32 (DAVID, http://david.abcc.ncifcrf.gov/home.jsp, version 6.7) and based on the pathway data obtained from the Kyoto Encyclopedia of Genes and Genomes database33 (KEGG, http://www.genome.jp/kegg/, updated on April 18, 2016).

Molecular docking simulation

To verify the binding affinity of candidate targets to constituent compounds of YHQFC, a molecular docking simulation was performed using the program LibDock implemented in Discovery Studio 2.5 (DS 2.5). All the crystal structures of the candidate targets, including interleukin (IL)-3, IL-5, FCER1G, and EPX, were directly downloaded from the RCSB protein data bank (http://www.pdb.org/, updated on October 21, 2016) and were carefully examined for their resolutions. To evaluate the binding affinity of each candidate target to the corresponding compound, a docking score was calculated by the customizable scoring function of LibDock. The docking scores >72 (the median value of all the docking scores) indicated a strong binding affinity of candidate targets to their corresponding compounds.

Results and discussion

Putative targets of YHQFC

A total of 593 constituent compounds of the 14 herbs contained in YHQFC were collected, and 9,199 genes were predicted as the putative targets of YHQFC, including 464 of SCP, 901 of XJYZH, 607 of ZS, 53 of CSL, 1135 of MMH, 491 of PPY, 513 of ZBM, 1535 of YXY, 950 of GC, 197 of SSG, 881 of DQY, 119 of BTLZ, 463 of WWZ, and 890 of KXR. Detailed information about the putative targets of YHQFC is provided in Table S4.

Notably, there were different numbers of common putative targets among the herbs contained in YHQFC, as shown in Table 1, suggesting that these herbs might have several interactions in the course of treatment. Moreover, >50% of the putative targets of PPY were common to those of the other herbs, implying a close relationship between PPY and other herbs in the formula.

Table 1

Potential target overlaps among 14 herbs of YHQFC

Herbs	SCP(464)	XJYZH(901)	ZS(607)	CSL(53)	MMH(1135)	PPY(491)	ZBM(513)	YXY(1535)	GC(950)	SSG(197)	DQY(881)	BTLZ(119)	WWZ(463)	KXR(890)
SCP (464)	–	224	305	21	244	129	139	385	279	57	280	44	293	275
XJYZH (901)	224	–	375	25	661	412	294	808	482	69	321	118	299	510
ZS (607)	305	375	–	14	392	255	230	518	417	70	202	110	393	327
CSL (53)	21	25	14	–	27	11	53	34	53	8	53	1	9	33
MMH (1135)	244	661	392	27	–	346	241	908	535	91	514	73	274	405
PPY (491)	129	412	255	11	346	–	190	426	287	48	165	112	219	270
ZBM (513)	139	294	230	53	241	190	–	443	297	54	231	112	178	348
YXY (1535)	385	808	518	34	908	426	443	–	719	118	622	119	418	690
GC (950)	279	482	417	53	535	287	297	719	–	88	330	104	309	451
SSG (197)	57	69	70	8	91	48	54	118	88	–	44	13	57	64
DQY (881)	280	321	202	53	514	165	231	622	330	44	–	50	163	315
BTLZ (119)	44	118	110	1	73	112	112	119	104	13	50	–	107	115
WWZ (463)	293	299	393	9	274	219	178	418	309	57	163	107	–	251
KXR (890)	275	510	327	33	405	270	348	690	451	64	315	115	251	–

Note: Numbers shown within parentheses show total putative targets of each herb.

Abbreviations: YHQFC, Yin–Huang–Qing–Fei capsule; SCP, Shi Chang Pu (Acorus tatarinowii Schott); XJYZH, Xin Jiang Yi Zhi Hao (Artemisia rupestris L.); ZS, Zhi Shi (Citrus aurantium L.); CSL, Chuan Shan Long (Dioscorea nipponica Makino); MMH, Mi Ma Huang (Ephedra sinica Stapf); PPY, Pi Pa Ye (Eriobotrya japonica (Thunb.) Lindl.); ZBM, Zhe Bei Mu (Fritillaria thunbergii Miq.); YXY, Yin Xing Ye (Ginkgo biloba L); GC, Gan Cao (Glycyrrhiza uralensis Fisch.); SSG, Sheng Shi Gao (Gypsum Fibrosum); DQY, Da Qing Ye (Isatis indigotica Fort.); BTLZ, Bei Ting Li Zi (Lepidium apetalum Willd.); WWZ, Wu Wei Zi (Schisandra chinensis [Turcz.] Baill.); KXR, Ku Xing Ren (Semen Armeniacae Amarum).

Pharmacological mechanisms of YHQFC acting on chronic bronchitis

The putative YHQFC target-known therapeutic target of the chronic bronchitis network consisted of 10,828 nodes and 73,659 edges. Using the twofold median value of node degree (4) in this network as a cutoff point, 3,667 nodes were identified as hubs. Detailed information on these hubs is provided in Table S5.

Following the construction of a hub interaction network and the calculation of the four topological features, ”degree”, “betweenness”, “closeness”, and “coreness”, of each hub, 475 major hubs, including 395 putative YHQFC targets and 80 known therapeutic targets of chronic bronchitis, were screened; detailed information on these targets is provided in Table S6. To investigate the biological functions of these major hubs, a pathway enrichment analysis was performed. Consequently, the major hubs were significantly associated with various physiological processes, including inflammation processes and immune suppression (such as asthma pathways, cytokine–cytokine receptor interactions, and T cell receptor signaling pathways), vascular smooth muscle contraction, and neuroactive ligand–receptor interaction (Table 2). Cytokines are soluble mediators (polypeptides) that act as messengers for the immune system, among other biologic functions, and are important for leukocyte growth and differentiation, as well as activation and migration. The overexpression of inflammatory cytokines or decreased levels of anti-inflammatory cytokines can lead to inflammatory disorders.34 In human beings, allergic asthma is primarily initiated by a type I hypersensitivity reaction, which represents increased susceptibility to producing immunoglobulin E (IgE) in response to external allergens, the secretion and differentiation of which are dependent on CD4+ helper T cells (Th2 type).35,36 Activated Th2 cells may produce a number of cytokines (eg, IL-3, IL-4, IL-5, and IL-13), which in turn play roles in various processes, such as promoting production of IgE cells by B cells (IL-4 and IL-13), growth of mast cells (IL-4), and survival of eosinophils (IL-3 and IL-5).37 IL-4 and IL-13 also stimulate epithelial cells to produce transforming growth factor alpha, which may cause mucosal metaplasia and fibroblast proliferation.38 Proinflammatory mediators, including tumor necrosis factors and granulocyte macrophage colony-stimulating factors, stimulate the expression of vascular adhesion molecules on endothelial cells, which results in an increased influx of inflammatory leukocytes into the bronchial tree.39 The pathology of inflammatory processes in chronic airway diseases such as chronic bronchitis is directly associated with increased generation of reactive oxygen species (ROS) in the lungs.40 Various enzymatic and nonenzymatic antioxidants are available in the lungs to counteract ROS-mediated damage, which may lead to increased vascular permeability, mucus hypersecretion, smooth muscle contraction, and epithelial shedding.41–45 The regulation of vascular smooth muscle contractions contributes to the control of blood pressure and airway remodeling.46,47 It is reported that increase in the local concentration of inflammatory neuropeptides may stimulate bronchial fibers, which rapidly adapts receptors and provokes the cough reflex.48 The neuroactive ligand–receptor pathway includes many G-protein-coupled receptors, which are often located on the cell surface and involved in the transduction of an endogenous signal into an intracellular response.49 G-protein-coupled receptors regulate vital processes by controlling the expression and function of seven transmembrane receptors, which play a particularly important role in some widespread neurological disorders.19 For example, the activation of the A1 adenosine receptor in allergic rabbit airway smooth muscle causes the production of inositol 1,4,5-trisphosphate [Ins(1,4,5)P3] via a pertussis toxin-sensitive G protein-coupled phospholipase C. This signaling mechanism may be involved in the generation of contractile responses. It is hypothesized that this process may contribute to adenosine-induced bronchoconstriction in allergic asthma.50

Table 2

Main KEGG pathways significantly relating to major hubs

KEGG pathways	Counts	P-value
hsa04664:Fc epsilon RI signaling pathway	34	2.36×10−16
hsa04662:B cell receptor signaling pathway	31	3.96×10−14
hsa04660:T cell receptor signaling pathway	34	1.41×10−11
hsa04620:Toll-like receptor signaling pathway	29	6.16×10−9
hsa04080:Neuroactive ligand–receptor interaction	46	8.02×10−7
hsa04621:NOD-like receptor signaling pathway	19	1.62×10−6
hsa04270:Vascular smooth muscle contraction	26	3.58×10−6
hsa04630:Jak-STAT signaling pathway	31	8.62×10−6
hsa04650:Natural killer cell-mediated cytotoxicity	27	2.89×10−5
hsa04060:Cytokine–cytokine receptor interaction	36	0.003305
hsa05310:Asthma	8	0.008380
hsa05130:Pathogenic Escherichia coli infection	10	0.045010

Abbreviation: KEGG, Kyoto Encyclopedia of Genes and Genomes.

A herb-major hub-pathway network was constructed to visualize the relationships among herbs, major targets, and disease (Figure 2). This network illustrated that the major putative targets of YHQFC are mostly involved in pathways related to anti-inflammation processes and the immune system, which may be crucial elements during the attenuation of chronic bronchitis. Interestingly, we also found that YHQFC may exert its therapeutic effects on chronic bronchitis by regulating the pathways related to the contraction of airway vascular smooth muscle. The relevant nerve center regulation induced by the major putative targets of YHQFC was another pathway to calm cough and asthma.

Figure 2

Herb-major hub-pathway network.

Notes: Green diamonds represent each herb in YHQFC; round nodes represent putative targets of constituent compounds of YHQFC; nodes marked with red rings represent known therapeutic targets in the treatment of chronic bronchitis. Blue hexagons represent main pathways from enrichment analysis of major targets. Edges represent interactions among herbs, putative targets, and known therapeutic targets for the treatment of chronic bronchitis.

Abbreviations: YHQFC, Yin–Huang–Qing–Fei capsule; SCP, Shi Chang Pu (Acorus tatarinowii Schott); XJYZH, Xin Jiang Yi Zhi Hao (Artemisia rupestris L.); ZS, Zhi Shi (Citrus aurantium L.); CSL, Chuan Shan Long (Dioscorea nipponica Makino); MMH, Mi Ma Huang (Ephedra sinica Stapf); PPY, Pi Pa Ye (Eriobotrya japonica (Thunb.) Lindl.); ZBM, Zhe Bei Mu (Fritillaria thunbergii Miq.); YXY, Yin Xing Ye (Ginkgo biloba L); GC, Gan Cao (Glycyrrhiza uralensis Fisch.); SSG, Sheng Shi Gao (Gypsum Fibrosum); DQY, Da Qing Ye (Isatis indigotica Fort.); BTLZ, Bei Ting Li Zi (Lepidium apetalum Willd.); WWZ, Wu Wei Zi (Schisandra chinensis [Turcz.] Baill.); KXR, Ku Xing Ren (Semen Armeniacae Amarum).

YHQFC attenuates the pathological changes of chronic bronchitis via regulating the asthma pathway

Four major putative targets of YHQFC (IL-3, IL-5, the FcεRI receptors [FCER1G], and eosinophil peroxidase [EPX]) and four known therapeutic targets of chronic bronchitis (IL-4, IL-10, IL-13, and a member of the eotaxin subfamily of CC chemokines [CCL11]) were demonstrated to participate in the asthma pathway, which has been shown to be involved in the occurrence and progression of chronic bronchitis. As shown in Figure 3, the asthma pathway is mainly composed of inflammatory processes. According to the clusters of major hubs contained in the asthma pathway, YHQFC may be of vital importance for regulating these cytokines and chemokines.

Figure 3

Illustration of the asthma pathway induced by the major putative YHQFC targets and known therapeutic targets of chronic bronchitis.

Abbreviations: YHQFC, Yin–Huang–Qing–Fei capsule; IL, interleukin; APC, antigen-presenting cell; MHC, major histocompatibility complex; Th, T helper; TCR, T cell receptor; BCR, B cell receptor; IgE, immunoglobulin E; LTD4, leukotriene D4; PAF, platelet-activating factor; ECP, eosinophil cationic protein; MBP, myelin basic protein.

Among the major putative targets of YHQFC, IL-4 and IL-13 both play important roles in the class switching of B cells to IgE synthesis, as well as in the recruitment of mast cells, while IL-4 alone is crucial for maintaining the Th2 cell phenotype. Several animal studies have shown that blocking the production or inhibiting the effects of IL-4 has profound effects on the allergic phenotype.51 It is well established that the crosslinking of FcεRI receptors by allergen-bound IgE induces the release of inflammatory mediators (eg, histamine), which, in turn, causes the immediate symptoms of inflammatory reaction within 30 min.52 The suppression of FCER1G and down-regulation of the expression of FcεRI receptors may have effects on relieving the immediate reaction of chronic bronchitis. IL-10 is the main inhibitory cytokine produced by regulatory T (TReg) cells that inhibits the expression of many proinflammatory cytokines and chemokines, as well as proinflammatory enzymes.53 It has been reported that IL-10-deficient mice had increased airway inflammation after allergen challenges.54 The presence of IL-13, as well as epidermal, neural, vascular, and fibroblast growth factors, leads to mucous metaplasia and remodeling of the airway walls.55 The numerous functions of IL-13 in regulating IgE production, eosinophilic inflammation, airway smooth muscle hyperplasia, the induction of goblet cell hyperplasia with mucus production, and the recruitment of monocytes, macrophages, and T cells into the airway spaces make it a key therapeutic target in chronic bronchitis.56 Selective blocking of IL-13 has been achieved in mice using a soluble form of IL-13Rα2, which has a high affinity with IL-13, and has led to the reversal of airway hyper-responsiveness and mucus production in allergen-exposed sensitized mice.57 IL-3 and IL-5 can be released by Th2 cells and mast cells when they are stimulated.58,59 The function of IL-3 and IL-5 is to promote the maturation of eosinophils, and it is evident that the activation of eosinophils can lead to the production of ROS, including EPX, which are implicated in protein oxidation in the bronchus airway and the late asthmatic response.60–62 CCL11 (eotaxin) is a member of the eotaxin subfamily of CC chemokines and is a potent chemoattractant for eosinophils.63,64 It acts predominantly through the eotaxin receptor, CCR3, which has been identified on eosinophils, mast cells, Th2 cells, and airway epithelial cells.65 It is evident that CCL11 can be upregulated in a variety of diseases associated with eosinophilic infiltrates, including chronic bronchitis.66 Down-regulation of CCL11 and CCR3 receptor antagonists may be useful for blocking eosinophil entry and degranulation in chronic bronchitis.67

Importantly, among the four major putative targets involved in the asthma pathway, IL-3 is the putative target of 2-methyl-1,3,6-trihydroxyanthraquinone (constituent compound from GC) and angeloylgomisin Q (constituent compound from WWZ). FCER1G is the putative target of phenylalanine (constituent compound from YXY), mandelic acid (constituent compound from KXR), choline (constituent compound from ZBM), and isolobelanine (constituent compound from GC). EPX is the putative target of euasarone (constituent compound from SCP) and N-methyltyramine (constituent compound from ZS). IL-5 is the putative target of multiple constituent compounds, including naringenin (from ZS) and isoginkgetin (from YXY). These data may provide more evidence for the identification of the active compounds of YHQFC that act on chronic bronchitis.

Molecular docking validation

It is evident that the computational docking simulation is an invaluable tool in drug discovery and design. This technique is structure-based and could help researchers discover the interactions between the components of TCM and network targets.68 DS LibDock is an efficient and high-throughput module in which the ligand will be aligned in the active site. LibDock consensus score can be used to evaluate the binding efficiency of a ligand, which is based on the geometry of a candidate ligand docked into a target receptor structure using empirical functions.69 In this study, a docking score >72 indicated a strong binding capacity between the constituent compounds of YHQFC and the candidate targets. And 17 pairs of compound–candidate effector molecule interactions were of strong binding affinity. Detailed information about the docking scores of the candidate targets binding to their corresponding compounds contained in YHQFC is described in Table S7. To investigate the regulation of these compounds on the corresponding candidate targets, known drugs were collected with similar structures to these compounds and their experimentally validated interactions with the corresponding candidate targets using Medchem Studio. As drugs with similar structures often possess similar functions, the regulatory effects of compounds on candidate effectors of YHQFC may be similar to those of the corresponding similar known drugs. As a result, 29 pairs of compound–candidate effector molecule interactions have positive effects and 4 pairs have negative effects. Detailed information about the regulatory effects of similar known drugs on candidate targets of YHQFC is described in Table S8. These findings should be validated by further biological experiments.

Conclusion

Chinese herbal formulas are one of the most important types of complementary and alternative medicine and contribute greatly to respiratory disease therapy. This study proposed and applied network pharmacology-based analysis to suggest that YHQFC may attenuate the inflammatory response of chronic bronchitis by regulating its candidate targets implicated into the asthma pathway. Further experimental validations of these prediction results are required in our future studies.