Abstract

Background. Idiopathic pulmonary fibrosis (IPF) is a progressive disease with high mortality and poor prognosis. The prognostic signatures related to conventional therapy response remain limited. The Wenfei Buqi Tongluo (WBT) formula, a traditional Chinese medicine (TCM) formula, has been widely utilized to treat respiratory diseases in China, which is particularly effective in promoting inflammatory absorption. In this study, we aim to explore the mechanism of the WBT formula in the inhibition of inflammatory response during IPF, based on network pharmacology and in vivo experiments. Methods. Network pharmacology was applied to predict the changes of biological processes and potential pathways for the WBT formula against IPF. Histopathological changes, inflammatory factors (IL-6, IL-1β, and TNF-α), and the proteins of the TLR4/MyD88/NF-κB pathway in bleomycin- (BLM-) induced mice model were examined by hematoxylin-eosin (H&E), Masson or immunohistochemistry staining, Western blot, and enzyme-linked immunosorbent assay analysis. Results. A total of 163 possible components and 167 potential targets between the WBT formula and IPF were obtained. The enrichments of network pharmacology showed that inflammation response, TNF, and NF-κB pathways were involved in the treatment of WBT against IPF. The in vivo experiments indicated that the WBT formula could ameliorate inflammatory exudation and collagen deposition at a histopathology level in the BLM-induced mice model. The levels of IL-6, IL-1β, and TNF-α were reduced after the WBT formula treatment. Moreover, the expressions of phosphorylated-NF-κB p65, TLR4, and MyD88 were significantly downregulated by the WBT formula, compared with the BLM-induced group. Conclusion. These results indicated that the WBT formula can suppress BLM-induced IPF in a mouse model by inhibiting the inflammation via the TLR4/MyD88/NF-κB pathway. This study provides a new insight into the molecular mechanisms of the WBT formula in the application at the clinic.

1. Introduction

Idiopathic pulmonary fibrosis (IPF) is a chronic, progressive disease with a median survival time of three to five years since diagnosis [1]. IPF, the most common interstitial lung disease, is characterized by interstitial inflammation, fibrocyte proliferation of the alveolar wall, and fibrosis [2, 3]. Inflammatory and oxidative injuries, shortened telomeres, epithelial-mesenchymal transition (EMT), and endoplasmic reticulum stress lead to increased secretion of fibrotic factors [4]. Repeated inflammation and lung tissue injury appear to be an early-stage phenotype of pulmonary fibrosis (PF). The aberrant repair of inflammatory cells subsequently leads to the cross talk among epithelial cells, the extracellular matrix, and nearby mesenchymal cells [1]. Chronic inflammation triggered the secretion of transforming growth factor-β (TGF-β) in the alveolar compartments, which was the primary mechanism driving fibroblast activation and proliferation [5]. Inflammatory response, wounding, and oxidative stress play important roles in the prior phase of IPF [6]. It has been proven that mice deficient for the toll-like receptor 4 (TLR4) developed worse fibrosis phenotype and downregulation of cell surface hyaluronan, which has been evidenced in patients with IPF [7]. TLR4 activates the canonical NF-κB pathway through MyD88 as an innate immune response, which is related to IPF [8, 9]. Importantly, the NF-κB pathway as a key mediator for inflammatory response modulates the secretion of numerous cytokines and plays a key role in the inflammatory stage of PF [10, 11]. Therefore, the suppression of inflammation through the TLR4/MyD88/NF-κB pathway in an early stage of the fibrosis process can effectively improve lung fibrosis.

