A network pharmacology-based approach to analyse potential targets of traditional herbal formulas: An example of Yu Ping Feng decoction

Article | OPEN | Published: 30 July 2018 Huali Zuo, Qianru Zhang, Shibing Su, Qilong Chen, Fengqing Yang & Yuanjia Hu Scientific Reportsvolume 8, Article number: 11418 (2018) | Download Citation Abstract Herbal formulas from traditional Chinese medicines (TCMs) have been extensively used in clinics as effective therapies, but it is still a great challenge to demonstrate the scientific basis for their therapeutic effects at the level of molecular biology. By taking a classic herbal formula (Yu Ping Feng decoction, YPF) as an example, this study developed a novel network pharmacology based method to identify its potential therapeutic targets. First, this study constructed a “targets–(pathways)–targets” (TPT) network in which targets of YPF were connected by relevant pathways; then, this network was decomposed into separate modules with strong internal connections; lastly, the propensity of each module toward different diseases was assessed by a contribution score. On the basis of a significant association between network modules and therapeutic diseases validated by chi-square test (p-value < 0.001), this study identified the network module with the strongest propensity toward therapeutic diseases of YPF. Further, the targets with the highest centrality in this module are recommended as YPF’s potential therapeutic targets. By integrating the complicated “multi-targets–multi-pathways–multi-diseases” relationship of herbal formulas, the method shows promise for identifying its potential therapeutic targets, which could contribute to the modern scientific illustration of TCMs’ traditional clinical applications. Introduction Understanding the scientific basis and therapeutic mechanisms of herbal formulas from traditional Chinese medicines (TCMs) is essential since they play a vital role in complementary and alternative therapies in clinical practice to cure diseases or restore the whole-body balance; additionally, they provide a rich source for drug discovery. Owing to the complexity of the constituents in TCMs, reductionism, which individually considers each ingredient, cannot adequately explain the effects generated by an entire formula, and therefore faces several challenges comprehending and illustrating the mechanisms of herbal formulas. With the rapid development of bioinformatics, a combined attempt based on system/network biology and polypharmacology has created a new paradigm called network pharmacology (NP). In recent years, NP has been proposed as a promising approach for understanding herbal formulas1,2 and predicting potential new drugs or targets3,4,5. To some extent, NP-based research on the underlying therapeutic mechanisms of herbal formulas has been well-accepted6,7 and accounts for the strong compatibility of different scales of systems. A distinct advantage of network analysis is that it can excavate underlying information from multilevel interactions, especially when dealing with big data. Yu Ping Feng (YPF) decoction, a classic herbal formula, has been used in clinics to treat various diseases, such as respiratory tract infection8,9, nosocomial pneumonia10, chronic bronchitis11, and perennial allergic rhinitis12,13. As a safe and effective formula, it has been reported to relieve idiopathic sweating and augment appetite in end-stage cancer patients14, prevent viral infections, such as severe acute respiratory syndrome15,16, and have immunoregulatory or beneficial effects in asthma17,18, hepatic-fibrosis19, pulmonary-fibrosis20,21,22, and dermatitis models23. Since explanations for the clinical application of YPF decoction have not yet been illustrated by the conventional method (such as the in vitro activity assays of selected compounds) due to difficulties in the experimental design of complex systems, its underlying mechanism remains to be explored. Thus, this study aimed to develop an NP-based method for identifying potential therapeutic targets of traditional herbal formulas by taking YPF decoction as an example. This method gathered relevant bioinformatics information about YPF decoction to identify its underlying therapeutic targets and shed light on the NP-based mechanism of herbal formulas’ traditional clinical applications. The workflow of the NP-based approach and its application in YPF decoction is shown in below (Fig. 1). Figure 1 Figure 1 Workflow of the novel NP-based method. Full size image In this study, first, the putative targets of compounds from YPF decoction were enriched to pathways to create the targets–pathways interaction, and the diseases related to these pathways were collected to create the pathways–diseases interactions. Next, the targets–pathways interactions were used to generate a “targets–(pathways)–targets” (TPT) network with partitioned modules (clusters), in the TPT network, targets act as nodes, and they are connected by relevant pathways (edges). Based on the TPT network, we hypothesised that if a group of protein targets were gathered to form a module naturally, they might share a functional similarity. Herein, a contribution score (CS) was calculated to evaluate the contribution of each module to different diseases based on the network modules’ partition and pathways–diseases interactions. Further, the module with the highest CS to relevant diseases of YPF as well as the notable targets within the module were identified to explore potential therapeutic targets of YPF. The detailed procedures are more fully described in the methods section. Results and Discussion The targets–pathways–diseases interaction profile of YPF decoction YPF decoction consists of three herbs and a total of 352 compounds, which were collected with PubChem CID and they belonged to flavonoids, saponins, saccharides, and volatile components (including phenylpropanoids, monoterpenes, sesquiterpenes, diterpenes, and triterpenes). The rapid development of bioinformatics and computational sciences has provided a wide range of approaches for targets prediction of small molecules24. In this study, the putative targets of compounds were identified based on the similarity ensemble approach (SEA: http://sea.bkslab.org/)25. The potential targets of 250 compounds were successfully predicted by the SEA, with a total of 968 protein targets. (A worth-mentioning issue is about the high amount of 968 potential “targets” of YPF decoction, especially in terms of the relative small target pool of SEA. We tried to narrow the scope of potential “targets” of YPF decoction by setting up a pre-enrichment filter based on Tc threshold value. In this way, the number of “targets” for enrichment analysis sharply decreases from 968 to 304. However, the results of the most-important potential therapeutic targets and potential active compounds of YPF have no obvious changes. Thus, we still chose the post-enrichment filtering strategy, considering the advantage of this strategy with more informative reference to alternative targets. Throughout this paper, the term “targets” refers to the putative proteins to which the ingredients from YPF decoction directly bind, according to SEA predictions; however, it is not known whether targets are responsible for the therapeutic effect of YPF decoction. The term “therapeutic targets” refers to partial proteins from targets that have been identified to be responsible for the therapeutic effect of YPF decoction by our novel NP-based method.) Then, the target–pathway interactions were generated to elucidate the biological pathways that YPF decoction might impact; moreover, the overrepresented KEGG pathways of the predicted 968 protein targets of YPF were identified from enrichment analysis. A total of 132 highly related pathways with 549 relevant protein targets were obtained. Please refer to supplementary dataset for more detailed information about the compounds (Supplementary Table S1), putative targets (Supplementary Table S2) and relevant pathways (Supplementary Table S3). To understand what role the enriched pathways might play from a therapeutic angle, the Kyoto Encyclopedia of Genes and Genomes26 (KEGG; http://www.genome.jp/kegg/pathway.html) was used to identify the diseases with which each pathway might be involved. Since each pathway might participate in various bioprocesses in vivo, it is also possible that each pathway could be relevant to multi-diseases, and that each disease might be contributed to by multi-pathways. Next, a total of 43 diseases that the 132 pathways might impact were retrieved from the KEGG database to create the pathways–diseases interaction. Finally, the targets–pathways–diseases interaction profile of YPF decoction was generated. (Throughout this paper, the term “diseases” refers to the disease category of the diseases relevant to the pathways. Please refer to Supplementary Table S4 for detailed information about the pathways-diseases relationship.) Targets–(pathways) –targets (TPT) network and module identification Network modules (clusters or communities) refer to subnetworks whose nodes are more strongly connected to one another than to the rest of the network. That is to say, modules are sets of highly interconnected nodes. The identification of these modules is of crucial importance as they might help to uncover the hidden structural information within a network. Network modules have been utilised in NP-based TCM studies to uncover the combinations rule of herbal medicines27, chemical modules with similar structures28, and proteins with similar functions29,30,31. Especially in the context of TCMs as a typical complex system of multi-components and multi-targets, network modularity based of the Law of Like Attracts Like was considered to be a powerful way to sketch the complexity of TCMs. The complicated targets–pathways–diseases interaction of YPF in Fig. 2(a) showed a preliminary NP-based bioinformatics profile of YPF decoction, by employing the most frequently used network analysis method in research on traditional remedies for therapeutic targets’ prediction. The nodes represent putative targets of YPF (orange), relevant pathways (light blue), and diseases (deep purple). From Fig. 2(a), it is difficult to determine which targets might be most relevant to the therapeutic effect that YPF exerts; therefore, the novel TPT network was established to realise the recognition of therapeutic-based targets of YPF decoction, as shown in Fig. 2(b). Eight modules (modules 1–8) were detected and divided by the Louvain algorithm incorporated in the Gephi software at a resolution of 1.032,33. According to graph theory, the TPT network shown in Fig. 2(b) suggests some significant structural information, especially when considering the relatively separated modules with dense internal connections. Then, the inner interactions (edges representing the common pathways between every two nodes) of each target module were further investigated. Figure 2 Figure 2 (a) Targets–pathways–diseases network, and (b) Targets–(pathways)–targets network with modularity partition. Nodes in (a) represent targets of YPF decoction (orange), relevant pathways (light blue), and diseases (deep purple), while edges represent the interaction of targets–pathways–diseases. Nodes in (b) represent the targets of YPF decoction, while edges refer to relevant pathways of these targets. If two nodes are connected by an edge, this indicates that they have been enriched to participate in at least one of the same pathways. Eight modules (modules 1–8) presented in different colours were detected and partitioned by the Louvain algorithm incorporated in the Gephi software. Full size image The targets in module 1, module 2, and module 3 were connected by the “proteasome” (hsa03050), “lysine degradation” (hsa00310), and “aminoacyl-tRNA biosynthesis” (hsa00970) pathways, respectively; the inner connection of the targets in module 6 was “protein digestion and absorption” (hsa04974) and in module 8, it was the “insulin signaling pathway” (hsa04910). Three major modules (module 4, module 5 and module 7) could be easily identified from the network in Fig. 2(b). In addition, the inner connection of the targets in module 4 was the “neuroactive ligand-receptor interaction” (hsa04080) pathway. The connections between targets in module 5 were almost metabolism-related pathways, and the pathways in module 7 were pathways that immune-related, which was the most complicated. The results are shown in Table 1. From the results, we inferred that the target modules in the TPT network might represent the module shares’ different functional propensities. Thus, a contribution score (CS) was calculated from the newly established contribution-scoring algorithm to assess the relationship between each module and relevant diseases in the following study by integrating the modularity partition of the TPT network, target–pathway interactions, and pathway–disease interactions. The contribution-scoring algorithm is presented in the Method section in detail. Table 1 The inner interactions (edges) of each target module. Full size table Contribution score (CS) of modules toward diseases Following the partition of the modules, the CS was calculated to evaluate the contribution of each module to each disease by the contribution scoring algorithm, the results of which are shown in a heat map (Fig. 3); the modules are listed in a row and the diseases in a column. The gradually changing colours, from deep blue to white, represent the highest to the lowest CS, respectively. Figure 3 Figure 3 The contribution scores (CS) of targets in each module to different diseases. The deep blue and white refer to the highest to lowest CS, respectively. Full size image The result which was presented in Fig. 3 implies that each module might contribute differently to the 43 diseases, while the statistical results of the chi-square (χ2) test (p-value < 0.001) further indicate the significant association of modules 1–8 to these diseases. The contribution propensity of each module to diseases was further discussed as follows. It could be found from Fig. 3 that module 1, module 2, module 3, module 6, and module 8 commonly show a slight propensity toward diseases, but module 1 shows a notable contribution (CS = 1.000) to autism, compared with the other modules. In addition, it has been mentioned above that the targets in module 2 were connected by the “lysine degradation” pathway (hsa00310). Herein, the targets in module 2 might interfere with autism through lysine degradation, which has been already implicated in autism by previous studies34,35,36. Compared to other modules, the targets in module 4, which were connected by the neuro-relevant “neuroactive ligand-receptor interaction” pathway (hsa04080), contribute to hypopituitarism (CS = 1.000) and mental disorders (CS = 0.667) most. Though hypopituitarism is a chronic endocrine illness, it is still regulated by the nervous system37,38. Thus, we inferred that the targets in module 4 would exert neuro-relevant function when affected by compounds. The targets in module 5 were mainly connected by metabolism-relevant pathways, and they contribute to congenital developmental disorders (CS = 1.000), mental retardation (CS = 0.625), and inherited metabolic disease (CS = 0.503). Module 7 contributes to various diseases, including infectious disease (CS = 0.856), immune system disease (CS = 0.844), primary immunodeficiency (CS = 0.803), cancer (CS = 0.780), gastrointestinal disease (CS = 0.759), and autoimmune disease (CS = 0.750). Thus, it was inferred that the targets in module 7 would possess a higher possibility of exerting an immune-relevant function when affected by compounds, compared with other modules. To ensure the reliability of the results, targets in the network that are associated with the drugs approved for the market are used to validate the association between diseases and network modules. The approved drug information and drug–target interactions were extracted from the Drugbank database39 and DrugCentral database40. The therapeutic area of the drugs was determined by the Anatomical Therapeutic Chemical Classification System code (ATC code), which is attributed to a drug by the WHO Collaborating Centre (WHOCC) for Drug Statistics Methodology41. The ATC code classifies drugs according to the following five levels: level 1, the organ or anatomical system on which they act; level 2, the pharmacological action; levels 3 and 4, the chemical, pharmacological, and therapeutic subgroups; and level 5, the specific single drug or drug combination. In this study, the ATC code of the approved drugs was abbreviated to level 1 to understand the organ or anatomical system on which they act. For example, the ABL1 (tyrosine-protein kinase ABL1) has been reported to be targeted by the drug “dasatinib”, whose ATC code is L01XE06, and the connotation is shown as follows: Level 1 (L): antineoplastic and immunomodulating agents Level 2 (L01): antineoplastic agents Level 3 (L01X): other antineoplastic agents Level 4 (L01XE): protein kinase inhibitors Level 5 (L01XE06): dasatinib The level 1 abbreviation of the ATC code for dasatinib is L, which denotes that the therapeutic area of the drug dasatinib is antineoplastic and immunomodulating. Then, it was further inferred that the target ABL1 might be bound to exert antineoplastic and immunomodulating functions. The verification result showed that, in the TPT network, 40 targets were identified to be targeted by approved drugs to exert therapeutic effects on the nervous system, among which 30 targets (~75%) are located in module 4, seven in module 5, and only two in module 6 and module 7, respectively. This indicates that the targets in module 4 contribute most to the nervous system, which is consistent with the contribution-scoring result; in total, 54 targets were identified to be targeted by approved drugs to exert antineoplastic and immunomodulating activity. Among these, 34 targets (~65%) are located in module 7, 15 (~25%) in module 5, two in module 1 and module 4, respectively, and only one in module 8. This indicates that the targets in module 7 contribute most to the immune system, which is consistent with the contribution-scoring result. The verification results can also be found in Fig. 4. (Please refer to Supplementary Table S5 for detailed information of the targets with relevant approved drugs.) The chi-square (χ2) test (p-value < 0.001) of the results further indicates the significant association between network modules and diseases in terms of the targets of marketed drugs and the consistency in the main modules between predicted functions and drug indications in the real world. Figure 4 Figure 4 The count of approved drug targets in each module that are targeted to exert neuro-relevant function (N: Red) and antineoplastic or immunomodulating activity (L: Blue). Full size image The module that contributes most to YPF-relevant diseases To understand which network module contributes most to the therapeutic effects of YPF decoction, the CS of each module toward relevant diseases of YPF was selected for the following analysis. The relevant diseases of YPF decoction are immune system disease16,42,43,44, autoimmune diseases16,45,46, infectious diseases16,45,46, respiratory diseases20,46,47,48, and cancer49,50, according to the information retrieved from the publications, which is listed in Table 2. (The term “relevant diseases” refers to the diseases that are affected by YPF, while the term “irrelevant diseases” denotes those diseases that are relevant to the enriched pathways of YPF targets but have not been reported to be affected by YPF). The ratio of the CS of modules 1–8 to relevant diseases of YPF is presented in Fig. 5, which indicates that targets in module 7 contribute most to immune system diseases, infectious diseases, autoimmune diseases, and cancer; respiratory diseases are mainly contributed to by targets in module 6 and module 7. From a general perspective, YPF decoction might interfere with the targets in module 7 with a higher possibility of exerting therapeutic activity. Table 2 Reported diseases that are associated with YPF decoction. Full size table Figure 5 Figure 5 The contribution ratio of each module to relevant diseases of YPF decoction. Full size image YPF decoction is a classic formula with immune-regulatory effects, and module 7 showed the highest potential for curing its relevant disorders compared to other modules. Thus, the 199 targets (nodes) were considered to be the potential therapeutic targets of YPF decoction. Please refer to Supplementary Table S6 for more detailed information about the targets. Although YPF decoction is less known to exhibit anti-neoplastic activity directly in traditional applications, modern research has shown its potential in the treatment of idiopathic sweating of unknown cause in end-stage cancer patients14 and Lewis lung cancer in vitro50. There is extensive evidence showing a close relationship between inflammation and cancer51,52,53,54; that is, chronic inflammation contributes to cancer. Thus, the anti-neoplastic activity of YPF decoction might mainly occur via the regulation of the immune system. Identification of potential therapeutic targets in module 7 The above identification and CS of modules combined with statistical analysis were carried out to determine the therapeutic targets module that was highly associated with the therapeutic activity of YPF decoction. The following centrality analysis was conducted to evaluate the importance of the targets (nodes) within the therapeutic targets module from the perspective of the topological analysis of the TPT network. Four centrality indicators, i.e., degree centrality, closeness centrality, betweenness centrality, and eigenvector centrality, were used to measure the importance of the targets (nodes) from the different perspectives55,56,57. Then, an integrated indicator, namely target importance (TI), was calculated accordingly. Next, the top 10% therapeutic targets by TI were selected and recommended for further pharmacology study in the context of the mechanism illustration of YPF decoction. From a holistic perspective, YPF might interfere with targets in module 7 in a synergistic way, some of which (with higher centrality), as illustrated in Table 3, can be considered to be the most potential therapeutic targets of YPF. Moreover, the compounds relevant to these targets according to the SEA prediction can be recommended as potential active ingredients in YPF decoction in Table 3. Table 3 Top 10% targets in module 7 ranking by TI and relevant compounds. Full size table Previous studies have shown that YPF potentially attenuates allergic inflammation via the regulation of NF-kappa B activation42. Thus, it was inferred that YPF might also be relevant to NFKB1 according to the results of this study. Moreover, the top 10% potential therapeutic targets have been reported as a pleiotropic locus that is associated with cancer and a critical factor for acute inflammation associated with asthma, such as PRKCA58,59, EGFR60, AKT161, RAC162,63, RHAO64,65, HDAC166. In the activation by GC, NR3C1 can activate MAPK8 (a seed gene) and several MAPK family members, including MAPK3, which could trigger GAB1 to interfere with asthma disease67; PRKACA shares a close relationship with cardiovascular disease and cancer68; ADCY5 has been found to play a critical role in sensory neuropeptide release and neurogenic inflammation and it is critical in the inflammatory process69,70; CAMK2A71 is a calcium signaling molecule involved in cell growth and stress-signals integration, and it has been reported to play a tumor-supportive role in osteosarcomas; PLCG2 is an enzyme responsible for ligand-mediated signaling in cells of the hematopoietic system, and it plays a key role in the regulation of immune responses72, CREB173,74 and CAMK2B75 are also important for the regulation of immune responses; HRAS and NRAS, which belong to the p21 RAS subfamily of small GTPases, regulates cell proliferation, cytoskeletal organization, and other signaling networks, and they are among the most frequent targets of activating mutations activations in cancer76. Moreover, six of these top-ranking targets have been reported to be targeted by approved drugs: EGFR, HDAC1, MAPK3, NFKB1, PRKCA, and RXRA; except for NFKB1, the others are targeted to exert antineoplastic or immunoregulating effects. The drugs cetuximab, gefitinib, afatinib, olmutinib, and brigatinib act on EGFR, and they have been used in cancer therapy; vorinostat and romidepsin act on HDA1, and they are used in primary cutaneous T-cell lymphoma; arsenic trioxide acts on MAPK3, and it has been used to treat acute promyelocytic leukemia; midostaurin acts on PRKCA to treat gastrointestinal stromal tumors; and RXRA is targeted by alitretinoin and bexarotene with the indication of Kaposi’s sarcoma skin lesions and primary cutaneous T-cell lymphoma, respectively. All these previous studies provided validation for our novel TPT network-based study to some extent when using YPF as an example. As a result, targets with the highest centrality in this module are recommended as potential therapeutic targets of YPF decoction. Moreover, among them, there are eight most-important potential therapeutic targets, i.e., PRKCA, EGFR, NFKB1, AKT, CALM1, HDAC1, CREB1, and CAMK2B, which are identified by the relevance to YPF’s compounds with Tc value higher than the threshold of SEA method. The compounds associated with the top ranking targets The compounds that were predicted to bind to the top-ranking targets (potential therapeutic targets) by SEA are presented in Table 3. Thus, these compounds were recommended as the potential active ingredients of YPF decoction. Among them, the compounds with a Tc value greater than the threshold of SEA are mainly flavonoids, including oroxylin A, psoralen, quercetin, tectochrysin, wogonin, chrysin, and isorhamnetin. This result is consistent with a previous report that the flavonoid in YPF shares immunomodulatory effects9,42. Methods Data acquisition and processing YPF decoction is composed of Astragali Radix (Huang qi in Chinese, AR), Atractylodis Macrocephalae Rhizoma (Bai zhu in Chinese, AMR), and Saposhnikoviae Radix (Fang feng in Chinese, SR). First, the compound information of AR, AMR, and SR were collected from chemical databases, including the traditional Chinese medicines integrated database (TCMID, http://www.megabionet.org/tcmid/, which contains information from the TCM@Taiwan, traditional Chinese medicine information database (TCM-ID), herbal ingredients’ targets (HIT; accessed in May 2016) database, traditional Chinese medicine systems pharmacology (TCMSP; http://lsp.nwsuaf.edu.cn/tcmsp.php; accessed in May 2016) database, the cardiovascular disease herbal database (CVDHD, http://pkuxxj.pku.edu.cn/CVDHD/; accessed in May 2016), and the Chinese Academy of Sciences chemistry database (http://www.organchem.csdb.cn/scdb/; accessed in May 2016). Literature mining also supplemented the information from these databases. The common amino acids and compounds with high molecular weight, such as polysaccharides, were not included in this study. The names of compounds were standardised through SciFinder (https://scifinder.cas.org) according to the Chemical Abstracts Service (CAS) number. Additionally, these compounds’ PubChem CIDs were obtained from PubChem (https://pubchem.ncbi.nlm.nih.gov/) based on chemical names or structures. The protein targets of compounds from YPF decoction were predicted by similarity ensemble approach (SEA: http://sea.bkslab.org/)25, which is a chemical similarity searching-based prediction method recognised worldwide for its accuracy77,78. Though the SEA method is capable of accounting only ~2,800 potential active proteins as alternative binding targets, the scale is consistent with the recent recognition of the druggable genome (~3000)79,80,81. In addition, the limitation of a chemical similarity searching-based prediction method, including the SEA, would be that the prediction is based on the two-dimensional (2D) structure (SMILE information) of compounds, which might produce the same potential targets of the isomeric compounds. Herein, it is not enough to predict solely using SEA or any other chemical similarity searching-based prediction method, when isomeric compound domains in the research. The NP-based study of YPF decoction involved 13 pairs of isomeric compounds in 352 compounds, which would interfere with the result slightly, from a holistic perspective. Then, enrichment analysis is used to acquire target-pathway interactions which are essential for the further network construction and analysis. The targets were further enriched to KEGG pathways to create the targets–pathways relationship using the online tool STRING82 (http://string-db.org/). The enriched KEGG pathways with false rate discovery (FDR) < 0.01 were used in the subsequent research. The Kyoto Encyclopedia of Genes and Genomes26 (KEGG; http://www.genome.jp/kegg/pathway.html) was used to identify the diseases with which each pathway might be involved, the disease category information was used to represent the disease. The diseases that are affected by YPF decoction were retrieved from previous research. First, publications about the topic of YPF were extracted from the Web of Science, and the topic keywords were set as “yupingfeng”, “Yu-Ping-Feng”, or “gyokuheifusan” (the Japanese name of YPF decoction). Second, the relevant diseases’ information was identified, where available. Then, the diseases were classified according to the KEGG database. Network construction and module identification From a perspective of topology, the pathway-based associations between targets by enrichment analysis may be described as a network. A familiar representation is obtained by letting N be a set of nodes representing targets, and E be a set of edges where elements of E are unordered pairs of distinct nodes ni, nj representing a pathway-based link between a target pair {ni, nj}. The two sets together are called a simple network G = (N, E). Specifically speaking, a TPT network was established on the one-mode targets–targets interaction, which was transferred by Pajek software83 from the two-mode “targets–pathways” relationship; then, the TPT network was visualised and analysed by Gephi software84. In the TPT network, the nodes refer to protein targets, and each set of two nodes is connected by an edge, which indicates that they are both involved in at least one of the same pathways. Next, the Louvain algorithm for module identification incorporated in Gephi32 was adopted to explore the modularity structure of the TPT network. Since the targets in a module are connected by the same pathways or a functionally similar pathway, it was hypothesised that each module would exert a kind of specific function, and different modules would have effects on different diseases. Herein, an evaluation algorithm, named the contribution scoring algorithm, was established to weigh the contribution of each network module to various diseases and, further, to illustrate the TPT network. Contribution scoring algorithm On the basis of an existing algorithm for measuring the associated intensity of individual elements between interactive sets85,86, an indicator measuring the intensity of associating a specific network module with a specific disease, named a “contribution score (CS)” was established in this study by integrating the modularity partition of the TPT network, the targets–pathways interactions, and the pathways–diseases interactions. In this algorithm, the higher the value (contribution score), the greater the contribution made by the investigated module for a particular disease. Then, the targets in the module with a significant high value toward relevant diseases affected by YPF were considered to be its potential therapeutic targets. Figure 6 presents a diagram for illustrating the contribution scoring algorithm, where M refers to a set of modules, mi (i = 1, …, I); P refers to a set of pathways, pw (w = 1, …, W); and D refers to a set of diseases, dj (j = 1, …, J). In the TPT network, if uiw modules are connected by pw, each relevant mi contributes 1uiw to pw, namely the contribution of mi to pw, Cmipw=1uiw 85,86. If vwj pathways are relevant to dj, each relevant pw contributes 1vwj to dj, namely the contribution of pw to dj, Cpwdj=1vwj 85,86. Figure 6 Figure 6 A diagram illustrating the contribution scoring algorithm. Full size image The contribution scoring algorithm is defined as follows: Cmidj=∑pw∈XijCmipwCpwdj (1) where Xij is a subset of P and refers to the pathways that are relevant to mi and dj simultaneously, and where Cmidj refers to the CS of mi to dj, which is the sum of the contribution of mi to dj through all its relevant pw in Xij. The value of Cmidj varies from 0 to 1: the higher the value, the greater the contribution mi might make to dj; and all the modules contribute 1 to a particular disease. Below, an example is given to demonstrate the application of the contribution scoring algorithm. The Fig. 7(a) represents a TPT network with modules identified, and the Fig. 7(b) represents a pipeline which was abstracted from Fig. 7(a). Figure 7 Figure 7 An example demonstrating the contribution scoring algorithm. (a) a TPT network with modules identified: In the diagram, the nodes represent targets, and the edges represent pathways. Nodes T-A1-T-A9 in yellow refer to targets in module 1 (m1); nodes T-B1–T-B8 in blue refer to targets in module 2 (m2); nodes T-C1-T-C7 in green refer to targets in module 3 (m3), assuming that p1 − p5 are pathways that are relevant to d1. (b) a pipeline abstracted from (a). Full size image The TPT network was divided into three modules: module 1 (m1), module 2 (m2), and module 3 (m3). Since p1 and p2 connected the protein targets from only one module, m1, the contribution of m1 to p1 or p2 could be calculated as Cmipw=1uiw=1,(i=1;w=1,2) (2) As p3 and p4 were contributed to by the protein targets from m1 and m2, and from m2 and m3, separately, the contribution of m1 or m2 to m3, and from m1 or m3 to m4, could be calculated as Cmipw=1uiw=12,(i=1,2;w=3,ORi=1,3;w=4) (3) and since all three modules were related to p5, the contribution of m1, m2, or m3 to p5 could be calculated as Cmipw=1uiw=13,(i=1,2,3;w=5) (4) It was assumed that d1 was relevant to five pathways (p1 − p5) in the TPT network, herein, vwj=5,(w=1,2,3,4,5;j=1), and that each pathway contributed 15 to d1, namely: Cpwdj=1vwj=15,(w=1,2,3,4,5;j=1) (5) Ultimately, according to equation (1), the result of Cm1d1=∑pw∈X11Cm1pwCpwd1=23,X11={p1,p2,p3,p4,p5}, (6) Cm2d1=∑pw∈X21Cm2pwCpwd1=16,X21={p3,p5}, (7) Cm3d1=∑pw∈X31Cm3pwCpwd1=16,X31={p4,p5} (8) could be calculated: Cm1d1 shared the highest value, indicating that the targets in m1 contributed the most to d1 in the TPT network. That is, the targets in m1 might exert the most potential therapeutic proclivity toward d1 compared with the targets in other modules. Additionally, the total contribution of the three modules to d1 was 1. Statistical validation The CS, which reflected the contribution of each target module to relevant diseases, could be calculated from the established contribution scoring algorithm. To validate the hypothesis that different modules would have effects on different diseases, the association between network modules and therapeutic diseases was examined by chi-square (χ2) test. The verification of the contribution scoring algorithm based on the information of targets and approved drugs was also conducted by chi-square (χ2) test. Centrality analysis An integrated indicator, namely target importance (TI), was calculated accordingly as the equation below. TIi=14(DCimax(DC)+BCimax(BC)+CCimax(CC)+ECimax(EC)) (9) where TIi denotes the importance of target i, DCi, BCi, and CCi, and ECi denotes the degree, betweenness, closeness, and eigenvector centrality of target i, respectively; meanwhile, max(DC), max(BC), max(CC), and max(EC) denote the maximum degree, betweenness, closeness, and eigenvector centrality in the therapeutic target module. The value of TIi varied from 0 to 1, and the higher the value shared by a target, the more important the target was in the network from the topological perspective. Data availability The source data for analysis in this work are available in PubChem [https://pubchem.ncbi.nlm.nih.gov/], SEA [http://sea.bkslab.org/], KEGG [http://www.genome.jp/kegg/pathway.html]. Drugbank [https://www.drugbank.ca/], Drug Central [http://drugcentral.org/]. The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request. Conclusions A novel NP-based method that could help discover the core mechanisms of herbal formulas from complicated information and interactive relationships was presented in this study. YPF decoction was investigated by using the NP-based method, which integrated the “multi-targets–multi-pathways–multi-diseases” relationship into a TPT network and a contribution scoring algorithm. The modules in the TPT network with a similar function, especially coinciding with the already known application of YPF decoction, were inferred from the finding that the protein targets in module 7 were highly associated with the underlying mechanism of YPF decoction, especially those with highest centrality as presented in Table 3. Additionally, some of these targets have been reported to be relevant to the therapeutic activity that YPF might exert. Herein, the network-based method provided an easier and more reliable strategy for uncovering the potential therapeutic targets of YPF. In short, the proposed TPT network combined with the contribution scoring algorithm could be widely applied in the analysis of traditional medicines/formulas with robust effects in clinical care. And, importantly, the most therapeutically related modules coinciding with their clinical applications have a high chance of being associated with potential mechanisms worthy of further experimental investigation. Additional information Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. References 1. Liang, X. J., Li, H. Y. & Li, S. A novel network pharmacology approach to analyse traditional herbal formulae: the Liu-Wei-Di-Huang pill as a case study. Mol. Biosyst. 10, 1014–1022 (2014). Show context for reference 1 ArticlePubMed Google Scholar 2. Zhang, G. B., Li, Q. Y., Chen, Q. L. & Su, S. B. Network pharmacology: a new approach for Chinese herbal medicine research. Evid. Based Complement. Alternat. Med. 2013, 621423 (2013). Show context for reference 2 PubMedPubMed Central Google Scholar 3. Tao, W. Y. et al. Network pharmacology-based prediction of the active ingredients and potential targets of Chinese herbal Radix Curcumae formula for application to cardiovascular disease. J. Ethnopharmacol. 145, 1–10 (2013). Show context for reference 3 ArticlePubMed Google Scholar 4. Zhang, Q. R. et al. Natural formulas and the nature of formulas: Exploring potential therapeutic targets based on traditional Chinese herbal formulas. PLoS One 12, e0171628, https://doi.org/10.1371/journal.pone.0171628 (2017). Show context for reference 4 ArticlePubMedPubMed Central Google Scholar 5. Mao, Y. et al. Network pharmacology-based and clinically relevant prediction of the active ingredients and potential targets of Chinese herbs in metastatic breast cancer patients. Oncotarget 8, 27007–27021 (2017). Show context for reference 5 PubMedPubMed Central Google Scholar 6. Kibble, M. et al. Network pharmacology applications to map the unexplored target space and therapeutic potential of natural products. Nat. Prod. Rep. 32, 1249–1266 (2015). Show context for reference 6 ArticlePubMed Google Scholar 7. Xi, D. et al. State of the science: Cancer complementary and alternative medicine therapeutics research—NCI strategic workshop highlights of discussion report. J. Natl. Cancer Inst. Monogr. 2017, 62–67 (2017). Show context for reference 7 Article Google Scholar 8. Wang, S. J., Jiao, C. Y. & Sun, X. D. Clinical observation of effect of jiawei yupingfeng mixture for prevention and treatment of 100 children with repeated respiratory tract infection. Chinese Journal of Integrated Traditional and Western Medicine 29, 742–745 (2009). Show context for reference 8 PubMed Google Scholar 9. Du, C. Y. Q. et al. Yu Ping Feng San, an ancient Chinese herbal decoction containing Astragali Radix, Atractylodis Macrocephalae Rhizoma and Saposhnikoviae Radix, regulates the release of cytokines in murine macrophages. PLoS One 8, e78622, https://doi.org/10.1371/journal.pone.0078622 (2013). Show context for reference 9 ADSArticlePubMedPubMed Central Google Scholar 10. Yan, L. et al. Prevention of hospital-acquired pneumonia with Yupingfeng Powder in patients with acute cerebral vascular diseases: a randomized controlled trial. Journal of Chinese Integrative Medicine 8, 25–29 (2010). Show context for reference 10 ArticlePubMed Google Scholar 11. Huang, H. Q., Lu, F. & Lu, S. Y. Effect of jiawei yupingfeng powder on T-lymphocyte subsets in patients with senile chronic bronchitis in acute onset stage. Chinese Journal of Integrated Traditional and Western Medicine 24, 109–111 (2004). Show context for reference 11 PubMed Google Scholar 12. Shi, H. Y., Zhuang, Y. & Wang, X. Y. Effect of yupingfeng droppill in treatment of allergic rhinitis. China Journal of Chinese Materia Medica 39, 2364–2366 (2014). Show context for reference 12 PubMed Google Scholar 13. Li, M., Bai, H. & Li, S. L. Therapeutic effect of herbal medicines with the effects of Qi-boosting and exterior-securing on allergic rhinitis. Chinese Journal of Otorhinolaryngology in Integrative Medicine 15, 23–25 (2007). Show context for reference 13 Google Scholar 14. Chiu, S. C. et al. The therapeutic effect of modified Yu Ping Feng San on idiopathic sweating in end-stage cancer patients during hospice care. Phytother. Res. 23, 363–366 (2009). Show context for reference 14 ArticlePubMed Google Scholar 15. Lau, T. F. et al. Using herbal medicine as a means of prevention experience during the SARS crisis. Am. J. Chin. Med. 33, 345–356 (2005). Show context for reference 15 ArticlePubMed Google Scholar 16. Poon, P. M. K. et al. Immunomodulatory effects of a traditional Chinese medicine with potential antiviral activity: A self-control study. Am. J. Chin. Med. 34, 13–21 (2006). Show context for reference 16 ArticlePubMed Google Scholar 17. Wang, H. Z., Hong, M., Gui, L. L., Hua, Y. Q. & Xu, H. Q. Effect of Yupingfeng San against OVA-induced allergic asthma in mice. China Journal of Chinese Materia Medica 38, 1052–1055 (2013). Show context for reference 17 PubMed Google Scholar 18. Zhang, G. Q. & Liu, W. D. Add-on effect of Yu Ping Feng formula for childhood asthma: A meta-analysis of randomized controlled trials. Eur. J. Integr. Med. 9, 9–17 (2017). Show context for reference 18 Article Google Scholar 19. Wang, J. J. et al. Preventive effects of a fractioned polysaccharide from a traditional Chinese herbal medical formula (Yu Ping Feng San) on carbon tetrachloride-induced hepatic fibrosis. J. Pharm. Pharmacol. 62, 935–942 (2010). Show context for reference 19 ArticlePubMed Google Scholar 20. Li, L. C. et al. Total extract of Yupingfeng attenuates bleomycin-induced pulmonary fibrosis in rats. Phytomedicine 22, 111–119 (2015). Show context for reference 20 ArticlePubMed Google Scholar 21. Cui, W. et al. Total glycosides of Yupingfeng protects against bleomycin-induced pulmonary fibrosis in rats associated with reduced high mobility group box 1 activation and epithelial-mesenchymal transition. Inflammation Res. 64, 953–961 (2015). Show context for reference 21 Article Google Scholar 22. Li, L. C. & Kan, L. D. Traditional Chinese medicine for pulmonary fibrosis therapy: progress and future prospects. J. Ethnopharmacol. 198, 45–63 (2017). Show context for reference 22 ArticlePubMed Google Scholar 23. Nakatsukasa, H., Tago, F., Okamoto, T., Tsukimoto, M. & Kojima, S. Therapeutic effects of gyokuheifusan on NC/Nga mouse model of allergic dermatitis. J. Health Sci. 55, 516–524 (2009). Show context for reference 23 Article Google Scholar 24. Fang, J. S., Liu, A. L. & Du, G. H. Research advance in the drug target prediction based on chemoinformatics. Acta Pharmaceutica Sinica 49, 1357–1364 (2014). Show context for reference 24 PubMed Google Scholar 25. Keiser, M. J. et al. Relating protein pharmacology by ligand chemistry. Nat. Biotechnol. 25, 197–206 (2007). Show context for reference 25 ArticlePubMed Google Scholar 26. Kanehisa, M. & Goto, S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30 (2000). Show context for reference 26 ArticlePubMedPubMed Central Google Scholar 27. Li, S., Zhang, B., Jiang, D., Wei, Y. Y. & Zhang, N. B. Herb network construction and co-module analysis for uncovering the combination rule of traditional Chinese herbal formulae. BMC Bioinformatics 11, 1–12 (2010). Show context for reference 27 Google Scholar 28. Ding, F. et al. An analysis of chemical ingredients network of Chinese herbal formulae for the treatment of coronary heart disease. PLoS One 10, e0116441, https://doi.org/10.1371/journal.pone.0116441 (2015). Show context for reference 28 ArticlePubMedPubMed Central Google Scholar 29. Guimera, R. & Amaral, L. A. N. Functional cartography of complex metabolic networks. Nature 433, 895–900 (2005). Show context for reference 29 ADSArticlePubMedPubMed Central Google Scholar 30. Wilkinson, D. M. & Huberman, B. A. A method for finding communities of related genes. Proc. Natl. Acad. Sci. USA 101, 5241–5248 (2004). Show context for reference 30 ADSArticlePubMed Google Scholar 31. Spirin, V. & Mirny, L. A. Protein complexes and functional modules in molecular networks. Proc. Natl. Acad. Sci. USA 100, 12123–12128 (2003). Show context for reference 31 ADSArticlePubMed Google Scholar 32. Blondel, V. D., Guillaume, J. L., Lambiotte, R. & Lefebvre, E. Fast unfolding of communities in large networks. J. Stat. Mech.-Theory E. 2008, P10008 (2008). Show context for reference 32 Article Google Scholar 33. Lambiotte, R., Delvenne, J. C. & Barahona, M. Laplacian dynamics and multiscale modular structure in networks. arXiv.org, e-Print Arch., Phys. 8, 1770 (2008). Show context for reference 33 Google Scholar 34. Crawley, J. N., Heyer, W. D. & LaSalle, J. M. Autism and cancer share risk genes, pathways, and drug targets. Trends Genet. 32, 139–146 (2016). Show context for reference 34 ArticlePubMedPubMed Central Google Scholar 35. James, S., Shpyleva, S., Melnyk, S., Pavliv, O. & Pogribny, I. Complex epigenetic regulation of engrailed-2 (EN-2) homeobox gene in the autism cerebellum. Translational psychiatry 3, e232 (2013). Show context for reference 35 ArticlePubMedPubMed Central Google Scholar 36. Marrale, M., Albanese, N. N., Calì, F. & Romano, V. Assessing the impact of copy number variants on miRNA genes in autism by Monte Carlo simulation. PLoS One 9, e90947, https://doi.org/10.1371/journal.pone.0090947 (2014). Show context for reference 36 ADSArticlePubMedPubMed Central Google Scholar 37. Kornblum, R. N. & Fisher, R. S. Pituitary lesions in craniocerebral injuries. Arch. Pathol. 88, 242–248 (1969). Show context for reference 37 PubMed Google Scholar 38. Iovino, M. et al. Hypopituitarism in neurocritical patients: A case report. Endocr. metab. immune. 16, 28–31 (2016). Show context for reference 38 Google Scholar 39. Knox, C. et al. DrugBank 3.0: a comprehensive resource for ‘omics’ research on drugs. Nucleic Acids Res. 39, D1035–D1041 (2011). Show context for reference 39 ArticlePubMed Google Scholar 40. Ursu, O. et al. DrugCentral: online drug compendium. Nucleic Acids Res. 45, D932–D939 (2016). Show context for reference 40 ArticlePubMedPubMed Central Google Scholar 41. Santos, R. et al. A comprehensive map of molecular drug targets. Nat. Rev. Drug Discov. 16, 1–19 (2017). Show context for reference 41 Article Google Scholar 42. Shen, D. D. et al. Screening active components from Yu-Ping-Feng-San for regulating initiative key factors in allergic sensitization. PLoS One 9, e107279, https://doi.org/10.1371/journal.pone.0107279 (2014). Show context for reference 42 ADSArticlePubMedPubMed Central Google Scholar 43. Makino, T., Ito, Y., Sasaki, S., Fujimura, Y. & Kano, Y. Preventive and curative effects of Gyokuheifu-san, a formula of traditional Chinese medicine, on allergic rhinitis induced with Japanese cedar pollens in guinea pig. Biol. Pharm. Bull. 27, 554–558 (2004). Show context for reference 43 ArticlePubMed Google Scholar 44. Gao, J. et al. Antiinflammatory and immunoregulatory effects of total glucosides of Yupingfeng powder. Chin. Med. J. (Engl.) 122, 1636–1641 (2009). Show context for reference 44 Google Scholar 45. Du, C. Y. Q. et al. Yu Ping Feng San, an ancient Chinese herbal decoction, induces gene expression of anti-viral proteins and inhibits neuraminidase activity. Phytother. Res. 29, 656–661 (2015). Show context for reference 45 ArticlePubMed Google Scholar 46. Song, T. et al. Adjuvant treatment with yupingfeng formula for recurrent respiratory tract infections in children: A meta-analysis of randomized controlled trials. Phytother. Res. 30, 1095–1103 (2016). Show context for reference 46 ArticlePubMed Google Scholar 47. Wang, Z. G. et al. Yupingfeng pulvis regulates the balance of T cell subsets in asthma mice. Evid. Based Complement. Alternat. Med. 2016, 6916353, https://doi.org/10.1155/2016/6916353 (2016). Show context for reference 47 PubMedPubMed Central Google Scholar 48. Yang, Z. S. et al. Anti-inflammatory effect of Yu-Ping-Feng-San via TGF-beta 1 signaling suppression in rat model ofCOPD. Iran. J. Basic Med. Sci. 19, 993–1002 (2016). Show context for reference 48 PubMedPubMed Central Google Scholar 49. Huang, J. H. et al. Effect of Yupingfeng granules on HA and Foxp3(+) Treg expression in patients with nasopharyngeal carcinoma. Asian Pac. J. Trop. Med. 8, 662–664 (2015). Show context for reference 49 Google Scholar 50. Luo, Y. B. et al. NK Cell-Dependent Growth Inhibition of Lewis Lung Cancer by Yu-Ping-Feng, an Ancient Chinese Herbal Formula. Mediators Inflamm. 2016, 3541283, https://doi.org/10.1155/2016/3541283 (2016). Show context for reference 50 PubMedPubMed Central Google Scholar 51. Zhang, X. et al. Exosomes for immunoregulation and therapeutic intervention in cancer. J. Cancer 7, 1081–1087 (2016). Show context for reference 51 ArticlePubMedPubMed Central Google Scholar 52. Marx, J. Inflammation and cancer: the link grows stronger: research into a long-suspected association between chronic inflammation and cancer reveals how the immune system may be abetting tumors. Science 306, 966–969 (2004). Show context for reference 52 ArticlePubMed Google Scholar 53. Guo, Y. C. et al. Multiscale modeling of inflammation-induced tumorigenesis reveals competing oncogenic and onco-protective roles for inflammation. Cancer Res. 77, 6429–6441 (2017). Show context for reference 53 ArticlePubMed Google Scholar 54. Hofseth, L. J. & Ying, L. Identifying and defusing weapons of mass inflammation in carcinogenesis. BBA-rev. cancer 1765, 74–84 (2006). Show context for reference 54 Google Scholar 55. Wasserman, S. & Faust, K. Social network analysis: Methods and applications. Vol. 8 (Cambridge university press, 1994). Show context for reference 55 56. Hu, Y. J., Scherngell, T., Man, S. N. & Wang, Y. T. Is the United States still dominant in the global pharmaceutical innovation network? PLoS One 8, e77247, https://doi.org/10.1371/journal.pone.0077247 (2013). Show context for reference 56 ADSArticlePubMedPubMed Central Google Scholar 57. Heymann, S. In Encyclopedia of social network analysis and mining. 612–625 (Springer, 2014). Show context for reference 57 58. Murphy, A. et al. PRKCA: A positional candidate gene for body mass index and asthma. Am. J. Hum. Genet. 85, 87–96 (2009). Show context for reference 58 ArticlePubMedPubMed Central Google Scholar 59. Paraboschi, E. M. et al. miR-634 is a Pol III-dependent intronic microRNA regulating alternative-polyadenylated isoforms of its host gene PRKCA. Biochim. Biophys. Acta, Gen. Subj. 1861, 1046–1056 (2017). Show context for reference 59 Article Google Scholar 60. Hwang, S. et al. A protein interaction network associated with asthma. J. Theor. Biol. 252, 722–731 (2008). Show context for reference 60 ArticlePubMed Google Scholar 61. Di Lorenzo, A., Fernandez Hernando, C., Cirino, G. & Sessa, W. C. Akt1 is critical for acute inflammation and histamine-mediated vascular leakage. Proc. Natl. Acad. Sci. USA 106, 14552–14557 (2009). Show context for reference 61 ADSArticlePubMed Google Scholar 62. Leavy, O. Asthma and allergy: Apoptotic cell clearance lets you breathe easy. Nat. Rev. Immunol. 13, 157, https://doi.org/10.1038/nri3404 (2013). Show context for reference 62 ArticlePubMed Google Scholar 63. Henson, P. M. & Bratton, D. L. Allergy: airway epithelial Rac1 suppresses allergic inflammation. Curr. Biol. 23, R104 (2013). Show context for reference 63 ArticlePubMed Google Scholar 64. Kume, H. RhoA/Rho-kinase as a therapeutic target in asthma. Curr. Med. Chem. 15, 2876–2885 (2008). Show context for reference 64 ArticlePubMed Google Scholar 65. Chiba, Y., Matsusue, K. & Misawa, M. RhoA, a possible target for treatment of airway hyperresponsiveness in bronchial asthma. J. Pharmacol. Sci. 114, 239–247 (2010). Show context for reference 65 ArticlePubMed Google Scholar 66. van den Bosch, T., Kwiatkowski, M., Bischoff, R. & Dekker, F. J. Targeting transcription factor lysine acetylation in inflammatory airway diseases. Epigenomics 9, 1013–1028 (2017). Show context for reference 66 ArticlePubMed Google Scholar 67. Sharma, A. et al. A disease module in the interactome explains disease heterogeneity, drug response and captures novel pathways and genes in asthma. Hum. Mol. Genet. 24, 3005–3020 (2015). Show context for reference 67 ArticlePubMedPubMed Central Google Scholar 68. Turnham, R. E. & Scott, J. D. Protein kinase A catalytic subunit isoform PRKACA; History, function and physiology. Gene 577, 101–108 (2016). Show context for reference 68 ArticlePubMed Google Scholar 69. Nemeth, J. et al. Effect of pituitary adenylate cyclase activating polypeptide-38 on sensory neuropeptide release and neurogenic inflammation in rats and mice. Neuroscience 143, 223–230 (2006). Show context for reference 69 ArticlePubMed Google Scholar 70. Shen, X. et al. Elucidation of the anti-inflammatory mechanisms of Bupleuri and Scutellariae Radix using system pharmacological analyses. Mediators Inflamm. 2017, 1–10 (2017). Show context for reference 70 Article Google Scholar 71. Wang, S. et al. In Proceedings of the American Association for Cancer Research Annual Meeting, 4787 (2017). Show context for reference 71 72. Zhou, Q. et al. A hypermorphic missense mutation in PLCG2, encoding phospholipase Cgamma2, causes a dominantly inherited autoinflammatory disease with immunodeficiency. Am. J. Hum. Genet. 91, 713–720 (2012). Show context for reference 72 ArticlePubMedPubMed Central Google Scholar 73. Wen, A. Y., Sakamoto, K. M. & Miller, L. S. The role of the transcription factor CREB in immune function. The Journal of Immunology 185, 6413–6419 (2010). Show context for reference 73 ArticlePubMed Google Scholar 74. Bartel, S. et al. Pulmonary microRNA profiles identify involvement of Creb1 and Sec. 14l3 in bronchial epithelial changes in allergic asthma. Sci. Rep. 7, 46026, https://doi.org/10.1038/srep46026 (2017). Show context for reference 74 ADSArticlePubMedPubMed Central Google Scholar 75. Zimmerer, J. M., Lesinski, G. B., Radmacher, M. D., Ruppert, A. & Carson, W. E. STAT1-dependent and STAT1-independent gene expression in murine immune cells following stimulation with interferon-alpha. Cancer Immunol. Immunother. 56, 1845–1852 (2007). Show context for reference 75 ArticlePubMed Google Scholar 76. Oliveira, J. B. et al. NRAS mutation causes a human autoimmune lymphoproliferative syndrome. Proc. Natl. Acad. Sci. USA 104, 8953–8958 (2007). Show context for reference 76 ADSArticlePubMed Google Scholar 77. Keiser, M. J. et al. Predicting new molecular targets for known drugs. Nature 462, 175–182 (2009). Show context for reference 77 ADSArticlePubMedPubMed Central Google Scholar 78. Lounkine, E. et al. Large-scale prediction and testing of drug activity on side-effect targets. Nature 486, 361–367 (2012). Show context for reference 78 ADSArticlePubMedPubMed Central Google Scholar 79. Hopkins, A. L. & Groom, C. R. The druggable genome. Nat. Rev. Drug Discov. 1, 727–730 (2002). Show context for reference 79 ArticlePubMed Google Scholar 80. Russ, A. P. & Lampel, S. The druggable genome: an update. Drug Discov. Today 10, 1607–1610 (2005). Show context for reference 80 ArticlePubMed Google Scholar 81. Griffin, R. et al. Glimmers in illuminating the druggable genome. Nat. Rev. Drug Discov (2018). Show context for reference 81 82. Szklarczyk, D. et al. The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res. 39, D561–D568 (2011). Show context for reference 82 ArticlePubMed Google Scholar 83. Batagelj, V. & Mrvar, A. In Graph Drawing. (eds Mutzel, P., Jünger, M. & Leipert, S.). 477–478 (Springer, 2002). Show context for reference 83 84. Bastian, M., Heymann, S. & Jacomy, M. In International AAAI Conference on Weblogs and Social Media Vol. 8. 361–362 (2009). Show context for reference 84 85. Ahlswede, R., Cai, N., Robert Li, S. Y. & Yeung, R. W. Network information flow. IEEE Trans. Inf. Theory 46, 1204–1216 (2000). Show context for reference 85 MathSciNetArticleMATH Google Scholar 86. Suo, T. C. et al. Combining chemical profiling and network analysis to investigate the pharmacology of complex prescriptions in traditional Chinese medicine. Sci. Rep. 7, 40529, https://doi.org/10.1038/srep40529 (2017). Show context for reference 86 ADSArticlePubMedPubMed Central Google Scholar Download references Acknowledgements We sincerely acknowledge Ms. Xiangjun Kong for her technical assistance, and Prof. Jin Qi and Mr. Huiheng Lin for their valuable suggestions. We also thank Ms. Kunmeng Liu for her valuable advice on statistical analysis, and the Science and Technology Development Fund of Macao SAR and the University of Macau for financial support for projects 013-2015-A1, MYRG2015-00172-ICMS-QRCM, MYRG2015-00145-ICMS-QRCM, and MYRG2016-00144-ICMS-QRCM for this research. Author information Affiliations State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences, University of Macau, Macau SAR, China Huali Zuo, Qianru Zhang & Yuanjia Hu School of Pharmacy, Zunyi Medical University, Guizhou, China Qianru Zhang Research Center for Traditional Chinese Medicine Complexity System, Shanghai University of Traditional Chinese Medicine, Shanghai, China Shibing Su & Qilong Chen School of Chemistry and Chemical Engineering, Chongqing University, Chongqing, China Fengqing Yang Contributions H.Y.J. and Y.F.Q. conceived and designed the study; Z.H.L. and Z.Q.R. analysed the raw data; S.S.B. and C.Q.L. participated in the study design; H.Y.J. and Y.F.Q. reviewed the manuscript and carried out extensive revisions to the manuscript. All authors read and approved the final content. Competing Interests The authors declare no competing interests. Corresponding authors Correspondence to Fengqing Yang or Yuanjia Hu. Electronic supplementary material Supplementary Table S1 Supplementary Table S2 Supplementary Table S3 Supplementary Table S4 Supplementary Table S5 Supplementary Table S6 Rights and permissions Creative Commons BY Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. About this article Publication history Received 17 January 2018 Accepted 03 July 2018 Published 30 July 2018 DOI https://doi.org/10.1038/s41598-018-29764-1 Share this article Anyone you share the following link with will be able to read this content: https://rdcu.be/4mvLCopy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative Subjects AllergyAsthmaDrug development