Skip to main content

Systematic synergy modeling: understanding drug synergy from a systems biology perspective


Owing to drug synergy effects, drug combinations have become a new trend in combating complex diseases like cancer, HIV and cardiovascular diseases. However, conventional synergy quantification methods often depend on experimental dose–response data which are quite resource-demanding. In addition, these methods are unable to interpret the explicit synergy mechanism. In this review, we give representative examples of how systems biology modeling offers strategies toward better understanding of drug synergy, including the protein-protein interaction (PPI) network-based methods, pathway dynamic simulations, synergy network motif recognitions, integrative drug feature calculations, and “omic”-supported analyses. Although partially successful in drug synergy exploration and interpretation, more efforts should be put on a holistic understanding of drug-disease interactions, considering integrative pharmacology and toxicology factors. With a comprehensive and deep insight into the mechanism of drug synergy, systems biology opens a novel avenue for rational design of effective drug combinations.


Drug combinations have been widely used in treating complex diseases such as cancer, HIV and cardiovascular diseases. For example, cisplatin and sabarubicin are co-administrated in lung cancer [1], combination of aspirin and dipyridamole is used to reduce the risk of stroke [2], and a novel FDA-approved drug combination: Stribild TM [3] will been applied in the HIV-1 treatment. According to the Drug Combination Database (DCDB) [4], there have been 330 FDA-approved and 1033 investigational drug combinations until the second version. With more and more researches on this new form of therapy, drug combinations have been proven a promising way to combat complex diseases [59].

The most prominent benefit of a drug combination is the synergy effect among different drugs, where synergy means that the overall therapeutic effect of the combination is greater than the sum of effects caused by individual components [10]. In order to discover new combinations, quantification of drug synergy is one of the essential works. However, synergy quantification is not an easy task. Different dose–response methods, including Loewe additivity [11, 12], Bliss independence [13], and Chou-Talalay method [14] have been proposed. The main concept of these methods is to find out whether the observed combination effect or response departs from the expected effect, as reviewed in [10, 15, 16]. The dose–response methods play key roles in synergy quantification. However, no single method can describe synergy under all possible conditions. Besides, none of these methods has the overwhelming superiority under all conditions [10]; different forms of dose–response methods may even produce different results [17]. At present, it is still impractical to screen all possible drug combinations for different diseases. In addition, these dose–response methods cannot explain the potential synergy mechanism for drug combinations. Such problems suggest that the identification of synergistic drug combinations remains a non-trivial and challenging task. As “network-based drug discovery is taking the pharmaceutical industry into a new age” [18], systems biology and computational technologies have provided a powerful tool for multi-target drug discovery [1921]. Consequently, systems biology has become a compensative way to predict novel synergistic drug combinations and to illustrate underlying mechanism. Such predictions can help accelerate the generation of hypotheses on possible synergistic drug combinations and guide the experimental conductions.

Systems biology allows us to describe the interactions among different molecules from a comprehensive network perspective. Currently, some systems biology approaches have been applied in drug synergy researches to predict the complex human body responses to drug combinations, to recognize the specific biological network features which are more likely to produce synergistic outcomes, and to illustrate the complicated molecular interactions which may contribute to synergistic mechanism [2225]. Understanding the explicit relationships among different objects involved in drug combinations allows us to appreciate how drugs interact with each other to produce synergistic therapeutic effects.

In the future, systematic synergy modeling will provide an assistant tool independent of screening data to accelerate the generation of hypotheses and guide the experimental activities. Here we review the current status of this new field and describe how, in our view, systems biology models (as classified in Table 1) can help us to develop a predictive and mechanistic model of drug synergy.

Table 1 A classification of systematic drug synergy prediction models


Predicting synergy combinations based on the PPI network

The traditional “one drug, one target” therapeutic mode can only yield limited effects on complex diseases, because these diseases are mostly controlled by complex biological networks which involve compensative biological processes [26, 27]. On the contrary, drug combinations, with multiple targets belonging to interlinked processes, can combat the systematic pathological states through the cooperative mechanism. Researchers have recently applied network-based systems biology approaches to predict potential synergistic drug combinations, as well as to reveal the synergy mechanism from the perspective of biological molecular network. Most of such researches based on the hypothesis that synergistic drugs tend to have effects on targets with distinguishing network topological features [2831].

One of the most frequently used biological molecular networks is the protein-protein interaction (PPI) network which is based on intentional physical or functional associations between proteins. The PPIs are the core of the whole interactomics system of any living organism and are also the foundation of more highly specified networks such as a disease-specific molecular network. The disease-specific PPI network often refers to a protein network in which the proteins are related with one certain disease [3234]. Based on a disease-specific PPI network, a heuristic synergy score can be calculated based on the topology connection and centrality of drug targets to predict the possible synergy degree. Other associated scores can also be coupled with to adjust the topology-based scores. For example, two evaluation scores, topology score (TS) and agent score (AS) were created to evaluate drug synergy for two drug combination by an algorithm termed NIMS (network target-based identification of multicomponent synergy) [28] (as outlined in Fig. 1). The TS was based on the topology relationship among targets from different components in the context of a previously defined disease-specific PPI network. The TS for two drugs will be high if they target on proteins which are close to each other on the PPI network and with high network centralities. The AS was to measure the phenotype similarity between the corresponding phenotypes of drugs. This phenotype similarity quantifies the overlap of their OMIM descriptions and is calculated by a text mining method. At last, a predictive synergy score was calculated as the product of TS and AS. Some synergistic ingredients from anti-angiogenic traditional Chinese medicine were successfully recognized by this score. Likewise, in the background of a disease-specific PPI network integrated of different data sources, Vitali et al. [29] assigned a score called Topological Score of Drug Synergy (TSDS) to evaluate the multi-target combination effects based on the target topological features including node reachability, global effect and synergistic effect. This work has been applied to identify novel targets which may lead to synergistic outcomes for Type 2 Diabetes Mellitus. To cover more knowledge about drug behavior, Huang et al. [30] proposed another synergy evaluation tool: “DrugComboRanker”. In this approach, the synergistic score incorporated topological relatedness of targets in the signaling network, semantic similarity of gene ontologies, and dissimilarity of the gene expression profiles of different drugs. The method was assessed on lung adenocarcinoma and endocrine receptor positive breast cancer and most of the top-ranked synergistic combinations were confirmed by literature reports.

