Skip to main content

ORIGINAL RESEARCH article

Front. Physiol., 08 May 2018
Sec. Vascular Physiology
This article is part of the Research Topic Traditional Chinese Medicine: Traditional Medicine and Cardiovascular/Complex diseases - Volume I View all 29 articles

Prediction of Molecular Mechanisms for LianXia NingXin Formula: A Network Pharmacology Study

\r\nYang Yang&#x;Yang Yang1†Kuo Yang&#x;Kuo Yang2†Teng HaoTeng Hao3Guodong ZhuGuodong Zhu4Ruby LingRuby Ling1Xuezhong Zhou*Xuezhong Zhou2*Ping Li*Ping Li1*
  • 1The Third Affiliated Hospital, Beijing University of Chinese Medicine, Beijing, China
  • 2Beijing Key Lab of Traffic Data Analysis and Mining, School of Computer and Information Technology, Beijing Jiaotong University, Beijing, China
  • 3Department of Psychiatry, Beijing ChaoYang Hospital of Traditional Chinese Medicine, Beijing, China
  • 4Department of Cardiovascular, Beijing Chaoyang Integrative Medicine Emergency Medical Center, Beijing, China

Objectives: Network pharmacological methods were used to investigate the underlying molecular mechanisms of LianXia NingXin (LXNX) formula, a Chinese prescription, to treat coronary heart disease (CHD) and disease phenotypes (CHD related diseases and symptoms).

Methods: The different seed gene lists associated with the herbs of LXNX formula, the CHD co-morbid diseases and symptoms which were relieved by the LXNX formula (co-morbid diseases and symptoms) were curated manually from biomedical databases and published biomedical literatures. Module enrichment analysis was used to identify CHD-related disease modules in the protein–protein interaction (PPI) network which were also associated to the targets of LXNX formula (LXNX formula’s CHD modules). The molecular characteristics of LXNX formula’s CHD modules were investigated via functional enrichment analysis in terms of gene ontology and pathways. We performed shortest path analysis to explore the interactions between the drug targets of LXNX formula and CHD related disease phenotypes (e.g., co-morbid diseases and symptoms).

Results: We identified two significant CHD related disease modules (i.e., M146 and M203), which were targeted by the herbs of LXNX formula. Pathway and GO term functional analysis results indicated that G-protein coupled receptor signaling pathways (GPCR) of M146 and cellular protein metabolic process of M203 are important functional pathways for the respective module. This is further confirmed by the shortest path analysis between the drug targets of LXNX formula and the aforementioned disease modules. In addition, corticotropin releasing hormone (CRH) and natriuretic peptide precursor A (NPPA) are the only two LXNX formula target proteins with the low shortest path length (on average shorter than 3) to their respective CHD module and co-morbid disease and symptom gene groups.

Conclusion: G-protein coupled receptor signaling pathway and cellular protein metabolic process are the key LXNX formula’s pathways to treat CHD disease phenotypes, in which CRH and NPPA are the two key drug targets of LXNX formula. Further evidences from Chinese herb pharmacological databases indicate that Pinellia ternata (Banxia) has relatively strong adjustive functions on the two key targets.

Introduction

Coronary heart disease (CHD) is a cardiac disease caused by coronary artery arteriosclerosis leading to cardiac vascular stenosis or obstruction, and/or by functional change (or vasospasm) of coronary artery which leads to myocardial ischemia or necrosis (Ni et al., 2016). According to the data published on the World Health Organization website (World Health Organization, 2017), 7.4 million people die from CHD each year. One in seven people die of CHD in United State (Benjamin et al., 2018). CHD has become a serious threat to human life and health and has been a focus of medical research all over the world for a long time. More and more researchers have realized that the development of CHD is related not only to the circulatory system, but also others such as the immune (Cao et al., 2018), nervous (Alcántara and Davidson, 2014), and endocrine systems (Baxter et al., 2003). Therefore, many researchers are now investigating new treatments for CHD by looking for medical targets in systems other than the circulatory system. Around 2000 years ago, traditional Chinese medicine (TCM) had already started and documented research on CHD treatment, as noted in texts such as Synopsis of the Golden Chamber (Jin Gui Yao Lue) (Xu et al., 2017). Even today, TCM continues to play an irreplaceable role in CHD treatment. TCM treats CHD in accordance with its holistic principle, and thus, TCM may have already been treating CHD through systems aside from the circulatory system (Zhang et al., 2017).

LianXia NingXin formula (LXNX formula) was created by Professor Ping Li, a Chinese medicine expert in our research team, on the basis of the TCM theory of the heart controlling the blood vessels, nerves and spirit (Yang et al., 2016). It also embodies principles such as TCM master Zhizheng Lu view of “treating CHD via the spleen and stomach” (Song and Lu, 2011), and academician Yongyan Wang’s summation of CHD causes as “poison damaging the veins” (Zhang et al., 2013). LXNX formula targets CHD phlegm-fire syndrome in TCM through the ways to reduce the heat and phlegm in the heart. The results of our previous clinical study (Ma, 2015) showed that after 4 weeks of treatments with LXNX formula, physical activity levels, angina occurrence frequency and degree of pain for CHD patients noticeably improved, while co-morbid symptoms (e.g., palpitation, dizziness, chest pain, chest tightness, insomnia) were relieved significantly. These co-morbid symptoms seem to be related to autonomic nervous system (Cheng et al., 2017). Pre-clinical researches (Li et al., 2015a,b) indicated that LXNX formula decreased tyrosine hydroxylase (TH) mRNA, acetylcholinesterase (AchE) mRNA, nerve growth factor (NGF) mRNA, and tyrosine kinase receptor A (TrKA) mRNA in previous test subjects, mice with myocardial infarction. The results showed that LXNX formula had positive effects on heart sympathetic nerve remodeling and the autonomic nervous system by influencing inflammatory factors, such as NGF (Montoya et al., 2017), and endocrine proteins, such as AchE (Rotundo, 2017). These studies showed that LXNX formula could adjust the nervous, inflammatory, and endocrine systems; LXNX formula’s targets and pathways to treat CHD may also relate to systems other than the circulatory system. Therefore, further research into the exact pharmacological mechanism of action for LXNX formula to treat CHD, co-morbid diseases, and symptoms is helpful for discovering new drugs to improve CHD treatment.

