Peng Zeng1†
Yao Yi1†
Hong-Fei Su1
Chao-Yuan Ye1
Yi-Wen Sun1
Xin-Wen Zhou1
Youming Lu1
Anbing Shi2*
Qing Tian1*- 1Key Laboratory of Neurological Disease of National Education Ministry, Department of Pathology and Pathophysiology, School of Basic Medicine, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
- 2Department of Biochemistry and Molecular Biology, School of Basic Medicine, Tongji Medical College, Cell Architecture Research Institute, Huazhong University of Science and Technology, Wuhan, China
Presently, the treatment options for ischemic stroke (IS) are limited due to the complicated pathological process of the disease. Chuanxiong Rhizome (CR), also known as Conioselinum anthriscoides “Chuanxiong” (rhizome), is the most widely used traditional Chinese medicine for treating stroke. This study aimed to uncover the key phytochemicals and biological functions of CR against IS through a network pharmacology approach combining with IS pathophysiology analysis. We employed permanent unilateral common carotid artery ligation to construct a mouse model of global cerebral ischemia and found that cerebral ischemia injuries were improved after 7 days of gavage treatment of CR (1,300 mg/kg/day). CR exerts protective effects on neurons mainly by acting on targets related to synaptic structure, synaptic function, neuronal survival and neuronal growth. A total of 18 phytochemicals from CR based on UHPLC-MS/MS that corresponded to 85 anti-IS targets. Coniferyl ferulate, neocnidilide and ferulic acid were identified as the key phytochemicals of CR against IS. Its brain protective effects involve anti-inflammatory, anti-oxidative stress, and anti-cell death activities and improves blood circulation. Additionally, the two most important synergistic effects of CR phytochemicals in treating IS are prevention of infection and regulation of blood pressure. In brain samples of Sham mice, L-tryptophan and vanillin were detected, while L-tryptophan, gallic acid, vanillin and cryptochlorogenic acid were detected in IS mice by UHPLC–MS/MS. Our findings provide a pathophysiology relevant pharmacological basis for further researches on IS therapeutic drugs.
Introduction
Ischemic stroke (IS), which accounts for 87% of all cerebral strokes, is a leading cause of neurological morbidity and mortality worldwide (Benjamin et al., 2018). It can be induced by a variety of events, such as cardiac ischemia, occlusion of cerebral small blood vessels and carotid atherosclerosis. With the accelerated growth of the aging population (65 years and older), the incidence of IS is increasing yearly, and the overall burden is shifting to a younger population, especially in third world countries (Wang W. et al., 2017). The course of IS is broadly divided into the acute phase (from minutes to 6–12 h), sub-acute phase (from 6–12 h to days) and recovery phase (after more than 3 weeks). A variety of complex molecular and cellular interactions drive the pathophysiological effects of IS. Reduced blood flow to the brain causes stress, cell death (necrosis and apoptosis) and loss of neuronal function. Inflammation, oxidative stress, acidosis, overload of intracellular calcium, excitotoxicity, free radical injury, cytokine injury, complement activation, blood-brain barrier (BBB) impairment and activation of astrocytes and microglia are also key events contributing to IS pathology.
Treatment options for IS are currently limited. The recanalization of blood flow in acute IS relies on the administration of tissue plasminogen activator (t-PA) and mechanical thrombolytic technology in time, which have shown some benefits in helping patients recover from IS. However, the shortcoming of these treatments is that reperfusion can lead to the production of highly harmful reactive oxygen species (ROS) and oxidative stress. This is the major cause of reperfusion damage after ischemic injury. Oxidative stress causes apoptosis, autophagy, and necrosis in the brain. The ischemia time is a good biomarker for brain tissue viability. Clinically, the therapeutic window for tPA treatment is very limited, as tPA must be administered within 6 h of the onset of symptoms (van der Worp and van Gijn, 2007). Additionally, tPA may promote cerebral hemorrhage and anaphylaxis (Donnan et al., 2008). Therefore, there is still a great need to find therapeutic agents for IS during ischemia and recanalization, extend the therapeutic window and further improve clinical outcomes.
In clinical practice, Chinese herbal medicine has been shown to be effective in treating IS (Yuan et al., 2008). From nearly 1,000 anti-stroke Chinese medicine prescriptions, 192 anti-stroke herbs were identified (Liu et al., 2017). Chuanxiong Rhizome (CR), also known as Conioselinum anthriscoides “Chuanxiong” (rhizome), an herb first recorded in the Shennong’s Classic of Materia (simplified Chinese: 神农本草经), is the most frequently used herb for stroke (509 times) (Liu et al., 2017). CR is also used in treating cardiovascular diseases, respiratory diseases, pain, and trauma (Chen et al., 2018). Tetramethylpyrazine, a natural alkaloid extracted from CR, has been used extensively for the treatment of IS (Lin et al., 2021; Zhu et al., 2021). Z-ligustilide, a major component of CR, can significantly decrease the infarct volume, and mitigate neurological dysfunction in rats with middle cerebral artery occlusion (MCAO) (Li et al., 2017). CR contains numerous phytochemicals, such as a variety of aromatic acids, phthalides and alkaloids (Zhang et al., 2019; Wang et al., 2020). Analyzing and summarizing the effective phytochemicals of CR and their mechanism of action against stroke is of great significance for understanding the key intervention links in the treatment of stroke and the development of related drugs.
As a systematic research method (Hopkins, 2008; Zeng et al., 2021a; 2021b), network pharmacology can be used to study the effects of drugs in a systematic and comprehensive way. The aim of this research was to investigate the anti-IS effects of key active phytochemicals of CR and determine the anti-IS mechanism of CR based on network pharmacology research and pathophysiological changes in IS. To improve the quality of collected data, the bioactive phytochemicals of CR were obtained from a recent study based on ultra-HPLC with triple quadrupole MS (UHPLC-MS/MS) (Wang et al., 2020). Based on the pathophysiology of IS and 18 bioactive phytochemicals of CR corresponding to 85 anti-IS targets, CR was found exert neuroprotective, anti-inflammatory, anti-oxidative stress, anti-cell death and improved blood circulation effects to protect the brain. Coniferyl ferulate, neocnidilide and ferulic acid were suggested as the key phytochemicals of CR against IS. Moreover, CR phytochemicals synergistically prevent infection and regulate blood pressure in treating IS. In brain of CR treated IS mice, four components of CR, e.g., L-tryptophan, gallic acid (GA), vanillin and cryptochlorogenic acid (CA), were detected. The mechanisms by which these phytochemicals antagonize IS induced brain damage through intracerebral and extracerebral pathways require more in-depth and systematic studies.
Materials and Methods
Determination of the Main Phytochemicals of CR and ADME Evaluation
A total of 20 bioactive phytochemicals of CR were obtained from a recent study based on UHPLC-MS/MS (Wang et al., 2020). In addition, tetramethylpyrazine, a main bioactive phytochemical of CR, was also included (Lin et al., 2021; Zhu et al., 2021). Canonical SMILES of the 21 bioactive phytochemicals in CR were extracted from PubChem database (https://pubchem.ncbi.nlm.nih.gov/) (Kim et al., 2016). The 2-dimensional chemical structures were generated by ChemDraw Ultra 8.0 software. Lipinski’s rule of five (RO5), i.e., molecule weight (MW) < 500, number of hydrogen bond donors (Hdon) ≤ 5, number of hydrogen bond acceptors (Hacc) ≤ 10, lipid-water partition coefficient (LogP) ≤ 5 and number of rotatable bonds (Rbon) ≤ 10, has been extensively used to evaluate bioavailability based on the structures of compounds (Lipinski et al., 2001). Here, we employed the SwissADME web tool (www.swissadme.ch) (Daina et al., 2017; Zeng et al., 2021c) to evaluate the compounds according to RO5. The main phytochemicals of CR complied with RO5 were filtered out for follow-up studies. The contents of phytochemicals (mg/g) of CR were obtained from 36 CR dried rhizome samples from six different production origins. GraphPad Prism software (version 8.0, San Diego, CA, United States) was used for graphical visualization.
Target Fishing
Collection of CR-Related Targets
The targets of the bioactive phytochemicals of CR were obtained using SwissTargetPrediction (http://www.swisstargetprediction.ch/) (Daina et al., 2019) based on their structures. Specifically, canonical SMILES were input into the SwissTargetPrediction, and the target species was set as Homo sapiens. Subsequently, target information was collected and organized using Microsoft Excel software (version 2019, Redmond, WA, United States). All software used in this study was performed on Windows (version 10).
Screening of Targets of CR Against IS
The disease associated with targets were retrieved from the Pathway Assembly from Literature Mining-an Information Search Tool (PALM-IST) database (http://www.hpppi.iicb.res.in/ctm/) (Mandloi and Chakrabarti, 2015) and GeneCards database (https://www.genecards.org/) (Stelzer et al., 2016) using “ischemic stroke” as a key phrase. Specifically, top 200 relevance score target genes were selected for further study in the GeneCards database. After removing duplicate targets, Venn diagrams of overlapping CR targets and disease-related targets were generated using Venny 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/index.html). The common targets represented the targets of CR against IS. We further investigated the functional classification of the CR targets associated with IS using the Panther classification system (http://pantherdb.org/) (Mi et al., 2007). Sankey diagrams were plotted using OriginPro 2021 software (OriginLab Corporation, Northampton, MA, United States).
Gene Ontology (GO) and the KEGG Pathway Enrichment Analysis
KEGG pathway database is a collection of manually drawn pathway maps of molecular interactions. GO and KEGG pathway enrichment analyses were performed using the ClusterProfiler R package (version 3.12.0 for Windows) (Yu et al., 2012). GO terms were divided into three categories: biological process (BP), cellular components (CC), and molecular function (MF). The p values were adjusted using a BH approach. Statistical significance was denoted if adjust p value <0.05. The top 10 GO terms and top 20 KEGG pathways sorted by the p value were visualized using an online tool (http://www.bioinformatics.com.cn/).
Protein-Protein Interaction (PPI) Network Construction and Molecular Complex Detection (MCODE) Clustering Analysis
The PPI network of the targets was constructed using the latest version of STRING database (version 11.5, https://string-db.org/) (Szklarczyk et al., 2019) and visualized using Cytoscape software (version 3.7.1) (Shannon et al., 2003). By default, only PPIs with an interaction score exceeding the threshold of 0.4 were included. In the PPI network, degree refers to the number of other nodes directly connected to a node. The higher the degree is, the more important the node is in the PPI network. The degree values were calculated using Network Analysis (a Cytoscape plugin), and the top 10 targets ranked by degree were selected and identified as core targets.
The Cytoscape plugin MCODE (http://apps.cytoscape.org/apps/mcode) (Bader and Hogue, 2003) was employed to identify highly interconnected clusters in the PPI network. The following MCODE criteria were used for selection: degree cutoff = 2, node score cutoff = 0.2, k-core = 2, and max depth = 100. The node with the highest weighted vertex was defined as seed node (key target in the cluster) by MCODE. Moreover, a CR main component target-IS target network was constructed using Cytoscape software (version 3.7.1).
Molecular Docking Simulations
The 3D molecular structure of CR key phytochemicals was retrieved from the PubChem database and the structure files of target proteins were acquired from the RCSB Protein Data Bank (PDB database, http://www.rcsb.org/) (Westbrook et al., 2002). Crystal structures of PTGS2 (PDB ID: 5F19), IL1B (PDB ID: 5R85), CXCL8 (PDB ID: 1ILQ) and MMP9 (PDB ID: 2OW0) were employed for molecular docking using the SwissDock web service (http://www.swissdock.ch/docking) (Grosdidier et al., 2011).
Animals and Treatments
Male C57BL/6J mice (3-month old, n = 20) were from Hua Fukang Experimental Animal Center (Certificate number: SCXK (Jing) 2019-0022, Beijing, China). Ethics approval was received from the Animal Care and Use Committee of Huazhong University of Science and Technology. The animal house was kept at a constant relative humidity (55 ± 15%) and room temperature (22 ± 2°C). Mice were maintained in 12-hour-light/-dark cycle. Blood supply of brain comes from four arteries, namely two internal carotid arteries (ICAs) and two thinner vertebral arteries. Two ICAs bring about two-thirds of the blood supply into the brain. We ligated the right common carotid artery of mouse by the 5/0 surgical wire to establish a cerebral ischemia model (IS mice). Briefly, the animal was intraperitoneally anesthetized by a mixture of ketamine hydrochloride and xylazine (150 mg/ kg ketamine and 25 mg/ kg xylazine body weight). Carefully expose the right common carotid artery and permanently ligate with surgical wire. Avoid damaging the vagus nerve during the operation. Sham mouse received the same surgical exposure of right common carotid artery without ligating.
CR formula granules (ZGB 2021-111, 2 g/package) from Jiangyin Tian Jiang Pharmaceutical Co., Ltd (Jiangyin City, China) were dissolved in boiling water at a concentration of 65 mg/ ml. To evaluate the phytochemicals of CR in brain, IS mice and Sham mice were given 7-days’ CR (1,300 mg/ kg/ day) orally after surgery. The dose of CR in mice is calculated based on the clinical dose (Zhang et al., 2017). Brain samples were collected 1 h after the last gavage and stored at −80°C.
Rotarod Test and Statistical Analysis
The rotarod test was conducted to evaluate rodent motor coordination. The mice were placed on the rotating rod for adaptive training before surgery, and the training lasted for 3 days, four rounds per day, with the rotating rod speed of 5, 10, 20, and 30 revolutions per minute (rpm) for each round of training, and the mice were trained for 180 s at each speed. The mice were tested at 10, 20, and 30 rpm in the rotating rod experiment 3 days after surgery, and the time of the mice staying on the rotating rod was recorded. In addition, the changes in body weight were recorded. Statistical analysis was performed using SPSS 12.0 (SPSS Inc., Chicago, IL). The 1-way ANOVA procedure followed by Tukey’s multiple comparisons test was used to determine the differences between groups. The level of significance was defined as p < 0.05.
Triphenyltetrazolium Chloride Staining and HE Staining
TTC was obtained from Solarbio (T8170, Solarbio, Beijing, China). The mice were anesthetized and the integral brains were quickly obtained and cut into 2 mm tissue slices, stained with 2% TTC at 37°C for 20 min and soaked in 4% formaldehyde for fixation 6 h. The brain slices were arranged in order and photographed. HE staining was conducted according to routine protocols (Zeng et al., 2019). HE Staining Kit was obtained from Solarbio (G1120, Solarbio, Beijing, China). Images were acquired under a microscope (Olympus SV120, Tokyo, Japan). Cell nuclei were stained blue, while the cytoplasm and extracellular matrix were stained red in HE staining images.
Quantitative Analysis of Key Phytochemicals of CR in Brain
The contents of 21 main bioactive phytochemicals of CR in the brain samples was analyzed by Beijing Qingxi Technology Research Institute (Beijing, China). Place 100 mg brain sample in a 15 ml centrifuge tube containing 2 ml of acetonitrile, and take 1 ml of supernatant after ultrasonication for 5 min to centrifuge at 14,000 rpm for 5 min. After the materials in the supernatant were blown dry with nitrogen, they were dissolved in 500 μL of 50% methanol aqueous solution and filtered with a 0.22 μm microporous membrane. Separation was performed on the ACQUITY UPLC HSS T3 column (2.1 × 100 mm, 1.8 µm). The flow rate was 0.3 ml/ min, the injection volume was 10 μL, and the column temperature was 35°C. Deionized water with 0.1% formic acid (phase A) and acetonitrile with 0.1% formic acid (phase B) were used as mobile phases. Data were acquired in full MS-ddMS mode using a Q Exactive Orbitrap mass spectrometer (Thermo Fisher Scientific, Bremen, Germany). The full scan (MS1) range was 100–1,200 m/ z with a resolution of 70,000, and MS2 has a resolution of 17,500. The ion source electrospray voltage was 3.2 kV, the capillary temperature was 320°C, and aux gas heater was 350°C. The sheath gas flow rate was 40 L/ min, and auxiliary gas flow rate was 15 L/ min. AGC target was set at 1e6, TopN was set at 5, and stepped normalized collision energy (NCE) of 30, 40, and 50 was used for fragmentation. By using Compound discover 3.2 software, the characteristic peaks of Raw mass spectra were extracted, and the mass deviations of the characteristic peaks for element matching, molecular formula prediction and isotope distribution matching were all set within 5 ppm. The mzcloud online database (https://www.mzcloud.org/) and the local mzVault natural products database were used to identify the characteristic peaks. The screening criteria of the positive results were mass deviation <5 ppm, coincidence isotope distribution and mzVault best match score >70.
Results
Main Phytochemicals of CR and Their Potential Biological Functions
This study included 21 main bioactive phytochemicals of CR, namely 12 aromatic acids, eight phthalides, and tetramethylpyrazine (Wang et al., 2020; Yin et al., 2019). According to the UHPLC-MS/MS data (Wang et al., 2020), the component with the highest content is z-ligustilide (7.16 ± 0.36 mg/g CR), followed by senkyunolide A (5.2 ± 0.36 mg/g CR) and coniferyl ferulate (1.98 ± 0.08 mg/g CR) (Figure 1). The content of tetramethylpyrazine in CR is very low (0.60 to 11.75 μg/ g CR) (Yin et al., 2018). The chemical structures of these 21 phytochemicals of CR were downloaded from the PubChem database (Kim et al., 2016). The SwissADME database (Daina et al., 2017; Zeng et al., 2021c) was employed to evaluate the phytochemicals according to RO5, a classical set of criteria used to evaluate the drug-likeness of a compound (Lipinski, 2004). Since three out of 21 phytochemicals, namely chlorogenic acid, CA and 3,5-O-Dicaffeoylquinic acid, violated RO5 (Table 1), we focused on the remaining 18 CR phytochemicals, which were consistent with RO5, in the follow-up studies.
FIGURE 1. The contents of bioactive phytochemicals (mg/g) of Chuanxiong Rhizoma (CR). The data were retrieved from the reference (Wang et al., 2020), which used 36 CR dried rhizome samples from six different production origins in Sichuan province. The data are presented as means ± S.E.M.
TABLE 1. Pharmacological and molecular properties of the main phytochemicals in CR.
The SwissTargetPrediction database (Daina et al., 2019) was used to study the potential targets of the 18 CR phytochemicals based on their structures, and a total of 376 potential targets were obtained. To provide insight into the biological functions of these 376 targets, we performed GO and KEGG pathway enrichment analyses using the ClusterProfiler R package with the adjusted value of p < 0.05 (Yu et al., 2012). The GO terms included BP, CC and MF. Analysis of BP terms showed that most of the targets were related to peptidyl-tyrosine modification (GO:0018212), peptidyl-tyrosine phosphorylation (GO:0018108), response to lipopolysaccharide (LPS) (GO:0032496), G protein-coupled receptor signaling pathway (GO:0007187), cellular calcium ion homeostasis (GO:0006874) and so on. It is noteworthy that most of the CR targets were mainly localized in synaptic membranes, such as integral component of synaptic membrane (GO:0099699), intrinsic component of synaptic membrane (GO:0099240), integral component of presynaptic membrane (GO:0099056) and membrane region (GO:0098589). The top 10 enriched cell components were associated with 70 CR targets, which formed a PPI network with 69 nodes and 397 edges (Supplementary Figure S1A). Among these targets, amyloid precursor protein (APP), GRM5, SRC, GABBR1, GRM1, CNR1, HTR2A, CTNNB1, epidermal growth factor receptor (EGFR), and FYN have important roles in the PPI network (Supplementary Figure S1B). KEGG pathway analysis revealed that these targets were enriched in multiple synapse related KEGG pathways (Supplementary Figure S1C). According to analysis of MF terms, CR related targets are primarily involved in neurotransmitter receptor activity (GO:0030594), protein tyrosine kinase activity (GO:0004713), transmembrane receptor protein tyrosine kinase activity (GO:0004714), G protein-coupled amine receptor activity (GO:0008227), and transmembrane receptor protein kinase activity (GO:0019199) (Figure 2A). These data indicate that CR exerts protective effects on synaptic structure and synaptic function by targeting multiple synapse related targets.
FIGURE 2. Enrichment analysis of the CR potential targets. GO functional (A) and KEGG pathway enrichment analysis (B) of 376 targets of CR main phytochemicals were identified using the ClusterProfiler package, in which the significance threshold is set as p value < 0.05. For the bubble plot: X-axis, rich factor (the ratio of genes in the background pathways); bubble size, the number of genes enriched; bubble color, p value. BP: biological process; CC: cellular component; MF: molecular function.
The top 20 most significantly enriched KEGG pathways for the CR targets according to KEGG pathway enrichment analysis are shown in Figure 2B. The main KEGG pathways were the neuroactive ligand-receptor interaction (hsa04080), cAMP signaling pathway (hsa04024), AGE-RAGE signaling pathway in diabetic complications (hsa04933), EGFR tyrosine kinase inhibitor resistance (hsa01521), serotonergic synapse (hsa04726), nitrogen metabolism (hsa00910), endocrine resistance (hsa01522), calcium signaling pathway (hsa04020) and HIF-1 signaling pathway (hsa04066). These results revealed that CR targets are involved in neuronal survival and neuron growth related signals.
Anti-IS Targets of CR and Their Functional Classifications
Through retrieving the PALM-IST database (Mandloi and Chakrabarti, 2015) and GeneCards database (Stelzer et al., 2016), altogether, 683 IS-related targets were obtained. To further analyze the underlying mechanism of CR against IS, 85 targets (Table 2) were obtained based on the intersection of CR targets and IS-related targets (Figure 3A). Based on their cellular functions, these 85 targets involved in the anti-IS effect of CR were categorized into seven different classes (Figure 3B). Metabolite interconversion enzyme (PC00262, 27.8%), protein modifying enzyme (PC00260, 22.2%) and transmembrane signal receptor (PC00197, 19.4%) ranked among the top three in terms of the all classes. 20 targets involved in metabolite interconversion enzyme, and these targets formed a complex PPI network (20 nodes and 54 edges) (Figure 3C). Arachidonate 15-lipoxygenase (ALOX15), arachidonate 5-lipoxygenase (ALOX5), CYP2C19, CYP2C9, CYP3A4, IDO1, MPO, PTGS1 and PTGS2 are oxygenases. Among the 16 targets identified as protein modifying enzymes, ACE, ENPEP, MMP1, MMP2, MMP3, MMP8 and MMP9 are metalloproteases, DPP4, ELANE, F2 and F7 are serine proteases (Figure 3D). The above data suggested that CR can stabilize intracellular proteins and the extracellular matrix by regulating various enzymes.
TABLE 2. The target information of CR against IS.
FIGURE 3. Protein functional classification of targets of CR against IS. (A) Venn diagram showing overlapping targets between the CR and IS targets. (B) Protein functional classification of 85 common targets. Numbers above bars represent the percentage of the protein in the given functional class. (C) PPI network of 20 targets involved in metabolite interconversion enzyme (PC00262). Higher the degree of a node, the darker the color of the node. The edges thickness is proportional to the combined score. (D) Sankey diagram of 16 targets involved in protein modifying enzyme (PC00260). The bar chart below the sankey diagram represents the number of targets involved in each classification.
STRING database (version 11.5) based PPI analysis (Szklarczyk et al., 2019) was performed to investigate the relationship between the 85 anti-IS targets. The PPI network contained a total of 85 nodes and 920 edges, and the average node degree was 21.6 (Figure 4A). When the targets were sorted by degree, ALB, VEGFA, PTGS2, MAPK1, IL1B, CXCL8, SRC, MMP9, EGFR and MMP2 were identified as the core targets (Figure 4B). Furthermore, through MCODE analysis (Bader and Hogue, 2003), four functional clusters that may be the most relevant to the anti-IS mechanism of CR were identified from the 85 anti-IS targets (Figures 4C–F). Cluster one contained 28 nodes and 328 edges with a score of 24.3 (Figures 4C,G). The seed node of this cluster was EGFR (Figure 4C), which plays a major role in cell growth, proliferation, survival and differentiation. Cluster one involved targets that are mainly involved in the inflammatory response (GO:0006954), and MMP9, PTGS2, CXCL8, TLR4, SERPINE1, IL1B may be potential therapeutic targets in the effect of CR on the inflammatory response. CXCL8, IL1B, MMP9 and PTGS2 are also core targets, which indicates the importance of the inflammatory response in the therapeutic effect of CR on IS. Lowest energy docking model based molecular docking analysis was applied to validate the binding of the core targets (CXCL8, IL1B, MMP9, PTGS2) and CR phytochemicals. The results showed that all the core targets were well combined with the corresponding CR phytochemicals (Table 3). Among these targets, senkyunolide I showed the highest binding energy with PTGS2, with score values of −5.53 kcal/ mol; rulic acid showed the highest binding energy with PTGS2 and MMP9, with score values of −5.11 and −7.04 kcal/ mol, respectively; coniferyl ferulate showed the highest binding energy with MMP9, with score values of −7.91 kcal/ mol; nocnidilide and lvistilide A showed the highest binding energy with IL1B and CXCL8, respectively.
FIGURE 4. MCODE network clustering analysis of common targets between the CR and IS. (A) The PPI network was constructed for the 85 potential targets of CR against IS. The darker the node, the greater the degree. The edges thickness is proportional to the combined score. (B) The top 10 core targets were extracted from (A). (C–F) The four tightly connected network clusters were identified by MCODE algorithm. The seed node of each cluster is indicated by red font. (G) Comparison of MCODE scores for different clusters.
TABLE 3. Molecular docking of the core targets with its corresponding CR phytochemicals.
Cluster two contained eight nodes and 12 edges with a score of 3.4 (Figures 4D,G) and was primarily involved in blood circulation (GO:0008015). The seed node of this cluster was APP, which encodes a cell surface receptor and transmembrane precursor protein that is cleaved by secretases to form various peptides (Figure 4D). Cluster 3 contained 6 nodes and 8 edges with a score of 3.2 (Figures 4E,G). The seed node of this cluster was CREB binding protein (CREBBP) (Figure 4E). CREBBP (Gouveia et al., 2017), NFE2L2 (Kidani et al., 2020) and PARP1 (Liu et al., 2021) are known for their roles in the pathogenesis of stroke and post-stroke neurovascular remodeling and functional recovery. Cluster 4 contained 8 nodes and 11 edges with a score of 3.1 (Figures 4F,G) and ALOX5AP, ALOX15, CYP2C9, CYP3A4, TBXAS1, HTR1A, TTR and BCHE are oxidative stress related targets. ALOX15 is the seed node of this cluster (Figure 4F). Of note, ALOX15 is up-regulated following IS and contributes to neuronal cell death and hemorrhagic transformation (Gaberel et al., 2019). These data suggest that CR has anti-inflammation, anti-oxidative stress, anti-cell death activities and improves blood circulation to treat IS.
Key Phytochemicals of CR in Treating IS
To identify the key phytochemicals of CR in treating IS, a CR main component target-IS target network was constructed using Cytoscape software (version 3.7.1) (Shannon et al., 2003). The interaction network comprised 18 active phytochemicals of CR and 85 corresponding targets (shown as gene symbols) (Figure 5A). The average number of targets per CR component was 12.1, and the mean of degree of phytochemicals per target was 2.6. These results fully reflect the multicomponent and multitarget characteristics of CR in the treatment of IS. The results showed that coniferyl ferulate, neocnidilide and ferulic acid have the most targets (degree = 26), followed by caffeic acid (degree = 23), levistilide A (degree = 19), senkyunolide A (degree = 15) and z-ligustilide (degree = 14) (Figure 5A). The chemical structures of CR phytochemicals with top 10 anti-IS targets are presented in Supplementary Figure S2. Coniferyl ferulate acts on four core targets, including EGFR, MMP2, MMP9 and SRC. Neocnidilide acts on two core targets (IL1B and PTGS2) and ferulic acid acts on four core targets (EGFR, MMP2, MMP9 and PTGS2). Interestingly, coniferyl ferulate, neocnidilide and ferulic acid have many common targets (Figure 5B). Furthermore, they all act on the inflammatory response (GO:0006954) and oxidative stress (GO:0006979). Coniferyl ferulate and ferulic acid both regulate cytokine production (Figure 5C). Thus, coniferyl ferulate, neocnidilide and ferulic acid are the key phytochemicals of CR in treating IS.
FIGURE 5. Key phytochemicals of CR in treating IS. (A) CR main phytochemicals target-IS target network. Red nodes represent the main phytochemicals in CR, and red numbers below these nodes represent the number of CR targets against IS. Green circular nodes represent 85 common targets between the CR and IS. (B) Venn plot showing the number of overlapping potential targets among the three key ingredients. (C) Enriched GO biologic process for top three key phytochemicals -related targets. Numbers above bars are the number of targets enriched in terms.
Potential Synergistic Effects of CR Phytochemicals in Treating IS
As shown in Figure 5, the phytochemicals of CR have the same targets the different targets, suggesting that CR may act on different pathological processes of IS. Therefore, we used the clusterProfiler R package to determine the potential synergistic effects of CR phytochemicals in treating IS according to the 85 IS associated targets of CR by analyzing their enrichment for GO BP terms and KEGG pathways.
The 85 shared targets of CR and IS are significantly involved in response to LPS (GO:0032496), regulation of blood pressure (GO:0008217), vascular process in circulatory system (GO:0003018), extracellular structure organization (GO:0043062), regulation of inflammatory response (GO:0050727) and so on (Figure 6A). The endotoxin, LPS, triggers inflammatory cascade, complement activation, prooxidative stress, and tissue destruction. 22 IS associated CR targets are involved in the response to LPS (GO:0032496) and form a PPI network that includes 21 nodes and 137 edges (Figure 6B). Notably, ICAM1, VCAM1, MAPK1, CXCL8 and PTGS2 play a key role in this PPI network (Figure 6B). Infection is the leading cause of death in patients following stroke. Acute IS is followed by profound immunoreactions, including an inflammatory response and subsequent stroke-induced immunodepression syndrome (SIDS). SIDS further promotes infection complications, such as the most common pneumonia, and brain damage. The regulation of inflammatory response (GO:0050727), neuroinflammation response (GO:0150076), reactive oxygen species metabolic process (GO:0072593) and leukocyte migration (GO:0050900) are all important for the treatment of SIDS.
FIGURE 6. GO biological process (BP) enrichment analysis of targets of CR against IS. (A) Bubble plot of representative GO enrichment terms of common targets between the CR and IS. Only functional BP terms passing adjust p value threshold <0.05 were considered. X-axis, rich factor; bubble size, the number of genes enriched; bubble color, p value. (B, C) PPI network of targets involved in response to lipopolysaccharide (GO:0032496) and regulation of blood pressure (GO:0008217). Higher the degree of a node, the darker the color of the node. The edges thickness is proportional to the combined score.
An abrupt increase in blood pressure is the most common clinical symptom of acute IS. As blood pressure changes during acute IS vary with the time of onset and between patients, the optimal blood pressure management after acute IS remains a challenge. As shown (Figure 6A), regulation of blood pressure (GO:0008217) was the second most enriched term, with the PPI network including 18 nodes and 46 edges (Figure 6C). Additionally, the vascular process in circulatory system (GO:0003018), regulation of blood vessel diameter (GO:0097746) and regulation of blood vessel size (GO:0050880) (Figure 6A) all contribute to blood pressure management. According to these results, infection prevention and blood pressure regulation are the two important synergistic effects of CR phytochemicals in treating IS.
CR Improves Cerebral Ischemic Injury and Its Phytochemicals Entering the Brain
Twenty-four hours after establishing the global cerebral ischemia model, cerebral infarction was observed in IS mice. The infarcted brain tissue was pale, and normal brain tissue was red (Figure 7A). The cerebral ischemia injury was improved after 7 days of gavage treatment of CR (1,300 mg/ kg/ day) (Figure 7B). Intragastric treatment with CR can improve weight loss in IS mice (Figure 7C). To assess cerebral ischemia-induced impairment of motor coordination, the rotarod test was performed. Compared with Sham and IS + CR mice on day 3, the IS mice showed a significant reduction in retention time (Figure 7D). Neuronal morphology of the mice cerebral cortex was observed 24 h after cerebral ischemia surgery by HE staining (Figure 7E). Cortical neurons in Sham mice after CR treatment were arranged regularly, and the structures of neurons were clearly with round, large and regular nuclei. Most of the neurons of IS mice were arranged disorderly, and the nucleus was pyknotic or severely shrunken. These results indicated that treatment with CR may improve brain injuries caused by cerebral ischemia.
FIGURE 7. CR treatment significantly ameliorates cerebral ischemic injury. Representative images of TTC-stained brain slices obtained 24 h (A), or 7 days (B) after global cerebral ischemia surgery. (C) Body weight changes were monitored over a 7-day period. (D) The time on the rod in the rotarod test on day 3 (n = 6/group). ** p < 0.01, *** p < 0.001 vs. Sham; ## p < 0.01, ### p < 0.001 IS + CR vs. IS. (E) Representative images of HE staining of the cortex after 7 days of CR treatment. Normal cortical neurons revealed normal cells with rich cytoplasm and round and slightly stained nucleus with relatively large and clear nucleolus formation. Scale bar = 50 µm. (F) The phytochemicals identified by UHPLC/MS/MS in the brains of Sham and IS mice that were orally administered CR. (G) Four components detected only in the brains of CR treated IS mice.
After 7 days of CR treatment, less cellular damage was observed. Studying components of CR in the brain may help to understand its cerebral protections. Therefore, an UHPLC/MS/MS method was employed to detect the relative contents of CR components in the brain samples from CR treated Sham and IS mice as shown in Supplementary Figures S3A,B. A total of 148 chemical constituents were screened out. In brain samples of CR treated Sham mice, L-tryptophan (0.091%) and vanillin (0.001%) were detected, suggesting they enter the brain tissue under physiological condition. In CR treated IS mice L-tryptophan (0.092%), GA (0.001%), vanillin (0.001%) and CA (0.013%) were detected (Figure 7F). In the UHPLC/MS/MS data, we also found the increased contents of hypoxanthine (33.82%) and adenine (0.221%) in brains of IS mice (Figure 7F). Hypoxanthine acts an endogenous monoamine oxidase (MAO) inhibitor. The serum level of hypoxanthine in the patients with IS was significantly decreased (Wang D. et al., 2017). Adenosine is a key nucleoside and regulates energy homeostasis to affect the function of the brain. Under IS condition, GA and CA entered the brain (Figure 7F). GA, a well-known antioxidant compounds, has neuroprotective actions in different models of neurodegeneration, neurotoxicity, and oxidative stress (Daglia et al., 2014). CA could attenuate LPS-induced inflammatory symptoms and oxidative stress (Zhao et al., 2020). Additionally, it was found that linoleic acid, astragalin, hyperoside and (-)-Catechin Gallate (CG) detected only in the brains of IS mice (Figure 7G). Linoleic acid, an essential polyunsaturated fatty acid, is required for normal growth and development (Taha, 2020). Astragalin (Riaz et al., 2018), hyperoside (Chen et al., 2021) and CG (Lee et al., 2020) are well known for their activities such as anti-inflammation, anti-oxidant and neuroprotection. From above data, anti-oxidant and anti-inflammation may be two important aspects of CR’s cerebral protections in IS, which needs to be further verified.
Discussion
In this study we used a data processing method that combines network pharmacology analysis with the important pathophysiological processes of IS, and uncovered the key phytochemicals and biological functions of CR and its underlying mechanisms in the treatment of IS. A total of 18 bioactive phytochemicals of CR corresponding to 85 anti-IS targets were obtained, and coniferyl ferulate, neocnidilide and ferulic acid were identified as the key phytochemicals of CR against IS. CR exerts protective effects on neurons mainly by acting on targets related to synaptic structure, synaptic function, neuronal survival and neuronal growth. CR also has the activities of anti-inflammation, anti-oxidative stress, anti-cell death and improves blood circulation. Infection and an abrupt increase in blood pressure are common complications of IS. In this research, we found that infection prevention and blood pressure regulation are the two most important synergistic effects of CR phytochemicals in treating IS.
Although network pharmacology approach contributes to collect active phytochemicals and pharmacological actions of herbs or Chinese medicine prescriptions, this approach also has many limitations. One important issue is the reliability of the data source. Normally, vast majority of network pharmacology studies used public databases to screen the active phytochemicals in traditional Chinese medicine. This would result in a large number of non-specific phytochemicals to be screened out. For example, seven active phytochemicals in CR have an oral bioavailability (OB) ≥ 30% and drug likeness (DL) ≥ 0.18 according to the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP) (Ru et al., 2014). However, this screening method does not consider the content of each component. Among these phytochemicals, sitosterol exists widely in 179 kinds of herbs and mandenol exists in 35 kinds of herbs (data from HERB, http://herb.ac.cn/, date of search: October 20, 2021) (Fang et al., 2021). Obviously, compounds that are widely distributed in a variety of herbs and have good pharmacological properties are more likely to be identified as active phytochemicals. Unlike the conventional approach, the present study obtained data on active CR component, including aromatic acids, phthalides and tetramethylpyrazine from recent UHPLC-MS/MS-based studies (Yin et al., 2019; Wang et al., 2020). In stark contrast, butylidenephthalide, neocnidilide, senkyunolide A and z-ligustilide are only present in 3, 6, 10, and 13 herbs, respectively, including CR (data from HERB, date of search: October 20, 2021). There are differences in the phytochemicals of the same Chinese medicine with different harvest time and producing area. The contents of CR phytochemicals were analyzed using 36 CR samples from six production origins. It is necessary to design an algorithm to quantitatively evaluate the impact of component content on drug efficacy in a future study. Thus, the reliability of the source of CR component data and the follow-up analysis are significantly increased in this study. However, identification of key phytochemicals only through network pharmacological analysis has certain limitations, such as the complex dose-effect relationship and toxicological evaluation was not considered. Notably, most research based on network pharmacology is still static network analysis, which contradicts the dynamic changes of disease occurrence and development. Moreover, chemical changes that occur during various drug manufacturing processes are not considered in this study.
Among the phytochemicals of CR, 18 active phytochemicals meet with RO5 (Lipinski, 2004) and correspond to 376 targets. The intersection of drugs and disease targets can lead to inflated pathway enrichment analysis. Therefore, before analyzing the anti-IS effects of CR, it was necessary to obtain a more comprehensive understanding of the biological functions of CR-related targets. Unexpectedly, most CR targets are synaptic membrane proteins and involve multiple synapses. Synapses, as an essential unit for neuronal communication in the brain, and synaptic plasticity can accelerate functional recovery after IS (Clarkson et al., 2010; Clarkson et al., 2011). Neurons in the center of the ischemic region undergo liquefactive necrosis, during which the body and axons of neurons disappear (Park et al., 2013; Nguyen et al., 2016). IS can lead to neurological deficits (Kuklina et al., 2012; Esenwa and Gutierrez, 2015), and significantly enriched neuroactive ligand–receptor interaction signaling pathway (hsa04080) is directly related to neurological function. Mounting evidence suggests that the cAMP signaling pathway plays important biological roles in IS (Xin et al., 2020). Astaxanthin activates the cAMP/PKA signaling pathway to promote axonal regeneration around the infarction and improve motor function (Wang et al., 2019). Apoptosis is an important cause of neuronal death after IS, and hypoxia and calcium overload can induce cell apoptosis (Radak et al., 2017). KEGG pathway enrichment analysis revealed that CR can affect the apoptosis pathway by acting on the calcium signaling pathway and HIF-1 signaling pathway. In addition to exerting neuroprotective effects, CR affects multiple signaling pathways and structures associated with cerebral circulation, such as the HIF-1 signaling pathway, the VEGF signaling pathway and gap junctions.
The CR main component target-IS target network indicated that coniferyl ferulate, neocnidilide, ferulic acid, caffeic acid, levistilide A, senkyunolide A and z-ligustilide are likely to become key phytochemicals for CR against IS. Coniferyl ferulate possesses multiple pharmacological activities such as antibacterial (Chou et al., 2006), antioxidant (Li et al., 2007) and vasodilating effects (Dai et al., 2017). Treatment with coniferyl ferulate (50 or 100 mg/ kg for 14 days) exerts antidepressant effect by acting on NMDAR-CaMKII-MAPKs signaling pathway, oxidative stress, and mitochondrial apoptotic pathways (Gong et al., 2020). The content of z-ligustilide in CR was the highest and has vasodilatation (Cao et al., 2006). Neocindilide exhibits anti-inflammatory and antioxidant activities. The role of ferulic acid included anti-thrombosis (Hsieh et al., 2002). 5-lipoxygenase (ALOX5) is a key enzyme metabolizing arachidonic acid to produce leukotrienes. Caffeic acid is not only a specific inhibitor of ALOX5 but also possesses antioxidant and anti-inflammatory properties (Habtemariam, 2017). An animal study demonstrated that caffeic acid (30 or 50 mg/ kg) ameliorates global cerebral ischemia-reperfusion injury in rats by inhibiting ALOX5 (Liang et al., 2015). In addition, caffeic acid can inhibit platelet aggregation induced by collagen (Hsiao et al., 2007; Park, 2009). Currently, levistilide A is a potential P-gp modulator and used for treating cancer. Senkyunolide A and z-ligustilide are used in treating inflammation associated with cerebrovascular diseases (Or et al., 2011). In this study L-tryptophan and vanillin were detected in brains of Sham and IS mice, while linoleic acid, astragalin, hyperoside and CG detected only in the brains of IS mice. L-tryptophan, a precursor of serotonin, is used to improve the mood of healthy individuals (Kikuchi et al., 2021) and depressed individuals (Alvarez-Mon et al., 2021). The significantly decreased L-tryptophan in the brain was reported in acute ischemic stroke (Ormstad et al., 2013). Vanillin is reported as an anti-oxidative, anti-apoptotic, anti-inflammatory, and neuroprotective compound (Bezerra-Filho et al., 2019). Furthermore, vanillin has the ability of maintaining the integrity of the BBB (Lan et al., 2019). Linoleic acid was reported significantly increase the axonal outgrowth (Hennebelle et al., 2020). An inverse association between linoleic acid levels in circulation (Ye et al., 2020) and adipose tissue (Veno et al., 2018) with the risk of IS was reported. Astragalin (Riaz et al., 2018), hyperoside (Chen et al., 2021) and CG (Lee et al., 2020) are well known for their activities such as anti-inflammation, anti-oxidant and neuroprotection. Although these compounds have been detected in the mouse brain, we do not know how these components were produced from CR in the mouse body.
Inflammatory response is a key component in the pathophysiology of IS (Chen et al., 2017). Cluster one was identified to be involved in the inflammatory response, and its targets CXCL8, IL1B, MMP9 and PTGS2 are the core targets of CR in the treatment of IS. These core targets are involved in the NF-κB signaling pathway and interleukin-17 (IL-17) signaling pathway. The proinflammatory chemokine CXCL8, also known as interleukin-8 (IL-8), is involved in acute inflammation. Animal experiments have demonstrated that functional suppression of CXCL8 promotes neuroglial activation and inhibits neuroinflammation by regulating the PI3K/Akt/NF-κB signaling pathway (Lv et al., 2019). Functional analysis with DNA microarray also revealed that CXCL8 and TNF are the feature genes of IS and may be valid targets for IS treatment (Zhang et al., 2014). Moreover, the serum CXCL8 level may be a prognostic indicator for IS (He et al., 2018). MMP9 and MMP2 (belong to gelatinase) are activated during IS and can activate numerous proinflammatory agents as chemokine CXCL8, IL1B or TNFα (Kurzepa et al., 2014). As an important prognostic factor for IS, elevated levels of serum MMP9 in the acute phase of IS are associated with increased risk of mortality and major disability (Zhong et al., 2017). Many studies have shown that the level of matrix metalloproteinase (especially MMP9) increases after stroke, and is associated with BBB disruption (Turner and Sharp, 2016). A clinical study showed that the increased serum MMP9 level in acute phase of IS was associated with 3-months cognitive impairment (Zhong et al., 2018). Selective inhibition of gelatinase can be used to treat IS. Administration of MMPs inhibitors (BB-94) before rtPA treatment reduced mortality, suggesting that blocking MMPs activity reduces the risk associated with thrombolysis (Pfefferkorn and Rosenberg, 2003). In the acute phase of IS, IL1B and TNFα promote inflammatory injury and induce cell necrosis or apoptosis (Bi et al., 2015). Reduction of serum IL1B ultimately improve the clinical outcome of patients with IS (Kadri et al., 2020). A previous study showed that PTGS2 (also known as cyclooxygenase-2) produced prostaglandins and thromboxanes, which are important mediators of IS-induced inflammatory cascade (Dong et al., 2019). Protein expression of PTGS2 increased significantly after IS. The specific knockdown of PTGS2 can inhibit NF-κB signaling pathway, thus inhibit apoptosis, promote the proliferation, migration and angiogenesis of endothelial progenitor cells, and protect ischemic stroke mice (Zhou et al., 2019). Furthermore, the inflammatory response after IS can also lead to neuron apoptosis, forming a vicious circle, which in turn aggravates the inflammatory response (Mergenthaler et al., 2004). Overall, our findings indicate that CR can act on the above inflammatory response related core targets in the treatment of IS.
Conclusions
This study provides a preliminary exploration of the therapeutic effect of CR on ischemic brain injuries and its possible phytochemicals, and coniferyl ferulate, neocnidilide and ferulic acid are identified the key phytochemicals of CR against IS. Its brain protective effects involve anti-inflammation, anti-oxidant, anti-cell death and improving blood circulation. Simultaneously, preventing infection and regulating blood pressure are the two important manifestations of the synergistic effects of CR in treating IS. However, subsequent in vivo (oxygen-glucose deprivation model) and in vitro experiments need to be further proved. In addition, toxicity and safety studies should be performed before clinical translation. In summary, this study provides a pathophysiologically relevant pharmacological basis for further research on IS.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The animal study was reviewed and approved by The Animal Care and Use Committee of Huazhong University of Science and Technology.
Author Contributions
PZ conceived this research. PZ and H-FS contributed to the acquisition, analysis, and interpretation of data. C-YY and Y-WS helped in the data acquisition. PZ drafted this manuscript. YY performed the animal experiment. X-WZ and YL gave important suggestions on the writing and revised the draft. AS and QT contributed to the data analysis strategy and improving the quality of this paper. All authors agree to be accountable for all aspects of work ensuring integrity and accuracy.
Funding
This work was supported by the grants from the National Natural Science Foundation of China (82041478) and the Health Committee of Hubei Province (ZY 2021M041).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2021.758049/full#supplementary-material
References
Alvarez-Mon, M. A., Ortega, M. A., García-Montero, C., Fraile-Martinez, O., Monserrat, J., Lahera, G., et al. (2021). Exploring the Role of Nutraceuticals in Major Depressive Disorder (MDD): Rationale, State of the Art and Future Prospects. Pharmaceuticals 14 (8), 821. doi:10.3390/ph14080821
Bader, G. D., and Hogue, C. W. (2003). An Automated Method for Finding Molecular Complexes in Large Protein Interaction Networks. BMC Bioinformatics 4, 2. doi:10.1186/1471-2105-4-2
Benjamin, E. J., Virani, S. S., Callaway, C. W., Chamberlain, A. M., Chang, A. R., Cheng, S., et al. (2018). Heart Disease and Stroke Statistics-2018 Update: A Report from the American Heart Association. Circulation 137 (12), e67–e492. doi:10.1161/CIR.0000000000000558
Bezerra-Filho, C. S. M., Barboza, J. N., Souza, M. T. S., Sabry, P., Ismail, N. S. M., and de Sousa, D. P. (2019). Therapeutic Potential of Vanillin and its Main Metabolites to Regulate the Inflammatory Response and Oxidative Stress. Mini Rev. Med. Chem. 19 (20), 1681–1693. doi:10.2174/1389557519666190312164355
Bi, B. L., Wang, H. J., Bian, H., and Tian, Z. T. (2015). Identification of Therapeutic Targets of Ischemic Stroke with DNA Microarray. Eur. Rev. Med. Pharmacol. Sci. 19 (21), 4012–4019.
Cao, Y. X., Zhang, W., He, J. Y., He, L. C., and Xu, C. B. (2006). Ligustilide Induces Vasodilatation via Inhibiting Voltage Dependent Calcium Channel and Receptor-Mediated Ca2+ Influx and Release. Vascul Pharmacol. 45 (3), 171–176. doi:10.1016/j.vph.2006.05.004
Chen, S., Dong, Z., Cheng, M., Zhao, Y., Wang, M., Sai, N., et al. (2017). Homocysteine Exaggerates Microglia Activation and Neuroinflammation through Microglia Localized STAT3 Overactivation Following Ischemic Stroke. J. Neuroinflammation 14 (1), 187. doi:10.1186/s12974-017-0963-x
Chen, Z., Zhang, C., Gao, F., Fu, Q., Fu, C., He, Y., et al. (2018). A Systematic Review on the Rhizome of Ligusticum Chuanxiong Hort. (Chuanxiong). Food Chem. Toxicol. 119, 309–325. doi:10.1016/j.fct.2018.02.050
Chen, X., Famurewa, A. C., Tang, J., Olatunde, O. O., and Olatunji, O. J. (2021). Hyperoside Attenuates Neuroinflammation, Cognitive Impairment and Oxidative Stress via Suppressing TNF-α/NF-κB/caspase-3 Signaling in Type 2 Diabetes Rats. Nutr. Neurosci. [Epub ahead of print], 1–11. doi:10.1080/1028415X.2021.1901047
Chou, S. C., Everngam, M. C., Sturtz, G., and Beck, J. J. (2006). Antibacterial Activity of Components from Lomatium Californicum. Phytother. Res. 20 (2), 153–156. doi:10.1002/ptr.1828
Clarkson, A. N., Huang, B. S., Macisaac, S. E., Mody, I., and Carmichael, S. T. (2010). Reducing Excessive GABA-Mediated Tonic Inhibition Promotes Functional Recovery after Stroke. Nature 468 (7321), 305–309. doi:10.1038/nature09511
Clarkson, A. N., Overman, J. J., Zhong, S., Mueller, R., Lynch, G., and Carmichael, S. T. (2011). AMPA Receptor-Induced Local Brain-Derived Neurotrophic Factor Signaling Mediates Motor Recovery after Stroke. J. Neurosci. 31 (10), 3766–3775. doi:10.1523/JNEUROSCI.5780-10.2011
Daglia, M., Di Lorenzo, A., Nabavi, S. F., Talas, Z. S., and Nabavi, S. M. (2014). Polyphenols: Well beyond the Antioxidant Capacity: Gallic Acid and Related Compounds as Neuroprotective Agents: You Are what You Eat!. Curr. Pharm. Biotechnol. 15 (4), 362–372. doi:10.2174/138920101504140825120737
Dai, X., Pang, L., Zhang, Z., Yang, C., and Li, Y. (2017). Development of a Sensitive LC-MS/MS Method for Quantification of Coniferyl Ferulate and its Metabolite Coniferyl Alcohol in Rat Plasma: Application to a Pharmacokinetic Study. J. Pharm. Biomed. Anal. 146, 201–205. doi:10.1016/j.jpba.2017.08.033
Daina, A., Michielin, O., and Zoete, V. (2017). SwissADME: a Free Web Tool to Evaluate Pharmacokinetics, Drug-Likeness and Medicinal Chemistry Friendliness of Small Molecules. Sci. Rep. 7, 42717. doi:10.1038/srep42717
Daina, A., Michielin, O., and Zoete, V. (2019). SwissTargetPrediction: Updated Data and New Features for Efficient Prediction of Protein Targets of Small Molecules. Nucleic Acids Res. 47 (W1), W357–W364. doi:10.1093/nar/gkz382
Dong, H., Cui, B., and Hao, X. (2019). MicroRNA22 Alleviates Inflammation in Ischemic Stroke via P38 MAPK Pathways. Mol. Med. Rep. 20 (1), 735–744. doi:10.3892/mmr.2019.10269
Donnan, G. A., Fisher, M., Macleod, M., and Davis, S. M. (2008). Stroke. Lancet 371 (9624), 1612–1623. doi:10.1016/S0140-6736(08)60694-7
Esenwa, C., and Gutierrez, J. (2015). Secondary Stroke Prevention: Challenges and Solutions. Vasc. Health Risk Manag. 11, 437–450. doi:10.2147/VHRM.S63791
Fang, S., Dong, L., Liu, L., Guo, J., Zhao, L., Zhang, J., et al. (2021). HERB: a High-Throughput experiment- and Reference-Guided Database of Traditional Chinese Medicine. Nucleic Acids Res. 49 (D1), D1197–D1206. doi:10.1093/nar/gkaa1063
Gaberel, T., Gakuba, C., Zheng, Y., Lépine, M., Lo, E. H., and van Leyen, K. (2019). Impact of 12/15-Lipoxygenase on Brain Injury after Subarachnoid Hemorrhage. Stroke 50 (2), 520–523. doi:10.1161/STROKEAHA.118.022325
Gong, W., Zhou, Y., Gong, W., and Qin, X. (2020). Coniferyl Ferulate Exerts Antidepressant Effect via Inhibiting the Activation of NMDAR-CaMKII-MAPKs and Mitochondrial Apoptotic Pathways. J. Ethnopharmacol. 251, 112533. doi:10.1016/j.jep.2019.112533
Gouveia, A., Seegobin, M., Kannangara, T. S., He, L., Wondisford, F., Comin, C. H., et al. (2017). The aPKC-CBP Pathway Regulates Post-stroke Neurovascular Remodeling and Functional Recovery. Stem Cel Rep. 9 (6), 1735–1744. doi:10.1016/j.stemcr.2017.10.021
Grosdidier, A., Zoete, V., and Michielin, O. (2011). SwissDock, a Protein-Small Molecule Docking Web Service Based on EADock DSS. Nucleic Acids Res. 39 (Web Server issue), W270–W277. doi:10.1093/nar/gkr366
Habtemariam, S. (2017). Protective Effects of Caffeic Acid and the Alzheimer's Brain: An Update. Mini Rev. Med. Chem. 17 (8), 667–674. doi:10.2174/1389557516666161130100947
He, Q., Shi, X., Zhou, B., Teng, J., Zhang, C., Liu, S., et al. (2018). Interleukin 8 (CXCL8)-CXC Chemokine Receptor 2 (CXCR2) axis Contributes to MiR-4437-Associated Recruitment of Granulocytes and Natural Killer Cells in Ischemic Stroke. Mol. Immunol. 101, 440–449. doi:10.1016/j.molimm.2018.08.002
Hennebelle, M., Morgan, R. K., Sethi, S., Zhang, Z., Chen, H., Grodzki, A. C., et al. (2020). Linoleic Acid-Derived Metabolites Constitute the Majority of Oxylipins in the Rat Pup Brain and Stimulate Axonal Growth in Primary Rat Cortical Neuron-Glia Co-cultures in a Sex-dependent Manner. J. Neurochem. 152 (2), 195–207. doi:10.1111/jnc.14818
Hopkins, A. L. (2008). Network Pharmacology: the Next Paradigm in Drug Discovery. Nat. Chem. Biol. 4 (11), 682–690. doi:10.1038/nchembio.118
Hsiao, G., Lee, J. J., Lin, K. H., Shen, C. H., Fong, T. H., Chou, D. S., et al. (2007). Characterization of a Novel and Potent Collagen Antagonist, Caffeic Acid Phenethyl Ester, in Human Platelets: In Vitro and In Vivo Studies. Cardiovasc. Res. 75 (4), 782–792. doi:10.1016/j.cardiores.2007.05.005
Hsieh, M. T., Tsai, F. H., Lin, Y. C., Wang, W. H., and Wu, C. R. (2002). Effects of Ferulic Acid on the Impairment of Inhibitory Avoidance Performance in Rats. Planta Med. 68 (8), 754–756. doi:10.1055/s-2002-33800
Kadri, A., Sjahrir, H., Sembiring, R. J., and Ichwan, M. (2020). Combination of Vitamin A and D Supplementation for Ischemic Stroke: Effects on Interleukin-1ss and Clinical Outcome. Med. Glas (Zenica) 17 (2), 425–432. doi:10.17392/1137-20
Kidani, N., Hishikawa, T., Hiramatsu, M., Nishihiro, S., Kin, K., Takahashi, Y., et al. (2020). Cerebellar Blood Flow and Gene Expression in Crossed Cerebellar Diaschisis after Transient Middle Cerebral Artery Occlusion in Rats. Int. J. Mol. Sci. 21 (11), 4137. doi:10.3390/ijms21114137
Kikuchi, A. M., Tanabe, A., and Iwahori, Y. (2021). A Systematic Review of the Effect of L-Tryptophan Supplementation on Mood and Emotional Functioning. J. Diet. Suppl. 18 (3), 316–333. doi:10.1080/19390211.2020.1746725
Kim, S., Thiessen, P. A., Cheng, T., Yu, B., Shoemaker, B. A., Wang, J., et al. (2016). Literature Information in PubChem: Associations between PubChem Records and Scientific Articles. J. Cheminform 8, 32. doi:10.1186/s13321-016-0142-6
Kuklina, E. V., Tong, X., George, M. G., and Bansil, P. (2012). Epidemiology and Prevention of Stroke: a Worldwide Perspective. Expert Rev. Neurother. 12 (2), 199–208. doi:10.1586/ern.11.99
Kurzepa, J., Kurzepa, J., Golab, P., Czerska, S., and Bielewicz, J. (2014). The Significance of Matrix Metalloproteinase (MMP)-2 and MMP-9 in the Ischemic Stroke. Int. J. Neurosci. 124 (10), 707–716. doi:10.3109/00207454.2013.872102
Lan, X. B., Wang, Q., Yang, J. M., Ma, L., Zhang, W. J., Zheng, P., et al. (2019). Neuroprotective Effect of Vanillin on Hypoxic-Ischemic Brain Damage in Neonatal Rats. Biomed. Pharmacother. 118, 109196. doi:10.1016/j.biopha.2019.109196
Lee, S. B., Choi, E. H., Jeong, K. H., Kim, K. S., Shim, S. M., and Kim, G. H. (2020). Effect of Catechins and High-Temperature-Processed green tea Extract on Scavenging Reactive Oxygen Species and Preventing Aβ1-42 Fibrils' Formation in Brain Microvascular Endothelium. Nutr. Neurosci. 23 (5), 363–373. doi:10.1080/1028415X.2018.1507618
Li, S. Y., Yu, Y., and Li, S. P. (2007). Identification of Antioxidants in Essential Oil of Radix Angelicae Sinensis Using HPLC Coupled with DAD-MS and ABTS-Based Assay. J. Agric. Food Chem. 55 (9), 3358–3362. doi:10.1021/jf070140t
Li, J., Yu, J., Ma, H., Yang, N., Li, L., Zheng, D. D., et al. (2017). Intranasal Pretreatment with Z-Ligustilide, the Main Volatile Component of Rhizoma Chuanxiong, Confers Prophylaxis against Cerebral Ischemia via Nrf2 and HSP70 Signaling Pathways. J. Agric. Food Chem. 65 (8), 1533–1542. doi:10.1021/acs.jafc.6b04979
Liang, G., Shi, B., Luo, W., and Yang, J. (2015). The Protective Effect of Caffeic Acid on Global Cerebral Ischemia-Reperfusion Injury in Rats. Behav. Brain Funct. 11, 18. doi:10.1186/s12993-015-0064-x
Lin, J., Hao, C., Gong, Y., Zhang, Y., Li, Y., Feng, Z., et al. (2021). Effect of Tetramethylpyrazine on Neuroplasticity after Transient Focal Cerebral Ischemia Reperfusion in Rats. Evid. Based Complement. Alternat Med. 2021, 1587241. doi:10.1155/2021/1587241
Lipinski, C. A., Lombardo, F., Dominy, B. W., and Feeney, P. J. (2001). Experimental and Computational Approaches to Estimate Solubility and Permeability in Drug Discovery and Development Settings. Adv. Drug Deliv. Rev. 46 (1-3), 3–26. doi:10.1016/s0169-409x(00)00129-0
Lipinski, C. A. (2004). Lead- and Drug-like Compounds: the Rule-Of-Five Revolution. Drug Discov. Today Technol. 1 (4), 337–341. doi:10.1016/j.ddtec.2004.11.007
Liu, J. Q., Dai, S. X., Zheng, J. J., Guo, Y. C., Li, W. X., Li, G. H., et al. (2017). The Identification and Molecular Mechanism of Anti-stroke Traditional Chinese Medicinal Compounds. Sci. Rep. 7, 41406. doi:10.1038/srep41406
Liu, S., Luo, W., and Wang, Y. (2021). Emerging Role of PARP‐1 and PARthanatos in Ischemic Stroke. J. Neurochem. doi:10.1111/jnc.15464
Lv, H., Li, J., and Che, Y. Q. (2019). CXCL8 Gene Silencing Promotes Neuroglial Cells Activation while Inhibiting Neuroinflammation through the PI3K/Akt/NF-κB-Signaling Pathway in Mice with Ischemic Stroke. J. Cel. Physiol. 234 (5), 7341–7355. doi:10.1002/jcp.27493
Mandloi, S., and Chakrabarti, S. (2015). PALM-IST: Pathway Assembly from Literature Mining-Aan Information Search Tool. Sci. Rep. 5, 10021. doi:10.1038/srep10021
Mergenthaler, P., Dirnagl, U., and Meisel, A. (2004). Pathophysiology of Stroke: Lessons from Animal Models. Metab. Brain Dis. 19 (3-4), 151–167. doi:10.1023/b:mebr.0000043966.46964.e6
Mi, H., Guo, N., Kejariwal, A., and Thomas, P. D. (2007). PANTHER Version 6: Protein Sequence and Function Evolution Data with Expanded Representation of Biological Pathways. Nucleic Acids Res. 35 (Database issue), D247–D252. doi:10.1093/nar/gkl869
Nguyen, T. V., Frye, J. B., Zbesko, J. C., Stepanovic, K., Hayes, M., Urzua, A., et al. (2016). Erratum to: Multiplex Immunoassay Characterization and Species Comparison of Inflammation in Acute and Non-acute Ischemic Infarcts in Human and Mouse Brain Tissue. Acta Neuropathol. Commun. 4 (1), 104. doi:10.1186/s40478-016-0371-y10.1186/s40478-016-0376-6
Or, T. C., Yang, C. L., Law, A. H., Li, J. C., and Lau, A. S. (2011). Isolation and Identification of Anti-inflammatory Constituents from Ligusticum Chuanxiong and Their Underlying Mechanisms of Action on Microglia. Neuropharmacology 60 (6), 823–831. doi:10.1016/j.neuropharm.2010.12.002
Ormstad, H., Verkerk, R., Aass, H. C., Amthor, K. F., and Sandvik, L. (2013). Inflammation-induced Catabolism of Tryptophan and Tyrosine in Acute Ischemic Stroke. J. Mol. Neurosci. 51 (3), 893–902. doi:10.1007/s12031-013-0097-2
Park, J. Y., Byeon, J. H., Park, S. W., Eun, S. H., Chae, K. Y., and Eun, B. L. (2013). Neuroprotective Effect of Human Placental Extract on Hypoxic-Ischemic Brain Injury in Neonatal Rats. Brain Dev. 35 (1), 68–74. doi:10.1016/j.braindev.2012.01.009
Park, J. B. (2009). 5-Caffeoylquinic Acid and Caffeic Acid Orally Administered Suppress P-Selectin Expression on Mouse Platelets. J. Nutr. Biochem. 20 (10), 800–805. doi:10.1016/j.jnutbio.2008.07.009
Pfefferkorn, T., and Rosenberg, G. A. (2003). Closure of the Blood-Brain Barrier by Matrix Metalloproteinase Inhibition Reduces rtPA-Mediated Mortality in Cerebral Ischemia with Delayed Reperfusion. Stroke 34 (8), 2025–2030. doi:10.1161/01.STR.0000083051.93319.28
Radak, D., Katsiki, N., Resanovic, I., Jovanovic, A., Sudar-Milovanovic, E., Zafirovic, S., et al. (2017). Apoptosis and Acute Brain Ischemia in Ischemic Stroke. Curr. Vasc. Pharmacol. 15 (2), 115–122. doi:10.2174/1570161115666161104095522
Riaz, A., Rasul, A., Hussain, G., Zahoor, M. K., Jabeen, F., Subhani, Z., et al. (2018). Astragalin: A Bioactive Phytochemical with Potential Therapeutic Activities. Adv. Pharmacol. Sci. 2018, 9794625. doi:10.1155/2018/9794625
Ru, J., Li, P., Wang, J., Zhou, W., Li, B., Huang, C., et al. (2014). TCMSP: a Database of Systems Pharmacology for Drug Discovery from Herbal Medicines. J. Cheminform 6, 13. doi:10.1186/1758-2946-6-13
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 13 (11), 2498–2504. doi:10.1101/gr.1239303
Stelzer, G., Rosen, N., Plaschkes, I., Zimmerman, S., Twik, M., Fishilevich, S., et al. (2016). The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses. Curr. Protoc. Bioinf. 54, 1.30.1–1.30.33. doi:10.1002/cpbi.5
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING V11: Protein-Protein Association Networks with Increased Coverage, Supporting Functional Discovery in Genome-wide Experimental Datasets. Nucleic Acids Res. 47 (D1), D607–D613. doi:10.1093/nar/gky1131
Taha, A. Y. (2020). Linoleic Acid-Good or Bad for the Brain? NPJ Sci. Food 4, 1. doi:10.1038/s41538-019-0061-9
Turner, R. J., and Sharp, F. R. (2016). Implications of MMP9 for Blood Brain Barrier Disruption and Hemorrhagic Transformation Following Ischemic Stroke. Front. Cel. Neurosci. 10, 56. doi:10.3389/fncel.2016.00056
van der Worp, H. B., and van Gijn, J. (2007). Clinical Practice. Acute Ischemic Stroke. N. Engl. J. Med. 357 (6), 572–579. doi:10.1056/NEJMcp072057
Venø, S. K., Bork, C. S., Jakobsen, M. U., Lundbye-Christensen, S., Bach, F. W., Overvad, K., et al. (2018). Linoleic Acid in Adipose Tissue and Development of Ischemic Stroke: A Danish Case-Cohort Study. J. Am. Heart Assoc. 7 (13), e009820. doi:10.1161/JAHA.118.009820
Wang, D., Kong, J., Wu, J., Wang, X., and Lai, M. (2017). GC-MS-based Metabolomics Identifies an Amino Acid Signature of Acute Ischemic Stroke. Neurosci. Lett. 642, 7–13. doi:10.1016/j.neulet.2017.01.039
Wang, Y. L., Zhu, X. L., Sun, M. H., and Dang, Y. K. (2019). Effects of Astaxanthin Onaxonal Regeneration via cAMP/PKA Signaling Pathway in Mice with Focal Cerebral Infarction. Eur. Rev. Med. Pharmacol. Sci. 23 (3 Suppl. l), 135–143. doi:10.26355/eurrev_201908_18640
Wang, X., Yao, Y., An, C., Li, X., Xiang, F., Dong, Y., et al. (2020). Simultaneous Determination of 20 Bioactive Components in Chuanxiong Rhizoma from Different Production Origins in Sichuan Province by Ultra-high-performance Liquid Chromatography Coupled with Triple Quadrupole Mass Spectrometry Combined with Multivariate Statistical Analysis. Electrophoresis 41 (18-19), 1606–1616. doi:10.1002/elps.202000082
Wang, W., Jiang, B., Sun, H., Ru, X., Sun, D., Wang, L., et al. (2017). Prevalence, Incidence, and Mortality of Stroke in China: Results from a Nationwide Population-Based Survey of 480 687 Adults. Circulation 135 (8), 759–771. doi:10.1161/CIRCULATIONAHA.116.025250
Westbrook, J., Feng, Z., Jain, S., Bhat, T. N., Thanki, N., Ravichandran, V., et al. (2002). The Protein Data Bank: Unifying the Archive. Nucleic Acids Res. 30 (1), 245–248. doi:10.1093/nar/30.1.245
Xin, M., Feng, J., Hao, Y., You, J., Wang, X., Yin, X., et al. (2020). Cyclic Adenosine Monophosphate in Acute Ischemic Stroke: Some to Update, More to Explore. J. Neurol. Sci. 413, 116775. doi:10.1016/j.jns.2020.116775
Ye, D., Huang, H., Wu, D. J. H., Zhang, W., Zhou, F., Qian, Y., et al. (2020). Association between Circulating Linoleic Acid and Risk of Ischemic Stroke. Front. Genet. 11, 582623. doi:10.3389/fgene.2020.582623
Yin, D. D., Yang, M., Wang, Y. L., Yin, D. K., Liu, H. K., Zhou, M., et al. (2018). High Tetramethylpyrazine Production by the Endophytic Bacterial Bacillus Subtilis Isolated from the Traditional Medicinal Plant Ligusticum Chuanxiong Hort. AMB Express 8 (1), 193. doi:10.1186/s13568-018-0721-1
Yin, D. D., Wang, Y. L., Yang, M., Yin, D. K., Wang, G. K., and Xu, F. (2019). Analysis of Chuanxiong Rhizoma Substrate on Production of Ligustrazine in Endophytic Bacillus Subtilis by Ultra High Performance Liquid Chromatography with Quadrupole Time-Of-Flight Mass Spectrometry. J. Sep. Sci. 42 (19), 3067–3076. doi:10.1002/jssc.201900030
Yu, G., Wang, L. G., Han, Y., and He, Q. Y. (2012). clusterProfiler: an R Package for Comparing Biological Themes Among Gene Clusters. OMICS 16 (5), 284–287. doi:10.1089/omi.2011.0118
Yuan, Y., Zeng, X., Luo, Y., Li, Z., and Wu, T. (2008). Chuanxiong-type Preparations for Acute Ischemic Stroke. Cochrane Database Syst. Rev. 4, CD005569. doi:10.1002/14651858.CD005569.pub2
Zeng, P., Shi, Y., Wang, X. M., Lin, L., Du, Y. J., Tang, N., et al. (2019). Emodin Rescued Hyperhomocysteinemia-Induced Dementia and Alzheimer's Disease-Like Features in Rats. Int. J. Neuropsychopharmacol. 22 (1), 57–70. doi:10.1093/ijnp/pyy090
Zeng, P., Fang, M., Zhao, H., and Guo, J. (2021a). A Network Pharmacology Approach to Uncover the Key Ingredients in Ginkgo Folium and Their Anti-alzheimer's Disease Mechanisms. Aging (Albany NY) 13 (14), 18993–19012. doi:10.18632/aging.203348
Zeng, P., Wang, X. M., Ye, C. Y., Su, H. F., Fang, Y. Y., Zhang, T., et al. (2021b). Mechanistic Insights into the Anti-depressant Effect of Emodin: an Integrated Systems Pharmacology Study and Experimental Validation. Aging (Albany NY) 13 (11), 15078–15099. doi:10.18632/aging.203072
Zeng, P., Wang, X. M., Ye, C. Y., Su, H. F., and Tian, Q. (2021c). The Main Alkaloids in Uncaria Rhynchophylla and Their Anti-Alzheimer's Disease Mechanism Determined by a Network Pharmacology Approach. Int. J. Mol. Sci. 22 (7), 3612. doi:10.3390/ijms22073612
Zhang, Z. L., Wu, W. C., Liu, J. Q., Yao, Y. B., Pan, M. D., Yang, C. B., et al. (2014). Screening of Differentially Expressed Genes Related to Ischemic Stroke and Functional Analysis with DNA Microarray. Eur. Rev. Med. Pharmacol. Sci. 18 (8), 1181–1188.
Zhang, X., Ding, J., Gou, C., Wen, T., Li, L., Wang, X., et al. (2017). Qingchangligan Formula Attenuates the Inflammatory Response to Protect the Liver from Acute Failure Induced by D-Galactosamine/lipopolysaccharide in Mice. J. Ethnopharmacol. 201, 108–116. doi:10.1016/j.jep.2016.11.007
Zhang, X., Feng, Z. M., Yang, Y. N., Jiang, J. S., and Zhang, P. C. (2019). Bioactive Butylphthalide Derivatives from Ligusticum Chuanxiong. Bioorg. Chem. 84, 505–510. doi:10.1016/j.bioorg.2018.12.032
Zhao, X. L., Yu, L., Zhang, S. D., Ping, K., Ni, H. Y., Qin, X. Y., et al. (2020). Cryptochlorogenic Acid Attenuates LPS-Induced Inflammatory Response and Oxidative Stress via Upregulation of the Nrf2/HO-1 Signaling Pathway in RAW 264.7 Macrophages. Int. Immunopharmacol. 83, 106436. doi:10.1016/j.intimp.2020.106436
Zhong, C., Yang, J., Xu, T., Xu, T., Peng, Y., Wang, A., et al. (2017). Serum Matrix Metalloproteinase-9 Levels and Prognosis of Acute Ischemic Stroke. Neurology 89 (8), 805–812. doi:10.1212/WNL.0000000000004257
Zhong, C., Bu, X., Xu, T., Guo, L., Wang, X., Zhang, J., et al. (2018). Serum Matrix Metalloproteinase-9 and Cognitive Impairment after Acute Ischemic Stroke. J. Am. Heart Assoc. 7 (1), e007776. doi:10.1161/JAHA.117.007776
Zhou, Z., Lu, C., Meng, S., Dun, L., Yin, N., An, H., et al. (2019). Silencing of PTGS2 Exerts Promoting Effects on Angiogenesis Endothelial Progenitor Cells in Mice with Ischemic Stroke via Repression of the NF-κB Signaling Pathway. J. Cel. Physiol. 234 (12), 23448–23460. doi:10.1002/jcp.28914
Zhu, T., Wang, L., Feng, Y., Sun, G., and Sun, X. (2021). Classical Active Ingredients and Extracts of Chinese Herbal Medicines: Pharmacokinetics, Pharmacodynamics, and Molecular Mechanisms for Ischemic Stroke. Oxid. Med. Cel. Longev. 2021, 8868941. doi:10.1155/2021/8868941
Glossary
ALOX15 arachidonate 15-lipoxygenase
APP amyloid beta precursor protein
BBB blood-brain barrier
BP biological process
CC cellular component
CR Chuanxiong Rhizome
CREBBP CREB binding protein
DL drug likeness
EGFR epidermal growth factor receptor
GO Gene ontology
Hacc hydrogen bond acceptors
Hdon hydrogen bond donors
ICAs internal carotid arteries
IL-8 interleukin-8
IS Ischemic stroke
KEGG Kyoto Encyclopedia of Genes and Genomes
LogP lipid-water partition coefficient
MAO monoamine oxidase
MCAO middle cerebral artery occlusion
MCODE Molecular Complex Detection
MF molecular function
MW molecule weight
OB oral bioavailability
PALM-IST Pathway Assembly from Literature Mining-an Information Search Tool
PPI protein-protein interaction
Rbon rotatable bonds
RO5 Lipinski’s rule of five
ROS reactive oxygen species
TCMSP Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform
TTC Triphenyltetrazolium chloride
tPA tissue plasminogen activator
UHPLC-MS/MS ultra-HPLC with triple quadrupole MS
Keywords: Chuanxiong Rhizoma, ischemic stroke, network pharmacology, key phytochemicals, pathophysiology
Citation: Zeng P, Yi Y, Su H-F, Ye C-Y, Sun Y-W, Zhou X-W, Lu Y, Shi A and Tian Q (2021) Key Phytochemicals and Biological Functions of Chuanxiong Rhizoma Against Ischemic Stroke: A Network Pharmacology and Experimental Assessment. Front. Pharmacol. 12:758049. doi: 10.3389/fphar.2021.758049
Received: 16 August 2021; Accepted: 06 December 2021;
Published: 21 December 2021.
Edited by:
Zhong-Qiu Liu, Guangzhou University of Chinese Medicine, ChinaReviewed by:
Maoxiang Zhang, Guangzhou University of Chinese Medicine, ChinaPei-Wen Hsieh, Chang Gung University, Taiwan
Copyright © 2021 Zeng, Yi, Su, Ye, Sun, Zhou, Lu, Shi and Tian. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Anbing Shi, ashi@hust.edu.cn; Qing Tian, tianq@hust.edu.cn
†These authors have contributed equally to this work