Fig. 1
figure 1

A schematic example of the PPI network-based synergy score calculation. The calculation strategy in this example was from the method described in [28]. The green and red triangles represent two different drugs. Each drug is targeting on certain protein targets denoted as circles with the same color as the corresponding drug. All the targets are mapped into a disease-specific PPI network, in which purple circles are pathological proteins, while white circles are other high-related proteins, and different proteins are connected by PPIs. Based on the disease-specific PPI network (G), topology score TS(A, B) can be calculated by the topology connection and centrality for targets of drug A (TA) and targets of drug B (TB), where IP means a integrative topological parameter calculated by integrating betweenness, closeness and a variant of the eigenvector PageRank through PCA, IP A (i) is calculated for the target i of drug A and IP B (i) for the target j of drug B, and min(d i,j ) represents the minimum shortest path from the ith target of drug A to all targets of drug B. Another agent score (AS), in which P i,j represents the similarity score between the ith phenotype of drug A and the jth phenotype of drug B, can also be coupled with to adjust the topology-based score. At last, a synergy score can be calculated by the product of TS and AS

PPI network-based method is an effective and straightforward way to elucidate the inter-relationship between a complex disease and drug interventions. In a generous way, Wang et al. [31] have indicated that drug combinations tend to target proteins which are closer in the genetic interaction networks comparing against random combinations. In network-based synergy evaluation, target acquisition and disease network construction are two basic steps. By mapping drug targets onto disease-related molecular networks, it can not only help reveal the mechanisms of action for each drug, but also highlight how different drugs in a combination cooperate with each other to achieve synergy effects. The corresponding topology relationships of drug targets are of importance for estimating synergy among combinations. However, it is a fact that the behavior of a drug is not only determined by its targets; other factors, like drug side effects, drug-resistance and other pharmacokinetic characteristics, will also be the root of drug synergy. Therefore, a more comprehensive model which is not limited to PPIs should be constructed to model the complicated interactions involved in the mechanism of action for drug combination.

Exploring synergy mechanism from the perspective of pathway

A pathway refers to a series of chemical reactions or molecular interactions within a cell that generate a specific product, response or change in the cell [35]. Although the PPI network-based approaches described above can help discover potential synergistic drug combinations, pathways, which can represent specific parallel, cross-talk or feedback structures of molecular networks, are even more powerful for the explanation of synergy mechanism in more details. Different drugs can have influences on the same or different pathways, and induce synergy effects produced by targets aggregating at specific pathway structures, as revealed in [36]. Modeling the dynamic changes and network structures of pathways can help improve our understanding on the potential synergy mechanism for drug combinations.

Dynamic pathway simulation

Dynamic simulation of specific pathway under drug interferences can help us understand the dynamic behavior of drugs, and provide insights into the synergy mechanism of certain drug combinations. Such pathway simulation models are often represented in a network form: nodes stand for concentrations or activity levels of pathway components (e.g., gene, protein or metabolite) and edges reflect the interaction of one node on the time derivative of another (Fig. 2a). Moreover, dose–response data can be obtained from the simulation results (Fig. 2b). Based on these results, synergy can be quantified by traditional dose–response methods like Bliss independence and Loewe additivity.

Fig. 2
figure 2

An example of pathway simulation model. a. A network-based pathway simulation model of the EGFR signaling pathway. Circles represent different components in a pathway. Triangles are drugs. Edges represent the reactions among different elements of the pathway, and the quantitative dynamic changes for each element can be described by one ODE or other computational models. This pathway network is adapted from Klinger et al. [44]. b. Dose–response plots. The dose–response curves are based on the pathway simulation results for each individual drug (a, b) and the drug combination (ab) respectively, where the half maximal inhibitory concentration (IC50) can be utilized to quantitatively assess drug synergism by dose–response methods like Loewe additivity. If the ratios of two drugs in the combination are both 0.5, an interaction index can be calculated as: \( \frac{0.5\times IC{50}_{AB}}{IC{50}_A}+\frac{0.5\times IC{50}_{AB}}{IC{50}_B} \), where the denominators represent doses of each drug in the combination that yield a half inhibition effect, while the numerators represent the dose of each drug to produce the same effect when given alone. If this index is larger than one, two drugs are antagonistic, if it is less than one, the combination is synergistic. For the given example, this index is less than one, so drug A and B are synergistic