LianXia NingXin formula is composed of nine Chinese herbs including Coptis chinensis (Huang Lian), Pinellia Tuber (Ban Xia), Orange peels (Chen Pi), Porio coco (Fu Ling), etc. These nine Chinese herbs have more than one hundred drug ingredients which can affect more than one thousand human biological proteins. It will be extremely difficult for us to identify LXNX formula’s drug targets one by one from thousands of proteins through traditional method. In recent years, it is well recognized that human body is a complex system. A drug that only regulates one protein target may not be able to affect the whole human protein network (Lin, 2012). For such reason, the network medicine (Barabási et al., 2011; Zhou et al., 2018) approach for pharmacology was rapidly developed recently. Hopkins presented and elaborated the concept of network pharmacology for the first time (Hopkins, 2008), in accordance with the network pharmacological theory, all the human biological proteins are working as a protein network in which the proteins are the nodes and the interaction relationships between the proteins are the edges. This protein network has many characteristics as a complex network. We can find the interactions between the drug and specific node or module gene group and gain more understanding about the interaction between drug and human body through these characteristics. At present, the analysis of network pharmacology is used in western medicine widely to develop new drugs, to research on the multi-target drug treatment and drug actions, etc. (Tang et al., 2012; Peng et al., 2015). The systemic and synergetic approach of network pharmacology coincides well with the principles of traditional Chinese medicine (e.g., whole body treatment and the treatment mechanism characters of Chinese herbs formulas which had multi-ingredients and influenced multi-targets and multi-pathways) (Liu and Sun, 2012; Zhang Y.Q. and Li X., 2015). Using network pharmacology to analyze herbal prescriptions could result in rapid progress for TCM pharmacology (Li S. et al., 2014; Li, 2016), accelerating translation from an experience-based system into an evidence-based one. Through the preliminary works in this research (Xu et al., 2018), eight disease modules associated to CHD pathogenicity mechanism (CHD modules) have been identified. The identification methods were as follows: Firstly, dividing human protein–protein interaction network (PPI) into 314 topological modules (named from M1 to M314) by community detection algorithm (Dittrich et al., 2008; Wang et al., 2017), which had specific biological function and disease pathogenicity mechanism (Szklarczyk et al., 2011; Menche et al., 2015); secondly, finding the CHD seed gene (Wang et al., 2017; Xu et al., 2018) group composed of the CHD genes derived from public biomedical databases and literatures; thirdly, using the module enrichment analysis between topological modules and the CHD seed gene group to identify eight CHD modules. After that, on the basis of these CHD modules, this study planned to find CHD modules, drug targets, molecule pathways and ingredients related to LXNX formula action mechanism to treat CHD and its co-morbid diseases and symptoms by network pharmacological analysis methods such as disease module analysis, the shortest path analysis, pathway enrichment analysis and gene ontology enrichment analysis (GO term analysis).

Materials and Methods

Identifying the Terms Related to CHD Co-morbid Diseases and Symptoms

Coronary heart disease co-morbid diseases and symptoms related to LXNX formula were identified on the basis of the previous pre-clinical and clinical studies on using LXNX formula to treat CHD patients. To obtain the phenotype terms related to LXNX formula, we referred to the NLM’s Medical Subject Heading (MeSH) database1 to identify and confirm official terminology of co-morbid diseases and symptoms treatable by LXNX formula.

Gene Dataset Collection of Co-morbid Diseases and Symptoms

We used the disease terms to filter the associated genes for LXNX formula treated phenotypes from several well-established phenotype-genotype association databases. In particular, genes associated with the co-morbid diseases were identified via the OMIM2 and Malacards3 databases. Gene information for co-morbid symptoms was gathered from the Malacards, HPO4 and DiseaseConnect5 databases.

Construction of Drug Protein Dataset

The drug ingredients and protein targets information associated to nine herbs which constitute the LXNX formula were curated from HIT6, TCMID7 databases, and the biomedical literatures in PubMed8 and CNKI9 databases. The data was then organized into an herbs-ingredients-targets relationship dataset for LXNX formula.

Identifying CHD Modules Related to LXNX Formula

Module enrichment analysis was used between the eight CHD modules determined from the previous preliminary experiment and the following seed gene groups: each individual herb, each co-morbid disease, each co-morbid symptom, set of all herbs, set of all co-morbid diseases, and set of all co-morbid symptoms. Fisher’s exact test was used to calculate the enrichment degree between each seed gene group and the CHD modules. And then, the significant CHD modules related to LXNX formula were identified on the basis of the positive results (P < 0.05) of these enrichment results.

Enrichment Analysis and Shortest Path Analysis of LXNX Formula

Firstly, the functional pathways of CHD module related to LXNX formula were analyzed using pathway and gene ontology (GO) enrichment analysis based on the David platform (Huang et al., 2009). On the platform, the enrichment degree (p-value) between CHD modules related to LXNX formula was calculated by Fisher’s exact test. Gene ontology term analysis was based on Gene Ontology database, and was carried out by the plugin BiNGO 2.44 in the software Cytoscape 2.8.2. Pathway enrichment analysis was based on Reactome database and was carried out by online software KOBAS2.072. P-value was calculated in these two enrichment analyses, and P < 0.05 means the enrichment degree had statistically significance and the pathway results would be important functional mechanisms of disease modules. Secondly, the drug protein targets of LXNX formula in CHD modules related to LXNX formula were identified through shortest path analysis. The first step of the shortest path analysis is to calculate the shortest path between any node pair in the PPI network; subsequent calculations are then done to find the average shortest path length of a fixed source node to given node set. The shorter the average length is, the more important this source node is in connecting the entire node set. The Dijkstra algorithm (Dijkstra, 1959; Misa and Frana, 2010), which is the most classic algorithm to calculate shortest path of networks, was used to calculate the shortest path length of PPI network. In this study, we conducted the shortest path analysis between protein targets of LXNX formula and the genes of LXNX formula’s CHD gene modules, and the genes of all co-morbid diseases or symptoms. The average shortest path length for strongly connecting nodes was defined as less than 3. The nodes which were found by the shortest path analysis were the key drug targets of LXNX formula in CHD modules of LXNX formula or all co-morbid diseases or all co-morbid symptoms. Finally, we identified the CHD modules, key drug targets and molecular pathways related to LXNX formula for CHD treatment.