The Wenfei Buqi Tongluo (WBT) formula is an effective prescription for treating IPF, based on a traditional and classical Chinese medicine formula, Buyang Huanwu decoction, including Astragalus membranaceus (Fisch.) Bunge (Huang qi), Angelica sinensis (Oliv.) Diels (Dang gui), Prunus persica (L.) Batsch (Tao ren), Pheretima aspergillum (E. Perrier) (Di long), Ligusticum striatum DC. (Chuan xiong), Carthamus tinctorius L. (Hong hua), and Radix Paeoniae Rubra (Chi shao), which was created by Qing-Ren Wang for hundreds of years during the Qing Dynasty. Buyang Huanwu decoction ameliorates PF through inhibiting the PI3K/Akt signaling pathway [12, 13]. The WBT formula is composed of 13 Chinese medicines, containing the five medicines in Buyang Huanwu decoction and other eight herbs, such as Scutellaria baicalensis Georgi (Huang qin), Salvia miltiorrhiza Bunge (Dan shen), Polygonum cuspidatum Siebold & Zucc (Hu zhang), Aster tataricus L.f. (Zi wan), Tussilago farfara L. (Kuan donghua), Pinellia ternata (Thunb.) Makino (Ban xia), Clematis chinensis Osbeck (Wei lingxian), and Siegesbeckia orientalis L. (Xi xiancao) (Table 1). As we previously reported, the WBT formula inhibited cell proliferation, morphology, and EMT in the TGF-β1-induced A549 cell model [14]. Multiple components in the WBT formula have been proven that they are effective in the prevention and treatment of main pathological progresses of IPF. Isorhamnetin and kaempferol are bioactive compounds from Tussilago farfara L. and Pinellia ternata (Thunb.) Makino, which reduces inflammatory cytokines and oxidative stress to protect acute lung injury in mice model by regulating the NF-κB pathway [1517]. A bioactive compound from Ligusticum striatum DC., ligustilide, can suppress oxidative stress and inflammation to avoid BLM-induced PF in rat model through the inhibition of the TLR4/MyD88/NF-κB pathway [18]. However, the protective effects and possible mechanisms of the WBT formula in the early phase of IPF, especially the inflammatory process, have not been thoroughly explored. In this study, we first screened out IPF-related targets and potential chemical components of the WBT formula by network pharmacology, according to several common databases. Based on the prediction data, the pathway enrichment analysis was performed to elucidate the possible mechanisms of the WBT formula against IPF. Furthermore, the mice model after BLM induction for 7 days was employed to validate the anti-inflammatory effects of the WBT formula and investigate its regulatory roles on the TLR4/MyD88/NF-κB pathway. Our study might provide new insights into the molecular mechanism of the WBT formula for inhibiting inflammatory response during IPF.

2. Materials and Methods

2.1. Reagents

