Abstract

Aim. Stroke is the second significant cause for death, with ischemic stroke (IS) being the main type threatening human being’s health. Acorus tatarinowii (AT) is widely used in the treatment of Alzheimer disease, epilepsy, depression, and stroke, which leads to disorders of consciousness disease. However, the systemic mechanism of AT treating IS is unexplicit. This article is supposed to explain why AT has an effect on the treatment of IS in a comprehensive and systematic way by network pharmacology. Methods and Materials. ADME (absorbed, distributed, metabolized, and excreted) is an important property for screening-related compounds in AT, which were screening out of TCMSP, TCMID, Chemistry Database, and literature from CNKI. Then, these targets related to screened compounds were predicted via Swiss Targets, when AT-related targets database was established. The gene targets related to IS were collected from DisGeNET and GeneCards. IS-AT is a common protein interactive network established by STRING Database. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment were analysed by IS-AT common target genes. Cytoscape software was used to establish a visualized network for active compounds-core targets and core target proteins-proteins interactive network. Furthermore, we drew a signal pathway picture about its effect to reveal the basic mechanism of AT against IS systematically. Results. There were 53 active compounds screened from AT, inferring the main therapeutic substances as follows: bisasaricin, 3-cyclohexene-1-methanol-α,α,4-trimethyl,acetate, cis,cis,cis-7,10,13-hexadecatrienal, hydroxyacoronene, nerolidol, galgravin, veraguensin, 2′-o-methyl isoliquiritigenin, gamma-asarone, and alpha-asarone. We obtained 398 related targets, 63 of which were the same as the IS-related genes from targets prediction. Except for GRM2, remaining 62 target genes have an interactive relation, respectively. The top 10 degree core target genes were IL6, TNF, IL1B, TLR4, NOS3, MAPK1, PTGS2, VEGFA, JUN, and MMP9. There were more than 20 terms of biological process, 7 terms of cellular components, and 14 terms of molecular function through GO enrichment analysis and 13 terms of signal pathway from KEGG enrichment analysis based on . Conclusion. AT had a therapeutic effect for ischemic via multicomponent, multitarget, and multisignal pathway, which provided a novel research aspect for AT against IS.

1. Introduction

13.7 million people suffered from stroke every year in the world, the second largest reason causing death with around 5.5 million dead [1, 2]. It not only killed and enabled the elderly people [3] but also made a negative effect for the youth that morbidity rate and incidence had enhanced among them. Especially, mortality rate among young people suffering from stroke in developing countries was in a uptendency, which highly increased disease burden [4]. Over past 30 years, the same situation happened in China [5], where ischemic stroke is the major type of stroke, accounting for 70.8% [6]. The ratio of mortality rate to disablility rate among the youth in China was 34.5% to 37.1% [7, 8]. Thrombolytic therapy is a common and important method to battle against ischemic stroke. The major mechanism of thrombolytic therapy is that profibrolysin is dissociated with plasmin, which can be used to degrade fibrin to dissolve thrombus. Unluckily, its application is limited and is only applied within 4.5–6 hours after suffering from it. Even if patients could be saved in time via thrombolytic therapy, there were only one-third to half of them accomplishing ischemic reperfusion [9], which might probably lead to transfer hemorrhagic stroke as a side effect.

Traditional Chinese medicine is always used without side effect to reduce disability rates, improve the quality of life, and decrease the risk of recrudesce for cerebral ischemic patients. The major mechanism involved anti-inflammatory reaction, antioxidant, inhibiting neuron apoptosis, promoting vascular regeneration, and maintaining the BBB permeability. AT was firstly recorded in Sheng Nong’s herbal classic thousands years ago, and its dried rhizomes attributed to Araceae Juss, with the effect of reducing phlegm by resuscitation, improving intelligence, dehumidifying, and appetizing [10], which was used to treat tinnitus and epicophosis, Alzheimer disease, epilepsy, depression, and stroke. AT is one important component of Ditan decoction, widely used in clinical stroke treatment, which could decrease the level of inflammatory factors and oxidative stress after cerebral ischemic and inhibit the formation of thrombus to promote stroke patients recovery [11]. Many studies found that AT and its main ingredients had a positive effect on reducing the cerebral ischemic volume [12, 13], which was related to antioxidants [14], and regulating autophagy level [13, 15].

