Network Pharmacology Analysis and Experimental Verification Strategies Reveal the Action Mechanism of Danshen Decoction in Treating Ischemic Cardiomyopathy

Liu, Mengnan; Yuan, Gang; Luo, Gang; Guo, Xin; Chen, Mingtai; Yang, Huayi; He, Fan; Yang, Tingfu; Zhang, Xinyue; Wu, Qibiao; Zhou, Hua; Yang, Sijin

doi:https://doi.org/10.1155/2022/7578055

Evidence-Based Complementary and Alternative Medicine

On this page

Abstract Introduction Methods Results Discussion Conclusion Abbreviations Data Availability Disclosure Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 7578055 | https://doi.org/10.1155/2022/7578055

Network Pharmacology Analysis and Experimental Verification Strategies Reveal the Action Mechanism of Danshen Decoction in Treating Ischemic Cardiomyopathy

Mengnan Liu,^1,2,3Gang Yuan,¹Gang Luo,^2,3Xin Guo,¹Mingtai Chen,^1,4Huayi Yang,¹Fan He,⁵Tingfu Yang,^2,3Xinyue Zhang,⁶Qibiao Wu,^1,7Hua Zhou ,⁵and Sijin Yang^1,2

Academic Editor: Jing Zhao

Received17 Feb 2022

Accepted13 Apr 2022

Published02 May 2022

Abstract

Background. Danshen Decoction comprises Salvia miltiorrhiza, Santalum album, and Amomum villosum. It can promote blood circulation and remove blood stasis, and is commonly used in the treatment of gastric and duodenal ulcers, coronary heart disease, angina pectoris, etc. This research is based on network pharmacology and is experimentally verified to explore the potential mechanism of Danshen Decoction in the treatment of ischemic cardiomyopathy (ICM). Methods. The effective components and targets of Danshen Decoction were firstly extracted from Traditional Chinese Medicine Systems Pharmacology (TCMSP) Database and Analysis Platform, the drug-component-target-disease network was then constructed, the protein-protein interaction (PPI) network was constructed, the Gene Ontology (GO) enrichment analysis was carried out, and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway was analyzed in order to find the potential active components and therapeutic mechanisms. Finally, the in vitro hypoxia/reoxygenation model in H9c2 cells was established to verify the predicted active components and therapeutic mechanisms. Results. The results showed that Danshen Decoction has 67 potential active components and 109 therapeutic targets in treating ICM. These targets were rich in a variety of gene functions and different signaling pathways; the main gene targets include TP53, c-Jun, and Akt1. Go enrichment analysis showed that response to drug, membrane raft, and G protein-coupled amine receiver activity rank first in each process, and the main signaling pathways include PI3K-Akt signaling pathway. Through molecular docking and experimental verification of the major active components and core therapeutic targets, the active components of Danshen Decoction demonstrated an ability to reduce the cell damage caused by hypoxia/reoxygenation in H9c2 cells by regulating the core therapeutic target including Akt1, c-Jun, and TP53. Conclusion. Danshen Decoction has the effect of treating ICM in multiple ways, which is consistent with the results of network pharmacology. This laid a foundation for further study in exploring the active principles and pharmacological mechanism of Danshen Decoction.

1. Introduction

Ischemic cardiomyopathy (ICM) is diffuse myocardial fibrosis caused by long-term myocardial ischemia. It belongs to the advanced stage of coronary heart disease (CHD) [1]. Its pathophysiological basis is that coronary atherosclerosis causes myocardial ischemia and hypoxia, as well as myocardial cell reduction, necrosis, myocardial fibrosis, and myocardial scar formation; it produces clinical syndromes similar to primary dilated cardiomyopathy, often manifested as angina pectoris, heart failure, and arrhythmia [2]. It was reported that in North America and Europe, ischemic heart disease is the main cause of more than 50% of heart failure (HF) patients, and HF affects more than 26 million people worldwide [1, 3]. ICM is mainly treated with cardiotonic, diuretic, lipid-lowering, symptomatic treatment of heart failure combined with anti-atherosclerosis and surgery, if necessary. Patients are often hospitalized repeatedly due to heart failure or arrhythmia, and the prognosis is poor [4].

Traditional Chinese medicine (TCM) has unique advantages in the treatment of myocardial ischemia. Danshen Decoction, composed of Salvia miltiorrhiza Bunge, root, and rhizome (Salviae Miltiorrhizae Radix et Rhizoma) (Danshen), Santalum album L., lignum (Santali Albi Lignum) (Tanxiang) and Amomum villosum Lour, fruit (Amomi Fructus) (Sharen), is commonly used in TCM for the treatment of “different pains in the heart and abdomen” [5]. Danshen Decoction was first recorded in the Chinese Qing Dynasty medical book “Shi Fang Ge Kuo.” The ratio of these three herbs is 10 : 1 : 1, each dose of medicine is 44.64 g and decocted with water [6]. It can promote blood circulation and remove blood stasis, promote qi and relieve pain, and is currently often used in the treatment of gastric and duodenal ulcers, coronary heart disease, angina pectoris, etc. Clinical and basic studies have verified that it can inhibit blood blot, and promote fibrinolysis, anti-thrombosis, and anti-atherosclerosis [5, 7, 8]. However, there are many effective ingredients in TCM decoctions; the curative effect has been clear, and the understanding of the ingredients that exert pharmacological effects is still quite limited. Although there have been some reports on the study of Danshen Decoction or Salvia miltiorrhiza in the treatment of diabetic cardiomyopathy, angina pectoris, and CHD [5, 9], we still have the value of perfecting the research on the related mechanisms of Danshen Decoction in the treatment of ICM.

