{'fraction_non_alphanumeric': 0.033709520582255355, 'fraction_numerical': 0.016177033915785132, 'mean_word_length': 5.392540968167263, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 5, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Systems pharmacology-based approach for dissecting the active ingredients and potential targets of the Chinese herbal Bufei Jianpi formula for the treatment of COPD

Background The Chinese herbal Bufei Jianpi formula (BJF) provides an effective treatment option for chronic obstructive pulmonary disease (COPD). However, the systems-level mechanism underlying the clinical effects of BJF on COPD remains unknown. Methods In this study, a systems pharmacology model based on absorption filtering, network targeting, and systems analyses was applied specifically to clarify the active compounds and therapeutic mechanisms of BJF. Then, a rat model of cigarette smoke- and bacterial infection-induced COPD was used to investigate the therapeutic mechanisms of BJF on COPD and its comorbidity. Results The pharmacological system successfully identified 145 bioactive ingredients from BJF and revealed 175 potential targets. There was a significant target overlap between the herbal constituents of BJF. These results suggested that each herb of BJF connected with similar multitargets, indicating potential synergistic effects among them. The integrated target–disease network showed that BJF probably was efficient for the treatment of not only respiratory tract diseases but also other diseases, such as nervous system and cardiovascular diseases. The possible mechanisms of action of BJF were related to activation of inflammatory response, immune responses, and matrix metalloproteinases, among others. Furthermore, we demonstrated that BJF treatment could effectively prevent COPD and its comorbidities, such as ventricular hypertrophy, by inhibition of inflammatory cytokine production, matrix metalloproteinases expression, and other cytokine production in vivo. Conclusion This study using the systems pharmacology method, in combination with in vivo experiments, helped us successfully dissect the molecular mechanism of BJF for the treatment of COPD and predict the potential targets of the multicomponent BJF, which provides a new approach to illustrate the synergetic mechanism of the complex prescription and discover more effective drugs against COPD.


Introduction
Chronic obstructive pulmonary disease (COPD) is a serious health problem characterized primarily by irreversible airflow limitation due to an enhanced chronic inflammatory response in the airways and lung to noxious particles or gases. 1,2 Cigarette smoking is considered to be the major risk factor for COPD. Apart from smoking, bacterial infections accentuate airway inflammation, which accelerates irreversible disease progression. 3,4 To prevent the progression of COPD, millions of patients are treated by drug therapies that include bronchodilators, steroids, and phosphodiesterase inhibitors. These drugs offer therapeutic benefit but cause serious side effects and cannot effectively reduce the progression and mortality of this disease. For instance, aminophylline is used as a bronchodilator in COPD. The commonest side effects are headache, nausea and vomiting, abdominal discomfort, and restlessness. Convulsions, cardiac arrhythmias, and death may occur at high concentrations. 5 Hence, there is a pressing need for the development of more effective therapies for COPD. [6][7][8] Traditional Chinese medicine (TCM) is a comprehensive medicinal system that has provided effective relief of symptoms in COPD patients for thousands of years. Bufei Jianpi formula (BJF), a TCM, is composed of 12 medicinal herbs and has been proven effective for the treatment of COPD. Clinical studies revealed a variety of desirable pharmacological effects of BJF on COPD, such as alleviating the clinical symptoms of stable COPD patients, reducing the exacerbation frequency, delaying acute exacerbation, and improving pulmonary function and exercise capacity. 9 However, the molecular details regarding the effects of BJF on COPD are still unclear.
Studies on the pharmacological effects of TCM, such as BJF, face various difficulties for three main reasons. Firstly, identification of the complex active compounds contained in medicinal herbs is laborious and time-consuming. Secondly, there may be highly dynamic interactions among the active components, which also is one of the possible mechanisms of its treatment. Lastly, the multiple ingredients contained in most formulas and a lack of effective methods of study cause difficulties in identifying the targets of the herbal medicines. 10,11 Thus, a novel method that can identify the active compounds and the targets of herbal drugs, as well as clarify the mechanisms of action of Chinese medicine, is necessary.
Systems pharmacology, a new emerging research field of pharmacology, has provided promising avenues to help clarify the active compounds and therapeutic mechanisms of Chinese medicine. [11][12][13][14] This technological platform combines pharmacokinetics evaluation and drug-target-disease network analysis to get a global view of the multiple effects and mechanisms of Chinese herbal medicine. 15,16 Furthermore, the term "systems pharmacology" presently describes another promising new field of study that incorporates computational and experimental methods to explain the effects and molecular mechanisms of herbal medicines. In this work, we constructed a systems pharmacological model by integrating oral bioavailability predictions, target predictions, and network analysis to provide novel insights into the effect and mechanism of BJF. Finally, in order to supply the in vivo experimental evidence to validate the mechanisms of BJF action that were predicted in the computational experiment (systems pharmacology), we administered BJF to COPD model rats and examined the effect of BJF on cigarette smoke-and bacterial infection-induced pulmonary inflammatory responses, protease-antiprotease imbalance, and respiratory dysfunction.

Materials and methods
This trial was registered at the Chinese Clinical Trial Register Center: ChiCTR-TRC-11001406.

Dataset construction
All ingredients from the 12 herbs of BJF were collected mainly from the Chinese Academy of Sciences Chemistry Database (http://www.organchem.csdb.cn), Chinese Herbal Drug Database, and the literature. 17-20 For orally administered drugs, glucosides can be metabolized extensively to their deglycosylation products by enteric bacteria in the intestinal tract; 21 thus, both glucosides and deglycosylation products are considered to be the constituents of herbal drugs. Taken together, a total of 886 chemicals were included: 87 in

Drug-likeness prediction
The drug-likeness index was used to evaluate the structural similarity between the herbal ingredients and the drugs in the DrugBank database (http://www.drugbank.ca/) and help remove compounds that are considered to be chemically and pharmacologically unsuitable as drugs. 24 In this study, the database-dependent drug-likeness prediction approach was calculated as follows: where A represents the herbal compounds and B represents the average molecular drug-likeness index of all compounds in the DrugBank database. A drug-likeness index $0.18 (average value for all DrugBank molecules) was set as the threshold to select drug-like compounds. The compounds that overcame both the OB and drug-likeness screens were retained as candidate compounds. In addition, several compounds such as tangshenoside II, atractylenolide I, atractylenolide III, ergosterol, naringin, hesperidin, bergenin, icariside I, and anhydroicaritin initially were omitted according to these screening rules; however, these compounds were supported by literature evidence and, therefore, also were obtained as candidate compounds for further analysis. [21][22][23][24][25][26][27][28][29] Drug targeting analysis The targets related to the candidate compounds were predicted by the systematic drug targeting tool, 30 which efficiently integrates the chemical, genomic, and pharmacological information for drug targeting by RandomForest and Support Vector Machine methods. This model shows a good predictive performance for drug-target interactions, and the results are represented by a probability of interactions between each compound and a number of targets from the DrugBank database. In this work, the compound-target interactions with a Support Vector Machine score $0.8 and RandomForest score $0.7 were considered as potential targets of the candidate compounds. Then, the diseases related to the potential targets were collected from the Therapeutic Target, DrugBank, and PharmGkb databases, and the obtained diseases were classified further into different groups according to the MeSH Browser (2014 MeSH).

network construction
The herbs, candidate compounds, potential targets, and diseases related with the potential targets were used to build the compound-target and target-disease networks. These networks were generated by Cytoscape 3.2.1. 31 In the graphical networks, nodes encode the herbs, compounds, target proteins, or diseases, and edges represent the herbcompound, herb-target, or target-disease interactions. The degree of a node is the number of edges associated with it. Nodes with a high degree can be considered the key nodes in a network. The topological properties of the networks were analyzed using the Network Analysis plugin of Cytoscape. 32

COPD model and drug administration
Rats were placed into a chamber connected to a cigarette smoke-producing apparatus (volume 300 L) and exposed to the tobacco (Hongqi Canal ® Filter tip cigarette; tobacco type, tar: 10 mg; nicotine content: 1.0 mg; carbon monoxide: 12 mg, Henan Tobacco Industry, Zhengzhou, People's Republic of China) smoke of eight cigarettes for 30 minutes, twice per day with 3-hour smoke-free intervals during the first 2 weeks. Then, COPD rats were exposed to the smoke of 15 cigarettes for 30 minutes, thrice per day with 3-hour smoke-free intervals from weeks 3 to 12. Additionally, 100 μL of Klebsiella pneumonia suspension (6×10 8 CFU/ mL) was dropped slowly into the nasal cavities of smokeexposed rats every 5 days from weeks 1 to 8. On week 9, two COPD rats were sacrificed to obtain the lung tissues to confirm that this rat model was successful. 33 International Journal of COPD 2015:10 submit your manuscript | www.dovepress.com

2636
Zhao et al COPD rats were divided into three groups with 10 rats each (week 9). Then, COPD rats were intragastrically administrated normal saline (2 mL), BJF (4.84 g/kg), and aminophylline (2.3 mg/kg) every day for 12 weeks. The control rats also were administrated intragastrically normal saline (2 mL) every day for 12 weeks. All rats were sacrificed, and the heart and lung tissues were collected at week 20. The components of BJF were as follows: AR 15 g, PR 15 g, CR 15 g, AMR 12 g, Po 12 g, FTB 9 g, Pheretima 12 g, MOC 9 g, CRP 9 g, ATR 9 g, AJH 15 g, and EH 6 g. The herbs were identified and prepared in fluid extract according to the standard operating procedure of the Department of Pharmaceutics of the First Affiliated Hospital, Henan University of Traditional Chinese Medicine, Zhengzhou, People's Republic of China.

respiratory data collection
Respiratory data were collected by unrestrained pulmonary function testing plethysmography (Buxco Inc., Wilmington, NC, USA) conducted every 4th week from weeks 0 to 20. Buxco air flow transducers were connected to the chambers and to a reference chamber to compensate for pressure changes. Each chamber was calibrated to its respective transducer. We focused on three measures: tidal volume (TV), peak expiratory flow (PEF), and 50% TV expiratory flow.

histological analyses
Heart and lung tissue samples were fixed in 10% formalin neutral buffer solution for 24 hours and then embedded in paraffin before being sliced into 4 μm sections. For histological examination, sections were stained with Mayer's hematoxylin and then with 1% eosin alcohol solution (hematoxylin and eosin staining). Samples were inspected with the aid of an Olympus BX51 microscope (Tokyo, Japan).
Bronchia, lung injury, and bronchiole stenosis were observed under an optical microscope. Alveolar number, alveolar diameter, small pulmonary vessels, and bronchial wall thickness were measured using Image-Pro Plus ® 6.0 software (Media Cybernetics, Rockville, MD, USA). The morphometric analysis at the light microscopic level was conducted by an investigator blinded to the study protocol.
For immunohistochemical analysis, the right ventricle and lung tissue sections were transferred to gelatin-coated slides and then deparaffinized and rehydrated. The tissue slices were washed with 0.3% Triton X-100 in phosphate buffer (PB) and quenched with endogenous peroxidase (3% hydrogen peroxide). The sections were incubated with phosphate-buffered saline containing normal goat serum to reduce nonspecific reactions. Then, tissue sections were incubated overnight at 4°C with 1:300 dilutions of primary antibodies (Santa Cruz Biotechnology) against ET-1, TGF-β, VEGF, bFGF, MMP-2, MMP-9, TIMP-1, IL-6, IL-10, TNF-α, soluble TNF-α receptor 2, and collagens I, III, and IV in the presence of 2.5% bovine serum albumin. Slices were washed with phosphate buffer and incubated with the secondary antibody (Biocare Medical, Concord, CA, USA) for 2 hours. Immune complexes were visualized using the Catalyzed Signal Amplification System (Dako Denmark A/S, Glostrup, Denmark). Sections also were counter-stained with hematoxylin. Samples were inspected and measured using Image-Pro Plus ® 6.0 software (Media Cybernetics).

right ventricular hypertrophy index
After removing the arterial and adipose tissue on the epicardium, the right ventricle (RV), left ventricle (LV), and interventricular septum (S) were separated and weighed. The right ventricular hypertrophy index (RVHI) was calculated using the equation: RVHI = RV/(LV + S). 34 (2)

Myocardial ultrastructure
The right ventricle was fixed with 3% glutaraldehyde, then 1% osmium tetroxide, and dehydrated. Through epoxy resin embedding and deployment of a hardener, accelerator, and growth agent, ultrathin sections of 70 nm thickness were cut and stained with uranyl acetate. Changes in the muscular fibers and mitochondria of the myocardial ultrastructure were observed with transmission electron microscopy.

2637
Bufei Jianpi formula for the treatment of COPD USA) experiments. Specificity was confirmed by electrophoretic analysis of the reaction products and inclusion of template-or reverse transcriptase-free controls. Glyceraldehyde-3-phosphate dehydrogenase and hypoxanthine guanine phosphoribosyl transferase cDNA were used as an internal standard.

statistical analysis
Values are expressed as mean ± standard error of mean. Statistical differences between groups were analyzed by one-way analysis of variance with the SPSS 19.0 software package (IBM Corporation, Armonk, NY, USA).

Results and discussion
Candidate compounds screening In most cases, oral administration is the most predominant way to deliver drugs to systemic circulation for patients in TCM therapy; however, most compounds in herbs fail to reach to their cellular targets because they lack suitable pharmacokinetic properties, such as high oral bioavailability. 35 In addition, drug-like compounds contain functional groups and have physical properties that are consistent with clinical drugs. 20 Thus, a combination of OB screening and druglikeness property evaluation was used to determine whether a compound was pharmaceutically active in the complex herbal mixture. Pheretima, an herbal constituent of BJF, was excluded because all of its compounds showed poor OB or drug-likeness indices. Then, 131 potential ingredients with OB $30% and drug-likeness index $0.18 were obtained from the eleven herbs. Additionally, we found that 14 molecules with low OB (#30%) or drug-likeness indices (#0.18) were the most abundant constituents and possessed extensive pharmacological activities. [21][22][23][24][25][26][27][28][29] Thus, these 14 ingredients also were included in the candidate active compounds for further target prediction and network analysis. Finally, 145 compounds of the eleven herbs were regarded as "active compounds", including the 131 readily absorbed compounds and the 14 pharmacologically active compounds (  27,36,37 Similarly, luteolin and quercetin, which occur widely in herbs such as AR, ATR, AJH, and EH, have been shown to overcome OB barriers by inhibition of permeability glycoprotein and possess a broad spectrum of biological and pharmacological properties. 38,39 Target identification In general, TCM formula contains tens of thousands of active compounds, which provide bright prospects for prevention and treatment of complex diseases in a synergistic manner. To illuminate the detailed mechanisms of the effects of TCM formula, it is crucial to uncover the therapeutic targets of active compounds contained in the formulas. In this study, a pharmacophore modeling technique was applied to search potential targets based on the "candidate compounds". A total of 175 candidate targets were identified for the 131 candidate compounds; however, 14 candidate compounds had no candidate target under this criterion (Table S2). Finally, 131 candidate compounds yielded 175 potential targets, and the connections between them reached 1,849. The numbers of potential targets connected by AR, PR, CR, AMR, Po, FTB, MOC, CRP, ATR, AJH, and EH were 134, 70, 92, 33, 36, 56, 43, 56, 121, 99, and 140, respectively. Although the numbers of each herb-related target were different, there was a significant target overlap between the eleven herbs. These results indicated that each herb could connect with similar multitargets, suggesting the potential synergistic effects among them. For example, tangshenoside II (from CR), atractylenolide III (from AMR), honokiol (from MOC), and nobiletin (from CRP) have been shown to attenuate the function of nitric oxide synthase, which plays an important role in regulating blood pressure and, thus, ameliorate coronary artery vascular tone and local or systemic inflammatory disorders. 22,23,37,40 Compound-target network TCM formulas are regarded as multicomponent therapeutics, which means that the active compounds of the herbal drugs simultaneously interact with the same or different targets. To understand the drug action at a systems level, the network theory was used to uncover the underlying interactions between the candidate compounds and targets (Table S3). Here, we mapped the compounds from the eleven herbs and targets onto compound-target networks. Figure 1 shows a global view of the compound-target network with colorcoded nodes: herb (square, orange), candidate compounds (rhombus, blue), and candidate targets (ellipse, green). The network consisted of 1,849 compound-target interactions connecting the 131 candidate compounds to 175 targets, resulting in an average number of targets per active compound   1. This suggests that many compounds bind to multiple targets, thereby giving rise to polypharmacology. These compounds play important roles in the networks, which were selected to illustrate the pharmacological activities of BJF, such as the two highly connected nodes, quercetin (degree =87) and kaempferol (degree =55). For instance, quercetin was one of the largest hubs connecting to different targets. Quercetin had high affinities with peroxisome proliferator activated receptor gamma (PPARγ), nitric oxide synthase, inducible, prostaglandin G/H synthase 1/2 (also known as COX-1, -2), and estrogen receptor, which induce a complex pharmacological profile including anti-inflammatory, immunomodulatory, antioxidant, and antiproliferative effects, among others. Furthermore, these results also indicated that many target are cross-linked together in the compound-target network. For example, the androgen receptor, the most highly connected target, is linked to 79% (104 of 131) of the compounds, followed by the estrogen receptor (degree =93), prostaglandin G/H synthase 1/2 (degree =66), nitric oxide synthase, inducible (degree =56), dipeptidyl peptidase IV (degree =56), thrombin (degree =54), and PPARγ (degree =53). Interestingly, nitric oxide synthase, inducible, PPARγ, and prostaglandin G/H synthase 1/2 expression is correlated closely with inflammation and immune system activation. 41,42 In an effort to illuminate the therapeutic mechanism of the 175 selected targets, we applied ClueGO, a Cytoscape plugin, to explore the biological interpretation of the potential targets, 43 and divided the results into three stratums: molecular function, the immune system processes, and the reactome analysis. As shown in Figure 2A, the main molecular functions were classified into four categories: oxidoreductase activity, mitogen-activated protein kinase (MAPK) activity, neurotransmitter receptor activity, and G-protein-coupled amine receptor activity. The large majority of targets were related to oxidoreductase and MAPK activities, which are the activities mainly responsible for COPD. 44,45 In Figure 2B, the immune system processes of the targets are associated mainly with positive regulation of myeloid leukocyte differentiation, T-cell lineage commitment, positive regulation of immune effector processes, and the toll-like receptor 10 signaling pathway. Consequently, the reactome of the targets are related to activation of the activator protein-1 family of transcription factors, MAPK targets/nuclear events mediated by MAPK, and activation of MMPs ( Figure 2C).
Multiple intermediate risk factors, such as chronic bronchitis, small-airway disease, and emphysema, 3,6,46 which might induce differences in protein expression, are incriminated in COPD etiology. Thus, a single compound or herb is probably insufficient for COPD therapy. From our results, common targets with multiple functions shared by the active compounds suggest that the BJF could possess synergistic therapeutic effects for COPD and possibly could be more effective than a single herb or compound.

Target-disease network
All of the 175 targets were projected into the DrugBank, Therapeutic Target, and PharmGkb databases to collect their related diseases. Then, 348 diseases belonging to 16 groups derived from the MeSH Browser (2014 MeSH) were obtained. Finally, the target-disease interactions are shown in the target-disease network (Figure 3). Interestingly, among these 348 diseases, we found that most of them belong to neoplasms (87/348), nervous system diseases (65/348), cardiovascular diseases (50/348), and nutritional and metabolic diseases (21/348), apart from respiratory tract diseases (26/348). This illustrated that BJF probably has effectiveness not only on respiratory diseases but also on these diseases. This is consistent with the use of Chinese medicines for the treatment of various diseases. For instance, PPARγ connected with inflammatory response and immunomodulatory activities in atherosclerosis, type 2 diabetes, inflammatory bowel disease, arthritis, asthma, COPD, and other airway diseases and, as a result, has the potential to become a therapeutic target for COPD and other diseases. [47][48][49] In the target-disease network, we also found that PPARγ connected with many of these diseases. In addition, to clarify further the relationship among the compounds, the targets, and the correlated diseases, experimental studies were conducted.

effect of BJF on COPD and its comorbidity
In a previous study, BJF was confirmed to have beneficial clinical efficacy in COPD patients. 9 In addition, the targetdisease network predicts that BJF has efficiency not only for the treatment of respiratory diseases but also for other diseases, such as nervous system and cardiovascular diseases. To test further this prediction, we next examined the effect of BJF on COPD rats and its comorbidity, ventricular hypertrophy. [50][51][52] To determine the effect of BJF on COPD rats, the intragastric administration of BJF to rats was commenced 9 weeks after the administration of cigarette smoke, and lung mechanics and pulmonary histopathology were assessed. As shown in Figure 4, compared with the control rat, the TV, PEF, and 50% TV expiratory flow clearly decreased in the model rat from weeks 4 to 20, whereas, compared with the model rats, BJF had significantly increased the TV and PEF in the COPD rats at week 20. Furthermore, aminophylline, a classical bronchodilator, also increased the TV and PEF in COPD rats ( Figure 4A and B). Similarly, BJF and aminophylline also could increase the 50% TV expiratory flow at week 20 ( Figure 4C). In addition, lung injury scores, bronchiole wall thickness, small pulmonary vessels wall thickness, bronchiole stenosis, and alveolar diameter increased in the model rat, and this increase was clearly suppressed by the treatment of animals with BJF ( Figure 5A-F). Moreover, BJF treatment also markedly increased the alveolar number in COPD rats ( Figure 5G). These results suggested that BJF is an effective Chinese medicine for the treatment of COPD.
Many different studies have shown that comorbidities, such as cardiovascular, metabolic, muscular, and bone disorders, occur more frequently in patients with COPD than in smokers and never-smoked control subjects of similar age. 50,51 In this study, we also found that rats with COPD are at an increased risk of ventricular hypertrophy. Therefore, we examined the effects of BJF on ventricular hypertrophy, which is a common comorbidity of COPD. As shown in Figure 6, compared with the model rat, BJF significantly increased the myocardial sarcomere length and mitochondrial density of cardiocytes and reversed the RVHI at week 20. Furthermore, studies have demonstrated that hypertrophic stimuli, such as ET-1, TGF-β, VEGF, bFGF, and angiotensin II, can activate a range of hypertrophic signaling mediators and transcription factors, such as extracellular signal-regulated kinase 1/2 and nuclear factor-κB, and enhance myocardial angiogenesis. 53,54 In this work, we found that BJF clearly decreased the expression of VEGF, bFGF, TGF-β, and ET-1. Furthermore, the effects of BJF were higher than those of aminophylline (Figure 7).
These findings suggested that BJF treatment can effectively prevent COPD and ventricular hypertrophy, and these results are consistent with those obtained by the targetdisease network.

Effect of BJF on inflammatory responses in COPD rats
In the systems pharmacology study, we found that numerous prediction targets of BJF are related to activation of the activator protein-1 family of transcription factors and MAPK targets/nuclear events mediated by MAPK, which play important roles in inflammatory responses, such as inflammatory cytokines production. Many studies have shown that airway inflammation of COPD patients, both acute and chronic, is induced by cytokines, such as ILs (IL-1β, IL-4, IL-6, IL-10, and IL-13) and TNF-α, and other mediators, such as the cysteinyl leukotrienes and isoprostanes. 52 To verify this prediction, we examined the levels of inflammatory cytokines in the lung of COPD rats. As shown in Figure 8, levels of IL-6, IL-10, TNF-α, and soluble TNF-α receptor 2 increased after the cigarette smoke and bacterial infection exposures, and this increase was clearly suppressed by the treatment of animals with BJF. These results demonstrated that BJF achieves its anti-inflammatory activity by decreasing the levels of inflammatory cytokines in the lung, which also is consistent with the results obtained by target protein function analysis.

effect of BJF on collagen degradation and protease-antiprotease imbalance
A protease-antiprotease imbalance leads to the breakdown of connective tissue components and is the critical mechanism in the pathogenesis of COPD. Furthermore, there is increasing interest in the role of MMPs in COPD. 55 Increased concentrations of MMP-2 and MMP-9 were observed in bronchoalveolar lavage fluid from COPD patients, and there is increased activity of MMP-9 and decreased activity of TIMP-1 in the lung parenchyma. 56,57 In COPD, abnormally high collagen remodeling also occurs within the lung tissue.
In particular, the turnover rate of type I, III, and IV collagen is changed significantly, leading to excessive remodeling and accumulation of structural proteins. 58 In this work, the systems pharmacology study suggested that the potential targets were related to activation of MMPs. Therefore, we examined the effect of BJF on the expression of MMP-2, MMP-9, TIMP-1, and collagens I, III, and IV in lung tissues. As shown in Figure 9A, BJF treatment clearly decreased the levels of MMP-2 and MMP-9 and increased the level of TIMP-1. This alteration also was observed at the messenger RNA level ( Figure 9B). As shown in Figure 10, collagens I, III, and IV were increased significantly by the cigarette smoke and bacterial infection exposures, and this increase was markedly suppressed by the treatment of animals with BJF.

2649
Bufei Jianpi formula for the treatment of COPD These results suggested that BJF suppressed the cigarette smoke-and bacterial infection-induced collagen deposition and protease-antiprotease imbalance by inhibiting the expression of collagens I, III, and IV, MMP-2/-9, and increasing the expression of TIMP-1. From these results, systems pharmacology under the systems biology framework has been used successfully to investigate the active compounds of BJF and their targets and understand the multiple uses and molecular mechanisms of BJF. Furthermore, we experimentally validated that BJF suppressed cigarette smoke-and bacterial infection-induced pulmonary inflammation, collagen deposition, proteaseantiprotease imbalance, and hypertrophic factor production, suggesting this formula could be effective for the treatment of COPD and its comorbidity.

Conclusion
TCM, an important part of natural medicine, greatly improves the efficiency for the treatment of complex diseases, especially COPD. Even though TCM has accumulated much clinical experience and provided many effective applications, it always has been discredited as a mystery due to a lack of systems methods for identifying the bioactive ingredients and their therapeutic targets. In this work, we proposed a systems pharmacology approach to dissect the molecular mechanisms of the action of BJF by systematically incorporating OB prediction, target proteins prediction, and the target-related disease network analysis. In addition, a rat model of cigarette smoke-and bacterial infection-induced COPD was applied to validate the effects and mechanism of BJF on COPD and its comorbidity. The main findings are as follows. 1) The candidate ingredients and their potential targets in the BJF prescriptions were identified, which provided insights into the synergetic effects of BJF and clues to clarify the mechanisms of BJF for the prevention of COPD and other diseases. 2) The target-disease network showed that BJF has efficiency not only for the treatment of COPD but also for other diseases, such as neoplasms, nervous system diseases, cardiovascular disease, and musculoskeletal diseases, suggesting that different diseases can be treated with the same formula. 3) The experimental study found that BJF ameliorated COPD and its comorbidity in vivo via mechanisms that are dependent on its effects on inflammatory cytokine production, MMP expression, collagen deposition, Bufei Jianpi formula for the treatment of COPD  and other cytokine production, and these results were consistent with the predictions of the systems pharmacology studies. 4) The integrated systems pharmacology platform constructed in this work illuminated the effects and mechanisms of BJF on COPD and its comorbidity, which provides valuable insights to support the development of novel drugs.
In summary, systems pharmacology currently describes an innovative research field of study that incorporates experimental and computational approaches to explain both active ingredients and targets of TCM, which could help us to explore the complex TCM mechanism-of-action profiles and promote TCM modernization in the future study. The major limitation of this study is its direct demand for further experimental validation such as the proportion of every component, which will be resolved in the follow-up studies.
Tiao-Bu Fei-Shen Therapies in Rats with Chronic Obstructive Pulmonary Disease on Regulation of Multidimensional Molecular Network, 81130062).

Author contributions
JSL and YHW designed the outline of the study, and contributed toward data analysis and critically revising the paper. PZ performed experiments, conceived the study, wrote the draft and revised the manuscript. YGT and YL were involved in performing experiments, acquisition of data, and statistical analysis. CLZ contributed to the data analysis and interpretation. All authors contributed toward data analysis, drafting and critically revising the paper, gave final approval of the version to be published, and agree to be accountable for all aspects of the work.

Disclosure
The authors report no conflicts of interest in this work.

Publish your work in this journal
Submit your manuscript here: http://www.dovepress.com/international-journal-of-chronic-obstructive-pulmonary-disease-journal The International Journal of COPD is an international, peer-reviewed journal of therapeutics and pharmacology focusing on concise rapid reporting of clinical studies and reviews in COPD. Special focus is given to the pathophysiological processes underlying the disease, intervention programs, patient focused education, and self management protocols.
This journal is indexed on PubMed Central, MedLine and CAS. The manuscript management system is completely online and includes a very quick and fair peer-review system, which is all easy to use. Visit http://www.dovepress.com/testimonials.php to read real quotes from published authors.