It is unclear whether AT and AT-contained ingredients have other mechanisms against IS. Due to lack of systematic and comprehensive explanation about it, network pharmacology was used to elucidate the interactive relation among drug, targets, and disease. Potential target prediction helps us to understand more mechanisms for multiple pharmacological functions of AT and its role in biological network, which will be beneficial for the curative effect improvement [16, 17]. Many kinds of complex compounds are involved in Chinese medicine, so we applied the network pharmacology method to reveal the effect of TCM. Thus, the aim of this paper is to find out the mechanism of AT against IS in a comprehensive and systemic way.

2. Methods and Materials

The flowchart of this study is shown in Figure 1.

2.1. Active Compound Screening and Target Prediction

We searched active compounds AT contained via TCMSP Database [18] (version2.3, http://tcmspw.com/tcmsp.php), TCMID Database (v2.01, http://119.3.41.228:8000/tcmid/search/), Chemistry Database [19] (http://chemdb.sgst.cn/scdb/default.asp), and literature from CNKI to provide more useful and important compounds. Based on pharmacokinetics (absorbed, distributed, metabolized, and excreted), we set up a series of parameters, such as BBB absorption>−0.3 and DL ≥ 0.18 [20], to screen suitable compounds in TCMSP Database. Compounds could be obtained from TCMID Database and Chemistry Database without ADME parameters, so we need to get chemical formula of them from PubChem (https://pubchem.ncbi.nlm.nih.gov/) to finish Swiss ADME prediction [21], which was requested that BBB permeability was equal to YES and at least two terms of Drug-likeness were yes. All compounds were screened as above-needed target prediction through Swiss Target [22]. Probability is used to balance the connection between compounds and targets, which is more close to 1, and it is more connective. We screened targets by a median of probability to establish potential target database related to compounds in AT.

2.2. Obtaining Gene Targets Related to IS

We obtained IS-related gene targets by searching in Genecards (version 4.14, https://www.genecards.org/) and DisGeNET [2325] (version 7.0, https://www.disgenet.org/) with the keywords “ischemic stroke,” “cerebral ischemic,” “cerebral infarction.” These gene targets were put into IS-related targets database, respectively. Genes from Genecards in searching different keywords will be screened by its scores, which more than its median are constituted in IS-related targets database from Genecards.

2.3. AT-Compounds-Targets Network Construction

AT-compounds-targets visualized network was conducted by Cytoscape (version 3.7.2), where diversity of shape nodes was distinguished from AT, active compounds in AT, and targets, where edges represented the interactive relationship between AT and compounds and compounds and targets.

2.4. Extracting Core Target-Related AT-IS

By the aid of Microsoft Excel, we intersected potential target database related to compounds in AT- and IS-related target database to get the core target related to AT-IS, which was used to draw a Venn diagram online (http://www.bioinformatics.com.cn/).

2.5. Core Targets PPI Network Construction

Core targets related to AT-IS were imported into STRING (version 11.0, https://string-db.org/), selecting species “Homo sapiens”. Potential protein targets with a medium confidence score of 0.400 were designed as an interaction network, with hiding disconnect core targets to get PPI analysis results. All these results could be used in constructing AT-IS PPI visualized network, based on its degree value by analysing.

2.6. GO and KEGG Enrichment Analysis

GO and KEGG enrichment analysis of core targets related to AT-IS was accomplished by an online tool, Metascape [26] (https://metascape.org/), choosing species “Homo sapiens” and then clicking “custom analyse.” GO enrichment analysis is composed of biological process, molecular function, and cellular components. KEGG enrichment analysis is used to uncover the signal pathway showing potential therapeutic effect of some diseases. Before GO and KEGG enrichment for above core genes, “Min Overlap = 3,” “ value cutoff = 0.01,” and “Min enrichment = 1.5” should be set for passway and process enrichment and “Min Network Size = 3” and “Max Network Size = 500” for protein-protein interaction enrichment.

3. Results

3.1. Active Compounds in AT

There were 53 related active compounds satisfying ADME conditions, when filtering out repeat ones. MolID, molecular name, and Pubchem CID of AT active compounds are shown in Table 1. Structure, BBB permeability, and DL of these compounds are filled in Supplementary 3.

3.2. Compound Target Prediction Results and Its Network Construction

Probability scores were used to predict these compound targets, which median is 0.0604024588 (AT-14 (palmitic acid) and AT-22 (2-furfuraldehyde) were lower than this median; no target was related to AT-28 blumenol A; thus, these three compounds were out of the network). 1409 targets were predicted, 398 targets remained after deleting repeat ones, and specific information is shown in Supplementary Table S1. Compound-target network is presented in Figure 2. The network property of compounds in AT, degree, betweenness centrality, closeness centrality, and neighborhood connectivity is involved in Table 2. Nodes with highest connectivity have a high influence on the network [29]. We proposed that the higher the degree it had, the more significant the effect it exerted in therapy. The top 10 degree compounds were bisasaricin, 3-cyclohexene-1-methanol-α,α,4-trimethyl,acetate, cis,cis,cis-7,10,13-hexadecatrienal, hydroxyacoronene, nerolidol, galgravin, veraguensin, 2′-o-methyl isoliquiritigenin, gamma-asarone, and alpha-asarone.

3.3. AT against IS Core Targets Protein-Protein Interaction (PPI) Construction

We screened 2391 gene targets in Genecards Database, after filtering out based on median, when 560 gene targets were obtained from DisGeNET (Supplementary Table S2). Genecards Database was established based on gene-centric data from approximately 150 web sources, involved in genomic, transcriptomic, proteomic, genetic, clinical, and functional information. DisGeNET collected data from expert curated repositories, GWAS catalogues, animal models, and the scientific literature, which were identified well in relationship with disease. To avoid data redundancy, genes from Genecards and DisGeNET were intersected to improve targets relativity. In sum, 63 core gene targets remained after intersection with AT (Figure 3).

Target GRM2 was dissociated from the PPI network, without interactive relationship, which should be omitted from it. Thus, 62 remaining targets are presented in Table 3, which were used to construct the PPI visualized network in Figure 4. The lighter color and the bigger size in nodes represented the higher degree value, which made a great contribution to its effect against IS.

3.4. GO Enrichment Analysis

GO analysis of 62 potential core targets for AT against IS was performed by using the Metascape database to understand the relationship between functional units and their underlying significance in the biological system networks. The results were divided into three parts: biological processes (Figure 5(a)), cellular component (Figure 5(b)), and molecular function (Figure 5(c)), which were shown based on in statistics.

Firstly, the top 20 terms of biological process for enrichment analysis are shown in Figure 5(a), such as response to oxidative stress, postive regulation of cell migration, response to bacterium, positive regulation of response to external stimulus, coagulation, reactive oxygen species metabolic process, cellular response to nitrogen compound, fatty acid transport, neuroinflammatory response, icosanoid biosynthetic process, organophosphate biosynthetic process, cellular response to external stimulus, collagen metabolic process, positive regulation of cell death, regulation of DNA-binding transcription factor activity, and neuron death. Secondly, there were only 7 terms of cellular components enrichment analysis shown in Figure 5(c) containing membrane raft, extracellular matrix, transcription factor complex, receptor complex, ficolin-1-rich granule, perinuclear region of cytoplasm, and neuromuscular junction. Finally, 14 terms of molecular function are presented in Figure 5(b), which were lipid binding, steroid hormone receptor activity, phospholipase A2 activity, serine hydrolase activity, monocarboxylic acid binding, heme binding, oxidoreductase activity acting on NADPH, protein kinase binding, receptor regulator activity, lipopolysaccharide binding, steroid binding, phosphatase binding, protein kinase activity, and hsp90 protein bind.

3.5. KEGG Enrichment Analysis

There were 13 terms of signal pathway of AT against IS by KEGG enrichment analysis, which were ordered as ascending tendency according to -value() and contained fluid shear stress and atherosclerosis, HIF-1 signaling pathway, IL-17 signaling pathway, arachidonic acid metabolism, platelet activation, bladder cancer, inflammatory mediator regulation of TRP channels, transcriptional misregulation in cancer, PPAR signaling pathway, complement and coagulation cascades, GnRH signaling pathway, regulation of lipolysis in adipocytes, and serotonergic synapse (Figure 6). According to KEGG enrichment analysis results, we made a systematic summary for these core targets in the signal pathway in Figure 7. And we were informed for the specific role of these targets in a certain signal pathway (Figure 8).

4. Discussion

Absorbed, distributed, metabolized, and excreted are significant parameters to determine the transformation and delivery of one drug in the body. A CNS-active drug having a high BBB permeability is crucial for its antistoke effect for entering into the brain. It was evidenced than if BBB permeability was lower than −0.30, and it had no effect on entering the brain through BBB, so we set it for more than that [20]. Only 6 compounds (bisasaricin, veraguensin, lupeol, cycloartenol, asatone, and eudesmin) fitted this ADME screening in TCMSP. So it was hard to point out AT in treating IS systemically and comprehensively, and we searched on TCMID Database, Chemistry Database, and literature to supply its composition to find more potential targets. In total, we found 50 compounds in AT with effect of treating IS probably, which mainly composed of volatile oil. α-Asarone and β-asarone accounted for 95% in volatile oil, both of which could promote neuron differentiation via PI3K by aid of growth factors [30]. In our study, it was confirmed that α-asarone, with a high degree value in compound-target network (degree = 59), uncovered its core role in AT.

There were 62 core targets for AT against IS, IL6, TNF, IL1B, TLR4, NOS3, MAPK1, PTGS2, VEGFA, JUN, and MMP9 involved in the top 10 degree targets. IL6, TNF, IL1B, TLR4, and NOS3 are related to the HIF-1 signal pathway, fluid shear stress and atherosclerosis, IL-17 signal pathway, and inflammatory mediator regulation of TRP channels, which are called inflammatory factors as effect for enhancing the level of inflammatory reaction in cerebral ischemic [3133]. GO enrichment results contained some biological processes related to antistroke closely, such as response to oxidative stress, positive regulation of cell migration, coagulation, reactive oxygen species metabolic process, cellular response to nitrogen compound, fatty acid transport, and neuroinflammatory response. Similarly, some of KEGG enrichment pathways were related to secondary prevention of ischemic stroke by intervening atherosclerosis, regulation of lipolysis in adipocytes, coagulation cascades, arachidonic acid (AA) metabolism, and platelet activation. Others might have a potential therapeutic effect for IS, which were PPAR signal pathway, HIF-1 signal pathway, and inflammatory mediator regulation of TRP channels.

Atherosclerosis is the main pathological mechanism leading to cerebral ischemic, half of which has a connection with it. It was commonly known that it was beneficial to decrease the incidence or palindromia of cerebral ischemic by intervening atherosclerosis as secondary prevention [34]. In fact, fluid shear stress and atherosclerosis is a signal pathway with dual regulation for the formation of atherosclerosis, containing the biological process, antiatherosclerosis [35], and proatherosclerosis [36, 37], where nuclear transcripts Nrf2 and JUN (a part of AP-1 protein dimers) [38] are responsible for the process, respectively [39]. There was a evidence that lupeol, one component of AT, decreased lipid peroxidation level in the early stage of hypercholesterolemia artery atherosclerosis [40]. It was demonstrated that lupeol could exert a protective effect against cerebral I/R by activating Nrf2 and inhibiting p38-MAPK, with decreasing proinflammatory factors TNF-α, IL-1β, and IL-6, increasing anti-inflammatory factor IL-10, and suppressing oxidative stress level [41]. JNK and p38 MAPK promotes the expression of JUN to promote the proatherosclerosis process. Phosphorylation of p38-MAPK and JNK was inhibited by lupeol, when it suppresses the activation of microglia and astrocytes induced by LPS [42]. In addition, another component, bornyl acetate, which was found to promote HUVEC cell vitality recovery, insulted by ox-LDL, suppressing monocytes adhered to HUVEC cells, and decreasing TNF-α and IL-1β proinflammatory factors expression to antiatherosclerosis [43].

Disorder of fatty metabolism not only induces atherosclerosis but also is a high risk of stroke. It is necessary to regulate adipose lipolysis to avoid metabolism syndrome occurrence as possible to decrease the incidence of stroke indirectly. However, it required to further prove the role of adipose lipolysis in ischemic stroke for fat mice. It was confirmed that stroke rats with a high level of inflammatory factor in plasma and adipose tissue induced adipose lipolytic enzymes and free fatty acids expression increased [44]. It is a pity that amounts of compounds were conformed having a positive effect on fatty metabolism not through FABP4. Eudesmin could downregulate S6K1–H2BS36p to impair lipoblast differentiation [45]. Eugenol inhibited hepatic lipid accumulation by downregulating SREBP1 gene expression via increasing CAMKK, AMPK, and acetyl-CoA carboxylase (ACC) and suppressing phosphorylation of mammalian target of rapamycin (mTOR) and p70S6K [46]. Thymol could enhance PPARγ and PPARδ expression through overactivated p38MAPK, AMPK, and PKA to promote white adipose cell browning and increased lipid degeneration [47].

PPAR signal pathway has a close relationship with fatty metabolism, which also exerts a protective effect on cerebral ischemic. PPARα and PPARγ, nuclear transcription from nuclear receptor family, and PPARα/γ agonist decreased inflammatory reaction to induce the neuroprotective effect after cerebral ischemic [48]. α-Asarone conjugate structure activated by PPARα, increasing CPT1A gene expression related to degeneration and metabolism of fatty acid [49]. Linalool was identified as a PPARα ligand directly [50], which improved postischemic neurological scores and cognitive ability by decreasing the COX-2, IL-6, Nrf expression in cortex, and hippocampus of IS rat [51].

Antiplatelet aggregative activity and anticoagulation are common and major therapy to prevent emboli and thrombus formation for ischemic stroke. Previous studies confirmed AT suppressing platelet activation and thrombus formation [52, 53]. Platelet activation and accumulation is a cascade process, and TXA2 is a positive feedback protein, platelet agonist, a second wave mediater, which is not unnecessary for its accumulation. It can trigger platelet activated by G protein-coupled receptor [54]. In our study, we found MAPK, ERK, PLA, PTGS1, and TBXAS1 core targets involved in the process of platelet activation to regulate TXA2 produce. SQ29548, TXA2R antagonist, was found to inhibit microglia/macrophages activation and enrichment to reduce injury of cerebral ischemic [55]. It has been proved based on network pharmacology that β-asarone could be against coagulation via APP, PTGS2, and TBXAS1 [56]. Eugenol inhibited platelet activation better than elemicin, in a dose-dependent manner, which also preceded the effect of ASA-COX inhibitors such as asprin [57]. Paeonol increased NO and PGI2 expression and decreased ET-1 and TXA2 expression to inhibit platelet activation and accumulation to suppress thrombus formation [58]. 3,11-Eudesmadien-2-one [59] and α-asarone [53] also had a similar effect, which is indispensable to further identify its specific mechanism.

In addition, arachidonic acid (AA) metabolism promotes platelet accumulation. COX-2(PTGS2) makes AA degeneration, whereas PLA2 promotes AA synthesis, which takes participation in TXA2 formation leading to platelet accumulation. A series of evidence has proved that a high level of COX-2 and PLA2 was filled with MCAO animal models and stroke patients [6063]. Eugenol inhibited AA metabolism via cyclooxygenase and lipoxygenase pathways in human platelets [64]. In vitro, paeonol inhibited MAPKs activation in macrophages, decreasing iNOS, COX2, and IL6 expression to inhibit inflammatory reaction [65]. α-Iso-cubebene inhibited iNOS, COX2, and MMP9 expression and the phosphorylation of JNK and p38 to exert a neuroprotective effect, which was released by microglia amyloid beta induced [66]. However, whether paeonol or α-iso-cubebene has a neuroprotective effect for IS via inhibiting AA metabolism need to be further proved.

AT antithrombus was involved in an intrinsic way of coagulation cascades. FXa activated JNK pathway through PAR-1, which triggered F2 to induce neuron death [67]. In vitro, hippocampus clips insulted by OGD strengthened the activity of thrombus, which enhanced strength of synapsis by NAMDR; oppositely, if is inhibited thrombin/PAR-1, plasticity of synapsis was recovered [68]. It was a pity that there was no study proving AT or its components against thrombus by an intrinsic way of coagulation cascades so that we inferred it would be a novel research direction for neuroprotection of AT. Because of increasing Th-17 and decreasing Treg contributing to the pathology process of cerebral ischemic in clinic [69], we inferred that another new therapeutic point of view might be related to Th17 differentiation. In vivo, thymol decreased Th1/Treg and Th17/Treg in spleen for mice immune to Ova, which avoided overactivation of Th1 and Th17 [70].

Acute cerebral ischemic patients had a high level of HIF-1α in the early stage of clinic, indicating a serious expectation for 90 days [71]. TLR4 is one target in the HIF-1 signal pathway, where neutralizing HIF-1α weakened the increase in TLR4 in BV-2 cells in hypoxia condition with downregulation of TNF-α expression [72]. In fact, paeonol inhibited TLR4 expression, without influence of TNF-α [73]. β-Asarone had the ability of antioxidation, which decreased HIF-1α in cortex [74].

Both HRPA1 and TPR1 are members of TRP superfamily, which are structurally dependent nonselective cation channels and mediators of several signaling pathways. Methyl eugenol was proved to be hTPR1 agonist, selectively activating hTPRA1 [75]. Microglia was filled with TPRV1, which was activated to strengthen transmission of glutamate [76]. Previous study was identified that thymol improved spontaneous excitability transmission of spinal substantia gelatinosa neurons [77]. However, it is unknown whether methyl eugenol or thymol has a protective effect for IS via inflammatory mediator regulation of TRP channels, which needs more experiments to be proved. TPRV4 activated prompted expression of VEGFA and eNOS in cerebral ischemic mice, which was beneficial for proliferation and migration of neural stem cells (NPCs) and angiogenesis [78]. TPRV4 agonist impaired cerebral in IS rats by upregulating PI3K/AKT and downregulating p38MAPK (MAPK14) [79].

Taken together, these signal pathways, related to AT preventing and treating IS, could be approximately classified into three categories as follows: atherosclerosis, regulation of lipolysis in adipocytes, and PPAR signal pathway were associated with lipid metabolism; coagulation cascades, arachidonic acid (AA) metabolism, and platelet activation had a connection with the therapy of antiplatelet aggregative activity and anticoagulation; and HIF-1 signal pathway, inflammatory mediator regulation of TRP channels, and Th17 differentiation might be the potential therapeutic signal pathway. In addition, we inferred that an intrinsic way of coagulation cascades and Th17 differentiation was probably new therapeutic or preventive direction for the role of AT against IS.

5. Conclusion

In summary, we explored and explained multiple compounds, multiple pathways, and multiple targets of AT-regulated ischemic stroke treatment based on a network pharmacology method. Our data indicated that amounts of compounds contained in AT might prevent and treat ischemic stroke through some signal pathways, most of which need an elaborated detailed mechanism in the related signal pathway. Two novel research directions (IL-17 signaling pathway and complement and coagulation cascades) of AT therapeutic for preventing IS were proposed. In the future, more studies should focus on providing experimental evidence and enhancing the effect of AT in the cerebral ischemic on a comprehensive level and improving ability of AT targeting to the brain and crossing the BBB.

Data Availability

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

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

FZ L conceived and designed the study. FZ L and Q Z wrote the paper. SX L accomplished statistical data. YZ X, DR Z, and S W gathered active components and target information. FZ L, LQ Z, and Y G reviewed and edited the manuscript. All authors read and approved the manuscript.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (Grant no. 81573926) and National Key R&D Program of China (Grant no. 2018YFC1705000).

Supplementary Materials

Table S1: compound targets for each component in AT by prediction. Table S2: IS-related targets in GeneCards and DisGeNET. Table S3: specific property of compounds contained in AT. (Supplementary Materials)