BLM and pirfenidone (PFD) were purchased from MedChemExpress (Monmouth Junction, NJ, USA). Antibodies against TLR4 (95/120 kDa, sc-293072), MyD88 (33 kDa, sc-74532), p-NF-κB p65 (65 kDa, sc-166748), and fibronectin (FN, 220 kDa, sc-8422) were purchased from Santa Cruz Biotechnology (Rosemont, IL, USA). IκBα (39 kDa, #4814) and NF-κB p65 (65 kDa, #8242) were obtained from Cell Signaling Technology (Beverly, MA, USA). Tubulin (50 kDa, 11224-1-AP) was purchased from ProteinTech (Rosemont, IL, USA).

2.2. Network Pharmacology Analysis for the WBT Formula and IPF

The potential gene targets of the 13 Chinese medicines in the WBT formula were identified via network pharmacology analysis. The monomer chemicals in the WBT formula were retrieved from Traditional Chinese Medicine databases (TCMSP, http://tcmspw.com/tcmsp.php) [19] and Traditional Chinese Medicine Information Database (TCMID, http://www.megabionet.org/tcmid/) [20]. In the TCMSP database, oral bioavailability ≥30% and drug-likeness ≥0.18 were conditioned to the absorption parameters for the compound screening. To collect the potential gene targets of these components of the WBT formula, the following two databases, TCMSP and SwissTargetPrediction Database (STPD, http://old.swisstargetprediction.ch/), were used [21]. Then, the Universal Protein Resource Knowledgebase (UniProt Knowledgebase, http://www.uniprot.org/, entering at Jan 2021) was applied to unify the official gene symbol, which was a collection from the candidates above. The IPF-associated target genes were acquired from GeneCards (https://www.genecards.org/) [22], the screening parameter was “.” Protein-protein interaction (PPI) data were obtained from STRING (https://cn.string-db.org/) with parameter conditions filtered by “Homo sapiens” () and visualized using Cytoscape. Finally, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched by using the Metascape (http://metascape.org/) tool.

2.3. Preparation of the WBT Formula

The 13 Chinese medicines present in the WBT formula were purchased from the Hongjian Pharmacy (Changchun, China), which were deposited in the Jilin Ginseng Academy, Changchun University of Chinese Medicine (Changchun, China). Raw materials of the WBT formula were codecocted in the distilled water for 1 h twice using 10 times of the total weight of the medicine mixture to obtain the aqueous extract. After filtering and centrifuging at 3,500 rpm for 15 min, the supernatants of WBT were subjected to vacuum to obtain the powder for further experiments.

2.4. Establishment of BLM-Induced Pulmonary Inflammation and Fibrosis and Drug Treatment

All mice (the Vital River Laboratory Animal Technology Co. Ltd, Beijing, China) were housed under standard conditions at the Experimental Animal Center, Changchun University of Chinese Medicine (Changchun, China). All experiments were approved by the Experimental Animal Administration Committee of Changchun University of Chinese Medicine (Approval No. 2021230). After three days of acclimatization, a total of 60 male C57BL/6N mice with body weight about , aged 6–8 weeks, were randomly divided into six groups (): sham, BLM (3 mg/kg), WBT (3, 6, and 12 g/kg), and PFD (200 mg/kg) groups. The pulmonary inflammation and fibrosis model was established by intratracheal injection with BLM (3 mg/kg) under anesthesia by 0.3% pentobarbital sodium on day 0. The sham group was instilled with an equal volume of normal saline only. After BLM administration, mice were intragastrically administered with different doses of the WBT formula or PFD once a day for 7 consecutive days. On day 8, all the mice were sacrificed to collect serum, bronchoalveolar lavage fluid (BALF), and lung tissues for further analysis.

2.5. Interleukin-6 (IL-6), TNF-α, and Interleukin-1β (IL-1β) Measurement

According to the instructions, the IL-6 level in the serum and TNF-α and IL-1β in BALF from the mice of different groups were determined by the enzyme-linked immunosorbent assay kit (Sinobest Bio, Shanghai, China).

2.6. Histopathological Analysis

The mouse lung tissue was fixed using 4% formaldehyde, was embedded in paraffin, and then was cut into the slices at 4 μm thickness. The lung sections were deparaffinized with xylene after being rehydrated in water by graded concentrations of ethanol solution. The sections were stained to estimate lung inflammation and fibrotic changes using H&E and Masson’s trichrome stainings. For immunohistochemical (IHC) staining, the slides were submerged with a citric acid solution (pH 6.0, 0.01 mol/L; Servicebio, Wuhan, China) for 10 min for antigen retrieval. After incubation with 3% hydrogen peroxide for 25 min, the slides were incubated with primary antibodies against TLR4, MyD88, and p-NF-κB overnight at 4°C and secondary antibody for 1 h. The positive staining was determined with a 3,3-diaminobenzidine substrate and counterstained with hematoxylin. Histological images were acquired using an M8 microscope (PreciPoint, Thüringen, Germany).

2.7. Western Blot Analysis

Proteins from mouse lung tissues were extracted, separated on 12% SDS-PAGE gels, and transferred to the PVDF membranes (Roche, Basel, Switzerland). After blocking with 5% BSA or nonfat milk for 1.5 h, the membranes were incubated with primary antibodies against TLR4, MyD88, IκBα, NF-κB p65, p-NF-κB p65, and Tubulin overnight at 4°C. After the incubation with HRP-conjugated secondary antibody (ProteinTech, IL, USA) for 1 h at room temperature, the membranes were washed three times with Tris-buffered saline with 0.1% Tween 20. The protein brands were visualized and detected using the ECL Moon kit (Beyotime Biotechnology, Shanghai, China) by a chemiluminescence detection system (ChemiDoc XRS, Bio-Rad, CA, USA). The quantification of the bands was performed using the ImageJ software (National Institutes of Health, Bethesda, MD, USA).

2.8. Statistical Analysis

Data were analyzed on the Prism software version 8 (GraphPad, San Diego, CA, USA). All values are expressed as the . The significant differences among the groups were evaluated by one-way ANOVA with Tukey’s post hoc test. The results are accepted as the level of significance was set at value of < 0.05.

3. Results

3.1. The Activated Compounds and Potential Targets of the WBT Formula against IPF Are Predicted by Network Pharmacology

Using the TCMSP and TCMID databases, 243 candidate compounds were found in the WBT formula. Among these components, 20, 36, 65, 10, 2, 7, 23, 10, 19, 22, 13, 7, and 9 compounds were identified in Astragalus membranaceus (Fisch.) Bunge, Scutellaria baicalensis Georgi, Bunge Salvia miltiorrhiza Bunge, Polygonum cuspidatum Siebold & Zucc, Angelica sinensis (Oliv.) Diels, Ligusticum striatum DC., Prunus persica (L.) Batsch, Pheretima aspergillum (E. Perrier), Aster tataricus L.f., Tussilago farfara L., Pinellia ternata (Thunb.) Makino, Clematis chinensis Osbeck, and Siegesbeckia orientalis L., respectively (Table 2). Among these compounds, 40 compounds with the repeated names and 38 compounds without targets were deleted to obtain 165 compounds for further analysis. In the WBT formula, 352 targets of these 165 components (Supplementary Table 1) were predicted at TCMID and STPD websites. Then, we input 352 targets into the UniProt Knowledgebase website to unify the standard nomenclature. Meanwhile, 1,531 pathogenic genes-related with IPF were collected from the databases mentioned above (Supplementary Table 2). After the overlapping analysis, 167 intersection targets, as the potential targets of the WBT formula for treating IPF, were acquired and shown in a Venn diagram (Figure 1(a)). Furthermore, two compounds not for 167 targets of IPF were screened out to obtain 163 candidate compounds. In addition, those intersection targets were used to build the network diagram of the components of the WBT formula and their targets, which were composed of 330 nodes and 999 edges (Figure 1(b) and Supplementary Table 3). Collectively, 163 candidate compounds and 167 potential targets of the WBT formula against IPF were predicted by a network pharmacology-based method.

3.2. Potential Pathways and Gene Functional Enrichments for the WBT Formula against IPF

The GO annotation and KEGG pathway enrichments for 167 overlapping targets showed the top 20 significantly enriched terms in the biological process (BP), and cell component (CC) and molecular function (MF) categories are shown in Figure 2(a) and Supplementary Figure 1. As suggested from the results, 167 key targets showed tight relations to the major BP, such as inflammatory response, response to wounding, response to oxidative stress, and regulation of cell adhesion. The WBT formula was reported for integrating multiple signaling pathways in the fibrosis process, inflammatory response, and cancers. The enrichment of target-pathway was built to delve into the mechanisms of potential targets acting on their corresponding signal pathways (Figure 2(b)). The PPI network for therapeutic targets of the WBT formula against IPF was constructed (Figure 2(c)). Furthermore, the therapeutic effect of the WBT formula on IPF was most likely achieved through modulating multiple pathways such as the PI3K-Akt, TNF, and NF-κB pathways. Among these pathways, the NF-κB pathway might be the core potential pathway of the WBT formula for inhibiting the inflammatory response, wounding response, and subsequent fibrosis. The targets of the NF-κB pathway subnetwork were built and are shown in Figure 2(d). Based on the results of the network pharmacology enrichment analysis, we further examined whether the protective effect of the WBT formula against inflammation and subsequent fibrosis were dependent on the regulation of the NF-κB-mediated inflammatory pathways.

3.3. Anti-Inflammatory and Fibrotic Effects of the WBT Formula on BLM-Induced Mice

To further validate the mechanism of the WBT formula for preventing inflammatory response during lung fibrosis, the protocol for animal group, drug administration, and mechanism experiment were designed and are visualized in Figure 3(a). As shown in Figures 3(b) and 3(c), histopathological changes with inflammatory and fibrotic responses due to BLM induction showed the gathering of multiple inflammatory cells and collagen fibrogenesis in alveolar spaces after 7 days. Alveolar infiltration in the model group was suppressed in a dose-dependent manner after treatment with the WBT formula (Figure 3(b)). Masson’s trichrome staining reflected collagen accumulation, as an indicator for lung fibrosis. The newly formed collagen fibers were decreased by the treatment of the WBT formula and PFD. The WBT formula at the dose of 12 g/kg and PFD effectively inhibited the formation of collagen fiber in the alveoli (Figure 3(c)), compared to the BLM-induced group. Subsequently, we examined the content of the IL-6 level in the serum and the levels of IL-1β and TNF-α in BALF from different groups. We found that WBT or PFD treatment significantly inhibited serum IL-6 level, compared with the BLM-induced group (Figure 3(d)). Importantly, the contents of IL-1β and TNF-α in BALF induced by BLM were lower than that of WBT or PFD treatment (Figures 3(e) and 3(f)). In addition, Western blot analysis showed that the expression of FN, a key molecule involved in fibrosis, was upregulated by BLM induction, which was obviously suppressed by WBT or PFD treatment (Figures 3(g) and 3(h)). Together, these data indicated that the WBT formula can inhibit inflammatory and fibrotic responses in the BLM-induced mouse model.

3.4. The WBT Formula Attenuates BLM-Induced Inflammatory Response by Inhibiting the TLR4/MyD88/NF-κB Pathway to Protect Lung Tissue from Fibrosis

To further validate the mechanism of the WBT formula for inflammatory response from the prediction, Western blot and IHC staining methods were used to analyze the NF-κB inflammatory pathway for clarifying the mechanism of the WBT formula. Western blotting revealed that IκBα was reduced by BLM, while NF-κB p-p65/p65 was upregulated in the BLM-induced lung tissues (Figures 4(a)4(c)). WBT treatment promoted the expression of IκBα and inhibited the phosphorylation of NF-κB p65, compared with the BLM group (Figures 4(a)4(c)). Similarly, the increased phosphorylation of NF-κB p65 induced by BLM in lung tissues was inhibited by the WBT formula according to IHC staining (Figure 4(d)). Recent studies have shown that TLR4 activity is crucial for inflammation and PF [23], which is upstream of the NF-κB pathway. MyD88 can be recruited to mediate the release of proinflammatory cytokines when the TLR4 pathway is activated [24, 25]. Therefore, we used Western blot and IHC staining to further investigate the effect of the WBT formula on the levels of TLR4 and MyD88 in the BLM-induced lung tissues. As shown in Figures 5(a)5(c), Western blot analysis showed that BLM increased the expressions of TLR4 and MyD88 in the lung tissues, which were significantly downregulated after WBT formula intervention. For the IHC staining shown in Figure 5(d), similar results were found that WBT treatment significantly suppressed the BLM-induced increases of TLR4 and MyD88. Taken together, the WBT formula can inhibit the TLR4/MyD88/NF-κB signaling pathway to reduce inflammatory response in a BLM-induced mouse model with IPF.

4. Discussion

In the early stages following inflammation, fibroblasts and myofibroblasts provide a tissue scaffold for repairing the injured alveolar epithelial cell [26]. The inflammatory response is the initial pathogenesis of IPF, which stimulates the secretion of profibrotic cytokines to initiate fibroblast activation. Hence, the suppression of inflammation could slow down the IPF process [27]. The potential targets of the WBT formula for inhibiting inflammatory response were predicted using a network pharmacology approach. The WBT formula is a TCM and has extensively achieved efficacy in the Affiliated Hospital of Changchun University of Chinese Medicine (Changchun, China) for the prevention and treatment of IPF. However, its exact mechanism of anti-inflammatory action on IPF is still unclear. This study provides experimental evidence for the inhibitory effect of the WBT formula against inflammatory response by inactivating the TLR4/MyD88/NF-κB pathway, which could be a potential mechanism of the WBT formula treatment for IPF patients (Figure 6).

Some herbs and bioactive components have been shown to possess pharmacologic activities for the treatment of IPF. Astragalus membranaceus (Fisch.) Bunge is a kind of the most widely used traditional Chinese herbal medicines for antioxidant and anti-inflammatory effects [28]. Extracts of Scutellaria baicalensis Georgi and its major chemical constituents have been reported to possess antioxidant and anti-inflammatory functions [29], which promoted wound healing by suppressing the secretion of inflammatory cytokines [28]. After the prediction by network pharmacology, the degrees representing the connections of the nodes with direct neighbors were calculated to estimate the significance of the nodes in a network between the compounds and the targets [30]. In the present study, some targets with higher degree values are important for the network. The top ten degree-ranked compounds are quercetin (MOL000098), xanthinin, luteolin (MOL000006), kaempferol (MOL000422), wogonin (MOL000173), tanshinone IIA (MOL007154), baicalein (MOL002714), beta-sitosterol (MOL000358), 7-O-methylisomucronulatol (MOL000378), and dihydrotanshinlactone (MOL007100), indicating their essential roles of these compounds for the target network of the inflammatory response during IPF. Certainly, those components in the WBT formula need to require LC-MS methods to identify the active compounds in further experiments.

The early stage of IPF is manifested inflammation, especially in the mouse model after BLM for 7 days. According to recent studies and our results, the TLR4/MyD88/NF-κB p65 signaling pathway appears to have potential functions in IPF [18, 31]. In our study, some compounds in the WBT formula have been shown to inhibit the TLR4/MyD88/NF-κB pathway. Researchers have found that quercetin has anti-inflammatory effects and can reduce the inflammatory response in a variety of diseases and inhibit inflammation-related signaling pathways [3234]. Luteolin can inhibit TLR4 activation in multiple diseases [3537]. To be of worth, PI3K/Akt, AGE-RAGE, and other pathways were also predicted by network pharmacology, which is important in the multiple stages of IPF progression and requires further experimental validation.

5. Conclusion

This research is aimed at exploring the mechanism of the WBT formula using a network pharmacology analysis in the early stage of IPF. After the prediction, 163 candidate compounds and 167 key targets of the WBT formula against IPF were identified. The inflammatory response, TNF, and NF-κB pathways were considered the potential pathway of the WBT formula, which was validated using a mouse model of IPF. The results indicated that the WBT formula can suppress IPF progression by inhibiting the TLR4/MyD88/NF-κB pathway-mediated inflammation. This study provides a new insight to the molecular mechanisms of the WBT formula in the application at the clinic.

Data Availability

The datasets used and/or investigated during the current study are available from the corresponding authors upon reasonable request.

Ethical Approval

All animal experiments were approved by the Experimental Animal Administration Committee of Changchun University of Chinese Medicine (Approval No. 2021230) and carried out in accordance with the institutional guidelines.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

ZeW, XL, and DZ conceptualized the research topic and supervised its implementation; SS, JW, GL, LD, YL, HQ, LW, JZ, TC, MZ, ZiW, and YY designed the method, conducted the experiments, analyzed the data, and prepared figures/tables; SS and JW wrote the manuscript; XL and ZeW revised the manuscript. All authors have read and approved the final manuscript. SS and JW contributed equally to this work.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (No. 81804013), the Key Research and Development Project of Jilin Province (No. 20200404057YY), and the Science and Technology Project of the Education Department of Jilin Province (No. JJKH20210964KJ).

Supplementary Materials

Supplementary Figure 1: GO analysis for the 167 shared targets from WBT formula and IPF-related target genes. (a) Molecular function. (b) Cell component. Supplementary Table 1: the possible targets of 165 compounds from the WBT formula. Supplementary Table 2: the target genes of idiopathic pulmonary fibrosis. Supplementary Table 3: a total of 167 potential targets of 163 candidate compounds in the WBT formula were obtained by network pharmacology. (Supplementary Materials)