Classical chemical kinetics utilizes ordinary differential equations (ODE) to describe system dynamics [3739]. A system of ordinary differential equations is the most fundamental way to quantitatively simulate the dynamic response of each individual molecule in a pathway under different conditions like the interventions of drug combinations [4043]. Based on a kinetic ODEs model of the epidermal growth factor receptor (EGFR) signaling pathway, Araujo et al. [40] simulated the effects of combinatorial kinase inhibitors and discovered that a simultaneous inhibition of multiple nodes in a signaling cascade with small molecule kinase inhibitors will be a new promising combination strategy for cancer treatment. Modification of classical kinetic ODE models can improve the ability in reproducing cellular dynamics. Nelander et al. [41] derived a non-linear ODE model to capture the important phenomena which cannot be described by simple linear model, like the epistasis and saturation effects in a cellular system. Based on this model, the authors can predict the quantitative outcomes of combinatorial perturbations in breast cancer cells. Likewise, Miller et al. [43] also utilized a non-linear ODE model to quantitatively simulate the dynamic changes of signaling pathways in dedifferentiated liposarcoma cells. Integrating with the experimental drug combination screening, the simulation helps to explain that the synergy of CDK4 and IGF1R inhibitors may depend on the activity of AKT pathway.

Except for the ODE models, other computational methods have also been applied in pathway modeling for drug combination analysis. A modeling framework based on the modular response analysis (MRA) was developed to build a simulation model for the EGFR signaling pathway [44]. The MRA model can quantitatively simulate the dynamic changes triggered by feedback cycle, feedforward loop as well as cross talk structures in pathways. Based on the MRA model, Klinger et al. [44] deduced that downstream mutations do not necessarily invalidate upstream drug treatment if a suitable downstream inhibitor were combined. Different researches pay attention to different types of pathway structures. Yan et al. [45] proposed a mathematical model to simplify the serial and parallel structures of biological pathways. Besides, Asfar et al. [46] utilized the network modeling based on the Ingenuity Pathway Anallysis (IPA) Software which considers the microarray datasets at different time points to understand central synergy mechanisms for the combination of MI219 and oxaliplatin. Sometimes, simulations only based on the disease-specific pathway are not enough to describe the drug synergy mechanism, extra pathways which can indirectly affect the therapeutic process should also be considered. In the study of Li et al. [47], which aimed to investigate on the role of immune system response in mediating anti-influenza drug synergy, both of the influenza A virus life cycle pathway and two immune-related pathways were simulated by a series of delay differential equations.

Synergy-specific network motifs in pathways

The above-mentioned pathway simulation models provide a more concrete way to predict effective drug combinations. What’s more, they offer tools to discover which types of structures in pathways are more inclined to generate drug synergies, thus illustrating the mechanism of synergy at the pathway level. However, the pathway simulation models are often hampered by their incompleteness and complexity. To address this problem, reduced or abstract network models which can capture key dynamical properties of a network can be used. These reduced models, termed as network motifs, are composed of same number of elements and are used to describe distinct connectivity patterns that occur frequently in the whole network or different pathways [4850] (Fig. 3). Previous researches have begun to depict the dynamics of network motifs by kinetic ODE methods, and describe how variations in the kinetic parameters affect the human body responses to drug interventions [5052]. Different studies may produce distinctive forms of motifs considering their specific focuses. In the study of Zhang, et al. [53], based on ten three-node small network motifs covering main pathway sub-structures like positive feedback loop, negative feedback loop, positive auto-regulation, and feed-forward loop, they recognized which kind of concurrent variations could act synergistically (or antagonistically) to alter the responses of the motifs and which type of motifs were more likely to lead to synergistic (or antagonistic) responses. They discovered that combinations may be more probably to exhibit synergistic effects when targeting on a negative feedback loop motif or a mutual inhibition loop motif. As another example, Yin et al. [54] used three-node enzymatic motifs to study about the synergistic or antagonistic interactions within drug combinations considering different motif topologies and parameters. This study found that drug combination effects largely depend on network topology, and the kinetic parameter variations only yield limited influences on the drug combination effects. The specific synergistic or antagonistic network motifs for drug combinations were identified: synergistic motifs often encompass serial and parallel or mix-type structures, while antagonistic motifs are mostly with a positive feedback loop and a downstream link.

Fig. 3
figure 3

Examples of network motifs in pathways. a. Positive feedforward; b. Positive feedback; c. Negative feedforward; d. Negative feedback. S represents stimulus, P1, P2 are proteins in a pathway, and R represents the effector protein. Different types of the directed edges represent activation or inhibition effects between two connected elements. For example, P1 has an activation effect on P2 in Fig. A, while P1 has an inhibition effect on R in Fig. D. The study of Zhang et al. [53] discovered that combinations may be more probably to exhibit synergistic effects when targeting on a negative feedback motif. The study of Yin et al. [54] found antagonistic drug interactions are mostly related with a positive feedback loop. This figure is reproduced from [50], with permission from Elsevier (RightsLink license number: 3687391088768)

Although the pathway or network motif simulation results may not fully recreate the absolutely real biological phenomena, considering the network structures of the pathways which are more inclined to generate synergistic effects, these models have also come out with some hypotheses about the synergy mechanism underlying drug combinations: drugs targeting on certain types of structures, like cross-talk, parallel or mutual inhibition loop in pathways will lead to synergistic drug combination effects. Although these hypotheses may not hold true in all cases, they can still provide guidance for the rational design of potentially effective drug combinations. Some studies can reduce the combinatorial drug screening space by selecting drugs with targets on compensatory or inter-linked pathways [5557].

