Research Article | Open Access
Investigating the Mechanisms of Action of Depside Salt from Salvia miltiorrhiza Using Bioinformatic Analysis
Salvia miltiorrhiza is a traditional Chinese medicinal herb used for treating cardiovascular diseases. Depside salt from S. miltiorrhiza (DSSM) contains the following active components: magnesium lithospermate B, lithospermic acid, and rosmarinic acid. This study aimed to reveal the mechanisms of action of DSSM. After searching for DSSM-associated genes in GeneCards, Search Tool for Interacting Chemicals, SuperTarget, PubChem, and Comparative Toxicogenomics Database, they were subjected to enrichment analysis using Multifaceted Analysis Tool for Human Transcriptome. A protein-protein interaction (PPI) network was visualised; module analysis was conducted using the Cytoscape software. Finally, a transcriptional regulatory network was constructed using the TRRUST database and Cytoscape. Seventy-three DSSM-associated genes were identified. JUN, TNF, NFKB1, and FOS were hub nodes in the PPI network. Modules 1 and 2 were identified from the PPI network, with pathway enrichment analysis, showing that the presence of NFKB1 and BCL2 in module 1 was indicative of a particular association with the NF-B signalling pathway. JUN, TNF, NFKB1, FOS, and BCL2 exhibited notable interactions among themselves in the PPI network. Several regulatory relationships (such as JUN TNF/FOS, FOS NFKB1 and NFKB1 BCL2/TNF) were also found in the regulatory network. Thus, DSSM exerts effects against cardiovascular diseases by targeting JUN, TNF, NFKB1, FOS, and BCL2.
Salvia miltiorrhiza (also named Danshen) is widely used as a traditional Chinese medicinal herb to prevent and treat vascular diseases [1, 2]. Depside salt from S. miltiorrhiza (DSSM) is a new medicine that contains the active components magnesium lithospermate B (MLB), lithospermic acid (LA), and rosmarinic acid (RA) . MLB plays protective roles in relieving atherosclerosis and combating myocardial ischaemia-reperfusion injury . RA and LA also have beneficial effects on cardiovascular diseases, such as atherosclerosis and neointimal hyperplasia [5, 6]. Thus, understanding the mechanisms of action of DSSM is important for its better utilisation in a clinical setting.
In patients with unstable angina, DSSM can suppress platelet activation and aggregation as well as matrix metallopeptidase 9 (MMP-9) expression and secretion . MLB protects against diabetic atherosclerosis by inducing the nuclear factor erythroid 2-related factor-2-antioxidant responsive element-NAD(P)H: quinone oxidoreductase-1 pathway . Du et al. reported that MLB can be used to treat ischaemic heart diseases as it specifically inhibits transforming growth factor β-activated protein kinase 1-binding protein 1-p38 apoptosis signalling . Kim et al. also demonstrated that RA inhibits adriamycin-induced cardiotoxicity by suppressing reactive oxygen species generation as well as extracellular signal-regulated kinase and c-Jun N-terminal kinase (JNK) activation . Moreover, LA inhibits foetal bovine serum-induced vascular smooth muscle cell (VSMC) proliferation by arresting cell cycle progression and suppressing cyclin D1 expression and lipopolysaccharide-induced VSMC migration by downregulating MMP-9 expression; thus, LA may be used for preventing neointimal hyperplasia, restenosis, and atherosclerosis . However, no comprehensive survey of DSSM-associated genes has been reported.
In this study, we searched for DSSM-associated genes in several common databases. Subsequently, we applied multiple bioinformatic methods to identify further DSSM-associated key genes; these methods included enrichment analysis, protein-protein interaction (PPI) network and module analyses, and transcriptional regulatory network analysis. This study may provide a better understanding of the mechanisms of action of DSSM.
2.1. Search for DSSM-Associated Genes
Several databases were used in this study. The human genomic database GeneCards (http://www.genecards.org/) was used , which provides concise genome, transcriptome, proteome, and function data of all predicted and known human genes. Search Tool for Interacting Chemicals (STICH, version 5.0, http://stitch.embl.de/) was also sued , which is a database that integrates the interactions between chemicals and proteins. SuperTarget (http://insilico.charite.de/supertarget/)  was also used, which is a database that includes drug-associated information correlated with drug metabolism, adverse drug effects, medical indications, Gene Ontology (GO) terms, and pathways for target proteins. The PubChem (https://www.ncbi.nlm.nih.gov/pccompound/)  database was also used, which provides information on chemical substances and corresponding biological activities and is linked to the National Institutes of Health PubMed Entrez. Moreover, the study features the Comparative Toxicogenomics Database (CTD, http://ctdbase.org/) , which is a database associated with environmental chemical-gene product interactions as well as the transportation and accumulation of chemical substances in the human body. In GeneCards , STICH (parameters set as organism = human, score > 0.4) , SuperTarget , PubChem , and CTD  databases, “magnesium lithospermate B”, “lithospermic acid” and “rosmarinic acid” were used as keywords to search for genes associated with them (which were combined into a set of DSSM-associated genes).
2.2. Functional and Pathway Enrichment Analyses
The GO (http://www.geneontology.org) database describes the associations of gene products with the categories of biological process (BP), molecular function (MF), and cellular component (CC), as well as with more specific subcategories . The Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.ad.jp/kegg) is a database used for annotating the functions of genes or other molecules . The “BioCloud” online platform was developed for managing problems encountered in the analysis of high-throughput data. By applying the Multifaceted Analysis Tool for Human Transcriptome (MATHT, http://www.biocloudservice.com) in the “BioCloud” online platform, GO functional and KEGG pathway enrichment analyses were conducted for DSSM-associated genes, setting a threshold of <0.05 for the false discovery rate (FDR).
2.3. PPI Network and Module Analyses
On combining the Search Tool for the Retrieval of Interacting Genes (STRING, version 10.0, http://www.string-db.org/, combined score > 0.4) , Biological General Repository for Interaction Datasets (BioGRID, version 3.4, https://wiki.thebiogrid.org/)  and Human Protein Reference Database (release 9, http://www.hprd.org/)  interaction databases, PPI pairs among DSSM-associated genes were predicted. Subsequently, the PPI network was visualised for DSSM-associated genes using the Cytoscape software (http://www.cytoscape.org) . Using the CytoNCA plug-in  (version 2.1.6, http://apps.cytoscape.org/apps/cytonca) in Cytoscape, degree centrality (DC), betweenness centrality (BC), and closeness centrality of the nodes were analysed to obtain the hub proteins in the PPI network . The parameter was set as “without weight.”
Based on the MCODE plug-in  (version 1.4.2; http://apps.cytoscape.org/apps/mcode; parameters set as degree cut-off = 2, maximum depth = 100, node score cut-off = 0.2, and -core = 2) in Cytoscape, module analysis was conducted for the PPI network. Subsequently, KEGG pathway enrichment analysis was performed for the nodes of significant modules, with FDR < 0.05 as the cut-off criterion.
2.4. Transcriptional Regulatory Network Construction
Using the transcriptional regulatory relationships unravelled by a sentence-based text-mining (TRRUST, http://www.grnpedia.org/trrust/)  database, transcription factors (TFs) among DSSM-associated genes were searched and then their targets were screened. Finally, a transcriptional regulatory network was constructed using Cytoscape .
3.1. Search for DSSM-Associated Genes
The numbers of MLB-, RA-, and LA-associated genes identified from GeneCards, STICH, SuperTarget, PubChem, and CTD databases are listed in Table 1. The MLB-, RA-, and LA-associated genes were combined into a group of 73 DSSM-associated genes.
|MLB, magnesium lithospermate B; RA, rosmarinic acid; LA, lithospermic acid, STICH, Search Tool for Interacting Chemicals; CTD, Comparative Toxicogenomics Database.|
3.2. Functional and Pathway Enrichment Analyses
The main BP, CC, and MF terms, as well as KEGG pathways enriched for DSSM-associated genes, are shown in Figure 1. The significantly enriched terms included inflammatory response (BP), cytosol (CC), protein binding (MF), and tumour necrosis factor (TNF) signalling pathway (KEGG pathway).
3.3. PPI Network and Module Analyses
The PPI network for DSSM-associated genes comprised 65 nodes and 431 edges (Figure 2). The 65 nodes included 21 MLB-associated genes, 56 RA-associated genes, and eight LA-associated genes. Among them, two genes [nitric oxide synthase 2 (NOS2) and nitric oxide synthase 3 (NOS3)] were associated with all MLB, RA, and LA; four genes (interleukin 2, interleukin 1 beta, caspase 3, and v-rel avian reticuloendotheliosis viral oncogene homolog A) were associated with both MLB and RA; and four genes (xanthine dehydrogenase, aldo-keto reductase family 1 member B1, protein tyrosine phosphatase no-receptor type 1 and procollagen-lysine, and 2-oxoglutarate 5-dioxygenase) were associated with both MLB and LA. Jun protooncogene (JUN), TNF, nuclear factor kappa B subunit 1 (NFKB1), and Fos protooncogene (FOS) were noted as particular hub nodes according to their BC, closeness centrality, and DC scores (Table 2). Moreover, two modules were also identified from the PPI network: module 1 (21 nodes, score = 13.2) and module 2 (12 nodes, score = 7.8) (Figure 3). In addition, pathway enrichment analysis demonstrated that the nuclear factor kappa B (NF-B) signalling pathway [Figure 4; = 3.14E−06; involving NFKB1 and B-cell CLL/lymphoma 2 (BCL2)] and leishmaniasis ( = 4.72E−07) were particularly associated with the nodes in modules 1 and 2, respectively (Table 3 and Figure 5). JUN, TNF, NFKB1, FOS, and BCL2 also exhibited notable interactions among themselves in the PPI network.
|PPI, protein-protein interaction; DC, degree centrality; BC, betweenness centrality; CC, closeness centrality.|
3.4. Transcriptional Regulatory Network Analysis
Based on the TRRUST database, TFs among DSSM-associated genes were searched for and their targets were screened. The constructed transcriptional regulatory network had 37 nodes (including eight TFs and 29 target genes) and 108 relationship pairs (such as JUN TNF/FOS, FOS NFKB1, and NFKB1 BCL2/TNF) (Figure 6).
In this study, a total of 73 DSSM-associated genes were identified from various databases. In the PPI network constructed on the basis of these genes and their interactions, JUN, TNF, NFKB1, and FOS were established as hub nodes according to their BC, closeness centrality, and DC scores. JUN, TNF, NFKB1, FOS, and BCL2 were also shown to interact among themselves in the PPI network. Moreover, two distinct modules (modules 1 and 2) of the PPI network were identified. In addition, several regulatory relationships (such as JUN TNF/FOS, FOS NFKB1, and NFKB1 BCL2/TNF) were found to be involved in the transcriptional regulatory network.
Prolonged anti-TNF-α therapy has beneficial effects on the signs of subclinical cardiovascular disease in patients with severe psoriasis . Studies have shown that, compared with nonbiologic disease-modifying antirheumatic drugs, TNF-α-blocking agents may contribute to the reduction of the risk of cardiovascular events in patients with rheumatoid arthritis [28, 29]. TNF-α plays a critical role in vascular dysfunction in cardiovascular diseases, which may be exploited to treat inflammation in a clinical setting . TNF-α boosts atherosclerosis development by promoting the transcytosis of low-density lipoprotein (LDL) across endothelial cells, thus contributing to the reservoir of LDL in vascular walls . These findings indicate that TNF may be associated with the functions of DSSM in cardiovascular diseases.
Blocking the JNK pathway may provide a new strategy for treating the cardiomyocyte death induced by myocardial ischaemia/reperfusion [32, 33]. Dominant negative c-Jun (DN-c-Jun) gene transfer inhibits VSMC proliferation, and JUN is associated with the intimal hyperplasia induced by balloon injury . Saliques et al. considered that FOS is a novel factor that determines the severity and development of atherosclerosis and is thus involved in tobacco toxicity in coronary artery disease (CAD) patients . Palomer et al. also demonstrated that miR-146a-targeting FOS is a potential tool for treating enhanced inflammation-associated cardiac disorders . Thus, JUN and FOS may also be targets of DSSM in cardiovascular diseases.
Studies involving Uygur and Han women in China have demonstrated that the DD genotype of NFKB1 polymorphism (rs28362491) may be a genetic marker of CAD [37, 38]. The NFKB1-94ins/del ATTG polymorphism may reduce the susceptibility to myocardial infarction by decreasing activated NF-B, which is in turn correlated with the reduction of plasma inflammatory markers . By increasing BCL2 expression, miR-21 promotes heart failure progression with preserved left ventricular ejection fraction and subsequently inhibits cardiac fibrosis . Moreover, Tang et al. reported that miR-1 targets BCL2 functions in mediating cardiomyocyte apoptosis . In this study, pathway enrichment analysis showing the presence of NFKB1 and BCL2 in module 1 indicated a particular association with the NF-B signalling pathway, indicating that NFKB1 and BCL2 may be correlated with the effects of DSSM against cardiovascular diseases through this pathway.
In conclusion, a total of 73 DSSM-associated genes were identified by conducting a search of public databases. JUN, TNF, NFKB1, FOS, and BCL2 were also revealed as potential targets of DSSM for treating cardiovascular diseases. However, further experimental research should be performed to confirm these findings obtained by bioinformatic analysis.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
This work was supported by the National Nature Science Foundation of China (no. 81500229).
- X. Li, C. Yu, Y. Lu, and et al., “Pharmacokinetics, tissue distribution, metabolism, and excretion of depside salts from Salvia miltiorrhiza in rats,” Drug Metabolism and Disposition the Biological Fate of Chemicals, vol. 35, no. 2, pp. 234–239, 2006.
- B.-Q. Wang, “Salvia miltiorrhiza chemical and pharmacological review of a medicinal plant,” Journal of Medicinal Plants Research, vol. 425, no. 25, pp. 2813–2820, 2010.
- C. X. Zhou, H. W. Luo, and M. Niwa, “Studies on isolation and identification of water-soluble constituents of Salvia miltiorrhiza,” Journal of China Pharmaceutical University, vol. 30, no. 6, pp. 411–416, 1999.
- X. Li, C. Yu, L. Wang et al., “Simultaneous determination of lithospermic acid B and its three metabolites by liquid chromatography/tandem mass spectrometry,” Journal of Pharmaceutical and Biomedical Analysis, vol. 43, no. 5, pp. 1864–1868, 2007.
- S.-S. Huang and R.-L. Zheng, “Rosmarinic acid inhibits angiogenesis and its mechanism of action in vitro,” Cancer Letters, vol. 239, no. 2, pp. 271–280, 2006.
- L. Wang, Q. Zhang, X. Li et al., “Pharmacokinetics and metabolism of lithospermic acid by LC/MS/MS in rats,” International Journal of Pharmaceutics, vol. 350, no. 1-2, pp. 240–246, 2008.
- MU. Xinyu, “Influence study of depside salt from Salvia Miltiorrhiza on blood platelet function and MMP-9 in patients with unstable angina,” China Medical Herald, vol. 6, no. 24, pp. 59-60, 2009.
- K. Y. Hur, S. H. Kim, M.-A. Choi et al., “Protective effects of magnesium lithospermate B against diabetic atherosclerosis via Nrf2-ARE-NQO1 transcriptional pathway,” Atherosclerosis, vol. 211, no. 1, pp. 69–76, 2010.
- C.-S. Du, R.-F. Yang, S.-W. Song et al., “Magnesium lithospermate B protects cardiomyocytes from ischemic injury via inhibition of TAB1-p38 apoptosis signaling,” Frontiers in Pharmacology, vol. AUG, Article ID Article 119, 2010.
- D.-S. Kim, H.-R. Kim, E.-R. Woo, S.-T. Hong, H.-J. Chae, and S.-W. Chae, “Inhibitory effects of rosmarinic acid on adriamycin-induced apoptosis in H9c2 cardiac muscle cells by inhibiting reactive oxygen species and the activations of c-Jun N-terminal kinase and extracellular signal-regulated kinase,” Biochemical Pharmacology, vol. 70, no. 7, pp. 1066–1078, 2005.
- L. Chen, W.-Y. Wang, and Y.-P. Wang, “Inhibitory effects of lithospermic acid on proliferation and migration of rat vascular smooth muscle cells,” Acta Pharmacologica Sinica, vol. 30, no. 9, pp. 1245–1252, 2009.
- M. Safran, I. Dalah, J. Alexander et al., “GeneCards Version 3: the human gene integrator,” Database, vol. 2010, Article ID baq020, 2010.
- 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.
- N. Hecker, J. Ahmed, J. von Eichborn et al., “SuperTarget goes quantitative: update on drug-target interactions,” Nucleic Acids Research, vol. 40, no. 1, Article ID gkr912, pp. D1113–D1117, 2012.
- S. Kim, P. A. Thiessen, E. E. Bolton et al., “PubChem substance and compound databases,” Nucleic Acids Research, vol. 44, no. 1, pp. D1202–D1213, 2016.
- A. P. Davis, C. J. Grondin, R. J. Johnson et al., “The Comparative Toxicogenomics Database: Update 2017,” Nucleic Acids Research, vol. 45, no. 1, pp. D972–D978, 2017.
- S. Tweedie, M. Ashburner, K. Falls et al., “FlyBase: enhancing Drosophila Gene Ontology annotations,” Nucleic Acids Research, vol. 37, no. 1, pp. D555–D559, 2009.
- M. Kanehisa and S. Goto, “KEGG: kyoto encyclopedia of genes and genomes,” Nucleic Acids Research, vol. 28, no. 1, pp. 27–30, 2000.
- D. Szklarczyk, A. Franceschini, S. Wyder et al., “STRING v10: protein-protein interaction networks, integrated over the tree of life,” Nucleic Acids Research, vol. 43, pp. D447–D452, 2015.
- A. Chatr-Aryamontri, B.-J. Breitkreutz, R. Oughtred et al., “The BioGRID interaction database: 2015 update,” Nucleic Acids Research, vol. 43, no. 1, pp. D470–D478, 2015.
- T. S. Keshava Prasad, R. Goel, K. Kandasamy et al., “Human protein reference database—2009 update,” Nucleic Acids Research, vol. 37, no. 1, pp. D767–D772, 2009.
- R. Saito, M. E. Smoot, K. Ono et al., “A travel guide to Cytoscape plugins,” Nature Methods, vol. 9, no. 11, pp. 1069–1076, 2012.
- Y. Tang, M. Li, J. Wang, Y. Pan, and F.-X. Wu, “CytoNCA: a cytoscape plugin for centrality analysis and evaluation of protein interaction networks,” BioSystems, vol. 127, pp. 67–72, 2015.
- X. He and J. Zhang, “Why do hubs tend to be essential in protein networks?” PLoS Genetics, vol. 2, no. 6, p. e88, 2006.
- G. D. Bader and C. Hogue, “An automated method for finding molecular complexes in large protein interaction networks,” BMC Bioinformatics, vol. 4, 2003.
- H. Han, H. Shim, D. Shin et al., “TRRUST: A reference database of human transcriptional regulatory interactions,” Scientific Reports, vol. 5, Article ID 11432, 2015.
- E. Herédi, J. Végh, L. Pogácsás et al., “Subclinical cardiovascular disease and it's improvement after long-term TNF-α inhibitor therapy in severe psoriatic patients,” Journal of the European Academy of Dermatology and Venereology, vol. 30, no. 9, pp. 1531–1536, 2016.
- D. H. Solomon, J. R. Curtis, K. G. Saag et al., “Cardiovascular risk in rheumatoid arthritis: comparing tnf-α blockade with nonbiologic DMARDs,” American Journal of Medicine, vol. 126, no. 8, pp. 730.e9–730.e17, 2013.
- S. L. Westlake, A. N. Colebatch, J. Baird et al., “Tumour necrosis factor antagonists and the risk of cardiovascular disease in patients with rheumatoid arthritis: a systematic literature review,” Rheumatology, vol. 50, no. 3, pp. 518–531, 2011.
- H. Zhang, Y. Park, J. Wu et al., “Role of TNF-α in vascular dysfunction,” Clinical Science, vol. 116, no. 3, pp. 219–230, 2009.
- Y. Zhang, X. Yang, F. Bian, and et al., “TNF-α promotes early atherosclerosis by increasing transcytosis of LDL across endothelial cells: Crosstalk between NF-κB and PPAR-γ,” Journal of Molecular and Cellular Cardiology, vol. 72, pp. 85–94, 2014.
- C. Ferrandi, R. Ballerio, P. Gaillard et al., “Inhibition of c-Jun N-terminal kinase decreases cardiomyocyte apoptosis and infarct size after myocardial ischemia and reperfusion in anaesthetized rats,” British Journal of Pharmacology, vol. 142, no. 6, pp. 953–960, 2004.
- G. Milano, S. Morel, C. Bonny, and et al., “A peptide inhibitor of c-Jun NH2-terminal kinase reduces myocardial ischemia-reperfusion injury and infarct size in vivo,” American Lournal of Physiology, Heart and Circulatory Physiology, vol. 292, no. 4, pp. H1828–H1835, 2007.
- H. Yasumoto, S. Kim, Y. Zhan et al., “Dominant negative c-Jun gene transfer inhibits vascular smooth muscle cell proliferation and neointimal hyperplasia in rats,” Gene Therapy, vol. 8, no. 22, pp. 1682–1689, 2001.
- S. Saliques, J.-R. Teyssier, C. Vergely et al., “Smoking and FOS expression from blood leukocyte transcripts in patients with coronary artery disease,” Atherosclerosis, vol. 219, no. 2, pp. 931–936, 2011.
- X. Palomer, E. Capdevila-Busquets, G. Botteri, and et al., “miR-146a targets Fos expression in human cardiac cells,” Disease Models and Mechanisms, vol. 8, no. 9, pp. 1081–1091, 2015.
- Y. N. Yang, J. Y. Zhang, and Y. T. Ma, “GW25-e3074 94 ATTG Insertion/Deletion Polymorphism of the NFKB1 Gene Is Associated with Coronary Artery Disease in Han and Uygur Women in China,” Genetic Testing and Molecular Biomarkers, vol. 18, no. 6, pp. 430–438, 2014.
- H. Lai, Q. Chen, X. Li et al., “Association between genetic polymorphism in NFKB1 and NFKBIA and coronary artery disease in a Chinese Han population,” International Journal of Clinical and Experimental Medicine, vol. 8, no. 11, pp. 21487–21496, 2015.
- V. Boccardi, M. R. Rizzo, R. Marfella et al., “-94 ins/del ATTG NFKB1 gene variant is associated with lower susceptibility to myocardial infarction,” Nutrition, Metabolism and Cardiovascular Diseases, vol. 21, no. 9, pp. 679–684, 2011.
- S. Dong, W. Ma, B. Hao, and et al., “microRNA-21 promotes cardiac fibrosis and development of heart failure with preserved left ventricular ejection fraction by up-regulating Bcl-2,” International Journal of Clinical and Experimental Pathology, vol. 7, no. 2, pp. 565–574, 2014.
- Y. Tang, J. Zheng, Y. Sun, Z. Wu, Z. Liu, and G. Huang, “MicroRNA-1 regulates cardiomyocyte apoptosis by targeting Bcl-2,” International Heart Journal, vol. 50, no. 3, pp. 377–387, 2009.
Copyright © 2017 Hua Li and Hongying Liu. 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.