- 1Department of Pharmaceutical Analysis, School of Pharmacy, Xi’an Jiaotong University, Xi’an, China
- 2Department of Pharmacy, Xijing Hospital, Air Force Medical University, Xi’an, China
- 3Precision Pharmacy and Drug Development Center, Department of Pharmacy, Tangdu Hospital, Air Force Medical University, Xi’an, China
- 4Department of Pharmaceutical Analysis, College of Pharmacy, Shenzhen Technology University, Shenzhen, China
Ischemic stroke (IS), predominantly triggered by blockages in cerebral blood flow, is increasingly recognized as a critical public health issue. The combination of Salvia miltiorrhiza (SM) and Cortex moutan (CM), traditional herbs in Eastern medicine, are frequently used for managing heart and brain vascular conditions. However, the exact mechanisms by which this herb pair (SC) combats IS remain largely unexplored. This investigation focuses on pinpointing the active constituents in SC that contribute to its protective role and deciphering the mechanisms countering cerebral ischemia, particularly in a middle cerebral artery occlusion (MCAO) rat model. We employed UPLC-Q-TOF-MS/MS alongside network pharmacology for predicting SC’s target actions against IS. Key ingredients were examined for their interaction with principal targets using molecular docking. The therapeutic impact was gauged through H&E, TUNEL, and Nissl staining, complemented by transcriptomic and metabolomic integration for mechanistic insights, with vital genes confirmed via western blot. UPLC-Q-TOF-MS/MS analysis revealed that the main components of SC included benzoylpaeoniflorin, salvianolic acid B, oxypaeoniflora, salvianolic acid A, and others. Network pharmacology analysis indicated that SC’s mechanism in treating IS primarily involves inflammation, angiogenesis, and cell apoptosis-related pathways, potentially through targets such as AKT1, TNF, PTGS2, MMP9, PIK3CA, and VEGFA. Molecular docking underscored strong affinities between these constituents and their targets. Our empirical studies indicated SC’s significant role in enhancing neuroprotection in IS, with transcriptomics suggesting the involvement of the VEGFA/PI3K/AKT pathway and metabolomics revealing improvements in various metabolic processes, including amino acids, glycerophospholipids, sphingomyelin, and fatty acids metabolisms.
1 Introduction
Ischemic stroke (IS) represents a critical neurological condition caused by insufficient cerebral blood supply, resulting in sensory and motor deficits, cognitive decline, and potentially fatal outcomes (Mendelson and Prabhakaran, 2021). It makes up approximately 85% of all strokes and stands as the third most significant contributor to worldwide mortality and disability, imposing a significant economic and public health load (Virani et al., 2020). Current IS treatments primarily involve t-PA thrombolytic therapy, but surgical thrombectomy can paradoxically worsen ischemia-reperfusion injury, exacerbating brain damage (Henderson et al., 2018). Therefore, the exploration of novel therapeutic strategies for stroke is imperative.
The Salvia miltiorrhiza Bunge (Salvia miltiorrhiza, SM) and Paeonia × suffruticosa Andrews (Cortex moutan, CM) herb pair (SC) is a traditional combination in Eastern medicine known for its efficacy for its blood-flow enhancing properties. It has been historically applied to a spectrum of ailments, including cardiovascular and cerebrovascular diseases, liver diseases, pneumonia, renal diseases, and even cancers (Sun et al., 2005; Lee et al., 2006; Yang et al., 2010; Li et al., 2022; Zhang M. et al., 2022). SM contains bioactive compounds such as tanshinones and phenolic acids, known for their diverse pharmacological effects, including antioxidant, endothelial protective, anti-atherosclerotic, anti-inflammatory, and vasodilatory properties. For instance, salvianolate has been verified to protect against cerebral ischemic injury in diabetic rats by regulating inflammatory factors and the antioxidant pathway (Wang et al., 2017). Additionally, SM ethanol extract can alleviate cerebral ischemic injury in rats, inhibit thrombosis formation, and attenuate platelet aggregation (Fei et al., 2017). CM, on the other hand, contains paeonol, a prominent component with potent vasodilatory effects, potentially reducing intracellular calcium concentration by modulating Ca2+ influx and release (Li et al., 2010). However, due to the complexity and cooperativity of SC, the molecular mechanisms underlying its neuroprotective effects remain elusive.
The Phosphoinositide 3-Kinase (PI3K)/AKT pathway, crucial in cell survival and metabolism, plays a pivotal role in IS pathology (Liu et al., 2023b). Hypoxic conditions and inflammatory cytokines can stimulate the release of growth factors, many of which activate the AKT signal transduction pathway (Koh and Lo, 2015). These growth factors promote neuroprotection by activating AKT and influencing processes such as inflammation, apoptosis, and angiogenesis. Gomisin N, isolated from Schisandra chinensis (Turcz.) Baill, has been shown to attenuate IS-induced neuroinflammation and autophagy by regulating the PI3K/AKT/mTOR signaling pathway (Li R. et al., 2023). Recent research suggests that PI3K/AKT signaling can modulate arterial thrombosis (Liberale et al., 2023) and attenuate blood-brain barrier disruption in the penumbra, consequently improving post-ischemia neurological function. Vascular endothelial growth factor (VEGF) is a well-known angiogenic factor, that exerts its angiogenesis promotion and neuroprotective effects by binding to the VEGF receptor. SM has been demonstrated to promote angiogenesis in ischemic heart disease (Li X. et al., 2023).
The aforementioned findings suggest that SC may have therapeutic potential for IS; however, no relevant reports have been published to date. In this research, we aim to verify the neuroprotective effects of SC via the VEGFA/PI3K/AKT signaling pathway. Firstly, we will identify the components of SC and establish MCAO rat models to explore the effects of SC on neuroprotective development. Additionally, we will employ a multidisciplinary approach, including systems pharmacology, transcriptomics, and metabolomics, to unravel the mechanisms through which SC acts against cerebral ischemic injury.
2 Materials and methods
2.1 Preparation of the SC herb pair decoction
Salvia miltiorrhiza (SM) and Cortex moutan (CM) were procured from Xijing Hospital (Xi’an, China). SC herb pair decoction was prepared according to the method of the Chinese Pharmacopoeia (2020 Edition). First, Cortex moutan (100 g) was powdered and extracted by adding six times the volume of water for distillation. Then, Salvia miltiorrhiza (200 g) and the dregs of Cortex moutan were extracted by adding six times the volume of water boiled twice for 2 h, filtered, mixed, and then combined with the Cortex moutan distillate. The resulting aqueous extract was concentrated to 2 g/mL for Salvia miltiorrhiza and 1 g/mL for Cortex moutan crude drug.
2.2 Determination of SC components
The SC concoction underwent vacuum freeze-drying and was then analyzed through ultra-high-performance liquid chromatography-quadrupole time-of-flight tandem mass spectrometry (UPLC-Q-TOF-MS/MS, Waters Acquity I-CLASSTM UPLC, Waters Xevo G2-XS, United States). The SC powder was dissolved in 30% methanol (2 mg/mL), agitated via vortex, and centrifuged at 14,000 rpm for 10 min. The supernatant was then filtered through a 0.22-μm filter for analysis. An ACQUITY HSS T3 column (2.1 × 100 mm, 1.8 μm, United States) at 40°C with a 0.4 mL/min flow rate was used for chromatographic separation. The gradient elution program utilized a mobile phase containing acetonitrile (A) and 0.1% formic acid water (B). The program was as follows: 0–2 min, 95% B; 2–14 min, 90%–65% B; 14–15 min, 65%–60% B; 15–25 min, 60%–20% B; 25–25.1 min, 20%–95% B; 25.1–30 min, 95% B. Mass spectrometry was performed in both negative and positive modes, with data spanning m/z 50–1,500 Da. Analysis was facilitated by Waters UNIFI software (version 1.9.1).
2.3 Experimental design and animal treatments
Sixty adult male Sprague-Dawley rats (250 ± 20 g) were obtained from the Air Force Military Medical University’s Experimental Animal Center (Xi’an, China, license number for animal production: SCXK-2019-001). Animals were housed under a 12-h light-dark cycle at 23°C ± 2°C and 50%–70% humidity, following approval by the Institutional Animal Care and Utilization Committee (IACUC-202305177). The rats were divided into six groups for various treatments (n = 10): sham, MCAO, MCAO + SC-L (3.75 g/kg/day), MCAO + SC-M (7.5 g/kg/day), MCAO + SC-H (15 g/kg/day), and MCAO + ginaton (0.1 g/kg/day). The drug administration of SC and ginaton groups was conducted through oral gavage from 3 days pro-surgery to 14 days post-surgery. Sham and MCAO groups received saline.
2.4 MCAO/R model establishment
The middle cerebral artery occlusion/reperfusion (MCAO/R) model was induced as per established methods (Bai et al., 2023). The external carotid artery was exposed through a midline incision after positioning the rats on a warming platform following their anesthesia with pentobarbital sodium (30 mg/kg). A nylon filament was used to occlude the middle cerebral artery for 2 h, followed by reperfusion. The sham procedure was identical, excluding filament insertion.
2.5 Neurological and infarct volume assessment
Fourteen days after MCAO/R, neurological function was assessed using the Zea longa 5-grade scoring system (Longa et al., 1989; Liu et al., 2023a). Brain slices were stained with 2,3,5-triphenyl tetrazolium chloride (TTC) for infarct volume determination, and analyzed with Image-Pro Plus software.
2.6 Rotarod test and corner test
The rotarod test (LE8200 Panlab, Harvard Apparatus, United States) was conducted to evaluate the coordination and balance in the forelimbs and hindlimbs. The rats were subjected to the accelerating rotor mode, which varied from 5 to 40 rpm over a duration of 300 s. An investigator, blinded to the experimental groups, recorded the duration each animal stayed on the rod. For statistical analysis, the final score was represented as the average time a rat was able to maintain its position on the rod across two trials (Xia et al., 2023). The corner test serves as a commonly utilized functional evaluation for unilateral sensorimotor cortical impairment. In a home cage, two cardboard plates were positioned at a 30° angle. Each rat was positioned between these plates and permitted to navigate freely to the corner. Following the ischemic and reperfusion injury to the sensorimotor cortex, the rats exhibited a tendency to turn biasedly towards the side corresponding with their brain damage (Han et al., 2015). The frequency of left turns made by the rats during 10 trials was documented.
2.7 H&E, TUNEL, and Nissl staining
Histological analysis was conducted to assess histological damage, neuronal loss, and apoptosis in the ischemic hemisphere. Brain tissue underwent fixation, paraffin embedding, and sectioning for histological analyses. To evaluate neuronal apoptosis levels, the stained tissue sections were observed under a fluorescence microscope at a magnification of 200x to visualize TUNEL-positive cells. Additionally, observations of H&E staining and Nissl staining were made at a magnification of 400x. This allowed for the assessment of various aspects of the brain tissue.
2.8 Network pharmacology analyses
SC’s active components, identified via UPLC-Q-TOF-MS/MS, were assessed for Oral Bioavailability (OB) and Drug Similarity (DL) using the Traditional Chinese Medicine Systems Pharmacology (TCMSP) database. Only compounds with OB ≥ 30% and DL ≥ 0.18 were identified and preselected as active ingredients. Target predictions were made using Swiss Target Prediction, while ischemic stroke-related targets were collated from multiple databases, including Online Mendelian Inheritance in Man (OMIM), Genecards, DrugBank, and Therapeutic Target Database (TTD). Overlapping targets underwent protein-protein interaction (PPI) analysis using STRING, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted through DAVID, and KEGG results were visualized using the bioinformatics platform (http://www.bioinformatics.com.cn/). The final composite-target-pathway network was created using Cytoscape software (version 3.7.1).
2.9 Molecular docking analyses
Molecular structures of SC components were sourced from PubChem and prepared Mol2 files using Chem3D. Protein targets’ crystal structures with high resolution were retrieved from the RCSB Protein Data Bank (PDB), and edited using PyMOL (version 1.7.2). The optimal conformations of ingredients at the active sites of protein targets were predicted using Auto Dock Tools (version 1.5.7). The interaction and binding modes with the lowest binding energies were analyzed using Discovery Studio 2020 Client and visualized.
2.10 Molecular dynamics simulation
Based on the promising outcomes achieved from molecule docking, additional insights into the dynamic behavior and binding mechanisms were revealed through molecular dynamics simulation (MD). The MD simulations were performed using Gromacs 2022 software, where small molecules were modeled with the GAFF force field, proteins with the AMBER14SB force field, and a TIP3P water model was employed to create a simulation environment for the complexes. These simulations followed periodic boundary conditions under constant pressure and temperature. To maintain stability during MD simulations, all hydrogen bonds were constrained using the LINCS algorithm with a time step of 2 fs Electrostatic interactions were computed utilizing the Particle-Mesh Ewald (PME) method within a cutoff distance of 1.2 nm. Non-bonded interactions had a cutoff set at 10 Å and updated every 10 steps. Temperature was kept at 298 K using V-rescale method for temperature coupling while pressure was regulated at 1 bar via Berendsen technique. Equilibrium simulations including both NVT and NPT ensembles were conducted for 100 ps at 298 K, followed by an extensive MD simulation of the complex for 100 ns with conformational snapshots recorded every 10 ps Post-simulation analysis involved trajectory examination using Visual Molecular Dynamics (VMD) and PyMOL software tools, along with calculation of root mean square deviation (RMSD) and root mean square fluctuation (RMSF) between protein and small molecule ligands employing g_mmpbsa program.
2.11 Transcriptomic analyses
Transcriptome analysis was performed on ischemic penumbra brain tissues from the sham, MCAO, and SC-H groups (n = 3). The extraction of Total RNA was carried out utilizing the TRIzol reagent, followed by assessing its purity and integrity. The reverse transcription of mRNA resulted in the creation of cDNA, and subsequently, the construction of cDNA libraries commenced. Afterward, sequencing using the 2 × 150 bp paired-end method was performed on an Illumina Novaseq™ 6000 sequence platform (LC-Bio Technology Co., Ltd. Hangzhou, China). The transcriptome was aligned to the reference genome using HISAT2 (version 2.2.1), and differential expression analysis was performed using the DEGseqR package with a significance threshold (q-value ≤0.05, |log2FC|>1). GO enrichment and KEGG pathway analysis were performed to identify significantly enriched GO terms and pathways related to the effects of SC on ischemic stroke based on differentially expressed genes (DEGs).
2.12 Metabolomics analyses
Serum samples from the sham, MCAO, and SC-H groups were subjected to extraction using 50% methanol buffer and subsequent centrifugation. Ultra-high-performance liquid chromatography (ThermoFisher, Germany) coupled with a high-resolution tandem mass spectrometer TripleTOF6600plus (SCIEX, United States) was used for analysis (LC-Bio Technology Co., Ltd. Hangzhou, China). Mass spectrometry data were processed using XCMS software, and differential metabolites were identified based on variable importance in the projection (VIP) scores >1 and q-value ≤0.05. Differential metabolites were identified and analyzed for pathway implications using KEGG and HMDB databases.
2.13 Multi-omics integrated analysis
To investigate the relationship among various omics data, an integrated analysis was performed. Initially, homologous gene conversion was executed, converting all IDs to human gene equivalents, followed by standardization through the UniProt database. Subsequently, the predicted targets of SC were compared with IS-related genes and differentially expressed genes. The overlapping targets were deemed as high-confidence candidates and were analyzed utilizing the DAVID database. The KEGG outcomes were illustrated using a bioinformatics platform (http://www.bioinformatics.com.cn/). Ultimately, the Joint Pathway Analysis module in MetaboAnalyst 6.0 synthesized results from transcriptomic and metabolomic analyses for the identified treatment-related DEGs and differential metabolites. Scatterplots were generated to illustrate the identified KEGG pathways, with those exhibiting a p-value of less than 0.05 considered significantly enriched.
2.14 Western blotting
Standard western blotting was performed following the established methodology (Zhang et al., 2021). The brain tissue samples were homogenized and lysed before undergoing centrifugation. Protein concentrations were ascertained using a bicinchoninic acid assay. Proteins were then separated via SDS-PAGE and transferred onto PVDF membranes. To reduce non-specific interactions, membranes were blocked with 5% skim milk. Overnight incubation at 4°C with primary antibodies against PI3K, p-PI3K, AKT, p-AKT, HIF-1α, β-Actin, GAPDH (1:1,000, CST), VEGFA (1:1,000, Abcam), and MMP9 (1:1,000, Service Bio) was performed. This was followed by incubation with horseradish peroxidase-conjugated secondary antibody (goat anti-rabbit IgG, 1:3,000, Service Bio) and visualization. Protein band detection and quantification were conducted using ImageJ software (version 1.53).
2.15 Immunofluorescence assay
In the immunofluorescence assay, selected brain tissue sections underwent staining with primary antibodies against rabbit anti-VEGFA and anti- HIF-1α (1:500, CST). After appropriate blocking and washing, sections were incubated with fluorescently labeled secondary antibodies. Nuclei were counterstained with DAPI, and the slides were examined under a Nikon Eclipse C1 microscope, providing insights into the cellular and molecular changes induced by SC.
2.16 Statistical analysis
Data presentation utilized mean ± SD format, with group differences assessed through one-way ANOVA and Tukey’s post hoc test. A p-value <0.05 was deemed statistically significant. GraphPad Prism software (version 10.0.2) facilitated all statistical analyses.
3 Results
3.1 Chemical compounds analysis of SC herb pair decoction
Employing UPLC-Q-TOF-MS/MS, 46 principal compounds were identified in SC herb pair decoction, including phenolic acids, tanshinones, paeoniflorins, paeonol, and tannins, based on retention times and MS spectra from literature or reference standards (Figure 1; Table 1). These compounds were recognized as potentially active and used for subsequent target and pathway prediction via network pharmacology.
Figure 1. Comprehensive LC-MS/MS profiling of SC herb pair constituents. (A) Total ion chromatogram displaying the profile of SC herb pair in negative ion mode. (B) Total ion chromatogram displaying the profile of SC herb pair in positive ion mode.
Table 1. Forty-six key chemical constituents were identified from the SC herb pair using UPLC-Q-TOF-MS/MS in the positive and negative ion modes.
3.2 SC treatment alleviates MCAO-associated injury in rats
SC treatment markedly improved neurological outcomes in MCAO-induced rats, significantly reducing deficit scores (p < 0.01) (Figure 2A). Infarct volumes, assessed by TTC staining, were notably less in SC-H and ginaton groups (p < 0.001) (Figures 2B, C). Next, we performed the rotarod test and corner test to assess motor function and found that the SC-treated rats significantly increased latency to fall off the rod (Figure 2D), and increased the number of left turns in the corner test (Figure 2E). H&E staining demonstrated that the MCAO group exhibited nuclear deformation and neuronal loss, leading to serious cell structure damage in the penumbra, while reduced cellular damage in SC-treated groups compared to the MCAO group (Figure 2F). TUNEL staining (Figure 2G) showed a dose-dependent reduction in apoptosis in the SC-treated groups. Nissl staining (Figure 2H) revealed a significant reduction in the number of Nissl bodies and disrupted neuronal morphology in the MCAO group, while SC-treated groups exhibited an increase in Nissl bodies and improved neuronal morphology. These findings demonstrate that SC reduces brain damage, enhances neuroprotection, and restores brain function after ischemic stroke.
Figure 2. Therapeutic efficacy of SC herb pair in mitigating ischemic stroke induced by MCAO (n = 10). (A) Graphical representation of neurological deficit scores among different rat groups. (B) TTC staining images depicting cerebral infarction in various treatment groups. (C) Quantitative analysis of cerebral infarct area across different groups. (D) Latency time for rats to fall from the rotating rod. (E) The turning bias in the corner test. (F–H) Microscopic images of brain tissues stained with H&E, TUNEL, and Nissl (Scale bar = 50 μm). *Statistical significance indicated: *p < 0.05, **p < 0.01, ***p < 0.001 vs MCAO group; #p < 0.05, ##p < 0.01, ###p < 0.001 vs sham group.
3.3 Predictive network pharmacology
A total of 502 predicted targets for effective active compounds of SC were identified by using the Swiss Target Prediction, and 300 targets with high relevance scores were filtered from 4824 IS-related targets obtained in GeneCards (4,789), OMIM (20), TTD (17), and Drugbank (57) database using the keyword “Ischemic stroke” on 30 April 2023. Ultimately, 67 overlapping targets between SC and IS were identified as pivotal targets (Figure 3A). The PPI network highlighted important proteins including AKT1, TNF, PTGS2, MMP9, PIK3CA, VEGFA, and so on (Figure 3B). GO and KEGG pathway analyses connected these targets to positive regulation of angiogenesis, response to hypoxia, apoptosis regulation, and pathways like PI3K-AKT, TNF, VEGF, and HIF-1α signaling in IS (Figures 3C, D). A comprehensive network illustrating the interactions among SC components, targets, and pathways was constructed (Figure 3E).
Figure 3. Network pharmacology analysis elucidating SC’s mechanism of action in ischemic stroke. (A) Venn diagram depicting the overlap of targets between SC and ischemic stroke. (B) PPI network of shared targets. (C) GO enrichment analysis results of shared targets. (D) KEGG pathway analysis of core genes in the network. (E) Integrated network connecting ischemic stroke targets with SC ingredients.
3.4 Molecular docking and molecular dynamics simulation of key ingredients
Molecular docking analysis was conducted between six target proteins identified in the PPI analysis and six potentially active compounds. The results revealed good binding potential between key components in SC and target genes, including PIK3CA, TNF, VEGFA, AKT1, MMP9, and PTGS2 (Figures 4A–F). A total of 31 pairs were highlighted with significant binding activity with the affinity < −7 kcal/mol (Table 2). Next, salvianolic acid B, benzoylpaeoniflorin, and salvianolic Acid A were selected as core components for molecular dynamics simulations with PIK3CA, TNF, VEGFA, AKT1, MMP9, and PTGS2. The binding of salvianolic acid B to the ligands reached a steady state rapidly without excessive fluctuations. After a brief fluctuation, the binding of benzoylpaeoniflorin and salvianolic acid A to the ligand also formed a steady state. The results are shown in Figures 4G, H and Supplementary Figure S1.
Figure 4. Molecular docking interaction diagrams and molecular dynamics simulation between SC components and target proteins. Schematic diagrams (3D and 2D) of the docking models. (A) Salvianolic acid B with PIK3CA, (B) benzoylpaeoniflorin with TNF (C) salvianolic acid A with VEGFA, (D) benzoylpaeoniflorin with AKT1 (E) salvianolic acid B with MMP9, and (F) salvianolic acid B with PTGS2. (G) RMSD plot during molecular dynamic simulations of salvianolic acid B with PIK3CA. (H) RMSF plot during molecular dynamic simulations of salvianolic acid B with PIK3CA.
3.5 Transcriptomic insights
Transcriptomic analysis was performed on brain tissues from the sham, MCAO, and SC-H groups. Analysis of gene expression differences revealed that in the MCAO vs sham group, there were 360 genes found to be upregulated and 282 genes found to be downregulated. Similarly, in the SC-H vs MCAO group, 384 genes were upregulated and 433 genes were downregulated (Figures 5A, B). Meanwhile, we screened out 191 differentially expressed genes (DEGs) in the three groups by a Venn diagram (Figure 5C). The heatmap was used for mapping and clustering analysis of the 40 common DEGs, including VEGFA, SPP1, FGFR1, CNKN1A, and BDNF (Figure 5D). The identified DEGs predominantly correlated with biological processes, including cell migration, proliferation, ERK cascade, and cytokine response (Figure 5E). KEGG pathway analysis indicated involvement in pathways such as the PI3K-AKT, NF-kappa B, TNF, and MAPK (Figure 5F).
Figure 5. Transcriptomic revealed DEGs for the anti-ischemic stroke effect of SC (n = 3). (A, B) Volcano plots representing differentially expressed genes (DEGs) between MCAO vs sham and SC vs MCAO groups. (C) Overlapping DEGs illustrated via a Venn diagram. (D) Heatmap showcasing hierarchical clustering of DEGs between samples. (E, F) GO and KEGG pathway analyses highlight significant biological processes and pathways influenced by SC in MCAO models.
3.6 Metabolomics analysis of the changes induced by SC
Significant variations in serum metabolites among the three groups were observed through UPLC-Q-TOF-MS analysis and multivariate statistical analysis. The PCA and PLS-DA exhibited a distinct trend of separation among the three groups, indicating that the metabolites of SC differed significantly from those of the sham and MCAO groups, as presented in Figures 6A, B. A total of 444 metabolites were found to be differentially expressed in the MCAO vs sham group. In addition, 1,039 metabolites exhibited differential expression in the SC vs MCAO group, as demonstrated in Figure 6D. Sixty differential metabolites were selected and depicted in a heatmap (Figure 6C). Identification and annotation of these differential metabolites and trends in the serum of rats were shown in Supplementary Table S1. Metabolomic pathway analysis suggested associations with amino acids, glycosylphosphatidylinositol (GPI)-anchor, unsaturated fatty acids biosynthesis, alpha-linolenic acid, butanoate, glycerophospholipid, and inositol phosphate metabolism (Figures 6E, F).
Figure 6. Metabolomic landscape altered by SC in ischemic stroke (n = 3). (A) PCA of metabolite profiles in sham, MCAO, and SC groups. (B) PLS-DA scores contrasting MCAO vs sham and SC vs MCAO. (C) Volcano plots illustrating metabolic changes between MCAO vs sham and SC vs MCAO. (D) Heatmap delineating key metabolites across groups. (E, F) Pathway analyses of significant metabolites distinguishing MCAO from sham and SC-treated groups.
3.7 Integrated analysis of multiple omics
To obtain high-confidence candidates, we accessed 4824 IS-related genes from various databases using the keyword “Ischemic stroke” in previous studies. After comparing the predicted targets of SC, IS-related genes, and DEGs obtained from transcriptome data, with 40 intersecting targets were considered high-confidence candidates (Figure 7A), and shown in Supplementary Table S2. In the KEGG analysis, SC affected the TNF, HIF-1, VEGF, and PI3K-AKT signaling pathways, etc (Figure 7B). The correlation analysis was performed between the 817 DEGs and 1,039 metabolites using MetaboAnalyst 6.0 for Joint Pathway Analysis, and a total of 49 pathways associated with these genes or metabolites were identified. Among these metabolic pathways, including glycerophospholipid metabolism, glutathione metabolism, glycerolipid metabolism, phenylalanine metabolism, and inositol phosphate metabolism were notable (Figures 7C, D).
Figure 7. Integrated analysis of multi-omics on SC treating IS. (A) Venn diagram depicting the overlap of targets between SC, IS, and transcriptome. (B) Correlation analysis of network pharmacology and transcriptomics in the KEGG pathway. (C) The bubble diagram from the combined analysis of transcriptomics and untargeted metabolomics based on the MetaboAnalyst database. (D) Correlation analysis of untargeted metabolomics and transcriptomics in the KEGG pathway.
3.8 VEGFA/PI3K/AKT pathway modulation by SC
Western blot analysis indicated increases in PI3K and AKT phosphorylation levels in the MCAO group (p < 0.05), which were further enhanced in the SC-treated groups, with the SC-H group exhibiting the most significant increase (p < 0.001). Additionally, the protein expression of VEGFA, HIF-1α, and MMP9 followed a similar trend, and enhanced the VEGFA/PI3K/AKT signaling pathway by SC treatment, especially in the SC-H group (Figures 8A–F). Immunofluorescence studies, detailed in Figures 8G, H, corroborated these findings, showing an increase in VEGFA and HIF-1α expressions upon SC administration.
Figure 8. Impact of SC on key signaling proteins in the ischemic stroke pathway. (A) Western blot images showing expression levels of pathway proteins in ischemic brain tissue. (B, C) Quantitative analysis of protein expressions, specifically p-PI3K/PI3K and p-AKT/AKT ratios. (D–F) Quantitative evaluation of VEGFA, MMP9, and HIF-1α protein levels in brain tissues. Data presented as mean ± SD (n = 3). Statistical significance: *p < 0.05, **p < 0.01, ***p < 0.001 compared to MCAO group; #p < 0.05 compared to sham group. (G, H) Detailed immunofluorescence staining results for VEGFA and HIF-1α, respectively (Scale bar = 50 μm).
4 Discussion
Ischemic stroke is a condition characterized by high recurrence, disability, and mortality rates. Its underlying mechanisms involve multiple factors including oxidative stress, inflammation, apoptosis, and angiogenesis. In traditional Chinese medicine (TCM), herb pairs are frequently used to enhance the therapeutic effects or mitigate potential side effects in prescriptions (Wang et al., 2012; Pang et al., 2024). These herb pairs act as a bridge between single herbs and polyherbal medicines, making them a focus of research in the field of rehabilitation therapy for IS. For example, the efficacy of Danhong Injection (DHI) has been demonstrated in enhancing brain function scores and mitigating cell apoptosis in the brain tissues of individuals with IS. DHI is derived from a combination of S. miltiorrhiza Bunge and Carthamus tinctorius L. (CT) (Feng et al., 2020; Zhang et al., 2023a). Guhong injection (GHI), which is composed of CT extract and aceglutamide, has demonstrated its ability to have antiapoptotic, anti-inflammatory, and neuroprotective effects in myocardial ischemia/reperfusion injury (Zhang J. et al., 2022; Chen et al., 2023). While, because of the complexity of herbal compounds, understanding their synergistic mechanisms at the molecular level remains challenging.
Our study first confirmed that SC could have anti-cerebral ischemic injury effects, inhibit the apoptosis of neuronal cells, promote angiogenesis, and have neuroprotective effects. Based on neurobehavioral and pathological observations, treatment with mid-and high-dose SC (7.5 g/kg/d and 15 g/kg/d) for 14 days post-surgery significantly reduced neurological deficit scores, infarction area, neuronal loss, and improved neuronal morphology. We endeavored to elucidate the specific active ingredients of SC against cerebral ischemic injury. Through UPLC-Q-TOF-MS/MS analysis, we confirmed that oxypaeoniflora, benzoylpaeoniflorin, salvianolic acid B, salvianolic acid A, and tanshinone IIA had neuroprotective activity. The research found that oxypaeoniflora improves MI/R injury by activating the Sirt1/Foxo1 signaling-mediated anti-apoptotic pathway (Wang and Hu, 2020). Tanshinone IIA has been shown to have neuroprotective effects including anti-apoptotic, anti-inflammatory, and antioxidant properties (Sherawat and Mehan, 2023). In ischemic rats, salvianolic acid B has the potential to safeguard against neuronal apoptosis and mitigate neurological damage by promoting angiogenesis, a process that might be linked to the increase in VEGF and VEGFR2 through the targeting of STC1 (Bi et al., 2022). The systems pharmacology revealed that SC could regulate 67 common targets related to IS, including core targets like AKT1, TNF, PTGS2, MMP9, PIK3CA, and VEGFA, all closely linked to inflammation, cell apoptosis, and angiogenesis. The molecular docking and molecular dynamics simulation results indicated that the active ingredients mentioned above exhibited a strong affinity with core targets. This suggests that SC may have a therapeutic effect on IS by modulating these target-related pathways.
In recent years, integrated multi-omics analysis has emerged as the estimable approach to uncover the molecular mechanisms of drug action (Cheng et al., 2024; Li et al., 2024; Sun et al., 2023). We conducted a combined systems pharmacology, transcriptome, and metabolome analysis to elucidate the mechanism of action of SC against IS. Transcriptomics results revealed that SC treatment upregulated the expression of key genes, including VEGFA, secreted phosphoprotein 1 (SPP1), type I fibroblast growth factor receptor (FGFR1), cyclin-dependent kinase inhibitor 1A (CDKN1A), and brain-derived neurotrophic factor (BDNF), in pathways related to inflammation, angiogenesis, and apoptosis such as PI3K/AKT, IL-17, TNF, and NF-kappa B signaling. These findings suggest that SC’s protective effect against IS may be associated with the regulation of these pathways.
As previously reported, VEGFA is a key driver of angiogenesis and is considered indispensable in regulating blood vessel formation (Guo et al., 2021; Zhang et al., 2023b). VEGFA activates VEGF receptor-2 (VEGFR2), which, in turn, triggers downstream signaling pathways like PI3K/AKT and MEK/ERK, ultimately regulating cell development, growth, and differentiation (Kilic et al., 2006). The PI3K family of enzymes, which generates PIP3 upon activation, plays a vital role in controlling cellular proliferation and apoptosis (Koch et al., 2021; Goyal et al., 2023). The PI3K/AKT pathway has been shown to regulate angiogenesis, inflammation, and vascular endothelial homeostasis after IS (Koutsaliaris et al., 2022; Cheng et al., 2023). Additionally, PTGS2 is involved in the regulation of angiogenesis and is associated with the VEGF signaling pathway in various diseases (Jimenez-Luna et al., 2021). MMP9 is critical in degrading the basement membrane, marking the initiation of angiogenesis (Quintero-Fabián et al., 2019). HIF-1α is among the first genes to be upregulated by hypoxia in ischemic conditions and directly stimulates VEGF gene transcription by binding to the VEGF promoter (Fu et al., 2018). The integrated analysis of network pharmacology and transcriptomics revealed the VEGF pathway and PI3K/AKT signaling pathway as potential pivotal targets for SC treatment of IS. Subsequently, we investigated the potential impact of SC on the PI3K/AKT signaling pathway. Western blot and immunofluorescence staining partially validated that SC could activate the PI3K/AKT signaling pathway, induce AKT phosphorylation, inhibit degradation of HIF-1α, facilitate its accumulation, and subsequently translocate it to the nucleus for initiating transcriptional activation of VEGF. Consequently, enhanced expression of VEGF occurs which upon binding to its receptor VEGFR2 activates the downstream PI3K/AKT signaling cascade. These findings underscore the remarkable potential of SC in promoting recovery from IS and safeguarding neural tissue. There are still limitations in our study and more experiments are needed to explore the intrinsic mechanisms.
In addition, the transcriptomics analysis revealed that a substantial proportion of the genes displaying differential expression were significantly enriched in metabolic pathways. This observation highlights the potential impact of SC on regulating metabolic processes. Clinical studies have associated higher blood and cerebrospinal fluid concentrations of glutamate with neurological deficits in acute ischemic stroke (Castillo et al., 1997). Furthermore, branched-chain amino acids like valine, leucine, and isoleucine have been linked to various stages of stroke and may serve as diagnostic markers for IS (Sidorov et al., 2019). Altered levels of lyso-PI and lyso-PS in serum are indicative of impaired phospholipid metabolism and could potentially serve as biomarkers for IS (Sun et al., 2017). Acetylcarnitine, linked to fatty acid metabolism, increased during IS and is considered a potential target for IS therapy (Huang et al., 2023). Integrated analysis of metabolomics and transcriptomics revealed that SC treatment corrected disturbances in amino acids, glycerophospholipids, sphingomyelin, and fatty acid metabolism in MCAO rats. Glutamate, phenylalanine, leucine, lyso-phosphatidylinositol (lyso-PI), lyso-phosphatidylserine (lyso-PS), and acylcarnitine are key metabolites in the neuroprotective effect of SC.
One limitation of this study is the exclusive use of male rat models, which overlooks the physiological and hormonal disparities between males and females. It is crucial to acknowledge that the therapeutic effects of the SC herbal combination may vary based on sex. To obtain a more comprehensive understanding of its efficacy and safety, future studies should incorporate female rats as well. By including both sexes in the research, more generalized conclusions can be drawn, better reflecting the potential clinical applications of the SC herbal combination across diverse patient populations.
5 Conclusion
In summary, our research first time demonstrates the neuroprotective efficacy of SC in a rat MCAO model, mediated through the activation of the VEGFA/PI3K/AKT pathway, as evidenced by multi-omics analyses. These results suggest SC’s potential as a therapeutic agent for IS. It is worth mentioning that the scope of this investigation is restricted to research conducted at the animal level. Subsequently, additional exploration is required to assess the clinical potential of SC.
Data availability statement
The transcriptome library construction and sequencing data presented in the study was deposited in the NCBI SRA BioProject repository, and the accession number was PRJNA1110725.
Ethics statement
The animal study was approved by Institutional Animal Care and Utilization Committee of Air Force Military Medical University. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
CZ: Formal Analysis, Investigation, Writing–original draft. XB: Investigation, Validation, Visualization, Writing–original draft. YD: Conceptualization, Formal Analysis, Funding acquisition, Methodology, Writing–review and editing. AW: Conceptualization, Funding acquisition, Writing–review and editing. QF: Conceptualization, Funding acquisition, Methodology, Project administration, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was financially supported by the National Natural Science Foundation of China (Nos. 82074321, 81773689, and 82274313), Natural Science Foundation of Shaanxi Province (No. 2021SF-430), Key Construction Discipline Research Capacity Enhancement Project of Guangdong Province (2022ZDJS119), and Shaanxi Provincial Key R&D Program Social Development Project (2023YBSF-564).
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.2024.1431692/full#supplementary-material
Abbreviations
IS, Ischemic stroke; SM, Salvia miltiorrhiza; CM, Cortex moutan; SC, Salvia miltiorrhiza-Cortex moutan; MCAO, middle cerebral artery occlusion; STRING, String online platform; C-T-P, compound-target-pathway; SPP1, secreted phosphoprotein 1; FGFR1, type I fibroblast growth factor receptor.
References
Bai, M., Cui, N., Liao, Y., Guo, C., Li, L., Yin, Y., et al. (2023). Astrocytes and microglia-targeted Danshensu liposomes enhance the therapeutic effects on cerebral ischemia-reperfusion injury. J. Control Release 364, 473–489. doi:10.1016/j.jconrel.2023.11.002
Bi, S.-J., Dong, X.-Y., Wang, Z.-Y., Fu, S.-J., Li, C.-L., Wang, Z.-Y., et al. (2022). Salvianolic acid B alleviates neurological injury by upregulating stanniocalcin 1 expression. Ann. Transl. Med. 10, 739. doi:10.21037/atm-21-4779
Castillo, J., Dávalos, A., and Noya, M. (1997). Progression of ischaemic stroke and excitotoxic aminoacids. Lancet 349, 79–83. doi:10.1016/S0140-6736(96)04453-4
Chen, H., Zhou, H., Yang, J., Wan, H., and He, Y. (2023). Guhong injection mitigates myocardial ischemia/reperfusion injury by activating GST P to inhibit ASK1-JNK/p38 pathway. Phytomedicine 109, 154603. doi:10.1016/j.phymed.2022.154603
Cheng, M., Li, T., Hu, E., Yan, Q., Li, H., Wang, Y., et al. (2024). A novel strategy of integrating network pharmacology and transcriptome reveals antiapoptotic mechanisms of Buyang Huanwu Decoction in treating intracerebral hemorrhage. J. Ethnopharmacol. 319, 117123. doi:10.1016/j.jep.2023.117123
Cheng, X., Yang, X., Cui, H., Zhang, B., Chen, K., Yang, X., et al. (2023). Chuanxiong improves angiogenesis via the PI3K/AKT/Ras/MAPK pathway based on network pharmacology and DESI-MSI metabolomics. Front. Pharmacol. 14, 1135264. doi:10.3389/fphar.2023.1135264
Fei, Y.-X., Wang, S.-Q., Yang, L.-J., Qiu, Y.-Y., Li, Y.-Z., Liu, W.-Y., et al. (2017). Salvia miltiorrhiza Bunge (Danshen) extract attenuates permanent cerebral ischemia through inhibiting platelet activation in rats. J. Ethnopharmacol. 207, 57–66. doi:10.1016/j.jep.2017.06.023
Feng, C., Wan, H., Zhang, Y., Yu, L., Shao, C., He, Y., et al. (2020). Neuroprotective effect of Danhong injection on cerebral ischemia-reperfusion injury in rats by activation of the PI3K-Akt pathway. Front. Pharmacol. 11, 298. doi:10.3389/fphar.2020.00298
Fu, C., Tyagi, R., Chin, A. C., Rojas, T., Li, R.-J., Guha, P., et al. (2018). Inositol polyphosphate multikinase inhibits angiogenesis via inositol pentakisphosphate-induced HIF-1α degradation. Circ. Res. 122, 457–472. doi:10.1161/CIRCRESAHA.117.311983
Goyal, A., Agrawal, A., Verma, A., and Dubey, N. (2023). The PI3K-AKT pathway: a plausible therapeutic target in Parkinson’s disease. Exp. Mol. Pathol. 129, 104846. doi:10.1016/j.yexmp.2022.104846
Guo, X., Yi, H., Li, T. C., Wang, Y., Wang, H., and Chen, X. (2021). Role of vascular endothelial growth factor (VEGF) in human embryo implantation: clinical implications. Biomolecules 11, 253. doi:10.3390/biom11020253
Han, L., Cai, W., Mao, L., Liu, J., Li, P., Leak, R. K., et al. (2015). Rosiglitazone promotes white matter integrity and long-term functional recovery after focal cerebral ischemia. Stroke 46, 2628–2636. doi:10.1161/STROKEAHA.115.010091
Henderson, S. J., Weitz, J. I., and Kim, P. Y. (2018). Fibrinolysis: strategies to enhance the treatment of acute ischemic stroke. J. Thromb. Haemost. 16, 1932–1940. doi:10.1111/jth.14215
Huang, X.-X., Li, L., Jiang, R.-H., Yu, J.-B., Sun, Y.-Q., Shan, J., et al. (2023). Lipidomic analysis identifies long-chain acylcarnitine as a target for ischemic stroke. J. Adv. Res. S2090-1232 (23), 133–149. doi:10.1016/j.jare.2023.08.007
Jimenez-Luna, C., González-Flores, E., Ortiz, R., Martínez-González, L. J., Antúnez-Rodríguez, A., Expósito-Ruiz, M., et al. (2021). Circulating PTGS2, JAG1, GUCY2C and PGF mRNA in peripheral blood and serum as potential biomarkers for patients with metastatic colon cancer. J. Clin. Med. 10, 2248. doi:10.3390/jcm10112248
Kilic, E., Kilic, Ü., Wang, Y., Bassetti, C. L., Marti, H. H., Hermann, D. M., et al. (2006). The phosphatidylinositol-3 kinase/Akt pathway mediates VEGF’s neuroprotective activity and induces blood brain barrier permeability after focal cerebral ischemia. FASEB J. 20, 1185–1187. doi:10.1096/fj.05-4829fje
Koch, P. A., Dornan, G. L., Hessenberger, M., and Haucke, V. (2021). The molecular mechanisms mediating class II PI 3-kinase function in cell physiology. FEBS J. 288, 7025–7042. doi:10.1111/febs.15692
Koh, S. H., and Lo, E. H. (2015). The role of the PI3K pathway in the regeneration of the damaged brain by neural stem cells after cerebral infarction. J. Clin. Neurol. 11, 297–304. doi:10.3988/jcn.2015.11.4.297
Koutsaliaris, I. K., Moschonas, I. C., Pechlivani, L. M., Tsouka, A. N., and Tselepis, A. D. (2022). Inflammation, oxidative stress, vascular aging and atherosclerotic ischemic stroke. Curr. Med. Chem. 29, 5496–5509. doi:10.2174/0929867328666210921161711
Lee, T.-Y., Chang, H.-H., Wang, G.-J., Chiu, J.-H., Yang, Y.-Y., and Lin, H.-C. (2006). Water-soluble extract of Salvia miltiorrhiza ameliorates carbon tetrachloride-mediated hepatic apoptosis in rats. J. Pharm. Pharmacol. 58, 659–665. doi:10.1211/jpp.58.5.0011
Li, J., Ma, Z., Yang, Z., Yang, M., Li, C., Li, M., et al. (2024). Integrating transcriptomics and network pharmacology to reveal the mechanisms of total Rhizoma Coptidis alkaloids against nonalcoholic steatohepatitis. J. Ethnopharmacol. 322, 117600. doi:10.1016/j.jep.2023.117600
Li, R., Zheng, Y., Zhang, J., Zhou, Y., and Fan, X. (2023a). Gomisin N attenuated cerebral ischemia-reperfusion injury through inhibition of autophagy by activating the PI3K/AKT/mTOR pathway. Phytomedicine 110, 154644. doi:10.1016/j.phymed.2023.154644
Li, X., Liu, R., Liu, W., Liu, X., Fan, Z., Cui, J., et al. (2023b). Panax quinquefolium L. And Salvia miltiorrhiza Bunge. Enhances angiogenesis by regulating the miR-155-5p/HIF-1α/VEGF Axis in acute myocardial infarction. Drug Des. Devel Ther. 17, 3249–3267. doi:10.2147/DDDT.S426345
Li, X., Yang, Y., Song, Y., Nie, F., Fu, C., and Qin, Y. (2022). Effect of shuangdan mingmu capsule on diabetic retinopathy in rats via regulation of miRNAs. Diabetes Metab. Syndr. Obes. 15, 3181–3194. doi:10.2147/DMSO.S379611
Li, Y., Bao, J.-X., Xu, J., Murad, F., and Bian, K. (2010). Vascular dilation by paeonol--a mechanism study. Vasc. Pharmacol. 53, 169–176. doi:10.1016/j.vph.2010.07.001
Liberale, L., Puspitasari, Y. M., Ministrini, S., Akhmedov, A., Kraler, S., Bonetti, N. R., et al. (2023). JCAD promotes arterial thrombosis through PI3K/Akt modulation: a translational study. Eur. Heart J. 44, 1818–1833. doi:10.1093/eurheartj/ehac641
Liu, T., Bai, M., Liu, M., Li, T., Liao, Y., Zhao, C., et al. (2023a). Novel synergistic mechanism of 11-keto-β-boswellic acid and Z-Guggulsterone on ischemic stroke revealed by single-cell transcriptomics. Pharmacol. Res. 193, 106803. doi:10.1016/j.phrs.2023.106803
Liu, T., Wang, W., Li, X., Chen, Y., Mu, F., Wen, A., et al. (2023b). Advances of phytotherapy in ischemic stroke targeting PI3K/Akt signaling. Phytother. Res. 37, 5509–5528. doi:10.1002/ptr.7994
Longa, E. Z., Weinstein, P. R., Carlson, S., and Cummins, R. (1989). Reversible middle cerebral artery occlusion without craniectomy in rats. Stroke 20, 84–91. doi:10.1161/01.str.20.1.84
Mendelson, S. J., and Prabhakaran, S. (2021). Diagnosis and management of transient ischemic attack and acute ischemic stroke: a review. JAMA 325, 1088–1098. doi:10.1001/jama.2020.26867
Pang, H.-Q., Guo, J.-X., Yang, Y., Xu, L., Wang, J., Yang, F., et al. (2024). Elucidating the chemical interaction effects of herb pair Danshen-Chuanxiong and its anti-ischemic stroke activities evaluation. J. Ethnopharmacol. 318, 117058. doi:10.1016/j.jep.2023.117058
Quintero-Fabián, S., Arreola, R., Becerril-Villanueva, E., Torres-Romero, J. C., Arana-Argáez, V., Lara-Riegos, J., et al. (2019). Role of matrix metalloproteinases in angiogenesis and cancer. Front. Oncol. 9, 1370. doi:10.3389/fonc.2019.01370
Sherawat, K., and Mehan, S. (2023). Tanshinone-IIA mediated neuroprotection by modulating neuronal pathways. Naunyn Schmiedeb. Arch. Pharmacol. 396, 1647–1667. doi:10.1007/s00210-023-02476-8
Sidorov, E., Sanghera, D. K., and Vanamala, J. K. P. (2019). Biomarker for ischemic stroke using metabolome: a clinician perspective. J. Stroke 21, 31–41. doi:10.5853/jos.2018.03454
Sun, H., Zhao, J., Zhong, D., and Li, G. (2017). Potential serum biomarkers and metabonomic profiling of serum in ischemic stroke patients using UPLC/Q-TOF MS/MS. PLoS One 12, e0189009. doi:10.1371/journal.pone.0189009
Sun, J., Chen, X., Wang, Y., Song, Y., Pan, B., Fan, B., et al. (2023). Neuroprotective effects of Longxue Tongluo Capsule on ischemic stroke rats revealed by LC-MS/MS-based metabolomics approach. Chin. Herb. Med. 15, 430–438. doi:10.1016/j.chmed.2022.12.010
Sun, J., Huang, S. H., Tan, B. K.-H., Whiteman, M., Zhu, Y. C., Wu, Y. J., et al. (2005). Effects of purified herbal extract of Salvia miltiorrhiza on ischemic rat myocardium after acute myocardial infarction. Life Sci. 76, 2849–2860. doi:10.1016/j.lfs.2004.11.016
Virani, S. S., Alonso, A., Benjamin, E. J., Bittencourt, M. S., Callaway, C. W., Carson, A. P., et al. (2020). Heart disease and stroke statistics-2020 update: a report from the American heart association. Circulation 141, e139–e596. doi:10.1161/CIR.0000000000000757
Wang, F., He, Q., Wang, J., Yuan, Q., Guo, H., Chai, L., et al. (2017). Neuroprotective effect of salvianolate lyophilized injection against cerebral ischemia in type 1 diabetic rats. BMC Complement. Altern. Med. 17, 258. doi:10.1186/s12906-017-1738-8
Wang, K., and Hu, W. (2020). Oxypaeoniflorin improves myocardial ischemia/reperfusion injury by activating the Sirt1/Foxo1 signaling pathway. Acta Biochim. Pol. 67, 239–245. doi:10.18388/abp.2020_5206
Wang, S., Hu, Y., Tan, W., Wu, X., Chen, R., Cao, J., et al. (2012). Compatibility art of traditional Chinese medicine: from the perspective of herb pairs. J. Ethnopharmacol. 143, 412–423. doi:10.1016/j.jep.2012.07.033
Xia, Q., Zhang, X., Zhan, G., Zheng, L., Mao, M., Zhao, Y., et al. (2023). A cell-penetrating peptide exerts therapeutic effects against ischemic stroke by mediating the lysosomal degradation of sirtuin 5. MedComm 4, e436. doi:10.1002/mco2.436
Yang, W., Ju, J., Jeon, M. J., Han, X., and Shin, I. (2010). Danshen (Salvia miltiorrhiza) extract inhibits proliferation of breast cancer cells via modulation of Akt activity and p27 level. Phytother. Res. 24, 198–204. doi:10.1002/ptr.2945
Zhang, H., Wang, X., Chen, W., Yang, Y., Wang, Y., Wan, H., et al. (2023a). Danhong injection alleviates cerebral ischemia-reperfusion injury by inhibiting autophagy through miRNA-132-3p/ATG12 signal axis. J. Ethnopharmacol. 300, 115724. doi:10.1016/j.jep.2022.115724
Zhang, H., Zhou, J., Li, J., Wang, Z., Chen, Z., Lv, Z., et al. (2023b). N6-Methyladenosine promotes translation of VEGFA to accelerate angiogenesis in lung cancer. Cancer Res. 83, 2208–2225. doi:10.1158/0008-5472.CAN-22-2449
Zhang, J., Zhao, J., Ma, Y., Wang, W., Huang, S., Guo, C., et al. (2021). Investigation of the multi-target mechanism of guanxin-shutong capsule in cerebrovascular diseases: a systems pharmacology and experimental assessment. Front. Pharmacol. 12, 650770. doi:10.3389/fphar.2021.650770
Zhang, J., Zhou, R., Cao, G., Zhang, Y., Xu, H., and Yang, H. (2022a). Guhong injection prevents ischemic stroke-induced neuro-inflammation and neuron loss through regulation of C5ar1. Front. Pharmacol. 13, 818245. doi:10.3389/fphar.2022.818245
Keywords: Salvia miltiorrhiza, cortex moutan, ischemic stroke, transcriptomics, metabolomics, VEGFA/PI3K/AKT pathway
Citation: Zhao C, Bai X, Ding Y, Wen A and Fu Q (2024) Combining systems pharmacology, metabolomics, and transcriptomics to reveal the mechanism of Salvia miltiorrhiza-Cortex moutan herb pair for the treatment of ischemic stroke. Front. Pharmacol. 15:1431692. doi: 10.3389/fphar.2024.1431692
Received: 12 May 2024; Accepted: 27 August 2024;
Published: 09 September 2024.
Edited by:
Adolfo Andrade-Cetto, National Autonomous University of Mexico, MexicoReviewed by:
Sergei V. Fedorovich, Belarusian State University, BelarusZiyi Wang, Johns Hopkins University, United States
Yanzhi Liu, Guangdong Medical University, China
Copyright © 2024 Zhao, Bai, Ding, Wen and Fu. 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: Qiang Fu, ZnVxaWFuZ0BtYWlsLnhqdHUuZWR1LmNu
†These authors have contributed equally to this work