Network Pharmacology and Experimental Verification Revealed the Mechanism of Yiqi Jianpi Recipe on Chronic Obstructive Pulmonary Disease

Chen, Ke; Zhang, Min; Fan, Liming; Wang, Zhen

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

Evidence-Based Complementary and Alternative Medicine

On this page

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

Special Issue

The Role of Natural Products in Inflammatory Lung Diseases

View this Special Issue

Research Article | Open Access

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

Network Pharmacology and Experimental Verification Revealed the Mechanism of Yiqi Jianpi Recipe on Chronic Obstructive Pulmonary Disease

Ke Chen,¹Min Zhang,¹Liming Fan,¹and Zhen Wang^1,2

Academic Editor: Huantian Cui

Received16 May 2022

Revised19 Jul 2022

Accepted13 Aug 2022

Published07 Sept 2022

Abstract

Objective. The study aimed to explore the active ingredients, targets, and mechanism of action of Yiqi Jianpi recipe (YQJPR) in the treatment of COPD based on the network pharmacology and COPD rat models. Methods. The active ingredients and targets of YQJPR were collected by TCMSP. Disease-related protein targets were obtained from GeneCards. The Venn diagram was used to show the key therapeutic targets of COPD in YQJPR. The PPI network was established by STRING, and cytoHubba plug-in was used to screen the core targets within the network. GO functional enrichment and KEGG pathway enrichment analysis were performed to describe the functions and pathways of the core targets. Cytoscape software was used to construct the ingredient-target network and the core target-enrichment pathway network. The chemical constituents of YQJPR were analyzed by HPLC-MS/MS. Results. The network pharmacology showed 152 active ingredients and 225 targets in YQJPR for the treatment of COPD. The key active ingredients were quercetin, luteolin, kaempferol, tanshinone IIA, and baicalein. The contents of quercetin and luteolin in YQJPR were quantitatively measured by HPLC-MS/MS. 22 core genes were screened, including AKT1, IL-6, JUN, VEGFA, and CASP3, which were mainly involved in BPs such as cell proliferation and differentiation, oxidative/chemical stress, and regulation of DNA-binding transcription factor activity and regulated viral infection, tumor, HIF-1, MAPK, TNF, and IL-17 pathways. Animal experiments showed that YQJPR could significantly reduce the expression of p-ERK1/2, p-Akt, c-Myc, cleaved caspase-3, and p-Stat3 in lung tissue (). HE staining showed that, compared with the model group, YQJPR significantly improved lung tissue morphology and reduced lung inflammation in rats. Conclusion. The effects of YQJPR on COPD may involve multiple components, pathways, and targets. This study provides new ideas for further and more comprehensive exploration of the therapeutic effect of YQJPR on COPD in the future.

1. Introduction

Chronic obstructive pulmonary disease (COPD) is currently the third leading cause of death in the world after cardiovascular disease and stroke [1], which brings a huge economic and social burden to people and seriously affects the quality of life. The World Economic Forum estimates that the global cost of COPD will reach $50 trillion per year by 2030, more than the cost of cardiovascular disease [2]. The commonly used and effective treatment methods for COPD are mainly inhaled drugs, especially bronchodilators. Nondrug treatments, such as pulmonary rehabilitation, oxygen inhalation, and other adjuvant therapy, are also effective and are current research hotspots, but their costs are high. These treatments are underutilized worldwide due to insufficient funding, resources, and patient awareness [3]. To date, no therapy has been found to slow the progression of COPD and no drug has significantly reduced its mortality or complication rates [2]. Therefore, it is necessary to explore new treatment methods to improve the quality of life of patients.

Yiqi Jianpi recipe (YQJPR) is an empirical prescription for the clinical treatment of COPD by Professor Wang, a famous Chinese medicine doctor from the Zhejiang Provincial Hospital of Traditional Chinese Medicine. It is mainly composed of Huang Qi (Radix Astragalus), Dang Shen (Codonopsis), Bai Zhu (Atractylodes macrocephala Koidz), Fu Lin (Poria cocos), Zhu Li (Bamboo Juice), Ban Xia (Pinelliae Rhizoma), Chen Pi (Pericarpium Citri Reticulatae), Qian Hu (Peucedani Radix), Dan Pi (Cortex Moutan), and Gan Cao (Glycyrrhizae Radix). Although YQJPR has a curative effect on COPD, its composition is complex and its pharmacodynamic mechanism is not yet clear. Therefore, there is great interest in studying active compounds and molecular targets for the treatment of COPD in YQJPR.

Network pharmacology is a recently proposed new method to study traditional Chinese medicine (TCM) by combining systematic network analysis and pharmacology to elucidate the interactions between compounds, genes, proteins and diseases [4]. Therefore, the active ingredients, targets, and mechanism of action of YQJPR in the treatment of COPD through network pharmacology, which can provide a theoretical basis for future experimental research and clinical application, were explored.

2. Materials and Methods

2.1. Screening of Active Ingredients and Targets of YQJPR