However, one limitation of these pathway models is that we only partially understand the pathways of most diseases or drugs, as pointed out in [58]. With the accumulation of the “omic” data and various reverse engineering methods [59, 60], more and more disease-specific pathways can be reconstructed. These new pathways will provide more bases for drug synergy mechanism illustration. However, most of the existing researches focus on the target-pathway (or motifs) interactions, ignoring other drug-organism interactions (e.g., drug absorption or metabolism). To obtain a more comprehensive description of synergy mechanism, future studies should pay more attention to construction of comprehensive network models considering not only drug targeting pathways but also other drug behaviors like absorption, metabolism, transportation, and even toxicity-related pathways.

Computational models based on integrated drug similarities

In addition to the network-based synergy principles, other forms of synergy mechanism, which can be owed to certain chemical or pharmacological properties such as chemical structure, drug indication, and drug side effects, should also be characterized. In order to understand drug synergy from a more global perspective, different types of drug properties need to be considered simultaneously. Some computational methods have been proposed to integrate various similarities based on different drug properties to build classification models for drug combinations.

Drugs or drug combinations can be represented as their corresponding features described by multiple drug properties. Thereby different machine learning methods can be utilized in the prediction or classification for synergistic drug combinations. For example, Zhao et al. [61] developed a computational model for the prediction of effective drug combinations based on the correlations between candidate combinations and confirmed combinations. In this study, each drug pair was featured as the pairwise combinations of their comprehensive molecular and pharmacological elements, including drug targets and their corresponding downstream pathways, medical indication areas, Anatomical Therapeutic Chemical (ATC) codes, and side effects. Based on the enrichment levels of the candidate drug pair features among all the drug pair features from known drug combinations, Zhao et al. [61] have predicted novel pairwise drug combinations, among them 69 % can be validated by literature reports. Another research used associations of drugs in a combination to describe drug combination features [62]. Considering three different kinds of drug relationships including chemical interactions between drugs in the combination, protein interactions between the targets of drugs, and target enrichment of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, each drug combination was represented as a high-dimensional numeric vector. Then a random forest classification model was constructed based on the key features extracted from high-dimensional numeric vectors to predict synergistic drug combinations. A 5-fold crossover validation showed this classification model was well performed with a Matthews correlation coefficient of 0.673.

Prediction or classification models based on properly selected drug properties can help to identify potential synergistic drug combination. The most important part of such computational models is the feature calculation for drug combinations. Informative features will improve the prediction accuracy and help to reveal synergy mechanism. Besides, samples collection is another important step. Consequently, more efforts should be put in the elaborate design of effective features for identification of drug combinations and collection of drug combinations with or without confirmed synergistic effects.

“Omic”-based methods in synergy identification

With advances in bioinformatics, systems biology and molecular biology, different types of high-throughput “omic” data such as genomics, transcriptomics and metabolomics data have emerged as critical information in the biomedical and pharmaceutical fields [6365]. The “omic” data have been widely used in the establishment of interactions among small molecules, genes and diseases and the identification of molecules specifically associated with certain pathological processes [6669]. Owing to its benefits in reflecting biological responses, “omic” data has also made great contributions to drug combination studies. One of the pioneering studies was the CMAP project by Lamb et al. [69]. Based on this project, the researchers illustrated the possibility of rational drug combinations [33, 69].

Considering the aim of “omic” is at collectively quantifying pools of biological molecules [70], specific computational approaches has also been proposed to predict the human-body responses to drug combinations or constructing novel molecular networks based on “omic” data [33, 7174].

The most commonly used “omic” data are the genomic data. To reduce the burden in genomic data acquisition for all possible combinations, the corresponding genomic data for drug combination has been predicted by computational methods based on the data of individual ones. Based on this principle, Wu et al. [71] proposed a computational scheme to predict the gene expression profiles for drug combinations based on gene profiles treated with individual drugs. Then, according to the differential expression changes of genes, they applied a new integer programming model to identify sub-networks affected by single drugs or combinations from the background molecular network which included protein interactions, protein-DNA interactions, and signaling pathways. Finally, the best candidates of combination drugs were determined by taking into account efficacy and side-effect based on the identified sub-networks. The genomic data can also help reconstruct specific molecular network. As mentioned above, the synergy ranker: “DrugComboRanker” [33] utilized the genomic data for drug functional network construction. Moreover, the Petri-net model which is a graph-based mathematical method for the description of distributed systems has also been applied in drug combination analysis based on the genomic data. Jin et al. [72] proposed an enhanced Petri-net (EPN) model based on the gene expression data to simulate drug effects. According to the EPN model, the synergistic effects of pairwise drug combinations were identified when there was at least one molecule that shows an enhanced effect in the combination comparing with the summation effects of two individual drugs; and the researchers illustrated that the synergy mechanism was related with the feedback loops among molecules. Besides, some researchers adopted correlation-based strategies to calculate the similarities of “omic”-based responses to drug interventions [73, 74]. Based on the genome-wide expression profiles, Zhao et al. [73] have revealed the synergistic drug effects by quantitatively profiling cellular responses to drugs. According to the drug similarities based on a chemogenomic profile, Jansen et al. [74] identified antifungal synergies. They found that compound pairs that have correspondingly similar profiles are more likely to be synergistic when compared with randomly chosen compounds.

In addition to genomic data, proteomics is another type of “omic” data for drug synergy study. This type of data can help reconstruct the signaling pathways. As we have mentioned in “Dynamic pathway simulation”, the signaling pathways were reconstructed based on the phosphorylation data of key signaling proteins [4046] by different pathway simulation methods.