Network pharmacology is the use of computers, high-throughput omics data analysis, and other technologies to predict the pharmacological mechanism of Chinese medicine, which provides a new way for the study of action mechanism of Chinese medicine prescriptions, and the steps of our research also refer to the requirements of Network Pharmacology Evaluation Method Guidance [10–12]. This study uses a large database to screen the effective bioactive components and targets of Danshen Decoction, analyzes its key targets and signal pathways in the treatment of ICM, verifies it with molecular docking and in vitro experiments for further research on the Chinese medicine decoction, and provides guidance for disease treatment. The detailed workflow of this research is shown in Figure 1.

2. Methods

2.1. Network Pharmacology Study

2.1.1. Chemical Constituents Screening

All the effective ingredients of Chinese herbal medicine contained in Danshen Decoction, including Salvia miltiorrhiza, Santalum album, and Amomum villosum, were retrieved through the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP, https://tcmspw.com/index.php) [13, 14]; and then the candidate active compounds were screened by oral bioavailability (OB) ≥ 30% and drug similarity (DL) ≥ 0.18 as conditions, where OB refers to the relative amount and rate of absorption of the drug into the blood circulation by the body after oral administration and DL reflects the similarity between the specific functional gene in the compound and the known drug, and both are of great significance for the evaluation of the activity of the chemical components of TCM [15].

2.1.2. Genes Related to ICM Retrieval

