Abstract Objective Baolier Capsule (BLEC) is a Traditional Mongolian Medicine comprising fifteen herbs. This study aims to illustrate the synergistic mechanism of BLEC in the treatment of Coronary Artery Disease (CAD) by using network pharmacology method, molecular docking and experimental validation. Methods Searching and screening the active ingredients of different herbs in BLEC and target genes related to CAD in multiple databases. Subsequently, Protein–Protein Interactions Network (PPI-Net), gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment were used to identify the key targets. AutoDock was used to verify the binding ability between the active ingredient and key target through molecular docking. Reverse Transcription-Quantitative Real-Time Polymerase Chain Reaction (RT-qPCR) was used to verify the effect of active ingredient of BLEC on the key target gene. Finally, effect of BLEC on the degree of blood lipids and atherosclerosis was validated by animal experiment. Results There are 144 active components and 80 CAD-related targets that are identified in BLEC in the treatment of CAD. What is more, 8 core genes were obtained by clustering and topological analysis of PPI-Net. Further, GO and KEGG analysis showed that fluid shear stress and atherosclerosis are the key pathways for BLEC to treat CAD. These results were validated by molecular docking method. In vitro, active compounds of BLEC (Quercetin, luteolin, kaempferol, naringenin, tanshinone IIA, β-carotene, 7-O-methylisomucronulatol, piperine, isorhamnetin and Xyloidone) can inhibit 8 core gene (AKT1, EGFR, FOS, MAPK1, MAPK14, STAT3, TP53 and VEGFA) expression. Moreover, BLEC not only improve blood lipid levels but also inhibit the development of atherosclerosis in ApoE-knockout mice. Conclusion Our research first revealed the basic pharmacological effects and related mechanisms of in the treatment of CAD. The predicted results provide some theoretical support for BLEC or its important active ingredients to treat CAD. Keywords: coronary artery disease, Baolier capsule, traditional Mongolian medicine, network pharmacology, experimental validation Introduction At present, major progress has been made to prevent and treat cardiovascular disease (CVD) through lifestyle changes and drugs, but CVD is still the main cause of death in humans, especially in developed country and developing country.[32]1 It is worth noting that the residual risk of CVD in patients is still high, even though on the basis of the best lifestyle intervention and high-dose statin therapy, new lipid-lowering drugs such as ezetimibe, proprotein convertase subtilisin/kexin type 9 (PCSK9) inhibitors and cholesterol ester transfer protein (CETP) inhibitors, and even kanazumab are used.[33]2–6 Therefore, looking for drugs that can prevent and treat coronary artery disease (CAD) is of great significance to further decrease the incidence of major adverse cardiovascular events (MACEs) in patients who suffer from CAD. Traditional Chinese medicine has opened up a new path for the treatment of CAD. Among them, Mongolian medicine, as one of the outstanding cultural heritage of the Mongolian nationality, is a summary of the Mongolian experience in fighting against diseases for a long time, and gradually formed by absorbing the experience of traditional Chinese medicine and Tibetan medicine. The curative effect is remarkable and it is increasingly accepted by people. Yogurt, which is a popular Mongolian dairy product, contributed to prevent and manage hyperlipidemia, blood glucose and inflammatory marker.[34]7 Syringa pinnatifolia Hemsl., which is a Mongolian folk medicine in China, may protect against myocardial ischemia.[35]8 Eerdun Wurile (EW) is a well established traditional Mongolian medicine with proven clinical application for treatment of stroke recovery.[36]9,[37]10 Qingganjiuwei powder (QGJWS) is a commonly used Mongolian drug to treat patients with chronic hepatic disease through protecting liver fibrosis.[38]11 The use of Mongolian drugs to treat CAD has certain application prospect. Baolier capsule (BLEC), one of the Mongolian drugs, comes from the classic Mongolian medicine book “Coral Experience Prescription”, “The Set of Mongolian medicine prescription” and “Zhigao prescription”.[39]12 BLEC consists of fifteen herbs.[40]12 BLEC can protect vascular endothelial function, inhibit inflammation and significantly improve the symptoms of patients with CAD, including number and duration of episodes of angina and number of using for Nitroglycerin.[41]12 However, BLEC is formed of dozens of compounds from Chinese patent medicine and Chinese herbal compound formula, which leads to the difficulties in exploring BLEC mechanism. It is important clinically to discover the target genes and compounds of BLEC in discovering new drugs. Our study innovatively applies network pharmacology method, molecular docking and cellar or animal experiment to analyze key targets, and signal pathways of BLEC in treating CAD, to guide future researches of this formula. The detailed workflow is shown in [42]Figure 1. Figure 1. [43]Figure 1 [44]Open in a new tab Schematic of the network pharmacological analysis used to identify the potential mechanism of BLEC action on CAD. Abbreviations: BLEC, Baolier capsule; CAD, Coronary artery disease. Materials and Methods Network Pharmacology Prediction Collection and Screening of Candidate Compounds in Baolier Capsule The herbs in BLEC decoction include Guangzao (Choerospondiatis Fructus), Danshen (Radix Salviae), Roudoukou (Myristicae Semen), Zhizi (Gardeniae Fructus), Qiancao (Rubia Cordifolia), Honghua (Carthami Flos), Sanqi (Panax Notoginseng (Burk). F. H. Chen Ex C. Chow), Muxiang (Aucklandiae Radix), Tanxiang (Santalum Album L)., Niuhuang (Bovis Calculus), Dahuang (Radix Rhei Et Rhizome), Mutong (Caulis Akebiae), Huangqi (Hedysarum Multijugum Maxim). Biba (Piperis Longi Fructus) and Hezi (Chebulae Fructus). The active ingredients in BLEC are gained from the Traditional Chinese Medicine Systems Pharmacology (TCMSP) Database and Analysis Platform ([45]http://tcmspw.com/tcmsp.php), with the keywords Guangzao, Danshen, Roudoukou, Zhizi, Qiancao, Honghua, Sanqi, Muxiang, Tanxiang, Niuhuang, Dahuang, Mutong, Huangqi, Biba and Hezi respectively. TCMSP is a Chinese herbal medicine system pharmacology platform that includes 449 herbs, can catch the relations between diseases, targets and drugs. The database comprises targets, chemicals and corresponding drug-target networks. Oral Bioavailability (OB) is an important parameter in pharmacokinetics, which is directly related to bioavailability. Drug-like (DL) indicates active ingredient’s similarities to present drugs. We selected active ingredients with OB > 30% and DL > 0.18 for later research.[46]13 Identification of Potential Targets of Active Ingredients and CAD We obtain potential targets of active ingredients in BLEC from the database, retrieve related targets of CAD from Gene cards ([47]https://www.genecards.org), Online Mendelian Inheritance in Man (OMIM, [48]https://www.omim.org/), PharmGKB ([49]http://www.pharmgkb.org), Therapeutic Target Database (TTD, [50]http://db.idrblab.net/ttd/) and DrugBank ([51]https://go.drugbank.com/). In addition, genes related to CAD, which are obtained from Gene cards, were selected according to relevance score greater or equal to 30. UniProt ([52]https://www.uniprot.org) can help people find full name, ID and symbols of related gene. Construction of Ingredient-Target-Disease Network To further explore the relations between potential targets and active ingredients, we loaded them into Cytoscape 3.8.0 to form an ingredient-target-disease network. Construction of Protein–Protein Interactions Network (PPI-Net) STRING ([53]http://string-db.org) is a convenient online database to build PPI-Net.[54]14 Every interaction between target proteins will be scored. Then, following steps will be performed: selecting “multiple proteins”, inputting target genes that may interact with active ingredients in BLEC and related genes of CAD, searching in “Homo sapiens”, setting minimum required interaction score to high confidence (0.9) and hiding disconnected nodes in the network. Finally, target proteins will be put into Cytoscape to build the interaction network. Analysis of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) Target genes which are potential ones of active ingredients in BLEC and related ones in CAD are put into the R project for GO and KEGG analysis, in order to discover what role the selected targets play in signal pathways. Molecular Docking Construction of Receptor Structure The ID of core target proteins of PPI-Net was obtained from UniProt. According to the ID, receptor structure file was downloaded from RCSB PDB ([55]http://www1.rcsb.org/) as pdb format. We import the receptor structure file into pyMOL software and save it in pdb format after removing molecules and small molecule ligands. Construction of Ligand Structure The 2D structure of the ligand was downloaded from PubChem ([56]https://pubchem.ncbi.nlm.nih.gov/) as SDF format. The ligand structure file is imported into Chem3D software and converted into 3D structure then saved as mol2 format file after using the MM2 position to optimize the energy of all small molecules. Molecular Docking Calculation method of molecular docking Autodock vina ADT auto dock tool MGL tools Use ADT 1.5.6 (Auto Dock Tool) to convert the format of receptor and ligand structure files. Import The receptor structure pdb format file was imported into ADT 1.5.6 (Auto Dock Tool) and was saved as pdbqt format after adding hydrogens. The ligand structure pdb format file was converted to pdbqt format. The receptor structure was converted to 2D structure through ADT 1.5.6. set. The number of points in X-dimension, Y-dimension and Z-dimension are set as 40, respectively. Spacing (angstrom) was set as 1. Active pocket parameters is exported to GPF format file. The input file of ligand and receptor structure is modified according to active pocket parameters. The maximum energy difference is 5. The number of models is 20. The outputting the pdbqt file has been imported the pyMOL software together with the receptor file to save the optimal model. Experimental Verification Cell Culture The murine RAW264.7 macrophage cell line was purchased from the China Center for Type Culture Collection (CCTCC, China). The RAW264.7 macrophages were seeded in 60-mm cell culture dishes at 1.0×10^6 cells per dish in DMEM containing 10% fetal bovine serum and maintained at 37°C in a humidified atmosphere of 5% CO[2]. The cells were synchronized by replacing the medium with DMEM supplemented with 0.5% BSA for 12 h before the experiments. After 12h, the cells were transfected with PBS or active ingredient of drugs (Final concentration about 25 μmol/L) for 48 h. Quercetin, luteolin, kaempferol, naringenin, tanshinone IIA, β-carotene, 7-O-methylisomucronulatol, piperine, isorhamnetin and Xyloidone were purchased commercially (Solarbio, Beijing, China). Reverse Transcription-Quantitative Real-Time Polymerase Chain Reaction (RT-qPCR) Total mRNA was extracted from adenovirus-transfected RAW264.7 cells treated with PBS or active ingredient of drugs using TRIzol reagent (Invitrogen, Gaithersburg, MD) according to the manufacturer’s instructions. The mRNA was reverse transcribed using a PrimeScript™ RT Reagent Kit (TaKaRa, China). Real-time quantitative PCR (qPCR) was performed on a LightCycler 480 II (Roche, Switzerland) using a SYBR Premix Ex Taq Kit (TaKaRa, China). The mRNA expression level of GAPDH was used as an internal control. The quantitative measurements were determined using the ^ΔΔCt method.[57]15 All samples were assayed in triplicate, and the mean value was considered in the comparative analyses. All qPCR primer sequences were obtained from PrimerBank and are listed in [58]Supplementary Table 1. Animals Model ApoE-knockout (ApoE^−/−) mice are a classic model of atherosclerosis. All animal experiments were conducted in strict accordance with Southern Medical University guidelines for the Care and Use of Experimental Animals. ApoE^−/− mice were approved by the Animal Experiment Committee of Nanfang Hospital at Southern Medical University to be used in this study. The research ethical code for the study is 2020-BSGZ-07-01. The mice were fed with a high-fat diet (HFD, 15% fat and 1.2% cholesterol) and intragastrically administered BLEC once per day for 12 weeks to establish an AS model. According to previous studies and dose of clinical use about BLEC,[59]12 ApoE^−/− mice were randomly divided into 4 groups: Control group, Low-Dose BLEC group (37.5mg/kg), Middle-Dose BLEC group (75mg/kg), High-Dose BLEC group (150mg/kg). At the end of the experiment, blood samples and aorta of the mice were collected for further studies. ApoE^−/− mice on a C57BL/6 background were obtained from Vital River Laboratory Animal Technology Co. Ltd. (Beijing, China) and HFD was purchased from Guangdong Provincial Animal Medicine Experimental Center. During the experiment, all mice were weighed twice a week, and the drug doses were adjusted accordingly. The Detection of Lipids Serum was collected to measure the level of lipids, including triglyceride (TG), total cholesterol (TC), low-density lipoprotein cholesterol (LDL-C) and high-density lipoprotein cholesterol (HDL-C) using assay kits. Whole blood was collected into ethylenediaminetetraacetic acid anticoagulant tubes and stored at 4°C. Samples were centrifuged for 15 minutes at 1000×g at 4°C within 30 minutes of collection and then stored at −80°C until analysis. The quantification assay kits were purchased from Nanjing Jiancheng Biotech Co., Ltd (Nanjing, Jiangsu, China). Hematoxylin and Eosin (H&E) Staining When the mice were sacrificed, the tissue of the aortic root was collected, fixed with formalin and made into paraffin sections (5μm) and frozen sections (5μm). The above sections were subjected to H&E. Images are best obtained within 7 days of staining using a Japanese Olympus BX63 microscope. Statistical Analysis All statistical analyses were performed using SPSS 25.0 software. The measurement data are expressed as mean±SD. One-way analysis of variance (ANOVA) and Fisher’s least significant difference (LSD) were used. P value <0.05 was considered statistically significant. Results Active Ingredients in BLEC One thousand three hundred and thirty-one bioactive ingredients are selected from TCMSP database consisting of 43 types in Guangzao, 202 types in Danshen, 64 types in Roudoukou, 98 types in Zhizi, 54 types in Qiancao, 189 types in Honghua, 119 types in Sanqi, 106 types in Muxiang, 70 types in Tanxiang, 19 types in Niuhuang, 92 types in Dahuang, 43 types in Mutong, 87 types in Huangqi, 104 types in Biba and 41 types in Hezi. According to the standards OB and DL, 228 active ingredients are selected for further study, including 9 types in Guangzao, 65 types in Danshen, 9 types in Roudoukou, 15 types in Zhizi, 19 types in Qiancao, 22 types in Honghua, 8 types in Sanqi, 6 types in Muxiang, 3 types in Tanxiang, 5 types in Niuhuang, 16 types in Dahuang, 8 types in Mutong, 20 types in Huangqi, 15 types in Biba and 8 types in Hezi ([60]Supplementary Table 2). Potential Target Genes of Active Ingredients in BLEC After removal of duplicates, 280 target genes of BLEC are acquired from TCMSP database and stored for further study ([61]Supplementary Table 3). Related Genes of CAD As [62]Figure 2 shows, 618 target genes of CAD are attained from Gene cards, OMIM, PharmGKB, TTD and DrugBank. As [63]Figure 3 shows, there are a total of 80 target genes, which are potential ones of BLEC’s active ingredients and CAD’s related genes. Figure 2. [64]Figure 2 [65]Open in a new tab Identification of the CAD-related genes by taking a union of all the results from 5 database. Abbreviation: CAD, Coronary artery disease. Figure 3. [66]Figure 3 [67]Open in a new tab Identification of the drug-target disease-related genes by taking an intersection of drug target genes and CAD-related genes. Abbreviation: CAD, Coronary artery disease. Construction of Ingredient-Target-Disease Interactive Network To thoroughly reveal the potential ingredients and targets of BLEC against CAD, the “component-target-disease” interactive network is constructed. Moreover, this research applies network pharmacological analysis to conduct visual analysis with relevant parameters of each target. One hundred and ninety-two nodes and 543 edges with 144 ingredients and 80 targets were displayed in this network. The network shows that the nodes possessing more edges and greater degree value play a more important role in the regulation ([68]Figure 4). Accordingly, this research also establishes the network of ingredients with common targets of drugs and the disease, to directly display the interactions between potential targets and active ingredients. Degree refers to the number of routes that connect to one node and other in the network. The higher the value is, the more significant the corresponding ingredient or target is. After the Network Analyzer plugin in Cytoscape software done its analysis, ten essential active ingredients were listed in [69]Tables 1. The result tells us that active ingredients of BLEC mainly achieved its effect on CAD through acting on the essential targets above. Figure 4. [70]Figure 4 [71]Open in a new tab Construction of the drug-target pharmacology network. The drug-targets interaction pharmacology network. Circles represent the small molecule active compounds in BLEC. Each colour represents a traditional Chinese medicine ingredient. Rectangle represents the CAD-related target genes, and edges represent the interaction between the small molecule compounds and the target genes. Abbreviation: CAD, Coronary artery disease. Table 1. Active Ingredients of BLEC in the Treatment of CAD Molecular ID Ingredient Degree Source MOL000098 Quercetin 51 Zhizi/Sanqi/Guangzao/Huangqi/Honghua MOL000006 Luteolin 18 Tanxiang/Danshen/Honghua MOL000422 Kaempferol 18 Zhizi/Guangzao/Huangqi/Honghua MOL004328 Naringenin 16 Guangzao MOL007154 Tanshinone IIA 13 Danshen MOL002773 β-carotene 13 Honghua MOL000378 7-O-methylisomucronulatol 13 Huangqi MOL001592 Piperine 11 Biba MOL000354 Isorhamnetin 10 Tanxiang/Huangqi MOL006174 Xyloidone 10 Qiancao [72]Open in a new tab Outcomes of PPI-Net We use STRING to construct PPI-Net to study the mechanism of target genes of BLEC in treating CAD in accordance with protein-level high confidence (0.9). There are 541 edges and 73 nodes in this network ([73]Figure 5). In order to obtain core target proteins, PPI network was filtered using Cytoscape software. If the Betweenness (BC), Closeness (CC), Degree (DC), Eigenvector (EC), LAC and Network (NC) of the protein is greater than the median, the protein is retained. After two filtration (Filter 1: BC = 21.71263724, CC = 0.1518718315, DC = 4.5, EC = 0.0574660525, LAC = 1.6, NC = 2.7083333335. Filter 2: BC = 8.3583916085, CC = 0.605263158, DC = 8, EC = 0.189531028, LAC = 3.25, NC = 4.3714285715), a total of 8 core target proteins were retained ([74]Figure 6). Figure 5. [75]Figure 5 [76]Open in a new tab Protein–Protein interaction (PPI) network. Figure 6. [77]Figure 6 [78]Open in a new tab Identification of key subnetwork using Cytoscape. PPI network and the first filtration by CytoNca, the yellow nodes were screened with each score higher than median. Subnetwork constructed by a second filtration via CytoNca. The yellow nodes were screened with a score higher than the median. Final key subnetwork screened after two filtrations using CytoNca. Key subnetwork of top 9 nodes analysed by CytoHubba. Results of GO and KEGG Enrichment Analysis This research uses DAVID database for GO analysis of BLEC with potential targets for CAD. Two thousand one hundred and ninety-five GO entries are selected based on P value (P < 0.05). It includes 128 cellular components, 73 molecular function, and 1994 biological processes. In these biological processes, potential targets mainly focus on responding nutrient levels, oxidative stress, oxygen levels and so on. And the targets mainly relate to membrane raft, membrane microdomain, membrane region and so on in the cellular components. For molecular function, potential targets mainly work for phosphatase binding, protein phosphatase binding, cytokine receptor binding and so on ([79]Figure 7 and [80]Supplementary Table 4). Figure 7. [81]Figure 7 [82]Open in a new tab The top 10 of GO enrichment analysis. When we use DAVID to do the KEGG enrichment analysis of BLEC for CAD (P < 0.05), 162 pathways are acquired. According to P value, the top 20 pathways are selected ([83]Figure 8A and [84]Supplementary Table 5). The pathway map of the fluid shear stress and atherosclerosis is illustrated in [85]Figure 8B. Figure 8. [86]Figure 8 [87]Open in a new tab KEGG enrichment analysis and pathway map. (A) KEGG enrichment analysis of the target genes. Gene ratio refers to the ratio of enriched genes to all target genes. Counts refer to the number of the enriched genes. (B) Pathway map of fluid shear stress and atherosclerosis as the most significant enriched pathway. Molecular Docking of Active Compounds and Core Target Proteins This research use the 8 core target proteins from PPI-Net to perform molecular docking. Then, we get active compounds targeting core target proteins from the compound-target interaction network respectively. Subsequently, as it shows in molecular docking, all active compounds can enter and bind the active pocket of the protein easily ([88]Supplementary Figure 1) and the interactions with the biding site with the key amino acid ([89]Figure 9A–H). The scores of docking are displayed in [90]Table 2. Figure 9. [91]Figure 9 [92]Open in a new tab Screening of the key genes in the subnetwork and further molecular docking. (A) Molecular docking between the six small molecule ligands and protein 5WBL (encoded by AKT1). (B) Molecular docking between the two small molecule ligands and protein 2RGP (encoded by EGFR). (C) Molecular docking between the four small molecule ligands and protein 1FOS (encoded by FOS). (D) Molecular docking between the three small molecule ligands and protein 3W55 (encoded by MAPK1). (E) Molecular docking between the one small molecule ligands and protein 6NUQ (encoded by STAT3). (F) Molecular docking between the eight small molecule ligands and protein 3L8X (encoded by MAPK14). (G) Molecular docking between the five small molecule ligands and protein 2VPF (encoded by VEGFA). (H) Molecular docking between the six small molecule ligands and protein 7BWN (encoded by TP53). Abbreviations: AKT1, AKT Serine/Threonine Kinase 1; EGFR, Epidermal Growth Factor Receptor; FOS, Proto-oncogene c-Fos; MAPK1, Mitogen-activated protein kinase 1; MAPK14, Mitogen-activated protein kinase 14; STAT3, Signal transducer and activator of transcription 3; TP53, Cellular tumor antigen p53; VEGFA, Vascular endothelial growth factor A. Table 2. Molecular Docking Score Protein Name Molecule Name Docking Score (kcal/mol) AKT1 Baicalein −7.8 β-carotene −8.0 Kaempferol −7.7 Luteolin −8.0 Naringenin −8.0 Quercetin −8.1 EGFR Luteolin −6.0 Quercetin −5.6 FOS Baicalein −9.5 Piperine −8.6 Quercetin −10.2 Tanshinone IIA −9.4 MAPK1 Luteolin −9.0 Naringenin −8.7 Quercetin −8.7 MAPK14 3_Methylkempferol −7.7 7-O-methylisomucronulatol −7.5 Calycosin −7.4 Formononetin −7.5 Isoguaiacin −7.5 Isorhamnetin −7.5 Sudan III −7.4 Xyloidone −7.5 STAT3 Cryptotanshinone −7.4 TP53 Aloe-emodin −7.8 Baicalein −8.1 Ellipticine −8.8 Luteolin −8.6 Quercetin −8.3 Tanshinone IIA −9.0 VEGFA Baicalein −7.5 Β-carotene −7.8 Ellagic acid −7.5 Luteolin −7.5 Quercetin −7.3 [93]Open in a new tab Effects of Active Compounds of BLEC on Core Target Gene Expression in RAW264.7 Cells Cellular experiments were conducted to validate the effects of active compounds of BLEC on key genes. The mRNA expressions of the 8 core target gene in active compounds group were significantly reduced relative to control group ([94]Figure 10A–J). In conclusion, active compounds of BLEC inhibits target genes of BLEC in treating CAD expression in macrophages. Figure 10. [95]Figure 10 [96]Open in a new tab Effects of active compounds of BLEC on core target gene expression in RAW264.7 cells. (A) The relative mRNA expression level of akt1, egfr, fos, mapk1, tp53, vegfa in RAW264.7 cells with PBS or quercetin. (B) The relative mRNA expression level of akt1, egfr, mapk1, tp53, vegfa in RAW264.7 cells with PBS or luteolin. (C) The relative mRNA expression level of akt1 in RAW264.7 cells with PBS or kaempferol. (D) The relative mRNA expression level of akt1 and mapk1 in RAW264.7 cells with PBS or naringenin. (E) The relative mRNA expression level of fos and tp53 in RAW264.7 cells with PBS or tanshinone IIA. (F) The relative mRNA expression level of akt1 and vegfa in RAW264.7 cells with PBS or β-carotene. (G) The relative mRNA expression level of mapk14 in RAW264.7 cells with PBS or 7-O-methylisomucronulatol. (H) The relative mRNA expression level of fos in RAW264.7 cells with PBS or piperine. (I) The relative mRNA expression level of mapk14 in RAW264.7 cells with PBS or isorhamnetin. (J) The relative mRNA expression level of mapk14 in RAW264.7 cells with PBS or Xyloidone. N=6. Data are shown as means±SD. *p < 0.05 vs. The control group. BLEC Attenuates the Development of Atherosclerosis and Decreases TC, TG, LDL-C and HDL-C in vivo Cholesterol metabolic disorder is a key causative factor in atherosclerosis, and thus we further examined the effect of BLEC on serum lipid levels in ApoE^−/− mice fed an HFD. As shown in [97]Figure 11A, the mice treated with BLEC showed lower plasma levels of TC, TG, LDL-C and HDL-C than the control mice. The aorta of the above mice in four groups were stained with H&E staining to observe the occurrence of AS and lipid core area. It was suggested that compared with the control mice, mice treatment with BLEC had less lipid core area ([98]Figure 11B). These results show that BLEC inhibits the development of atherosclerosis by modulating plasma lipid levels. Figure 11. [99]Figure 11 [100]Open in a new tab BLEC attenuates the development of atherosclerosis and decreases TC, TG, LDL-C and HDL-C in vivo. (A) The serum lipid levels in ApoE^−/− mice treatment with LD, MD and HD BLEC and fed an HFD for 12 weeks was analyzed. (B) Hematoxylin and eosin (H&E)-stained cross-sections of aortic root tissue from ApoE^−/− mice treatment with LD, MD and HD BLEC and fed an HFD for 12 weeks (magnification, 10×). The circular region indicates the necrotic core area. Necrotic core areas were defined as anuclear, antifibrotic, and eosin-negative areas. Scale bar = 200 μm. All data are expressed as the mean ± SD. *P < 0.05 vs The control group. Abbreviations: TC, total cholesterol; TG, triglyceride; LDL-C: low-density lipoprotein cholesterol; HDL-C, high-density lipoprotein cholesterol. Discussion CAD is called “heart tingling” in Mongolian medicine. It is due to the dysfunction of the “three roots” and “seven elements”, which hinders the movement of Heyi and the blood. The three roots refer to Heyi Xila and Badagan. Seven elements refer to dietary essence, blood, muscle, fat, bone, bone marrow and semen. Combining modern Western medicine theories, Mongolian medicine currently divides CAD into Heyi heart tingling, bloody heart tingling, adhesive heart tingling and Huyang heart tingling. Mongolian medicine mainly treats CAD by promoting the differentiation of turbidity and improving Heyi and blood circulation.[101]16 BLEC is a prescription confirmed by famous Mongolian medicine experts repeatedly based on the above-mentioned etiology and pathogenesis and combined with Mongolian medicine theory and clinical practice experience.[102]12 It can treat CAD by adjusting functions of the body’s three root and seven element.[103]12 To the best of our knowledge, this is the first study integrating network pharmacology and molecular docking analyses to reveal the pharmacological mechanisms of BLEC for treating CAD. We identified a total of 228 active compounds in BLEC. Furthermore, a total of 80 potential target genes related to the action of BLEC on CAD were identified. PPI analysis revealed that the top-ranking genes, AKT Serine/Threonine Kinase 1 (AKT1), Epidermal Growth Factor Receptor (EGFR), Proto-oncogene c-Fos (FOS), Mitogen-activated protein kinase 1 (MAPK1), Mitogen-activated protein kinase 14 (MAPK14), Signal transducer and activator of transcription 3 (STAT3), Cellular tumor antigen p53 (TP53) and Vascular endothelial growth factor A (VEGFA), may be the crucial targets for BLEC treating CAD. Functional enrichment analysis revealed the over-represented GO terms and their functional domains. KEGG pathway enrichment analysis revealed that the 80 target proteins were significantly enriched in 162 related signal pathways. Among them, fluid shear stress and atherosclerosis appeared to be the most critical pathways involved in the treatment of CAD. Furthermore, molecular docking analysis demonstrated that the representative compounds could bind to the target protein binding site. In this study, systematic pharmacological method is applied to forecast and explain the molecular mechanisms of BLEC on CAD. In BLEC’s active components-targets network, we get 81 targets affected by 144 bioactive compounds. Several recent studies have found that the active compounds of plants play an important role in the treatment of diseases. Such as, polysaccharides, flavonoids and saponins of genus Polygonatum have function about controlling blood glucose and preventing obese diabetes.[104]17 Xanthohumol of Humulus lupulus L.has strong anti-cancer activity against various types of cancer cells.[105]18 A considerable phenolic content for both varieties in Phaseolus vulgari, and a strong antiobesity activity.[106]19 The 10 top-ranking compounds including quercetin, luteolin, kaempferol, naringin, tanshinone IIA, beta-carotene, 7-O-methylisomucronulatol, piperine, isorhamnetin, Xyloidone are identified as BLEC’s active ingredients. Their biological activities against CAD were mentioned above. Quercetin is a natural bioflavonoid found in vegetables and fruits including onions, grapes, tea, apples, and red wine. It can improve endothelial cell function,[107]20–22 inhibit dendritic cell activation,[108]23 inhibit foam cell formation by inhibiting oxidized low-density lipoprotein-induced apoptosis of macrophages[109]24 and regulates cholesterol metabolism and lipoprotein metabolism by regulating cholesterol metabolism-related gene expression.[110]25–27 Luteolin, is a natural flavone which is abundant in edible plants, including green chilies, broccoli, French beans, onion leaf, white radish, carrots, ragweed pollen and clover blossom.[111]28 Studies have shown that it can inhibit the occurrence of atherosclerosis by inhibiting macrophage inflammation.[112]29,[113]30 Kaempferol, is a common natural yellow flavonoid with a low molecular weight, which exists in a number of traditional Chinese medicine and plant-derived foods. Kaempferol and its glycosides have many pharmacological functions, for example, they can serve as anti-inflammatory, antioxidant, and anticancer, antidiabetic, neuroprotective, analgesic, antimicrobial, and anti-allergic drugs.[114]31 Studies have shown that kaempferol conducts anti-atherosclerotic effect by regulating the expression of gene and protein in inflammatory molecules.[115]32 Naringin is the bitter component in a kind of grapefruit which is often found in its flowering parts, peels, juice and vegetative parts.[116]33 Naringin can improve blood lipid levels,[117]34 metabolic imbalances,[118]35 and inhibit inflammatory reactions,[119]36 and inhibit atherosclerosis.[120]37 Tanshinone IIA is a main bioactive component isolated from the roots of the Chinese herb Salviae miltiorrhiza Bunge (Danshen),[121]38 can regulate the autophagy and polarization of macrophages,[122]39 and even the expression of non-coding RNA related to atherosclerosis and inflammation,[123]40 antioxidant,[124]41 improve endothelial cell function,[125]42 and regulate cholesterol metabolism-related pathways affect cholesterol metabolism.[126]43,[127]44 Epidemiological evidence suggests that beta-carotene helps decrease the risk for atherosclerosis by protecting LDL from oxidation.[128]45,[129]46 Isorhamnetin, a bioactive compound, exists in herbal plants, and possesses multiple biological properties. It can inhibit macrophage apoptosis and inhibit atherosclerosis.[130]47 According to PPI core network analysis and molecular docking results, it is shown that the main components of BLEC can prevent and treat CAD by affecting CAD-related top-ranking genes AKT1, EGFR, FOS, MAPK1, MAPK14, STAT3, TP53 and VEGFA. A large number of studies have confirmed that intervention of the above key genes can achieve the effect of preventing and treating atherosclerosis. The lack of AKT1 can promote the release of pro-inflammatory factors, and enhance the apoptosis of endothelial cells and macrophages, leading to the development of atherosclerosis.[131]48 The expression of FOS in monocytes in the blood of patients with CAD is significantly increased, and its expression level is related to the severity of atherosclerosis.[132]49 Inhibiting the expression of FOS of THP-1 cells can reduce the release of inflammatory factors.[133]49 The activation of STAT3 can promote the occurrence of atherosclerosis by affecting endothelial cell dysfunction, macrophage polarization, inflammation, and immunity. Inhibition of STAT3 can achieve the effect of inhibiting atherosclerosis.[134]50 P53 is a classic tumor suppressor gene. Inhibition of P53 can significantly increase the proportion of cell proliferation and apoptosis, and promote the occurrence of atherosclerosis.[135]51 The MAPK pathway is a classic inflammatory pathway. Inhibiting MAPK can reduce the release of macrophages MCP-1, inhibit the infiltration of macrophages into the plaque, and inhibit the occurrence and development of atherosclerosis.[136]52 The expression of EGFR and its ligands on arterial plaques of patients with CAD is significantly increased. By inhibiting EGFR, the uptake and inflammatory release of lipids by macrophages can be restricted, thereby inhibiting the occurrence and development of atherosclerosis.[137]53 VEGFA plays an important role in the formation of blood vessels, and can induce apolipoprotein to change in the direction of pro-inflammatory, and its expression increases in atherosclerosis model mice. Inhibiting the expression of VEGFA in vascular smooth muscle can reduce the size of atherosclerotic plaques.[138]54 Limitations However, the study also has its limitations. On the one hand, we ignored the actual dosage and performance of traditional Mongolian medicine, which can be significant in research. Besides, network pharmacology is based on the data in the existing literature. The prediction is limited. Conclusion In this pharmacology network–based study, active ingredients, core targets, and key signal pathways of BLEC in treating CAD are systematically investigated. In vitro, active compounds of BLEC (Quercetin, luteolin, kaempferol, naringenin, tanshinone IIA, β-carotene, 7-O-methylisomucronulatol, piperine, isorhamnetin and Xyloidone) can inhibit 8 core target gene (AKT1, EGFR, FOS, MAPK1, MAPK14, STAT3, TP53 and VEGFA) expression. While Animal experiments confirmed that BLEC is a potent prescription to improve the development of atherosclerosis and plasma lipid level. Acknowledgment