The chemical compositions of YQJPR were obtained from the TCMSP database (https://old.tcmsp-e.com/tcmsp.php). This database provides 12 pharmacokinetic properties [5], including oral bioavailability (OB), drug-likeness (DL), drug half-life (HL), and blood-brain barrier (BBB). OB is one of the most important pharmacokinetic parameters for studying oral drugs [6], and high OB is often a key indicator for determining the drug-like properties of bioactive molecules as therapeutics. It is generally considered that DL ≥ 0.18 of a compound is the selection criterion for “drug-like” compounds in traditional Chinese herbal medicine [7]. Therefore, compounds with OB ≥ 30% and DL ≥ 0.18 as active ingredients were selected for this study. Then, the targets of each herb were collected in the TCMSP database and converted into official genetic symbol through UniProt database (https://sparql.uniprot.org/).

2.2. Screening of Targets Related to COPD

COPD-related protein targets were provided by the GeneCard database (https://www.genecards.org/). Using “chronic obstructive pulmonary disease” as the keyword, the genes with protein-coding function and GeneCards Inferred Functionality Scores (GIFtS) ≥30 were selected.

2.3. Venn Diagram

The Venn diagram (made by Venny 2.1.0) was used to obtain the intersection genes of COPD and YQJPR, which were defined as “key therapeutic targets for COPD.”

2.4. Construction of the Ingredient-Target Network

Using the Cytoscape 3.9.1 software to link drug and disease intersection genes, the relationship between active ingredients and COPD-related targets in YQJPR was visualized.

2.5. Preparation of YQJPR

The preparation of YQJPR sample was performed as Dong et al. [8] and improved. Ten herbs were included in YQJPR, which were Huang Qi 30 g, Dang Shen 20 g, Bai Zhu 12 g, Fu Lin 15 g, Zhu Li Ban Xia 12 g, Chen Pi 9 g, Qian Hu 12 g, Dan Pi 12 g, and Gan Cao 3 g. All the above herbs (125 g in total) were mixed with boiling water and boiled twice (1.5 h for the first time, the amount of water was 8 times the quality of the herbs that was 1000 ml; 0.5 h for the second time, the amount of water was 6 times that was 750 ml). The two decoctions were filtered, mixed together, concentrated to 400 ml, and then frozen at −80°C overnight. After being placed in a freeze dryer for about 72 hours, the sample of YQJPR freeze-dried powder was obtained, and the powder yield was 30%. The lyophilized powder was dissolved in methanol as the test drug (namely, YQJPR) for further study. The contents of quercetin and luteolin in YQJPR were quantitatively measured by high-performance liquid chromatographic-mass spectrometry/mass spectrometry (HPLC-MS/MS). HPLC experiments were performed on an Agilent LC1290-QQQ-6470 instrument. The samples were filtered through 0.22 μm membrane before injection.

2.5.1. Chromatographic Conditions

HPLC separation was performed on an ACE UltraCore C18 column (2.1 × 75 mm, 2.5 μm) at a column temperature of 40°C. Gradient elution was conducted with 0.1% formic acid (A) and acetonitrile (B): 0 min, 90% A; 2.5 min, 50% A; and 5 min, 10% A. The solvent flow rate was 350 μL·min⁻¹, and the injection volume was 3 μL.

2.5.2. Mass Spectrometry Conditions

The MS data were acquired in negative ionization mode and multiple reaction monitoring scan mode. Instrumental parameters were optimized as follows: capillary: 3500 V (positive) and −3000 V (negative); nozzle voltage: −1000 V (negative); source gas (nebulizer gas), 30 psi; gas temperature, 150°C; gas flow, 12 L·min⁻¹; sheath gas temperature, 350°C; sheath gas flow, 12 L·min⁻¹.

2.6. Construction of the PPI Network

Using the STRING database (https://cn.string-db.org/), the protein-protein interaction (PPI) network was established for key therapeutic targets in COPD by setting “Homo sapiens” in “Organism” and “moderate confidence 0.4.” The obtained PPI network data were imported into Cytoscape 3.9.1, and then, the important nodes and subnetworks in the network were screened by the cytoHubba plugin. CytoHubba predicts and explores important nodes and subnetworks in the network through several topological analysis methods [9]. In this study, two analysis methods, degree and closeness, were selected to rank nodes, select core genes, and construct subnetworks.

2.7. GO and KEGG Pathway Enrichment Analysis

Gene ontology (GO) analysis describes the function of genes from three aspects: biological process (BPs), cellular components (CCs), and molecular functions (MFs), and group genes with similar functions into one group [10]. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis reflects the interactions of genes or gene products in pathways to predict cellular function [11].

The R packages including “colorspace,” “string,” and “ggplot2” and the Bioconductor package (https://www.bioconductor.org/) including “DOSE,” “clusterProfiler,” and “enrichplot” were used to perform GO and KEGG pathway analysis on the core genes of YQJPR and COPD [12]. According to “, q ≤ 0.05,” the top 20 GO terms and KEGG pathways were selected. Subsequently, the data of the top 20 KEGG pathways were imported into the Cytoscape software to construct a KEGG network.

2.8. In Vivo Experiments in Rats

2.8.1. Animals and Reagents

SD rats (male, 4–6 weeks old, 200 ± 20 g) were purchased from Zhejiang Chinese Medical University Animal Center (Grade SPF). All rats had free access to food and water and were tested in the SPF barrier laboratory of Zhejiang Chinese Medical University Animal Center. The whole protocol was approved by the Animal Ethics Committee of Zhejiang Chinese Medical University. Every effort was made to minimize the suffering of the animals during the experiment.

The main formula of YQJPR is as follows: Huang Qi, 30 g; Dang Shen, 20 g; Bai Zhu, 12 g; Fu Lin, 15 g; Zhu Li Ban Xia, 12 g; Chen Pi, 9 g; Qian Hu, 12 g; Dan Pi, 12 g; and Gan Cao 3 g. All herbs were provided by the TCM Pharmacy of the first affiliated hospital of Zhejiang Chinese Medical University (Hangzhou, China), and one dose was finally decocted to 125 ml of liquid (1 g/ml).

Hongtashan cigarettes (84-mm flue-cured cigarettes, smoke nicotine content: 0.9 mg, tar content: 10 mg, smoke carbon monoxide content: 11 mg) were provided by Hongta Tobacco Co., Ltd.

Antibody types are as follows: primary antibodies [Phospho-Akt (CST, 9271), Phospho-p44/42 MAPK (p-ERK1/2) (CST, 4370), c-Myc (CST, 5605), cleaved caspase-3 (CST, 9664), Phospho-Stat3 (CST, 9134), and GAPDH (BOSTER, BM1623)], and secondary antibodies [anti-rabbit IgG, HRP-linked antibody (CST, 7074) and anti-mouse IgG, HRP-linked Antibody (CST, 7076)].

2.8.2. Grouping and Intervention

After three days of adaptive feeding, ten of the 30 SD rats were randomly selected as the control group, and the rest were selected as the model establishment group (Figure 1). Referring to previous experience [13, 14], the model establishment group (n = 20) used cigarette smoke exposure (CSE) combined with intratracheal instillation of LPS (Sigma) to prepare a COPD model for 12 weeks. During this period, 20 rats passively inhaled smoke (1.5 cigarettes/rat, 1 h/day) in a self-made smoking box (approximately 46 cm × 33 cm × 26 cm, made of plastic), except for the intratracheal instillation of 0.2 ml (200 ug) LPS (no smoke exposure) on days 1 and 14. From Week 9 onward, the model establishment group was randomly divided into the model group (n = 10) and the YQJPR group (n = 10). The model group continued to passively smoke until the end of 12 weeks, and the YQJPR group started TCM intervention and continued to passively smoke 1 hour after gavage. The YQJPR group continued the TCM intervention after the 12th week of modelling, and the whole TCM intervention lasted for 12 weeks. The dose of YQJPR to rats was calculated from the dose conversion method between humans and rats [15]. The adult clinical dosage of YQJPR was 1.786 g/(kg·d), and then, the YQJPR group was given 11.25 g/(kg·d) by gavage, while the control and model groups were given the same amount of normal saline by gavage. After treatment, all rats were fasted overnight and anesthetized by intraperitoneal injection of 2.5% sodium pentobarbital (0.3 ml/100 g) the next day. The right part of the lung tissue was stored at −80°C for western blot, and the left part was fixed with 4% paraformaldehyde for HE staining to observe the morphological structure of the lung tissue.

2.8.3. Western Blot

The lung tissue was homogenized with RIPA lysis buffer (Boster, Wuhan, China) and centrifuged (, 10 min, 4°C) to collect total protein. After denaturation with SDS-PAGE Sample Loading Buffer (Beyotime, Shanghai, China), proteins were separated by 10% or 12% SDS-PAGE and transferred to PVDF membranes (Boster, Wuhan, China). The membranes were blocked with 5% nonfat dry milk (Dawen Biotec, Hangzhou, China) at room temperature. After 1.5 h, they were rinsed with TBST and incubated overnight at 4°C with primary antibodies against p-ERK1/2 (1 : 1000, anti-rabbit), p-Akt (1 : 1000, anti-rabbit), c-Myc (1 : 500, anti-rabbit), cleaved caspase-3 (1 : 1000, anti-rabbit), p-Stat3 (1 : 500, anti-rabbit), and GAPDH (1 : 5000, anti-mouse). Then, the blots were washed with TBST (10 min3). The blots were subsequently incubated with a secondary antibody (1 : 1000, goat anti-rabbit or goat anti-mouse antibody) at room temperature for 1.5 h. Finally, the blots were washed with TBST (10 min3) and visualized with enhanced chemiluminescence reagent (Fdbio science, Hangzhou, China) and analyzed with Bio-Rad Machine (USA).

ImageJ software (version 1.8.0) analyzed the gray value of the blots and took the ratio with the gray value of GAPDH as the relative expression level, and then, the ratio was normalized.

2.8.4. Statistical Analysis

GraphPad Prism 9 was used for statistics. All experimental data were obtained from at least three repeated independent trials and expressed as mean ± SD. Statistical differences between the groups were identified by one-way ANOVA. indicates a statistical difference.

3. Results

3.1. Active Ingredient-Target Analysis

Searched in the TCMSP database and screened with OB ≥ 30% and DL ≥ 0.18 as the condition, a total of 208 active ingredients of YQJPR were obtained, including 20 in Radix Astragalus, 21 in Codonopsis, 7 in Atractylodes macrocephala Koidz, 15 in Poria cocos, 13 in Pinelliae Rhizoma, 5 in Pericarpium Citri Reticulatae, 24 in Peucedani Radix, 11 in Cortex Moutan, and 92 in Glycyrrhizae Radix. Bamboo juice was filtered out during this screening process (Supplementary Table 1). The targets of each herb were collected, screened, and normalized by the TCMSP and UniProt databases. A total of 2926 drug targets were obtained after removing duplicate targets, which were summarized by drug-ingredient-target correspondence.

3.2. Acquisition and Selection of COPD-Related Targets

With “chronic obstructive pulmonary disease” as the keyword, the genes with protein-coding function and GIFtS ≥30 were selected as the targets in the GeneCards database, and 5057 targets were finally collected.

3.3. Venn Diagram

We constructed Venn diagrams (Figure 2) by Venny 2.1.0 to depict the intersection targets of COPD and YQJPR (the key therapeutic targets for COPD). The analysis revealed that there were 225 targets related to both COPD and YQJPR (Supplementary Table 2).

3.4. Ingredient-Target Network Analysis

The ingredient-target network (Figure 3) was constructed by the Cytoscape 3.9.1 software to visualize the relationship between 152 active ingredients and 225 targets. The network included 377 nodes and 1931 edges. These edges indicated interactions between ingredients and targets. The active ingredients and the targets with greater correlation were screened by the cytoHubba plug-in. The degree value represented the number of edges associated with the node. The greater the degree value, the greater the correlation between the ingredient and the target.

(a)

(b)

The results showed that in this network (Table 1), MOL000098 (quercetin), MOL000006 (luteolin), MOL000422 (kaempferol), MOL007154 (tanshinone IIA), MOL000378 (7-O-methylisomucronulatol), MOL003896 (7-methoxy-2-methyl isoflavone), and MOL002714 (baicalein) were associated with at least 30 targets, and PTGS2, ESR1, HSP90AA1, AR, NOS2, PPARG, NCOA2, PRSS1, SCN5A, CDK2, PTGS1, and GSK3B were regulated by at least 60 compounds.

3.5. Quality Evaluation of YQJPR

The HPLC-MS/MS parameters of the standards after optimization are shown in Table 2. Two key compounds (quercetin and luteolin) in YQJPR were well identified by comparing their retention times and mass spectra with the standards by using HPLC-MS/MS analysis (Figure 4), whose contents in YQJPR freeze-dried powder were 70.2 and 60.9 ng/g. This result was calculated by quantitative analysis.

(a)

(b)

3.6. PPI Network Construction

The PPI network (Figure 5) was constructed from the same 225 key therapeutic targets for COPD, including 225 nodes and 4224 edges, with the average node degree of 37.5 and the average local clustering coefficient of 0.594.

The nodes of the PPI network were calculated by the cytoHubba plug-in and were ranked according to degree and closeness. The top 22 targets with “degree > 122, closeness > 150” in the PPI network were selected (Figure 6(a)), which were AKT1, IL-6, JUN, VEGFA, CASP3, IL-1B, ESR1, PTGS2, EGFR, MYC, HIF1A, MAPK3, STAT3, etc. These targets were suggested to be the core genes of YQJPR in treating COPD, and then, a subnetwork (Figure 6(b)) was constructed for the core genes.

(a)

(b)

3.7. GO and KEGG Analysis

GO and KEGG enrichment analysis of 22 core genes was performed to determine the underlying molecular mechanism of YQJPR in the treatment of COPD. Results for “, q ≤ 0.05” were retained, and bubble charts were drawn for the top 20 clearly enriched entries.

3.7.1. GO Analysis

The top 20 BPs with the most enriched core genes were analyzed (Table 3 and Figure 7).

3.7.2. KEGG Pathway Analysis

KEGG analysis showed that the relevant YQJPR signaling pathways in the treatment of COPD included virus infection, TNF, IL-17, HIF-1, MAPK, AGE-RAGE, and cancer pathways (Supplementary Table 3). The top 20 pathways with the most enriched core genes for the treatment of COPD are selected in Figure 8.

3.7.3. KEGG Network Analysis

The KEGG network was constructed to visualize the relationship between the top 20 pathways and enriched genes through the Cytoscape software (Figure 9). The results showed that MAPK3, AKT1, EGFR, JUN, FOS, STAT3, CCND1, VEGF, CASP3, and EGF were associated with at least 10 different pathways.

3.8. Experimental Verification

3.8.1. Histopathological Examination

To determine whether the COPD model was established successfully, HE staining of rat lung tissue was performed, mainly to observe the structure of the trachea and alveoli. Figure 10 shows that the COPD rat model was successfully established. In the control group, the lung tissue structure was intact, and there was no obvious inflammatory cell infiltration in the alveolar cavity. Compared with the control group, the model group showed alveolar expansion and fusion and obvious inflammatory cell infiltration around the trachea. Compared with the model group, the lung tissue morphology of the YQJPR group was improved, and the peritracheal inflammatory infiltration was significantly reduced.

3.8.2. Western Blot of Lung Tissue

In order to identify the results of the network pharmacology analysis, five targets were selected: MAPK3, AKT1, MYC, STAT3, and CASP3 for protein validation (Figure 11).

(a)

(b)

The results showed that compared with the control group, the expressions of p-ERK1/2, p-Akt, c-Myc, cleaved caspase-3, and p-Stat3 in the lung tissue of the model group were significantly increased. Compared with the model group, YQJPR significantly inhibited the expressions of p-ERK1/2, p-Akt, c-Myc, cleaved caspase-3, and p-Stat3.

4. Discussion

COPD is a preventable and treatable progressive inflammatory disease of the airway, alveoli, and microvascular system characterized by respiratory symptoms with incomplete and persistent airflow restriction [16]. The pathogenesis of COPD is complex, and the existing treatments are difficult to meet expectations, so it is urgent to explore new treatment methods. YQJPR is an empirical prescription with decades of clinical experience by Professor Wang, a famous Chinese medicine doctor from Zhejiang. YQJPR has unique advantages in the treatment of COPD, but it is difficult to popularize and use due to the complex ingredients of the TCM compound. Network pharmacology provides a new approach for the research of the pharmacological mechanism of TCM, and we try to use it to predict the components, targets, and mechanisms of the TCM compounds.

In this study, 152 ingredients and 225 targets of YQJPR related to COPD were identified, indicating that YQJPR is a treatment of COPD through combined regulation of multiple ingredients and multiple targets. Among them, quercetin, luteolin, kaempferol, tanshinone IIA, and baicalein regulated at least 30 targets and were considered as key active ingredients. And it was verified that quercetin and luteolin did exist in YQJPR by HPLC-MS/MS. The above five compounds are all flavonoids with various pharmacological effects, which include antibacterial, anti-inflammatory, antioxidant, and anticancer. Studies have shown that a diet rich in quercetin can reduce the risk of asthma, bronchial hyperresponsiveness, and COPD [17, 18]. In vitro and in vivo experiments have shown that quercetin improves lung function in COPD mice and significantly reduces oxidative stress, lung inflammation, and goblet cell metaplasia induced by cigarette smoke exposure [19, 20]. Luteolin can modulate multiple signaling pathways (such as NF-κB, JAK-STAT, and TLR) and inhibit inflammatory factors (such as IL-1β, IL-6, IL-8, IL-17, IL-22, TNF-α, and COX-2) to exert anti-inflammatory effects [21]. Luteolin can prevent NF-κB pathway activation and inhibit apoptosis and inflammation in vitro and in vivo by downregulating miR-132 [22]. Kaempferol can induce apoptosis and downregulate the expression of AKT/PI3K phosphorylation, ERK pathway, and MAPK, and inhibit the growth of A549 cells in a concentration-dependent manner [23]. Tanshinone IIA inhibits inflammation and oxidative stress induced by cigarette smoke by blocking MAPK/HIF-1α pathway, effectively alleviating the progression of COPD [24, 25]. Baicalein mainly has anti-inflammatory and anticancer effects and can treat pneumonia caused by the streptococcus pneumoniae infection by reducing the permeability of the alveolar capillary membrane, inhibiting pulmonary inflammatory response and apoptosis [26]. Baicalein can also significantly reduce hemodynamic parameters, inhibit pulmonary remodeling, reduce apoptosis and inflammation, and effectively improve lung injury in pulmonary hypertension model [27].

Through constructing PPI network of intersection targets between YQJPR and COPD, 22 core genes, including AKT1, IL-6, JUN, VEGFA, CASP3, IL-1B, PTGS2, MYC, MAPK3, and STAT3, were identified. IL-1 and IL-6 are both secreted proteins associated with inflammation and immune responses, with IL-1 mainly being produced by pathogens that stimulate cells of the innate immune system, such as monocytes and macrophages. IL-1 is highly expressed in small airway epithelial cells, which inhibits CSE-induced cell proliferation and promotes apoptosis, but can also promote inflammation and aggravate cell damage [28, 29]. IL-6 is a multifunctional cytokine with various roles in immunity, development, metabolism, aging, and cancer and can stimulate B cells to produce antibodies [30]. Jun and Fos are the main components of transcriptional activator protein 1 (AP-1), which are closely related to cell proliferation, apoptosis, and inflammatory response [31–33]. Caspase-3 is a key executor of apoptosis [34], and STAT3 is a signaling protein closely related to cell growth, differentiation, and survival [35]. MAPKs are a family of serine/threonine protein kinases that can regulate various extracellular stimulus signals through four linear cascades composed of MAPKs to regulate corresponding cellular processes, in which ERK1/2 (MAPK1/3) cascade preferentially regulates proliferation, differentiation, and migration [36]. MYC is a multifunctional transcription factor that can regulate the expression of many genes involved in biological processes such as cell growth, proliferation, apoptosis, and metabolism [37]. COX2 (encoded by PTGS2) is a product of multiple stimuli, linked to inflammation and tumor cell growth [38]. It is hardly expressed in normal cells but is expressed in highly proliferating cells. COX2 can activate macrophages or other inflammatory cells to transform the inflammatory site into a precancerous microenvironment, which has carcinogenic effects [39].

To understand the molecular mechanism of YQJPR in the treatment of COPD, GO analysis and KEGG enrichment analysis of 22 core genes were performed. GO analysis revealed that core target genes were strongly associated with the following BPs: cell proliferation and differentiation, oxidative/chemical stress, glandular development, regulation of DNA-binding transcription factor activity, and response to metal ions. The analysis of the KEGG pathway was primarily related to TNF, IL-17, HIF-1, MAPK, AGE-RAGE, virus infection, and cancer signaling pathways. MAPK is mainly involved in cell proliferation, differentiation, and immune regulation. It has been shown that the p38 MAPK pathway can be activated by CSE to stimulate pro-inflammatory cytokine expression [40]. Studies have found that the TNF and IL-17 signaling pathways play important roles in the inflammatory response. TNF-α, the most important cytokine in the TNF pathway, is elevated in COPD and can induce the release of other pro-inflammatory factors [29]. IL-17A exacerbates inflammation and mucus production in the pathogenesis of COPD in an autocrine way [41]. IL-17, a signature cytokine secreted by Th17 cells, plays a major role in the IL-17 pathway in a pro-inflammatory manner. It controls inflammation by regulating the expression of inflammatory genes in cells in most nonhematopoietic organs, such as IL-1, IL-6, IL-8, TNF, CCL2, and MMP1/3/9/13 [42]. The pulmonary circulation is chronically hypoxic in COPD [43], leading to HIF-1 signaling pathway activation and upregulated HIF-1α expression. The HIF-1 pathway is a major regulator of angiogenesis, an adaptive response to tissue hypoxia [44, 45]. HIF-1α accumulation directly upregulates VEGF levels, which contribute to airway and vascular remodeling, to improve the oxygen supply to the pulmonary circulation [46]. The results suggest that YQJPR may act on the target genes to treat COPD through the above-mentioned pathways.

Western blotting was used to verify five key targets (MAPK3, AKT1, MYC, STAT3, and CASP3) that were predicted in the PPI network. According to the results of the BPs and the KEGG analysis, these five genes were mainly enriched in the virus infection, cancer, MAPK, HIF-1, TNF, and IL-17 signaling pathways. It is closely related to biological processes such as cell proliferation and differentiation, oxidative/chemical stress, response to metal ion, and regulation of DNA-binding transcription factor activity. The results showed that compared with the model group, YQJPR significantly downregulated the expression of p-ERK1/2, p-Akt, c-Myc, cleaved caspase-3, and p-Stat3 in lung tissue. Therefore, the molecular mechanism of YQJPR treatment of COPD has been more clearly verified.

5. Conclusions

In this study, the key active ingredients, key targets, and pathways of YQJPR in the treatment of COPD were selected by network pharmacology. The results show that YQJPR can interfere with the development of COPD through multiple ingredients, targets, and pathways. In conclusion, YQJPR effectively improved lung inflammation in COPD rats. Especially, quercetin, luteolin, kaempferol, tanshinone IIA, and baicalein were the main effectors in YQJPR, which primarily downregulated the expressions of five key targets (MAPK3, AKT1, MYC, STAT3, and CASP3), and subsequently inhibited the activation of virus infection, cancer, MAPK, HIF-1, TNF, and IL-17 pathways, thereby alleviating CS/LPS-induced lung inflammation.

The disadvantage of this study is that only one detection method was used in the COPD rat model, and the pathways identified above were not verified, so further in-depth and comprehensive studies are needed. This study provides ideas and a basis for exploring the therapeutic effect of YQJPR on COPD in the future.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Ke Chen and Min Zhang contributed equally to this work.

Acknowledgments

This work was supported by the National Key R&D Program of China (2018YFC2002500) and the National Natural Science Foundation of China (grant no. 82174302). All the authors of the manuscript are immensely grateful to the First Central Laboratory of Zhejiang Chinese Medical University for providing the experimental platform.

Supplementary Materials

Supplementary Table 1: information of collected YQJPR active ingredients. Supplementary Table 2: results of Venn diagram. Supplementary Table 3: data of KEGG enrichment analysis. (Supplementary Materials)

References

GBD Chronic Respiratory Disease Collaborators, “Prevalence and attributable health burden of chronic respiratory diseases, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017,” Lancet Respiratory Medicine, vol. 8, no. 6, pp. 585–596, 2020 Jun.
View at: Publisher Site | Google Scholar
P. J. Barnes, “Copd 2020: new directions needed,” American Journal of Physiology-Lung Cellular and Molecular Physiology, vol. 319, no. 5, pp. L884–L886, 2020 Nov 1.
View at: Publisher Site | Google Scholar
C. L. Rochester, I. Vogiatzis, A. E. Holland et al., “An official American thoracic society/European respiratory society policy statement: enhancing implementation, use, and delivery of pulmonary rehabilitation,” American Journal of Respiratory and Critical Care Medicine, vol. 192, no. 11, pp. 1373–1386, 2015.
View at: Publisher Site | Google Scholar
K. Wang, L. Lei, J. Cao et al., “Network pharmacology-based prediction of the active compounds and mechanism of Buyang Huanwu Decoction for ischemic stroke,” Experimental and Therapeutic Medicine, vol. 22, no. 4, p. 1050, 2021.
View at: Publisher Site | Google Scholar
J. Ru, P. Li, J. Wang et al., “TCMSP: a database of systems pharmacology for drug discovery from herbal medicines,” Journal of Cheminformatics, vol. 6, no. 1, p. 13, 2014.
View at: Publisher Site | Google Scholar
X. Xu, W. Zhang, C. Huang et al., “A novel chemometric method for the prediction of human oral bioavailability,” International Journal of Molecular Sciences, vol. 13, no. 6, pp. 6964–6982, 2012.
View at: Publisher Site | Google Scholar
W. Tao, X. Xu, X. Wang et al., “Network pharmacology-based prediction of the active ingredients and potential targets of Chinese herbal Radix Curcumae formula for application to cardiovascular disease,” Journal of Ethnopharmacology, vol. 145, no. 1, pp. 1–10, 2013.
View at: Publisher Site | Google Scholar
B. Dong, C. Peng, P. Ma, and X. Li, “An integrated strategy of MS-network-based offline 2DLC-QTOF-MS/MS coupled with UHPLC-QTRAP®-MS/MS for the characterization and quantification of the non-polysaccharides in Sijunzi decoction,” Analytical and Bioanalytical Chemistry, vol. 413, no. 13, pp. 3511–3527, 2021.
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, 2014.
View at: Publisher Site | Google Scholar
D. W. Huang, B. T. Sherman, and R. A. Lempicki, “Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources,” Nature Protocols, vol. 4, no. 1, pp. 44–57, 2009.
View at: Publisher Site | Google Scholar
M. Kanehisa and Y. Sato, “KEGG Mapper for inferring cellular functions from protein sequences,” Protein Science, vol. 29, no. 1, pp. 28–35, 2020.
View at: Publisher Site | Google Scholar
G. Yu, L. G. Wang, Y. Han, and Q. Y. He, “Clusterprofiler: an R package for comparing biological themes among gene clusters,” OMICS: A Journal of Integrative Biology, vol. 16, no. 5, pp. 284–287, 2012.
View at: Publisher Site | Google Scholar
D. Wang, Q. Du, H. Li, and S. Wang, “The isosteroid alkaloid imperialine from bulbs of fritillaria cirrhosa mitigates pulmonary functional and structural impairment and suppresses inflammatory response in a COPD-like rat model,” Mediators of Inflammation, vol. 2016, Article ID 4192483, 17 pages, 2016.
View at: Publisher Site | Google Scholar
J. Z. Weng, Y. Wang, and T. Y. Sun, “Cathelicidin LL-37 restoring glucocorticoid function in smoking and lipopolysaccharide-induced airway inflammation in rats,” Chinese Medical Journal, vol. 132, no. 5, pp. 569–576, 2019.
View at: Publisher Site | Google Scholar
Z. Yan, C. Huang, G. Huang et al., “The effect of Jiedu Huoxue decoction on rat model of experimental nonbacterial prostatitis via regulation of miRNAs,” Pharmaceutical Biology, vol. 58, no. 1, pp. 745–759, 2020.
View at: Publisher Site | Google Scholar
K. F. Rabe and H. Watz, “Chronic obstructive pulmonary disease,” The Lancet, vol. 389, no. 10082, pp. 1931–1940, 2017.
View at: Publisher Site | Google Scholar
P. Knekt, J. Kumpulainen, R. Jarvinen et al., “Flavonoid intake and risk of chronic diseases,” The American Journal of Clinical Nutrition, vol. 76, no. 3, pp. 560–568, 2002.
View at: Publisher Site | Google Scholar
C. Tabak, I. C. Arts, H. A. Smit, D. Heederik, and D. Kromhout, “Chronic obstructive pulmonary disease and intake of catechins, flavonols, and flavones: the MORGEN Study,” American Journal of Respiratory and Critical Care Medicine, vol. 164, no. 1, pp. 61–64, 2001.
View at: Publisher Site | Google Scholar
S. Ganesan, A. N. Faris, A. T. Comstock et al., “Quercetin prevents progression of disease in elastase/LPS-exposed mice by negatively regulating MMP expression,” Respiratory Research, vol. 11, no. 1, p. 131, 2010.
View at: Publisher Site | Google Scholar
N. P. da Silva Araujo, N. A. de Matos, S. Leticia Antunes Mota et al., “Quercetin attenuates acute lung injury caused by cigarette smoke both in vitro and in vivo,” COPD: Journal of Chronic Obstructive Pulmonary Disease, vol. 17, no. 2, pp. 205–214, 2020.
View at: Publisher Site | Google Scholar
F. Gendrisch, P. R. Esser, C. M. Schempp, and U. Wolfle, “Luteolin as a modulator of skin aging and inflammation,” BioFactors, vol. 47, no. 2, pp. 170–180, 2021.
View at: Publisher Site | Google Scholar
X. Liu and J. Meng, “Luteolin alleviates LPS-induced bronchopneumonia injury in vitro and in vivo by down-regulating microRNA-132 expression,” Biomedicine & Pharmacotherapy, vol. 106, pp. 1641–1649, 2018.
View at: Publisher Site | Google Scholar
M. Imran, B. Salehi, J. Sharifi-Rad et al., “Kaempferol: a key emphasis to its anticancer potential,” Molecules, vol. 24, no. 12, p. 2277, 2019.
View at: Publisher Site | Google Scholar
D. Li, D. Sun, L. Yuan et al., “Sodium tanshinone IIA sulfonate protects against acute exacerbation of cigarette smoke-induced chronic obstructive pulmonary disease in mice,” International Immunopharmacology, vol. 81, Article ID 106261, 2020.
View at: Publisher Site | Google Scholar
R. Guan, J. Wang, Z. Li et al., “Sodium tanshinone IIA sulfonate decreases cigarette smoke-induced inflammation and oxidative stress via blocking the activation of MAPK/HIF-1α signaling pathway,” Frontiers in Pharmacology, vol. 9, p. 263, 2018.
View at: Publisher Site | Google Scholar
G. Qinghe, L. Guangyu, L. Xiaobin, and H. Qingsong, “Baicalein exerts a protective role in pneumonia caused by Streptococcus pneumoniae,” Frontiers in Bioscience, vol. 24, no. 5, pp. 849–858, 2019.
View at: Publisher Site | Google Scholar
R. Shi, Z. Wei, D. Zhu et al., “Baicalein attenuates monocrotaline-induced pulmonary arterial hypertension by inhibiting vascular remodeling in rats,” Pulmonary Pharmacology & Therapeutics, vol. 48, pp. 124–135, 2018.
View at: Publisher Site | Google Scholar
G. Lopez-Castejon and D. Brough, “Understanding the mechanism of IL-1β secretion,” Cytokine & Growth Factor Reviews, vol. 22, no. 4, pp. 189–195, 2011.
View at: Publisher Site | Google Scholar
T. Liu, “miR-937 serves as an inflammatory inhibitor in cigarette smoke extract-induced human bronchial epithelial cells by targeting IL1B and regulating TNF-α/IL-17 signaling pathway,” Tobacco Induced Diseases, vol. 19, no. July, pp. 1–11, 2021.
View at: Publisher Site | Google Scholar
S. Kang, M. Narazaki, H. Metwally, and T. Kishimoto, “Historical overview of the interleukin-6 family cytokine,” Journal of Experimental Medicine, vol. 217, no. 5, Article ID e20190347, 2020.
View at: Publisher Site | Google Scholar
M. Castellazzi and A. Sergeant, “The C-Jun oncoprotein,” Bull Cancer, vol. 80, no. 9, pp. 757–759, 1993.
View at: Google Scholar
S. Leppa and D. Bohmann, “Diverse functions of JNK signaling and c-Jun in stress response and apoptosis,” Oncogene, vol. 18, no. 45, pp. 6158–6162, 1999.
View at: Publisher Site | Google Scholar
H. B. Schonthaler, J. Guinea-Viniegra, and E. F. Wagner, “Targeting inflammation by modulating the Jun/AP-1 pathway,” Annals of the Rheumatic Diseases, vol. 70, no. 1, pp. i109–i112, 2011.
View at: Publisher Site | Google Scholar
S. Shalini, L. Dorstyn, S. Dawar, and S. Kumar, “Old, new and emerging functions of caspases,” Cell Death & Differentiation, vol. 22, no. 4, pp. 526–539, 2015.
View at: Publisher Site | Google Scholar
A. C. Guanizo, C. D. Fernando, D. J. Garama, and D. J. Gough, “STAT3: a multifaceted oncoprotein,” Growth Factors, vol. 36, no. 1-2, pp. 1–14, 2018.
View at: Publisher Site | Google Scholar
G. Maik-Rachline, I. Wortzel, and R. Seger, “Alternative splicing of MAPKs in the regulation of signaling specificity,” Cells, vol. 10, no. 12, p. 3466, 2021.
View at: Publisher Site | Google Scholar
A. S. Farrell and R. C. Sears, “MYC degradation,” Cold Spring Harbor Perspectives in Medicine, vol. 4, no. 3, Article ID a014365, 2014.
View at: Publisher Site | Google Scholar
A. Alexanian and A. Sorokin, “Cyclooxygenase 2: protein-protein interactions and posttranslational modifications,” Physiological Genomics, vol. 49, no. 11, pp. 667–681, 2017.
View at: Publisher Site | Google Scholar
Y. Sun, P. Chen, B. Zhai et al., “The emerging role of ferroptosis in inflammation,” Biomedicine & Pharmacotherapy, vol. 127, Article ID 110108, 2020.
View at: Publisher Site | Google Scholar
M. A. Bewley, K. B. R. Belchamber, K. K. Chana et al., “Differential effects of p38, MAPK, PI3K or rho kinase inhibitors on bacterial phagocytosis and efferocytosis by macrophages in COPD,” PLoS One, vol. 11, no. 9, Article ID e0163139, 2016.
View at: Publisher Site | Google Scholar
M. Wu, T. Lai, D. Jing et al., “Epithelium-derived IL17A promotes cigarette smoke-induced inflammation and mucus hyperproduction,” American Journal of Respiratory Cell and Molecular Biology, vol. 65, no. 6, pp. 581–592, 2021.
View at: Publisher Site | Google Scholar
N. Amatya, A. V. Garg, and S. L. Gaffen, “IL-17 signaling: the Yin and the Yang,” Trends in Immunology, vol. 38, no. 5, pp. 310–322, 2017.
View at: Publisher Site | Google Scholar
R. M. Tuder, J. H. Yun, A. Bhunia, and I. Fijalkowska, “Hypoxia and chronic lung disease,” Journal of Molecular Medicine (Berlin), vol. 85, no. 12, pp. 1317–1324, 2007.
View at: Publisher Site | Google Scholar
G. H. Fong, “Mechanisms of adaptive angiogenesis to tissue hypoxia,” Angiogenesis, vol. 11, no. 2, pp. 121–140, 2008.
View at: Publisher Site | Google Scholar
A. Zimna and M. Kurpisz, “Hypoxia-inducible factor-1 in physiological and pathophysiological angiogenesis: applications and therapies,” BioMed Research International, vol. 2015, Article ID 549412, 13 pages, 2015.
View at: Publisher Site | Google Scholar
X. Fu and F. Zhang, “Role of the HIF-1 signaling pathway in chronic obstructive pulmonary disease,” Experimental and Therapeutic Medicine, vol. 16, no. 6, pp. 4553–4561, 2018.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Ke 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