Systems Pharmacology-Based Approach to Comparatively Study the Independent and Synergistic Mechanisms of Danhong Injection and Naoxintong Capsule in Ischemic Stroke Treatment

To provide evidence for the better clinical use of traditional Chinese medicine preparations (TCMPs), comparison of the pharmacological mechanisms between TCMPs with similar therapeutic effect is necessary. However, methodology for dealing with this issue is still scarce. Danhong injection (DHI) and Naoxintong capsule (NXT) are representative TCMPs for ischemic stroke (IS) treatment, which are also frequently used in combination. Here they were employed as research objects to demonstrate the feasibility of systems pharmacology approach in elucidation of the independent and combined effect of TCMPs. By incorporating chemical screening, target prediction, and network construction, a feasible systems pharmacology model has been established to systematically uncover the underlying action mechanisms of DHI, NXT, or their pair in IS treatment. Systematic analysis of the created TCMP-Compound-Target-Disease network revealed that DHI and NXT shared common targets such as PTGS2, F2, ADRB1, IL6, ALDH2, and CCL2, which were involved in the vasomotor system regulation, blood-brain barrier disruption, redox imbalance, neurotrophin activity, and brain inflammation. In comparative mechanism study, the merged DHI/NXT-IS PPI network and pathway enrichment analysis indicated that DHI and NXT exerted the therapeutic effects mainly through immune system and VEGF signaling pathways. Meanwhile, they had their own unique pathways, e.g., calcium signaling pathway for DHI and gap junction for NXT. While for their synergistic mechanism, DHI and NXT participated in chemokine signaling pathway, T cell receptor signaling pathway, VEGF signaling pathway, gap junction, and so on. Our study provided an optimized strategy for dissecting the different and combined effect of TCMPs with similar actions.


Introduction
Ischemic stroke (IS) is a common disease caused by intracranial or neck blood vessel occlusion that impairs blood flow to a portion of the brain and leads to subsequent cerebral necrosis [1]. Due to its high morbidity and mortality, as well as high-risk of permanent disability in surviving individuals, IS represents a substantial economic and social burden throughout the world [2,3]. Currently, pharmacological or mechanical thrombolytic therapy for restoration of blood flow to ischemic tissue is the main treatment to improve patient outcome [4,5]. However, their unique benefits are shadowed by limited therapeutic time window and serious complications [6]. Moreover, cellular consequences of ischemic brain injury are very complex and chemical drugs are hardly to achieve favorable outcomes because of their single therapeutic target [7][8][9]. Thus, novel therapeutic drugs and strategies are urgently needed to improve the efficiency in treating IS.
Traditional Chinese medicine (TCM), characterized as multiple compounds and targets, differs in substance, methodology, and philosophy to chemical drugs [10]. It has been practiced in China for thousands of years and is now becoming more frequently used in the West. In recent years, TCM preparations (TCMPs) have been proposed as important complementary and alternative medicines for IS treatment [11,12]. Danhong injection (DHI) and Naoxintong 2 Evidence-Based Complementary and Alternative Medicine capsule (NXT) are both commonly used TCMPs for the prevention and treatment of cardiovascular and cerebrovascular diseases [13][14][15][16]. DHI is extracted from Radix Salviae Miltiorrhizae and Flos Carthami, while NXT consists of 16 kinds of herbs including Radix Astragali, Radix Paeoniae Rubra, Radix Salviae Miltiorrhizae, and so on. Their efficiency in protecting against IS has been widely validated, whether by used alone or in combination [17][18][19][20][21]. According to a recent report, DHI-NXT combination therapy could achieve better therapeutic effect than monotherapy in cerebral ischemiareperfusion model of rats [22]. However, the comparative and synergistic effect of the action mechanisms between DHI and NXT still remain vague and warrant further investigation.
As the multicomponent and multitarget nature of TCMPs, it is very difficult to systematically study their mechanisms using routine methods [23]. Therefore, systems pharmacology method is now widely perceived as an integral and efficient tool to study the complex molecular mechanisms of herbal medicines [24][25][26][27][28][29][30]. This approach usually includes three steps, including active compounds screening, drug targets prediction, and network/pathway analysis [31,32]. It could provide global active compounds-therapeutic targets interactions based on the network pharmacology framework, which matches well with the integrity and systemic conception of TCM. For example, a network of 31 active compounds interacting with 42 disease targets was constructed by Yue et al. to elucidate the polypharmacological mechanisms underlying the efficiency of herb pair Danggui-Honghua for blood stasis syndrome treatment [33]. Still, methodology for comparative study of the mechanisms between TCMPs is scarce.
To provide evidence for the better clinical use of DHI and NXT, here we proposed a feasible systems pharmacologybased method to explore the similarities and differences of the pharmacological mechanisms between DHI and NXT. Their synergistic mechanism in treatment of IS was also investigated. To our knowledge, no work has been reported about the mechanism comparison and synergistic effect of DHI and NXT.

