Evidence Construction of Chuankezhi Injection Against Chronic Obstructive Pulmonary Disease: A Systematic Review and Network Pharmacology

Abstract

Objective

Chronic obstructive pulmonary disease (COPD) is a chronic respiratory disease with high prevalence, morbidity, and mortality. Chuankezhi (CKZ) injection, a Chinese patent medicine, has been commonly used for treating COPD. This study evaluated the clinical efficacy of CKZ injections in COPD patients and explored potential underlying mechanisms by integrating meta-analysis and network pharmacology.

Research Methods

Randomized controlled trials (RCTs) were search in database by Web of Science, Cochrane Library and PubMed as of November 2022 for literature collection, and the Review Manager 5.4 was used to analyze the data. Through the network pharmacology method, the chemical components and their targets, as well as the disease targets were further analyzed.

Results

A total of 15 RCTs including 1212 patients were included. The results of meta-analysis showed that CKZ injection can significantly improve the clinical effective rate (RR = 1.25, 95% CI: 1.14 to 1.36), and the clinical advantage was that it can significantly reduced acute exacerbation rate (RR = 0.29, 95% CI: 0.12 to 0.70) and COPD assessment test (CAT) scores (MD =−4.62, 95% CI:-8.966 to-0.28). A total of 31 chemical compounds and 178 potential targets for CKZ injection were obtained from the online databases. Molecular docking revealed that most key components and targets could form stable structure.

Conclusion

This systematic review with meta-analysis and network pharmacology demonstrates that CKZ could effectively improve the clinical efficacy and safety in the treatment of COPD. Such efficacy may be related to an anti-inflammatory effect and immunoregulation of CKZ via multiple components, multiple targets and multiple pathways.

Graphical Abstract

Keywords:

• meta-analysis
• network pharmacology
• chronic obstructive pulmonary disease
• Chuankezhi injection

Introduction

Chronic obstructive pulmonary disease (COPD) is a common obstructive pulmonary disease, which can be prevented and treated.¹ It is characterized by persistent respiratory symptoms and airway obstruction induced by abnormalities in the airways alveoli as a result of prolonged exposure to noxious substances, particularly tobacco smoke and pollutants.² More than 300 million individuals worldwide suffer from COPD, which is now a significant global burden and is responsible for over 3.2 million annual deaths.³ This disease will become the third biggest cause of mortality worldwide by 2030, according to the World Health Organization (WHO).⁴ Importantly, reduplicative acute exacerbations can lead to a loss of pulmonary function and quality of life, culminating in death.⁵ Thus, it is necessary to manage and prevent the progression of COPD⁶ since inhaled corticosteroids and bronchodilators, which are the front-line therapies, are oftentimes ineffective.⁷

Due to its relative safety and distinct superiority, Traditional Chinese medicine (TCM) has begun to gain international recognition during the past few years.⁸ TCM is widely used in China to treat COPD patients and is gradually being accepted by the public due to its clinical efficacy and low incidence of side effects.⁹ The Chuankezhi injection (CKZ) is made up of two traditional Chinese herbal medicines, Epimedium sagittatum (Sieb. et Zucc.) Maxim. and Morinda officinalis F.C.How.¹⁰ These medicines are originated from the Shen Nong Ben Cao Jing, in which Epimedium sagittatum (Sieb. et Zucc.) Maxim.is thought to “warm the kidneys” and “help yang to stop coughing”.¹¹ It is said in “Ben Cao Jing Shu” that Morinda officinalis F.C.How. can subsidize Yuan Yang, while “Ben Cao Qi Yuan” states that it can resolve phlegm and stop coughing and wheezing.¹² According to a previous study, Epimedium sagittatum (Sieb. et Zucc.) Maxim. was shown to alleviate lipopolysaccharide (LPS)-induced inflammatory response by inhibiting the NF-κB signaling pathway.¹³ Meanwhile, Epimedium sagittatum (Sieb. et Zucc.) Maxim. can attenuate oxidative stress in human lung epithelial cells via quenching reactive oxygen species (ROS), up-regulating GSH and modulating PI3K/Akt/Nrf2 signaling.¹⁴

Numerous studies have demonstrated that CKZ is able to cure cough and asthma through its anti-allergy, anti-inflammatory, stress response, and immunoregulatory functions.^15,16 CKZ has also been used in the adjuvant treatment of COPD patients with its with its remarkable effects being demonstrated in previous studies in China.^17,18 Nevertheless, systematical assessment of efficacy, safety and underlying mechanism of CKZ against COPD is also necessary to be further processed. Combination of evidenced-based medicine study and network pharmacology have been regarded as an effective approach to better clarify above issues. Hence, we conducted a meta-analysis to estimate the efficacy and safety of CKZ against COPD. Furthermore, network pharmacology and molecular docking were utilized to explore potential pharmacological mechanisms of CKZ against COPD.

Methods

The Preferred Reporting Items for Systematic Reviews and Meta-Analyses checklist was used to guide this study.¹⁹ The systematic review protocol was registered at PROSPERO(CRD:42022370971).

