Evidence-Based Complementary and Alternative Medicine
Volume 2015 (2015), Article ID 497314, 12 pages
http://dx.doi.org/10.1155/2015/497314

Research Article

The Mechanism Research of Qishen Yiqi Formula by Module-Network Analysis

Shichao Zheng, Yanling Zhang, and Yanjiang Qiao

Key Laboratory of TCM-Information Engineer of State Administration of TCM, School of Chinese Pharmacy, Beijing University of Chinese Medicine, Beijing 100102, China

Received 28 May 2015; Accepted 4 August 2015

Academic Editor: Lixing Lao

Copyright © 2015 Shichao Zheng 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.

Abstract.

Qishen Yiqi formula (QSYQ) has the effect of tonifying Qi and promoting blood circulation, which is widely used to treat the cardiovascular diseases with Qi deficiency and blood stasis syndrome. However, the mechanism of QSYQ to tonify Qi and promote blood circulation is rarely reported at molecular or systems level. This study aimed to elucidate the mechanism of QSYQ based on the protein interaction network (PIN) analysis. The targets’ information of the active components was obtained from ChEMBL and STITCH databases and was further used to search against protein-protein interactions by String database. Next, the PINs of QSYQ were constructed by Cytoscape and were analyzed by gene ontology enrichment analysis based on Markov Cluster algorithm. Finally, based on the topological parameters, the properties of scale-free, small world, and modularity of the QSYQ’s PINs were analyzed. And based on function modules, the mechanism of QSYQ was elucidated. The results indicated that Qi-tonifying efficacy of QSYQ may be partly attributed to the regulation of amino acid metabolism, carbohydrate metabolism, lipid metabolism, and cAMP metabolism, while QSYQ improves the blood stasis through the regulation of blood coagulation and cardiac muscle contraction. Meanwhile, the “synergy” of formula compatibility was also illuminated.

1. Introduction

Qishen Yiqi formula (QSYQ), consisting of Radix Salvia miltiorrhiza, Panax notoginseng, Dalbergia odorifera, and Astragalus membranaceus, has the effect of tonifying Qi, promoting blood circulation and relieving pain, and hence it has been widely used to treat the cardiovascular diseases with Qi deficiency and blood stasis syndrome [1, 2]. Pharmacological researches have shown that the mechanism of QSYQ is related to improve myocardial function, inhibit platelet aggregation, prevent enlargement of end-diastolic diameter, and slow down ventricular remodeling [3–7]. However, the mechanism of QSYQ was mostly elucidated macroscopically by pharmacology indexes of animal experiments or clinical trials. For example, Tong et al. [6] had reported that the myocardial protection function of QSYQ may relate to the reduction of myocardial cell apoptosis in adriamycin-induced cardiomyopathy animal model. The research of Cui et al. [7] had shown that, by detecting clinical indexes which included the right ventricular end-diastolic volume (RVEDV), end-systolic volume (RVESV), stroke volume (SV), and right ventricular ejection fraction (RVEF), QSYQ could significantly improve the right heart function on patients undergoing valve replacement. The findings of these studies explain the action mechanism of QSYQ to some extent, but the further study of QSYQ is still to be done. Up to now, the mechanism of QSYQ to tonify Qi and promote blood circulation is rarely reported at molecular or systems level. In this study, the mechanism of QSYQ was illuminated by the network analysis approach which has the advantage of evaluating TCM’s pharmacological effect as a whole unity at molecular level [8, 9].

Proteins are vital macromolecules, at both cellular and systematic levels, but they rarely act alone. And protein-protein interactions (PPIs) are major bearers of the biological process. So, protein interaction network (PIN) could provide the basis of understanding cellular organization and processes. The GO [10] project is a collaborative effort to construct ontologies which facilitate biologically meaningful annotation of gene products. It provides a collection of well-defined biological terms, spanning biological processes, molecular functions, and cellular components. GO enrichment is a common statistical method used to identify shared associations between proteins and annotations to GO. Module-network and GO analysis may provide an efficient way to illustrate the molecular mechanism of QSYQ.

In this study, a network analysis approach based on functional modules is applied to systematically illuminate the mechanism of QSYQ. The PINs of QSYQ were constructed by Cytoscape, while properties of scale-free, small word, and modularity were analyzed based on topological parameters. Then, the functional modules were identified by gene ontology (GO) enrichment analysis based on Markov Cluster (MCL) algorithm. This study aimed to provide an efficient way to elucidate the mechanism of QSYQ based on functional modules at the molecular level.

2. Materials and Methods

2.1. Targets Mining of Main Active Components of QSYQ

The main active components of QSYQ were used to study the mechanism of QSYQ. By literature retrieval from PubMed and CNKI database, the main active components of QSYQ were obtained based on the principles that components are the main efficacy components, have rich content, and can be absorbed into the blood. The information of the main active components of QSYQ is shown in Table 1.

Table 1: The information of the main active components of QSYQ.

Herbs
Active components
Reference

Salvia miltiorrhiza
Tanshinone IIA, cryptotanshinone, salvianolic acid A, salvianolic acid B, tanshinol, and protocatechuic aldehyde
[46–51]

Panax notoginseng
Dencichine, ginsenoside Rb1, ginsenoside Rg1, and notoginsenoside R1
[52–54]

Dalbergia odorifera
Butein, formononetin, isoliquiritigenin, nerolidol
[55–58]

Astragalus membranaceus
Calycosin, astragaloside Ι, formononetin, and astragaloside IV
[59, 60]