Candidate Targets Prediction of LXNX Formula

Similarly, we identified the drug targets of LXNX formula in regards to gene sets of co-morbid diseases and symptoms using shortest path analysis, respectively. Based on these drug targets and the targets for LXNX formula’ s CHD modules, key candidate drugs targets of LXNX formula to treat both CHD and its co-morbid diseases and symptoms were identified. Figure 1 showed the flowchart of the pathway analysis and target identification of LXNX formula to treat CHD.

FIGURE 1
www.frontiersin.org

FIGURE 1. A flowchart of the pathway analysis and targets identification of LXNX formula. To obtain reliable results, the following steps were taken: (A) The genetic information of LXNX formula herb ingredients, co-morbid diseases and symptoms was acquired from open biological gene databases. (B) For determining the key CHD modules related to LXNX formula, module enrichment analysis was used between each CHD module determined from the previous preliminary experiment and the following seed gene groups: each individual herb, each co-morbid disease, each co-morbid symptom, set of all herbs, set of all co-morbid diseases, and set of all co-morbid symptoms. (C) The function pathways of the key CHD modules of LXNX formula were determined using pathway and GO term enrichment analysis. (D) Using shortest path analysis, the key drug targets of LXNX formula in the key modules, respectively, all co-morbid diseases and symptoms were determined. Green lines showed the process of shortest path analysis. (E) LXNX formula key pathways to treat CHD were determined based on the results of module functional pathway analysis and LXNX formula key targets in the modules. (F) LXNX formula key drug targets were determined based on the results of shortest path analysis.

Results

CHD Related Disease Phenotype MeSH Terms

The results of previous pre-clinical studies (Li J. et al., 2014; Li et al., 2015a,b; Zhang H.N. and Li P., 2015) on LXNX formula indicated that LXNX formula had good effects on heart sympathetic nerve remodeling, arrhythmia and autonomic nerve system. By comparing the symptoms before and after the patients were treated with LXNX formula, clinical research (Ma, 2015) showed that LXNX formula can relieve CHD co-morbid symptoms significantly (i.e., chest pain, shortness of breath, palpitation, thirst, bitter taste in mouth, dizziness, chest tightness, and dream-disturbed sleep). It can also significantly relieve anxiety and insomnia. Based on the searching results from MeSH database, eight MeSH headings for five co-morbid diseases and ten Mesh headings for nine co-morbid symptoms were found in total (see Table 1).

TABLE 1
www.frontiersin.org

TABLE 1. The MeSH terms and number of genes related to LXNX formula co-morbid diseases and symptoms.

CHD-Co-morbid Diseases and Symptoms Related Genes

In accordance with the searching results from gene databases, there were 3230 genes related to LXNX formula co-morbid diseases and symptoms. These genes had 4605 interactions with co-morbid diseases and symptoms. For the symptoms of chest tightness, bitter taste in mouth, and dream-disturbed sleep, no related genes were found (see Table 1 and Supplementary Tables S1, S2).

The Herb Targets of LXNX Formula

LianXia NingXin formula ingredients and their corresponding targets information had been curated from various databases and documents with the metadata like data sources, curated time and keywords are kept as well. Only one record had been kept if there were any records with the same information about Chinese herb ingredients and targets. Finally, 134 chemical ingredients were identified from the nine herbs included in LXNX formula. These chemical ingredients were noted to have influence on 1860 human biological proteins (with totally 3424 protein interactions) (see Supplementary Table S3).

Key CHD Modules Related to LXNX Formula

Module enrichment analysis between the preliminary eight CHD modules and various seed gene groups for co-morbid diseases, co-morbid symptoms, and Chinese herbs showed that M146 and M203 had statistically significant results (p < 0.05). Specifically, M146 had 12 positive results, while M203 had 8 positive results (see Table 2).

TABLE 2
www.frontiersin.org

TABLE 2. The positive result summary of module enrichment analysis.

Pathway and GO Term Enrichment Analysis

There were 46 gene pathways related to M146 and 2 gene pathways related to M203 according to pathway enrichment analysis (P < 0.05) (see Supplementary Tables S7, S9). There were 119 gene pathways related to M146 and 29 gene pathways related to M203 according to GO term enrichment analysis (see Supplementary Tables S6, S8). Based on the results of shortest path analysis between each LXNX formula’s CHD modules and LXNX formula’s protein targets, there were 30 key drug targets of LXNX formula in M146 (see Supplementary Table S10), and 3 key drug targets of LXNX formula in M203 (see Supplementary Table S11). The G-protein coupled receptor (GPCRs) signaling pathway that M146 enriches contained more key drug targets of LXNX formula in the M146 than other enriched pathways. Cellular protein metabolic process, which was one of the M203 GO enrichment, contained all key drug targets of LXNX formula in the M203. Other functional pathway enrichment results in M203 did not contain as many key drug targets. Table 3 shows the functional pathway results of GPCRs signaling pathway related to M146, and specific cellular protein metabolic process related to M203 based on the pathway enrichment analysis and GO term enrichment analysis.

TABLE 3
www.frontiersin.org

TABLE 3. The functional pathway analysis of CHD modules related to LXNX formula.

Molecular Interactions Between Herb Targets and Modules