The keyword “ischemic cardiomyopathy” was entered into Genecards database (https://www.genecards.org/), OMIM database (https://omim.org/), and DrugBank database (https://www.drugbank.ca/) to collect disease-related genes with PharmGkb database (https://www.pharmgkb.org/), and all collected genes were combined and deduplicated.

2.1.3. Danshen Decoction-ICM Co-Action Targets Screening and Regulatory Network Visualization

The Uniprot database (https://www.uniprot.org/) was used to search for the human gene name corresponding to the therapeutic targets of the compounds contained in the Danshen Decoction, which were then compared with the protein target gene of ICM screened in the previous step. Then, the common chemical components and common protein target genes of Danshen Decoction-ICM were introduced into Cytoscape 3.8.2 software for visualization. The targets were presented in the form of a network diagram, in which compounds and targets were represented by nodes, and the interaction relationship between nodes were represented by edges to screen core compounds. The STRING data platform (https://string-db.org/) was used to construct the protein-protein interaction (PPI) network of the target protein of Danshen Decoction in the treatment of ICM [10]. The obtained intersection targets were imported into the STRING database, the species was set to “Same Homo sapiens,” the minimum interaction threshold was set to “0.9 medium confidence,” and “hide disconnected nodes” was set, and the resulting network was saved in TSV format for future use [16]. The TSV file was imported into Cytoscape 3.8.2 software, CytoNCA plugin was used and select “Betweenness,” “Closeness,” “Degree,” “Eigenvector,” “Local Average Connectivity-Based Method,” and “Network” as the scoring basis. The larger these values, the more important the position of the node in the network [17]. After obtaining the scoring file, all genes greater than the median value were selected to obtain the core network, and the screening was repeated 3 times to obtain the core gene.

2.1.4. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes Enrichment Analysis

The Bioconductor database (https://www.bioconductor.org/) was used to query the gene IDs of common target genes for drugs-diseases. The R language was used to install the relevant installation package of the Bioconductor platform, P value = 0.05, q value = 0.05 were set for Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, respectively, and the relevant pathway diagrams have been queried [18].

2.1.5. Molecular Docking Verification

The obtained core gene was used as a protein receptor. The 2D structure of the relevant small molecule ligand was searched in the PubChem database (https://pubchem.ncbi.nlm.nih.gov/), and the ChemOffice software was used to convert it into a 3D structure; protein receptors were searched from PDB database (http://www.rcsb.org/); the Auto Dock Vina docking model was used. The lower the binding energy of the ligand and the receptor, the more stable the binding between the two. The binding energy ≤−5.0 kcal/mol⁻¹ is usually used as the standard to predict the binding between the core active ingredient and the core target [19].

2.2. Experimental Verification

2.2.1. Reagents Used in Validation

Tanshinone IIA, luteolin, and β-sitosterol were purchased from Chengdu Must Bio-Technology Co., Ltd. (Chengdu, China), with a purity of >98%, and the catalog numbers were A0057, A0108, and A0197, respectively. Propranolol hydrochloride was purchased from Shanghai Macklin Biochemical Co., Ltd. (Shanghai, China), with a purity of >98%, and the catalog number was P854704. TP53, c-Jun, p-c-Jun (Ser73), Akt1, and p-Akt1 (Thr308) antibodies were purchased from Affinity Biosciences Co., Ltd. (USA), the catalog numbers were DF7238, AF6090, AF3095, AF4718, and AF0832, respectively. Lactate dehydrogenase (LDH) and Malondialdehyde (MDA) kits were purchased from Beijing Solarbio Science & Technology Co., Ltd. (Beijing, China), and the catalog numbers were BC0685 and CA1410, respectively.

2.2.2. Cells Culture

H9c2 cells were purchased from ATCC (Rockville, MD, United States) and cultivated under conventional conditions (37°C, 95% air/5% CO₂) using high-sugar DMEM medium with 100 μg/mL streptomycin, 100 U/mL penicillin, and 10% fatal bovine serum (FBS). The cells were passaged every 2–3 days according to cell status.

2.2.3. Cell Models

The hypoxia chamber (Stem Cell Technology, USA) was used to establish a hypoxia/reoxygenation (H/R) model to create an in vitro hypoxic environment. The method refers to our previously published literature [20, 21]. 1 × 10⁵ cells were seeded in a six-well plate and cultured in a normal culture environment (95% air/5% CO₂) until the cell confluence reached 80%, and then washed twice with KRB buffer (composition in mM: NaCl 115, KCl 4.7, CaCl 2.5, KH₂PO 4.2, MgSO₄ 1.2, NaHCO₃ 24, HEPES 10, pH 7.4), pre-equilibrate with N₂ overnight at 4°C. Then, 100 μL KRB and 0.01% (w/v) bovine serum albumin (BSA) were added to each well and the cells were placed in the hypoxia chamber and replaced the gas with N₂ at a flow rate of 30 mL/min for 5 min, then all the connectors of the hypoxia chamber were closed, and the chamber was moved into the 37°C incubator for 2 h. After that, the KRB buffer in the well was replaced with fresh complete medium and incubated in a normal culture environment for 1 h. Dosing intervention should be started 12 h before modeling and continue until the end of modeling.

2.2.4. Cell Viability

The MTT method was used to detect the effects of the selected compounds on the proliferation of H9c2 cells. H9c2 cells in logarithmic growth phase were seeded in 96-well plates at a concentration of 5 × 10³ cells. After 24 h, the medium was changed and standards were added with concentrations of 1, 10, and 100 μmol/L and incubated for 12 h, repeated 5 times for each group. After incubating for modeling, 20 μL of MTT solution was added and incubated at 37°C for 4 h. The supernatant was aspirated and 150 μL SDS hydrochloric acid was added to dissolve the crystals; microplate reader was used to detect the OD value at a wavelength of 490 nm. Cell proliferation rate = average value of experimental group/average value of control group × 100%. More accurate effective concentration was further explored based on the above experimental results. Optical microscope was used to observe cell morphology observation.

2.2.5. Investigation of Predicted Compounds in Protecting from H/R Injury

H9c2 cells were processed according to the method in the section “Cell models.” They were divided into control group (cultivating under normal conditions), H/R model group, H/R + compound group (adding different compounds to incubate for 12 h during the establishment of H/R model), and H/R + propranolol group (incubating with positive control drug for 12 h during the establishment of H/R model). The process was repeated three times for each group, the supernatant was collected, and the detection of LDH and MDA were performed in accordance with the kit instructions.

2.2.6. Western Blot

Protein extraction was completed on ice after modeling. The concentration of the protein sample was determined by the Bradford method, and the protein loading buffer was added in proportion to the denaturation in a metal bath at 100°C for 10 min. The proteins were separated by SDS-bisacrylamide gel electrophoresis and then transferred onto the PVDF membrane. The TBST solution containing 5% BSA was used for blocking, then TP53, c-Jun, p-c-Jun, Akt1, and p-Akt1 (1 : 1000) antibodies were added and incubated at 4°C overnight. After washing the next day, the secondary antibody (1 : 5000) was added and incubated at room temperature for 2 h. After washing again, the Amersham ImageQuant 800 system (Global Life Sciences Solutions USA LLC, Marlborough, MA, USA) was used to acquire images and analyze the gray values of the bands.

2.3. PCR

The Trizol method was used to extract total RNA from cells. After reverse transcription and amplification, fluorescent quantitative PCR instrument (Roche, LightCycler 480II) was used for routine melting curve analysis to determine the Ct value, and GAPDH was used as an internal reference for control. The primers needed for the experiment were synthesized by Shanghai Shenggong Company. The primer sequences used in the present study were listed as follows: Akt1-q-F: CACAGGTCGCTACTATGCCATGAAG, Akt1-q-R: GCAGGACACGGTTCTCAGTAAGC, TP53-q-F: GTACCGTATGAGCCACCTGAG,

TP53-q-R: TCCAGCGTGATGATGGTAAG, c-Jun-q-F: GTCCTCCATAAATGCCTGTTCC, c-Jun-q-R: GATGCAACCCACTGACCAGAT, GAPDH-q-F: GGACCTCATGGCCTACATGG, GAPDH-q-R:TAGGGCCT CTC TTGCTCAGT.

2.4. Statistical Analysis

GraphPad Prism (Ver 8.2.1, GraphPad Software, San Diego, CA, USA) was used to perform statistical analysis on the results. The data were expressed as mean ± standard deviation. One-way ANOVA analysis of variance was used in statistics. The LSD-t test was used for pairwise comparisons between groups, and P < 0.05 was considered as statistically significant.

3. Results

3.1. Identification of Putative Ingredient Targets

The compounds contained in Salvia miltiorrhiza, Santalum album, and Amomum villosum were 202, 70, and 165, respectively. According to the characteristics of the OB and DL of the compounds, 58, 3, and 10 compounds were screened out; these compounds correspond to 932, 78, and 100 targets, respectively; the full name of targets were converted to Gene Symbol through analysis and conversion in the Uniprot database, the compounds contained in Salvia miltiorrhiza, Santalum album, and Amomum villosum correspond to 824, 89, and 63 genes, respectively. The active compounds of each component herb involved in this study are shown in Table 1, and the detailed information of the hypothetical components of the target genes are shown in Table S1.

3.1.1. Identification of Disease-Related Genes

There are 2092, 8, 99, and 115 ICM-related gene targets found in the four databases of GeneCards, OMIM, PharmGkb, and DrugBank, respectively, and the targets obtained from the four databases were combined and deduplicated to determine them as targets related to ICM, and a total of 2223 ICM targets were obtained. The genes of these databases were integrated in Table S2. The Venn diagram of overlapping genes is shown in Supplementary Figure S1.

3.1.2. Network Visualization

The Venny 2.1 online tool was used to find the interconnection of 136 therapeutic targets contained in Danshen Decoction and 2223 ICM-related targets, and a total of 109 overlapping genes were obtained; the Venn diagram is shown in Figure S2. The active compounds of Danshen Decoction and their corresponding disease-drug overlapping targets were introduced into Cytoscape 3.8.2 to form a network diagram of Danshen Decoction-compounds-ICM targets as shown in Figure 2. The network consists of 169 nodes and 737 edges. The nodes represent the active compounds and target proteins that meet the relevant targets, and the edges represent the interaction between the active compounds and the target proteins. Among these targets, they were ranked by the number of interconnections as shown in Table S3. Luteolin involves 43 targets, tanshinone IIA involves 35 targets, and β-sitosterol involves 24 targets. The top 5 of these compounds were used as candidate compounds to verify their relationship with the targets.

Subsequently, 109 drug-disease overlapping genes were submitted to the STRING database, and 96 PPI network target genes were obtained. The network relationship consisting of 96 nodes and 388 edges is shown in Figure 3(a). The edges represent the degree of correlation between proteins. The TSV file was downloaded in STRING and imported into Cytoscape 3.8.2 software, CytoNCA was used to perform PPI network analysis on the TSV file to obtain the core gene network, then 3 core genes were finally screened out. The screening process is shown in Figure 3(b), and the scoring data are shown in Table S4. These genes include TP53, Akt1, and Jun. In addition, some nodes with high degree such as EGFR, BCL2, VEGFA, FOS, CASP3, and CASP8 were also considered to be important targets for the treatment of ICM.

(a)

(b)

3.1.3. GO and KEGG Enrichment Analyses

To further clarify the possible effects of the 109 candidate targets, GO enrichment analysis and KEGG biological pathway analysis were carried out on the 109 candidate targets of Danshen Decoction in the treatment of ICM. There are 2005 entries in biological process (BP), 930 entries in cellular component (CC), and 664 entries in molecular function (MF). They rank first in response to drug, membrane raft, and G protein-coupled amine receiver activity. The top 10 entries in each process are shown in Figure 4. In BP, the response to drugs, the response to metal ion ranked first; the regulation of cell membranes and cell synapses are ranked first in CC; in MF, it is closely related to DNA-binding transcription factor binding.

(a)

(b)

KEGG biological pathway analyzed a total of 159 related signal pathways, the top three related signal pathways include PI3K-Akt signaling pathway, human cytomegalovirus infection, and hepatitis B, and the top 30 entries are shown in Figure 5. The results showed that the mechanism of Danshen Decoction in the treatment of ICM mainly relates to multiple tumor-related pathways. The PI3K-Akt signal pathway with the highest degree of enrichment was selected to draw a pathway diagram, as shown in Figure S3.

(a)

(b)

3.1.4. Molecular Docking Verification

The TP53, Akt1, and c-Jun screened out in the PPI network were used as core receptors for molecular docking with corresponding compounds. Among them, TP53 corresponds to tanshinone IIA and luteolin; Akt1 corresponds to luteolin; c-Jun corresponds to tanshinone IIA, luteolin, and β-sitosterol. Next, Pymol 2.3.2 software was used for visual analysis of molecular docking. When the three compounds were docked with the corresponding targets, the analysis found that multiple active sites could form hydrogen bonds with different amino acids, and the average affinity of this docking was less than −5.0 kcal/mol⁻¹. However, although tanshinone IIA and β-sitosterol cannot bind to JUN at the hydrogen bonding site, the potential binding of other sites were still shown in the simulated docking. Therefore, they all showed good binding activity to the target, which may be the therapeutic agent for the corresponding target, as shown in Table 2.

3.1.5. Intervention Concentration Screening

Three potential effective compounds including tanshinone IIA, luteolin, and β-sitosterol (Table 2) were selected for in vitro experiments to verify the virtual screening, and propranolol was used as the positive control drug. The concentrations of 1, 10, and 100 μM were firstly used and the results showed that when the concentration reached 100 μM, the cell activity decreased significantly, so 10 μM was selected for further experiment (Figure 6(a)). Then the improvement rate under different concentrations was studied to obtain the best dosage. After H/R intervention, we got similar results, the optimal dosages of tanshinone IIA, luteolin, β-sitosterol and propranolol were 10 μM respectively, when the concentration exceeds 100 μM, the relative cell viability of each test decreases to varying degrees (Figure 6(b)). In terms of cell morphology, the control group has a large number of cells, complete morphology, and monolayer cluster growth; the model group has a large number of nuclear shrinkages, the cells become rounded, and floated in the medium; after the intervention of the compounds, the cell morphology has improved to varying degrees with fewer floating cells (Figure 6(c)).

(a)

(b)

(c)

3.1.6. Predicted Effective Compounds Protect H9c2 Cells from H/R Injury

Based on the findings of this study for predicting therapeutic targets and related signal pathways, the effect of the core compounds of Danshen Decoction on the H/R model was verified. After H/R treatment, the content of LDH and MDA in the cell supernatant was reduced by these compounds (P < 0.05), which reflects the different degrees of protection of these compounds on cells (Figure 7(a)). Then, the indicators of TP53, c-Jun, and Akt1 were used to detect protein expression, and the immunoblotting bands were shown in Figure 7(b). Compared with the H/R group, these compounds all had a significant inhibitory effect on the expression of p-c-Jun protein (P < 0.05), and had a significant enhancement effect on the expression of p-Akt1 protein (P < 0.05); however, these compounds did not show a significant effect on TP53 compared with the model group (P > 0.05) (Figure 7(c)). However, at the level of gene expression, we found that these compounds significantly increased Akt and decreased c-Jun, TP53 mRNA expression (P > 0.05) (Figure 7(d)). In short, the active compounds of Danshen Decoction can regulate Jun and Akt1 as protein targets to improve myocardial hypoxia.

(a)

(b)

(c)

(d)

Figure 7

The effect of predicted compounds on improving cell hypoxia in multiple ways. (a) The expressions of LDH and MDA from different groups (n = 3). (b) Western blot bands. (c) The expressions of TP53, p-c-Jun, and p-Akt from different groups were detected by Western blot assay (n = 3). (d) The expression of Akt, c-Jun, and TP53 mRNA in different groups was detected by PCR (n = 6). The results were presented as mean ± SD. ^∗P < 0.05, ^∗∗P < 0.01 vs. Model group.

4. Discussion

Network pharmacology is designed to study drug treatments with unclear targets or multiple targets to predict possible therapeutic effects [15]. Danshen Decoction is a mixture of 3 kinds of herbs with multiple compounds and targets, the network pharmacology study of Danshen Decoction in the treatment of ICM conforms to the above point of view, and the compounds we studied were closed to the protective effect of the positive control drug propranolol on the H/R mode; this provides a reference for future in-depth research. TCM with Danshen as the main ingredient has been verified in some studies to have the effect of treating or preventing cardiovascular diseases. For example, Danshen injection can significantly prevent myocardial fibrosis, cardiac hypertrophy, hemodynamic deterioration, and systolic and diastolic dysfunction characterized by failed hearts [8]. Danshen-Gegen Decoction can inhibit cell apoptosis induced by H/R by inhibiting the transition of mitochondrial permeability [22]. Many studies published in Chinese journals also showed that Danshen Decoction has a definite effect on myocardial ischemia [23]. Therefore, it is not necessary for us to carry out repeated in vitro experimental verification. We focused on these active ingredients and the most closely interacting targets screened by network pharmacology in order to illustrate the material basis of Danshen Decoction against myocardial ischemia.

Three effective ingredients were screened out in this study. Luteolin is a widely distributed flavonoid compound that can be found in a variety of plants. A large number of studies have found that luteolin is a potential cardioprotective agent, and it also has many evidences of cardioprotection in epidemiology [24]. β-sitosterol can also be found in a variety of plants. The current research showed that the content of sitosterol in the circulating blood may have a certain relationship with the cardiovascular disease [25]. An in vitro study by Huang et al. showed that β-sitosterol regulates the cell cycle by promoting cellular glutathione production and protects H9c2 cells from apoptosis induced by H/R [26]. Tanshinone IIA is obviously the most specific compound among these three herbs. It is the main fat-soluble component of Salvia miltiorrhiza. Modern pharmacological studies have shown that tanshinone IIA has anti-inflammatory and antioxidant activities and has a wide range of treatments for cardiovascular diseases [27]. In this experimental verification, we were able to determine the role of the compound itself, which on the other hand also shows that the method of network pharmacology has a certain degree of reliability. However, it cannot be ignored that although these three compounds have been confirmed to have the effect of saving cell death in our experimental verification, the dose-response relationship cannot be well reflected due to the large drug dose gap set in our study. But there is no doubt that the therapeutic effects of these compounds have been demonstrated.

GO function analysis showed that the core targets of Danshen Decoction in the treatment of ICM were TP53, Akt1, Jun, etc. TP53 is a human tumor suppressor gene, and it has been confirmed that the expression of this gene is closely related to the suppression of tumor cells [28]. Regulating the expression of TP53 protein can promote the repair process of heart endothelial cells and help vascular remodeling to reduce cardiovascular damage [29]. In our study, the protective effect of TP53 on cardiomyocytes was not significant, but its trend was obvious. We considered that this was related to our short-term modeling method, and the protein expression did not reach the peak at this time. Akt involves in the occurrence and development of inflammation, cancer, diabetes, and cardiovascular diseases; the functions include the regulation of cell cycle and transcription [30]. Akt1 plays an important role in multiple signaling pathways such as PI3K-Akt, NF-κB, and Toll-like receptors. Research has shown that the expression and activation of Akt1 are involved in the progression of myocardial fibrosis [31]; Jun can serve as an important transcription factor in the cell and widely involves in tumor regulation, smooth muscle proliferation, and apoptosis [32, 33]. According to the enrichment of targets, we also found that some other important targets, such as CASP3, play a central role in the execution of apoptosis, and can regulate cell proliferation, survival, migration, and other physiological functions to protect damaged myocardium [34]. NF-κB participates in the response of cells to external stimuli, such as cytokines, radiation, heavy metals, viruses, etc., and plays a key role in the process of cellular inflammation and immune response [35]. IL-1β, as a downstream inflammatory factor of inflammation-related signal pathways, has been proved to have extensive regulation on the inflammatory response of blood vessels and myocardium [35, 36]. MAPK14 and MAPK1 can regulate various physiological and pathological processes and the inflammatory responses of cells through the MAPK signaling pathway [37]. VEGFA is the central factor in the regulation of neovascularization, which can accelerate the proliferation and differentiation of vascular endothelium [38]. From the verification of molecular docking, it was found that the core active ingredients predicted by Danshen Decoction have good binding activity with the core target for the treatment of ICM, which is consistent with the results of network pharmacology screening and has potential therapeutic significance. Our experiments also verified that TP53, Akt1, and Jun were involved in the repair of cardiomyocyte H/R model, which is consistent with the trend of research reported in the literature [5, 39–42]. Reactive oxygen species (ROS) are highly active molecules in vivo. In the human body, ROS acts on lipid to produce peroxidation reaction, the end product of oxidation is MDA, which has cytotoxicity and causes cross-linking polymerization of proteins, nucleic acids, and other life macromolecules. In clinical practice, LDH is often used as an auxiliary diagnostic index of acute myocardial infarction. The determination of three indexes was used to investigate the protective effect of core compounds on the H/R model. The experiment showed that these indexes significantly decreased compared with the model group after using the compounds intervention, indicating that the predicted core compounds to protect heart is likely the active ingredient of Danshen Decoction.

The enrichment of the KEGG pathway of core genes showed that the core targets were involved in multiple pathways related to tumors. Others include PI3K-Akt, fluid shear stress and atherosclerosis, IL-17 signaling pathway, apoptosis, HIF-1, tumor necrosis factor, EGFR, etc. These signaling pathways have been studied in depth in cardiovascular diseases [39, 41]. PI3K-Akt signaling pathway enriched most genes in this study, which involves in the regulation of cell proliferation and plays an important role in combating myocardial H/R injury and inhibiting cell apoptosis. Regarding the preliminary verification of this pathway, we have found that Akt1 is involved in the repair process of cell H/R injury, demonstrating the reliability of this network pharmacology study.

Although this study has been partially verified, the following problems still exist due to the limitations of the network pharmacology method: First, the components directly extracted from TCMSP may be inconsistent with the exact compounds absorbed by the patient; second, the compounds included in this study need to meet the OB and DL limits, and there is excessive screening; third, due to the diversity of the corresponding targets of the compound, there may be errors in the enrichment analysis of GO and KEGG; fourth, with the continuous update of the data of each database, the results of future research may be more abundant. However, there must be some errors in the current research. Such problems have been exposed in this study; for example, β-sitosterol in this study, the corresponding therapeutic targets were Akt1 only, and the corresponding targets of tanshinone IIA include TP53 and c-Jun. However, after experimental verification, we found that these compounds have regulatory effects on Akt1, TP53, and c-Jun, which may be bound up with the delay in updating the network pharmacology–related database. Furthermore, the interaction between multiple compounds was not considered in this study, and the most attractive part of traditional Chinese medicine itself is that many compounds play a wide range of regulatory roles, which should be paid attention to in future research. In addition, we did not detect the contents of various predicted compounds in Danshen Decoction, and the data of these substances were obtained from the database. Unfortunately, the proportion of these three compounds verified by our experiments has not been obtained through literature review, which makes it impossible for us to mix these three compounds in proper proportion to observe their improvement on cell damage. Last but not least, verification was carried out at the level of in vitro experiments, and verification in animal models will also be required in future studies with higher levels of evidence. In future research, we will continue to reveal the mechanism of Danshen Decoction in the treatment of related diseases through in vitro and in vivo experiments.

5. Conclusion

Our research provides molecular evidence for reference to prove that Danshen Decoction has the effect of treating ischemic cardiomyopathy, and the mechanism of such treatment is the participation of multiple compounds, pathways, and targets. Its core components can regulate cytokine levels and improve hypoxia from multiple aspects, so as to protect the cardiomyocytes of ICM patients. Our network pharmacology research and experimental verification have laid the foundation for further study on the pharmacological mechanism of Danshen Decoction.

Abbreviations

Akt:	Protein kinase B
BCL:	B-cell lymphoma
BP:	Biological process
CASP:	Cysteinyl aspartate specific proteinase
CC:	Cellular component
CHD:	Coronary heart disease
DL:	Drug similarity
EGFR:	Epidermal growth factor receptor
GO:	Gene ontology
KEGG:	Kyoto encyclopedia of genes and genomes
HF:	Heart failure
ICM:	Ischemic cardiomyopathy
IL:	Interleukin
Jun:	Jun activation domain-binding protein
LDH:	Lactate dehydrogenase
MDA:	Malondialdehyde
MF:	Molecular function
OB:	Oral bioavailability
ROS:	Reactive oxygen species
TCM:	Traditional Chinese medicine
TCMSP:	Traditional Chinese medicine systems pharmacology database and analysis platform
TP53:	Tumor protein P53
VEGFA:	Vascular endothelial growth factor.

Data Availability

The datasets used or analyzed during the current study are available from the corresponding author.

Disclosure

The funder had no role in the study design, data analysis, or decision to publish. Mengnan Liu, GangYuan, Gang Luo, and Xin Guo are co-first authors.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Mengnan Liu, GangYuan, Gang Luo, and Xin Guo are equal contributors. Mengnan Liu and Gang Yuan were involved in conceptualization, methodology, data analysis, and manuscript writing. Gang Luo and Xing Guo performed experiment and data analysis. Mingtai Chen, Huayi Yang, and Tingfu Yang did experimental design. Fan He and Xinyue Zhang did investigation. Qibiao Wu, Hua Zhou, and Sijin Yang were in charge of funding acquisition, conceptualization, project design, and manuscript revision. All data were generated in-house, and no paper mill was used. All authors agree to be accountable for all aspects of work ensuring integrity and accuracy.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (82074378), Science and Technology Development Fund, Macao SAR (0001/2021/AKP, 0098/2021/A2), the Science and Technology Planning Project of Guangdong Province (2020B1212030008), the National Traditional Chinese Medicine Inheritance and Innovation Team (ZYYCXTD-C-202207), and the Southwest Medical University Project (2021ZKQN143).

Supplementary Materials

Figure S1. Venn diagram showing gene targets related to ICM. Figure S2. Venn diagram showing the overlapped genes between Danshen Decoction and ICM. Figure S3. PI3K-Akt signal pathway. Table S1. Eligible and putative ingredients with targeted genes. Table S2. ICM-related genes. Table S3. Sequence of compounds linkage gene targets. Table S4. Core gene target score. (Supplementary Materials)

References

B. R. Pagliaro, F. Cannata, G. G. Stefanini, and L. Bolognese, “Myocardial ischemia and coronary disease in heart failure,” Heart Failure Reviews, vol. 25, no. 1, pp. 53–65, 2020.
View at: Publisher Site | Google Scholar
N. Briceno, A. Schuster, M. Lumley, and D. Perera, “Ischaemic cardiomyopathy: pathophysiology, assessment and the role of revascularisation,” Heart, vol. 102, no. 5, pp. 397–406, 2016.
View at: Publisher Site | Google Scholar
G. Savarese and L. H. Lund, “Global public health burden of heart failure,” Cardiac Failure Review, vol. 3, no. 1, pp. 7–11, 2017.
View at: Publisher Site | Google Scholar
D. Gajanana, M. Shah, P. Junpapart, A. Romero-Corral, V. M. Figueredo, and B. Bozorgnia, “Mortality in systolic heart failure revisited: ischemic versus non-ischemic cardiomyopathy,” International Journal of Cardiology, vol. 224, pp. 15–17, 2016.
View at: Publisher Site | Google Scholar
K.-P. Yan, Y. Guo, Z. Xing et al., “Dan-Shen-Yin protects the heart against inflammation and oxidative stress induced by acute ischemic myocardial injury in rats,” Experimental and Therapeutic Medicine, vol. 3, no. 2, pp. 314–318, 2012.
View at: Publisher Site | Google Scholar
Q. Dong, H.-T. Yu, Y.-T. Dai et al., “Establish quality evaluation method based on standard decoction of Danshen extract,” Zhongguo Zhong Yao Za Zhi, vol. 42, no. 5, pp. 817–822, 2017.
View at: Publisher Site | Google Scholar
X.-D. MEIm, Y.-F. Cao, and Y.-Y. Che, “Danshen: a phytochemical and pharmacological overview,” Chinese Journal of Natural Medicines, vol. 17, pp. 59–80, 2019.
View at: Google Scholar
L. Wang, J. Yu, P. A. Fordjour et al., “Danshen injection prevents heart failure by attenuating post-infarct remodeling,” Journal of Ethnopharmacology, vol. 205, pp. 22–32, 2017.
View at: Publisher Site | Google Scholar
Z.-M. Li, S.-W. Xu, and P.-Q. Liu, “Salvia miltiorrhizaBurge (Danshen): a golden herbal medicine in cardiovascular therapeutics,” Acta Pharmacologica Sinica, vol. 39, no. 5, pp. 802–824, 2018.
View at: Publisher Site | Google Scholar
C. Huang, C. Zheng, Y. Li, Y. Wang, A. Lu, and L. Yang, “Systems pharmacology in drug discovery and therapeutic insight for herbal medicines,” Briefings in Bioinformatics, vol. 15, no. 5, pp. 710–733, 2014.
View at: Publisher Site | Google Scholar
S. Li, Y.-T. Chen, and Q.-Y. Ding, “Network pharmacology evaluation method guidance - draft,” WJTCM, vol. 7, 2021.
View at: Google Scholar
S. Li and B. Zhang, “Traditional Chinese medicine network pharmacology: theory, methodology and application,” Chinese Journal of Natural Medicines, vol. 11, no. 2, pp. 110–120, 2013.
View at: Publisher Site | Google Scholar
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, 2014.
View at: Publisher Site | Google Scholar
H. Yu, J. Chen, X. Xu et al., “A systematic prediction of multiple drug-target interactions from chemical, genomic, and pharmacological data,” PLoS One, vol. 7, no. 5, Article ID e37608, 2012.
View at: Publisher Site | Google Scholar
T.-T. Luo, Y. Lu, S.-K. Yan, X. Xiao, X. Rong, and J. Guo, “Network pharmacology in research of Chinese medicine formula: methodology, application and prospective,” Chinese Journal of Integrative Medicine, vol. 26, no. 1, pp. 72–80, 2020.
View at: Publisher Site | Google Scholar
N. T. Doncheva, J. H. Morris, J. Gorodkin, and L. J. Jensen, “Cytoscape StringApp: network analysis and visualization of proteomics data,” Journal of Proteome Research, vol. 18, no. 2, pp. 623–632, 2019.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
L. Chen, Y.-H. Zhang, S. Wang, Y. Zhang, T. Huang, and Y. D. Cai, “Prediction and analysis of essential genes using the enrichments of gene ontology and KEGG pathways,” PLoS One, vol. 12, no. 9, Article ID e0184129, 2017.
View at: Publisher Site | Google Scholar
P. Yang, H. He, S. Xu, P. Liu, and X. Bai, “Potential molecular target prediction and docking verification of hua-feng-dan in stroke based on network pharmacology,” Evidence-Based Complementary and Alternative Medicine, vol. 2020, Article ID 8872593, 12 pages, 2020.
View at: Publisher Site | Google Scholar
G. Dong, T. Chen, X. Ren et al., “Rg1 prevents myocardial hypoxia/reoxygenation injury by regulating mitochondrial dynamics imbalance via modulation of glutamate dehydrogenase and mitofusin 2,” Mitochondrion, vol. 26, pp. 7–18, 2016.
View at: Publisher Site | Google Scholar
Y.-H. Zuo, Q.-B. Han, G.-T. Dong et al., “Panax ginseng polysaccharide protected H9c2 cardiomyocyte from hypoxia/reoxygenation injury through regulating mitochondrial metabolism and RISK pathway,” Frontiers in Physiology, vol. 9, 2018.
View at: Publisher Site | Google Scholar
P. Y. Chiu, H. Y. Leung, P. K. Leong et al., “Danshen-Gegen decoction protects against hypoxia/reoxygenation-induced apoptosis by inhibiting mitochondrial permeability transition via the redox-sensitive ERK/Nrf2 and PKCε/mKATP pathways in H9c2 cardiomyocytes,” Phytomedicine, vol. 19, no. 2, pp. 99–110, 2012.
View at: Publisher Site | Google Scholar
Y. Wang, L. Xiao, W. Mu et al., “A summary and evaluation of current evidence for myocardial infarction with Chinese medicine,” Chinese Journal of Integrative Medicine, vol. 23, no. 12, pp. 948–955, 2017.
View at: Publisher Site | Google Scholar
D. Pan and D. Li, “At the crossroads from bench to bedside: luteolin is a promising pharmacological agent against myocardial ischemia reperfusion injury,” Annals of Translational Medicine, vol. 4, no. 23, p. 475, 2016.
View at: Publisher Site | Google Scholar
H. Tada, A. Nohara, A. Inazu, N. Sakuma, H. Mabuchi, and M. Kawashiri, “Sitosterolemia, hypercholesterolemia, and coronary artery disease,” Journal of Atherosclerosis and Thrombosis, vol. 25, no. 9, pp. 783–789, 2018.
View at: Publisher Site | Google Scholar
H. S. Wong, N. Chen, P. K. Leong, and K. M. Ko, “β-Sitosterol enhances cellular glutathione redox cycling by reactive oxygen species generated from mitochondrial respiration: protection against oxidant injury in H9c2 cells and rat hearts,” Phytotherapy Research, vol. 28, no. 7, pp. 999–1006, 2014.
View at: Publisher Site | Google Scholar
R. Guo, L. Li, J. Su et al., “Pharmacological activity and mechanism of tanshinone IIA in related diseases,” Drug Design, Development and Therapy, vol. 14, pp. 4735–4748, 2020.
View at: Publisher Site | Google Scholar
T. Soussi and K. G. Wiman, “TP53: an oncogene in disguise,” Cell Death & Differentiation, vol. 22, no. 8, pp. 1239–1249, 2015.
View at: Publisher Site | Google Scholar
A. A. Tzovaras, F. Gentimi, and M. Nikolaou, “Tumor protein p53 gene and cardiovascular disease,” Angiology, vol. 69, no. 8, pp. 736-737, 2018.
View at: Publisher Site | Google Scholar
P. Abeyrathna and Y. Su, “The critical role of Akt in cardiovascular function,” Vascular Pharmacology, vol. 74, pp. 38–48, 2015.
View at: Publisher Site | Google Scholar
B. T. O’Neill and E. D. Abel, “Akt1 in the cardiovascular system: friend or foe?” Journal of Clinical Investigation, vol. 115, no. 8, pp. 2059–2064, 2005.
View at: Publisher Site | Google Scholar
R. Eferl, M. Sibilia, F. Hilberg et al., “Functions of c-Jun in liver and heart development,” Journal of Cell Biology, vol. 145, no. 5, pp. 1049–1061, 1999.
View at: Publisher Site | Google Scholar
F. Wu, F. Gao, S. He, Y. Jiang, G. Luo, and Y. Xiao, “c-Jun promotes the survival of H9c2 cells under hypoxia via PTEN/Akt signaling pathway,” J Physiol Biochem, vol. 75, no. 4, pp. 433–441, 2019.
View at: Publisher Site | Google Scholar
B. Yang, D. Ye, and Y. Wang, “Caspase-3 as a therapeutic target for heart failure,” Expert Opinion on Therapeutic Targets, vol. 17, no. 3, pp. 255–263, 2013.
View at: Publisher Site | Google Scholar
M. Bartekova, J. Radosinska, M. Jelemensky, and N. S. Dhalla, “Role of cytokines and inflammation in heart function during health and disease,” Heart Failure Reviews, vol. 23, no. 5, pp. 733–758, 2018.
View at: Publisher Site | Google Scholar
S. Toldo and A. Abbate, “The NLRP3 inflammasome in acute myocardial infarction,” Nature Reviews Cardiology, vol. 15, no. 4, pp. 203–214, 2018.
View at: Publisher Site | Google Scholar
T. Yokota and Y. Wang, “p38 MAP kinases in the heart,” Gene, vol. 575, no. 2, pp. 369–376, 2016.
View at: Publisher Site | Google Scholar
M. Braile, S. Marcella, L. Cristinziano et al., “VEGF-A in cardiomyocytes and heart diseases,” International Journal of Molecular Sciences, vol. 21, no. 15, Article ID E5294, 2020.
View at: Publisher Site | Google Scholar
G.-M. Li, C.-L. Zhang, R.-P. Rui, B. Sun, and W. Guo, “Bioinformatics analysis of common differential genes of coronary artery disease and ischemic cardiomyopathy,” European Review for Medical and Pharmacological Sciences, vol. 22, no. 11, pp. 3553–3569, 2018.
View at: Publisher Site | Google Scholar
C.-Y. Su, Q.-L. Ming, K. Rahman, T. Han, and L.-P. Qin, “Salvia miltiorrhiza: traditional medicinal uses, chemistry, and pharmacology,” Chinese Journal of Natural Medicines, vol. 13, no. 3, pp. 163–182, 2015.
View at: Publisher Site | Google Scholar
Y.-S. Weng, W.-W. Kuo, Y.-M. Lin et al., “Danshen mediates through estrogen receptors to activate Akt and inhibit apoptosis effect of Leu27IGF-II-induced IGF-II receptor signaling activation in cardiomyoblasts,” Food and Chemical Toxicology, vol. 56, pp. 28–39, 2013.
View at: Publisher Site | Google Scholar
L. Zhang, L. Han, and X. Wang, “Exploring the mechanisms underlying the therapeutic effect of Salvia miltiorrhiza on diabetic nephropathy using network pharmacology and molecular docking,” Bioscience Reports, vol. 41, 2021.
View at: Google Scholar

Copyright

Copyright © 2022 Mengnan Liu 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.

PDF Download Citation

Download other formats

Order printed copies