Evidence-Based Complementary and Alternative Medicine

Evidence-Based Complementary and Alternative Medicine / 2019 / Article

Research Article | Open Access

Volume 2019 |Article ID 6503752 | https://doi.org/10.1155/2019/6503752

Yan Dong, Lian Duan, Heng-wen Chen, Yong-mei Liu, Yun Zhang, Jie Wang, "Network Pharmacology-Based Prediction and Verification of the Targets and Mechanism for Panax Notoginseng Saponins against Coronary Heart Disease", Evidence-Based Complementary and Alternative Medicine, vol. 2019, Article ID 6503752, 13 pages, 2019. https://doi.org/10.1155/2019/6503752

Network Pharmacology-Based Prediction and Verification of the Targets and Mechanism for Panax Notoginseng Saponins against Coronary Heart Disease

Academic Editor: Jae Youl Cho
Received27 Feb 2019
Revised08 May 2019
Accepted03 Jun 2019
Published03 Jul 2019


Coronary heart disease (CHD) is the worldwide leading cause for cardiovascular death. Panax notoginseng saponin (PNS), which is the main bioactive compound of panax notoginseng, has been generally accepted to exert a remarkable effect on CHD for a long time. However, to reveal the underlying treatment target and corresponding mechanism of PNS against CHD is still a substantial challenge. In this work, the targets and mechanism of PNS against CHD were successfully achieved by pharmacology-based prediction and experimental verification. 36 common targets were screened out through integrating the gene expression profile of CHD and the chemical-protein data of PNS. Then, two key nodes were further selected for verification by experiment after analyzing GO function, KEGG pathway, coexpression, and topology analysis. Results showed that PNS has protected the human umbilical vein endothelial cells from H2O2-induced oxidative stress by inhibiting early cell apoptosis via upregulating VEGFA mRNA expression. Therefore, our research has successfully pointed out one treatment target and apoptotic inhibition caused by PNS with method of integrating bioinformatics prediction and experimental verification, which has partially explained the pharmacological mechanism of PNS against CHD.

1. Introduction

Cardiovascular disease (CVD), as the leading reason in death worldwide, would increasingly result in over 23.6 million deaths by 2030 [1]. Coronary heart disease (CHD) is the most common cause of cardiovascular death, which accounts for about 44% [2]. 2018 AHA statistical data indicated that 16.5 million people over the age of 20 had CHD, and myocardial infarction (MI) occurred per 40 seconds approximately in America [2]. Moreover, silent MI or undiagnosed MI may account for almost 50% of the whole MIs [3]. This nonmanifestation MI is also an important risk factor for development of future heart failure [4]. Therefore, it could be speculated that the actual incidence of CHD is possibly higher than the reported data. Currently, the treatment strategy of CHD has achieved much progress involving new oral anticoagulant/antiplatelet drugs and drug eluting stent. However, recent treatment methods such as antiplatelet therapy, lipid regulation, and percutaneous coronary intervention (PCI) are still unsatisfactory for lots of patients worldwide due to difficulties in overcoming adverse effects. For example, as one cornerstone of CHD treatment, antiplatelet therapy encountered challenges of drug resistance [5, 6] and inappropriate use [7, 8] in clinical practice. Statins for lipid regulation sometimes induced the damage of liver and muscular cells [9]. Regarding PCI, there was increasing death rate observed from 2004 to 2014 [2]. As a result, it is noted that more effective and safer treatments for CHD are still urgently needed and continue to be challenging.

Panax notoginseng, as one of the traditional Chinese medicines (TCM) with effects of promoting blood circulation and stopping bleeding, has attracted enormous researches in the treatment of CVD. Panax notoginseng saponin (PNS) is the main bioactive compound in panax notoginseng, exerting cardiovascular protection. PNS mainly contains five compounds, namely, notoginsenoside R1, ginsenoside Rg1 (also called sanchinoside C1), ginsenoside Re, ginsenoside Rb1, and ginsenoside Rd (Figure 1), according to Pharmacopoeia of the People’s Republic of China [10] and Hong Kong Chinese Materia Medica Standards [11]. Currently, PNS has been made into several preparations for CHD therapy. A meta-analysis revealed that panax notoginseng preparations are beneficial for unstable angina pectoris patients [12], possibly owing to the vital importance of PNS in antiplatelet aggregation [13], anti-inflammation [14], antioxidation [15], antiatherosclerotic [16], antiapoptosis [17], and lipid regulation [16]. However, essential compounds and specific targets of PNS are not clearly determined and the underlying mechanism against CHD remains poorly defined. Therefore, to reveal and understand the treatment targets and mechanism of PNS against CHD is a subject which is worthy of study.

Since network pharmacology is an effective approach for drug target prediction, it has been increasingly growing in popularity [1820]. It could offer complementary support for new drug development and better understanding of complex mechanisms for drugs. Moreover, it can also provide systematic and holistic methods to deeply reveal the drug-gene-disease interactions which are especially suitable for understanding the multitarget, complementary, and synergic essence of TCM in molecular networks. Nevertheless, analysis of drug target based on network pharmacology is still in its infancy stage [21]. Predicted targets still need further confirmation from experiments. In this study, we predicted and verified the targets and mechanism for PNS against CHD through coupling the network pharmacology approach and experimental method. By integrating bioinformation from databases involving TCMSP, PharmMapper, STITCH, Gene Expression Omnibus (GEO) DataSets, WebGestalt, and GeneMANIA, key targets and corresponding mechanism of PNS against CHD were predicted. Then, cell experiment was applied to verify two key genes. As a result, one verified gene was identified as the treatment target of PNS against CHD, and the potential mechanism was regulation of apoptosis. The whole flow chart of our network prediction and experimental verification was shown in Figure 2.