Results and Discussion
In this study, a network pharmacology approach based on chemical, pharmacokinetic, and pharmacological data was employed to systematically investigate the comparative and synergistic effects of DHI and NXT against IS. The overall procedure involved five steps ( Figure 1): (1) the chemical ingredients of DHI and NXT were collected and evaluated by in silico ADME (absorption, distribution, metabolism, and excretion) system to screen candidate compounds; (2) the putative targets of the candidate compounds and IS were data-mined from various databases and literatures; (3) a basal TCMP-Compound-Target-Disease network was constructed to illustrate the interaction among DHI/NXT, ingredients, targets and IS; (4) for comparative study of the action mechanisms, two protein-protein interaction (PPI) networks named DHI-IS and NXT-IS PPI network were merged and pathway enrichment analysis was separately performed to explore the common and unique mechanisms of DHI and NXT; (5) for the study of their synergistic effect, DHI-NXT-IS PPI network was established among DHIingredient targets, NXT-ingredient targets, and IS targets. The enriched pathways were analyzed to explain the synergistic mechanisms of the drug combination.
. . Candidate Ingredients in DHI and NXT. In most systems pharmacology studies, the components in TCMPs were mechanically collected from chemical databases by simple summation of the compounds in each herb. As a matter of fact, the chemical substances and composition ratios of TCMPs differ dramatically due to different preparation processes. For example, tanshinones, as the main components from ethanol extract of Radix Salviae Miltiorrhizae, could not be experimentally detected in DHI because of its water extraction process [34,35]. Hence, the chemical ingredients in DHI and NXT were collected from the high-resolution LC-MS profile data reported by related literatures to guarantee the accuracy and objectivity of this work [34][35][36][37][38]. Since the glycosides in NXT might be deglycosylated by the glycosidase in the intestinal tract, 14 aglycones were also considered as chemical ingredients labeled by qt. As a result, a total of 293 ingredients were retrieved for DHI (65) and NXT (228). The detailed information of these ingredients was displayed in Table S1.
Although TCMPs usually contain tens or hundreds of chemical ingredients, only a few of them with favorable pharmacodynamic and pharmacokinetic properties are responsible for their therapeutic effects. Drug-likeness (DL) is an important drug research parameter that used to assess whether a compound functions as a drug. Drug-like compounds are those which 'contain functional groups and/or have physical properties consistent with the majority of known drugs' [39]. Therefore, an optimal DL model based on molecules in the DrugBank database (https://www.drugbank.ca/) was employed to prescreen potential active ingredients in DHI and NXT. In this study, DL ≥ 0.18 was used to select candidate compounds for further research. Since NXT is consumed orally, oral bioavailability (OB ≥ 30%) and Caco-2 cell permeability (Caco-2 > -0.4) were also used to screen the candidate ingredients from NXT [31][32][33]. Further, a few compounds that do not meet these criteria but possessed high amounts and extensive pharmacological activities were also adopted. Consequently, 35 and 82 candidate ingredients were selected from DHI and NXT, respectively. Among them, 10 ingredients were common for both DHI and NXT, such as chlorogenic acid, hydroxysafflor yellow A, ferulic acid, rosmarinic acid, and salvianolic acid B (Table 1).          text-mining method [41,42]. The detailed information of all the targets was displayed in Table S2. A total of 179 targets were found to be associated with IS. Among them, 71 targets could be targeted by the candidate ingredients of DHI and NXT. Of note, 37 targets were shared by both of them; 40 and 68 targets were regarded as direct anti-IS targets of DHI and NXT, respectively (Figure 2(b)). These results showed that the targets overlapped dramatically in these two TCMPs, which also indicated different ingredients in DHI and NXT shared common or similar targets with synergistic effects.

. . Target Collection and Analysis of TCMP-Compound
For the purpose of interpreting the complex relationships among TCMPs, candidate ingredients, putative targets, and IS, a TCMP-Compound-Target-Disease network (Figure 2(a)) was established using Cytoscape (version 3.6.0). This network consisted of 693 nodes (2 TCMPs, 107 candidate ingredients and 584 targets) and 2561 compound-target interactions. The IS targets targeted by the ingredients of DHI and NXT were labeled by yellow circle. The degree of a node, defined the number of edges connected to it, was analyzed by Network Analysis plugin. The average degree of candidate ingredients was 25.03 and 29 ingredients possessed degrees greater than 25, indicating that these ingredients were more likely to exert therapeutic effects by acting on multiple targets. Among them, the one with a higher degree was more important in the network [43], especially quercetin (degree = 291), apigenin (degree = 143), kaempferol (degree = 128), luteolin (degree = 126), rutin (degree = 116), genistein (degree = 116), caffeic acid (degree = 108), rosmarinic acid (degree = 90), ferulic acid (degree = 76), and biochanin A (degree = 67). For example, in focal cerebral ischemia rats, quercetin can decrease cell apoptosis, prevent free radicals associated with oxidative damage, and reduce the elevated MMP-9 activity, suggesting that quercetin might have a potential role in the treatment of IS [44][45][46]. Rosmarinic acid, a representative compound in both DHI and NXT, showed protective effects against IS induced brain injury by synaptogenic activity and anti-inflammatory action [47,48].
As shown in the TCMP-Compound-Target-Disease network (Figure 2(a)), the ingredients of DHI and NXT could not only modulate the crucial targets involved in the regulation of the vasomotor system and the disruption of bloodbrain barrier (PTGS2, F2, ADRB1, and EGFR), but also could concentrate on the redox imbalance, neurotrophin activity, and brain inflammation (IL6, ALDH2, CCL2, STAT3, and AKT1). These processes were closely associated with IS pathogenesis. For example, increased F2 generation was associated with the risk of IS in the elderly [49,50]. The protein F2 was found to be connected with 33 ingredients, among which danshensu, salvianolic acid C and kaempferol were reported to inhibit the activity of F2 [51][52][53]. In addition, the neurotoxic proinflammatory mediator PTGS2 also had interactions with multiple ingredients. Specifically, ferulic acid, tanshinone IIA, and (Z)-ligustilide could exert antiinflammatory effects by attenuating the PTGS2 expression [54][55][56].
. . Comparative Study of the Anti-IS Mechanisms of DHI and NXT. The interaction between putative targets and IS targets were derived by mapping them onto a PPI network. In this section, the IS targets and ingredient targets from DHI or NXT were integrated into STRING database (https://string-db.org/) to acquire PPI data, respectively. The species was limited to "Homo sapiens" and the confidence score was set as 0.7. Irrelevant proteins, including the putative targets and IS targets that could not directly or indirectly interact with each other, were removed. Two PPI networks, named DHI-IS and NXT-IS PPI network, were generated by Cytoscape and then merged to identify common and unique targets of DHI and NXT (Figure 3). A hub node should be two times greater than the average node degree in the network [57]; 47 and 75 hub nodes were identified for DHI-IS and NXT-IS PPI network, respectively. In the DHI-IS PPI network, 29 hub nodes were IS targets and 29 were DHI-ingredient targets. Besides the 11 overlapping targets, the remaining 18 DHI-ingredient targets were regarded as indirect anti-IS targets of DHI, illustrating the important indirect interactions with the IS targets. While in NXT-IS PPI network, 23 out of 75 hub nodes were IS targets, 67 were NXT-ingredient targets, and 15 were overlapping targets. From them, 52 targets were considered as indirect anti-IS targets of NXT.
To comparatively dissect the underlying mechanisms of DHI and NXT, 58 putative targets of DHI (combining 40 direct anti-IS targets with 18 indirect anti-IS targets) and 120 putative targets of NXT (combining 68 direct anti-IS targets with 52 indirect anti-IS targets) were employed to conduct pathway enrichment analysis, respectively. The pathway information was achieved from Database for Annotation, Visualization and Integrated Discovery (DAVID, https://david.ncifcrf.gov/). The top 20 pathways (Benjamini Hochberg corrected p value < 0.05) for DHI and NXT were shown in Figure 4. DHI and NXT shared 11 common pathways, including 4 immune system pathways (T cell receptor signaling pathway, toll-like receptor signaling pathway, chemokine signaling pathway, and NOD-like receptor signaling pathway), 2 signal transduction pathways (ErbB signaling pathway and VEGF signaling pathway), and 5 cancer-related pathways (colorectal cancer, pancreatic cancer, bladder cancer, prostate cancer, and pathways in cancer). Among them, T cell receptor signaling pathway, toll-like receptor signaling pathway, chemokine signaling pathway, NOD-like receptor signaling pathway, and VEGF signaling pathway were closely associated with cardiocerebrovascular diseases [58]. For example, angiogenesis was observed in stroke patients, which contained extracellular matrix degradation, endothelial cell proliferation, and new vessel formation [41]. Activation of VEGF signaling pathway played a vital role in angiogenesis [59], which had also been regarded as a promising strategy for IS treatment [60]. Consistent with Wan's report [22], DHI and NXT could both activate VEGF signaling pathway and promote angiogenesis.
Besides these common pathways, DHI was also involved in 5 IS-related pathways, including cytokine-cytokine receptor interaction, calcium signaling pathway, complement and coagulation cascades, neuroactive ligand-receptor interaction, and MAPK signaling pathway. Among them, the calcium signaling pathway participated in the regulation of agonist-stimulated contraction and myogenic tone in vascular smooth muscle cells [61]. Previous study suggested that modification of calcium activity may be a novel therapeutic strategy for cardiocerebrovascular diseases, such as IS [62]. While for NXT, the 120 putative targets were also enriched in gap junction, endometrial cancer, glioma, progesterone-mediated oocyte maturation, acute myeloid leukemia, chronic myeloid leukemia, melanoma, Fc epsilon RI signaling pathway, and non-small-cell lung cancer. Besides the 6 cancer-related pathways, gap junction may play a crucial role during the development of IS due to its vital function in neuroprotection [63,64].
Overall, it could be concluded that DHI and NXT exerted the therapeutic effects in the treatment of IS mainly through immune system pathways and VEGF signaling pathway. Furthermore, DHI could also focus on some signaling pathways, while NXT involved in gap junction. These results demonstrated that the anti-IS mechanisms of DHI and NXT had similarities and differences, which may provide important evidence for clinical choice of these two TCMPs.

. . Elucidation of the Synergistic Mechanism of DHI and NXT.
To reveal the synergistic mechanism of DHI and NXT, in this section, the 179 IS targets and 584 ingredient targets (including 201 DHI-ingredient targets and 560 NXTingredient targets) were entirely integrated into STRING database to acquire PPI data. The DHI-NXT-IS PPI network was shown in Figure 5. The average degree of all the nodes was 23.25 and 76 nodes had degree bigger than 46. Among these hub nodes, 24 were IS targets and 68 were ingredient targets. Besides the 16 overlapping targets, the remaining 52 targets were considered as indirect anti-IS targets of this TCMP pair.
The 123 putative targets of DHI and NXT (combining 71 direct anti-IS targets with 52 indirect anti-IS targets) were employed to elucidate the synergistic mechanism of DHI and NXT against IS using DAVID database. The information of top 20 pathways (Benjamini Hochberg corrected p value < 0.05) was summarized in Figure 6. Six pathways had high correlations with cardiocerebrovascular diseases, including four immune system pathways (chemokine signaling pathway, T cell receptor signaling pathway, toll-like receptor signaling pathway, and NOD-like receptor signaling pathway), one signaling pathways (VEGF signaling pathway), and one cell-related pathway (gap junction) [58]. Studies have reported that focal cerebral ischemia could result in the activation of inflammatory cytokines, chemokines, and  chemokine receptors in the immune system. These processes may be caused by the formation of thrombi in blood vessels, particularly inside microvessels [65,66]. Meanwhile, NODlike receptor signaling pathway was also closely associated with the putative targets of DHI and NXT. Recent studies have reported that NOD-like receptor played an important role in central nervous diseases, such as IS. The abnormal changes of NOD-like receptors can cause redox imbalance in the brain, which had a marked impact in the pathogenesis of IS [67]. Of note, DHI and NXT showed synergistic effect in activation of the VEGF signaling pathway, as proved experimentally by Wan's study [22]. Fortuitously, 11 cancerrelated pathways were also discerned in Figure 6. Danshensu from DHI was validated to exhibit antitumor activity by affecting on tumor angiogenesis and tumor cell invasion [68]. Additionally, butylidenephthalide, senkyunolide A, and (Z)-ligustilide from NXT showed selective cytotoxic and antiproliferative effects in human colon cancer cells [69].
In summary, DHI and NXT exerted the synergistic effect against IS by regulating 6 meaningful pathways, which indicated that the combination of DHI and NXT may provide better protective effect than monotherapy in immune system, angiogenesis, and neuroprotection.   were manually gathered from relevant researches. In order to explore the potential bioactive compounds from DHI and NXT, three important parameters were employed using in silico ADME models, including OB, Caco-2, and DL. Those ingredients with OB ≥30%, Caco-2 > −0.4 and DL ≥0.18 were preserved and adopted as candidate ingredients for further analysis. The details of all the ingredients and candidate ingredients were described in Tables S1 and 1, respectively.

. . Chemical Ingredients Collection and ADME
. . Target Identification. With names and/or chemical structures (denoted by the canonical SMILES format) as key words, the putative targets of the candidate ingredients were derived from STITCH, TCMSP, and Binding DB. The obtained targets were normalized to the official gene symbols using UniProt database (https://www.uniprot.org/) with the species limited to "Homo sapiens".
The targets related to IS were derived from five public databases, including TTD, DrugBank, PharmGKB, MalaCards, CTD with keywords "ischemic stroke", "cerebral infarction", "cerebral ischemia", "brain ischemia", and "cerebrovascular ischaemia". Then the targets were manually supplemented through a wide-scale text-mining method and annotated using the official gene symbols as described above.  The T-P network for elucidation of the synergistic mechanism of DHI and NXT. The lower 11 pathways were associated with cancers and the upper 9 pathways were related to other processes. compounds, targets, diseases, and pathways, several networks were constructed.
(1) TCMP-Compound-Target-Disease network. TCMPs, candidate ingredients of DHI and NXT and their putative targets, were employed to build this network. The IS targets targeted by the ingredients of DHI and NXT were labeled by yellow circle.
(2) Merged DHI-IS and NXT-IS PPI network and corresponding Target-Pathway (T-P) network. Based on the PPI data derived from STRING, with the species limited to "Homo sapiens" and a confidence score >0.7, two PPI networks named DHI-IS and NXT-IS PPI network were constructed for the ingredient targets of DHI and NXT interacting with IS targets, respectively. These two networks were merged to identify overlapping and distinct targets of DHI and NXT. Then the anti-IS targets of DHI and NXT were employed to conduct pathway enrichment analysis, respectively. The enriched pathway information of these targets was obtained from DAVID (used for pathway analysis according to Kyoto Encyclopedia of Genes and Genomes, KEGG), and then T-P networks composed of anti-IS targets of DHI or NXT and their corresponding pathways were created, respectively. These two T-P networks were merged to explore common and unique mechanisms of DHI and NXT for treating IS.
(3) DHI-NXT-IS PPI network and corresponding T-P network. A DHI-NXT-IS PPI network was established among DHI-ingredient targets, NXT-ingredient targets, and IS targets. Then the anti-IS targets were used to performed pathway enrichment analysis using DAVID database. The T-P network was generated to elucidate the synergistic mechanism of DHI and NXT in the treatment of IS.
All these networks were created by Cytoscape 3.6.0, open software for visualizing, integrating, modeling, and analyzing the biological networks.

Conclusions
To provide evidence for the better clinical use of TCMPs, comparison of the pharmacological mechanisms between TCMPs with similar therapeutic effect is necessary. DHI and NXT are representative TCMPs for IS treatment, which are also frequently used in combination. However, their independent and combined mechanisms have not been clearly clarified. In the present study, a feasible network pharmacology strategy based on chemical, pharmacokinetic, and pharmacological data was employed to elucidate the comparative and synergistic effects of DHI and NXT against IS. By systematic analysis of the constructed networks, we could learn that (1) the candidate ingredients of DHI and NXT not only could target the vasomotor systemrelated and blood-brain barrier-related proteins (PTGS2, F2 ADRB1, and ADRB2), but also could involve in the redox imbalance, neurotrophin activity, and brain inflammation (IL6, ALDH2, CCL2, STAT3, and AKT1). (2) DHI and NXT exerted the protective effects mainly through immune system and VEGF signaling pathways. But DHI and NXT also had their own unique pathways for treating IS, e.g., calcium signaling pathway for DHI and gap junction for NXT. (3) The synergistic mechanism of DHI and NXT was mainly associated with chemokine signaling pathway, T cell receptor signaling pathway, VEGF signaling pathway, gap junction, and so on. Overall, this study provided an effective approach for dissecting the comparative and synergistic effects of DHI and NXT in the treatment of IS.

Data Availability
The data of our research can be acquired from the Supplementary Materials uploaded with this article.