A literature search was conducted by screening four English databases (PubMed, Cochrane Library, Web of science and Embase) and four Chinese databases (China National Knowledge Infrastructure (CNKI), the Wan fang database, the China Biology Medicine disc (SinoMed) database and CQVIP Database (VIP) from their inceptions to November 2022. Additionally, the Chinese clinical trial registration website (http://www.chictr) and the international clinical trial register website (http://clinicaltrials.gov/) were used as resources. We also looked through the publication of meeting minutes, grey literature, and reference lists. English and Chinese were considered the primary search languages. The terms were entered alone or in different combinations and included: “Chronic Obstructive Pulmonary Disease”, “Chronic Obstructive Lung Disease”, “Chronic Obstructive Airway Disease”, “Chuankezhi” and “Chuankezhi Injection”. Details on search terms and search strategy are listed in Supplementary Table 1.

Inclusion Criteria

We incorporated randomized controlled trials (RCTs) that examined the effectiveness and safety of CKZ in patients with stable COPD. Adult patients (>18 years) were diagnosed according to the Global Initiative for Chronic Obstructive Lung Disease (GOLD)⁴ with a postbronchodilator fixed ratio of forced expiratory volume in first second (FEV₁)/forced vital capacity (FVC) < 0.70 measured by spirometry without restrictions in terms of gender, race, educational attainment, or socioeconomic situation. Stable COPD is defined as stable or mild symptoms of cough, expectoration, and dyspnea, with no acute exacerbation within 4 weeks.²

Conventional treatments of control groups included: oxygen inhalation therapy, anti-inflammatories, cough alleviation drugs, sputum reduction drugs, and inhaled bronchodilators. The intervention of the experimental group was a combination of CKZ injection and conventional treatments.

The following outcomes were collected: clinical total effective rate, COPD assessment test (CAT) score, pulmonary function index (FEV₁%pred, FEV₁/FVC), St George’s Respiratory Questionnaire (SGRQ) score, six-minute walk distance (6MWD), Quality of Life Instruments for Chronic Diseases-COPD (QLICD-COPD), modified Medical Research Council (mMRC) dyspnea score, TCM symptom score and adverse events.

Exclusion Criteria

The exclusion criteria included: (1) non-RCTs or uncontrolled clinical trials; (2) studies that used CKZ in conjunction other TCM therapies such as acupuncture, herbal paste applied to acupoints and decoction therapy; (3) patients with major disorders such as diffuse bronchiolitis, tuberculosis, congestive heart failure, asthma, bronchiectasis, or congestive heart failure; (4) the trial’s specifics were unclear or the data it produced were insufficient.

Study Selection

Two review authors (Xuan Wei and Yu Zhong) identified possibly suitable studies after individually examining titles and abstracts that were returned by the literature search. The same review writers then obtained the full text of each study and independently evaluated them in accordance with the inclusion and exclusion criteria. In cases of disagreement, a third reviewer (Xiaofei Yi) served as the arbitrator.

Data Extraction

Two review authors (Xuan Wei and Xiaofei Yi) extracted data and evaluated the literature quality for each qualifying study. Then, the above two reviewers checked the data extracted by the other party to ensure accuracy. A third researcher (Yu Zhong) was invited to make a final decision in cases of disagreement.

The following information and data were extracted according to a standard data extraction form that was created prior to analyzing the studies: title, first author name, year of publication, country, study design, participant characteristics, treatment duration, outcomes, and adverse events. To retrieve relevant information that could not be found in the published report, we e-mailed the corresponding author.

Risk of Bias Assessment

The Cochrane Risk of Bias Assessment Tool, Version 5.1.0 was adopted by two reviewers to assess the methodological quality of included studies.²⁰ The items used in this analysis were as follows: (1) random sequence generation (selection bias); (2) allocation concealment (selection bias); (3) blinding of participants and personnel (performance bias); (4) blinding of outcome assessment (detection bias); (5) incomplete outcome data (attrition bias); (6) selective reporting (reporting bias), and other biases. Each type of bias was ranked as: high risk, unclear risk or low risk. Reviewers settled disagreements by dialogue and, when that was impracticable, through adjudication by a third party.

Statistical Analysis

The meta-analysis was conducted by the Review Manager (RevMan) 5.4 software. To appraise dichotomous and continuous variables outcomes, Relative risk (RR) and its 95% confidence intervals (CIs), mean difference (MD) and its 95% CIs were respectively used. Data that could not be included in a meta-analysis are demonstrated in the final table I² was used to assess statistical heterogeneity among included studies. A fixed-effect model was chosen for when I²<50%; otherwise, a random-effect model was applied. We also performed a subgroup analysis to identify potential sources of heterogeneity. In order to perform the subgroup analysis, characteristics and differences between included studies were examined including intervention and follow-up times, including the treatment philosophies associated with TCM syndromes, disease severity, location and age of participants.

This study used secondary data from published RCTs, in which there was no personal information disclaimed. Hence, it was not necessary to obtain ethical approval and consent from participants of included studies.

CKZ Active Ingredients and Targets

The Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP; https://old.tcmsp-e.com/index.php) was used to identify active components of CKZ. The Canonical Simplified molecular input line entry system (Canonical SMILES) of each ingredient was searched at PubChem (https://pubchem.ncbi.nlm.nih.gov/). After that, the Swiss Target Prediction database (http://swisstargetprediction.ch/) was applied to forecast the therapeutic targets of CKZ and the UniProt database (https://www.uniprot.org/) was utilized to acquire canonical gene names of each target. CKZ-related targets were obtained after removing duplicates.

Screening for COPD-Related Targets

GeneCards (https://www.genecards.org/), Online Mendelian Inheritance in Man (OMIM; https://omim.org/), Drugbank (https://www.drugbank.com/), DisGeNET (https://www.disgenet.org/), and Target Database (TTD; http://db.idrblab.net/ttd/) were used to search COPD-related targets with “Chronic Obstructive Pulmonary Disease” as the keyword.

Network Topology Analysis of Protein-Protein Interaction

Potential targets of CKZ relevant to COPD were defined as the common intersection of drug and disease targets. Owing to the multi-component and multi-target characteristic of traditional Chinese herbal prescription, a protein-protein interaction (PPI) network was constructed to demonstrate the complex therapeutic mechanisms of CKZ. First, intersection target data were imported into the STRING database (https://string-db.org/) to build the PPI network.²¹ The conditions were set as “Homo sapiens” and the confidence score was defined as >0.4. Then, the PPI network was visualized in Cytoscape 3.9.0. The CytoHubba plugin was used to obtain hub genes. A node’s size corresponds to its degree value in the network, and there is a positive correlation between the two.

Target Function and Pathway Enrichment Analysis

The Database for Annotation, Visualization, and Integrated Discovery (DAVID) 6.8 (https://david.ncifcrf.gov/) was used to analyze the enrichment of GO and KEGG pathways to investigate the biological pathways of key genes. An adjusted P value<0.05, was deemed statistically significant.

Molecular Docking

To confirm the binding force of hub genes and key ingredients from the PPI network, we performed molecular docking. Specifically, we used CB-Dock, which is an online website providing an online protein-ligand docking approach to ascertain binding sites, evaluate center and size, and demonstrate molecular docking.

The Protein Data Bank database (PDB; https://www.rcsb.org) was utilized to obtain the 3D structure of target proteins, while the chemical structures of active components were acquired from PubChem (https://pubchem.ncbi.nlm.nih.gov). The Vina score represents the binding energy and the stability of the interconnection module has a positive correlation with the absolute Vina score value. Finally, PyMOL software was used to visually represent molecular docking results.

Results

Study Characteristics

A total of 15 RCTs published in Chinese were included. The study flowchart is outlined in Figure 1. These 15 trials included a total of 1, 212 participants and were published between 2009 and 2022. Further details are presented in Table 1.

Table 1 Characteristics of Included Studies

Study (Year)	Sample Size			Average Age (Year)		Male/Female		Treatment Duration (Week)	Interventions		Outcomes
	C	E	Total	C	E	C	E		C	E
Huang 2009^Citation22	29	27	56	68.52±8.48	66.41±9.95	2/27	1/26	12	Bronchodilators, corticosteroids	CT+CKZ 4mL im, qod	1,2,4,6,7,8,10
Jiang 2015^Citation23	57	59	116	57.80±5.30	58.20±5.20	47/10	46/13	4	Budesonide, formoterol inhalant	CT+CKZ 4mL im, bid	1,2,5,6,10,
Li 2020^Citation24	25	25	50	57.20±13.80	56.90±14.10	13/12	14/11	12	Theophylline, sustained-release tablets	CT+CKZ 4mL acupoint injection, biw	3,7,10,11,12,14,15
Li 2022^Citation25	39	39	78	62.14±9.74	62.75±9.83	23/16	22/17	2	Budesonide, oxygen	CT+CKZ 4mL iv, qd	1,2,3,10,12
Xiao 2020^Citation26	44	44	88	68.71±4.38	68.03±4.51	27/17	31/13	4	Budesonide aerosol	CT+CKZ 4mL im, bid	6,10,11
Pan 2016^Citation17	58	60	118	72.49± 5.69	74. 05±5.28	41/17	41/19	12	Bronchodilators, corticosteroids, oxygen	CT+CKZ 2mL acupoint injection, qd	1,9,10,11,16
Shan 2019^Citation18	80	80	160	58. 40±4.20	59. 20±3. 90	46/34	42/38	12	Oxygen, tiotropium bromide	CT+CKZ 4mL acupoint injection, biw	2,3,4,6,7,10 12,13,14,15
Wang 2014^Citation27	30	30	60	66.80±3.80	67.60±6.30	17/13	16/14	2	Salmeterol xinafoate and fluticasone	CT+CKZ 4mL im, bid	1,2,10,14,15
Wang 2017^Citation28	30	30	60	72.50±8.50	71.50±8.50	20/10	18/12	2	Salmeterol xinafoate and fluticasone	CT+CKZ 4mL iv, bid	7,8,10
Xu 2015^Citation29	58	60	118	72.49±5.69	74.05±5.28	41/17	41/19	12	Bronchodilators, corticosteroids	CT+CKZ 2mL acupoint injection, qd	1,9,10,11,16
Yu 2010^Citation30	30	30	60	83.00 ±4.50	82.00±5.30	23/7	22/8	4	Oxygen, tiotropium, Bromide	CT+CKZ 4mL im, qd	2,3,10,11,13
Yuan 2021^Citation31	49	49	98	66.12±4.15	66.75±4.03	27/22	28/21	12	Oxygen, tiotropium, Bromide	CT+CKZ 4mL im, bid	10,11
Zhan 2016^Citation32	15	15	30	56.27±8.47	58.40±7.71	7/8	9/6	12	Conventional therapy	CT+CKZ 4mL im, qd	2,3,4,7,10,11
Zhang 2021^Citation33	32	32	64	64.91±7.47	68.63±8.93	9/23	12/20	12	Conventional therapy	CT+CKZ 2mL acupoint injection, qd	1,2,3,4,7
Zheng 2010^Citation34	24	24	48	59	60	15/9	16/8	4	Conventional therapy	CT+CKZ 2mL aerosol inhalation, bid	3,10,11

Notes: Outcomes: 1. FEV₁/FVC; 2. FEV₁%pre; 3. clinical effective rates; 4. CAT scores/mMRc scores; 5. acute exacerbations; 6. 6MWD; 7. TCM syndrome scores; 8. SGRQ scores; 9. QLICD-COPD scores; 10. adverse events; 11. FEV₁(L)/FVC(L); 12. immunoglobulin; 13. T cell; 14. PaO₂; 15. PaCO₂; 16. Barthel scores.

Abbreviations: E, experiment group; C, control group; CT, conventional treatment group.

Figure 1 Flow chart of meta-analysis.

Bias Assessment

The Cochrane Risk of Bias Assessment Tool was used to estimate the quality of included studies. Four studies reported the use of randomization (eg using a random number table). Hence, these RCTs were judged as having a “low risk” for selection bias. The remainder were unclear in regard to their randomization techniques and were therefore classified as having an “unclear risk” for selection bias. None of the studies mentioned missing values as a result of withdrawals or dropouts. In all studies, there was a low risk of bias in selective result reporting. None of the research made reference to study design, sample size estimation, or randomization scheme concealment. Figure 2 depicts that it was unclear whether there were other biases in the RCTs herein analyzed.

Figure 2 Risk bias of included studies.

Notes: (A) Overall risk; (B) Detailed risk.

Efficacy Assessment

Clinical Effective Rate

This meta-analysis included six studies that provided categorical data for clinical effective rates. According to the data analysis, the heterogeneity was not statistically significant (x² = 3.92, P = 0.56; I² = 0%); hence, a fixed-effect model was chosen. The analysis results revealed that the clinical effective rate of the experimental group was higher than that of the control group (RR = 1.25; 95% CI: 1.14 to 1.36; Z = 4.91, P<0.00001) (Figure 3).

Figure 3 Forest plot of clinical effective rate.

Acute Exacerbation Rate

Two trials that provided categorical data for the acute exacerbation rate were included. Since there was no heterogeneity (x² =0.00, P = 1.00; I² = 0%), a fixed-effect model was used. The data analysis demonstrated that the experimental group had a significantly reduced acute exacerbation rate compared with the control group (RR = 0.29; 95% CI: 0.12 to 0.70; Z = 2.77, P = 0.006) (Figure 4).

Figure 4 Forest plot of acute exacerbation rate.

Pulmonary Function

Eight studies reported FEV₁%pred as a measure of pulmonary function. Again, there was no heterogeneity (x² = 8.34, P = 0.30; I² = 16%), hence, a fixed-effect model was chosen. The combined data demonstrated that the improvement of FEV₁%pre in the experiment group was superior to the improvement in the control group (MD = 7.35;95% CI:5.87 to 8.83; Z = 9.72, P<0.00001) (Figure 5). Meanwhile, eight studies provided numerical data for FEV₁/FVC and a random effects model was utilized due to considerable heterogeneity. In a subgroup analysis investigating different administration routes, one study found no statistically significant change in FEV₁/FVC between groups that used aerosol inhalation (MD = 3.62; 95% CI: −2.66 to 9.90; Z = 1.13, P = 0.96) (Figure 6), while the outcomes for studies using intramuscular injections of CKZ combined with conventional treatment were comparable in the subgroup analysis (MD = 1.55; 95% CI:-1.57 to 4.67; Z = 0.97, P = 0.33) (Figure 6). One study showed that a combination of intravenous injection of CKZ and conventional treatment resulted in better FEV₁/FVC values than the control group (MD = 8.35; 95% CI: 5.09 to 11.61; Z = 5.02, P<0.00001) (Figure 6), similar to when CKZ was administered through acupoint injection (MD = 7.95; 95% CI: 6.70 to 9.20; Z = 12.42, P<0.00001) (Figure 6).

Figure 5 Forest plot of FEV₁% pre.

Figure 6 Forest plot of FEV₁/FVC.

CAT Score

Three of the fifteen included studies reported variations in CAT scores. The pooled analysis revealed that the experimental group displayed a considerable elevation in their CAT score (MD =−4.62; 95% CI:-8.966 to-0.28; Z = 2.08, P = 0.04) compared to control patients. However, there was significant heterogeneity between the groups (x² = 22.36, P<0.0001; I² = 91%). To pinpoint the source of variability, a subgroup analysis of different CKZ administration routes was carried out. The results demonstrated that CAT scores decreased in both studies that used acupoint injections (MD =−5.34; 95% CI:-10.73 to 0.05; Z = 1.94, P<0.0001) and intramuscular injections of CKZ (MD =−2.74; 95% CI:-6.98 to 1.50; Z = 1.27, P = 0.20). However, heterogeneity was still present in this subgroup analysis (Figure 7A).

Figure 7 Forest plot of quality of life score. (A) CAT score. (B) QLICD-COPD score. (C) SGRQ score.

SGRQ Score

With 118 participants, two of the fifteen included studies mentioned changes in the SGRQ score. Indeed, CKZ combined with conventional treatment improved the SGRQ score of COPD patients to a greater extent than what was observed with conventional treatment alone (MD =−6.42; 95% CI:-10.71 to-2.13; Z = 2.93, P = 0.003). There was no statistical heterogeneity between studies (x² = 0.00, P = 1.00; I² = 0%) (Figure 7B).

QLICD-COPD

The QLICD-COPD score was reported in 2 of 15 studies, comprising 232 participants. The pooled analysis demonstrated that when CKZ was added to conventional treatment, the QLICD-COPD score significantly improved compared to conventional treatment alone (MD = 6.07; 95% CI: 3.79 to 8.35; Z = 5.22, P<0.0001). In addition, there was no statistical heterogeneity among groups (x² = 0.00, P = 1.00; I² = 0%) (Figure 7C).

6MWD

Four trials with 420 participants reported 6MWD data. Data analysis indicated that there was significant heterogeneity (x² = 13.21, P = 0.004; I² = 77%); hence, a random-effect model was adopted. The results revealed a significant improvement in 6MWD for patients in the experimental group compared to that of the control group (MD = 32.17; 95% CI: 17.76 to 46.95; Z = 4.37, P < 0.00001) (Figure 8).

Figure 8 Forest plot of 6MWD.

TCM Syndrome Score

Five trials with a total of 390 participants reported data for the TCM Syndrome Score. We identified significant heterogeneity (x² = 26.90, P < 0.0001; I² = 85%); hence, a random-effect model was chosen. The results demonstrated that the TCM syndrome score of the experimental group was much lower than that of the control group, suggesting a lower occurrence of severe clinical symptoms in patients treated with CKZ (MD =−3.39; 95% CI:-5.06 to-1.72; Z = 3.98, P < 0.0001) (Figure 9).

Figure 9 Forest plot of TCM Syndrome Score.

Other Outcomes

In addition, two studies compared immunoglobulin M, immunoglobulin A and immunoglobulin G levels and one study compared BODE scores between groups (Table 2). It was clear that CKZ treatment could improve circulating levels of immunoglobulins and reduce the BODE score compared with outcomes of the control group.

Table 2 Other Outcomes

Study (Year)	N(C/E)	Outcome Indicator	C		E		P
			Before	After	Before	After	P1	P2
Li 2020^Citation24	25/25	IgM	0.2±0.1	0.3±0.1	0.2±0.1	0.6±0.2	>0.05	<0.05
		IgA	0.5±0.1	0.8±0.2	0.5±0.2	1.1±0.4	>0.05	<0.05
		IgG	7.2±0.5	7.7±0.9	7.2±0.6	8.8±0.7	>0.05	<0.05
Li 2022^Citation25	39/39	IgM	0.86±0.08	1.34±0.12	0.87±0.08	1.65±0.15	0.583	<0.05
		IgA	1.38±0.13	2.34±0.21	1.40±0.13	3.17±0.30	0.499	<0.05
		IgG	6.41±0.63	8.92±0.88	6.37±0.62	11.52±1.13	0.778	<0.05
Zhang 2021^Citation33	32/32	BODE score	3.72±2.89	3.72±1.67	3.81±2.39	2.53±2.02	0.888	0.008

Notes: P1: p-value of the comparison between the control and the experimental groups before intervention; P2: p-value of the comparison between control and experimental groups after intervention.

Abbreviations: E, experiment group; C, control group;

Adverse Reactions

Four studies reported adverse reactions, such as palpitation, cough, hoarseness, and dizziness, among others. All adverse reactions were temporary and had no long-term negative repercussions. There were no adverse reactions or incidents reported by other studies.

Active Ingredients and Targets Screening of CKZ

According to pharmacokinetic parameters (drug-likeness (DL) ≥0.18, drug half-life (HL) ≥4, and relative molecular weight (MW)≤500), 31 active ingredients of CKZ were identified at the TCMSP database (Table 3). Through PubChem and Swiss Target Prediction databases, related targets of the 31 active ingredients were also identified. A total of 217 CKZ-related targets were obtained after eliminating redundant ones and standardizing them through UniProt.

Table 3 List of CKZ Ingredients

Mol ID	Molecule Name	MW	DL	HL
MOL000006	Luteolin	286.25	0.25	15.94
MOL000098	Quercetin	302.25	0.28	14.40
MOL000359	Sitosterol	414.79	0.75	5.37
MOL000422	Kaempferol	286.25	0.24	14.74
MOL001510	24-epicampesterol	400.76	0.71	4.5
MOL001645	Linoleyl acetate	308.56	0.2	7.48
MOL001771	Poriferast-5-en-3beta-ol	414.79	0.75	5.07
MOL001792	DFV	256.27	0.18	17.89
MOL003044	Chryseriol	300.28	0.27	16.31
MOL003542	8-Isopentenyl-kaempferol	354.38	0.39	15.37
MOL004373	Anhydroicaritin	368.41	0.44	15.01
MOL004380	C-Homoerythrinan,1,6-didehydro-3,15,16-trimethoxy- (3.beta)	329.48	0.49	6.58
MOL004382	Yinyanghuo A	420.49	0.77	14.44
MOL004384	Yinyanghuo C	336.36	0.5	15.74
MOL004386	Yinyanghuo E	352.36	0.55	15.47
MOL004391	8-(3-methylbut-2-enyl)-2-phenyl-chromone	290.38	0.25	18.73
MOL002883	Ethyl oleate (NF)	310.58	0.19	4.85
MOL000358	Beta-sitosterol	414.79	0.75	5.36
MOL000359	Sitosterol	414.79	0.75	5.37
MOL006147	Alizarin-2-methylether	254.25	0.21	30.82
MOL009495	2-hydroxy-1,5-dimethoxy-6-(methoxymethyl)-9,10-anthraquinone	328.34	0.37	29.98
MOL009496	1,5,7-trihydroxy-6-methoxy-2-methoxymethylanthracenequinone	330.31	0.38	28.25
MOL009500	1,6-dihydroxy-5-methoxy-2-(methoxymethyl)-9,10-anthraquinone	314.31	0.34	29.02
MOL009503	1-hydroxy-3-methoxy-9,10-anthraquinone	254.25	0.21	30.78
MOL009504	1-hydroxy-6-hydroxymethylanthracenequinone	254.25	0.21	30.8
MOL009513	2-hydroxy-1,8-dimethoxy-7-methoxymethylanthracenequinone	328.34	0.37	29.44
MOL009519	(2R,3S) - (+)-3’,5-Dihydroxy-4,7-dimethoxydihydroflavonol	332.33	0.33	14.13
MOL009524	3beta,20(R),5-alkenyl-stigmastol	414.79	0.75	5.22
MOL009525	3beta-24S(R)-butyl-5-alkenyl-cholestol	456.88	0.82	7.04
MOL009558	2-hydroxyethyl-5-hydroxy-2-(2-hydroxybenzoyl)-4-(hydroxymethyl)benzoate	332.33	0.26	17.61
MOL009562	Ohioensin-A	372.39	0.76	15.09

Abbreviations: MW, molecular weight; DL, drug-likeness; HL, drug half-life.

COPD-Related Target Gene Collection

COPD-related genes were collected from GeneCards, OMIM, DisGeNET, DrugBank and TTD databases. The score in the Genecards database has a positive correlation with how prevalent the condition is. The maximum score of the target obtained by GeneCards was 142.24, while the minimum score was 0.28, and the median was 8.11. Hence, the target with a score higher than the median was classified as a target associated with COPD. Overall, 3215 targets were identified from GeneCards, 544 targets from OMIM, 73 targets from TTD, 83 targets from DrugBank and 180 targets from DisGeNET. Finally, a total of 3778 targets associated with COPD were collected after eliminating duplicate genes.

Screening of Hub Genes and PPI Network Construction

We discovered that COPD and CKZ shared 178 targets by comparing the two sets of genes (Figure 10A). The common targets were imported into STRING online database to establish the PPI network. Hub genes were chosen using the CytoHubba plugin and the Maximum Cross-Correlation (MCC) computation method. The top 10 hub genes include: JUN, Casp3, VEGFA, AKT1, MMP9, PTGS2, TP53, EGF, NF-κB1A and HIF1A (Figure 10B). An active ingredient-target network was created in Cytoscape, in which green represents the targets, and orange represents the active ingredients of CKZ (Figure 11). The importance of each protein and component in the network was examined in accordance with the node degree, and a network diagram was created for visualization. The top five most significant elements in CKZ, as determined by the Analyze Network plugin and based on node degree, were: quercetin, luteolin, kaempferol, Anhydroicaritin and β-sitosterol.

Figure 10 Common targets and PPI network of hub genes.

Notes: (A) Common targets between CKZ and COPD; (B)PPI network of 10 hub genes.

Figure 11 Active ingredient-target network.

Notes: Green oval represent targets, orange rhombus represent active ingredients of CKZ.

GO Enrichment Analysis

To systematically identify the treatment pathways of CKZ in COPD, a GO enrichment analysis was carried out. We identified 110 items for cellular components (CC), 172 items for molecular functions (MF), and 836 items for biological processes (BP). Separately, we chose the top 10 BP, CC, and MF data to depict in a bubble diagram (Figure 12A).

Figure 12 GO and KEGG enrichment bubble diagram.

Notes: (A) GO enrichment bubble diagram; (B) KEGG enrichment bubble diagram.

KEGG Pathway Enrichment Analysis

The top 20 pathways were chosen for presentation and visualization out of a total of 181 pathway items that were enriched in the KEGG enrichment analysis. These items primarily included the AGE-RAGE signaling pathway, Cellular senescence, and PI3K/Akt signaling pathway. Based on the results of the KEGG enrichment, the first 20 associated enrichment results were visualized via a bubble diagram (Figure 12B).

Molecular Docking

According to the ingredient-target network diagram, the top 6 core ingredients were selected for molecular docking with their respective hub genes relevant to COPD. The core target protein ID codes were obtained through the PDB database for JUN (6Y3V), Casp3 (1RHJ), VEGFA (5DN2), AKT1 (4EJN), MMP9 (5I12), and PTGS2 (5F19). The results of docking scores of quercetin (MOL000098), luteolin (MOL000006), kaempferol (MOL000422), Anhydroicaritin (MOL004373) and β-sitosterol (MOL000358) with key targets are listed in Figure 13A. The binding energies of all ingredients were all lower than −5 kcal/mol, indicating a stable interaction. The top 5 interconnection modules with the lowest vina score were chosen for visual demonstration (Figure 13B).

Figure 13 Molecular docking results.

Notes: (A) Heat map of molecule docking affinity; (B) Molecular docking visualization results.

Discussion

COPD is a complex pulmonary inflammatory disease, with airway and parenchymal lung damage, which is now the third leading cause of death in the world and is estimated to cause the greatest economic burden of disease in China.^35,36 It is generally considered that the intake of pollutants, particularly cigarette smoke, is the primary cause of COPD.³⁷

Reiterative airway inflammation and airflow limitation can trigger an imbalance of immune function, which is strongly associated with the continuous development of COPD.³⁸ The TCM theory states that COPD falls within the categories of “cough”, “phlegm” and “lung distension”.³⁹ Lung-kidney Qi deficiency is the most common syndrome in stable COPD.⁴⁰ Hence, the treatment of COPD with TCM focuses on nourishing the kidney and warming yang.⁴¹ The two herbs in CKZ exert these properties. Therefore, this research first evaluated the efficacy and safety of CKZ for COPD using an evidence-based analytical approach. Our systematic review included 15 RCTs with 1212 individuals and suggested that CKZ was beneficial to COPD patients since it alleviated the risk of acute exacerbation, improved quality of life, pulmonary function, clinical effectiveness and activity endurance.

Exacerbation is a recurring issue in the management of COPD since it might impair pulmonary function if it occurs frequently.³⁶ Clinical effective rate is widely considered to be a subjective outcome measure in TCM treatment, which reflects symptom change.⁴² The SGRQ is recommended as the golden standard to evaluate the symptoms of COPD patients, given that it assesses the psychological state of COPD patients.⁴³ It also has an intimate correlation with FEV₁%pre, FEV₁/FVC, and clinical symptoms.⁴⁴ Meanwhile, the CAT score is a multi-dimensional assessment of a patient’s health, based on a comprehensive symptom score, and the mMRC is mainly used to evaluate the degree of dyspnea in COPD patients.⁴⁵ The QLICD-COPD, which includes modules on social, psychological, and physical functions, was created to assess the health condition of individuals with COPD.⁴⁶ Another significant feature of COPD is lower exercise tolerance, which may be caused by muscle atrophy and impaired cardiovascular function. Hence, the 6MWD test measures this decline in activity tolerance.⁴⁷ Based on subgroup analyses, we suggest that CKZ combined with conventional treatment improves the: SGRQ score, CAT score, 6MWD, QLICD-COPD score and clinical effective rate when compared to participants using conventional treatment alone. Both FEV₁/FVC and FEV₁%pre play a significant role in diagnostic standards as well as in evaluating disease progression.⁴⁸ Furthermore, our analyses suggest that acupoint injection of CKZ leads to a significant elevation in FEV₁/FVC compared to other administration routes. Therefore, the effectiveness of CKZ may depend on the route of administration.

Only four of the included research documented precise randomization methods, such as random number tables, while the remaining studies just referenced randomization without further details. Double-blinding was not stated in any of the included RCTs, and its absence may result in measurement bias. Only four of the 15 studies that were included reported adverse events. Even though the adverse events were tolerated by the patients, did not interfere with the course of the treatment, or disappeared following drug withdrawal, a personalized treatment plan must be created to account for individual variances. Due to the incompleteness of raw and comprehensive data, the baseline characteristics of different study participants (such as disease severity and progression) may also vary. This could be one of the reasons for the high heterogeneity of some analyses. While conducting the meta-analysis, we attempted to avert and prevent compromises due to high heterogeneity. Nevertheless, a high-quality RCT with large sample size is necessary to corroborate our findings, owing to the small sample size, multiple confounding variables, and poor quality of studies herein analyzed. Based on the aforementioned research findings, it is clear that CKZ can benefit COPD patients by improving clinical efficacy. In addition, the pooled analysis indicated that serum IgM, IgA and IgG were increased in the experimental group, suggesting that the therapeutic effectiveness of CKZ for COPD might be related to an enhancement in host immunity. These findings lay the groundwork for additional research into the pharmacological mechanisms of action of CKZ.

Therefore, we explored the pharmacological mechanisms of CKZ in COPD by network pharmacology and molecule docking. A total of 217 targets related to CKZ and 31 active components were found and PPI networks were established by integrating 178 intersecting targets of CKZ and COPD. The top 6 ingredients of the drug-active ingredient-target network were: quercetin, luteolin, kaempferol, Anhydroicaritin and β-sitosterol. It was found that JUN, Casp3, VEGFA, AKT1, MMP9 and PTGS2 were the most critical potential targets of CKZ relevant to COPD.

Quercetin is a natural compound found in a variety of plants that exerts anti-inflammatory, antioxidant, anti-cancer effects.⁴⁹ It was found that quercetin can reduce inflammatory storms and oxidative stress induced by LPS.⁵⁰ Besides, quercetin could inhibit oxidative stress, reduce inflammatory factors and the expression of MMP9, MMP12 and MUC5AC to decelerate emphysema in elastase-induced mice.⁵¹ Luteolin is found in a variety of plants and it is a flavonoid with anti-inflammatory, antibacterial, antioxidant, and anti-cancer properties.⁵² Indeed, luteolin considerably reduces the production of superoxide anion, and Raf1 activity to alleviate inflammation storms and oxidative stress induced by neutrophils in inflammatory arthritis.⁵³ Interestingly, luteolin can regulate the levels of Treg and Tregs-derived IL-10, and alleviate caspase-11-dependent pyroptosis to inhibit inflammation responses in the CLP-induced acute lung injury mouse model.⁵⁴ Kaempferol is an antimicrobial, anti-inflammatory, antioxidant, and antitumor flavonoid aglycone widely distributed in numerous plants.⁵⁵ Kaempferol was shown to mitigate allergic response and airway inflammation induced by OVA via the NF-κB signal pathway and by inhibiting the expression of otaxin-1 in eosinophils.⁵⁵ Anhydroicaritin, a natural component originated from Epimedii Herba, was shown to inhibit cancer cell multiplication, cell differentiation, invasion, and metastasis.⁵⁶ In addition, Anhydroicaritin can regulate the LKB1/AMPK/mTOR pathway to ameliorate diet-induced obesity and alleviate insulin resistance.⁵⁷ β-sitosterol is a primary component from common plants and vegetables that exerts anti-hyperlipidemia, anti-inflammation and anti-tumor effects.⁵⁸

Analysis of the disease target-ingredient network and PPI network suggested that JUN, Casp3, VEGFA, AKT1, MMP9, PTGS2, TP53 and EGF may be the primary targets of CKZ in COPD. The GO BP functional enrichment analysis revealed that the co-targets of CKZ and COPD were mainly enriched in the positive regulation of transcription from RNA polymerase II promoter, positive regulation of transcription, positive regulation of gene expression, response to drug, and other functions. The GO CC functional analysis mostly focused on the extracellular space, cytosol, nucleoplasm, plasma membrane, and other processes, while the GO MF functional analysis was primarily enriched in protein binding, enzyme binding, protein homodimerization activity, protein kinase binding, and other functions. The primary KEGG pathway included Pathways in cancer, lipid and atherosclerosis, AGE-RAGE signaling pathway in diabetic complications, PI3K/Akt signaling pathway, chemical carcinogenesis-receptor activation, human cytomegalovirus infection, IL-17 signaling pathway, and fluid shear stress and atherosclerosis. These functions and pathways are associated with immunity and inflammatory response during the onset and progression of COPD. Hence, these data are indicative that CKZ mediates airway and pulmonary inflammation via regulating inflammatory factors, which is crucial for the treatment of COPD. Finally, molecular docking was conducted to further explore the mechanism of CKZ in COPD. Our data demonstrate that most ingredients could stably bind to their target proteins. PTGS2, a cyclooxygenase, inhibits cellular and humoral immune responses in endometriosis via participating in the angiogenesis and regulation of T cell differentiation.⁵⁹ It has been be proved that Luteolin could suppress protein expression such as VEGF, MMP-2 and MMP-9, which are related to cell migration and invasion, and regulate PI3K/Akt and MAPK signaling pathways to prevent tumorigenesis of non-small cell lung cancer.⁶⁰ Quercetin was shown to accelerate the release of CASP8 and decrease the expression of EGFR, VEGFA, CCND1, and ERBB2 in cervical cancer cells, which improved the anticancer effects of cisplatin.⁶¹ Thus, the above molecular docking results reiterate that the active ingredients of CKZ are likely to bind to proteins implicated in COPD.

This study integrated meta-analysis and network pharmacology to evaluate the clinical efficacy of CKZ in COPD and explore potential mechanisms. Here we show how TCM can be used to treat COPD. Nevertheless, there were still some limitations. First, since CKZ is a TCM, all research included were published in Chinese and had limited sample sizes. More research should be conducted outside China to validate our findings. Second, a large-scale, multicenter, randomized controlled trials are needed to validate our data. Finally, since this study is primarily based on literature and database research, experiments are required to verify identified mechanisms.

Conclusion

In summary, based on meta-analysis, network pharmacology, and molecular docking, this study revealed that CKZ is an effective adjuvant therapy for COPD, and its mechanisms may be related to reducing inflammatory and responses.

Abbreviations

COPD, chronic obstructive pulmonary disease; FEV₁, forced expiratory volume in one second; FVC, forced vital capacity; HE, hematoxylin and eosin; CAT, COPD Assessment Test; mMRC, modified Medical Research Council scale; SGRQ, St George’s Respiratory Questionnaire; 6MWD, six-minute walk distance; QLICD-COPD, Quality of Life Instruments for Chronic Diseases-COPD; BC, Betweenness centrality; CC, Closeness centrality; DL, drug-like; DAVID, Database for Annotation, Visualization, and Integrated Discovery; TCM, traditional Chinese medicine; TCMSP, Traditional Chinese Medicine Systems Pharmacology; OB, oral bioavailability; TTD, Therapeutic Target Database; OMIM, Online Mendelian Inheritance in Man database; PPI, protein‒protein interaction; DC, Degree centrality; SD, standard deviation; BP, biological processes; MF, molecular functions; CC, cellular components; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; PDB, Protein Data Bank.

Ethics Approval and Consent to Participate

Ethical approval was not necessary for the current study because our study gathered data only from Embase, PubMed, the Cochrane Library, the China National Knowledge Infrastructure Database (CNKI), the Wanfang Database, the CQVIP Database (VIP), the China Biology Medicine disc (SinoMed) database, the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP), PubChem, UniProt database, Therapeutic Target Database (TTD), Online Mendelian Inheritance in Man (OMIM), DisGeNET, and STRING database, and this procedure did not address any patient’s personal data or harm any patient.

Disclosure

The authors report no conflicts of interest in this work.

Additional information

Funding

This study was supported by the National Natural Science Foundation of China (Grant No. 81860010) and Guangxi Medical High-level Backbone Personnel Training “139” Plan Training Project (Grant No. G201903008, G202002002).