Skip to main content
Download PDF

Exploring the protective mechanism of baicalin in treatment of atherosclerosis using endothelial cells deregulation model and network pharmacology

BMC Complementary Medicine and Therapies volume 22, Article number: 257 (2022) Cite this article

Abstract

Background

Baicalin is a generally available flavonoid with potent biological activity. The present study aimed to assess the underlying mechanism of baicalin in treatment of atherosclerosis (AS) with the help of network pharmacology, molecular docking and experimental validation.

Methods

The target genes of baicalin and AS were identified from public databases, and the overlapping results were considered to be baicalin-AS targets. Core target genes of baicalin were obtained through the PPI network and validated by a clinical microarray dataset (GSE132651). Human aortic endothelial cells (HAECs) were treated with Lipopolysaccharide (LPS) to construct an endothelial injury model. The expression of NOX4 was examined by real-time qPCR and western blot. Flow cytometry was used to detect intracellular levels of reactive oxygen species (ROS). Furthermore, HAECs were transfected with NOX4-specific siRNA and then co-stimulated with baicalin and LPS to investigate whether NOX4 was involved in the anti-oxidative stress effects of baicalin.

Results

In this study, baicalin had 45 biological targets against AS. Functional enrichment analysis demonstrated that most targets were involved in oxidative stress. Using the CytoHubba plug-in, we obtained the top 10 genes in the PPI network ranked by the EPC algorithm. Molecular docking and microarray dataset validation indicated that NOX4 may be an essential target of baicalin, and its expression was significantly suppressed in AS samples compared to controls. In endothelial injury model, intervention of HAECs with baicalin increased the expression levels of NOX4 and NOS3 (eNOS), and decreased LPS-induced ROS generation. After inhibition of NOX4, the anti-ROS-generating effect of baicalin was abolished.

Conclusion

Collectively, we combined network pharmacology and endothelial injury models to investigate the anti-AS mechanism of baicalin. The results demonstrate that baicalin may exert anti-oxidative stress effects by targeting NOX4, providing new mechanisms and insights to baicalin for the treatment of AS.

Peer Review reports

Background

Atherosclerosis (AS), the leading cause of cardiovascular disease (CVD), is difficult to detect in early stages until lipid-burdened plaque ruptures or erodes, eventually leading to thrombosis and tissue damage [1]. Endothelial cells (ECs) are the first line of defense for vascular health. Both activation of ECs and infiltration of circulating monocytes are crucial factors in occurrence and progression of AS. Inflammatory ECs recruit leukocytes by expressing adhesion molecules; circulating monocytes infiltrate into damaged endothelium and differentiate into macrophages, which are further transformed into foam cells in response to excessive lipid load. This process promotes the release of pro-inflammatory factors and generation of reactive oxygen species (ROS) [2]. ROS, mainly originating from the NADPH oxidase (NOX) complex, mediates oxidative stress that promotes the progression of AS by inducing inflammation, apoptosis and alterations in vascular tone [3]. Increased expression of NOX1, NOX2 and NOX5 has been reported as contributing to AS and vascular diseases [4, 5]. However, NOX4 is a protective gene in ECs that exerts anti-oxidative effects by increasing the release of NO through AKT-dependent phosphorylation [6, 7].

Chinese herbal medicine is the most widely used form of complementary and alternative medicine, which has been practiced for centuries in Asia. In treatment of coronavirus disease 2019, the effects were remarkable [8,9,10,11]. Studies on the components of herbal medicine have shown that baicalin is mainly found in the roots of Scutellaria baicalensis, but also in Salvia miltiorrhiza, Platycodon grandiflorum, Paeonia lactiflora and Sempervivum [12,13,14,15]. As a representative of the flavonoids, baicalin has been shown to have a variety of pharmacological functions, such as reducing oxidative stress and inflammation through the NF-κB pathway, blocking the growth cycle of vascular smooth muscle cells, and inhibiting intimal proliferation after injury [16].

Network pharmacology, first proposed by Andrewl Hopkins in 2007 [17], combines pharmacology with bioinformatics to reveal specific targets for drug interventions in disease process, helping to advance the development of precision medicine [18]. It also can be used to analyze the mechanism of herbal formula, such as BuYangHuanWu decoction [19] and Chaihu Lizhong Tang [20]. In current study, we used multiple public databases to predict baicalin-AS targets and set up pharmacological network. Through molecular docking and microarray dataset (GSE132651) validation, we hypothesized that NOX4 might be an essential target of baicalin in treatment of AS. Thereafter, endothelial injury model and NOX4-specific siRNA were constructed to investigate whether NOX4 was involved in the anti-oxidative stress effects of baicalin. The flow chart was exhibited in Fig. 1.

Fig. 1
figure 1

The flowchart for exploring the molecular mechanisms of baicalin in treatment of AS using network pharmacology and experimental validation. AS: atherosclerosis

Full size image

Methods

Predicting target genes of baicalin