Thirty proteins in LXNX formula were closely related to M146 genes (see Supplementary Table S10), while 3 proteins were related to M203 genes (see Supplementary Table S11), 144 proteins were related to all co-morbid symptom genes (see Supplementary Table S5), 29 proteins were related to whole co-morbid disease genes (see Supplementary Table S4). There was no single protein in LXNX formula associated with both M146 and M203. Among LXNX formula proteins, corticotropin releasing hormone (CRH) was the key drug target simultaneously in M146 and all co-morbid diseases and symptoms. In addition, natriuretic peptide precursor A (NPPA) was the key drug target simultaneously in M203 and all co-morbid diseases and symptoms (see Table 4). Besides CRH and NPPA, there were no other proteins which served as key drug targets simultaneously in LXNX formula’s CHD modules and all co-morbid diseases and symptoms.

TABLE 4
www.frontiersin.org

TABLE 4. The shortest path analysis results about CRH and NPPA.

Figure 2 showed the gene relationships between M146 and M203 in the human protein–protein interaction network, and it also reflected the key targets (those big nodes in the figure) and pathways (those bold lines between key targets) of LXNX formula to treat CHD based on the analysis results of the enriched functional pathway and the shortest path.

FIGURE 2
www.frontiersin.org

FIGURE 2. The protein interactions of M146 and M203. The red and blue nodes represent the proteins of M146 and M203, respectively. The red and blue lines represent the protein interactions of M146 and M203, respectively.

Discussion

Presently, the incompleteness of knowledge regarding human protein interaction and genes related to diseases resulted in the hindering of research into complex systemic diseases such as CHD. Recently, it has been recognized that proteins associated with diseases did not exist randomly in a scattered network, but rather were connected into functional subgroups called disease modules. Thus, disease module analysis can now be used to identify previously unknown genes and mechanisms associated with complex diseases (Vinayagam et al., 2016; Mukund and Subramaniam, 2017). For this study, protein–protein interaction (PPI) networks, the most intensely studied biology network (Zhou et al., 2014), were used to determine disease modules. By consulting open bioinformatics databases and using module enrichment analysis, disease modules related to different diseases, symptoms, and herbs were identified based on the gene groups of these different networks. Results with significant enrichment (p < 0.05) indicated that the biological mechanism of the disease module was related to the functional pathway of the disease, symptom, or herb seed gene groups. Disease module analysis results showed that the functional pathways in M146 and M203 were strongly associated with these aforementioned seed gene groups, both individually and collectively. Other CHD modules did not have such lots of positive results. Therefore, these results indicate that the treatment effect of LXNX formula for CHD and its co-morbid diseases and symptoms is mainly accomplished by modulating M146 and M203 functional pathways.

Based on the results of pathway and GO term enrichment analysis and of the shortest path analysis, the signaling pathways of G-protein coupled receptors (GPCRs) related to adrenaline, adenylate cyclase, dopamine, and calcium ion are significant in M146, and the drug targets of LXNX formula in M146 are all related to these pathways; the pathways of specific cellular protein metabolic process related to the apolipoprotein (APO) and natriuretic peptide are important in M203, and the drug targets in M203 are all related to these pathways. GPCRs signaling transduction pathway constitute a large protein family of cellular surface receptors which are responsible for various ligand responses, signal transduction from extracellular space to intracellular space, and activating cellular responses. GPCRs are largely expressed in human body, and play a key role in regulation. For example, β-adrenergic receptor (β-AR) is a typical ligand activated GPCR. Afterβ-AR is activated, adenylate cyclase is then activated, and the concentration of cyclic adenosine monophosphate (cAMP) is subsequently increased, inducing a series of biological effects. β-AR is widely distributed in the heart and thus plays an important role in regulating cardiac functions (Mahoney and Sunahara, 2016; Gurevich and Gurevich, 2017). The pathways of cellular protein metabolic process may control cell apoptosis or proliferation by adjusting the synthesis or decomposition of the proteins in cells, thus regulating the repair and necrosis of nerves and vascular smooth muscle. For example, apolipoprotein A1 (APOA1), apolipoprotein A4 (APOA4) may impact low-density lipoprotein (LDL) particle density (Gomez et al., 2010). LDL is a protein in the cell. LDL particle density is an important factor in atherosclerosis. In accordance with the module function analysis, the shortest path analysis results between LXNX formula protein targets and LXNX formula’s CHD module gene networks, respectively, and the mechanism characters of the pathway, the two pathogenic pathways, GPCRs signaling pathway and cellular protein metabolic process, both have three same characters. Firstly, they are both the important pathogenic pathways of LXNX formula’s CHD module, respectively. Secondly, they both are composed of a lot of LXNX formula’s important drug targets in, respectively, CHD modules of LXNX formula. Thirdly, they can influence CHD and CHD co-morbid diseases which can be treated by LXNX formula. Hence, GPCRs signaling pathway and specific cellular protein metabolic process both are important pathways of LXNX formula to treat CHD and its co-morbid diseases and symptoms based on the three characters of these two pathway mentioned above.

The proteins CRH and NPPA are critical drug targets for LXNX formula’s CHD module and all co-morbid diseases and symptoms; no other protein has shown such simultaneous results. CRH is synthesized in the hypothalamus and acts upon the hypothalamus-pituitary-adrenal gland axis (HPA axis). HPA axis is a complex system which includes direct effects and negative feedback mechanisms. When stressed (a physical reaction that occurs in imbalanced circumstances or when facing threats to internal stability), the hypothalamus secretes CRH (Zhou et al., 2018), which proceeds to combine with CRH receptor 1 (CRHR1) to promote the release of adreno-corticotropic-hormone (ACTH). ACTH then acts upon the adrenal gland to stimulate the release of glucocorticoid (GC), which has a negative feedback effect on the hypothalamus and pituitary glands (Asalgoo et al., 2017; Paragliola et al., 2017). Thus, HPA axis is an important part of the neuroendocrine system and is essential to controlling the effects of stress. The stimulation of chronic stress will lead to constant HPA axis excitation. It will cause disorders of the endocrine system, imbalanced physiological functions, and mental disorders (Contoreggi, 2015). The increase of GC synthesis is a key factor in HPA axis mechanism; the constant increase of GC promotes immune system activity and inflammatory response, injures vascular endothelial cells, and accelerates arteriolosclerosis (AS) (Deuchar et al., 2011; Cho et al., 2015). CRH causes the HPA axis to release GC by controlling ACTH release; therefore, CRH can affect CHD. There are three types of CRH receptors, CRH-R1, CRH-R2, and CRH-R3, which all belong to G-protein coupled receptors. After CRH is combined with CRH receptors, G protein activates adenylate cyclase, resulting in increased cAMP cell concentration and protein kinases A (PKA) activation. The CRH function that affects CHD is achieved through GPCRs signaling pathway (Inda et al., 2017). ACTH is a proopiomelanocortin (POMC) derivative (Parvin et al., 2017). CRH is the key target of LXNX formula in M146. CRHR1 and POMC are genes not only belonging to M146, but also belonging to GPCRs signaling transduction pathway based on the results of the GO term and pathway functional enrichment analysis. Therefore, LXNX formula’s mechanism of treating CHD and co-morbid diseases and symptoms could be as such: LXNX formula acts upon CRH, and then influences the HPA axis through GPCRs signaling transduction pathway related to CRHR1 and POMC.

