Mechanism of Fei-Xian Formula in the Treatment of Pulmonary Fibrosis on the Basis of Network Pharmacology Analysis Combined with Molecular Docking Validation

Chen, Xiao-Li; Tang, Cheng; Xiao, Qing-Ling; Pang, Zhong-Hua; Zhou, Dan-Dan; Xu, Jin; Wang, Qi; Zhao, Ya-xi; Zhu, Qi-Yong

doi:https://doi.org/10.1155/2021/6658395

Evidence-Based Complementary and Alternative Medicine

On this page

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

Special Issue

Innovative Methods and Technologies for Efficacy Evaluation of Traditional Medicine

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 6658395 | https://doi.org/10.1155/2021/6658395

Mechanism of Fei-Xian Formula in the Treatment of Pulmonary Fibrosis on the Basis of Network Pharmacology Analysis Combined with Molecular Docking Validation

Xiao-Li Chen,¹Cheng Tang,¹Qing-Ling Xiao,¹Zhong-Hua Pang,²Dan-Dan Zhou,³Jin Xu,¹Qi Wang,¹Ya-xi Zhao,¹and Qi-Yong Zhu¹

Academic Editor: Hongcai Shang

Received16 Nov 2020

Revised21 Jun 2021

Accepted16 Jul 2021

Published04 Aug 2021

Abstract

Objective. This study aimed to clarify the mechanism of Fei-Xian formula (FXF) in the treatment of pulmonary fibrosis based on network pharmacology analysis combined with molecular docking validation. Methods. Firstly, ingredients in FXF with pharmacological activities, together with specific targets, were identified based on the BATMA-TCM and TCMSP databases. Then, targets associated with pulmonary fibrosis, which included pathogenic targets as well as those known therapeutic targets, were screened against the CTD, TTD, GeneCards, and DisGeNet databases. Later, Cytoscape was employed to construct a candidate component-target network of FXF for treating pulmonary fibrosis. In addition, for nodes within the as-constructed network, topological parameters were calculated using CytoHubba plug-in, and the degree value (twice as high as the median degree value for all the nodes) was adopted to select core components as well as core targets of FXF for treating pulmonary fibrosis, which were subsequently utilized for constructing the core network. Furthermore, molecular docking study was carried out on those core active ingredients together with the core targets using AutoDock Vina for verifying results of network pharmacology analysis. At last, OmicShare was employed for enrichment analysis of the core targets. Results. Altogether 12 active ingredients along with 13 core targets were identified from our constructed core component-target network of FXF for the treatment of pulmonary fibrosis. As revealed by enrichment analysis, the 13 core targets mostly concentrated in regulating biological functions, like response to external stimulus (from oxidative stress, radiation, UV, chemical substances, and virus infection), apoptosis, cell cycle, aging, immune process, and protein metabolism. In addition, several pathways, like IL-17, AGE-RAGE, TNF, HIF-1, PI3K-AKT, NOD-like receptor, T/B cell receptor, and virus infection-related pathways, exerted vital parts in FXF in the treatment of pulmonary fibrosis. Conclusions. FXF can treat pulmonary fibrosis through a “multicomponent, multitarget, and multipathway” mean. Findings in this work lay foundation for further exploration of the FXF mechanism in the treatment of pulmonary fibrosis.

1. Introduction

Pulmonary fibrosis may take place in a variety of clinical settings and may threaten human life [1, 2]. It is featured by the disturbed cell homeostasis and changed cell composition in peripheral lungs, thus resulting in excess extracellular matrix (ECM) accumulation and lung dysfunction [3]. Over the last 10 years, investigators have identified that some molecular and cellular signal transduction pathways are related to the pulmonary fibrosis pathogenesis, which promotes to identify more novel therapeutic targets [4, 5]. Numerous treatments have been utilized to treat pulmonary fibrosis, but few of them elevate patient survival or improve their quality of life [6]. Nintedanib and pirfenidone have been approved to be used to treat pulmonary fibrosis, yet they have serious diverse reactions (such as gastrointestinal symptoms, photosensitivity, and abnormal liver function) [7–10]. Such phenomenon may be ascribed to the complicated regulatory networks involved in pulmonary fibrosis, which suppress or promote associated target genes or pathway expression [11–13]. It is increasingly suggested that epigenetic, genetic, or proteomic factors play important roles in these regulatory networks in pulmonary fibrosis [14–16]. Besides, the networks modulate fibrosis-related gene expression through activating or inactivating the relevant pathways, but not through one single pathway [17].

Traditional Chinese medicine (TCM) exerts multi-target effects and has been used to treat pulmonary fibrosis by regulating oxidant stress, inflammation, and so on. The TCM pathogenesis of pulmonary fibrosis is Yin deficiency and blood stasis. Thus, the treatment principle is dominated by promoting blood circulation to remove blood stasis, enriching Yin and nourishing the lung [18]. Fei-Xian formula (FXF), the cipher prescription used to treat pulmonary fibrosis at our institution, is made up of 7 herbs, Salviae Miltiorrhizae Radix et Rhizoma (Danshen, DS), Persicae Semen (Taoren, TR), Hirudo (Shuizhi, SZ), Paeoniae Radix Rubra (Chishao, CS), Asparagi Radix (Tiandong, TD), Canarii Fructus (Qingguo, QG), and Oroxyli Semen (Muhudie, MHD). FXF is effective in relieving the pulmonary fibrosis symptoms, improving lung function (evaluated by vital capacity max (VC Max), total lung capacity (TLC), diffusion capacity for carbon monoxide-single breath method (DLCO-SB), and diffusion capacity for carbon monoxide/alveolar volume (DLCO/VA)), and enhancing the exercise tolerance as well as quality of life of patients [19, 20]. Nonetheless, the multi-gene and multi-pathway mechanism of FXF in the treatment of pulmonary fibrosis remains largely unclear.