2. Materials and Methods

2.1. Evaluation of ADME-Related Properties for PNS

PNS mainly contains notoginsenoside R1, ginsenoside Rg1, ginsenoside Re, ginsenoside Rb1, and ginsenoside Rd. The last compound is often exhibited as the metabolite, ginsenoside Rd qt. In order to assess druggable characters of the above five compounds, ADME (abbreviation of absorption, distribution, metabolism, and excretion) properties are evaluated by using TCMSP server (http://lsp.nwu.edu.cn/tcmsp.php). Particularly, important parameters of oral bioavailability (OB), drug-likeness (DL), molecular weight (MW), blood-brain barrier (BBB), intestinal epithelial permeability (Caco-2), fractional negative surface area (FASA-), the polar surface area (TPSA), and RBN [22] are used. Among all the descriptive parameters, OB and DL are the most two important indicators [23, 24] to assess and identify candidate compounds or drugs. OB is closely related to drug absorption and delivery of an orally administered dose. Drugs with high OB value mean they are more likely to reach the potential targets and can be therapeutic agents. DL can be used to estimate the druggability of a drug qualitatively, which helps to optimize pharmacokinetic and pharmaceutical properties. Commonly, OB30% and DL0.18 are suggested as drug screening criteria in the traditional Chinese herbs [22]. In the TCMSP server, OB value is calculated by an in-house model OBioavail1.1 [22]. Besides, the Tanimoto coefficient [25] is used to assess DL indexes of each compound in PNS. The formula is shown as follows [26]:

2.2. Targets Prediction of PNS against CHD

Pharmmapper (http://lilab-ecust.cn/pharmmapper/) and STITCH (http://stitch.embl.de/) were used to identify targets related to PNS in the study. Pharmmapper is a reverse docking server which identifies potential protein targets for small molecule compounds via a pharmacophore mapping approach. By calculating the fit scores and z’-scores (normalized fit scores) of a compound with different targets respectively, the correlation degree between the query compound and its targets are predicted. Generally, large positive z’-score indicates high significance [27], while STITCH provides a collection of protein-chemical interactions that integrates many sources of experimental and manually curated lines of evidence with text-mining information and interaction predictions [28, 29]. It calculates a score for each pair of protein-chemical interaction. Usually, interactions with score 0.4 are considered as medium confidence and those with the score0.7 are high confidence [30].

In addition, gene expression series of CHD were downloaded from GEO DataSets [31] (http://www.ncbi.nlm.nih.gov/geo/). Four researches with GSE20686, GSE71226, GSE18612, and GSE120774 accession numbers were offered in the platform of GPL4133, GPL570, and GPL6244, respectively. The expression profiles were all received by array. After unifying different gene names as Gene Symbol by UniProtKB (http://www.uniprot.org), the differential expression genes (DEGs) of CHD were obtained. Furthermore, Venn diagram (http://bioinformatics.psb.ugent.be/webtools/Venn/) was drawn to identify the intersection nodes of PNS and CHD. Finally, 36 common targets were predicted.

2.3. GO Function, KEGG Pathway, and Key Targets Analysis

GO function and KEGG pathway of 36 common targets were analyzed with WebGestalt (http://www.webgestalt.org/). Then, coexpression network of common targets was established by GeneMANIA (http://pages.genemania.org/) which was usually applied to generate hypotheses about gene function, analyze gene lists, and prioritize genes for functional assays [32]. In addition, topology analysis was conducted by CentiScape and MCODE to calculate the degree centrality and extract subnetwork, in order to identify key nodes in the coexpression network. Ultimately, two key targets which participate in apoptotic signaling pathway were selected to perform further cell experiment.

2.4. Experimental Materials and Cell Culture

Human umbilical vein endothelial cells (HUVECs) were purchased from Bei Na Chuanglian Biotechnology Research Institute, Beijing, China. PNS was purchased from National Institutes for Food and Drug Control, Beijing, China. PNS was dissolved with Phosphate Buffered Saline (1X) (Hyclone, Lot AC12557265) before experiment. Dulbecco’s modified Eagle’s medium (DMEM, Hyclone, Lot AD19822332) supplemented with 10% fetal bovine serum (Gibco, Lot 1739463) and 1% penicillin-streptomycin (Hyclone, Lot J160016) was used to culture HUVECs. The cultures were maintained in a humidified atmosphere containing 5% CO2 at 37°C. Then cells were inoculated onto 6-well plates to perform experiments and each group had 3 samples. Treatment groups were pretreated with 30 μg/mL PNS for 21 h and then added into 100 μmol·L−1 H2O2 for 3 h. In contrast, model groups were only treated with 100 μmol·L−1 H2O2 in the last 3 h. Cells in control groups were cultured in medium added with equal volume of DMEM.

2.5. RNA Isolation and Real-Time PCR Analysis

Total RNAs were extracted by using TRIzol reagent (Ambion, Lot 124006). Reverse transcript PCR was performed by FastQuant RT Kit (TianGen Biotech co., ltd, N2828, Beijing, China). Real-time quantitative PCR was conducted by PowerUP SYBR Green PCR Master Mix (Applied Biosystems, Lot 1710043) on 7900HT fluorescence quantitative PCR machine (Applied Biosystems). The primers of the two target genes and reference gene were as follows: vascular endothelial growth factor A (VEGFA) (forward: 5′-ACGAACGTACTTGCAGATGTG-3′; reverse: 5′-TTCTGTCGATGGTGATGGTGT-3′); B-cell lymphoma 2-related protein A1 (BCL2A1) (forward: 5′-GCACAATCACACACCTATGCT-3′; reverse: 5′-TGTGTTGGCAATCGTTTCCAT-3′); β-actin (forward: 5′-GAGACCTTCAACACCCCAGCC-3′; reverse: 5′-AATGTCACGCACGATTTCCC-3′). The relative mRNA expressions of VEGFA and BCL2A1 were calculated by   (ΔCt = Ct value of VEGFA or BCL2A1 minus Ct value of β-actin; ΔΔCt = ΔCt of the treatment group minus ΔCt of the treatment group) and fold change (FC). In order to ensure the quantitative accuracy, the experiments were performed in triplicate.

2.6. Cell Apoptosis Analysis

HUVECs were seeded in 6-well plates (1×106 cells/well) and incubated with 30 μg/mL PNS and/or 100 μmol·L−1 H2O2. After the culture for 24 h, cells were harvested and washed with cold PBS twice. Then cells were resuspended in 1×Binding Buffer and stained with FITC Annexin V and Propidium Iodide (PI) supplied in the FITC Annexin V Apoptosis Detection kit I (BD Pharmingen, USA). Next, the cells were incubated for 15 min at RT (25°C) in the dark. Finally, 1×Binding Buffer was added to cells and those cells were analyzed by a flow cytometer (BD FACSCalibur, USA) within 1 h. All experiments were performed in triplicate.

2.7. Statistical Analysis

The quantitative data were presented as mean±standard deviation (±SD). One-way analysis of variance (ANOVA) was applied to compare multiple means. When the variance was equal, the least-significant difference (LSD) method was used. Dunnett’s T3 test was conducted when the variance was not equal. Significance was defined as P<0.05. All statistical analysis was performed with SPSS 17.0 and data were presented as two decimals. In addition, all statistical column charts were drawn with Origin 8.0.

3. Results

3.1. ADME-Related Properties of PNS

All compounds of PNS showed superior DL value, except for ginsenoside Rb1. Particularly, ginsenoside Rg1 and ginsenoside Rd qt had the DL value above 0.18, which indicated that they presented excellent druggability such as better solubility and chemical stability [26]. Similarly, their OB indices which were estimated to be 10.04% and 12.23%, respectively, demonstrated that they can be efficiently absorbed to systemic circulation, compared with other three compounds. As for MW, TPSA, and RBN parameters, ginsenoside Rd qt showed better properties than others as well (see Table 1). Therefore, it could be seen that ginsenoside Rd qt was more druggable and may be the main bioactive compound in oral management.


MOL007487Notoginsenoside R15.430.13933.27-4.15-3.010.22298.1412
MOL005341Ginsenoside Rg110.040.28801.14-3.50-2.270.24239.2210
MOL005338Ginsenoside Re4.270.12947.3-4.39-3.20.25298.1412
MOL007476Ginsenoside Rb16.290.041,109.46-4.95-3.720.23377.2916
MOL007480Ginsenoside Rd qt12.230.77460.82-0.010.650.2460.694

3.2. Identification of Potential Targets for PNS

Potential targets were predicted by combination of Pharmmapper and STITCH servers. Firstly, the Tripos mol2 files of 5 compounds from TCMSP server were put into Pharmmapper, respectively. According to the z’-score, the top 300 potential targets of each compound were obtained (see attachment 1). Then, chemical names of notoginsenoside R1, ginsenoside Rg1, ginsenoside Re, ginsenoside Rb1, and ginsenoside Rd were put into STITCH to match their potential targets, with the minimum required interaction score being 0.4 (medium confidence). Correspondingly, 24 different targets involving NOS3, NR3C1, PRKAA1, STK11, MMP9, RB1, VEGFA, CHRNA2, CHRNA4, CHRNA7, RHAG, CHRFAM7A, RBL2, BDNF, AKT1, CASP3, CREB1, MAPK11, MAPK14, HBEGF, NFE2L2, CREM, CYCS, and TRPM7 were thereby predicted. Furthermore, after unifying different gene names as Gene Symbol by UniProtKB, 679 different targets were identified totally, among which 46 common targets were shared by all the compounds (Figure 3).

3.3. Common Targets of PNS and CHD

CHD-related genes were searched in GEO database and the series (GSE20686, GSE71226, GSE18612, and GSE120774) were selected. The former two contained DEGs of 399 samples totally from peripheral blood, while the latter two involved DEGs of 49 samples totally from epicardial and subcutaneous adipose tissue which were closely related to coronary atherosclerosis [3335]. CHD is diagnosed when the luminal stenosis50% in at least one major vessel. Therefore, patients who met the above diagnostic criteria were assigned to the CHD group, and the rest belonged to non-CHD group, when analyzing GSE20686 and GSE71226 data. Regarding GSE18612 and GSE120774, patients who underwent cardiac surgery for coronary artery disease were assigned to the CHD group; and those undergoing cardiac surgery for other reasons were considered as non-CHD group. By comparing CHD and non-CHD groups with GEO2R, there were 1775 DEGs finally identified in the above four series with selection criteria of FC1.5 and P<0.05. Moreover, further comparison between 1775 CHD-related genes and 679 PNS targets identified 36 common genes that may become the predicted treatment targets of PNS against CHD (Figure 4).

3.4. GO Function and KEGG Pathway Analysis of 36 Common Targets

To investigate the 36 common genes of PNS and CHD, GO function and KEGG pathway were analyzed with WebGestalt. GO enrichment was analyzed from the perspective of biological process, cellular component, and molecular function. It was noted that the common genes were mainly responsible for protein binding, biological regulation, response to stimulus, and metabolic process (Figure 5). Importantly, two apoptotic signaling pathways have been enriched in the top 10 GO functions (as shown in Table 2).


GO:0002446neutrophil mediated immunity70.0000740.0247786.5312
GO:0036230granulocyte activation70.0000770.0247786.4789
GO:0030099myeloid cell differentiation60.0001400.0247787.3264
GO:2001233regulation of apoptotic signaling pathway60.0001520.0247787.2122
GO:0006979response to oxidative stress60.0002800.0276156.4424
GO:0048872homeostasis of number of cells50.0001660.0247789.4445
GO:0072593reactive oxygen species metabolic process50.0002040.0247789.0387
GO:0090130tissue migration50.0003250.0276158.1763
GO:0097193intrinsic apoptotic signaling pathway50.0003350.0276158.119
GO:0042737drug catabolic process40.0002000.02477813.611

Additionally, KEGG pathways of the common genes were analyzed and shown in Figure 6. Among all the pathways, hsa04926, hsa05219, hsa04915, hsa05418, hsa05205, hsa04370, hsa04211, hsa04657, hsa01522, and hsa04933 enriched the most genes and ranked the top 10 levels (Table 3). Particularly, relaxin signaling pathway, fluid shear stress and atherosclerosis, VEGF signaling pathway, and IL-17 signaling pathway were reportedly related to CHD [3638].


hsa04926Relaxin signaling pathway60.0000010.00018218.143
hsa05219Bladder cancer40.0000030.00046538.352
hsa04915Estrogen signaling pathway40.0003000.02703911.477
hsa05418Fluid shear stress and atherosclerosis40.0004000.02703911.312
hsa05205Proteoglycans in cancer40.0014000.0611657.823
hsa04370VEGF signaling pathway30.0004000.02703919.988
hsa04211Longevity regulating pathway30.0014000.06116513.251
hsa04657IL-17 signaling pathway30.0016000.06116512.681
hsa01522Endocrine resistance30.0018000.06116512.034
hsa04933AGE-RAGE signaling pathway in diabetic complications30.0019000.06116511.912

3.5. Coexpression and Topology Analysis of 36 Common Targets

To further evaluate coexpression characteristics of 36 common targets, GeneMANIA was used to analyze their interactions among each other and with other related genes. The results showed that 34 common targets involving PDK2, DLG1, ACHE, TASP1, TIMM9, ARHGAP5, MMP9, USP14, FKBP8, RHAG, SOD2, VEGFA, HBEGF, PRDX6, CREB1, CAT, YARS, SENP7, BCL2A1, FCGR2A, MED7, RBP7, MSX1, S100A12, HMGB2, STRBP, RAB31, REPS2, SULT1B1, CNTN2, MMP1, TPSB2, MAPK14, and NOS3 were recognized in GeneMANIA. Through complex interactions involving coexpression, physical interactions, pathways, and colocalization, those 34 nodes established a coexpression network with other 20 nodes in GeneMANIA (Figure 7).

In addition, topology analysis was performed to identify key nodes in the coexpression network, involving calculation of degree centrality and k-core decomposition by CentiScape and MCODE, respectively, in Cytoscape 3.6.1. Degree centrality refers to the degree to which a node in a network is associated with all other nodes. Node with higher degree usually means more interaction edges and more significance in the network. The k-core of a network is defined as the largest subnetwork where every node has at least k links. Commonly, the extracted subnetwork is considered to be more stable and often used for hub nodes screening [39]. In this study, we set the parameters of degree>15 and k-core=2 to help identify key nodes in the coexpression network. Results indicated that common genes involving BCL2A1, MMP9, MAPK14, S100A12, and VEGFA have degree above 15 (Figure 8(a)). Furthermore, three subnetworks were extracted from the coexpression network by k-core decomposition (Figure 8(b)). Interestingly, it is shown in our analysis that VEGFA also constituted one part of a subnetwork, showing more possibility as a key node in the whole network. More importantly, VEGFA participates in regulation of apoptosis and BCL2A1 is a well-known antiapoptotic gene [40, 41]. As a result, we selected the two key nodes of VEGFA and BCL2A1 as potential targets for next experimental verification. Additionally, cell apoptosis was also observed to reveal VEGFA and BCL2A1-involved pathway in the treatment.

3.6. PNS Upregulated VEGFA mRNA Expression in H2O2-Treated HUVECs

To verify the predicted targets of VEGFA and BCL2A1 in PNS against CHD, 100 μmol·L−1 H2O2-induced HUVECs injury was performed as the oxidative stress model in cardiovascular diseases [42, 43]. In our experiment, the mRNA expressions of VEGFA and BCL2A1 were detected by real-time PCR in control, model, and treatment groups, respectively. According to the analyzing method of real-time PCR data [44], we performed the calculation of values and analyzed corresponding FC and expression level. The detail values were described in Table 4. Furthermore, statistical differences of relative expression of VEGFA and BCL2A1 in each group were also analyzed, respectively. Results showed that H2O2 significantly reduced VEGFA expression compared with the control group (P=0.041); moreover, PNS promoted VEGFA markedly compared to model group (P=0.007) (Figure 9(a)). With respect to BCL2A1, the mRNA expression of BCL2A1 had the trend to be downregulated by H2O2 and upregulated by PNS, but current results showed no statistical significance (P>0.05) (Figure 9(b)).

GeneGroup comparison valueFC valueP valueExpression

VEGFAModel vs Control0.62-1.610.041Downregulation
Treatment vs Model2.282.280.007Upregulation

BCL2A1Model vs Control0.81-1.23>0.05Downregulation
Treatment vs Model2.402.40>0.05Upregulation

3.7. PNS Inhibited Early Apoptosis in H2O2-Treated HUVECs

Previous reports have pointed out that VEGFA can participate in suppression of cell apoptosis [45, 46]. Thus, we further carried out apoptosis analysis to find out whether VEGFA overexpression induced by PNS could inhibit cell apoptosis in HUVECs. Our results showed that HUVECs in the model group were obviously reduced and disorderly arranged, while treatment group cells grew well compared with the model group (Figure 10). Apoptosis analysis further indicated that H2O2 induced cell apoptosis, but PNS pretreatment could protect cell from apoptosis (Figure 11). Particularly, the early apoptosis rate in model group was significantly higher than that in control group (P=0.004, Figure 12(a)). In contrast, PNS pretreatment markedly relieved H2O2-induced damages and protected HUVECs from early apoptosis when compared with model group (P=0.009, Figure 12(a)). However, no significant difference existed in the late apoptosis, although declined rates were also observed in treatment group (Figure 12(b)). These results suggested that PNS inhibited HUVECs apoptosis mainly ascribing to the suppression of early phase.

4. Discussion

In this study, we revealed that Ginsenoside Rd qt had better pharmacokinetics properties than the other four compounds in PNS. Thus, Ginsenoside Rd qt may be the main bioactive compound that exerted protective effect in the treatment of CHD. Notably, all compounds had relatively low OB value compared with the recommended screening criteria by TCMSP server [22], which indicated that PNS was not suitable for oral administration directly. More preparations involving injection [12, 47, 48] and spray [49] were therefore widely applied in practice or studied for new drug development. Although PNS has shown positive druggability and efficacy in clinical application, current researches about its specific targets and mechanism against CHD are still insufficient. Targets prediction based on network pharmacology have been reported [50, 51], but those predicted results are not convincing enough. Experimental verification is further needed to confirm the authenticity and reliability. As a result, we combined network pharmacology and cell experiment to predict and verify target and underlying mechanism of PNS in the treatment of CHD in our study. The results showed that VEGFA and BCL2A1 had great potential as target genes. However, only VEGFA was proved to be able to act as the target and participate in inhibition of apoptosis by cell experiment. Previous studies reported consistent results, indicating that VEGFA was an important protective factor in apoptosis of cardiomyocytes [45] and microvascular endothelial cells [46]. Moreover, VEGFA also took part in lncRNA NEAT1/miR-377/VEGFA axis to gain the protective effect in cell apoptosis [52]. It shows that VEGFA overexpression induced by PNS which inhibits cell apoptosis is reasonable. Notably, the inhibitory effect focused on the early apoptosis, rather than late phase. Therefore, our study revealed that PNS could upregulate the mRNA expression of VEGFA to participate in the suppression of apoptosis in early phase.

VEGFA was identified as a treatment target of PNS against CHD by integrating bioinformatics prediction and experimental verification. But beyond that, some other nodes in the coexpression network may also play important roles in the treatment of CHD, especially for MMP9 and MAPK14. A meta-analysis indicated that MMP9-1562C/T polymorphism is markedly associated with CHD [53]. Thus, the circulating level of MMP9 was upregulated in coronary artery disease, which may act as a risk factor for evaluating the development and severity of disease [54]. A mechanistic study further revealed that MMP-9 activation could promote cell migration to stimulate angiogenesis and preserve cardiac function [55]. As for MAPK14, also called p38MAPK , it played negative role in expression of vascular smooth muscle cell (VSMC) contractile genes and differentiation of VSMC [56]. Besides, it promoted platelet activation and served as a potential indicator of highly thrombotic lesions and no-reflow after myocardial infarction [57]. As a result, experimental verification for more potential targets needs to be performed in the future. In addition, because of the fact that the expression level of mRNA and protein may be inconsistent after treatment [58], the detection of the protein expression from VEGFA should be carried out by Western blot assay in subsequent research. Furthermore, there have been some reports [46, 52] about the modulation of VEGFA by miRNA and lncRNA. Therefore, to reveal the regulatory mechanism of PNS against CHD based on lncRNA-miRNA-VEGFA axis in the apoptosis pathway is also a subject which is worthy of exploration.

5. Conclusion

In conclusion, our study predicted and verified that PNS can promote mRNA expression of VEGFA for the suppression of early apoptosis in vascular endothelial cells to exert curative effect for CHD. The results helped partially explain the pharmacological mechanism of PNS against CHD. Moreover, our research provided an integrated method for discovery of drug targets and corresponding mechanism.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare no conflicts of interest.


This study was supported by the National Natural Science Foundation of China (No. 81673847; No. 81473561).

Supplementary Materials

The supplementary materials contain 5 files which are the top 300 potential targets of each compound predicted from Pharmmapper Server. (Supplementary Materials)


  1. S. Mendis, S. Davis, and B. Norrving, “Organizational update: the World Health Organization global status report on noncommunicable diseases 2014; one more landmark step in the combat against stroke and vascular disease,” Stroke, vol. 46, no. 5, pp. e121–e122, 2015. View at: Publisher Site | Google Scholar
  2. E. J. Benjamin, S. S. Virani, C. W. Callaway et al., “Heart disease and stroke statistics-2018 update: a report from the american heart association,” Circulation, vol. 137, no. 12, pp. e67–e492, 2018. View at: Google Scholar
  3. T. K. Paul and D. Mukherjee, “Silent myocardial infarction and risk of heart failure,” Annals of Translational Medicine, vol. 6, Suppl 1, p. S35, 2018. View at: Google Scholar
  4. W. T. Qureshi, Z. Zhang, P. P. Chang et al., “Silent myocardial infarction and long-term risk of heart failure,” Journal of the American College of Cardiology, vol. 71, no. 1, pp. 1–8, 2018. View at: Publisher Site | Google Scholar
  5. H. Zhang, X. Chen, L. Liu et al., “High prevalence of aspirin resistance in elderly patients with cardiovascular disease (CVD) and hyperhomocysteinaemia,” Archives of Gerontology and Geriatrics, vol. 59, no. 2, pp. 491–495, 2014. View at: Publisher Site | Google Scholar
  6. J. Siller-Matula, K. Schrör, J. Wojta, and K. Huber, “Thienopyridines in cardiovascular disease: Focus on Clopidogrel resistance,” Thrombosis and Haemostasis, vol. 97, no. 03, pp. 385–393, 2017. View at: Publisher Site | Google Scholar
  7. R. S. Hira, K. Kennedy, V. Nambi et al., “Frequency and practice-level variation in inappropriate aspirin use for the primary prevention of cardiovascular disease,” Journal of the American College of Cardiology, vol. 65, no. 2, pp. 111–121, 2015. View at: Publisher Site | Google Scholar
  8. R. S. Hira, K. Kennedy, H. Jneid et al., “Frequency and practice-level variation in inappropriate and nonrecommended prasugrel prescribing: Insights from the NCDR PINNACLE registry,” Journal of the American College of Cardiology, vol. 63, no. 25, Pt A, pp. 2876-2877, 2014. View at: Publisher Site | Google Scholar
  9. P. D. Thompson, G. Panza, A. Zaleski, and B. Taylor, “Statin-associated side effects,” Journal of the American College of Cardiology, vol. 67, no. 20, pp. 2395–2410, 2016. View at: Publisher Site | Google Scholar
  10. National Pharmacopoeia Commission, Pharmacopoeia of the People's Republic of China (2015 edition A), China Medical Science and Technology Press, Beijing, China, 2015.
  11. HKCMMS, Hong Kong Chinese Materia Medica Standards (HKCMMS), 2017, http://www.cmd.gov.hk/html/b5/hkcmms/cmmlist.html.
  12. H. Song, P. Wang, J. Liu, and C. Wang, “Panax notoginseng preparations for unstable angina pectoris: A systematic review and meta-analysis,” Phytotherapy Research, vol. 31, no. 8, pp. 1162–1172, 2017. View at: Publisher Site | Google Scholar
  13. Q. Shen, J. Li, C. Zhang et al., “Panax notoginseng saponins reduce high-risk factors for thrombosis through peroxisome proliferator-activated receptor-gamma pathway,” Biomedicine & Pharmacotherapy, vol. 96, pp. 1163–1169, 2017. View at: Publisher Site | Google Scholar
  14. J. Wei, X. Wen, P. W. Bible, Z. Li, R. B. Nussenblatt, and L. Wei, “Panax notoginseng saponin controls IL-17 Expression in helper T cells,” Journal of Ocular Pharmacology and Therapeutics, vol. 33, no. 4, pp. 285–289, 2017. View at: Publisher Site | Google Scholar
  15. J. Fan, D. Liu, C. He, X. Li, and F. He, “Inhibiting adhesion events by Panax notoginseng saponins and Ginsenoside Rb1 protecting arteries via activation of Nrf2 and suppression of p38 – VCAM-1 signal pathway,” Journal of Ethnopharmacology, vol. 192, pp. 423–430, 2016. View at: Publisher Site | Google Scholar
  16. J.-S. Fan, D.-N. Liu, G. Huang et al., “Panax notoginseng saponins attenuate atherosclerosis via reciprocal regulation of lipid metabolism and inflammation by inducing liver X receptor alpha expression,” Journal of Ethnopharmacology, vol. 142, no. 3, pp. 732–738, 2012. View at: Publisher Site | Google Scholar
  17. Z. Zhou, J. Wang, Y. Song et al., “Panax notoginseng saponins attenuate cardiomyocyte apoptosis through mitochondrial pathway in natural aging rats,” Phytotherapy Research, vol. 32, no. 2, pp. 243–250, 2018. View at: Publisher Site | Google Scholar
  18. B. Boezio, K. Audouze, P. Ducrot, and O. Taboureau, “Network-based approaches in pharmacology,” Molecular Informatics, vol. 36, no. 10, 2017. View at: Google Scholar
  19. S. Chen and M. C. Cui, “Systematic understanding of the mechanism of salvianolic acid a via computational target fishing,” Molecules, vol. 22, no. 4, p. E644, 2017. View at: Publisher Site | Google Scholar
  20. M. Huo, Z. Wang, D. Wu, Y. Zhang, and Y. Qiao, “Using coexpression protein interaction network analysis to identify mechanisms of danshensu affecting patients with coronary heart disease,” International Journal of Molecular Sciences, vol. 18, no. 6, p. E1298, 2017. View at: Publisher Site | Google Scholar
  21. J. Zhao, P. Jiang, and W. Zhang, “Molecular networks for the study of TCM pharmacology,” Briefings in Bioinformatics, vol. 11, no. 4, pp. 417–430, 2009. View at: Publisher Site | Google Scholar
  22. J. Ru, P. Li, J. Wang et al., “TCMSP: a database of systems pharmacology for drug discovery from herbal medicines,” Journal of Cheminformatics, vol. 6, no. 1, article no. 13, 2014. View at: Publisher Site | Google Scholar
  23. X. Xu, W. Zhang, C. Huang et al., “A novel chemometric method for the prediction of human oral bioavailability,” International Journal of Molecular Sciences, vol. 13, no. 6, pp. 6964–6982, 2012. View at: Publisher Site | Google Scholar
  24. W. Tao, X. Xu, X. Wang et al., “Network pharmacology-based prediction of the active ingredients and potential targets of Chinese herbal Radix Curcumae formula for application to cardiovascular disease,” Journal of Ethnopharmacology, vol. 145, no. 1, pp. 1–10, 2013. View at: Publisher Site | Google Scholar
  25. Y. Yamanishi, M. Kotera, M. Kanehisa, and S. Goto, “Drug-target interaction prediction from chemical, genomic and pharmacological data in an integrated framework,” Bioinformatics, vol. 26, no. 12, Article ID btq176, pp. i246–i254, 2010. View at: Publisher Site | Google Scholar
  26. T. Pei, C. Zheng, C. Huang et al., “Systematic understanding the mechanisms of vitiligo pathogenesis and its treatment by Qubaibabuqi formula,” Journal of Ethnopharmacology, vol. 190, pp. 272–287, 2016. View at: Publisher Site | Google Scholar
  27. X. Liu, S. Ouyang, B. Yu et al., “PharmMapper server: a web server for potential drug target identification using pharmacophore mapping approach,” Nucleic Acids Research, vol. 38, no. 2, pp. W609–W614, 2010. View at: Publisher Site | Google Scholar
  28. C. Von Mering, L. J. Jensen, B. Snel et al., “STRING: known and predicted protein-protein associations, integrated and transferred across organisms,” Nucleic Acids Research, vol. 33, supplement 1, pp. D433–D437, 2005. View at: Publisher Site | Google Scholar
  29. M. Kuhn, D. Szklarczyk, S. Pletscher-Frankild et al., “STITCH 4: integration of protein-chemical interactions with user data,” Nucleic Acids Research, vol. 42, no. 1, pp. D401–D407, 2014. View at: Publisher Site | Google Scholar
  30. D. Szklarczyk, A. Santos, C. Von Mering, L. J. Jensen, P. Bork, and M. Kuhn, “STITCH 5: Augmenting protein-chemical interaction networks with tissue and affinity data,” Nucleic Acids Research, vol. 44, no. 1, pp. D380–D384, 2016. View at: Publisher Site | Google Scholar
  31. T. Barrett, D. B. Troup, S. E. Wilhite et al., “NCBI GEO: mining tens of millions of expression profiles—database and tools update,” Nucleic Acids Research, vol. 35, supplement 1, pp. D760–D765, 2007. View at: Publisher Site | Google Scholar
  32. D. W. Farley, S. L. Donaldson, O. Comes et al., “The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function,” Nucleic Acids Research, vol. 38, no. 2, pp. W214–W220, 2010. View at: Publisher Site | Google Scholar
  33. M. Shimabukuro, Y. Hirata, M. Tabata et al., “Epicardial adipose tissue volume and adipocytokine imbalance are strongly linked to human coronary atherosclerosis,” Arteriosclerosis, Thrombosis, and Vascular Biology, vol. 33, no. 5, pp. 1077–1084, 2013. View at: Publisher Site | Google Scholar
  34. A. Salgado-Somoza, E. Teijeira-Fernández, Á. L. Fernández, J. R. González-Juanatey, and S. Eiras, “Changes in lipid transport-involved proteins of epicardial adipose tissue associated with coronary artery disease,” Atherosclerosis, vol. 224, no. 2, pp. 492–499, 2012. View at: Publisher Site | Google Scholar
  35. O. V. Gruzdeva, O. E. Akbasheva, Y. A. Dyleva et al., “Adipokine and cytokine profiles of epicardial and subcutaneous adipose tissue in patients with coronary heart disease,” Bulletin of Experimental Biology and Medicine, vol. 163, no. 5, pp. 608–611, 2017. View at: Publisher Site | Google Scholar
  36. J. M. Raleigh, S. Toldo, A. Das, A. Abbate, and F. N. Salloum, “Relaxin' the heart: a novel therapeutic modality,” Journal of Cardiovascular Pharmacology and Therapeutics, vol. 21, no. 4, pp. 353–362, 2015. View at: Publisher Site | Google Scholar
  37. Z. Podemska-Jedrzejczak, A. Malinska, P. Sujka-Kordowska et al., “Vascular restenosis in coronary artery bypass grafting might be associated with VEGF-C/VEGFR-3 signaling pathway,” Heart and Vessels, vol. 33, no. 9, pp. 1106–1120, 2018. View at: Publisher Site | Google Scholar
  38. S. P. Barry, S. Ounzain, J. McCormick et al., “Enhanced IL-17 signalling following myocardial ischaemia/reperfusion injury,” International Journal of Cardiology, vol. 163, no. 3, pp. 326–334, 2013. View at: Publisher Site | Google Scholar
  39. Y. Zhang, Z. Li, M. Yang et al., “Identification of GRB2 and GAB1 coexpression as an unfavorable prognostic factor for hepatocellular carcinoma by a combination of expression profile and network analysis,” Plos One, vol. 8, no. 12, Article ID e85170, 2013. View at: Publisher Site | Google Scholar
  40. H. Zhang, Y. Wang, and Z. He, “Glycine-Histidine-Lysine (GHK) alleviates neuronal apoptosis due to intracerebral hemorrhage via the mir-339-5p/vegfa pathway,” Frontiers in Neuroscience, vol. 12, p. 644, 2018. View at: Google Scholar
  41. J. Kim, H. J. Cho, M. Cho et al., “Prenylated Rab acceptor RABAC1 inhibits anti-apoptotic protein BCL2A1 and induces apoptosis,” Biochemical and Biophysical Research Communications, vol. 513, no. 4, pp. 940–946, 2019. View at: Publisher Site | Google Scholar
  42. L. Ni, T. Li, B. Liu et al., “The protective effect of Bcl-xl overexpression against oxidative stress-induced vascular endothelial cell injury and the role of the Akt/eNOS pathway,” International Journal of Molecular Sciences, vol. 14, no. 11, pp. 22149–22162, 2013. View at: Publisher Site | Google Scholar
  43. S. Jiao, B. Tang, Y. Wang et al., “Pro-angiogenic role of danqi pill through activating fatty acids oxidation pathway against coronary artery disease,” Frontiers in Pharmacology, vol. 9, article no. 1414, 2018. View at: Publisher Site | Google Scholar
  44. T. D. Schmittgen and K. J. Livak, “Analyzing real-time PCR data by the comparative CT method,” Nature Protocols, vol. 3, no. 6, pp. 1101–1108, 2008. View at: Publisher Site | Google Scholar
  45. B. Guo, J. Xu, M. Xiao, M. Ding, and L. Duan, “Puerarin reduces ischemia/reperfusion-induced myocardial injury in diabetic rats via upregulation of vascular endothelial growth factor A/angiotensin-1 and suppression of apoptosis,” Molecular Medicine Reports, vol. 17, no. 5, pp. 7421–7427, 2018. View at: Publisher Site | Google Scholar
  46. W. Liu, “MicroRNA‐9 inhibits retinal neovascularization in rats with diabetic retinopathy by targeting vascular endothelial growth factor A,” Journal of Cellular Biochemistry, vol. 120, no. 5, pp. 8032–8043, 2018. View at: Publisher Site | Google Scholar
  47. X. Zhao, F. Zhang, and Y. Wang, “Proteomic analysis reveals Xuesaitong injection attenuates myocardial ischemia/reperfusion injury by elevating pyruvate dehydrogenase-mediated aerobic metabolism,” Molecular BioSystems, vol. 13, no. 8, pp. 1504–1511, 2017. View at: Publisher Site | Google Scholar
  48. J. Zhang, Q. He, T. Huang, and B. Zhang, “Effects of panax notoginseng saponins on homing of C-kit+ bone mesenchymal stem cells to the infarction heart in rats,” Journal of Traditional Chinese Medicine, vol. 31, no. 3, pp. 203–208, 2011. View at: Publisher Site | Google Scholar
  49. H.-M. Wang, T.-M. Fu, and L.-W. Guo, “Preparation of Panax notoginseng saponins-tanshinone IIA composite method for pulmonary delivery with spray-drying mothod and its characterization,” Zhongguo Zhongyao Zazhi, vol. 38, no. 4, pp. 559–563, 2013. View at: Google Scholar
  50. B. Zhu, W. Zhang, Y. Lu et al., “Network pharmacology-based identification of protective mechanism of Panax Notoginseng Saponins on aspirin induced gastrointestinal injury,” Biomedicine & Pharmacotherapy, vol. 105, pp. 159–166, 2018. View at: Publisher Site | Google Scholar
  51. W. Xie, X. Meng, Y. Zhai et al., “Panax notoginseng saponins: A review of its mechanisms of antidepressant or anxiolytic effects and network analysis on phytochemistry and Pharmacology,” Molecules, vol. 23, no. 4, article no. E940, 2018. View at: Google Scholar
  52. Z. Zhou, L. Zheng, X. Ren, A. Li, and W. Zhou, “LncRNA NEAT1 facilitates survival and angiogenesis in oxygen-glucose deprivation (OGD)-induced brain microvascular endothelial cells (BMECs) via targeting miR-377 and upregulating SIRT1, VEGFA, and BCL-XL,” Brain Research, vol. 1707, pp. 90–98, 2019. View at: Publisher Site | Google Scholar
  53. J. Li, H. Lu, F. Tao et al., “Meta-analysis of MMP9-562C/T and the risk of coronary heart disease,” Cardiology, vol. 124, no. 1, pp. 53–59, 2013. View at: Publisher Site | Google Scholar
  54. K. Mahmoodi, K. Kamali, E. Karami, and M. S. Soltanpour, “Plasma concentration, genetic variation, and gene expression levels of matrix metalloproteinase 9 in Iranian patients with coronary artery disease,” Journal of Research in Medical Sciences, vol. 22, no. 1, p. 8, 2017. View at: Google Scholar
  55. R. Chen, X. Cai, J. Liu, B. Bai, and X. Li, “Sphingosine 1-phosphate promotes mesenchymal stem cell-mediated cardioprotection against myocardial infarction via ERK1/2-MMP-9 and Akt signaling axis,” Life Sciences, vol. 215, pp. 31–42, 2018. View at: Publisher Site | Google Scholar
  56. X. Long, S. L. Cowan, and J. M. Miano, “Mitogen-activated protein kinase 14 is a novel negative regulatory switch for the vascular smooth muscle cell contractile gene program,” Arteriosclerosis, Thrombosis, and Vascular Biology, vol. 33, no. 2, pp. 378–386, 2013. View at: Publisher Site | Google Scholar
  57. P. Shi, L. Zhang, M. Zhang et al., “Platelet-specific p38αdeficiency improved cardiac function after myocardial infarction in mice,” Arteriosclerosis, Thrombosis, and Vascular Biology, vol. 37, no. 12, pp. e185–e196, 2017. View at: Publisher Site | Google Scholar
  58. M. Huang, Y. Chou, M. Li et al., “Epigenetic targeting DNMT1 of pancreatic ductal adenocarcinoma using interstitial control release biodegrading polymer reduced tumor growth through hedgehog pathway inhibition,” Pharmacological Research, vol. 139, pp. 50–61, 2019. View at: Publisher Site | Google Scholar

Copyright © 2019 Yan Dong et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.