Natriuretic peptide precursor A is the precursor to atrial natriuretic peptide (ANP) before synthesis; NPPA polymorphism may induce myocardial necrosis, hypertension, atrial fibrillation, water and sodium disequilibrium, and heart remodeling (Lynch et al., 2009; Francia et al., 2013). After ANP combines with the natriuretic peptide receptor A (NPR-A), cGMP in cell is adjusted, thus influencing myocardial contractions and hindering the proliferation of vascular smooth muscle cells (Song et al., 2015). ANP takes part in the metabolism of the endocardium macrophage inflammatory proteins, such as tumor necrosis factor-α(TNF-α), and interleukin-1 (IL-1). These inflammatory factors can lead to AS plaque occurrence and development, thus affecting the pathological development of AS diseases, including coronary heart disease (Clarke et al., 2010; Kleinbongard et al., 2010). ANP has two pathways to adjust macrophage protein metabolism. First, ANP adjusts the metabolism and synthesis of macrophage inflammatory factors directly through cyclic guanosine monophosphate (cGMP) (Kiemer et al., 2005; Howlett and Moore, 2006). Second, ANP interaction with amyloid factors, such as APOA1, APOA4, and transthyretin (TTR), adjusts macrophage protein metabolism by affecting amyloidosis deposits in coronary artery (Maioli et al., 2000; Torricelli et al., 2004a,b). NPPA, APOA1, APOA4, and TTR are the genes belonging to both M203 and the cell protein metabolic pathway which is one significant result of the M203 GO term analysis. Therefore, LXNX formula mechanism to treat CHD is related to NPPA and the pathway of macrophage protein metabolism containing APOA1, APOA4 and TTR.

In summary, CRH and NPPA have strong relationship with LXNX formula for the reasons that they are both the most important drug targets of LXNX formula in LXNX formula’s CHD modules, respectively, co-morbid diseases and symptoms, they are both relevant to the important pathways, respectively, of LXNX formula, and they both can have influence on CHD and CHD co-morbid diseases which LXNX formula can treat. Therefore, CRH and NPPA are key drug targets of LXNX formula to treat CHD and its co-morbid diseases and symptoms.

Based on the herb-ingredient-target relationship database for LXNX formula, gamma-aminobutyric acid (GABA) found in Pinellia ternata (Ban Xia) has a direct influence on CRH; adenine in Porio coco (Fu Ling), undecanal in dried orange peels (Chen Pi), guanosine and L-arginine (L-Arg) in Pinellia ternata (Ban Xia) have a strong impact on NPPA. Evidently, Pinellia ternata (Ban Xia), Porio coco (Fu Ling), and dried orange peels (Chen Pi) are important herbs in LXNX formula. GABA is an important inhibitory neurotransmitter in the brain. By combining with GABA receptors, it inhibits CRH neuron excitability to reduce the release of CRH (Kakizawa et al., 2016). L-Arg is the precursor of EDRF (Endothelium derived relaxing factors)/NO. Shortage of EDRF/NO is essential to the occurrence of cardiovascular diseases including CHD (Charpie et al., 1997). On the other hand, L-Arg decreased ANP gene expression of mice with myocardial hypertrophy caused by isoprenaline (Shen et al., 2004).

Based on the analysis mentioned above about previous documentation, GPCRs signaling transduction pathway related to CRH and cellular protein metabolic process pathway related to NPPA both can influence immune and nervous systems to treat CHD, this character is accordance with the mechanism of LXNX formula to treat CHD. Thus, the forecast results of network medicine concerning LXNX formula’s mechanism of action for treating CHD is reliable. In future research, further experimentation can be done to test the influence of LXNX on CRH and NPPA pathway by in vitro and in vivo experiments on mice with myocardial ischemia.

Conclusion

In accordance with this research, the important pharmacological mechanisms of LXNX formula to treat CHD disease phenotypes are that the ingredients such as L-arginine and gamma-aminobutyric acid, etc. in Pinellia Tuber, Poria coco, affect CRH in GPCRs signaling pathway and NPPA in specific cellular protein metabolic process.

Author Contributions

PL, XZ, YY, and KY conceived and designed the study. YY, TH, and GZ extracted the data. KY analyzed the data. PL, XZ, YY, KY, and RL wrote and revised the manuscript. All authors reviewed and approved the submitted version of the manuscript.

Funding

This work was supported by grants from youth teacher project, Beijing University of Chinese Medicine, China (2017-JYB-JS-124), National Key R&D Project (2017YFC1703506) and the Special Programs of Traditional Chinese Medicine (No. 201407001).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2018.00489/full#supplementary-material

Footnotes

  1. ^ http://www.nlm.nih.gov/mesh/MBrowser.html
  2. ^ http://www.omim.org
  3. ^ http://www.malacards.org
  4. ^ http://human-phenotype-ontology.github.io/
  5. ^ http://disease-connect.org/
  6. ^ http://lifecenter.sgst.cn/hit
  7. ^ http://www.megabionet.org/tcmid/search/
  8. ^ www.pubmed.com
  9. ^ http://www.cnki.net