In the treatment of disease, the traditional Chinese medicine (TCM) herbs are less toxic. Nonetheless, the mechanism of TCM remains unknown so far, which may thus spark criticism as it becomes increasingly popular nowadays. Therefore, it is of great importance to investigate those target genes and pathways for the sake of modernizing TCM. Network pharmacology, a critical technical approach for investigating TCM formulas in systems biology, can visualize the herb and formula pharmacological effects and the underlying molecular mechanisms through multidisciplinary integration, including computer technology (CT), high-throughput omics, network database retrieval, and pharmacology [21, 22]. In the present work, network pharmacology analysis was adopted to predict pharmacodynamic material basis together with the underlying molecular mechanism of FXF for the treatment of pulmonary fibrosis, and the results were then validated by molecular docking study. Findings in this work lay theoretical basis to conduct experimental study and to apply FXF in clinic.

2. Materials and Methods

2.1. Mining the Potential Pharmacodynamic Compounds and Related Targets of FXF

The names of 7 herbs were imported into the BATMAN-TCM, Traditional Chinese Medicine Systems Pharmacology (TCMSP) database, and the analysis platform successively, to acquire the chemical compounds with corresponding information [23]. Thereafter, the drug likeness (DL) ≥ 0.18 and bioavailability (OB) ≥ 30% were used as the thresholds to screen those possible pharmacodynamic compounds in FXF by the method in literature [24, 25]. Additionally, certain compounds that did not conform to the selection thresholds but had great levels in each herb or had wide pharmacological actions, or were utilized to identify single medicine in the pharmacopeia, were incorporated as the possible pharmacodynamic compounds as well, including hirudin, salvianolic acid A, amygdalin, and oroxin [26–29]. Moreover, for those active components, the possible targets were discovered against the TCSMP database based on drug structural similarity assessment as well as reverse molecular docking, which were later used for constructing a potential FXF pharmacodynamic target set.

2.2. Pulmonary Fibrosis-Associated Target Screening

