- 1Department of Neurology, Xuanwu Hospital, Capital Medical University, Beijing, China
- 2Capital Medical University, Beijing, China
- 3Institute of Chinese Materia Medica, China Academy of Chinese Medical Sciences, Beijing, China
Tanhuo formula (THF), a traditional Chinese medicinal formula, has been demonstrated to be effective in the clinical treatment of acute ischemic stroke (AIS). However, its active ingredients, potential targets, and molecular mechanisms remain unknown. Based on the validation of active ingredient concentrations, our study attempted to elucidate the possible mechanisms of THF based on network pharmacological analysis and experimental validation. Components of THF were screened using network pharmacological analysis, and a compound–target network and protein–protein interaction (PPI) network were constructed. In total, 42 bioactive compounds and 159 THF targets related to AIS were identified. The PPI network identified AKT1, TNF, IL6, IL1B, and CASP3 as key targets. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis demonstrated that the inflammation and apoptotic pathways were enriched by multiple targets. The main components of THF were identified via high-performance liquid chromatography. Subsequently, a validation experiment was conducted, and the expressions of GFAP, C3, TNF-α, and IL-6 were detected via immunofluorescence staining, confirming the inflammatory response at 30 min and 3 days post injury. Immunohistochemical staining for caspase-3 and TUNEL was also performed to assess apoptosis at the same time points. These results indicate that THF can effectively decrease neural cell apoptosis through the caspase-3 pathway and restrain excessive abnormal activation of astrocytes and the release of TNF-α and IL-6, which might be accompanied by the recovery of motor function. Thus, THF may serve as a promising therapeutic strategy for AIS through multiple targets, components, and pathways.
Introduction
Acute ischemic stroke (AIS), characterized by high morbidity, disability, and mortality, is one of the leading causes of human death and disability worldwide; it, thereby, imposes an immense financial burden on society and severe psychological stress on individuals (Herpich and Rincon, 2020; Gao et al., 2020). According to the World Health Organization, approximately 15 million people experience an ischemic stroke each year, of which approximately 5 million die (Esquiva et al., 2018; Krishnamurthi et al., 2020; Morovatdar et al., 2021). Recombinant tissue–type plasminogen activator, a unique, FDA-approved drug for the treatment of AIS, can only benefit approximately 5% of patients owing to the narrow treatment window for thrombolysis and the risk of cerebral hemorrhage induced by the drug (Saberi et al., 2021). Therefore, it is essential to develop more effective and safe drugs with minimal side effects for AIS treatment.
Increasing evidence has demonstrated that intervening inflammatory reactions are important targets in the clinical treatment of AIS. Astrocytes are important regulators of the neuroinflammatory response and have also been found to be effective targets for therapeutic intervention in ischemic stroke (Patabendige et al., 2021). Owing to the widespread distribution and intrinsic features of the central nervous system (CNS), astrocytes are rapidly activated and proliferated after ischemic stroke (Stary and Giffard, 2015). Activated astrocytes can immediately differentiate into two polarization phenotypes, A1 and A2. The former has pro-inflammatory and neurotoxic effects, whereas the latter has anti-inflammatory and neural restoration effects. Overactivated A1 astrocytes secrete various types of inflammatory cytokines and a large number of cytotoxic substances, such as ROS, iNOS, and glutamate. These substances amplify the inflammatory cascade reaction, aggravate the redistribution or degradation of key tight junctions (such as claudin-5 and occludin), and deteriorate the leakage of the blood–brain barrier, ultimately causing severe tissue damage and neuronal death (Xu et al., 2020; Qiu et al., 2021; Sun et al., 2021). Apoptosis is the major type of cell death at the location of cerebral infarction (Radak et al., 2017). Neurons are sensitive to the deterioration of the microenvironment and undergo massive apoptosis after a series of ischemia, hypoxia, oxidative stress, and inflammatory responses following AIS (Xu et al., 2020). Thus, inhibition of apoptosis and activation of A1 astrocytes are an essential neuroprotective strategy for AIS.
For many years, Tanhuo formula (THF), a traditional Chinese medicinal formula, has been generally acknowledged as an effective prescription for AIS caused by phlegm-heat syndrome in Xuanwu Hospital. Although our previous study confirmed that THF reduced infarct volume and attenuated neuronal damage in cerebral ischemic rats, the mechanism of antiapoptosis and astrocyte regulation has not yet been elucidated (Cong et al., 2021). THF comprises five herbal medicines: dried rhizome of Coptis chinensis Franch. [Ranunculaceae; Coptidis rhizoma; Huang-Lian (HL)], dried roots and rhizomes of Rheum officinale Baill., [Polygonaceae, Rhei radix et rhizoma, Da-Huang (DH)], dried fruits of Forsythia suspensa (Thunb.) Vahl [Oleaceae; Forsythiae fructus; Lian-Qiao (LQ)], dried leaves of Lophatherum gracile Brongn. [Poaceae; Lophatheri herba; Dan-Zhu-Ye (DZY)], and Rhizoma Arisaematis Cum Bile [Araceae; Arisaema cum bile; Dan-Nan-Xing (DNX)]. The herbal components of THF, such as aloe emodin in DH, phillyrin in LQ, berberine, and palmatine in HL, have been shown to inhibit inflammatory responses and neuronal apoptosis in the brain (Zhu et al., 2019; Guo X. et al., 2021; Jiang et al., 2021; Tang et al., 2021; Xian et al., 2021; Zhao et al., 2021). We determined the relatively high concentrations of these key components in THF using high-performance liquid chromatography (HPLC). Based on the required concentrations to exert pharmacological effects, it is reasonable to infer that these THF components can exert pharmacological effects in the human body and promote recovery by inhibiting inflammatory response and reducing apoptosis in AIS. However, THF consists of several herbs and contains a variety of components that act via multiple pathways and may influence each other. Therefore, we comprehensively investigated the pharmacodynamics of THF on AIS by considering a variety of components and using network pharmacology methods. Based on bioinformatics and systems biology, network pharmacology is used to study the multilevel network relationship of “drug–target–disease” by screening potential pharmacodynamic components, action targets, and pathways in a database. Although the formula for THF has multiple ingredients and targets, network pharmacology can help address this challenge by precisely determining the research target of THF and revealing the underlying mechanisms for AIS effectively. However, current network pharmacology studies focus on the chemical composition of herbal medicines instead of the dose–effect relationship, and the content of compounds also affects the effect of herbal medicines. Therefore, the prediction results of network pharmacology may be biased and need to be experimentally verified (Zhang et al., 2019). Herein, combined with content determination of the known active components, a comparative experiment of two doses was designed in this study to enhance the reliability of network pharmacology and further analyze the mechanisms of THF in AIS. The framework of this study is shown in Figure 1.
Materials and Methods
Materials and Reagents
LC-20AT HPLC (Shimadzu Corporation, Japan), Welch Ultimate XB-C18 column (250 mm × 4.6 mm, 5 μm), methanol and acetonitrile (chromatographic purity, Thermo Fisher Scientific, Shanghai, China), ultrapure water (Hangzhou Wahaha Beverage Co., Ltd., Hangzhou, China), ethanol, and other reagents were used as analytical reagents. LQ (batch number 140610003, purchased from Beijing Qiancao Chinese Medicine Tablet Co.), DNX (batch number 409051, purchased from Sinopharm Group Beijing Huamiao Pharmaceutical Co. Ltd.). HL (batch number 1402001), DZY (batch number 1406001), and DH (batch number 1402001) were purchased from Anguo Changda Chinese herbal medicine Pieces Co. Ltd. (Baoding, China). The identity of all the herbs was confirmed by a botanical expert, Professor Hui Li (Institute of Chinese Materia Medica, China Academy of Chinese Medical Sciences, Beijing, China), and the specimens were archived in the institute. Berberine hydrochloride (lot:110713-201212, purity: 86.7%) and phillyrin (lot:110821-201213, purity: 95.3%), obtained from the National Institutes for Food and Drug Control (Beijing, China), were used as reference standards.
Preparation of THF
The composition of THF was HL (9 g), DH (5 g), LQ (10 g), DZY (9 g), and DNX (9 g). LQ (166.7 g) and DH (83.3 g), ratio: 2:1, were extracted three times by reflux extraction with 65% ethanol (eight times the volume) for 2 h at a time and filtered and concentrated into a filtrate with a relative density of 1.10–1.20 (50–60°C). HL (150 g), DNX (150 g, wrap-boiling), and DZY (150 g), ratio: 1:1:1, were extracted with water three times (2 h each time): the first time with 10 times the volume of water and the second and third times with eight times the volume of water. The extracts were filtered and concentrated into a filtrate with a relative density of 1.20–1.30 (60–70°C). The two filtrates were mixed, dried under reduced pressure, and finally crushed into a powder.
Sample and Reference Solution Preparation
Preparation of standard solutions: berberine hydrochloride and phillyrin were prepared in methanol to yield 0.15 and 0.20 mg/ml control solution, respectively, and then filtered through a 0.45-μm membrane filter before HPLC analysis. To prepare a sample solution for berberine, epiberberine, coptisine, and palmatine: THF dried powder (0.5 g) was added to 50 ml of methanol-hydrochloric acid (100: 1) and then sonicated for 30 min. Methanol was added to restore the solution weight to the initial. After filtration, 2 ml of the filtrate was placed in a 10-ml volumetric flask and methanol was added till it reached the mark. The sample solution was obtained by filtration through a 0.45-μm membrane. To prepare a sample solution for phillyrin, 1.0 g of THF dried powder was weighed and 15 ml 70% methanol was added, follow by sonication for 30 min, 70% methanol was used to replace the lost weight, and ultimately the sample solution of phillyrin was obtained by filtration through a 0.45-µm membrane.
HPLC Analysis
The concentrations of THF compounds were analyzed by HPLC. HPLC analysis was performed using a Welch Ultimate XB-C18 column (250 mm × 4.6 mm, 5 μm).
HPLC conditions for berberine, epiberberine, coptisine, and palmatine were as follows: octadecylsilane-bonded silica gel was used as the filler for the HPLC column. Acetonitrile (0.05 mol/L): potassium dihydrogen phosphate (50:50) ( 0.4 g sodium dodecyl sulfate per 100 ml is added and then the pH is adjusted to 4.0, with phosphoric acid) were used as the mobile phase. Chromatographic analysis was performed at 30°C with a flow rate of 1 ml/min, detection wavelength of 345 nm, and an injection volume of 10 µl.
To perform HPLC for phillyrin, the mobile phase consisted of water (A) and acetonitrile (B), and the gradient elution program was (time/B%): 0–30 min, 23–24%; 30–31 min, 24–95%; 31–39 min, 95%; 39–40 min, 95–23%; 10–15 min, 50–80%; and 40–50 min, 23%. Chromatographic analysis was performed at 30°C with a flow rate of 1 ml/min, detection wavelength of 277 nm, and an injection volume of 10 µl.
Screening Conditions for Candidate Active Ingredients
The Traditional Chinese Medicine Database and Analysis Platform (TCMSP) was used to screen the active ingredients of THF. The TCMSP is a recognized Chinese herbal medicine platform that provides detailed pharmacokinetic data of herbal ingredients, such as oral bioavailability (OB), drug-likeness (DL), aqueous solubility, and the blood–brain barrier (BBB) etc.(Veber et al., 2002; Ru et al., 2014). Due to the breakdown of the neurovascular unit and BBB after cerebral infarction, huge amounts of herbal components directly enter the CNS without obstacles. Thus, the screening conditions on the platform were set as OB ≥30% and DL ≥0.18 (Huang et al., 2014).
Candidate Active Compounds and Ingredient-Related Target Prediction
THF compounds were clarified from the TCMSP (https://old.tcmsp-e.com/tcmsp.php), ETCM (http://www.tcmip.cn/ETCM/index.php/Home/Index/), TCMID (http://www.megabionet.org/tcmid/), and BATMAN (http://bionet.ncpsb.org/batman-tcm/index.php) databases and literature search. Subsequently, the TCMSP platform was used to identify target genes. In addition, the targets were predicted using the Swiss Target Prediction database (http://www.swisstargetprediction.ch/) for the quality marker (Q-maker) of each herb in THF. Subsequently, different ID types of the targets were transformed into UniProt IDs using the UniProt database (https://www.uniprot.org/).
Screening of Disease-Associated Targets
AIS-related targets were predicted using GeneCards (https://www.genecards.org) and OMIM (https://omim.org). “Cerebral infarction,” “stroke,” and “cerebral ischemia” were used as the keywords for target searching. The relevance score was set as greater than or equal to 5, duplicate targets were excluded, and cerebral infarction–related targets were screened out. Intersection of THF-AIS–related target genes was performed to verify common target genes.
Protein–Protein Interaction (PPI) and Component–Target Network
Common targets of THF and AIS were uploaded into the STRING database (https://string-db.org/) for PPI network mapping by setting species as “Homo sapiens” and moderate confidence. Subsequently, the PPI network and component–target network were created and visualized using Cytoscape 3.8.0 software, and the degree centrality, betweenness centrality, and closeness centrality were calculated. The size and color of a node represent the degree value, which further reflects its importance and the number of connections with other nodes.
Enrichment Analysis
The ClusterProfiler software package and R software were used to perform Gene Ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. GO function enrichment analysis involves biological processes (BP), molecular function (MF), and cellular component (CC) biology.
Molecular Docking
The target receptors were retrieved by downloading them from the PDB database (PDB: https://www.rcsb.org/). The compounds were downloaded from PubChem (https://pubchem.ncbi.nlm.nih.gov/). Prior to docking, all receptor proteins were treated with PyMol 2.5, including the removal of water molecules, salt ions, and small molecules. Molecular docking was then accomplished using AutoDock Vina 1.1.2 software, and docking was performed using the default docking parameters. The docking results were visualized using PyMol 2.5.
MCAO Models
All experimental procedures complied with the national guidelines for care and the use of laboratory animals were approved by the Experimental Animal Ethics Committee of Capital Medical University (Ethics Number: aeei-2019-086). All male 8-week-old Sprague–Dawley rats were obtained from Beijing Vital River Laboratory Animal Technology Co., Ltd. (License Number: scxk[Beijing]2016-0006) and raised in the experimental animal center of Capital Medical University in a specific pathogen-free SPF environment (18–22°C, 40–60% humidity, 12 h day/night). Twelve male rats were used to avoid the protective effects of estrogen. The rats were anesthetized using 5% isoflurane (Beyotime, China), and the right common carotid artery (CCA), internal carotid artery (ICA), and external carotid artery (ECA) were gently dissected and exposed. After permanent ligation of the CCA and ECA in the proximal portion, the distal end of the ICA was clamped using hemostatic forceps, and a small incision was made at the bifurcation of the CCA. The preprocessed nylon suture was then inserted from the CCA into the ICA for 18–22 mm until resistance was probed in the middle cerebral artery. Finally, a nylon suture (Beijing Shadong Biotechnology Co., Ltd., Beijing, China) was fixed to prevent spontaneous withdrawal, and the incision was sutured and disinfected. The sham group was treated using the same procedure but without MCAO.
Drug Treatment
The THF dried powder (purity: 5.34 g crude drug/kg dried powder; batch number: 201809-jgf) was provided by the China Academy of Chinese Medical Sciences. The daily consumption of THF is 4.9 g crude drug/kg in rats, which is equivalent to 0.7 g crude drug/kg in patients. We set low dose of THF (THF-L) and high dose of THF (THF-H) groups in rats, approximately one and two times the dosage of clinical patients, respectively: 0.92 g dried powder/kg; 1.84 g dried powder/kg.
Grouping
The animals were randomly assigned to four groups: sham group, model group (MCAO rats), THF-L group (0.92 g dried powder/kg), and THF-H group (1.84 g dried powder/kg). For the THF-L group, 9.2 g of dry powder was weighed and dissolved in 100 ml of deionized water, meaning 10 ml of the solution contained 0.92 g of the drug, and the volume administered to rats was 10 ml/kg/d. For the THF-H group, 18.4 g of dry powder was weighed and dissolved in 100 ml of deionized water, meaning 10 ml of the solution contained 1.84 g of the drug, and the volume administered to rats was 10 ml/kg/d. The rats in the administration group were intragastrically administered 6, 24, 48, and 72 h after cerebral ischemia surgery. The rats in the sham and model groups were administered the same volume of saline (10 ml/kg/day) at the same time points.
Immunofluorescence
For immunofluorescence (Shah et al., 2017), the rats in each group were euthanized, rinsed with warm saline via rapid left ventricular perfusion, and fixed with 4% paraformaldehyde in PBS solution for 60 min. This process was carried out at 3 days and 30 min post injury (dpi) for immunofluorescence studies. The striatal level of the brain was then dissected and fixed with 4% paraformaldehyde in PBS solution for 1 week. After that, the tissues were dehydrated with gradient ethanol, made transparent with xylene, embedded in paraffin, and cut into 0.5-μm sections using a biological tissue microtome. After dewaxing and hydration, the slides were immersed in 5% BSA for 1 h to block nonspecific antibodies. All paraffin sections were then incubated with primary antibodies overnight at 4°C, namely, GFAP (1:1000, rabbit source, Servicebio, Wuhan, China) and C3(1:400, mouse source, Santa Cruz, United States). These sections were then incubated with Alexa Fluor 594 and Alexa Fluor 488 (1:300, goat anti-mouse and goat anti-rabbit source, Abcam, United Kingdom) for 2 h.
Immunohistochemistry
For immunohistochemistry (Kusakawa et al., 2015), the slides were incubated with antibodies against caspase-3 (1:750, rabbit source, Servicebio, Wuhan, China), TNF-α (1:500, rabbit source, Servicebio, Wuhan, China), and IL-6 (1:1000, rabbit source, Servicebio, Wuhan, China) overnight at 4°C. Subsequently, paraffin sections of the tissues were washed with PBS, and the secondary antibody was added and incubated for 1 h.
TUNEL Staining
TUNEL staining (Beyotime Biotechnology, Shanghai, China) was used to detect apoptosis in tissues (Lv et al., 2021). First, the slides were fixed with 4% paraformaldehyde for 30 min at room temperature. The tissue sections were then treated with permeabilization solution for 5 min, followed by TUNEL solution for 1 h at 37°C. Finally, the tissues were observed under a fluorescence microscope (Nikon Eclipse E600, Boston, United States).
Statistical Analysis
Statistical analyses were performed using SPSS 26.0 (IBM Corp., Armonk, NY, United States). All data are expressed as the mean ± standard deviation. Homogeneity of variance was tested; one-way analysis of variance was used to compare multiple groups; and subsequently, the LSD test was performed for pairwise comparison. P < 0.05 was considered to be statistically significant.
Results
Quantitative Analysis of Key Compounds
HL is a monarch medicine, and LQ has the highest dosage in THF. A qualitative study of HL and LQ is important for controlling the quality of the entire prescription. Therefore, berberine, epiberberine, coptisine, palmatine of HL, and phillyrin of LQ were selected as quality control indices according to the standards of the Chinese Pharmacopoeia 2015 edition. HPLC was used to determine the concentrations of the key compounds. The results are shown in Figure 2 and Table 1. The average berberine, epiberberine, coptisine, palmatine, and phillyrin contents were 34.455, 4.266, 7.205, 9.291, and 4.22 mg/g, respectively.
FIGURE 2. Quantitative analysis of key compounds of THF by HPLC analysis. (A,B) HPLC analysis and chemical structures of four major ingredients of HL. 1. Epiberberine (C20H17NO4). 2. Coptisine(C19H13NO4). 3. Palmatine (C21H21NO4). 4. Berberine(C20H17NO4). (C,D) HPLC analysis and chemical structure of phillyrin. 5. Phillyrin(C20H34O11).
Candidate Active Compound of THF
Traditional Chinese medicine is a multitarget, multicomponent therapy model. By screening the TCMSP, ETCM, TCMID, and BATMAN databases, a total of 64 active ingredients were retrieved, of which 42 had regulatory effects on cerebral infarction. DH, HL, LQ, DNX, and DZY contained eight, 10, 16, five, and one chemical compounds/compound, respectively. Beta-sitosterol is a common component of DH and DNX, whereas quercetin is a common component of HL and LQ (Supplementary Table S1).
Common Targets of THF and AIS
A total of 518 targets were obtained by eliminating duplicate targets (Supplementary Table S2). These drug targets were further converted into Uniport ID. A total of 986 AIS-related genes were obtained from the GeneCards and OMIM databases. Finally, 159 common target genes of THF and AIS were identified (Figure 3A).
FIGURE 3. (A) Venn diagram for common genes of THF and AIS. (B) PPI network of the potential targets of THF. (C) Compound–target network of THF on AIS. The regular hexagon, rhombus, and ellipse represent the herbs, herb compounds, and genes, respectively. The number of connected lines indicates the importance of the node.
Protein–Protein Interaction Network
To identify the predominant common targets of THF and AIS, a network analysis of PPI was performed using the STRING database (Figure 3B). A total of 159 common target genes were imported into the STRING database platform, and the PPI network was constructed using Cytoscape 3.8.0 software. The network included 159 nodes and 2,779 edges, with an average degree value of 34.95 (Supplementary Table S3). The larger the node, the more important the role it plays in the network, and the top five targets contain AKT1, TNF, IL6, IL1B, and CASP3 (Table 2), which are not only important genes of THF for treating AIS but are also essential targets for experimental verification.
Compound-Target Network Analysis
To establish the correlation between THF and potential targets, the “Compound–Target Network” was conducted using Cytoscape 3.8.0 (Figure 3C). This network includes 206 nodes and 604 edges, among which the regular hexagon, rhombus, and ellipse represent herbs, herb compounds, and genes, respectively. The degree values, which reflect the number of nodes directly linked to the compound, were used to evaluate the importance of the compound in the network. We selected the top 15 compounds based on their degree values, as shown in Table 3. The results showed that quercetin (MOL000098, degree = 70), coptisine (MOL001458, degree = 38), palmatine (MOL000785, degree = 38), berberine (MOL001454, degree = 36), epiberberine (MOL002897, degree = 36), physcion (MOL000476, degree = 35), emodin (MOL000472, degree = 30), and chrysophanol (MOL001729, degree = 29) were identified. These compounds exhibited high degree values, demonstrating that one herb can target multiple genes and play a crucial role in the treatment of cerebral infarction.
GO and KEGG Pathway Enrichment Analysis
To determine the biological characteristics of potential targets of THF on AIS, GO and KEGG pathway enrichment analyses on 159 common targets were performed using the ClusterProfiler software package of the R platform. A total of 2681 GO functional processes were obtained and annotated in terms of BP, MF, and CC (p <0.05). The top 10 significant terms in the BP, CC, and MF categories are shown in Figure 4A. Among them, the top three enriched GO functions for target genes included cellular response to chemical stress, neuron death, response to oxidative stress in BP, membrane rafts, membrane microdomains, membrane regions in CC, phosphatase binding, protein phosphatase binding, and serine hydrolase activity in MF. KEGG pathway enrichment analysis showed that 168 pathways were associated with cerebral infarction. Detailed pathways related to cerebral infarction are shown in Figure 4B and Supplementary Table S4. Enrichment analysis suggested that the inflammation and apoptotic pathways were enriched by multiple targets. Collectively, these results suggest that restraining inflammation and apoptosis may be a promising therapeutic strategy for AIS treated with THF.
FIGURE 4. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. (A) GO functional enrichment analysis. (B) KEGG pathway enrichment analysis.
Molecular Docking
Berberine (the most abundant component in HL) and phillyrin were used for molecular docking. The key targets IL-6 and TNF-α were also selected for analysis. It is generally accepted that the lower the Vina score, the more stable the binding of the compound is to the crucial target. The results indicated that berberine and phillyrin had good binding activity for IL-6 and TNF-α (Figures 5A–E).
FIGURE 5. Molecular docking. (A) Berberine-(IL-6). (B) Phillyrin-(IL-6). (C) Berberine-(TNF-α). (D) Phillyrin-(TNF-α). (E) Absolute value of vina score (kcal/mol).
Inhibition of THF Activity in Apoptosis in Cerebral Infarction in the MCAO Model
As network pharmacology analysis suggested that THF may protect against neuronal damage by inhibiting the apoptotic pathway, TUNEL staining was used to estimate the number of apoptotic cells at the infarction locus. A marked increase in apoptotic cells was observed in the cortex of MCAO rats, whereas THF-L and THF-H remarkably reduced apoptosis (Figures 6A,C). Immunochemistry was further performed to determine the protective effect of THF on neural cells. The number of caspase-3–positive cells in the cortex was significantly lower than that in the MCAO group. In contrast to MCAO rats, the number of cells with positive expression of caspase-3 in the THF-L and THF-H groups was significantly decreased (Figures 6B,D). These results indicate that THF treatment significantly suppressed the expression of caspase-3 and protected neural cells from apoptosis in MCAO rats.
FIGURE 6. Effects of THF on the neural cell apoptosis in infarction focus at 3 days and 30 min. (A) TUNEL staining of apoptotic cells in the cortex. Scale bar = 50 μm. (B) Representative immunohistochemistry images showing caspase-3–positive cells in the cortex. Scale bar = 100 μm. (C) Quantitative analysis of apoptotic cells in the cortex. (D) Quantification of the number of caspase-3–positive cells in each group. n = 3. (**p < 0.01, compared to the MCAO group).
Inhibitory Activity of THF on the Pro-inflammatory Activation of A1 Astrocyte and Release of IL-6 and TNF-α
In the network pharmacology analysis, IL-6 and TNF were the hub genes of the PPI network, and the TNF signaling pathway was the key pathway in the KEGG enrichment, suggesting that the inflammatory response may play a pivotal role in treating AIS by THF. To evaluate the inflammatory response after cerebral infarction, immunofluorescence was performed to determine the activation of A1 astrocytes that induce an evident inflammatory response in many CNS diseases. First, double immunofluorescence staining of GFAP and C3 was performed to label astrocytes, in which positive cells labeled with C3/GFAP represented A1-type pro-inflammatory astrocytes, whereas positive cells labeled with GFAP represented all types of astrocytes. In this study, a large number of C3-/GFAP-labeled positive cells were detected at the edge of the infarct foci in MCAO rats, but not in the ipsilateral brains of rats in the sham group. Furthermore, the statistical results showed that the cell density of GFAP and C3/GFAP positive cells around infarct lesions in the THF-L and THF-H groups was lower than that in MCAO rats (Figures 7A,C,D). Simultaneously, immunohistochemical staining was performed to evaluate the release of IL-6 and TNF-α. Significant differences were observed between the sham and model groups. Accordingly, IL-6 expression in the cortex was decreased in both the THF-L and THF-H groups compared with that in the MCAO group, whereas THF-H markedly reduced TNF-α expression, suggesting that THF inhibits the release of inflammatory factors (Figures 7B,E,F). Collectively, our study demonstrated that THF curbs the inflammatory response and thereby protects neural cells from apoptosis during cerebral infarction.
FIGURE 7. THF inhibits the activation of A1 astrocytes and expression of IL-6 and TNF-α. (A) Representative immunofluorescence images represent GFAP (red) and C3 (green) staining in the cortex. GFAP, glial fibrillary acidic protein. Scale bar = 50 μm. (B) Immunohistochemistry of IL-6 and TNF-α in the cortex. Scale bar = 100 μm. (C) Quantification of the number of GFAP positive cells in each group. (D) Statistical diagram of GFAP/C3 positive cells. (E) Statistical diagram of TNF-α–positive cells. (F) Quantification of the density of IL-6–immunostained cells. IL-6, interleukin-6; TNF-α, tumor necrosis factor-α. n = 3. (*p < 0.05, **p < 0.01, ***p < 0.001compared to the MCAO group).
Discussion
AIS is a common disease in elderly people, with a stepwise sequence of pathophysiological processes that are strongly correlated to inflammatory responses and apoptosis during development and progression (Ao et al., 2018; Uzdensky, 2019). THF is a traditional herbal medicine formula with multiple active compounds, which has been used in treating AIS for more than 10 years in China owing to its excellent efficacy in “heating–clearing, fire–purgingg, and detoxification.” In this study, network pharmacology and experimental verification were applied to explore the underlying mechanism by which THF improves AIS outcomes.
Network pharmacological analysis was first used to predict and determine the active ingredients and potential targets of THF in AIS. First, 42 potential compounds associated with the therapeutic effects of cerebral infarction were identified. By constructing a PPI network for the common targets of THF and AIS, AKT1, TNF, IL6, IL1B, and CASP3 were predicted as pivotal targets. Finally, KEGG pathway enrichment analysis showed that inflammation- and apoptosis-related pathways were significantly enriched, and we focused on the TNF pathway. IL-6 is a pleiotropic cytokine primarily produced by activated microglia and astrocytes. Upregulation of IL-6 may act on microvascular endothelial cells and increase the expression of adhesion molecules and chemokines, thus mediating the inflammatory cascade and aggravating cerebral ischemic injury (Huang et al., 2006). CASP3 is generally recognized as a crucial molecule in apoptosis after cerebral ischemia. Thus, inhibition of the CASP3 pathway can simultaneously block neural cell apoptosis (Burguillos et al., 2011). TNF-α (tumor necrosis factor–α) is a cytokine associated with systemic inflammation, immune regulation, inhibition of tumor cell growth, and rejection of organ transplantation. There is considerable evidence that TNF-α is significantly associated with ischemic stroke, which further increases sharply in the acute stage of stroke, resulting in serious tissue damage (Shishkova et al., 2018; Zuo et al., 2019). Based on the HPLC analysis and molecular docking results, we speculated that berberine and phillyrin are likely to play a pivotal role in inflammation. Previous studies have shown that berberine alleviates cerebral ischemic injury by attenuating apoptosis, inhibiting M1 glial cell polarization, and inhibiting the expression of pro-inflammatory cytokines (Simões Pires et al., 2014; Zhu et al., 2019). Similarly, in vivo experiments showed that phillyrin inhibits microglia-mediated inflammation by promoting the polarization of “M2” microglia. In addition, phillyrin has been reported to alleviate apoptosis by inhibiting apoptotic pathways such as caspase-3 (Du et al., 2020; Jiang et al., 2021). Altogether, the network pharmacology results showed that TNF, IL6, and CASP3 are hub genes in the PPI network. In addition, KEGG enrichment analysis showed that the TNF inflammatory pathway was significantly enriched through multiple targets (including the hub targets CASP3 and IL6). These targets are actively involved in apoptotic and inflammatory pathways, suggesting that THF may exert anti-inflammatory and antiapoptotic effects on AIS by primarily regulating the TNF signaling pathway and targeting certain genes, such as CASP-3, IL-6, and TNF-α. Therefore, corresponding animal experiments were subsequently designed to validate the predictions of network pharmacology. Specifically, we emphasize that the CASP3 signaling pathway induces apoptosis and detects the activation of A1 astrocytes and the release of IL6 and TNF-α. CASP3, a marker of apoptosis, is located in the core of the apoptotic cascade pathway. Nerve cells are very sensitive to severe ischemia and hypoxia and die rapidly after ischemic stroke. In particular, the number of CASP3-positive neurons in the cerebral cortex and hippocampal CA1 area increased at 3 h, reached a peak at 48 h, and decreased to a stable level 14 days after MCAO (Liu et al., 2013). Therefore, inhibiting the expression of CASP3 and reducing neural cell apoptosis are important therapeutic targets for acute ischemic stroke (Akpan and Troy, 2013). Our findings confirm that THF alleviates apoptosis of neural cells by inhibiting the CASP3 pathway.
Numerous studies have demonstrated that abnormally activated A1 astrocytes exhibit a detrimental effect on the transmission of inflammatory mediators and glutamic acids, initiating immune cascade reactions and leading to increased tissue damage, such as the destruction of the BBB, brain edema, and neuronal degeneration and death (Hawkins and Davis, 2005; Buffo et al., 2010). TNF-α is an important mediator of the inflammatory response, which can activate microglia and astrocytes, leading to leakage of the BBB and aggravating the progression of stroke (Berti et al., 2002; Tuttolomondo et al., 2014). IL-6 is a pro-inflammatory cytokine that promotes the occurrence and development of neurovascular unit dysfunction and neuroinflammation (Pawluk et al., 2020). Therefore, restraining the abnormal activation of A1 astrocytes and production of TNF-α and IL-6 is conducive to alleviating neuronal damage after cerebral ischemia, which is also indispensable for functional rehabilitation. In our study, A1 astrocytes were significantly activated, and the levels of TNF-α and IL-6 were upregulated in MCAO rats, consequently causing neuronal death. However, this inflammation and tissue injury could be reversed by THF treatment. These data suggest a potential mechanism whereby THF can curb overactivation of A1 astrocytes and the release of TNF-α and IL-6 at an early stage after AIS.
Finally, although our study demonstrates that THF has an impact on the inflammatory response and apoptosis by targeting IL-6, TNF-α, and caspase-3 in MCAO rats, the upstream pathway has not been elucidated. In addition, a previous study showed that THF can improve clinical outcomes by altering the gut microbiota (Guo Q. et al., 2021). However, there is currently no evidence of the efficacy of THF in animals. Therefore, the effect of THF on improving animal prognosis and more in-depth mechanisms remain to be elucidated in the future.
Conclusion
To the best of our knowledge, ischemic stroke is a complex pathological process mediated by multiple targets and multiple pathways, and the complicated pathophysiological mechanism of THF in treating AIS was explored based on network pharmacology analysis and experimental validation. According to network pharmacology analysis, 42 bioactive compounds and 159 targets of THF related to AIS were obtained. Simultaneously, AKT1, TNF, IL6, IL1B, and CASP3 were identified as hub targets. Inflammation and apoptosis were significantly enriched in multiple signaling pathways and multiple targets. Therefore, THF is likely to be involved in modulating AIS by primarily regulating the TNF signaling pathway and targeting certain genes such as CASP3, IL-6, and TNF-α. Based on experimental validation, this study provides new insights into the mechanisms by which THF can effectively decrease neural cell apoptosis through the caspase-3 pathway and restrain the excessive activation of A1 astrocytes. Therefore, THF may serve as a prospective treatment strategy for AIS through multiple targets, components, and pathways.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics Statement
The animal study was reviewed and approved by the Ethics Committee of Capital Medical University.
Author Contributions
MQ and LG designed this research; YN and LW analyzed the network pharmacology data, conducted the experiments, and wrote the manuscript; LH performed the preparation of herbal prescriptions; JS, NW and LH contributed to data analysis and interpretation of the results; LG and MQ reviewed the manuscript. All authors have read and approved the manuscript.
Funding
This research was funded by the National Natural Science Foundation of China (81973759, 81573905), the National Key Research and Development Project (2020YFC2003100, 2020YFC2003103), and Capital’s Funds for Health Improvement and Research (SF2020–1–2011).
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.2022.859244/full#supplementary-material
References
Akpan, N., and Troy, C. M. (2013). Caspase Inhibitors: Prospective Therapies for Stroke. Neuroscientist 19 (2), 129–136. doi:10.1177/1073858412447875
Ao, L. Y., Yan, Y. Y., Zhou, L., Li, C. Y., Li, W. T., Fang, W. R., et al. (2018). Immune Cells after Ischemic Stroke Onset: Roles, Migration, and Target Intervention. J. Mol. Neurosci. 66 (3), 342–355. doi:10.1007/s12031-018-1173-4
Berti, R., Williams, A. J., Moffett, J. R., Hale, S. L., Velarde, L. C., Elliott, P. J., et al. (2002). Quantitative Real-Time RT-PCR Analysis of Inflammatory Gene Expression Associated with Ischemia-Reperfusion Brain Injury. J. Cereb. Blood Flow Metab. 22 (9), 1068–1079. doi:10.1097/00004647-200209000-00004
Buffo, A., Rolando, C., and Ceruti, S. (2010). Astrocytes in the Damaged Brain: Molecular and Cellular Insights into Their Reactive Response and Healing Potential. Biochem. Pharmacol. 79 (2), 77–89. doi:10.1016/j.bcp.2009.09.014
Burguillos, M. A., Deierborg, T., Kavanagh, E., Persson, A., Hajji, N., Garcia-Quintanilla, A., et al. (2011). Caspase Signalling Controls Microglia Activation and Neurotoxicity. Nature 472 (7343), 319–324. doi:10.1038/nature09788
Cong, L., Zhao, T., Hui, Z., Xuefen, F., Yun, L., Le, Y., et al. (2021). Protective Effects of TanhuoFormula on Neuroglial Cells and Cerebral Per-Fusion in Rats with Cerebral Ischemia. J. Beijing Univ. Tradit. Chin. Med. 44 (7), 605–614. doi:10.3969/j.issn.1006-2157.2021.07.005
Du, Y., You, L., Ni, B., Sai, N., Wang, W., Sun, M., et al. (2020). Phillyrin Mitigates Apoptosis and Oxidative Stress in Hydrogen Peroxide-Treated RPE Cells through Activation of the Nrf2 Signaling Pathway. Oxid. Med. Cel Longev. 2020, 2684672. doi:10.1155/2020/2684672
Esquiva, G., Grayston, A., and Rosell, A. (2018). Revascularization and Endothelial Progenitor Cells in Stroke. Am. J. Physiol. Cel Physiol. 315 (5), C664–c674. doi:10.1152/ajpcell.00200.2018
Gao, Y., Xie, Y. M., Cai, Y. F., Shen, X. M., Zhao, D. X., and Xie, Y. Z. (2018). Risk factors associated with recurrence within 90 days of ischemic stroke onset in chinese medicine hospital: A national cross-sectional study in China. World J. Trad. Chin. Med. 6 (4), 441–447. doi:10.4103/wjtcm.wjtcm_32_20
Guo, Q., Jiang, X., Ni, C., Li, L., Chen, L., Wang, Y., et al. (2021a). Gut Microbiota-Related Effects of Tanhuo Decoction in Acute Ischemic Stroke. Oxid. Med. Cel Longev. 2021, 5596924. doi:10.1155/2021/5596924
Guo, X., Cao, P., Lian, X., Hu, M., Zhao, J., Shen, W., et al. (2021b). The Neuroprotective Effect of Phillyrin in Intracerebral Hemorrhagic Mice is Produced by Activation of the Nrf2 Signaling Pathway. Eur. J. Pharmacol. 909, 174439. doi:10.1016/j.ejphar.2021.174439
Hawkins, B. T., and Davis, T. P. (2005). The Blood-Brain Barrier/neurovascular Unit in Health and Disease. Pharmacol. Rev. 57 (2), 173–185. doi:10.1124/pr.57.2.4
Herpich, F., and Rincon, F. (2020). Management of Acute Ischemic Stroke. Crit. Care Med. 48 (11), 1654–1663. doi:10.1097/ccm.0000000000004597
Huang, J., Upadhyay, U. M., and Tamargo, R. J. (2006). Inflammation in Stroke and Focal Cerebral Ischemia. Surg. Neurol. 66 (3), 232–245. doi:10.1016/j.surneu.2005.12.028
Huang, C., Zheng, C., Li, Y., Wang, Y., Lu, A., and Yang, L. (2014). Systems Pharmacology in Drug Discovery and Therapeutic Insight for Herbal Medicines. Brief Bioinform. 15 (5), 710–733. doi:10.1093/bib/bbt035
Jiang, Q., Wei, D., He, X., Gan, C., Long, X., and Zhang, H. (2021). Phillyrin Prevents Neuroinflammation-Induced Blood-Brain Barrier Damage Following Traumatic Brain Injury via Altering Microglial Polarization. Front. Pharmacol. 12, 719823. doi:10.3389/fphar.2021.719823
Krishnamurthi, R. V., Ikeda, T., and Feigin, V. L. (2020). Global, Regional and Country-Specific Burden of Ischaemic Stroke, Intracerebral Haemorrhage and Subarachnoid Haemorrhage: A Systematic Analysis of the Global Burden of Disease Study 2017. Neuroepidemiology 54 (2), 171–179. doi:10.1159/000506396
Kusakawa, Y., Mikawa, S., and Sato, K. (2015). BMP5 Expression in the Adult Rat Brain. Neuroscience 284, 972–987. doi:10.1016/j.neuroscience.2014.07.057
Liu, G., Wang, T., Wang, T., Song, J., and Zhou, Z. (2013). Effects of Apoptosis-Related Proteins Caspase-3, Bax and Bcl-2 on Cerebral Ischemia Rats. Biomed. Rep. 1 (6), 861–867. doi:10.3892/br.2013.153
Lv, H., Li, Y., Cheng, Q., Chen, J., and Chen, W. (2021). Neuroprotective Effects against Cerebral Ischemic Injury Exerted by Dexmedetomidine via the HDAC5/NPAS4/MDM2/PSD-95 Axis. Mol. Neurobiol. 58 (5), 1990–2004. doi:10.1007/s12035-020-02223-7
Morovatdar, N., Avan, A., Azarpazhooh, M. R., Di Napoli, M., Stranges, S., Kapral, M. K., et al. (2021). Secular Trends of Ischaemic Heart Disease, Stroke, and Dementia in High-Income Countries from 1990 to 2017: the Global Burden of Disease Study 2017. Neurol. Sci. 43, 255–264. doi:10.1007/s10072-021-05259-2
Patabendige, A., Singh, A., Jenkins, S., Sen, J., and Chen, R. (2021). Astrocyte Activation in Neurovascular Damage and Repair Following Ischaemic Stroke. Int. J. Mol. Sci. 22 (8), 4280. doi:10.3390/ijms22084280
Pawluk, H., Woźniak, A., Grześk, G., Kołodziejska, R., Kozakiewicz, M., Kopkowska, E., et al. (2020). The Role of Selected Pro-inflammatory Cytokines in Pathogenesis of Ischemic Stroke. Clin. Interv. Aging 15, 469–484. doi:10.2147/cia.S233909
Qiu, Y. M., Zhang, C. L., Chen, A. Q., Wang, H. L., Zhou, Y. F., Li, Y. N., et al. (2021). Immune Cells in the BBB Disruption after Acute Ischemic Stroke: Targets for Immune Therapy? Front. Immunol. 12, 678744. doi:10.3389/fimmu.2021.678744
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
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
Saberi, A., Ghayegran, A., Ghayegran, A., Abbasalizade, M., Ehtiatkar, Z., Ghorbani Shirkouhi, S., et al. (2021). The Outcome of Treatment with Recombinant Tissue Plasminogen Activator in Acute Ischemic Stroke. Caspian J. Neurol. Sci. 7 (3), 148–156. doi:10.32598/CJNS.7.26.1
Shah, S. A., Khan, M., Jo, M. H., Jo, M. G., Amin, F. U., and Kim, M. O. (2017). Melatonin Stimulates the SIRT1/Nrf2 Signaling Pathway Counteracting Lipopolysaccharide (LPS)-Induced Oxidative Stress to Rescue Postnatal Rat Brain. CNS Neurosci. Ther. 23 (1), 33–44. doi:10.1111/cns.12588
Shishkova, V. N., Adasheva, T. V., Remenik, A. Y., Valyaeva, V. N., and Shklovsky, V. M. (2018). Prognostic Significance of Clinical-Anthropometric, Biochemical, Metabolic, Vascular-Inflammatory and Molecular-Genetic Markers in the Development of the First Ischemic Stroke. Zh Nevrol Psikhiatr Im S S Korsakova 118 (2), 4–11. doi:10.17116/jnevro2018118214-11
Simões Pires, E. N., Frozza, R. L., Hoppe, J. B., and Menezes, B. d. M. (2014). Berberine was Neuroprotective against an In Vitro Model of Brain Ischemia: Survival and Apoptosis Pathways Involved. Brain Res. 1557, 26–33. doi:10.1016/j.brainres.2014.02.021
Stary, C. M., and Giffard, R. G. (2015). Advances in Astrocyte-Targeted Approaches for Stroke Therapy: an Emerging Role for Mitochondria and microRNAS. Neurochem. Res. 40 (2), 301–307. doi:10.1007/s11064-014-1373-4
Sun, Q., Xu, X., Wang, T., Xu, Z., Lu, X., Li, X., et al. (2021). Neurovascular Units and Neural-Glia Networks in Intracerebral Hemorrhage: from Mechanisms to Translation. Transl. Stroke Res. 12 (3), 447–460. doi:10.1007/s12975-021-00897-2
Tang, C., Hong, J., Hu, C., Huang, C., Gao, J., Huang, J., et al. (2021). Palmatine Protects against Cerebral Ischemia/Reperfusion Injury by Activation of the AMPK/Nrf2 Pathway. Oxid. Med. Cel Longev. 2021, 6660193. doi:10.1155/2021/6660193
Tuttolomondo, A., Pecoraro, R., and Pinto, A. (2014). Studies of Selective TNF Inhibitors in the Treatment of Brain Injury from Stroke and Trauma: a Review of the Evidence to Date. Drug Des. Devel. Ther. 8, 2221–2238. doi:10.2147/dddt.S67655
Uzdensky, A. B. (2019). Apoptosis Regulation in the Penumbra after Ischemic Stroke: Expression of Pro- and Antiapoptotic Proteins. Apoptosis 24 (9-10), 687–702. doi:10.1007/s10495-019-01556-6
Veber, D. F., Johnson, S. R., Cheng, H. Y., Smith, B. R., Ward, K. W., and Kopple, K. D. (2002). Molecular Properties that Influence the Oral Bioavailability of Drug Candidates. J. Med. Chem. 45 (12), 2615–2623. doi:10.1021/jm020017n
Xian, M., Cai, J., Zheng, K., Liu, Q., Liu, Y., Lin, H., et al. (2021). Aloe-emodin Prevents Nerve Injury and Neuroinflammation Caused by Ischemic Stroke via the PI3K/AKT/mTOR and NF-κB Pathway. Food Funct. 12 (17), 8056–8067. doi:10.1039/d1fo01144h
Xu, S., Lu, J., Shao, A., Zhang, J. H., and Zhang, J. (2020). Glial Cells: Role of the Immune Response in Ischemic Stroke. Front. Immunol. 11, 294. doi:10.3389/fimmu.2020.00294
Zhang, R., Zhu, X., Bai, H., and Ning, K. (2019). Network Pharmacology Databases for Traditional Chinese Medicine: Review and Assessment. Front. Pharmacol. 10, 123. doi:10.3389/fphar.2019.00123
Zhao, L., Li, H., Gao, Q., Xu, J., Zhu, Y., Zhai, M., et al. (2021). Berberine Attenuates Cerebral Ischemia-Reperfusion Injury Induced Neuronal Apoptosis by Down-Regulating the CNPY2 Signaling Pathway. Front. Pharmacol. 12, 609693. doi:10.3389/fphar.2021.609693
Zhu, J., Cao, D., Guo, C., Liu, M., Tao, Y., Zhou, J., et al. (2019). Berberine Facilitates Angiogenesis against Ischemic Stroke through Modulating Microglial Polarization via AMPK Signaling. Cell Mol. Neurobiol. 39 (6), 751–768. doi:10.1007/s10571-019-00675-7
Keywords: Tanhuo formula (THF), inflammation, astrocyte, apoptosis, network pharmacology
Citation: Nie Y, Wen L, Li H, Song J, Wang N, Huang L, Gao L and Qu M (2022) Tanhuo Formula Inhibits Astrocyte Activation and Apoptosis in Acute Ischemic Stroke. Front. Pharmacol. 13:859244. doi: 10.3389/fphar.2022.859244
Received: 21 January 2022; Accepted: 15 March 2022;
Published: 26 April 2022.
Edited by:
Fushun Wang, Nanjing University of Chinese Medicine, ChinaReviewed by:
Wenqiang Chen, Georgetown University, United StatesXiaoyan Gao, Beijing University of Chinese Medicine, China
Copyright © 2022 Nie, Wen, Li, Song, Wang, Huang, Gao and Qu. 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: Miao Qu, cXVtaWFvdGNtQDEyNi5jb20=; Li Gao, eHVhbnd1Z2FvbGlAMTI2LmNvbQ==
†These authors have contributed equally to this work