Firstly, we obtained the pharmacokinetics parameters of baicalin by TCMSP database [21]. The parameters (oral bioavailability = 40.12%, druglike = 0.75) indicated that baicalin has positive pharmacodynamic effects [22]. Following this, further predictions including absorption, distribution, metabolism, excretion and toxicity (ADMET) were performed using the admetSAR 2.0 server regression model [23]. The prediction results were presented in Supplementary Files (Table S1). Subsequently, we queried the PubChem database (https://pubchem.ncbi.nlm.nih.gov/) for the 2D structures, molecular formulas, and canonical SMILES strings. SDF file of baicalin was entered into three public databases for predicting potential targets including Similarity ensemble approach (SEA) [24], Swiss Target Prediction [25]. and Gene Cards [26]. In the SEA database, the top 100 results ranked by P-value were retained; for the Swiss Target Prediction, results with probability greater than 0.1 were selected; while among the Gene cards, those with score greater than 0.1 were reserved. For each database, species were restricted to ‘Homo sapiens’. After pooling the results, only genes predicted by two or three databases at the same time were screened as high confidence target genes.

Collection of target genes for AS

Similarly, to identify key genes associated with the development of AS, we searched for the keyword “atherosclerosis” in multiple databases(i.e., DisGeNET [27], Gene Cards and DrugBank [28]) and collected overlapping results as vital pathogenic targets of AS for further study.

Identification of baicalin-AS targets

The target genes of baicalin were matched with AS-related targets, and overlapping data were considered to be potential targets of baicalin for AS treatment. Venn diagrams were generated in R software (version 3.6.2) using “VennDiagram” package to visualize the number of intersecting genes.

Functional enrichment analysis of baicalin-AS targets

Gene Ontology (GO) and KEGG can identify enriched biological functions and pathways [29,30,31,32]. To obtain functional annotation information of baicalin-AS targets, we performed Gene Ontology (GO) and KEGG analysis by R packages including “clusterProfiler” [33] and “ggplot2” [34]. The results of enrichment analysis with statistical significance (adjusted p-value < 0.05) were screened and then visualized by bubble plots and multi-level interaction networks.

Protein–protein interaction (PPI) analysis and central module identification

The STRING database (https://string-db.org/, version 10.5) [35, 36] was used to build PPI network of baicalin-AS targets. To make the network more concise and with higher confidence, we kept nodes and edges with interaction score greater than 0.4 and hid the disconnected nodes. Thereafter, the PPI information was aggregated into a “TSV” format file and imported into Cytoscape software (version 3.8) [37] for visualization. The Edge Percolation Component (EPC) algorithm of the Cytohubba plug-in in Cytoscape was applied to calculate the score of each node and determine the top 10 hub genes based on their ranking.

Molecular docking

The SDF file for baicalin was downloaded from the PubChem database and converted to PDB format by OpenBabel. The 3D crystal structures of the five hub target genes were downloaded from the RCSB Protein Data Bank (http://www.pdb.org/). The ligand and receptor files were imported into AutoDock Tools 1.5.6 and prepared for docking, including dehydration and hydrogenation, which were exported and saved in PDBQT format. Then, the appropriate parameters were set and the ‘Grid box’ was set to maximum for blind docking. A binding energy of less than -5 kJ∙mol-1 indicates that the ligand can spontaneously bind to the receptor. Docking results were visualized using PyMOL.

Probing the expression of hub genes in endothelial injury model

To examine the gene expression patterns (up- or down-regulation) of hub genes in endothelial injury model, we selected GSE132651 dataset, a microarray study that quantifies gene expression of ECs in living subjects, for further analysis. This dataset consisted of 19 samples, including 6 normal controls and 13 subjects with abnormal endothelial function (showing a risk for early AS). We retrieved original matrix files from GEO repository (https://www.ncbi.nlm.nih.gov/geo/), filtered the mRNA expression levels of hub genes between two groups using “Limma” package [38] and visualized by GraphPad Prism 7.0.

Reagents

Baicalin (purity > 98%) was supplied by Yuan-ye biotechnology (Shanghai, China). Lipopolysaccharide (LPS), DMSO, and cell counting kit-8 (CCK-8) were obtained from Meilun biology (Dalian, China). EGM-2 Bullet Kit were obtained from Lonza (#CC-3162, Basel, Switzerland), including endothelial cell basal medium, FBS, VEGF, etc. RNA-quick purification kit was purchased from Yi-shan biotechnology (Shanghai, China). First strand cDNA synthesis kit and SYBR Green master mix were obtained from Thermo Scientific (Shanghai, China). PCR primers were synthesized by Sangon biotechnology (Shanghai, China). DCFDA/H2DCFDA cellular ROS assay kit was obtained from Abcam (Massachusetts, USA).

Cell culture

Human aortic endothelial cells (HAECs) were purchased from Chinese Tissue Culture Collections (CTCC). Cells were cultured in MEM medium (supplemented with 2% FBS, 2 mM L-glutamine, 1000 units/mL penicillin, 100 μg/mL streptomycin, and 0.25 μg/mL amphotericin B) and maintained in a humidified incubator at 37℃ with 5% CO2. For experimental group, cells were pretreated with optimal concentrations of baicalin for 12 h, and then stimulated with LPS (1 μg/mL) for 6 h.

Cytotoxicity assay of baicalin

The cytotoxicity of baicalin was tested by CCK-8 assay according to the manufacturer’s instructions. Briefly, to obtain optimal dose of intervention, HAECs (8 × 104/ml) were seeded in 96-well culture plates and exposed to gradually increasing concentrations (0-400 μM) of baicalin for 12 h. The absorbance value of each well was measured at 450 nm using EnSpire multifunctional microplate reader (PerkinElmer, USA). Cell viability at different concentrations was calculated from the absorbance of the blank (without HAECs) and control (baicalin = 0 μM) wells.

Reverse transcription quantitative PCR (RT-qPCR)

Total RNA was extracted from HAECs and its purity was checked using Nanodrop (Thermo Scientific, USA). After that, mRNA was transcribed into the first-strand cDNA and subjected to polymerase chain reaction on 7500 real-time PCR system (Applied Biosystems, USA). We chose GAPDH as endogenous control and compared gene expression based on the 2—△△CT algorithm. The primer sequences used in current study were listed in Table 1.

Table 1 RT-qPCR primers used in this study
Full size table

Western blot

Total proteins were extracted from HAECs using RIPA lysis buffer (Beyotime Biotechnology, China), followed by concentration detection with the BCA protein assay kit (Sangon Biotech, China). Equal amounts of protein were separated by 10% SDS-PAGE and transferred onto PVDF membranes. Afterwards, PVDF membranes were blocked with 5% skim milk for 1 h and then incubated with primary antibodies including β-actin (1:1000, Cell Signaling Technology (CST), #3700), NOX4(1:1000, abcam, #ab133303) at 4 °C overnight. The membranes were incubated with secondary antibodies (Beyotime Biotechnology, Shanghai, China) for 2 h. Differences in gray values were calculated using the Image Lab software (Bio-Rad Laboratories Inc). The experiment was repeated three times independently.

Transfection

For small interfering RNA (siRNA) transfection, HAECs were transfected with NOX4-specific siRNA (Tsingke Biotechnology, Beijing, China) or negative control siRNA (Table S2) for 6 h by Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA) with Opti-MEM Reduced-Serum Medium (Gibco). Then, the medium was replaced with regular culture medium containing serum and the cells were cultured for 12 h. Next, the cells were exposed to LPS (1 μg/mL) with or without baicalin. RT-qPCR and western blot analysis were used to evaluate the efficiency of siRNA.

Intracellular ROS measurement

Briefly, HAECs pretreated with baicalin for 12 h were harvested and stained with 20 μM DCFDA in culture medium. After incubation at 37℃ for 30 min (protected from light), the cells were washed with 1X buffer and treated with 1 μg/ml LPS for 6 h. Then, cells were gently pipetted to make single cell suspension and analyzed on flow cytometer (BD Accuri C6, USA). Data processing was carried out on FlowJo software (version 10, USA).

Statistical analysis

Statistical analysis was performed in R software or GraphPad Prism 7.0 (California, USA). When analyzed using GraphPad Prism, the data were shown as mean ± SEM. The t-test was used for comparisons between two groups; one-way analysis of variance was taken for comparisons between three or more groups. Results with p-value < 0.05 were considered significant.

Results

Collection target genes of baicalin

The skeletal formula of baicalin was created using ACD/ChemSketch software (version 12.0, Advanced Chemistry Development, Toronto, Canada) (Fig. 2A). We collated the target genes for baicalin predicted by the SEA, Swiss Target Prediction and Gene Cards databases. Based on the screening criteria mentioned in “Methods”, we obtained 120 targets with high confidence for baicalin.

Fig. 2
figure 2

Identification of baicalin-AS target genes. A The chemical structure of baicalin. B Venn diagram showing the overlap among target genes of baicalin (blue) and AS related target genes (yellow). A total of 45 gene products were identified as common

Full size image

Identification of AS-related target genes

By integrating disease targets derived from the DisGeNET, Gene Cards, and DrugBank databases, we ultimately retained 497 genes that appeared in at least two databases. Admittedly, the pathophysiology of AS is complicated, with an extensive number of genes involved and interacting with each other.

Baicalin-AS target genes

We obtained intersection data by matching targes of baicalin with AS-related target genes. As shown in Venn diagram (Fig. 2B and Table S3), 45 baicalin-AS target genes were generated, which were potential targets of baicalin for AS treatment.

Functional enrichment analysis of baicalin-AS targets

The potential biological processes (BP) and signaling pathways related to baicalin-AS targets were explored through functional enrichment analysis. The top 10 possible BP items ordered by “gene ratio” were shown in Fig. 3A, including responses to oxidative stress, cellular responses to oxidative stress, metabolic processes of reactive oxygen species (ROS), etc. On the other hand, KEGG pathway analysis showed that baicalin-AS targets were mainly involved in human cytomegalovirus infection, IL-17 signaling pathway, fluid shear stress and atherosclerosis, etc. (Fig. 3B). Overall, enrichment analysis showed that most of the target genes were associated with inflammation, suggesting that baicalin may act against AS through anti-inflammatory effects. We sorted out the relationships between baicalin, target genes and signaling pathways in Excel, and then imported the information into Cytoscape software to construct a multi-level interaction network. As shown in the Fig. 4, numerous genes were involved in inflammatory and immune-related pathways. Genes involved in fluid shear stress and atherosclerosis include IL1β, TNF, MMP9, AKT1, and NOX4; genes associated with AGE-RAGE signaling pathway in diabetic complications include AKT1, IL1β, and NOX4; while genes associated with HIF signaling pathway contain STAT3, NOX4, NFκB, and MAPK1.

Fig. 3
figure 3

Functional enrichment analysis of baicalin-AS target genes. A Top 10 terms of BP enrichment analysis for baicalin-AS target genes. B Top 10 terms of KEGG enrichment analysis for baicalin-AS target genes

Full size image
Fig. 4
figure 4

Compound-target-pathway networks(C-P–T). Green squares represented compounds and pink circles represented targets. The blue V shapes represented the pathways

Full size image

PPI network construction and central module identification

To further explore the association between the 45 target genes, we constructed a PPI network, containing 44 nodes and 542 edges, via STRING database (Fig. 5A). In parallel, the top 10 target genes (IL6, NOX4, TNF, STAT3, MMP9, IL1B, AKT1, CASP3, PTGS2 and HIF1A) were calculated using the cytoHubba plug-in and shown in Fig. 5B, which were considered to be hub target genes for baicalin in treatment of AS.

Fig. 5
figure 5

NOX4 was the core target gene of baicalin. A PPI network topological analysis. Based on gene intersection, this study acquired altogether 44 protein nodes and 10 hub target genes from topological analysis. As shown, the node color represented the score, with deeper and redder nodes meaning higher degree value. B The top 10 hub genes were screened from (A) using CytoHubba. As shown, the node color represented the score, with deeper and redder nodes meaning higher EPC scores. C Molecular docking of NOX4 with baicalin, affinity =  − 9.31 kcal/mol. D NOX4 mRNA expression in GSE132651 based on probe intensity. (*P < 0.05)

Full size image

Verification of ligand–receptor interaction

In order to select the most powerful binding molecules, we selected the top 5 target genes in the hub module for molecular docking with baicalin. It was confirmed that baicalin had a high binding activity with the hub genes, and the binding energy to NOX4 is the strongest (Table 2). Thus, it was hypothesized that NOX4 might be an essential target of baicalin in the treatment of AS (Fig. 5C).

Table 2 Results of molecular docking of hub target genes and baicalin
Full size table

NOX4 expression was down-regulated in endothelial injury condition

To further investigate the role of NOX4 in AS, we selected the GSE132651 dataset for validation. The details of gene expression matrix were shown in Supplementary Files. As illustrated in Fig. 5D, inter-group comparison showed that NOX4 expression was significantly suppressed in AS sample group (probe intensity, 176.8 ± 19.42) compared to controls (probe intensity, 285.1 ± 56.56) (P < 0.05).

Cytotoxic effects of baicalin on HAECs

The cytotoxic effects of multiple concentration gradients of baicalin on HAECs were examined by CCK-8 assay. The results showed that baicalin had no obvious cytotoxicity to HAECs when the treatment concentration was less than 30 μM (Fig. 6A). Therefore, we chose the maximum safe concentration (30 μM) for subsequent experiments.

Fig. 6
figure 6

Baicalin increased the expression levels of NOX4 and decreased the LPS-induced ROS generation. A Effect of different concentrations of baicalin on viability of HAECs assessed by CCK8. B Baicalin modulated mRNA expression level of NOX4 and eNOS. GAPDH was used as the control. C HAECs were cultured with LPS supplemented with baicalin or not. Cell lysates were then analyzed using Western blotting assay. Each histogram value is the ratio of the scan value for NOX4 divided by the value for β-actin. The values represent mean ± SEM of three independent experiments. D Baicalin inhibited ROS release. HAECs were exposed to baicalin for 12 h and stained with 20 μM DCFDA for 30 min, followed by 1 μg/ml LPS for 6 h. ROS levels were detected by flow cytometry analysis. Data were shown as mean ± SEM. (n = 3, *P < 0.05, **P < 0.01, ***P < 0.001, ***** P < 0.0001)

Full size image

Baicalin increased NOX4 and eNOS expression in LPS-stimulated HAEC

In the GO analysis section, we noted that the enrichment of “oxidative stress” and “ROS production” was very high, suggesting that baicalin may have a modulating effect on oxidative stress and ROS. To further probe whether baicalin has a protective effect on AS, we established an LPS-induced endothelial injury model for in vitro experiments [39, 40].

NOX4 is one of the hub targets of baicalin, which can inhibit oxidative stress by promoting the production of NO [41]. NOS3 (eNOS) can interact with L-arginine to promote the physiological synthesis of nitric oxide (NO) [42]. In the LPS-induced endothelial cell dysfunction model, mRNA expression of NOX4 and eNOS was decreased compared to the control group, which was notably reversed after 12 h of treatment with baicalin. (Fig. 6B). To further demonstrate the relationship between baicalin and NOX4, we examined NOX4 protein levels by western blot, which showed that baicalin also elevated NOX4 protein levels in HAECs (Fig. 6C). These data suggested that baicalin may exert protective effects on LPS-induced HAECs by targeting NOX4 and eNOS, thereby alleviating oxidative stress levels.

Baicalin attenuated ROS release from LPS-stimulated HAECs

Excessive intracellular ROS cause endothelial dysfunction and oxidative stress injury. In this study, intracellular ROS levels were measured by DCFH-DA staining to assess whether baicalin was involved in regulating ROS in HAECs. Flow analysis demonstrated that ROS levels were increased significantly in the LPS-stimulated group compared to the control group. And the ROS levels were down-regulated after baicalin treatment compared to the LPS-stimulated group. The results showed that baicalin could attenuate ROS release in HAECs (Fig. 6D).

NOX4 plays an essential role in the anti-oxidative stress of baicalin

To further confirm whether the anti-oxidative stress effects of baicalin was regulated by NOX4, three pairs of NOX4 siRNAs were separately transfected into HAECs to knock down endogenous NOX4. The reduction of NOX4 mRNA was 67.24%, 60.04% and 63.14% after transfection with siRNA-1, siRNA-2 or siRNA-3, respectively (Fig. 7A). Since NOX4 protein expression was significantly weaker in HAECs transduced with siRNA-3, we selected siRNA-3 for follow-up experiments (Fig. 7B). As shown in Fig. 7C, the ROS content in NOX4-siRNA did not change significantly compared to the control-siRNA. After LPS stimulation, ROS levels were significantly higher in the NOX4-siRNA group whereas the baicalin slightly reduced ROS but not statistically (p = 0.2647).

Fig. 7
figure 7

Low expression of NXO4 attenuated the inhibition of LPS-induced ROS production by baicalin. A NOX4 mRNA expression in HAECs infected with control-siRNA, NOX4-siRNA-1, NOX4-siRNA-2 and NOX4-siRNA-3. All NOX4-siRNAs remarkably downregulated NOX4 mRNA. B Western blotting analysis of the expression of NOX4 after the suppression of NOX4 expression. C The expression of ROS content in HAECs transfected with NOX4-siRNA or control-siRNA after stimulation with or without LPS (1 μg/ml). ROS levels were detected by flow cytometry analysis. Data were shown as mean ± SEM. (n = 3, *P < 0.05, ***P < 0.001)

Full size image

Discussion

Chinese medicine is widely used in clinical treatment in Asia. It has gradually shown its advantages owing to multi-target characteristics and positive therapeutic effects. However, the application of TCM is obstructed by its complex composition and unrecognized therapeutic mechanisms. With the flourishing development of bioinformatics and network pharmacology, it is possible to elucidate active ingredients and pharmacological mechanisms of Chinese herbal medicine. The preventive and therapeutic effects of herbs, such as Salvia miltiorrhiza, Rhizoma chuanxiong and Scutellaria baicalensis, on AS have been widely confirmed [43,44,45]. Scutellaria baicalensis has immunomodulatory effects on macrophages and dendritic cells, regulating the shift in cell polarization toward anti-inflammation and producing a preventive effect on foam cell formation, thereby inhibiting the development of AS [46]. Baicalin is a flavonoid derived from Scutellaria baicalensis Georgi and has antibacterial, anti-inflammatory, antioxidant, anticancer, and other functions [47, 48]. Li et al. [49] observed that baicalin inhibited cigarette smoke-induced inflammation through suppressing the phosphorylation of HDAC2. Sun et al. [50] found that baicalin could reduce the lipopolysaccharide (LPS)-induced periodontitis in rats by curbing the TLR/MYD88/NF-κB signaling pathway. However, the target genes and mechanisms of baicalin in alleviating oxidative stress and thus anti-AS remains unclear.

Oxidative modifications of lipids have been detected in vascular lesions and the degree of oxidation positively correlates with the severity of disease, suggesting a role of oxidative stress in AS [51], while the NADPH oxidase (NOX) family is an important mediator of ROS production in the vascular wall. NOX1 and NOX2 have been shown to be deleterious to vascular disease development, particularly in the context of AS. In contrast, NOX4, the major NADPH oxidase in endothelial cells, has been shown to protect against AS in several mouse models [52,53,54]. NOX4 mainly releases H2O2 to stimulate the expression of NOS and the formation of NO, which inhibits oxidative stress by scavenging oxygen free radicals [55, 56]. Schürmann et al. demonstrated that NOX4 is an endogenous anti-atherosclerotic enzyme and that genetic deletion of NOX4 accelerates the development of atherosclerosis in ApoE-/- mice. NOX4 deficiency in endothelial cells increases leukocyte adhesion and exacerbates endothelial inflammatory activation [57]. eNOS promotes the formation of NO in the vascular endothelium. ROS are spatially close to eNOS and can directly provoke the uncoupling of eNOS, reducing NO bioavailability and activating other redox-sensitive systems [58]. Oxidative stress was considered as the major contributor to eNOS uncoupling and endothelial dysfunction. It was also known that ROS impair endothelial NO-mediated coronary micro-vessel dilation [59].

In our study, the core gene of baicalin-AS targets was identified through network pharmacology, molecular docking, and microarray dataset analysis. As a result, NOX4 was shown to be down-regulated under endothelial damage conditions. In addition, functional enrichment analysis showed that baicalin-AS targets were associated with inflammation and oxidative stress pathways, suggesting that baicalin may act against AS via anti-inflammatory and anti-oxidative stress effects. Through literatures, we found that NOX4 was closely associated with oxidative stress. To further investigate whether baicalin is associated with reducing oxidative stress by targeting NOX4, we conducted HAECs-based molecular experiments, including RT-qPCR, western blot and flow cytometry. Interestingly, we found that baicalin elevated the expression of NOX4 and eNOS, and reduced intracellular ROS levels in LPS-induced endothelial injury model. Nevertheless, the inhibitory effect of baicalin on ROS production was diminished after down-regulation of NOX4 expression. Consequently, we conclude that baicalin may exert its pharmacological function by targeting NOX4 to alleviate oxidative stress.

Some limitations of this study should be noted. First, three different databases were used in order to obtain adequate targets for baicalin. However, some bias may exist due to the different computer algorithms used in each database. Secondly, our conclusions were only validated on one kind of endothelial cells and not tested on other endothelial cells and animals. In the next studies, other types of endothelial cells and Nox4−/− mice should be added to scientifically and rigorously demonstrate the in-depth mechanism of baicalin in the alleviation of AS.

Conclusion

Collectively, we combined network pharmacology and endothelial injury models to investigate the anti-AS mechanism of baicalin. The results demonstrate that baicalin may exert anti-oxidative stress effects by targeting NOX4, providing new mechanisms and insights to baicalin for the treatment of AS.

Availability of data and materials

The dataset analyzed during the current study are available in the GEO repository, GSE132651. Other data used to support the findings of this study are available from the corresponding author upon request.

References

  1. Benjamin EJ, Muntner P, Alonso A, Bittencourt MS, Callaway CW, Carson AP, et al. Heart Disease and Stroke Statistics-2019 update: a report from the American Heart Association. Circulation. 2019;139(10):e56–528.

    Article  PubMed  Google Scholar 

  2. Frostegard J. Immunity, atherosclerosis and cardiovascular disease. BMC Med. 2013;11:117.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Bedard K, Krause KH. The NOX family of ROS-generating NADPH oxidases: physiology and pathophysiology. Physiol Rev. 2007;87(1):245–313.

    Article  CAS  PubMed  Google Scholar 

  4. Wind S, Beuerlein K, Armitage ME, Taye A, Kumar AH, Janowitz D, et al. Oxidative stress and endothelial dysfunction in aortas of aged spontaneously hypertensive rats by NOX1/2 is reversed by NADPH oxidase inhibition. Hypertension. 2010;56(3):490–7.

    Article  CAS  PubMed  Google Scholar 

  5. Guzik TJ, Chen W, Gongora MC, Guzik B, Lob HE, Mangalat D, et al. Calcium-dependent NOX5 nicotinamide adenine dinucleotide phosphate oxidase contributes to vascular oxidative stress in human coronary artery disease. J Am Coll Cardiol. 2008;52(22):1803–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Takac I, Schröder K, Zhang L, Lardy B, Anilkumar N, Lambeth JD, et al. The E-loop is involved in hydrogen peroxide formation by the NADPH oxidase Nox4. J Biol Chem. 2011;286(15):13304–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Miura H, Bosnjak JJ, Ning G, Saito T, Miura M, Gutterman DD. Role for hydrogen peroxide in flow-induced dilation of human coronary arterioles. Circ Res. 2003;92(2):e31-40.

    Article  CAS  PubMed  Google Scholar 

  8. Zhang L, Yu J, Zhou Y, Shen M, Sun L. Becoming a Faithful Defender: Traditional Chinese Medicine against Coronavirus Disease 2019 (COVID-19). Am J Chin Med. 2019;2020:1–15.

    Google Scholar 

  9. Luo E, Zhang D, Luo H, Liu B, Zhao K, Zhao Y, et al. Treatment efficacy analysis of traditional Chinese medicine for novel coronavirus pneumonia (COVID-19): an empirical study from Wuhan, Hubei Province, China. Chin Med. 2020;15:34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Serafini M, Peluso I, Raguzzini A. Flavonoids as anti-inflammatory agents. Proc Nutr Soc. 2010;69(3):273–8.

    Article  CAS  PubMed  Google Scholar 

  11. Wang ZY, Jiang ZM, Xiao PT, Jiang YQ, Liu WJ, Liu EH. The mechanisms of baicalin ameliorate obesity and hyperlipidemia through a network pharmacology approach. Eur J Pharmacol. 2020;878:173103.

    Article  CAS  PubMed  Google Scholar 

  12. Ma X, Shao S, Xiao F, Zhang H, Zhang R, Wang M, et al. Platycodon grandiflorum extract: chemical composition and whitening, antioxidant, and anti-inflammatory effects. RSC Adv. 2021;11(18):10814–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Gharari Z, Bagheri K, Derakhshani B, Sharafi A. HPLC-DAD-ESI/MS(n) analysis of phenolic components of Scutellaria araxensis, S. bornmuelleri and S. orientalis. Nat Prod Res. 2022;36(9):2440–5.

    Article  CAS  PubMed  Google Scholar 

  14. Gao Q, Yu G, Yu M, Wang X. Bupleuri radix prevents the recurrences of resected colonic polyps by affecting Angiogenin-2-Induced Protein Kinase B/Akt signaling. J Oncol. 2020;2020:3531652.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Shi JW, Li ZZ, Wu JS, Jin WY, Chang XY, Sun H, et al. Identification of the bioactive components of Banxia Xiexin Decoction that protect against CPT-11-induced intestinal toxicity via UPLC-based spectrum-effect relationship analyses. J Ethnopharmacol. 2021;266:113421.

    Article  CAS  PubMed  Google Scholar 

  16. Dong LH, Wen JK, Miao SB, Jia Z, Hu HJ, Sun RH, et al. Baicalin inhibits PDGF-BB-stimulated vascular smooth muscle cell proliferation through suppressing PDGFRβ-ERK signaling and increase in p27 accumulation and prevents injury-induced neointimal hyperplasia. Cell Res. 2010;20(11):1252–62.

    Article  CAS  PubMed  Google Scholar 

  17. Hopkins AL. Network pharmacology. Nat Biotechnol. 2007;25(10):1110–1.

    Article  CAS  PubMed  Google Scholar 

  18. Boezio B, Audouze K, Ducrot P, Taboureau O. Network-based approaches in pharmacology. Mol Inform. 2017;36(10):1–10.

    Article  Google Scholar 

  19. Liu B, Song Z, Yu J, Li P, Tang Y, Ge J. The atherosclerosis-ameliorating effects and molecular mechanisms of BuYangHuanWu decoction. Biomed Pharmacother. 2020;123:109664.

    Article  CAS  PubMed  Google Scholar 

  20. Zhang M, Yuan Y, Zhou W, Qin Y, Xu K, Men J, et al. Network pharmacology analysis of Chaihu Lizhong Tang treating non-alcoholic fatty liver disease. Comput Biol Chem. 2020;86:107248.

    Article  CAS  PubMed  Google Scholar 

  21. Ru J, Li P, Wang J, Zhou W, Li B, Huang C, et al. TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J Cheminform. 2014;6:13.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Ertl P, Rohde B, Selzer P. Fast calculation of molecular polar surface area as a sum of fragment-based contributions and its application to the prediction of drug transport properties. J Med Chem. 2000;43(20):3714–7.

    Article  CAS  PubMed  Google Scholar 

  23. Yang H, Lou C, Sun L, Li J, Cai Y, Wang Z, et al. admetSAR 2.0: web-service for prediction and optimization of chemical ADMET properties. Bioinformatics. 2019;35(6):1067–9.

    Article  CAS  PubMed  Google Scholar 

  24. Keiser MJ, Roth BL, Armbruster BN, Ernsberger P, Irwin JJ, Shoichet BK. Relating protein pharmacology by ligand chemistry. Nat Biotechnol. 2007;25(2):197–206.

    Article  CAS  PubMed  Google Scholar 

  25. Daina A, Michielin O, Zoete V. SwissTargetPrediction: updated data and new features for efficient prediction of protein targets of small molecules. Nucleic Acids Res. 2019;47(W1):W357–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Stelzer G, Rosen N, Plaschkes I, Zimmerman S, Twik M, Fishilevich S, et al. The genecards suite: from gene data mining to disease genome sequence analyses. Curr Protoc Bioinformatics. 2016;54:1.30.1-1..3.

    Article  Google Scholar 

  27. Pinero J, Ramirez-Anguita JM, Sauch-Pitarch J, Ronzano F, Centeno E, Sanz F, et al. The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 2020;48(D1):D845–55.

    CAS  PubMed  Google Scholar 

  28. Wishart DS, Feunang YD, Guo AC, Lo EJ, Marcu A, Grant JR, et al. DrugBank 5.0: a major update to the DrugBank database for 2018. Nucleic Acids Res. 2018;46(D1):D1074–82.

    Article  CAS  PubMed  Google Scholar 

  29. Gene Ontology Consortium. Gene Ontology Consortium: going forward. Nucleic Acids Res. 2015;43(Database issue):D1049–56.

  30. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51.

    Article  CAS  PubMed  Google Scholar 

  33. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Wickham H. Elegant graphics for data analysis. Media. 2009;35(211):10.1007.

    Google Scholar 

  35. Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, et al. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45(D1):D362–8.

    Article  CAS  PubMed  Google Scholar 

  36. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.

    Article  CAS  PubMed  Google Scholar 

  37. Otasek D, Morris JH, Bouças J, Pico AR, Demchak B. Cytoscape Automation: empowering workflow-based network analysis. Genome Biol. 2019;20(1):185.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Jiang J, Huang K, Xu S, Garcia JGN, Wang C, Cai H. Targeting NOX4 alleviates sepsis-induced acute lung injury via attenuation of redox-sensitive activation of CaMKII/ERK1/2/MLCK and endothelial cell barrier dysfunction. Redox Biol. 2020;36:101638.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Cen M, Ouyang W, Zhang W, Yang L, Lin X, Dai M, et al. MitoQ protects against hyperpermeability of endothelium barrier in acute lung injury via a Nrf2-dependent mechanism. Redox Biol. 2021;41:101936.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Yu W, Li S, Wu H, Hu P, Chen L, Zeng C, et al. Endothelial Nox4 dysfunction aggravates atherosclerosis by inducing endoplasmic reticulum stress and soluble epoxide hydrolase. Free Radic Biol Med. 2021;164:44–57.

    Article  CAS  PubMed  Google Scholar 

  42. Wu RF, Ma Z, Liu Z, Terada LS. Nox4-derived H2O2 mediates endoplasmic reticulum signaling through local Ras activation. Mol Cell Biol. 2010;30(14):3553–68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Jiang F, Qian J, Chen S, Zhang W, Liu C. Ligustrazine improves atherosclerosis in rat via attenuation of oxidative stress. Pharm Biol. 2011;49(8):856–63.

    Article  CAS  PubMed  Google Scholar 

  44. Ko M, Oh GT, Park J, Kwon HJ. Extract of high hydrostatic pressure-treated danshen (Salvia miltiorrhiza) ameliorates atherosclerosis via autophagy induction. BMB Rep. 2020;53(12):652–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Song KH, Lee SH, Kim BY, Park AY, Kim JY. Extracts of Scutellaria baicalensis reduced body weight and blood triglyceride in db/db Mice. Phytother Res. 2013;27(2):244–50.

    Article  PubMed  Google Scholar 

  46. BaradaranRahimi V, Askari VR, Hosseinzadeh H. Promising influences of Scutellaria baicalensis and its two active constituents, baicalin, and baicalein, against metabolic syndrome: A review. Phytother Res. 2021;35(7):3558–74.

    Article  CAS  Google Scholar 

  47. Hou J, Wang J, Zhang P, Li D, Zhang C, Zhao H, et al. Baicalin attenuates proinflammatory cytokine production in oxygen-glucose deprived challenged rat microglial cells by inhibiting TLR4 signaling pathway. Int Immunopharmacol. 2012;14(4):749–57.

    Article  CAS  PubMed  Google Scholar 

  48. Yu WY, Gao CX, Zhang HH, Wu YG, Yu CH. Herbal active ingredients: potential for the prevention and treatment of acute lung injury. Biomed Res Int. 2021;2021:5543185.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Li L, Bao H, Wu J, Duan X, Liu B, Sun J, et al. Baicalin is anti-inflammatory in cigarette smoke-induced inflammatory models in vivo and in vitro: a possible role for HDAC2 activity. Int Immunopharmacol. 2012;13(1):15–22.

    Article  PubMed  Google Scholar 

  50. Sun JY, Li DL, Dong Y, Zhu CH, Liu J, Li JD, et al. Baicalin inhibits toll-like receptor 2/4 expression and downstream signaling in rat experimental periodontitis. Int Immunopharmacol. 2016;36:86–93.

    Article  CAS  PubMed  Google Scholar 

  51. Förstermann U, Xia N, Li H. Roles of vascular oxidative stress and nitric oxide in the pathogenesis of atherosclerosis. Circ Res. 2017;120(4):713–35.

    Article  PubMed  Google Scholar 

  52. Gray SP, Di Marco E, Kennedy K, Chew P, Okabe J, El-Osta A, et al. Reactive oxygen species can provide atheroprotection via NOX4-dependent inhibition of inflammation and vascular remodeling. Arterioscler Thromb Vasc Biol. 2016;36(2):295–307.

    Article  CAS  PubMed  Google Scholar 

  53. Craige SM, Kant S, Reif M, Chen K, Pei Y, Angoff R, et al. Endothelial NADPH oxidase 4 protects ApoE-/- mice from atherosclerotic lesions. Free Radic Biol Med. 2015;89:1–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Schröder K, Zhang M, Benkhoff S, Mieth A, Pliquett R, Kosowski J, et al. Nox4 is a protective reactive oxygen species generating vascular NADPH oxidase. Circ Res. 2012;110(9):1217–25.

    Article  PubMed  Google Scholar 

  55. Ochoa CD, Wu RF, Terada LS. ROS signaling and ER stress in cardiovascular disease. Mol Aspects Med. 2018;63:18–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Salazar G. NADPH oxidases and mitochondria in vascular senescence. Int J Mol Sci. 2018;19(5):1–21.

    Article  Google Scholar 

  57. Schürmann C, Rezende F, Kruse C, Yasar Y, Löwe O, Fork C, et al. The NADPH oxidase Nox4 has anti-atherosclerotic functions. Eur Heart J. 2015;36(48):3447–56.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Araujo JA, Barajas B, Kleinman M, Wang X, Bennett BJ, Gong KW, et al. Ambient particulate pollutants in the ultrafine range promote early atherosclerosis and systemic oxidative stress. Circ Res. 2008;102(5):589–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Thengchaisri N, Hein TW, Wang W, Xu X, Li Z, Fossum TW, et al. Upregulation of arginase by H2O2 impairs endothelium-dependent nitric oxide-mediated dilation of coronary arterioles. Arterioscler Thromb Vasc Biol. 2006;26(9):2035–42.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China. (Grant no.81973579).

Author information

Authors and Affiliations

  1. Taizhou Hospital, Shanghai University of Traditional Chinese Medicine, Taizhou, Zhejiang, China

    Mingshuang Li & Conglin Ren

  2. Taizhou Hospital of Traditional Chinese Medicine, Taizhou, Zhejiang, China

    Mingshuang Li & Conglin Ren

Authors
  1. Mingshuang Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Conglin Ren
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Mingshuang Li and Conglin Ren contributed to the experimental design for this study. Mingshuang Li contributed to carry out the experimental design and data analysis and wrote the paper. Conglin Ren performed the molecule docking analysis. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Conglin Ren.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no conflicts of interest in this work.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, M., Ren, C. Exploring the protective mechanism of baicalin in treatment of atherosclerosis using endothelial cells deregulation model and network pharmacology. BMC Complement Med Ther 22, 257 (2022). https://doi.org/10.1186/s12906-022-03738-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12906-022-03738-3

Keyword

BMC Complementary Medicine and Therapies

ISSN: 2662-7671

Contact us