Targets associated with the pathogenesis or treatment of pulmonary fibrosis were acquired based on DisGeNet (https://www.disgenet.org/) [30], TTD (http://db.idrblab.net/ttd/) [31], CTD (http://ctdbase.org/) [32], and GeneCards (https://www.genecards.org/) [33] databases with the following keyword: “pulmonary fibrosis.” Targets were sorted out according to disease specificity index (DSI) from high to low against DisGeNet database, while targets larger than median were chosen. The abnormal drugs with related targets in TTD database were ruled out. For CTD database, those 200 most significant genes with the highest inference score were chosen. With regard to GeneCards, targets of Score ≥10 were chosen. Afterwards, targets retrieved based on these 4 databases were combined for the construction of a target set related to pulmonary fibrosis.

2.3. Establishment of the Core Component-Target Network of FXF for the Treatment of Pulmonary Fibrosis

Firstly, we chose the Uniprot database [34] in this work and chose species as “Homo sapiens”; then, target names acquired from these two steps were normalized for the acquisition of distinct Uniprot IDs and gene names. Later, both target sets were incorporated into Venny tool (https://bioinfogp.cnb.csic.es/tools/venny/), respectively, for obtaining overlapping targets, and they were regarded as the candidate FXF pharmacodynamic targets for the treatment of pulmonary fibrosis. Moreover, we applied Cytoscape in constructing the FXF candidate component-target network for the treatment of pulmonary fibrosis, where “component” was set as square, whereas “target” was set as circle. Thereafter, the node degree in candidate component-target network was calculated and ranked by CytoHubba plug-in [35], and the value (twice as high as the median degree value for the nodes) was adopted to select core components as well as core targets of FXF for treating pulmonary fibrosis, which were subsequently utilized for constructing the core network.

2.4. Compound-Target Interaction Validation

AutoDock Vina [36] (version: 1.1.2) was used to validate the relationships of compounds with targets. The compound mol2 structures were obtained based on TCMSP database, and target protein 3D structures were acquired from RCSB PDB database (http://www.rcsb.org/) [37]. Prior to molecular docking, proteins and ligands were prepared using AutoDockTools [38] (version: 1.5.6) and the MOE method. With regard to target proteins, their crystal structures were pre-treated, including removal of water molecules, 3D hydrogenation, protonation, correction of protein structure, energy optimization, and retention of target active region. Additionally, the ligand structures must conform to the low-energy conformation. Moreover, the box size and coordinates in molecular docking were finally determined based on ligand position. To achieve higher calculation accuracy, the exhaustiveness parameter was set at 20. The remaining parameters were the defaults except as otherwise noted. Thereafter, the components were integrated to the target proteins in a semiflexible manner, which produced altogether 9 conformations. The most affinal conformation was chosen to be the eventual docking conformation.

2.5. GO and KEGG Enrichment Analysis for Core Targets

For better exploring those biological processes (BPs), molecular functions (MFs), cellular components (CCs), and regulatory signal transduction pathways related to the FXF core targets for the treatment of pulmonary fibrosis, we input the information of selected core targets. Thereafter, GO-BP along with KEGG enrichment analysis was carried out on these targets using OmicShare online software [39] (https://www.omicshare.com/) at a threshold.

3. Results

3.1. Selection of Potential Pharmacodynamic Components and Targets of FXF

By searching against the TCMSP database, altogether 201 chemical components were obtained from 7 herbs in FXF and 126 were selected as the core components upon the DL ≥ 0.18 and OB ≥ 30% thresholds. Additionally, compounds with lower OB or DL value that had great contents in each herb were also collected as the potential pharmacodynamic components of FXF. At last, DS, TR, SZ, CS, TD, QG, and MHD had 61, 23, 1, 27, 9, 6, and 19 potential pharmacodynamic components, separately, of which, baicalin, beta-sitosterol, ellagic acid, stigmasterol, quercetin, and 2,5-dihydroxy-6,7-dimethoxyflavone were widely distributed in several herbs. Table S1 presents the information for the possible pharmacodynamic components of FXF.

Thereafter, the TCMSP database was searched to identify targets for the 130 potential pharmacodynamic components, and finally 286 were found (Table S2). DS, TR, SZ, CS, TD, QG, and MHD possessed 140, 71, 2, 114, 203, 207, and 131 potential targets, respectively. The 7 herbs obviously overlapped, regardless of the different target numbers in each herb, suggesting that those diverse herbs of FXF possibly exerted different actions by modulating similar targets.

To systemically and holistically understand the component-target network of FXF, Cytoscape was utilized to establish a network map, including 2969 edges and 416 nodes (Figure 1(a)). The node degree in the map represented the target or edge number related to nodes based on topological analysis (Figure 1(b)). There were 65 components discovered in the as-constructed network upon the threshold of median ≥14 degrees. Quercetin, beta-sitosterol, stigmasterol, baicalein, and luteolin of these components functioned in 326, 285, 132, 100, and 67 targets, separately, and they might serve as the obvious core components of FXF.

(a)

(b)

3.2. Selection of Core Targets in FXF for the Treatment of Pulmonary Fibrosis

Pulmonary fibrosis represents one of the polygenic genetic diseases, and its pathogenesis may be studied by examining the interactions between genes or between genes and the environment. In this study, the databases including CTD, TTD, DisGeNet, and GeneCards were retrieved, which discovered 318 targets related to pulmonary fibrosis (Table S3). Besides, 87 of these discovered targets were also recognized as the potential targets for the pharmacodynamic components in FXF, which suggested that FXF had some therapeutic effect (Figure 2(a) and Table S4). These 87 targets may be the candidate targets of FXF in treating pulmonary fibrosis, while the 88 pharmacodynamic active components that exert regulatory effect on these targets are the candidate components.

(a)

(b)

(c)

Additionally, the candidate component-target network of FXF for the treatment of pulmonary fibrosis was constructed using Cytoscape (Figure 2(b)). For better selecting the FXF core components and core targets for the treatment of pulmonary fibrosis, we adopted the CytoHubba plug-in for calculating and sorting those node topological parameters (degrees) in the as-constructed network (Table S5). Thereafter, nodes with degree value greater than or equal to twice as high as the median degree value (=5) for all the nodes were selected as the core components and core targets in FXF for the treatment of pulmonary fibrosis. Later, the core component-core target network (including 12 core components and 13 core targets) was established (Figure 2(c)). Of these core components, quercetin, beta-sitosterol, baicalein, luteolin, ellagic acid, wogonin, naringenin, tanshinone IIa, stigmasterol, and cryptotanshinone from DS, TR, CS, TD, QG, and MHD in FXF were the most significant in line with the degree values. Some of them were suggested previously to postpone pulmonary fibrosis course and improve the symptoms [40–43]. Besides, PTGS2, NOS2, CDK2, GSK3B, CCNA2, PPARG, MAPK14, CASP3, BCL2, and RELA ranked the top 10 places based on the degree value.

3.3. Compound-Target Network Validation

The relationships between components and targets were evaluated by molecular docking analysis, which helped to reduce network complexity while improving its accuracy. Thus, 10 core targets as well as 10 core compounds were discovered by molecular docking (Table 1 and Figure 3).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

Thereafter, PTGS2, NOS2, CDK2, GSK3B, CCNA2, PPARG, MAPK14, CASP3, BCL2, and RELA were searched in the PDB protein database, respectively, to acquire 3D structures. It was observed from the results of binding free energy in Table 1 that most of the core compounds could tightly combine with core targets. Typically, stigmasterol is most closely bounded to BCL2 (Figure 3(i)), CASP3 (Figure 3(h)), and RELA (Figure 3(j)). In addition, cryptotanshinone was the compound with the most tight correlation with CCNA2 (Figure 3(e)), GSK3B (Figure 3(d)), PPARG (Figure 3(f)), and PTGS2 (Figure 3(a)). Meanwhile, tanshinone IIa was the compound with the lowest binding energy to CDK2 (Figure 3(c)), MAPK14 (Figure 3(g)), and NOS2 (Figure 3(b)) in this molecular docking experiment. In virtual docking (Figure 3(a) and Figure 3(d)-3(f)), cryptotanshinone formed the hydrogen bonds with ARG₁₂₀, ARG₁₄₁, MET₂₁₀, and SER₃₄₂ of COX-2, GSK-3 beta, Cyclin-A2, and PPAR-gamma, respectively. Meanwhile, tanshinone IIa formed the potent hydrophobic interaction with iNOS (π-sigma bonds with TRP₁₉₄, Figure 3(b)), CDK2 (alkyl bonds with ALA₃₁, Figure 3(c)), and MAP kinase 14 (alkyl bonds with MET₇₈, Figure 3(g)). Typically, the stigmasterol-caspase-3 complex became stable via forming the π-sigma and π-alkyl with the PHE₂₅₂ and PHE₂₅₆ (Figure 3(h)). Correspondingly, the stigmasterol–Bcl-2 complex and stigmasterol-p65 complex got steady severally via forming hydrophobic interaction with MET₁₁₅ residues (Figure 3(i)) and the hydrogen bond with PHE₂₃₉ residues (Figure 3(j)).

3.4. Enrichment Analysis of FXF Core Targets for the Treatment of Pulmonary Fibrosis

For better understanding the multi-target and multi-pathway mechanism by which FXF treated pulmonary fibrosis, the OmicShare online approach was used to carry out GO as well as KEGG analysis on the thirteen core targets, to identify the biological processes/molecular functions along with signal transduction pathways in FXF for the treatment of pulmonary fibrosis (, FDR < 0.05). As shown, external stimulus (from oxidative stress, radiation, UV, chemical substances, and virus infection), apoptosis, cell cycle, aging, immune process, and protein metabolism were the mostly concentrated GO-BP terms (Figures 4(a)–4(c)). These core targets were mostly concentrated on several KEGG pathways, like IL-17, AGE-RAGE, TNF, HIF-1, PI3K-AKT, NOD-like receptor, T/B cell receptor, and virus infection-related pathways, indicating their important parts in the treatment of pulmonary fibrosis.

(a)

(b)

(c)

(d)

Figure 4

Enrichment analysis of core targets for FXF on the treatment of pulmonary fibrosis through OmicShare: (a) The top fifteen enriched GO-biological process; (b) GO-molecular functions; (c) GO-cellular components terms; and (d) KEGG pathways. The abscissa shows the enrichment factor, and the ordinate shows the GO terms or KEGG pathways. The color of the dot represents the adjusted p-value/q-value, and the size of the dot represents the number of core targets mapped to the reference GO terms or pathways.

4. Discussion

Pulmonary fibrosis, an interstitial lung disease featured by the fibrotic, chronic, and progressive nature, is an uncommon disorder generally seen among the old people [44]. Its incidence in the North American and European countries is reported to be 2.8–9.3/100,000 persons. At present, few epidemiological data are available in China, yet the incidence of pulmonary fibrosis shows an increasing trend recently [45]. Usual interstitial pneumonia (UIP) is characteristic of pulmonary fibrosis on high-resolution chest CT or in patients with a pulmonary fibrosis history. As for the major clinical manifestations, pulmonary fibrosis mainly manifests as progressive dyspnea, hypoxia, restricted ventilation impairment, gas exchange disturbance, and respiratory failure [46]. For the time being, there is no radical treatment for pulmonary fibrosis, and treatment is mainly conducted for delaying disease progression, improving patient quality of life, and prolonging patient survival. Pirfenidone and nintedanib are the anti-fibrotic drugs adopted to treat pulmonary fibrosis, which are effective to some extent; however, their high prices and adverse reactions have restricted their application. Recently, Chinese herbs have exerted vital parts in treating pulmonary fibrosis. It is reported in numerous articles that Chinese herbs alleviate clinical symptoms, improve patient quality of life, and boost the exercise tolerance while delaying lung function decline among pulmonary fibrosis patients [47, 48].

Chinese herbs represent an important TCM means to treat diseases, which has been utilized due to the complicated pulmonary fibrosis pathogenesis as well as the common patient syndrome [49]. Pulmonary fibrosis can be classified as Yin deficiency and internal heat syndrome, lung and kidney qi deficiency syndrome, phlegm turbid/heat obstruction of lung syndrome, and blood stasis syndrome, based on the diffuse interstitial lung disease diagnostic criteria in TCM (2012 Edition). Typically, clinical medication mainly aims to promote qi, accelerate blood circulation, eliminate blood stasis, and reinforce Yin. Therefore, FXF, which has the effect of nourishing Yin and removing blood stasis, is used as a cipher prescription in our hospital for the treatment of pulmonary fibrosis.

Herbal medicine has different and complicated components. Previous research on herbal medicine is still encountered with certain problems, like the complicated composition and the unknown mechanism. Network pharmacological analysis has been extensively adopted to investigate the pharmacological mechanisms of TCM currently. As a result, the holistic view in TCM was used in combination with network pharmacology analysis and the syndrome differentiation science system in this study to identify core pharmacodynamic material basis and targets. Then, molecular docking study was carried out on those core active ingredients together with the core targets using AutoDock Vina for verifying results of network pharmacology analysis. Finally, these core targets were subjected to GO and pathway enrichment analysis.

As revealed in this work, quercetin, beta-sitosterol, baicalein, luteolin, ellagic acid, wogonin, naringenin, tanshinone IIa, stigmasterol, and cryptotanshinone were identified as the core compounds for treating pulmonary fibrosis. Some studies demonstrate the effect of quercetin, baicalein, tanshinone IIA, and cryptotanshinone on suppressing experimental pulmonary fibrosis induced by bleomycin or silica through promoting lung tissue self-healing, modulating disturbed redox-balance in lung tissue, alleviating local inflammation, reducing collagen deposition, promoting ECM degradation, inducing pulmonary fibrosis cell apoptosis, reversing angiotensin production, and regulating multiple signaling molecules and pathways (such as TGF-beta, sphingosine kinase, SMAD, Nrf2, STAT3) [41, 42, 50–54]. Our research results were consistent with those in previous study, indicating that these compounds were the core effective components of FXF to treat pulmonary fibrosis.

According to the node degrees in our constructed core component-target network, PTGS2, NOS2, CDK2, GSK3B, CCNA2, PPARG, MAPK14, CASP3, BCL2, and RELA were suggested as core targets of FXF for the treatment of pulmonary fibrosis. In addition, molecular docking was conducted to validate the associations of core components with core targets. Upon GO-BP annotation analysis, the core targets exerted vital parts in external stimulus (from oxidative stress, radiation, UV, chemical substances, and virus infection), apoptosis, cell cycle, aging, immune process, and protein metabolism. As revealed by KEGG pathway analysis, numerous pathways were tightly associated with the pulmonary fibrosis pathogenic mechanism. Typically, IL-17, AGE-RAGE, TNF, HIF-1, PI3K-AKT, NOD-like receptor, T/B cell receptor, and virus (EBV and hepatitis virus) infection-related pathways were the most significant pathways involved. The associations of diverse pathways, including EBV and hepatitis virus infections, with pulmonary fibrosis are identified through fundamental clinical studies to induce pulmonary fibrosis [55–57]. Meanwhile, epithelial-mesenchymal transition (EMT) in pulmonary fibrosis may also be related to activating the NOD-like receptor protein (NLRP) 1 and NLRP3 inflammatory pathways [58]. The advanced glycation end products (AGE) may be obtained from non-enzymatic reaction of proteins and lipids with various oxidants in the process of aging. In addition, the receptor for AGE (RAGE) has been demonstrated to be associated with alveolar homeostasis as well as pulmonary fibrosis. The aberrant epithelial stromal repair capacity in pulmonary fibrosis is linked with this aging process. The elevated AGE-RAGE proportion in pulmonary fibrosis is correlated with epithelial stromal repair during pulmonary fibrosis [59, 60]. In the pathological process of pulmonary fibrosis, the infiltration of various immune cells and the release of immuno-inflammatory factors play an important role, which is also considered as a new therapeutic target to reverse the disease, such as IL-17 [61] and T cell receptor [62]. In addition, the mitochondria-mediated apoptosis also exerts a vital part during pulmonary fibrosis genesis and progression mediated by both mitochondria and ER stress (triggered by HIF-1alpha) [63–65].

Nonetheless, certain limitations should be noted in the present work [66], like the data-based network pharmacology, which made it impossible to examine whether all data were comprehensively collected and whether the criteria to select core components were completely accurate. Besides, the efficacy of additional components and effects of dosage should be further examined.

5. Conclusion

To sum up, the mechanism of action by which FXF treats pulmonary fibrosis may be related to the regulation of several pathways, like IL-17, AGE-RAGE, TNF, HIF-1, PI3K-AKT, NOD-like receptor, T/B cell receptor, and virus infection-related pathways.

Data Availability

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Xiao-Li Chen and Cheng Tang contributed equally to this study.

Acknowledgments

This study was financially supported by the National Natural Science Foundation of China (81704040), Foundation for Jiangsu Administration of Traditional Chinese Medicine (YB2017029), and Foundation for Beijing Medical and Health (F2031E). The authors thank the members of their laboratory and their collaborators.

Supplementary Materials

Table S1: all the pharmacodynamic ingredients of FXF. Table S2: all the potential pharmacodynamic targets of FXF. Table S3: known pulmonary fibrosis-related targets. Table S4: FXF shared 87 potential pharmacodynamic targets with known pulmonary fibrosis-related targets. Table S5: degree values of nodes in the candidate active ingredient-target network of FXF in treating pulmonary fibrosis. (Supplementary Materials)

References

A. K. Agarwal and N. Huda, Interstitial Pulmonary Fibrosis, StatPearls, Treasure Island, FL, USA, 2020.
L. Richeldi, A. S. Rubin, S. Avdeev, Z. F. Udwadia, and Z. J. Xu, “Idiopathic pulmonary fibrosis in BRIC countries: the cases of Brazil, Russia, India, and China,” BMC Medicine, vol. 13, no. 1, 237 pages, 2015.
View at: Publisher Site | Google Scholar
C. Shimbori, J. Gauldie, and M. Kolb, “Extracellular matrix microenvironment contributes actively to pulmonary fibrosis,” Current Opinion in Pulmonary Medicine, vol. 19, no. 5, pp. 446–452, 2013.
View at: Publisher Site | Google Scholar
M. A. T. Freeberg, A. Perelas, J. K. Rebman, R. P. Phipps, T. H. Thatcher, and P. J. Sime, “Mechanical feed-forward loops contribute to idiopathic pulmonary fibrosis,” American Journal of Pathology, vol. 191, pp. 18–25, 2020.
View at: Google Scholar
L. A. Murtha, M. J. Schuliga, N. S. Mabotuwana et al., “The processes and mechanisms of cardiac and pulmonary fibrosis,” Frontiers in Physiology, vol. 8, 777 pages, 2017.
View at: Publisher Site | Google Scholar
R. J. Hewitt and T. M. Maher, “Idiopathic pulmonary fibrosis: new and emerging treatment options,” Drugs & Aging, vol. 36, no. 6, pp. 485–492, 2019.
View at: Publisher Site | Google Scholar
Y. Uchida, S. Ikeda, A. Sekine et al., “Tolerability and safety of nintedanib in elderly patients with idiopathic pulmonary fibrosis,” Respiratory Investigation, vol. 59, pp. 99–105, 2020.
View at: Google Scholar
L. Richeldi, M. Kolb, S. Jouneau et al., “Efficacy and safety of nintedanib in patients with advanced idiopathic pulmonary fibrosis,” BMC Pulmonary Medicine, vol. 20, no. 1, 3 pages, 2020.
View at: Publisher Site | Google Scholar
B. Marcos Ribes, J. N. Sancho-Chust, A. Talens et al., “Effectiveness and safety of pirfenidone for idiopathic pulmonary fibrosis,” European Journal of Hospital Pharmacy, vol. 27, no. 6, pp. 350–354, 2020.
View at: Publisher Site | Google Scholar
M. K. Glassberg, S. D. Nathan, C.-Y. Lin et al., “Cardiovascular risks, bleeding risks, and clinical events from 3 phase III trials of pirfenidone in patients with idiopathic pulmonary fibrosis,” Advances in Therapy, vol. 36, no. 10, pp. 2910–2926, 2019.
View at: Publisher Site | Google Scholar
L. H. Lancaster, J. A. de Andrade, J. D. Zibrak et al., “Pirfenidone safety and adverse event management in idiopathic pulmonary fibrosis,” European Respiratory Review: An Official Journal of the European Respiratory Society, vol. 26, 2017.
View at: Publisher Site | Google Scholar
M. L. Wencel, T. Haselkorn, S. L. Limb et al., “Real-world practice patterns for prevention and management of potential adverse events with pirfenidone in patients with idiopathic pulmonary fibrosis,” Pulmonary Therapy, vol. 4, no. 1, pp. 103–114, 2018.
View at: Publisher Site | Google Scholar
J. A. Galli, A. Pandya, M. Vega-Olivo, C. Dass, H. Zhao, and G. J. Criner, “Pirfenidone and nintedanib for pulmonary fibrosis in clinical practice: tolerability and adverse drug reactions,” Respirology, vol. 22, no. 6, pp. 1171–1178, 2017.
View at: Publisher Site | Google Scholar
J. Guiot, I. Struman, V. Chavez et al., “Altered epigenetic features in circulating nucleosomes in idiopathic pulmonary fibrosis,” Clinical Epigenetics, vol. 9, no. 1, 84 pages, 2017.
View at: Publisher Site | Google Scholar
A. Carleo, C. Landi, A. Prasse et al., “Proteomic characterization of idiopathic pulmonary fibrosis patients: stable versus acute exacerbation,” Monaldi Archives for Chest Disease, vol. 90, 2020.
View at: Publisher Site | Google Scholar
Y. Tian, H. Li, Y. Gao et al., “Quantitative proteomic characterization of lung tissue in idiopathic pulmonary fibrosis,” Clinical Proteomics, vol. 16, no. 1, 6 pages, 2019.
View at: Publisher Site | Google Scholar
Z. Xu, L. Mo, X. Feng, M. Huang, and L. Li, “Using bioinformatics approach identifies key genes and pathways in idiopathic pulmonary fibrosis,” Medicine, vol. 99, no. 36, Article ID e22099, 2020.
View at: Publisher Site | Google Scholar
Y. Feng, D. Zhou, Q. Xiao et al., “Treatment of pulmonary fibrosis based on Yin deficiency and blood stasis,” Acta Chinese Medicine and Pharmacology, vol. 3, pp. 64–66, 2020, in Chinese.
View at: Google Scholar
Q. Xiao, X. Chen, H. Sun et al., “An analysis of professor Zhu Qiyong’s experience in treating pulmonary fibrosis,” Chinese Journal of Ethnomedicine and Ethnopharmacy, vol. 28, pp. 61–63, 2019, in Chinese.
View at: Google Scholar
Y. Liu, Q. Zhu, X. Chen, Q. Xiao, and J. Ju, “Optimization of extraction technology of Fei-Xian formula,” Chinese Herbal Medicine, vol. 42, pp. 1581–1584, 2020, in Chinese.
View at: Google Scholar
Z. Zhou, B. Chen, S. Chen et al., “Applications of network pharmacology in traditional Chinese medicine research,” Evidence-Based Complementary and Alternative Medicine, vol. 2020, Article ID 1646905, 7 pages, 2020.
View at: Publisher Site | Google Scholar
C. Liang, W. Niu, F. Wu et al., “Network pharmacology for the identification of phytochemicals in traditional Chinese medicine for COVID-19 that may regulate interleukin-6,” Bioscience Reports, vol. 41, 2020.
View at: Publisher Site | Google Scholar
R. Zhang, X. Zhu, H. Bai, and K. Ning, “Network pharmacology databases for traditional Chinese medicine: review and assessment,” Frontiers in Pharmacology, vol. 10, 123 pages, 2019.
View at: Publisher Site | Google Scholar
J. Li, M. Ye, J. Gao, Y. Zhang, Q. Zhu, and W. Liang, “Systematic understanding of mechanism of Yi-Qi-Huo-Xue decoction against intracerebral hemorrhagic stroke using a network pharmacology approach,” Medical Science Monitor: International Medical Journal of Experimental and Clinical Research, vol. 26, Article ID e921849, 2020.
View at: Publisher Site | Google Scholar
S. Wang, H. Wang, and Y. Lu, “Tianfoshen oral liquid: a CFDA approved clinical traditional Chinese medicine, normalizes major cellular pathways disordered during colorectal carcinogenesis,” Oncotarget, vol. 8, no. 9, pp. 14549–14569, 2017.
View at: Publisher Site | Google Scholar
Y. Pan, H. Fu, Q. Kong et al., “Prevention of pulmonary fibrosis with salvianolic acid a by inducing fibroblast cell cycle arrest and promoting apoptosis,” Journal of Ethnopharmacology, vol. 155, no. 3, pp. 1589–1596, 2014.
View at: Publisher Site | Google Scholar
H. K. Du, F. C. Song, X. Zhou, H. Li, and J. P. Zhang, “Effect of amygdalin on serum proteinic biomarker in pulmonary fibrosis of bleomycin-induced rat,” Zhonghua Lao Dong Wei Sheng Zhi Ye Bing Za Zhi, vol. 28, pp. 260–263, 2010.
View at: Google Scholar
N. Li, W. Men, Y. Zheng, H. Wang, and X. Meng, “Oroxin B induces apoptosis by down-regulating microRNA-221 resulting in the inactivation of the PTEN/PI3K/AKT pathway in liver cancer,” Molecules, vol. 24, 4384 pages, 2019.
View at: Publisher Site | Google Scholar
L. Shen, S. Lei, L. Huang et al., “Therapeutic effects of the rhSOD2-Hirudin fusion protein on bleomycin-induced pulmonary fibrosis in mice,” European Journal of Pharmacology, vol. 852, pp. 77–89, 2019.
View at: Publisher Site | Google Scholar
J. Piñero, À. Bravo, N. Queralt-Rosinach et al., “DisGeNET: a comprehensive platform integrating information on human disease-associated genes and variants,” Nucleic Acids Research, vol. 45, no. D1, pp. D833–D839, 2017.
View at: Publisher Site | Google Scholar
F. Zhu, B. Han, P. Kumar et al., “Update of TTD: therapeutic target database,” Nucleic Acids Research, vol. 38, no. 1, pp. D787–D791, 2010.
View at: Publisher Site | Google Scholar
C. J. Mattingly, M. C. Rosenstein, G. T. Colby, J. N. Forrest Jr., and J. L. Boyer, “The Comparative toxicogenomics database (CTD): a resource for comparative toxicological studies,” Journal of Experimental Zoology Part A: Comparative Experimental Biology, vol. 305A, no. 9, pp. 689–692, 2006.
View at: Publisher Site | Google Scholar
N. Rappaport, S. Fishilevich, R. Nudel et al., “Rational confederation of genes and diseases: NGS interpretation via GeneCards, MalaCards and VarElect,” BioMedical Engineering OnLine, vol. 16, no. S1, 72 pages, 2017.
View at: Publisher Site | Google Scholar
S. Pundir, M. J. Martin, and C. O’Donovan, “UniProt protein knowledgebase,” Protein Bioinformatics, vol. 1558, pp. 41–55, 2017.
View at: Publisher Site | Google Scholar
C. H. Chin, S. H. Chen, H. H. Wu, C. W. Ho, M. T. Ko, and C. Y. Lin, “cytoHubba: identifying hub objects and sub-networks from complex interactome,” BMC Systems Biology, vol. 8, no. 4, S11 pages, 2014.
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, pp. 455–461, 2010.
View at: Google Scholar
S. K. Burley, H. M. Berman, C. Bhikadiya et al., “RCSB protein data bank: biological macromolecular structures enabling research and education in fundamental biology, biomedicine, biotechnology and energy,” Nucleic Acids Research, vol. 47, no. D1, pp. D464–D474, 2019.
View at: Publisher Site | Google Scholar
N. El-Hachem, B. Haibe-Kains, A. Khalil, F. H. Kobeissy, and G. Nemer, “AutoDock and AutoDockTools for protein-ligand docking: beta-site amyloid precursor protein cleaving enzyme 1 (BACE1) as a case study,” Methods in Molecular Biology, vol. 1598, pp. 391–403, 2017.
View at: Publisher Site | Google Scholar
M. Ni, X. Liu, J. Wu et al., “Identification of candidate biomarkers correlated with the pathogenesis and prognosis of non-small cell lung cancer via integrated bioinformatics analysis,” Frontiers in Genetics, vol. 9, 469 pages, 2018.
View at: Publisher Site | Google Scholar
S. Mehrzadi, P. Hosseini, M. Mehrabani et al., “Attenuation of bleomycin-induced pulmonary fibrosis in Wistar rats by combination treatment of two natural phenolic compounds: quercetin and gallic acid,” Nutrition and Cancer, pp. 1–11, 2020.
View at: Publisher Site | Google Scholar
Y. Gao, J. Lu, Y. Zhang, Y. Chen, Z. Gu, and X. Jiang, “Baicalein attenuates bleomycin-induced pulmonary fibrosis in rats through inhibition of miR-21,” Pulmonary Pharmacology & Therapeutics, vol. 26, no. 6, pp. 649–654, 2013.
View at: Publisher Site | Google Scholar
F. Feng, N. Li, P. Cheng et al., “Tanshinone IIA attenuates silica-induced pulmonary fibrosis via inhibition of TGF-β1-Smad signaling pathway,” Biomedicine & Pharmacotherapy, vol. 121, Article ID 109586, 2020.
View at: Publisher Site | Google Scholar
Q. Zhang, C. Gan, H. Liu et al., “Cryptotanshinone reverses the epithelial‐mesenchymal transformation process and attenuates bleomycin‐induced pulmonary fibrosis,” Phytotherapy Research, vol. 34, no. 10, pp. 2685–2696, 2020.
View at: Publisher Site | Google Scholar
M.-S. Cao, J. Sheng, T.-Z. Wang et al., “Acute exacerbation of idiopathic pulmonary fibrosis: usual interstitial pneumonitis vs. possible usual interstitial pneumonitis pattern,” Chinese Medical Journal, vol. 132, no. 18, pp. 2177–2184, 2019.
View at: Publisher Site | Google Scholar
J. Hutchinson, A. Fogarty, R. Hubbard, and T. McKeever, “Global incidence and mortality of idiopathic pulmonary fibrosis: a systematic review,” European Respiratory Journal, vol. 46, no. 3, pp. 795–806, 2015.
View at: Publisher Site | Google Scholar
G. Raghu, M. Remy-Jardin, J. L. Myers et al., “Diagnosis of idiopathic pulmonary fibrosis. An official ATS/ERS/JRS/ALAT clinical practice guideline,” American Journal of Respiratory and Critical Care Medicine, vol. 198, no. 5, pp. e44–e68, 2018.
View at: Publisher Site | Google Scholar
H.-Y. Zhang, L.-J. Pang, X.-D. Lv, C. Liu, and M.-H. Nan, “Multiple traditional Chinese medicine interventions for idiopathic pulmonary fibrosis: a protocol for systematic review and meta-analysis of overview,” Medicine, vol. 99, no. 39, Article ID e22396, 2020.
View at: Publisher Site | Google Scholar
K. Ji, J. Ma, L. Wang, N. Li, S. Dong, and L. Shi, “Efficacy and safety of traditional Chinese medicine in idiopathic pulmonary fibrosis: a meta-analysis,” Evidence-Based Complementary and Alternative Medicine, vol. 2020, Article ID 1752387, 11 pages, 2020.
View at: Publisher Site | Google Scholar
X.-q. Yu, S.-g. Yang, Y. Xie, and J.-s. Li, “Traditional Chinese medicine in the treatment of idiopathic pulmonary fibrosis based on syndrome differentiation: study protocol of an exploratory trial,” Journal of Integrative Medicine, vol. 18, no. 2, pp. 163–168, 2020.
View at: Publisher Site | Google Scholar
C. Veith, M. Drent, A. Bast, F. J. van Schooten, and A. W. Boots, “The disturbed redox-balance in pulmonary fibrosis is modulated by the plant flavonoid quercetin,” Toxicology and Applied Pharmacology, vol. 336, pp. 40–48, 2017.
View at: Publisher Site | Google Scholar
X. Zhang, Y. Cai, W. Zhang, and X. Chen, “Quercetin ameliorates pulmonary fibrosis by inhibiting SphK1/S1P signaling,” Biochemistry and Cell Biology, vol. 96, no. 6, pp. 742–751, 2018.
View at: Publisher Site | Google Scholar
H. Wu, Y. Li, Y. Wang et al., “Tanshinone IIA attenuates bleomycin-induced pulmonary fibrosis via modulating angiotensin-converting enzyme 2/angiotensin-(1-7) axis in rats,” International Journal of Medical Sciences, vol. 11, no. 6, pp. 578–586, 2014.
View at: Publisher Site | Google Scholar
Y. Zhang, W. Lu, X. Zhang et al., “Cryptotanshinone protects against pulmonary fibrosis through inhibiting Smad and STAT3 signaling pathways,” Pharmacological Research, vol. 147, Article ID 104307, 2019.
View at: Publisher Site | Google Scholar
Z. Zhu, Q. Li, C. Xu et al., “Sodium tanshinone IIA sulfonate attenuates silica-induced pulmonary fibrosis in rats via activation of the Nrf2 and thioredoxin system,” Environmental Toxicology and Pharmacology, vol. 80, Article ID 103461, 2020.
View at: Publisher Site | Google Scholar
K. Manika, S. Alexiou-Daniel, D. Papakosta et al., “Epstein-Barr virus DNA in bronchoalveolar lavage fluid from patients with idiopathic pulmonary fibrosis,” Sarcoidosis, Vasculitis, and Diffuse Lung Diseases: Official Journal of WASOG, vol. 24, pp. 134–140, 2007.
View at: Google Scholar
A. H. Jafarian, N. Mohammadian Roshan, H. Ayatollahi, A. A. Omidi, M. Ghaznavi, and M. Gharib, “Epstein-Barr virus and human herpesvirus 8 in idiopathic pulmonary fibrosis,” Iranian Journal of Pathology, vol. 15, no. 1, pp. 30–33, 2020.
View at: Publisher Site | Google Scholar
O. Okutan and Ö. Ayten, “Hepatitis C and pulmonary fibrosis,” Tuberkuloz ve Toraks, vol. 65, no. 2, pp. 131–137, 2017.
View at: Publisher Site | Google Scholar
R. Tian, Y. Zhu, J. Yao et al., “NLRP3 participates in the regulation of EMT in bleomycin-induced pulmonary fibrosis,” Experimental Cell Research, vol. 357, no. 2, pp. 328–334, 2017.
View at: Publisher Site | Google Scholar
S. Y. Kyung, K. H. Byun, J. Y. Yoon et al., “Advanced glycation end-products and receptor for advanced glycation end-products expression in patients with idiopathic pulmonary fibrosis and NSIP,” International Journal of Clinical and Experimental Pathology, vol. 7, pp. 221–228, 2014.
View at: Google Scholar
K. Yamaguchi, H. Iwamoto, W. Mazur et al., “Reduced endogenous secretory RAGE in blood and bronchoalveolar lavage fluid is associated with poor prognosis in idiopathic pulmonary fibrosis,” Respiratory Research, vol. 21, no. 1, 145 pages, 2020.
View at: Publisher Site | Google Scholar
W. Wang, Z. Liu, J. Niu et al., “Feibi recipe reduced pulmonary fibrosis induced by bleomycin in mice by regulating BRP39/IL-17 and TGFβ1/Smad3 signal pathways,” Evidence-Based Complementary and Alternative Medicine, vol. 2020, Article ID 5814658, 13 pages, 2020.
View at: Publisher Site | Google Scholar
S. Segawa, D. Goto, A. Iizuka et al., “The regulatory role of interferon-γ producing gamma delta T cells via the suppression of T helper 17 cell activity in bleomycin-induced pulmonary fibrosis,” Clinical & Experimental Immunology, vol. 185, no. 3, pp. 348–360, 2016.
View at: Publisher Site | Google Scholar
E. Delbrel, A. Soumare, A. Naguez et al., “HIF-1α triggers ER stress and CHOP-mediated apoptosis in alveolar epithelial cells, a key event in pulmonary fibrosis,” Scientific Reports, vol. 8, no. 1, Article ID 17939, 2018.
View at: Publisher Site | Google Scholar
Z. Shi, L. Xu, H. Xie et al., “Attenuation of intermittent hypoxia-induced apoptosis and fibrosis in pulmonary tissues via suppression of ER stress activation,” BMC Pulmonary Medicine, vol. 20, no. 1, 92 pages, 2020.
View at: Publisher Site | Google Scholar
C. L. Fattman, “Apoptosis in pulmonary fibrosis: too much or not enough?” Antioxidants & Redox Signaling, vol. 10, no. 2, pp. 379–386, 2008.
View at: Publisher Site | Google Scholar
L. Wu, Y. Chen, J. Yi, Y. Zhuang, L. Cui, and C. Ye, “Mechanism of action of Bu-Fei-Yi-Shen formula in treating chronic obstructive pulmonary disease based on network pharmacology analysis and molecular docking validation,” BioMed Research International, vol. 2020, Article ID 9105972, 12 pages, 2020.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Xiao-Li Chen 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

659

Downloads

846

Citations