The targets’ information of main active components of QSYQ was obtained from two parts: pharmacophore virtual screening and the component-protein interaction database including ChEMBL (https://www.ebi.ac.uk/chembl/#) [11] and STITCH 3.1 (http://stitch.embl.de/) [12]. The 27 pharmacophore models which were applied to virtual screen were constructed by our laboratory team [13, 14]. ChEMBL is a manually curated chemical database which contains compound bioactivity data against drug targets. STITCH is a database in which every interaction has a confidence score, and the interactions with a confidence score > 0.7 were selected.

2.2. Network Construction of Single Herb and Formula

The PPIs information of targets was obtained from the online updated database of String 9.1 (http://string-db.org/) which has a confidence score for every protein interaction [15]. PPIs with a confidence score > 0.7 were applied to construct PIN using Cytoscape which is one of the most popular open-source software tools for the visual exploration of biomedical networks composed of protein, gene, and other types of interactions [16]. Every single herb network is formed only by PPIs involving proteins of this herb, and the formula network is formed only by PPIs involving proteins of this formula.

2.3. Network Analysis

The analysis of topological properties based on topological parameters has become very popular for gaining insight into the organization and structure of the resultant large complex networks [17–19]. Therefore, the topological parameters such as degree distribution, average shortest path, and clustering coefficient were analyzed by Network Analyzer [20] in Cytoscape. Properties of scale-free, small word, and modularity of the QSYQ’s PIN were also investigated.

Functional modules of the network were explored by the MCL [21] which simulates a flow on the graph by calculating successive powers of the associated adjacency matrix and the value of the inflation parameter strongly influences the number of clusters. Compared to the other algorithms, for example, RNSC [22], MCODE [23], and SPC [24], the MCL is superior with highlighting the robustness to graph alterations [25]. Based on the identified modules, GO enrichment analysis was utilized to predict possible biological roles of the modules by evaluating the involved biological processes, using the BinGO [26] plugin for Cytoscape.

3. Results and Discussion

3.1. The Analysis of the Main Active Components of QSYQ

The main active components of QSYQ are all related to the effect of tonifying Qi or promoting blood stasis.

Tanshinone IIA, cryptotanshinone, salvianolic acid A, salvianolic acid B, tanshinol, and protocatechuic aldehyde are from Salvia miltiorrhiza which is a classical traditional Chinese medicine (TCM) which can promote blood circulation and remove blood stasis with 1000 years of clinical application [27]. It has been demonstrated that Salvia miltiorrhiza can reduce the area of cerebral infarct of ischemia-reperfusion injury rats which results from blood stasis [28]. The chemical components of Salvia miltiorrhiza are divided into water-soluble and liposoluble components. Among the liposoluble components, tanshinone IIA [29] has been reported to improve blood stasis syndrome of patients with coronary heart diseases by inhibiting the circulating inflammatory markers (including IL-6, TNF α, VCAM-1, CD40, sCD40L, MCP-1, and MMP9). Cryptotanshinone [30] has good pharmacological effects on atherosclerosis, while atherosclerosis is one of the diseases resulting from blood stasis. Salvianolic acids, as the main effective components of water-soluble components including salvianolic acid A, salvianolic acid B, tanshinol, and protocatechuic aldehyde, can inhibit thrombosis, thromboxane B2 formation, and platelet aggregation [31]. This indicated that the main active components from Salvia miltiorrhiza are all associated with blood stasis.

Dencichine, ginsenoside Rb1, ginsenoside Rg1, and notoginsenoside R1 are from Panax notoginseng, which is a highly-valued herb and is able to modulate vascular tone such as the activation of blood circulation, removal of blood stasis, and inhibition of platelet aggregation [32]. The main active components of Panax notoginseng include two types of bioactive molecules: one has been reported to have good hemostatic and antithrombotic effects, such as dencichine [33]. In addition, saponins, as the main blood-activating components, which include ginsenoside Rb1, ginsenoside Rg1, and notoginsenoside R1, have showed significant effectiveness on treating cardiovascular diseases [34, 35].

Butein, formononetin, isoliquiritigenin, and nerolidol are from Dalbergia odorifera. Dalbergia odorifera, as blood-activating and stasis-removing TCM, is widely used for promoting blood circulation, relieving pain, and removing blood stasis, which has the effects on antithrombosis, antiplatelet aggregation, antioxidant, antitumor, and anti-inflammation [36]. Volatile oil and flavonoid compounds are two main chemical components of Dalbergia odorifera. According to Guo et al. [37], the ethyl acetate part of Dalbergia odorifera can significantly shorten the bleeding time and clotting time of mice, and it indicated that volatile oil is the material basis of blood-activation in Dalbergia odorifera, while Nerolidol, as a main active component, accounts for 45.23~69.13% of the volatile oil. Butein, formononetin, and isoliquiritigenin, as flavonoid components, show antioxidant activity, antiplatelet aggregation, anti-inflammatory properties, and the capacity for treating cardiovascular diseases [38–40].

Calycosin, astragaloside Ι, formononetin, and astragaloside IV are from Astragalus membranaceus which is a popular Qi-tonifying herb with multiple biological functions, such as antioxidative, antihypertensive, antiaging, and immunomodulatory activities [41]. The main bioactive components including isoflavonoids and triterpene saponins are associated with effects on human health [42]. Isoflavonoids, which are considered “marker components” for the quality control of Astragalus membranaceus including calycosin and formononetin, show strong antioxidant activity, immunoregulation, anti-inflammatory properties, and the capacity for treating cardiovascular diseases [43]. Astragaloside, including astragaloside Ι and astragaloside IV, is the main effective component of astragalus polysaccharides and exerts significant effects on myocardial protection and immunity enhancement [44, 45].

3.2. Targets Information of Active Components of QSYQ

75 targets were obtained from pharmacophore virtual screening. 174 and 65 targets were, respectively, extracted from the ChEMBL and STITCH 3.1. The targets’ number of each active component is listed in Table 2, and the additional targets’ information is shown in Table S1 in Supplementary Material available online at http://dx.doi.org/10.1155/2015/497314.

Table 2: The targets’ number of each active component from QSYQ.

Active components
total

Tanshinone IIA
45

Cryptotanshinone
34

Salvianolic acid A
17

Salvianolic acid B
21

Tanshinol
8

Protocatechuic aldehyde
9

Dencichine
7

Ginsenoside Rb1
15

Ginsenoside Rg1
18

Notoginsenoside R1
11

Butein
23

Formononetin
27

Isoliquiritigenin
51

Nerolidol
3

Calycosin
8

Astragaloside Ι
7

Astragaloside IV
10

3.3. Construction of Network

PPIs information of the targets from String 9.1 with their confidence score > 0.7 was imported in Cytoscape 2.8.3, and then union calculation was carried out, followed by the removal of duplicated edges of PPIs using Advanced Network Merge [20] of Plugins. The structural information of constructed networks was listed in Table 3.

Table 3: The structural information of networks on herbs and formula.

Networks
Nodes
Edges

Salvia miltiorrhiza
604
2362

Panax notoginseng
264
963

Dalbergia odorifera
588
2379

Astragalus membranaceus
399
1294

QSYQ
993
4215

3.4. Network Analysis

3.4.1. Topological Analysis

All the topological parameters of QSYQ were calculated and they are shown in Table 4.

Table 4: The simple parameters of protein interaction network of QSYQ.

Parameters
PIN of QSYQ

Clustering coefficient
0.673

Network diameter (radius)
11 (1)

Network centralization
0.104

Shortest path
804676

Mean path length
4.455

Network heterogeneity
0.955

Notes. The network diameter is the longest distance between any pair of vertices and the radius of a graph is the minimum eccentricity of any vertex. Network centralization is a network index that measures the degree of dispersion of all node centrality scores in a network. And network heterogeneity quantifies the degree of uneven distribution of the network.

Biological networks have been proposed to have scale-free topology whose degree distribution follows a power law distribution () [61]. As shown in Figure 1(a), the degree distribution of the PIN of QSYQ followed the power law distribution and the equation is . So, the PIN of QSYQ was a scale-free network.

(a)

(b)

(c)

(d)

Figure 1: Topological properties of network. (a) The degree distribution of QSYQ network; (b) shortest path length distribution of QSYQ network; (c) average clustering coefficient of QSYQ network; (d) average clustering coefficient of random network.

Small world networks have a property that mean path length is short [62]. The shortest path length between any two proteins was calculated, and it turned out to be 4.455. As shown in Figure 1(b), network path length was mostly concentrated in 3–5 steps, which meant that most proteins were closely linked and the PIN of QSYQ was a small world network.

In graph theory, a clustering coefficient is a measure of the degree to which nodes in a graph tend to cluster together. As shown in Figures 1(c) and 1(d), compared with random network whose numbers of nodes and edges are the same as PIN of QSYQ, the clustering coefficient of PIN was higher. It meant the PIN of QSYQ was more modular. These results suggested that the network exhibited the properties of scale-free, small word, and modularity.

3.4.2. Clustering and GO Enrichment Analysis

With the MCL algorithm, 57, 24, 49, 29, and 85 modules were, respectively, identified from salvia miltiorrhiza, Panax notoginseng, Dalbergia odorifera, Astragalus membranaceus, and QSYQ. The modules of QSYQ are shown in Figure 2, and the others are shown in Figures S1–S4.

Figure 2: Modules in the PIN of QSYQ. With the MCL algorithm, 85 modules are extracted from the network.

The results of functional enrichment analysis of QSYQ using BinGO are shown in Table 5, and they show that QSYQ played a pharmacodynamics with the biological processes, such as DNA metabolic process, regulation of cAMP metabolic process, lipid metabolic process, and the regulation of blood coagulation. The results of functional enrichment analysis of salvia miltiorrhiza, Panax notoginseng, Dalbergia odorifera, and Astragalus membranaceus are shown in Tables S2–S5.

Table 5: GO biological process terms of the modules of QSYQ.

Modules
value
GO terms

1
1.97E − 14
Regulation of protein metabolic process

2
1.32E − 32
DNA metabolic process

3
2.01E − 28
Regulation of cAMP metabolic process

4
4.16E − 24
G-protein coupled receptor signaling pathway

5
5.28E − 23
DNA-dependent transcription, initiation

6
3.03E − 25
Transmembrane receptor protein tyrosine kinase signaling pathway

7
4.74E − 23
Cellular lipid metabolic process

8
6.43E − 15
Apoptotic process

9
1.80E − 18
Tricarboxylic acid cycle

10
9.96E − 27
G-protein coupled receptor signaling pathway

11
1.32E − 32
Xenobiotic metabolic process

12
2.74E − 15
Toll-like receptor signaling pathway

13
5.96E − 32
Potassium ion transport

14
1.50E − 20
Lipid metabolic process

15
2.57E − 20
Xenobiotic metabolic process

16
8.66E − 12
Positive regulation of RNA metabolic process

17
2.78E − 27
Regulation of blood coagulation

18
3.33E − 14
Inflammatory response

19
1.94E − 16
Immune response-activating signal transduction

20
3.65E − 19
Apoptotic process

21
1.04E − 11
Regulation of blood coagulation

22
2.64E − 16
Nucleotide metabolic process

23
4.14E − 13
Transmembrane receptor protein tyrosine kinase signaling pathway

24
3.85E − 12
Interferon-gamma-mediated signaling pathway

25
5.59E − 07
Regulation of cellular protein metabolic process

26
4.34E − 06
RNA processing

27
1.19E − 22
Cell cycle phase transition

28
2.71E − 06
Regulation of RNA metabolic process

29
9.03E − 14
Regulation of systemic arterial blood pressure by renin-angiotensin

30
5.05E − 18
Cardiac muscle contraction

31
2.53E − 06
Regulation of RNA splicing

32
1.13E − 09
Carbohydrate metabolic process

33
9.21E − 14
Insulin receptor signaling pathway

34
3.77E − 08
Lipid catabolic process

35
2.50E − 22
Cellular amino acid catabolic process

36
1.93E − 17
Regulation of cell cycle

37
2.42E − 14
Lipid metabolic process

38
3.52E − 17
mRNA metabolic process

39
2.31E − 06
Execution phase of apoptosis

40
1.16E − 15
Toll-like receptor signaling pathway

41
3.70E − 11
RNA biosynthetic process

42
3.65E − 12
Cell cycle phase

43
2.80E − 12
Cellular amino acid metabolic process

44
6.84E − 05
Regulation of cell proliferation

45
1.17E − 07
Inflammatory response

46
4.51E − 09
RNA biosynthetic process

47
3.74E − 11
Regulation of interleukin-1 secretion

48
1.49E − 13
DNA repair

49
1.67E − 18
Insulin receptor signaling pathway

50
3.18E − 06
Negative regulation of RNA metabolic process

51
1.24E − 04
Transport

52
8.34E − 14
Regulation of transforming growth factor beta receptor signaling pathway

53
4.47E − 06
Xenobiotic metabolic process

54
1.16E − 08
Negative regulation of inflammatory response

55
1.37E − 07
Positive regulation of RNA metabolic process

56
5.56E − 07
Transmission of nerve impulse

57
5.79E − 12
Mitotic cell cycle

58
6.69E − 06
Negative regulation of protein metabolic process

59
2.66E − 12
Regulation of apoptotic signaling pathway

60
2.00E − 09
Regulation of apoptotic signaling pathway

61
2.47E − 11
NIK/NF-kappaB cascade

62
4.69E − 15
Vascular endothelial growth factor receptor signaling pathway

63
1.98E − 06
DNA metabolic process

64
2.04E − 04
Chromatin organization

65
1.11E − 05
DNA packaging

66
9.39E − 07
Inflammatory response

67
3.22E − 09
Xenobiotic metabolic process

68
8.09E − 12
Insulin receptor signaling pathway

69
8.52E − 04
DNA-dependent transcription, initiation

70
1.56E − 05
DNA replication

71
5.51E − 07
Toll-like receptor signaling pathway

72
2.16E − 08
Copper ion import

73
2.90E − 08
Transmission of nerve impulse

74
4.63E − 07
Regulation of blood coagulation

76
6.90E − 04
Alanine catabolic process

77
6.76E − 05
RNA biosynthetic process

78
7.60E − 04
Lipid metabolic process

79
2.86E − 05
TRIF-dependent toll-like receptor signaling pathway

81
4.18E − 05
Regulation of type I interferon production

83
5.21E − 03
Cellular lipid metabolic process

84
1.27E − 11
Neural crest cell migration

85
2.10E − 03
AMP catabolic process

Notes. value is the probability of obtaining the observed effect under the null hypothesis; a very small value indicates that the observed effect is very unlikely to have arisen purely by chance.

(1) Modules Related to Qi Deficiency. In TCM, Qi refers to the energy which flows within our body, to support a variety of biological functions such as movement, digesting food, and fight against diseases [43]. Qi deficiency is reflected in the lack of energy. Therefore, the regulation of energy metabolism would improve the Qi deficiency. As shown in Table 5, QSYQ participated in the amino acid metabolic process, carbohydrate metabolic process, lipid metabolic process, and cAMP metabolic process which are related to the energy metabolism and also have been demonstrated to play critical roles in cardiovascular diseases [63–66]. Among them, amino acid metabolism, carbohydrate metabolism, and lipid metabolism are the main energy source of the body.

Amino acid metabolism (module 35) contained proteins such as GLUD2, GLUD1, and GLS. Glutamate dehydrogenase (GLUD) is an enzyme central to the glutamate and energy metabolism of the cell [67]. GLUD activity is raised in order to increase the amount of α-ketoglutarate produced, which can be used to provide energy by being used in the citric acid cycle to ultimately produce ATP. GLUD2 and GLUD1 are the GLUD’s isozymes that differ in amino acid sequence but catalyze the same chemical reaction. Glutaminase (GLS) is a multifunctional enzyme involved in energy metabolism [68]. And GLS is the GLS2’s isozyme, which regulates cellular energy metabolism by increasing production of glutamate and alpha-ketoglutarate and in turn results in enhanced mitochondrial respiration and ATP generation [69]. This shows that proteins in amino acid metabolism are all involved in the energy metabolism, and QSYQ can improve the Qi deficiency by regulating the amino acid metabolism.

Carbohydrate metabolism is the basis of the body to produce energy. Carbohydrate metabolism (module 32) contained proteins such as GALK1, SORD, and DCXR. Galactokinase 1 (GALK1) is an enzyme (phosphotransferase) that facilitates the phosphorylation of α-D-galactose to galactose 1-phosphate at the expense of one molecule of ATP. Sorbitol dehydrogenase (SORD) is an enzyme in carbohydrate metabolism converting sorbitol, the sugar alcohol form of glucose, into fructose [70]. Dicarbonyl/L-xylulose reductase (DCXR) is involved in carbohydrate metabolism and glucose metabolism which is a highly conserved and phylogenetically widespread enzyme converting L-xylulose into xylitol [71]. This shows that proteins in carbohydrate metabolism make contribution to provide energy for the body by participating in carbohydrate metabolism.

Lipid metabolism (module 83) contained proteins such as ACOT8, AACS. Acyl-coenzyme A thioesterase 8 (ACOT8) is a peroxisomal thioesterase involved more in the oxidation of fatty acids which are in order to generate acetyl-CoA, the entry molecule for the citric acid cycle, the main energy supply of animals [72]. Acetoacetyl-CoA synthetase (AACS) can directly activate ketone bodies for the synthesis of physiologically important lipidic substances such as cholesterol and fatty acid [73]. So, AACS can provide basic substances for energy metabolism. This shows that proteins in lipid metabolism are all related to energy metabolism, and QSYQ can improve the Qi deficiency by regulating the lipid metabolism.

cAMP metabolism (module 3) contained proteins such as GCG, ADCY7, and ADCYAP1. Glucagon (GCG) is a peptide hormone of cAMP metabolic process, which generally elevates the concentration of glucose in the blood by promoting gluconeogenesis and glycogenolysis [74]. Adenylate cyclase type 7 (ADCY7) is a membrane-bound adenylate cyclase that catalyses the formation of cyclic AMP from ATP [75]. ADCYAP1 is also known as pituitary adenylate cyclase-activating polypeptide (PACAP), which stimulates adenylate cyclase and subsequently increases the cAMP level and plays crucial roles in energy metabolism, including lipid metabolism [76].

This indicated that the QSYQ reinforced Qi efficacy by the regulation of the cAMP metabolism, amino acid metabolism, carbohydrate metabolism, and lipid metabolism. And Qi deficiency may be associate with the modules including amino acid metabolism, carbohydrate metabolism, lipid metabolism, and the cAMP metabolism.

(2) Modules Related to Blood Stasis. Blood stasis is caused by disturbance of blood circulation and is reflected in microcirculation relating to vessel and cell function, such as blood viscosity and blood cell adhesion [77]. As shown in Table 5, QSYQ took part in the regulation of blood coagulation and cardiac muscle contraction which can promote blood circulation.

The regulation of blood coagulation (module 17) contained proteins such as GGCX, F2, and SERPIND1. Gamma-glutamyl carboxylase (GGCX) catalyzes the posttranslational modification of vitamin K-dependent proteins which are involved in coagulation [78]. F2 is also known as thrombin (IIa) acts as a serine protease that converts soluble fibrinogen into insoluble strands of fibrin and activation of thrombin is crucial in physiological and pathological coagulation [79]. SERPIND1, known as heparin cofactor II, is a coagulation factor which rapidly inhibits thrombin in the presence of dermatan sulfate or heparin. SERPIND1 deficiency can lead to increased thrombin generation and a hypercoagulable state [80]. This shows that proteins in this module are all involved in the blood coagulation, and QSYQ can improve the blood stasis by the regulation of blood coagulation.

The cardiac muscle contraction (module 30) contained proteins such as MYL2, TNNC1, and TNNI3. MYL2 is also known as myosin regulatory light chain 2, ventricular/cardiac muscle isoform (MLC-2v) which plays a key role in the regulation of cardiac muscle contraction, through its interactions with myosin [81]. TNNC1 is also known as troponin C which is a protein that resides in the troponin complex on actin thin filaments of striated muscle (cardiac) and is responsible for binding calcium to activate muscle contraction [82]. Troponin I (TNNI3) has been shown to interact with TNNC1 [83] and has been reported to have a special role in the control of cardiac contractility [84]. This shows that proteins in this module are all participated in the cardiac muscle contraction. The mechanism of QSYQ has been reported to be related to improve myocardial function [7]. So, QSYQ can promote blood circulation and hence can improve the blood stasis by regulating the cardiac muscle contraction.

This indicated that QSYQ improved the blood stasis through the regulation of blood coagulation and cardiac muscle contraction. And blood stasis may be associated with the modules including the regulation of blood coagulation and cardiac muscle contraction.

3.4.3. The Synergetic Effects of QSYQ

Synergetic effects occur when the efficacy of herbs are combined. The scientific interpretation of these properties is a benefit to the explanation of the compatibility rule and it is further beneficial to the action mechanism of formulae. Synergy refers to the efficacy of combinations of herbs that is greater than the summed responses of each individual herb. As shown in Figure 3, Salvia miltiorrhiza, Panax notoginseng, Dalbergia odorifera, and Astragalus membranaceus all participate in the energy metabolism process, including cAMP metabolic process, carbohydrate metabolic process, and lipid metabolic process, and they hence have the synergetic effect on enhancing the Qi efficacy of QSYQ. The regulation of blood coagulation is involved by four herbs which reinforce the efficacy of promoting the blood circulation of QSYQ. This indicated that the synergy of formula can be illustrated based on the functional modules.

Figure 3: The schematic diagram of biological processes QSYQ and its herbs involved in. The hexagon represents formula. The triangle represents the herbs of QSYQ. The circle represents the biological processes.

4. Conclusion

In this paper, the PIN of QSYQ exhibited the properties of scale-free, small world, and modularity based on the analysis of topological parameters. A module-based network analysis approach was proposed to expound the mechanism of QSYQ. Qi-tonifying efficacy of QSYQ may be partly attributed to the regulation of amino acid metabolic process, carbohydrate metabolic process, lipid metabolic process, and the cAMP metabolic process, while QSYQ improves the blood stasis through the regulation of blood coagulation and cardiac muscle contraction. A systematic exploration of mechanism of QSYQ based on module-network analysis may bring out the best between research on drug molecules and TCM phenotypic information, so as to facilitate the therapy for the disease.

Further experiments are needed to confirm the conclusions. However, despite the lack of validation of wet experiments, this study provides an efficient way to understand the mechanisms of QSYQ faster and better considering the complexity of TCM analogous formulae. What is more, the scientific intension of “synergy” of TCM can be also illustrated based on the functional modules at the molecular level.

Conflict of Interests

The authors declared that there is no conflict of interests.

Acknowledgments

This research is financially supported by the National Natural Science Foundation of China (nos. 81430094 and 81173522) and the National Key Technology R&D Program (2008BAI51B01) in Beijing University of Chinese Medicine.

References

Y. B. Ge and J. Y. Mao, “Study on the distribution characteristics of the TCM syndromes of 7512 coronary artery disease patients,” Shandong Journal of Traditional Chinese Medicine, vol. 4, pp. 227–229, 2011.

Y. Z. Hou, S. Wang, Z. Q. Zhao et al., “Clinical assessment of complementary treatment with Qishen Yiqi dripping pills on ischemic heart failure: Study protocol for a randomized, double-blind, multicenter, placebo-controlled trial (CACT-IHF),” Trials, vol. 14, article 138, 2013.

D. X. Xie, X. Z. Huang, and B. Y. Mao, “Mechanisms of Qishen Yiqi dropping pills on ventricular remodeling after myocardial infarction of rat model,” Chinese Journal of Experimental Traditional Medical Formulae, no. 6, pp. 180–183, 2011.

Y. Z. Song, L. P. Guo, H. Shang, and J. Wang, “Effects of Qishenyiqi dripping pills on lipid metabolism in experimental hypercholesterolemia rabbits,” Jilin Journal of Traditional Chinese Medicine, vol. 1, no. 31, pp. 71–73, 2011.

Y. Wang, J. Wang, L. Guo, and X. Gao, “Antiplatelet effects of Qishen Yiqi Dropping Pill in platelets aggregation in hyperlipidemic rabbits,” Evidence-Based Complementary and Alternative Medicine, vol. 2012, Article ID 205451, 5 pages, 2012.

J.-Y. Tong, Y.-J. Xu, Y.-P. Bian, X.-B. Shen, L. Yan, and X.-Y. Zhu, “Effect and mechanism of Qishen Yiqi Pills on adriamycin-induced cardiomyopathy in mice,” Chinese Journal of Natural Medicines, vol. 11, no. 5, pp. 514–518, 2013.

Z.-T. Cui, W.-L. Wei, M. Liu, and W.-J. Wang, “Effect of pretreatment with Qishen Yiqi dropping pills on right cardiac function of patients undergoing valve replacement,” Zhongguo Zhongyao Zazhi, vol. 39, no. 5, pp. 916–919, 2014.

A.-L. Barabási, “Scale-free networks: a decade and beyond,” Science, vol. 325, no. 5939, pp. 412–413, 2009.

S. C. Janga and A. Tzakos, “Structure and organization of drug-target networks: insights from genomic approaches for drug discovery,” Molecular BioSystems, vol. 5, no. 12, pp. 1536–1548, 2009.

D. Pal, “On gene ontology and function annotation,” Bioinformation, vol. 1, no. 3, pp. 97–98, 2006.

A. P. Bento, A. Gaulton, A. Hersey et al., “The ChEMBL bioactivity database: an update,” Nucleic Acids Research, vol. 42, no. 1, pp. D1083–D1090, 2014.

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.

X. Wang, Y. Xiang, Z. Ren, Y. Zhang, and Y. Qiao, “Rational questing for inhibitors of endothelin converting enzyme-1 from Salvia miltiorrhiza by combining ligand-and structure-based virtual screening,” Canadian Journal of Chemistry, vol. 91, no. 6, pp. 448–456, 2013.

X. Wang, Z. Ren, Y. He, Y. Xiang, Y. Zhang, and Y. Qiao, “A combination of pharmacophore modeling, molecular docking and virtual screening for iNOS inhibitors from Chinese herbs,” Bio-Medical Materials and Engineering, vol. 24, no. 1, pp. 1315–1322, 2014.

A. Franceschini, D. Szklarczyk, S. Frankild et al., “STRING v9.1: protein-protein interaction networks, with increased coverage and integration,” Nucleic Acids Research, vol. 41, no. 1, pp. D808–D815, 2013.

G. Su, J. H. Morris, B. Demchak, and G. D. Bader, “Biological network exploration with cytoscape 3,” Current Protocols in Bioinformatics, vol. 47, pp. 8.13.1–8.13.24, 2014.

A.-L. Barabási and Z. N. Oltvai, “Network biology: understanding the cell's functional organization,” Nature Reviews Genetics, vol. 5, no. 2, pp. 101–113, 2004.

E. Almaas, “Biological impacts and context of network theory,” The Journal of Experimental Biology, vol. 210, part 9, pp. 1548–1558, 2007.

X. Zhu, M. Gerstein, and M. Snyder, “Getting connected: analysis and principles of biological networks,” Genes and Development, vol. 21, no. 9, pp. 1010–1024, 2007.

Y. Assenov, F. Ramírez, S.-E. S.-E. Schelhorn, T. Lengauer, and M. Albrecht, “Computing topological parameters of biological networks,” Bioinformatics, vol. 24, no. 2, pp. 282–284, 2008.

A. J. Enright, S. Van Dongen, and C. A. Ouzounis, “An efficient algorithm for large-scale detection of protein families,” Nucleic Acids Research, vol. 30, no. 7, pp. 1575–1584, 2002.

A. D. King, N. Pržulj, and I. Jurisica, “Protein complex prediction via cost-based clustering,” Bioinformatics, vol. 20, no. 17, pp. 3013–3020, 2004.

G. D. Bader and C. W. V. Hogue, “An automated method for finding molecular complexes in large protein interaction networks,” BMC Bioinformatics, vol. 4, no. 1, article 2, 2003.

M. Blatt, S. Wiseman, and E. Domany, “Superparamagnetic clustering of data,” Physical Review Letters, vol. 76, no. 18, pp. 3251–3254, 1996.

S. Brohée and J. van Helden, “Evaluation of clustering algorithms for protein-protein interaction networks,” BMC Bioinformatics, vol. 7, article 488, 2006.

S. Maere, K. Heymans, and M. Kuiper, “BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks,” Bioinformatics, vol. 21, no. 16, pp. 3448–3449, 2005.

X. Chen, J. Guo, J. Bao, J. Lu, and Y. Wang, “The anticancer properties of Salvia miltiorrhiza Bunge (Danshen): a systematic review,” Medicinal Research Reviews, vol. 34, no. 4, pp. 768–794, 2014.

C.-J. Lao, J.-G. Lin, J.-S. Kuo et al., “Effect of salvia miltiorrhiza bunge on cerebral infarct in ischemia-reperfusion injured rats,” The American Journal of Chinese Medicine, vol. 31, no. 2, pp. 191–200, 2003.

Q. Shang, H. Wang, S. Li, and H. Xu, “The effect of sodium tanshinone IIA sulfate and simvastatin on elevated serum levels of inflammatory markers in patients with coronary heart disease: a study protocol for a randomized controlled trial,” Evidence-based Complementary and Alternative Medicine, vol. 2013, Article ID 756519, 8 pages, 2013.

Z. Liu, S. Xu, X. Huang et al., “Cryptotanshinone, an orally bioactive herbal compound from Danshen, attenuates atherosclerosis in apolipoprotein E-deficient mice: role of LOX-1,” British Journal of Pharmacology, 2015.

R.-W. Jiang, K.-M. Lau, P.-M. Hon, T. C. W. Mak, K.-S. Woo, and K.-P. Fung, “Chemistry and biological activities of caffeic acid derivatives from Salvia miltiorrhiza,” Current Medicinal Chemistry, vol. 12, no. 2, pp. 237–246, 2005.

L. Y. Ma and P. G. Xiao, “Effects of Panax notoginseng saponins on platelet aggregation in rats with middle cerebral artery occlusion or in vitro and on lipid fluidity of platelet membrane,” Phytotherapy Research, vol. 12, no. 2, pp. 138–140, 1998.

J. V. Formica and W. Regelson, “Review of the biology of Quercetin and related bioflavonoids,” Food and Chemical Toxicology, vol. 33, no. 12, pp. 1061–1080, 1995.

X. J. Zhang, C. He, K. Tian, P. Li, H. Su, and J. Wan, “Ginsenoside Rb1 attenuates angiotensin II-induced abdominal aortic aneurysm through inactivation of the JNK and p38 signaling pathways,” Vascular Pharmacology, 2015.

C. Jia, M. Xiong, P. Wang et al., “Notoginsenoside R1 attenuates atherosclerotic lesions in ApoE deficient mouse model,” PLoS ONE, vol. 9, no. 6, Article ID e99849, 2014.

Z.-H. Yang, C. Mei, X.-H. He, and X.-B. Sun, “Advance in studies on chemical constitutions, pharmacological mechanism and pharmacokinetic profile of dalbergiae odoriferae lignum,” Zhongguo Zhong Yao Za Zhi, vol. 38, no. 11, pp. 1679–1683, 2013.

L. B. Guo, L. R. Huang, L. H. Zhao, and C. P. Tang, “Pharmacological screening of effective parts on Xingqi Zhitong and Huoxue Zhixue in Dalbergia odorifera,” Journal of Chinese Medicinal Materials, vol. 30, no. 6, pp. 696–698, 2007.

Z.-J. Cheng, S.-C. Kuo, S.-C. Chan, F.-N. Ko, and C.-M. Teng, “Antioxidant properties of butein isolated from Dalbergia odorifera,” Biochimica et Biophysica Acta, vol. 1392, no. 2-3, pp. 291–299, 1998.

C. Zhan and J. Yang, “Protective effects of isoliquiritigenin in transient middle cerebral artery occlusion-induced focal cerebral ischemia in rats,” Pharmacological Research, vol. 53, no. 3, pp. 303–309, 2006.

J. Sung and J. Lee, “Anti-inflammatory activity of butein and luteolin through suppression of NFκB activation and induction of heme oxygenase-1,” Journal of Medicinal Food, vol. 18, no. 5, pp. 557–564, 2015.

T.-S. Yeh, H.-L. Chuang, W.-C. Huang, Y.-M. Chen, C.-C. Huang, and M.-C. Hsu, “Astragalus membranaceus improves exercise performance and ameliorates exercise-induced fatigue in trained mice,” Molecules, vol. 19, no. 3, pp. 2793–2807, 2014.

L.-Z. Lin, X.-G. He, M. Lindenmaier et al., “Liquid chromatography-electrospray ionization mass spectrometry study of the flavonoids of the roots of Astragalus mongholicus and A. membranaceus,” Journal of Chromatography A, vol. 876, no. 1-2, pp. 87–95, 2000.

W. Zhou and Y. Wang, “A network-based analysis of the types of coronary artery disease from traditional Chinese medicine perspective: potential for therapeutics and drug discovery,” Journal of Ethnopharmacology, vol. 151, no. 1, pp. 66–77, 2014.

S. Ren, H. Zhang, Y. Mu, M. Sun, and P. Liu, “Pharmacological effects of Astragaloside IV: a literature review,” Journal of Traditional Chinese Medicine, vol. 33, no. 3, pp. 413–416, 2013.

M. Lu, F. Tang, J. Zhang et al., “Astragaloside IV attenuates injury caused by myocardial ischemia/reperfusion in rats via regulation of toll-like receptor 4/nuclear factor-kappaB signaling pathway,” Phytotherapy Research, vol. 29, no. 4, pp. 599–606, 2015.

G. Ye, C.-S. Wang, Y.-Y. Li, H. Ren, and D.-A. Guo, “Simultaneous determination and pharmacokinetic studies on (3,4-dihydroxyphenyl)-lactic acid and protocatechuic aldehyde in rat serum after oral administration of Radix Salviae miltiorrhizae extract,” Journal of Chromatographic Science, vol. 41, no. 6, pp. 327–330, 2003.

T. Lu, J. Yang, X. Gao et al., “Plasma and urinary tanshinol from Salvia miltiorrhiza (Danshen) can be used as pharmacokinetic markers for cardiotonic pills, a cardiovascular herbal medicine,” Drug Metabolism and Disposition, vol. 36, no. 8, pp. 1578–1586, 2008.

R. Mahesh, H. W. Jung, G. W. Kim, Y. S. Kim, and Y.-K. Park, “Cryptotanshinone from Salviae miltiorrhizae radix inhibits sodium-nitroprusside-induced apoptosis in neuro-2 cells,” Phytotherapy Research, vol. 26, no. 8, pp. 1211–1219, 2012.

S.-C. Chiu, S.-Y. Huang, S.-F. Chang et al., “Potential therapeutic roles of tanshinone IIA in human bladder cancer cells,” International Journal of Molecular Sciences, vol. 15, no. 9, pp. 15622–15637, 2014.

T. Zhang, J. Xu, D. Li et al., “Salvianolic acid A, a matrix metalloproteinase-9 inhibitor of Salvia miltiorrhiza, attenuates aortic aneurysm formation in apolipoprotein E-deficient mice,” Phytomedicine, vol. 21, no. 10, pp. 1137–1145, 2014.

J. Fu, H.-B. Fan, Z. Guo et al., “Salvianolic acid B attenuates spinal cord ischemia-reperfusion-induced neuronal injury and oxidative stress by activating the extracellular signal-regulated kinase pathway in rats,” Journal of Surgical Research, vol. 188, no. 1, pp. 222–230, 2014.

Y. Wei, P. Li, H. Fan et al., “Metabolism study of notoginsenoside R₁, ginsenoside Rg₁ and ginsenoside Rb₁ of radix panax notoginseng in zebrafish,” Molecules, vol. 16, no. 8, pp. 6621–6633, 2011.

J. Huang, L. Ding, D. Shi et al., “Transient receptor potential vanilloid-1 participates in the inhibitory effect of ginsenoside Rg1 on capsaicin-induced interleukin-8 and prostaglandin E₂ production in HaCaT cells,” Journal of Pharmacy and Pharmacology, vol. 64, no. 2, pp. 252–258, 2012.

C.-F. Qiao, X.-M. Liu, X.-M. Cui et al., “High-performance anion-exchange chromatography coupled with diode array detection for the determination of dencichine in Panax notoginseng and related species,” Journal of Separation Science, vol. 36, no. 15, pp. 2401–2406, 2013.

Z.-J. Cheng, S.-C. Kuo, S.-C. Chan, F.-N. Ko, and C.-M. Teng, “Antioxidant properties of butein isolated from Dalbergia odorifera,” Biochimica et Biophysica Acta, vol. 1392, no. 2-3, pp. 291–299, 1998.

S. H. Lee, J. Y. Kim, G. S. Seo, Y.-C. Kim, and D. H. Sohn, “Isoliquiritigenin, from Dalbergia odorifera, up-regulates anti-inflammatory heme oxygenase-1 expression in RAW264.7 macrophages,” Inflammation Research, vol. 58, no. 5, pp. 257–262, 2009.

L. Xu, H. Shi, T. Liang et al., “Selective separation of flavonoid glycosides in Dalbergia odorifera by matrix solid-phase dispersion using titania,” Journal of Separation Science, vol. 34, no. 11, pp. 1347–1354, 2011.

M. P. N. Silva, G. L. S. Oliveira, R. B. F. De Carvalho et al., “Antischistosomal activity of the terpene nerolidol,” Molecules, vol. 19, no. 3, pp. 3793–3803, 2014.

H.-J. Kwon, J. Hwang, S.-K. Lee, and Y.-D. Park, “Astragaloside content in the periderm, cortex, and xylem of Astragalus membranaceus root,” Journal of Natural Medicines, vol. 67, no. 4, pp. 850–855, 2013.

W. Li, Y. N. Sun, X. T. Yan et al., “Flavonoids from Astragalus membranaceus and their inhibitory effects on LPS-stimulated pro-inflammatory cytokine production in bone marrow-derived dendritic cells,” Archives of Pharmacal Research, vol. 37, no. 2, pp. 186–192, 2014.

A.-L. Barabási and R. Albert, “Emergence of scaling in random networks,” Science, vol. 286, no. 5439, pp. 509–512, 1999.

S. H. Strogatz, “Exploring complex networks,” Nature, vol. 410, no. 6825, pp. 268–276, 2001.

Z. Wang and F. Li, “Effects of ilexonin A on cAMP metabolism in platelets,” Chinese Medical Sciences Journal, vol. 8, no. 4, pp. 215–217, 1993.

J. M. Gaziano, H. D. Sesso, J. L. Breslow, C. H. Hennekens, and J. E. Buring, “Relation between systemic hypertension and blood lipids on the risk of myocardial infarction,” American Journal of Cardiology, vol. 84, no. 7, pp. 768–773, 1999.

J. Mann, “Dietary carbohydrate: relationship to cardiovascular disease and disorders of carbohydrate metabolism,” European Journal of Clinical Nutrition, vol. 61, supplement 1, pp. S100–S111, 2007.

K. U. Yun, C. S. Ryu, J. M. Oh et al., “Plasma homocysteine level and hepatic sulfur amino acid metabolism in mice fed a high-fat diet,” European Journal of Nutrition, vol. 52, no. 1, pp. 127–134, 2013.

L. Rosso, A. C. Marques, A. S. Reichert, and H. Kaessmann, “Mitochondrial targeting adaptation of the hominoid-specific glutamate dehydrogenase driven by positive darwinian selection,” PLoS Genetics, vol. 4, no. 8, Article ID e1000150, 2008.

N. Bae, Y. Wang, L. Li, S. Rayport, and G. Lubec, “Network of brain protein level changes in glutaminase deficient fetal mice,” Journal of Proteomics, vol. 80, pp. 236–249, 2013.

W. Hu, C. Zhang, R. Wu, Y. Sun, A. Levine, and Z. Feng, “Glutaminase 2, a novel p53 target gene regulating energy metabolism and antioxidant function,” Proceedings of the National Academy of Sciences of the United States of America, vol. 107, no. 16, pp. 7455–7460, 2010.

O. El-Kabbani, C. Darmanin, and R. P.-T. Chung, “Sorbitol dehydrogenase: structure, function and ligand design,” Current Medicinal Chemistry, vol. 11, no. 4, pp. 465–476, 2004.

S.-K. Lee, L. T. Son, H.-J. Choi, and J. Ahnn, “Dicarbonyl/l-xylulose reductase (DCXR): the multifunctional pentosuria enzyme,” International Journal of Biochemistry and Cell Biology, vol. 45, no. 11, pp. 2563–2567, 2013.

M. C. Hunt, M. I. Siponen, and S. E. H. Alexson, “The emerging role of acyl-CoA thioesterases and acyltransferases in regulating peroxisomal lipid metabolism,” Biochimica et Biophysica Acta, vol. 1822, no. 9, pp. 1397–1410, 2012.

S. Hasegawa, Y. Ikeda, M. Yamasaki, and T. Fukui, “The role of acetoacetyl-CoA synthetase, a ketone body-utilizing enzyme, in 3T3-L1 adipocyte differentiation,” Biological and Pharmaceutical Bulletin, vol. 35, no. 11, pp. 1980–1985, 2012.

B. Ahrén, “Glucagon—early breakthroughs and recent discoveries,” Peptides, vol. 67, pp. 74–81, 2015.

K. Hellevuo, R. Berry, J. M. Sikela, and B. Tabakoff, “Localization of the gene for a novel human adenylyl cyclase (ADCY7) to chromosome 16,” Human Genetics, vol. 95, no. 2, pp. 197–200, 1995.

S. Tomimoto, T. Ojika, N. Shintani et al., “Markedly reduced white adipose tissue and increased insulin sensitivity in Adcyap1-deficient mice,” Journal of Pharmacological Sciences, vol. 107, no. 1, pp. 41–48, 2008.

J. Hou, J. Wang, C. Lin et al., “Circulating microRNA profiles differ between Qi-stagnation and Qi-deficiency in coronary heart disease patients with blood stasis syndrome,” Evidence-Based Complementary and Alternative Medicine, vol. 2014, Article ID 926962, 9 pages, 2014.

C. R. King, E. Deych, P. Milligan et al., “Gamma-glutamyl carboxylase and its influence on warfarin dose,” Thrombosis and Haemostasis, vol. 104, no. 4, pp. 750–754, 2010.

K. Guria and G. T. Guria, “Spatial aspects of blood coagulation: two decades of research on the self-sustained traveling wave of thrombin,” Thrombosis Research, vol. 135, no. 3, pp. 423–433, 2015.

K. A. Parker and D. M. Tollefsen, “The protease specificity of heparin cofactor II. Inhibition of thrombin generated during coagulation,” The Journal of Biological Chemistry, vol. 260, no. 6, pp. 3501–3505, 1985.

F. Sheikh, R. C. Lyon, and J. Chen, “Getting the skinny on thick filament regulation in cardiac muscle biology and disease,” Trends in Cardiovascular Medicine, vol. 24, no. 4, pp. 133–141, 2014.

P. J. Townsend, M. H. Yacoub, and P. J. R. Barton, “Assignment of the human cardiac/slow skeletal muscle troponin C gene (TNNC1) between D3S3118 and GCT4B10 on the short arm of chromosome 3 by somatic cell hybrid analysis,” Annals of Human Genetics, vol. 61, no. 4, pp. 375–377, 1997.

D. G. Ward, S. M. Brewer, C. E. Gallon, Y. Gao, B. A. Levine, and I. P. Trayer, “NMR and mutagenesis studies on the phosphorylation region of human cardiac troponin I,” Biochemistry, vol. 43, no. 19, pp. 5772–5781, 2004.

R. J. Solaro, P. Rosevear, and T. Kobayashi, “The unique functions of cardiac troponin I in the control of cardiac muscle contraction and relaxation,” Biochemical and Biophysical Research Communications, vol. 369, no. 1, pp. 82–87, 2008.