Integration of different types of “omic” data can also help explore synergy mechanism. By incorporating phosphoproteomics, transcriptomics and chemical proteomics data, Winter et al. [75] discovered that the underlying synergy mechanism for the combination of danusertib and bosutinib in the treatment of imatinib-resistant Chronic myelogenous leukemia (CML) might be related with the nonobvious off target effects implicated in the Mitogen-activated protein kinases (MAPK) signaling pathway.

With accumulation of “omic” data, the field of systems biology has progressed at an even higher speed. The “omic”-based drug synergy studies are presently mainly conducted on the genomics and proteomics levels. To understand the complex interplay of drug combinations against the disease process, integrative “omic”-based models covering across various “omic” data like genetic, transcriptomic, proteomic and metabolomics data need to be developed.


Rational design of drug combinations with synergistic effects remains a challenging task despite tremendous experimental and clinical efforts. Examples presented in this review illustrate how systems biology approaches are making encouraging contributions to explore synergistic drug combinations. Given the complexity of synergy mechanism and the difficulty in drug combination screening, systematic synergy modeling will become a very important component in the field of drug combinations beyond question. However, there are still great challenges to overcome.

Network or pathway-based methods have already facilitated the prediction of synergistic drug combinations and promoted the understanding of synergy mechanism. Based on varied pathway modeling methods and the increasing “omic” data, pathways for different diseases can be reconstructed, thus providing firm foundations for the exploration of synergy. However, mostly focusing on drug targets to represent drug action, these models for drug synergy mainly depend on the interaction between targets and disease-specific network or pathways. Such methods neglect the fact that synergy can also arise when one drug increases the effects of other drugs by influencing on their certain properties like bioavailability as revealed in [76]. Drug behaviors are not simply determined by targets, other factors like the absorption, distribution, metabolism, and excretion (ADME) processes and off-target activities can also have significant impacts on drug therapeutic outcomes. Even though some classification models based on integrated drug features considering comprehensive chemical and pharmacological properties have overcome this “target-only” issue; however, the current number of effective drug combinations and their obtainable properties will limit the performance. Consequently, more attention should be paid on the acquisition and integration of extensive knowledge to model the drug-organism interplay.

Systems biology models for drug synergy evaluation are still in their infancy, but have shown prominent advantages in exploring synergistic drug combinations. The potential synergy mechanism has also been illustrated based on these methods. To understand drug synergy from a more comprehensive perspective, future efforts should be put on novel approaches which can model the synergistic effects from a more holistic perspective covering across pharmacology-considering both pharmacodynamics and pharmacokinetics elements-and toxicology spaces. With a deep insight into drug synergy, systems biology approaches will serve as a compensative tool for the rational design of effective drug combinations.

However, it should also be noted that systems biology modeling is only an aided tool but not the judge. The drug synergy predictions can provide valuable references for experimental drug combination screening, but they cannot substitute the experimental methods to determine the final statement on drug synergism or antagonism. The systematic synergy modeling methods can provided assistances for drug synergy researches, but the effectiveness and accuracy of these predictive models should still rely on experimental evaluations.



Protein-protein interactions


Food and Drug Administration


Drug Combination Database


Topology Score


Agent Score


Network target-based Identification of Multicomponent Synergy


Topological Score of Drug Synergy


Connectivity Map


Ordinary Differential Equations


Epidermal Growth Factor Receptor


Modular Response Analysis


half maximal Inhibitory Concentration


Ingenuity Pathway Anallysis


Anatomical Therapeutic Chemical


Kyoto Encyclopedia of Genes and Genomes


Enhanced Petri-Net


Chronic myelogenous leukemia


Mitogen-activated protein kinases


Absorption, Distribution, Metabolism, and Excretion