References

Alcántara, C., and Davidson, K. W. (2014). Mental disorders and coronary heart disease risk: could the evidence elude us while we sleep? Circulation 129, 139–141. doi: 10.1161/CIRCULATIONAHA.113.006515

PubMed Abstract | CrossRef Full Text | Google Scholar

Asalgoo, S., Tat, M., Sahraei, H., and Pirzad Jahromi, G. (2017). The psychoactive agent crocin can regulate hypothalamic-pituitary-adrenal axis activity. Front. Neurosci. 11:668. doi: 10.3389/fnins.2017.00668

PubMed Abstract | CrossRef Full Text | Google Scholar

Barabási, A. L., Gulbahce, N., and Loscalzo, J. (2011). Network medicine: a network-based approach to human disease. Nat. Rev. 12, 56–68. doi: 10.1038/nrg2918

PubMed Abstract | CrossRef Full Text | Google Scholar

Baxter, J. D., Young, W. F. Jr., and Webb, P. (2003). Cardiovascular endocrinology: introduction. Endocr. Rev. 24, 253–260. doi: 10.1210/er.2003-0099

PubMed Abstract | CrossRef Full Text | Google Scholar

Benjamin, E. J., Virani, S. S., Callaway, C. W., Chamberlain, A. M., Chang, A. R., Cheng, S., et al. (2018). Heart disease and stroke statistics-2018 update: a report from the American heart association. Circulation 137, e67–e492. doi: 10.1161/CIR.0000000000000558

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, D., Schiattarella, G. G., Villalobos, E., Jiang, N., May, H. I., Li, T., et al. (2018). Cytosolic DNA sensing promotes macrophage transformation and governs myocardial ischemic injury. Circulation doi: 10.1161/CIRCULATIONAHA.117.031046 [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Charpie, J. R., Charpie, P. M., Goud, C., Pitt, B., and Webb, R. C. (1997). Quinapril prevents hypertension and enhanced vascular reactivity in nitroarginine-treated rats. Blood Press. 6, 117–124. doi: 10.3109/08037059709061810

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, Z. Q., Liu, X. M., and Shi, C. J. (2017). Research development of autonomic dysfunction. Asia Pac. Tradit. Med. 13, 71–73.

Cho, W., Kang, J. L., and Park, Y. M. (2015). Corticotropin-Releasing Hormone(CRH) promotes macrophage foam cell formation via reduced expression of ATP Binding Cassette Transporter-1(ABCA1). PLoS One 10:e0130587. doi: 10.1371/journal.pone.0130587

PubMed Abstract | CrossRef Full Text | Google Scholar

Clarke, M. C., Talib, S., Figg, N. L., and Bennett, M. R. (2010). Vascular smooth muscle cell apoptosis induces interleukin-1-directed inflammation effects of hyperlipidemia-mediated inhibition of phagocytosis. Circ. Res. 106, 363–372. doi: 10.1161/CIRCRESAHA.109.208389

PubMed Abstract | CrossRef Full Text | Google Scholar

Contoreggi, C. (2015). Corticotropin releasing hormone and imaging, rethinking the stress axis. Nucl. Med. Biol. 42, 323–339. doi: 10.1016/j.nucmedbio.2014.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Deuchar, G. A., McLean, D., Hadoke, P. W. F., Brownstein, D. G., Webb, D. J., Mullins, J. J., et al. (2011). 11β-hydroxysteroid dehydrogenase type 2 deficiency accelerates atherogenesis and causes proinflammatory changes in the endothelium in apoe-/- mice. Endocrinology 152, 236–246. doi: 10.1210/en.2010-0925

PubMed Abstract | CrossRef Full Text | Google Scholar

Dijkstra, E. W. (1959). A note on two problems in connection with graphs. Numer. Math. 1, 269–271. doi: 10.1007/BF01386390

CrossRef Full Text | Google Scholar

Dittrich, M. T., Klau, G. W., Rosenwald, A., Dandekar, T., and Müller, T. (2008). Identifying functional modules in protein-protein interaction networks: an integrated exact approach. Bioinformatics 24, i223–i231. doi: 10.1093/bioinformatics/btn161

PubMed Abstract | CrossRef Full Text | Google Scholar

Francia, P., Ricotta, A., Frattari, A., Stanzione, R., Modestino, A., Mercanti, F., et al. (2013). Atrial natriuretic peptide single nucleotide polymorphisms in patients with nonfamilial structural atrial fibrillation. Clin. Med. Insights Cardiol. 7, 153–159. doi: 10.4137/CMC.S12239

PubMed Abstract | CrossRef Full Text | Google Scholar

Gomez, P., Perez-Martinez, P., Marin, C., Camargo, A., Yubero-Serrano, E. M., Garcia-Rios, A., et al. (2010). APOA1 and APOA4 gene polymorphisms influence the effects of dietary fat on LDL particle size and oxidation in healthy young adults. J. Nutr. 140, 773–778. doi: 10.3945/jn.109.115964

PubMed Abstract | CrossRef Full Text | Google Scholar

Gurevich, V. V., and Gurevich, E. V. (2017). Molecular mechanisms of GPCR signaling: a structural perspective. Int. J. Mol. Sci. 18:E2519. doi: 10.3390/ijms18122519.18

CrossRef Full Text | Google Scholar

Hopkins, A. L. (2008). Network pharmacology: the next paradigm in drug discovery. Nat. Chem. Biol. 4, 682–690. doi: 10.1038/nchembio.118

PubMed Abstract | CrossRef Full Text | Google Scholar

Howlett, G. J., and Moore, K. J. (2006). Untangling the role of amyloid in atherosclerosis. Curr. Opin. Lipidol. 17, 541–547. doi: 10.1097/01.mol.0000245260.63505.4f

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, D. W., Sherman, B. T., and Lempicki, R. A. (2009). Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nat. Protoc. 4, 44–57. doi: 10.1038/nprot.2008.211

PubMed Abstract | CrossRef Full Text | Google Scholar

Inda, C., Bonfiglio, J. J., Dos Santos Claro, P. A., Senin, S. A., Armando, N. G., Deussing, J. M., et al. (2017). cAMP-dependent cell differentiation triggered by activated CRHR1 in hippocampal neuronal cells. Sci. Rep. 7:1944. doi: 10.1038/s41598-017-02021-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Kakizawa, K., Watanabe, M., Mutoh, H., Okawa, Y., Yamashita, M., Yanagawa, Y., et al. (2016). A novel GABA-mediated corticotropin-releasing hormone secretory mechanism in the median eminence. Sci. Adv. 2:e1501723. doi: 10.1126/sciadv.1501723

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiemer, A. K., Fürst, R., and Vollmar, A. M. (2005). Vasoprotective actions of the atrial natriuretic peptide. Curr. Med. Chem. Cardiovasc. Hematol. Agents 3, 11–21. doi: 10.2174/1568016052773333

CrossRef Full Text | Google Scholar

Kleinbongard, P., Heusch, G., and Schulz, R. (2010). TNFalpha in atherosclerosis, myocardial ischemia/reperfusion and heart failure. Pharmacol. Ther. 127, 295–314. doi: 10.1016/j.pharmthera.2010.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S. (2016). Exploring traditional chinese medicine by a novel therapeutic concept of network target. Chin. J. Integr. Med. 22, 647–652. doi: 10.1007/s11655-016-2499-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S., Fan, T. P., Jia, W., Lu, A., and Zhang, W. (2014). Network pharmacology in traditional chinese medicine. Evid. Based Complement. Alternat. Med. 138460. doi: 10.1155/2014/138460

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J. L., Li, B. J., Li, P., Bai, F., Ma, H., Liu, W., et al. (2014). Effects of lianxia formula granule on experimental arrhythmia of rats induced by adrenaline. JETCM 23, 2153–2154, 2158.

Li, J. L., Li, P., and Ma, H. H. (2015a). Effects of lianxia formula granule on NGF mRNA and TrKA mRNA expression level at infarcted border zone of rats after myocardial infarction. Global Tradit. Chin. Med. 8, 303–306.

Li, J. L., Li, P., and Ma, H. H. (2015b). Effects of lianxia formula granule on TH mRNA and AchE mRNA expression level in rats after myocardial infarction. Global Tradit. Chin. Med. 8, 162–165.

Lin, B. Y. (2012). Systematic Biology. Zhejiang: Zhejiang university press, 16.

Liu, Z. H., and Sun, X. B. (2012). Network pharmacology: new opportunity for the modernization of traditional Chinese medicine. Acta Pharm. Sin. 47, 696–703.

PubMed Abstract | Google Scholar

Lynch, A. I., Claas, S. A., and Arnett, D. K. (2009). A review of the role of atrial natriuretic peptide gene polymorphisms in hypertension and its sequelae. Curr. Hypertens. Rep. 11, 35–42. doi: 10.1007/s11906-009-0008-7

CrossRef Full Text | Google Scholar

Ma, N. N. (2015). The case series study of coronary heart disease treatment using eliminating phlegm-fire syndrome. Beijing Univ. Chin. Med. 12–33.

Mahoney, J. P., and Sunahara, R. K. (2016). Mechanistic insights into GPCR-G protein interactions. Curr. Opin. Struct. Biol. 41, 247–254. doi: 10.1016/j.sbi.2016.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Maioli, E., Torricelli, C., Santucci, A., and Pacini, A. (2000). Molecular assembly of endogenous and synthetic big atrial natriuretic peptide (ANP) and its amyloidogenic implications. Biochim. Biophys. Acta 1500, 31–40. doi: 10.1016/S0925-4439(99)00083-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Menche, J., Sharma, A., Kitsak, M., Ghiassian, S. D., Vidal, M., Loscalzo, J., et al. (2015). Uncovering disease-disease relationships through the incomplete interactome. Science 347:1257601. doi: 10.1126/science.1257601

PubMed Abstract | CrossRef Full Text | Google Scholar

Misa, T. J., and Frana, P. L. (2010). An interview with Edsger Wybe Dijkstra. Commun. ACM 53, 41–47. doi: 10.1145/1787234.1787249

CrossRef Full Text

Montoya, J. G., Holmes, T. H., Anderson, J. N., Maecker, H. T., Rosenberg-Hasson, Y., Valencia, I. J., et al. (2017). Cytokine signature associated with disease severity in chronic fatigue syndrome patients. Proc. Natl. Acad. Sci. U.S.A. 114, E7150–E7158. doi: 10.1073/pnas.1710519114

PubMed Abstract | CrossRef Full Text | Google Scholar

Mukund, K., and Subramaniam, S. (2017). Co-expression network approach reveals functional similarities among diseases affecting human skeletal muscle. Front. Physiol. 8:980. doi: 10.3389/fphys.2017.00980

PubMed Abstract | CrossRef Full Text | Google Scholar

Ni, W., Chen, X. Y., and Zhang, Q. (2016). Internal medicine. China Press Tradit. Chin. Med. 4:129.

Paragliola, R. M., Papi, G., Pontecorvi, A., and Corsello, S. M. (2017). Treatment with synthetic glucocorticoids and the hypothalamus-pituitary-adrenal axis. Int. J. Mol. Sci. 18:E2201. doi: 10.3390/ijms18102201

PubMed Abstract | CrossRef Full Text | Google Scholar

Parvin, R., Saito-Hakoda, A., Shimada, H., Shimizu, K., Noro, E., Iwasaki, Y., et al. (2017). Role of NeuroD1 on the negative regulation of Pomc expression by glucocorticoid. PLoS One 12:e0175435. doi: 10.1371/journal.pone.0175435

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, X. F., Wang, L., and Bian, Y. H. (2015). The application and technology related to network pharmacology research. J. Tianjin Univ. Tradit. Chin. Med. 34, 121–124.

Rotundo, R. L. (2017). Biogenesis, assembly and trafficking of acetylcholinesterase. J. Neurochem. 142(Suppl. 2), 52–58. doi: 10.1111/jnc.13982

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, Y. H., Xiao, W., and Jin, L. (2004). Research of inhibitory effects of L-arginine on pathological cardiac hypertrophy. China Med. Herald 11, 15–18.

Song, J., and Lu, Z. Z. (2011). Study on medication law of Prof. LU Zhi-zheng in treatment of coronary heart disease in terms of adjusting spleen and stomach functions. World J. Integr. Tradit. West. Med. 6, 801–803.

Google Scholar

Song, W., Wang, H., and Wu, Q. Y. (2015). Atrial natriuretic peptide in cardiovascular biology and disease(NPPA). Gene 569, 1–6. doi: 10.1016/j.gene.2015.06.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Szklarczyk, D., Franceschini, A., Kuhn, M., Simonovic, M., Roth, A., Minguez, P., et al. (2011). The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res. 39, D561–D568. doi: 10.1093/nar/gkq973

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, P. P., Bai, M., and Miao, M. S. (2012). Traditional Chinese medicine research and network pharmacology. China J. Chin. Med. 27, 1112–1115.

Torricelli, C., Capurro, E., Santucci, A., Paffetti, A., D’Ambrosio, C., Scaloni, A., et al. (2004a). Multiple plasma proteins control atrial natriuretic peptide (ANP) aggregation. J. Mol. Endocrinol. 33, 335–341.

PubMed Abstract | Google Scholar

Torricelli, C., Capurro, E., Santucci, A., Paffetti, A., D’Ambrosio, C., Scaloni, A., et al. (2004b). Small HDL form via apo A-I a complex with atrial natriuretic peptide. Biochem. Biophys. Res. Commun. 315, 16–21.

Google Scholar

Vinayagam, A., Gibson, T. E., Lee, H. J., Yilmazel, B., Roesel, C., Hu, Y., et al. (2016). Controllability analysis of the directed human protein interaction network identifies disease genes and drug targets. Proc. Natl. Acad. Sci. U.S.A. 113, 4976–4981. doi: 10.1073/pnas.1603992113

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y. W., Liu, H. L., Lin, Y. Z., Liu, G. M., Chu, H. W., Zhao, P. Y., et al. (2017). Network-based approach to identify potential targets and drugs that promote neuroprotection and neurorepair in acute ischemic stroke. Sci. Rep. 7:40137. doi: 10.1038/srep40137

PubMed Abstract | CrossRef Full Text | Google Scholar

World Health Organization (2017). Cardiovascular Diseases (CVDs). Available at: http://www.who.int/mediacentre/factsheets/fs317/en/ [Accessed Mar 6, 2018].

Xu, W. M., Hu, J. Q., and Jiang, L. J. (2017). Analysis of the original meaning of Gualou Xiebai Three Decoctions in treating thoracic obstruction in Jin Gui Yao Lue. China J. Tradit. Chin. Med. Pharm. 32, 2288–2290.

Xu, W. M., Yang, K., Jiang, L. J., Hu, J. Q., and Zhou, X. Z. (2018). Integrated modules analysis to explore the molecular mechanisms of phlegm-stasis cementation syndrome with ischemic heart disease. Front. Physiol. 9:7. doi: 10.3389/fphys.2018.00007

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, T., Zhao, M. J., and Wang, L. (2016). Connotation and modern scientific evidence for “heart controlling mental activities”. J. Beijing Univ. Tradit. Chin. Med. 39, 811–814.

PubMed Abstract

Zhang, H. N., and Li, P. (2015). Effects of Ningxin Granule on arrhythmia after an acute myocardial infarction. China J. Tradit. Chin. Med. Pharm. 30, 1194–1196.

Zhang, J., Meng, H., Zhang, Y., Zhang, X., Shao, M., Li, C., et al. (2017). The therapeutical effect of chinese medicine for the treatment of atherosclerotic coronary heart disease. Curr. Pharm. Des. 23, 5086–5096. doi: 10.2174/1381612823666170803101602

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Zhang, Y. L., Guo, R. J., Chen, Z. G., and Wang, Y. Y. (2013). From toxin damaging brain collaterals to toxin damaging collaterals. Beijing J. Tradit. Chin. Med. 32, 483–486.

Google Scholar

Zhang, Y. Q., and Li, X. (2015). The present situation of Chinese medicine modern researches using network pharmacology. Chin. J. Pharmacol. Toxicol. 29, 883–892.

Zhou, J. J., Gao, Y., Zhang, X., Kosten, T. A., and Li, D. P. (2018). Enhanced hypothalamic NMDA receptor activity contributes to hyperactivity of HPA axis in chronic stress in male rats. Endocrinology doi: 10.1210/en.2017-03176 [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, X., Lei, L., Liu, J., Halu, A., Zhang, Y., Li, B., et al. (2018). A systems approach to refine disease taxonomy by integrating phenotypic and molecular networks. EBioMedicine (in press). doi: 10.1016/j.ebiom.2018.04.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, X., Menche, J., Barabási, A. L., and Sharma, A. (2014). Human symptoms-disease network. Nat. Commun. 5:4212. doi: 10.1038/ncomms5212

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: NPPA, CRH, network pharmacology, coronary heart disease, LianXia NingXin formula, pharmacological target

Citation: Yang Y, Yang K, Hao T, Zhu G, Ling R, Zhou X and Li P (2018) Prediction of Molecular Mechanisms for LianXia NingXin Formula: A Network Pharmacology Study. Front. Physiol. 9:489. doi: 10.3389/fphys.2018.00489

Received: 04 January 2018; Accepted: 17 April 2018;
Published: 08 May 2018.

Edited by:

Jing-Yan Han, Peking University, China

Reviewed by:

Suowen Xu, University of Rochester, United States
Aiping Lu, Hong Kong Baptist University, Hong Kong

Copyright © 2018 Yang, Yang, Hao, Zhu, Ling, Zhou and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Xuezhong Zhou, xzzhou@bjtu.edu.cn Ping Li, pearll2008@126.com

These authors have contributed equally to this work.

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.