Principal Component Analysis


  1. Bigioni M, Benzo A, Irrissuto C, Lopez G, Curatella B, Maggi CA, et al. Antitumour effect of combination treatment with Sabarubicin (MEN 10755) and cis-platin (DDP) in human lung tumour xenograft. Cancer Chemother Pharmacol. 2008;62:621–9.

    Article  CAS  PubMed  Google Scholar 

  2. Lenz TL, Hilleman DE. Aggrenox: a fixed-dose combination of aspirin and dipyridamole. Ann Pharmacother. 2000;34:1283–90.

    Article  CAS  PubMed  Google Scholar 

  3. Breeze S. Novel HIV-1 treatment Stribild gains regulatory approval. Expert Rev Clin Pharmacol. 2012;5:613.

    Article  PubMed  Google Scholar 

  4. Liu Y, Hu B, Fu C, Chen X. DCDB: drug combination database. Bioinformatics. 2010;26:587–8.

    Article  CAS  PubMed  Google Scholar 

  5. Albers GW, Amarenco P. Combination Therapy With Clopidogrel and Aspirin. Stroke. 2001;32:2948–9.

    Article  CAS  PubMed  Google Scholar 

  6. Bezin J, Pariente A, Lassalle R, Dureau-Pournin C, Abouelfath A, Robinson P, et al. Use of the recommended drug combination for secondary prevention after a first occurrence of acute coronary syndrome in France. Eur J Clin Pharmacol. 2014;70:429–36.

    Article  CAS  PubMed  Google Scholar 

  7. Singh S, Bouzinbi N, Chaturvedi V, Godreuil S, Kremer L. In vitro evaluation of a new drug combination against clinical isolates belonging to the Mycobacterium abscessus complex. Clin Microbiol Infect. 2014;20:O1124–7.

    Article  CAS  PubMed  Google Scholar 

  8. Lee HZ, Miller BW, Kwitkowski VE, Ricci S, DelValle P, Saber H, et al. U.S. Food and drug administration approval: obinutuzumab in combination with chlorambucil for the treatment of previously untreated chronic lymphocytic leukemia. Clin Cancer Res. 2014;20:3902–7.

    Article  CAS  PubMed  Google Scholar 

  9. Patankar NA, Pritchard J, van Grinsven M, Osooly M, Bally MB. Topotecan and doxorubicin combination to treat recurrent ovarian cancer: the influence of drug exposure time and delivery systems to achieve optimum therapeutic activity. Clin Cancer Res. 2013;19:865–77.

    Article  CAS  PubMed  Google Scholar 

  10. Sucher NJ. Searching for synergy in silico, in vitro and in vivo. Synergy. 2014;1:30–43.

    Article  Google Scholar 

  11. Loewe S. Die Mischiarnei. Klin Wochenschr. 1927;6:1077–85.

    Article  Google Scholar 

  12. Loewe S. The problem of synergism and antagonism of combined drugs. Arzneimittelforschung. 1953;3:285–90.

    CAS  PubMed  Google Scholar 

  13. Bliss CI. The toxicity of poison applied jointly. Ann Appl Biol. 1939;28:585–615.

    Article  Google Scholar 

  14. Chou TC. Drug combination studies and their synergy quantification using the Chou-Talalay method. Cancer Res. 2010;70:440–6.

    Article  CAS  PubMed  Google Scholar 

  15. Chou TC. Theoretical basis, experimental design, and computerized simulation of synergism and antagonism in drug combination studies. Pharmacol Rev. 2006;58:621–81.

    Article  CAS  PubMed  Google Scholar 

  16. Tallarida RJ. Quantitative methods for assessing drug synergism. Genes & cancer. 2011;2:1003–8.

    Article  Google Scholar 

  17. Doern CD. When does 2 plus 2 equal 5? A review of antimicrobial synergy testing. J Clin Microbiol. 2014;52:4124–8.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Leung EL, Cao ZW, Jiang ZH, Zhou H, Liu L. Network-based drug discovery by integrating systems biology and computational technologies. Brief Bioinform. 2013;14:491–505.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Csermely P, Agoston V, Pongor S. The efficiency of multi-target drugs: the network approach might help drug design. Trends Pharmacol Sci. 2005;26:178–82.

    Article  CAS  PubMed  Google Scholar 

  20. Lu JJ, Pan W, Hu YJ, Wang YT. Multi-target drugs: the trend of drug research and development. PLoS One. 2012;7:e40262.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Schrattenholz A, Groebe K, Soskic V. Systems biology approaches and tools for analysis of interactomes and multi-target drugs. Methods Mol Biol. 2010;662:29–58.

    Article  CAS  PubMed  Google Scholar 

  22. Feala JD, Cortes J, Duxbury PM, Piermarocchi C, McCulloch AD, Paternostro G. Systems approaches and algorithms for discovery of combinatorial therapies. Wiley Interdiscip Rev Syst Biol Med. 2010;2:181–93.

    Article  PubMed  Google Scholar 

  23. Galizzi JP, Lockhart BP, Bril A. Applying systems biology in drug discovery and development. Drug Metabol Drug Interact. 2013;28:67–78.

    Article  CAS  PubMed  Google Scholar 

  24. Li P, Huang C, Fu YX, Wang JA, Wu ZY, Ru JL, et al. Large-scale exploration and analysis of drug combinations. Bioinformatics. 2015;31:2007–16.

  25. Berg EL. Systems biology in drug discovery and development. Drug Discov Today. 2014;19:113–25.

    Article  CAS  PubMed  Google Scholar 

  26. Zimmermann GR, Lehar J, Keith CT. Multi-target therapeutics: when the whole is greater than the sum of the parts. Drug Discov Today. 2007;12:34–42.

    Article  CAS  PubMed  Google Scholar 

  27. Li Y, Agarwal P. A Pathway-Based View of Human Diseases and Disease Relationships. Plos ONE. 2009;4:e4346.

  28. Li S, Zhang B, Zhang N. Network target for screening synergistic drug combinations with application to traditional Chinese medicine. BMC Syst Biol. 2011;5 Suppl 1:S10.

    Article  Google Scholar 

  29. Vitali F, Mulas F, Marini P, Bellazzi R. Network-based target ranking for polypharmacological therapies. J Biomed Inform. 2013;46:876–81.

    Article  PubMed  Google Scholar 

  30. Huang L, Li F, Sheng J, Xia X, Ma J, Zhan M, et al. DrugComboRanker: drug combination discovery based on target network analysis. Bioinformatics. 2014;30:i228–36.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Wang YY, Xu KJ, Song J, Zhao XM. Exploring drug combinations in genetic interaction network. BMC bioinformatics. 2012;13 Suppl 7:S7.

    Article  CAS  Google Scholar 

  32. Ren G, Liu Z. NetCAD: a network analysis tool for coronary artery disease-associated PPI network. Bioinformatics. 2013;29:279–80.

    Article  CAS  PubMed  Google Scholar 

  33. Wan FC, Cui YP, Wu JT, Wang JM, -Z Liu Q, Gao ZL. The PPI network and cluster ONE analysis to explain the mechanism of bladder cancer. Eur Rev Med Pharmacol Sci. 2013;17:618–23.

    PubMed  Google Scholar 

  34. Wu B, Xie J, Du Z, Wu J, Zhang P, Xu L, et al. PPI network analysis of mRNA expression profile of ezrin knockdown in esophageal squamous cell carcinoma. Biomed Res Int. 2014;2014:651954.

    PubMed Central  PubMed  Google Scholar 

  35. Viswanathan GA, Seto J, Patil S, Nudelman G, Sealfon SC. Getting started in biological pathway construction and analysis. PLoS Comput Biol. 2008;4, e16.

    Article  PubMed Central  PubMed  Google Scholar 

  36. Jia J, Zhu F, Ma X, Cao Z, Li Y, Chen YZ. Mechanisms of drug combinations: interaction and network perspectives. Nat Rev Drug Discov. 2009;8:111–28.

    Article  CAS  PubMed  Google Scholar 

  37. Polynikis A, Hogan SJ, di Bernardo M. Comparing different ODE modelling approaches for gene regulatory networks. J Theor Biol. 2009;261:511–30.

    Article  CAS  PubMed  Google Scholar 

  38. Elias J, Dimitrio L, Clairambault J, Natalini R. The dynamics of p53 in single cells: physiologically based ODE and reaction–diffusion PDE models. Phys Biol. 2014;11:045001.

    Article  PubMed  Google Scholar 

  39. Khan FM, Schmitz U, Nikolov S, Engelmann D, Putzer BM, Wolkenhauer O, et al. Hybrid modeling of the crosstalk between signaling and transcriptional networks using ordinary differential equations and multi-valued logic. Biochim Biophys Acta. 1844;2014:289–98.

    Google Scholar 

  40. Araujo RP, Petricoin EF, Liotta LA. A mathematical model of combination therapy using the EGFR signaling network. Biosystems. 2005;80:57–69.

    Article  CAS  PubMed  Google Scholar 

  41. Nelander S, Wang W, Nilsson B, She QB, Pratilas C, Rosen N, et al. Models from experiments: combinatorial drug perturbations of cancer cells. Mol Syst Biol. 2008;4:216.

    Article  PubMed Central  PubMed  Google Scholar 

  42. Sun X, Bao J, Nelson KC, Li KC, Kulik G, Zhou X. Systems Modeling of Anti-apoptotic Pathways in Prostate Cancer: Psychological Stress Triggers a Synergism Pattern Switch in Drug Combination Therapy. PLoS Comput Biol. 2013;9:e1003358.

    Article  PubMed Central  PubMed  Google Scholar 

  43. Miller ML, Molinelli EJ, Nair JS, Sheikh T, Samy R, Jing X, et al. Drug synergy screen and network modeling in dedifferentiated liposarcoma identifies CDK4 and IGF1R as synergistic drug targets. Sci signal. 2013;6:ra85.

    PubMed Central  PubMed  Google Scholar 

  44. Klinger B, Sieber A, Fritsche-Guenther R, Witzel F, Berry L, Schumacher D, et al. Network quantification of EGFR signaling unveils potential for targeted combination therapy. Mol Syst Biol. 2013;9:673.

    Article  PubMed Central  PubMed  Google Scholar 

  45. Yan H, Zhang B, Li S, Zhao Q. A formal model for analyzing drug combination effects and its application in TNF-alpha-induced NFkappaB pathway. BMC Syst Biol. 2010;4:50.

    Article  PubMed Central  PubMed  Google Scholar 

  46. Azmi AS, Banerjee S, Ali S, Wang Z, Bao B, Beck FWJ, et al. Network Modeling of MDM2 Inhibitor-Oxaliplatin Combination Reveals Biological Synergy in wt-p53 solid tumors. Oncotarget. 2011;2:15.

    Google Scholar 

  47. Li Z, Zhou H, Lu Y, Colatsky T. A Critical Role for Immune System Response in Mediating Anti-influenza Drug Synergies Assessed by Mechanistic Modeling. CPT Pharmacometrics Syst Pharmacol. 2014;3, e135.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. Milo R, Shen-Orr S, Itzkovitz S, Kashtan N, Chklovskii D, Alon U. Network motifs: simple building blocks of complex networks. Science. 2002;298:824–7.

    Article  CAS  PubMed  Google Scholar 

  49. Wong E, Baur B, Quader S, Huang CH. Biological network motif detection: principles and practice. Brief Bioinform. 2012;13:202–15.

    Article  PubMed Central  PubMed  Google Scholar 

  50. Han Z, Vondriska TM, Yang L, Robb MacLellan W, Weiss JN, Qu Z. Signal transduction network motifs and biological memory. J Theor Biol. 2007;246:755–61.

    Article  PubMed Central  PubMed  Google Scholar 

  51. Tyson JJ, Novak B. Functional motifs in biochemical reaction networks. Annu Rev Phys Chem. 2010;61:219–40.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  52. Fioravanti F, Helmer-Citterich M, Nardelli E. Modeling gene regulatory network motifs using Statecharts. BMC bioinformatics. 2012;13 Suppl 4:S20.

    Article  PubMed Central  PubMed  Google Scholar 

  53. Zhang Y, Smolen P, Baxter DA, Byrne JH. Computational analyses of synergism in small molecular network motifs. PLoS Comput Biol. 2014;10, e1003524.

    Article  PubMed Central  PubMed  Google Scholar 

  54. Yin N, Ma W, Pei J, Ouyang Q, Tang C, Lai L. Synergistic and Antagonistic Drug Combinations Depend on Network Topology. PLoS ONE. 2014;9, e93960.

    Article  PubMed Central  PubMed  Google Scholar 

  55. Axelrod M, Gordon VL, Conaway M, Tarcsafalvi A, Neitzke DJ, Gioeli D, et al. Combinatorial drug screening identifies compensatory pathway interactions and adaptive resistance mechanisms. Oncotarget. 2013;4:622–35.

    PubMed Central  PubMed  Google Scholar 

  56. Dent P, Curiel DT, Fisher PB, Grant S. Synergistic combinations of signaling pathway inhibitors: Mechanisms for improved cancer therapy. Drug Resist Updat. 2009;12:65–73.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  57. Azmi AS, Wang Z, Philip PA, Mohammad RM, Sarkar FH. Proof of concept: network and systems biology approaches aid in the discovery of potent anticancer drug combinations. Mol Cancer Ther.2010;9:3137–44.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  58. Reddy AS, Zhang S. Polypharmacology: drug discovery for the future. Expert Rev Clin Pharmacol. 2013;6:41–7.

    Article  CAS  PubMed  Google Scholar 

  59. Parikh JR, Klinger B, Xia Y, Marto JA, Bluthgen N. Discovering causal signaling pathways through gene-expression patterns. Nucleic Acids Res. 2010;38:W109–17.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  60. Janes KA, Lauffenburger DA. Models of signalling networks - what cell biologists can gain from them and give to them. J Cell Sci. 2013;126:1913–21.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  61. Zhao XM, Iskar M, Zeller G, Kuhn M, van Noort V, Bork P. Prediction of drug combinations by integrating molecular and pharmacological data. PLoS Comput Biol. 2011;7, e1002323.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  62. Chen L, Li BQ, Zheng MY, Zhang J, Feng KY, Cai YD. Prediction of effective drug combinations by chemical interaction, protein interaction and target enrichment of KEGG pathways. Biomed Res Int. 2013;2013:723780.

    PubMed Central  PubMed  Google Scholar 

  63. Cui Y, Paules RS. Use of transcriptomics in understanding mechanisms of drug-induced toxicity. Pharmacogenomics. 2010;11:573–85.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  64. Zhao S, Iyengar R. Systems pharmacology: network analysis to identify multiscale mechanisms of drug action. Annu Rev Pharmacol Toxicol. 2012;52:505–21.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  65. Kaur G, Rajput B. Comparative analysis of the omics technologies used to study antimonial, amphotericin B, and pentamidine resistance in leishmania. J Parasitol Res. 2014;2014:726328.

    Article  PubMed Central  PubMed  Google Scholar 

  66. Lee KJ, Yin W, Arafat D, Tang Y, Uppal K, Tran V, et al. Comparative transcriptomics and metabolomics in a rhesus macaque drug administration study. Front Cell Dev Biol. 2014;2:54.

    Article  PubMed Central  PubMed  Google Scholar 

  67. Stratton MR. Exploring the genomes of cancer cells: progress and promise. Science. 2011;331:1553–8.

    Article  CAS  PubMed  Google Scholar 

  68. Tieri P, Zhou X, Zhu L, Nardini C. Multi-omic landscape of rheumatoid arthritis: re-evaluation of drug adverse effects. Front Cell Dev Biol. 2014;2:59.

    Article  PubMed Central  PubMed  Google Scholar 

  69. Lamb J, Crawford ED, Peck D, Modell JW, Blat IC, Wrobel MJ, et al. The Connectivity Map: using gene-expression signatures to connect small molecules, genes, and disease. Science. 2006;313:1929–35.

    Article  CAS  PubMed  Google Scholar 

  70. Coughlin SS. Toward a road map for global -omics: a primer on -omic technologies. AM J Epidemiol. 2004;180:1188–95.

    Article  Google Scholar 

  71. Wu Z, Zhao XM, Chen L. A systems biology approach to identify effective cocktail drugs. BMC Syst Biol. 2010;4 Suppl 2:S7.

    Article  PubMed Central  PubMed  Google Scholar 

  72. Jin G, Zhao H, Zhou X, Wong ST. An enhanced Petri-net model to predict synergistic effects of pairwise drug combinations from gene microarray data. Bioinformatics. 2011;27:i310–6.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  73. Zhao J, Zhang XS, Zhang S. Predicting cooperative drug effects through the quantitative cellular profiling of response to individual drugs. CPT Pharmacometrics Syst Pharmacol. 2014;3, e102.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  74. Jansen G, Lee AY, Epp E, Fredette A, Surprenant J, Harcus D, et al. Chemogenomic profiling predicts antifungal synergies. Mol Syst Biol. 2009;5:338.

    Article  PubMed Central  PubMed  Google Scholar 

  75. Winter GE, Rix U, Carlson SM, Gleixner KV, Grebien F, Gridling M, et al. Systems-pharmacology dissection of a drug synergy in imatinib-resistant CML. Nat Chem Biol. 2012;8:905–12.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  76. Cokol M, Chua HN, Tasan M, Mutlu B, Weinstein ZB, Suzuki Y, et al. Systematic exploration of synergistic drug pairs. Mol Syst Biol. 2011;7:544.

    Article  PubMed Central  PubMed  Google Scholar 

Download references


The writing of the review was supported by the National Natural Science Foundation of China (Grant No. 81303152), and the State Key Program of National Natural Science Foundation of China (Grant No. 81330086).

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Hongjun Yang or Peng Lu.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Di Chen and Xi Liu conceived and wrote manuscript; Yiping Yang, Hongjun Yang and Peng Lu drafted and revised the manuscript. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided 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 Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, D., Liu, X., Yang, Y. et al. Systematic synergy modeling: understanding drug synergy from a systems biology perspective. BMC Syst Biol 9, 56 (2015).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: