 Research Article
 Open Access
 Published:
Drug combinatorics and side effect estimation on the signed human drugtarget network
BMC Systems Biology volume 10, Article number: 74 (2016)
Abstract
Background
The mode of action of a drug on its targets can often be classified as being positive (activator, potentiator, agonist, etc.) or negative (inhibitor, blocker, antagonist, etc.). The signed edges of a drugtarget network can be used to investigate the combined mechanisms of action of multiple drugs on the ensemble of common targets.
Results
In this paper it is shown that for the signed human drugtarget network the majority of drug pairs tend to have synergistic effects on the common targets, i.e., drug pairs tend to have modes of action with the same sign on most of the shared targets, especially for the principal pharmacological targets of a drug. Methods are proposed to compute this synergism, as well as to estimate the influence of the drugs on the side effect of another drug.
Conclusions
Enriching a drugtarget network with information of functional nature like the sign of the interactions allows to explore in a systematic way a series of network properties of key importance in the context of computational drug combinatorics.
Background
In drug discovery research, moving beyond the “magic bullet” (onedrug, onetarget) paradigm is a trend that has become quite popular in recent times, motivated by the discovery of the value of polypharmacology in treating complex diseases [1–3], and by the promises of network pharmacology [4–8], as a “system” way to understand the effects (and side effects) of drugs on an organism. A network viewpoint can help in understanding the effect of a drug on system properties such as robustness, resilience and redundancy, and can be used for systemdriven drug discovery [9]. The approach is particularly meaningful to deal with the multifactor nature of many complex diseases, like cancer, asthma, diabetes, neurodegenerative disorders and cardiovascular diseases.
Network representations of drugtarget interactions have already been used for several tasks, like to extrapolate information of functional nature on the action of the drugs, to predict novel putative drugtarget interactions, and to provide strategies for an efficient use of multidrug therapies, see [10, 11] for an overview. For instance the topological structure and organization of a drugtarget network obtained from FDAapproved drugs is investigated in [12], as well as the clustering of its drugs according to the functional categories of the Anatomical Therapeutic Chemical classification. Chemical and genomic information is used in [13] to build a drugtarget network for various classes of human targets, like enzymes, ion channels and receptors, using supervised learning. Systematic attempts to identify new drugtarget interactions are very frequent in the literature, based on sequence [14], structure [15], pathwaypathway interactions [16], but also on chemical and phenotypical similarity among drugs and sets of ligands [17, 18], often leading to a large number of lowaffinity interactions of limited significance. Methods for predicting new potential drugtarget interactions for known drugs, using drugbased, targetbased and networkbased similarity scores are investigated in [19], inspired by algorithms used in the field of recommendation systems, or by random walking the network itself [20]. Drug repositioning is also one of those problems that can be investigated through drugtarget networks. Machine learning techniques aimed at extending the target space of already approved drugs are reviewed in [21], while in [22] it is shown how to use constraint based computational methods for metabolic drug repurposing. Another line of research deals with predicting drug combinatorics via gene profiling and gene network reconstruction [23]. As surveyed in [10], the number of datasets and tools specifically dedicated to the analysis of drugtarget networks has grown rapidly in the last years.
When investigating a given drugtarget network (for instance reconstructed from a database such as DrugBank [24], which is also our database of choice), most computational methods rely principally on the topological information that can be obtained from the bipartite drugtarget graph and on the ontology of its constituent components. In this work we aim to add another element of functional nature in the networkbased investigation of drugtarget interactions, namely the information on the mode of action of the interactions. When browsing DrugBank, many are the possible mechanisms of action of a drug on its targets: it can activate or inhibit the target, it can act as an agonist or an antagonist, as a potentiator or as a blocker, as an inducer or as a suppressor, and so on. Although qualitatively different and applicable to different categories of targets (proteins, macromolecules, nucleic acids, small molecules, etc.), these modes of action can be reasonably classified as positive or negative. If these and a few more modes of action admit such a signed classification, several more, such as for instance “modulators”, “binder”, “cleavage” (a complete list of categories will be given below) are instead impossible to classify with a sign, and hence cannot be included in the analysis we are proposing in this paper. On the drugtarget network, characterizing the modes of action as positive or negative corresponds to associating a sign to the edges of the bipartite graph. Signed graphs have been frequently used in Systems Biology [25–28], and we can draw inspiration from this literature to formulate and solve problems which are meaningful and insightful also for our signed drugtarget network.
One such problem, important in the context of multicomponent therapies, is to understand the joint effect of two or more drugs acting on the same targets [29]. To illustrate the point, let us look at two drugs having a common target: if their modes of action have the same sign then it is likely that the combined effect of the two drugs on the target is reinforced. On the contrary, if the modes of action have opposite signs then it is more plausible to assume that the two drugs tend to compensate each other’s action, and hence the overall effect on the target tends to be mitigated. In this paper two actions that have a common sign on a target (regardless of what that sign is) are called coherent. They are called incoherent when they have opposite signs. Understanding to what extent drug pairs act coherently or incoherently on common targets is an important aspect for instance of drug combinatorics. It is shown in the paper that coherent drug pairs are much more frequent than incoherent, and that the same observation is true when we count drug pairs acting simultaneously on pairs of targets (i.e., we look at length4 undirected cycles). These basic examples of coherent/incoherent motif counting hint at the signed human drugtarget network having a topological and functional (for what can be deduced from the edge signs) organization which is far from random. Such an organization can be investigated in a more systematic way if, rather than common actions on one or two targets, we consider the combined action of drug pairs on the ensemble of common targets. We show in the paper that the distribution of drug pairs acting simultaneously on multiple targets is significantly skewed both towards actions that are coherent on many targets but also towards actions that are incoherent on multiple targets. Neither of these two tails is present in null models. Their meaning is that current drugs tend to interact largely on overlapping targets, often exerting a similar action on all of them (hence leading to coherent drug pairs on multiple targets) but sometimes also exerting antithetical actions on all of them (hence leading to incoherent drug pairs on multiple targets). This is related to the observation that most drugs have a “monochromatic” mode of action (i.e. with same sign) on multiple targets.
By attaching signs to the mechanisms of action, we are also able to quantify the amount of synergism (i.e., when coherence prevails over incoherence) in a drug pair, and to classify all drug pairs accordingly. In the paper we do this systematically for all drugtarget interactions that are classified as pharmacological in DrugBank, which we take as principal actions (ontarget actions) of a drug. The vast majority of drug pairs has a beneficial synergistic action on their common pharmacological drug targets.
Failure of a drug to act properly in clinical trials is often due to unexpected side effects, in the form of drug activity on targets other than the principal targets for which the pharmacological action is intended. Methods to predict in a systematic way such offtarget interactions have attracted a considerable attention in recent years [17, 30–33]. In some of these works the side effect is defined phenotypically, in some others based on the offtargets of a drug. If we consider as offtargets of a drug the DrugBank interactions that are not classifiable as pharmacological, then exploiting the signs of the interactions we can try to quantify the action of a drug on the offtargets of a second drug, and in particular seek for drug pairs likely to have a mutualistic benefit on each other’s side effect, for instance drug pairs behaving synergistically on their common principal (pharmacological) targets but having incoherent actions on their common offtargets. An exhaustive classification of all these pairs is carried out for our signed drugtarget network.
Results
Consider the human drugtarget network reconstructed from the DrugBank database (see Methods and Table 1), and associate to its edges a sign according to the procedure described in Methods and in Table 2. Table 3 reports the resulting sign distribution. The fraction of edges which cannot be assigned a sign is around 63 %. The fraction decreases drastically if we restrict to actions which can be classified as pharmacological (as defined in the Methods), to less than 11 %.
Looking at the drug degree distribution of Fig. 1 a, one can see that the drugs that lack most edge signs are those with higher connectivity, while the situation is much improved for drugs with lower connectivity. Comparing the drug and target connectivity analysis of Fig. 1, the latter has a significant difference in the fact that also targets with high connectivity tend to have a certain fraction (even higher that 50 %) of edges having a known sign. In the drug connectivity histogram, on the other hand, when a drug has signed edges, then they are almost all of the same sign, see inset of Fig. 1 a.
Additional file 1: Figure S2 shows the equivalent drug and target connectivity analysis restricted to edges representing pharmacological actions. In both cases the connectivity decreases considerably when we compare it with Fig. 1, and the already mentioned fact that most pharmacological actions have a known sign (see Table 3) allows us to have a sharper picture of the sign distributions. For both drugs and targets, in fact, the sign distribution tends to be “monochromatic”, i.e., most edges adjacent to a drug or to a target tend to be positive or negative but rarely both, see inset plots in Additional file 1: Figure S2.
When restricting the drugtarget network to signed edges, we obtain a bipartite signed graph involving 1315 drugs and 820 targets, containing 1417 positive and 2711 negative edges, see Table 3. The drugtarget network described in Table 3 has one very large connected component involving 81 % of the drugs and 71 % of the targets, plus a number of other smaller connected components, see Additional file 1: Figure S1.
Combinatorics of signed drugs. The classification of drug actions into positive and negative modes of action allows us to characterize the effect of multiple drugs acting on the same target. Two drugs sharing the same target tend to reinforce their effect on the target if their modes of action have the same sign, while they tend to mitigate the overall effect if the signs are different. The 3 possible sign combinations of a drug pair on a common target, (+, +), (−, −), and (+, −). are shown in Fig. 2 a. The two combinations (+, +) and (−, −) will be referred to as coherent, as the action of a drug is reinforced by the presence of a second drug. In the remaining combination (+, −), instead, the presence of a second drug tends to counteract (and hence in general reduce) the action of the first drug. We will call this configuration incoherent.
Some statistics on pairs of drugs sharing at least a common target are shown in Fig. 2 b. Summing up the data of the top left subpanel of Fig. 2 b, the total number of drug pairs incident to at least a common target is 25461. The number of drug pairs having exactly 1 target in common is 15947 (highest bar in top left subpanel of Fig. 2 b). The statistics for the three pairs of edge sign combinations (+, +), (−, −) and (+, −) are reported in the remaining 3 subpanels of Fig. 2 b. As can be seen in Fig. 2 b, incoherent drug actions are not rare in the human drugtarget network. The total number of drug pairs acting incoherently on at least one target is 6431. These drug pairs act incoherently on a total of 10594 targets (sum of all bars in the bottom right subpanel of Fig. 2 b). The corresponding numbers for coherent actions are however much higher: 20017 drug pairs act coherently on at least one target, for a total number of coherent actions of 37439.
The shortcoming of the analysis carried out so far is that even after restricting to signed actions, 57 % of the 1315 drugs have 2 or more targets, hence when we look at multiple drugs applied simultaneously it is necessary to investigate their effect on all of their common targets. As can be seen in Fig. 2 b (top left subpanel), around 5000 drug pairs share 2 targets, and more than 2000 have 3 targets in common, while other 2400 pairs have more than 3 targets in common. The simplest possible approach to tackle this more complex problem consists in looking at all cycles of length 4 formed by two drugs having at least 2 targets in common. Cycles are here obviously intended as undirected, i.e., arrows on the edges are dropped. Three qualitatively different classes of length4 cycles can be identified, which we call (fully) coherent, mixed and (fully) incoherent. They are shown in Fig. 3. Coherent cycles are those for which both drugs act with the same sign on each of the targets, hence reinforcing each other’s action on both targets. In incoherent cycles, instead, the signs of the drug actions are conflicting on both targets, leading to mitigation of the effect on both targets. Mixed cycles occur when the action of a drug pair has the same sign on one of the targets but conflicting signs on the other target (i.e., coherent on one target but incoherent on the other). As can be seen on Table 4, the fraction of coherent length4 cycles is around 81 % of the total, while that of incoherent length4 cycles is 13.5 %. Mixed cycles are the most rare (around 5.5 %). Notice how fully coherent and fully incoherent length4 cycles are positive, i.e., they have an even number of negative edges, while mixed cycles are negative (odd number of negative edges). Hence, overall, the fraction of positive length4 cycles is around 94.5 % of the total, see Table 4.
Connectivity analysis of drug pairs: distribution of coherent/incoherent actions. With a total of 48033 targets being acted upon by at least two drugs, a natural question to ask is if such pairwise connectivity is high or low with respect to a null model having the same drug edge distribution. Comparing with a null model, obtained maintaining the same edge distribution (sign included) at the drug side, but reassigning randomly the edges to the targets, then the drugtarget topological properties one obtains are drastically different. On a null model, the number of drug pairs sharing at least a target reduces from 25461 to 10126.6 (average over 100 realizations of the null model) and so does the total number of drug pairs incident to a common target, from 48033 to 10388.7. We deduce that the signed drugtarget network is highly organized, and in particular highly redundant in its coverage of the targets. For instance, the original target connectivity is such that only near half (416) of the 820 targets are hit by more than one drug, while in a null model this number is around 95 %. However, the target edge distribution for the null model completely lacks the highly connected nodes which can instead be seen in panel (b) of Fig. 1. As can be seen in Fig. 4, the consequence is that in null models a pair of drugs very seldom exceeds a total of 3 common targets, counting both coherent and incoherent actions (green bars). In the true drugtarget network, instead, edges tend to concentrate on fewer targets, leading to abundance of targets shared by more than two drugs (in Fig. 4 green bars for the null model should be compare to all the other bars taken together, representing the original network).
Also the sign distribution in the true drugtarget network is highly nonrandom. If 77.9 % of the drug pairs hitting a common target have a coherent action, such number is drastically less in our null models, around 54.8 %, meaning that in the real drugtarget network not only known drugs tend to hit always the same targets, but they tend to do so with the same sign, leading to an overabundance of coherent joint actions. In spite of the limited frequency, also the fraction of incoherent drug pair actions is far from being negligible. As can be seen on Fig. 4, the histogram with counts of coherent and incoherent actions of drug pairs is highly skewed and has two tails: just like many drug pairs act coherently on many targets, there exists a considerable number of drug pairs having an incoherent action simultaneously on multiple targets, for instance pairs of drugs in which one of the two activates all its targets while the other inhibit them. The blue bars in Fig. 4 correspond to pairs whose incoherent action is statistically significant, given the total amount of pairs and its coherent/incoherent partition (cumulative binomial test, see “Methods”). Since the number of coherent action pairs is much higher, statistical significance for them requires a pair to hit coherently a larger number of targets (red bars).
Synergistic/compensatory pharmacological effect of drug pairs. Most drugs are designed to produce a specific action on selected targets, here denoted pharmacological targets. Other targets of the same drug (here called offtargets, see Methods) are often present but they are normally undesired. We refer to these offtargets as the side effect of the drug.
Let t _{ i } be the number of targets associated to the ith drug for which the sign of the action is available and p _{ i }, p _{ i }≤t _{ i }, the number of signed pharmacological targets of drug i. Let further be c _{ ij } the number of common (signed) targets of the drug pair (i, j), and p _{ ij }, p _{ ij }≤c _{ ij }, the number of (signed) pharmacological targets shared by (i, j). Denote s _{ i }=t _{ i }−p _{ i } the number of (signed) offtargets of drug i and s _{ ij }=c _{ ij }−p _{ ij } the number of (signed) offtargets in common between ith and jth drugs. To specify the sign of the action we will use upper indices: \( p_{ij}^{++}\), \( p_{ij}^{}\) and \( p_{ij}^{+}\) (respectively \( s_{ij}^{++}\), \( s_{ij}^{}\) and \( s_{ij}^{+}\)). In this paper we are assuming that the coherent pairs \( p_{ij}^{++} \) and \( p_{ij}^{} \) enhance the action of each single drug on the common pharmacological targets (i.e., act synergistically), while in an incoherent pair \( p_{ij}^{+}\) the actions of the two drugs tend to counteract each other (i.e., the action is compensatory). Consider the case p _{ ij }>0 and s _{ ij }=0, i.e., the two drugs share only common pharmacological targets. Additional file 1: Figure S4 tells us that for pharmacological actions the distribution of coherent/incoherent interaction pairs has a longer tail towards the coherent end, but does not tell us what fraction of pharmacological targets of a drug i receives a benefit from the presence of a second drug j sharing a certain number of pharmacological targets with it. In order to quantify this, we use the synergistic score coefficient a _{ i }(j) defined in Methods. This coefficient is by construction between −1 and 1. a _{ i }(j)>0 means that more coherent than incoherent actions are established on the pharmacological targets of drug i by the presence of a second drug j. Since a _{ i }(j)>0 if and only if a _{ j }(i)>0, this synergistic effect is always mutual, although a _{ i }(j)≠a _{ j }(i) when the two drugs have a different number of pharmacological targets. Negative synergistic score means that for a drug pair the incoherent actions are more numerous than the coherent ones. Clearly a _{ i }(j)<0 if and only if a _{ j }(i)<0.
Of the 14380 drug pairs with common pharmacological targets, 10611 have p _{ ij }>0 and s _{ ij }=0, hence for them we can compute a _{ i }(j). The synergistic score coefficient is positive in 81.6 % of the cases. In 45.2 % of the cases it is a _{ i }(j)≥0.5 and a _{ j }(i)≥0.5, i.e., a significant mutual synergistic score is produced, see Additional file 2, while in 29.8 % of the cases a _{ i }(j)≥0.5 but not a _{ j }(i), see Additional file 3. The 70 cases in which a _{ i }(j)=a _{ j }(i)=0 are listed in the Additional file 4, and the 1882 cases in which both a _{ i }(j) and a _{ j }(i) are <0 are in the Additional file 5. The histogram of the synergistic score coefficients is shown in Fig. 5.
Side effect improvement/aggravation through drug combinations. When two drugs share pharmacological targets but also offtargets, it becomes interesting to understand if it is possible to combine drugs so as to reduce their side effect. The principle that we follow is to look for drug pairs having a positive synergism on common pharmacological targets and incoherent actions on common offtargets. As incoherent actions correspond to opposite edge signs, they tend to compensate each other, hence the side effect tends to be mitigated.
Let us consider drug pairs such that p _{ ij }>0 and s _{ ij }>0, i.e., having some pharmacological target and some offtarget in common. In particular, we look for drug combinations in which

1.
the actions on the pharmacological targets are predominantly coherent: \( p_{ij}^{++} + p_{ij}^{} > p_{ij}^{+} \);

2.
the actions on offtargets tend to cancel each other: \( s_{ij}^{+}> s_{ij}^{++} +s_{ij}^{} \).
To quantify the benefit of adding a second drug j on the side effect of drug i, we consider the side effect score coefficient b _{ i }(j) (see Methods). By construction −1≤b _{ i }(j)≤1, with b _{ i }(j)<0 meaning an aggravating effect of the drug j on the side effect of i, b _{ i }(j)=0 meaning neutral effect in the side effect (\( s_{ij}^{+} = s_{ij}^{++}+s_{ij}^{}\), see Fig. 6 a), and b _{ i }(j)>0 meaning an improvement of the side effect of drug i due to drug j, see Fig. 6 b. In particular, b _{ i }(j)=1 means that the drug j hits all offtargets of drug i with an action which is opposite in sign to that of drug i. Notice that by construction b _{ i }(j) and b _{ j }(i) must have the same sign, meaning that for drug pairs with a positive side effect score the benefits are always mutualistic. Screening the 14380 drug pairs having pharmacological targets in common, we obtain 325 pairs. Among these, 205 pairs are beneficial i.e., have b _{ i }(j)>0 and b _{ j }(i)>0; 39 are neutral, i.e., have b _{ i }(j)=b _{ j }(i)=0 and 81 are aggravating, i.e., b _{ i }(j)<0 and b _{ j }(i)<0. The distribution of their side effect score coefficients b _{ i }(j) is shown in Fig. 6 c. For a significant fraction of drug pairs (56 out of 205), both coefficients b _{ i }(j) and b _{ j }(i) are ≥0.5, meaning that the drugs i and j have a significant reciprocal beneficial effect on each other’s side effect. These cases are reported in the Additional file 6. Several more cases (81) are obtained in which b _{ i }(j)≥0.5 for one of the 2 drugs, although not for both. These cases with “unilateral” side effect improvement are listed in the Additional file 7. They often correspond to cases in which one drug has a large side effect (in terms of number of offtargets) and the other drug a small side effect, see Fig. 6 b. Hence when the latter drug is added to the former it can only improve on a few of its many offtargets. The drug pairs giving b _{ i }(j)=b _{ j }(i)=0 and those aggravating the side effect are reported in Additional files 8 and 9.
Discussion
A drugtarget network contains valuable information for researchers interested in network pharmacology and drug combinatorics. This has to do not only with the biochemical classification of the drug compounds and with the ontological/therapeutical classification of the molecular targets, but also with the functional classification of the modes of action of the drugs. In order to explore systematically this aspect, we have to introduce a coarsegrained classification of the mechanisms of action. Such a “binary” classification covers a significant fraction of the known categories for drugtarget modes of action.
The fact that many human drugtarget interactions cannot be classified in terms of action signs is certainly a limitation of the present study and of the approach in general. The fact that the number of such unsigned edges decreases so drastically when we zoom on pharmacological targets (passing from 63 to 11 %) suggests that our results should give reliable predictions on the coherence/incoherence of drug pairs on the “primary” pharmacological targets (and hence on their potential synergies), but could have a limited predictability power on the effective side effect score, when many other unsigned interactions are shared by a pair of drugs alongside the coherent/incoherent offtargets we consider here.
Having so few unsigned pharmacological drugtarget interactions also suggests that the modes of action that can be characterized with signs are more “valuable” in terms of describing the therapeutical effects of a drug. Hence also the investigation of the amount of coherence/incoherence encoded in this subset of signed pharmacological interactions is important. An analysis like the one carried out in Figs. 1, 2 b and 4 is repeated in the Additional file 1 for the subnetwork of signed pharmacological targets (see Additional file 1: Figures S2, S3 and S4). It turns out that some of the properties mentioned for the entire signed drugtarget network still hold for the signed subnetwork of pharmacological actions. For instance, the long tail of coherent action pairs observed in the distribution of Fig. 4 is still visible in Additional file 1: Figure S4, much more than the corresponding incoherent tail, in accordance with the high synergistic score we found.
Moving beyond drug pairs (to triplets, quadruplets, etc.) is a more challenging problem. A possible approach consists in computing the socalled level of structural balance [34], i.e., the amount of “disorder” that characterizes the drugtarget network as a whole, intending with that the amount of contradictory “orders” that multiple drugs send to their common targets [25, 26]. Such information, which is particularly important in a network perspective, is carried out in the Additional file 1. Its limitation is that it does not allow to distinguish between coherent and incoherent edge pairs, but only between positive (undirected) cycles and negative (undirected) cycles. Only the latter contribute to the disorder of the network.
When classifying the common targets of drug pairs according to principal targets and offtargets, it can be observed that in the vast majority of cases in which drug pairs share common pharmacological targets they have no common offtarget, meaning that the synergistic score coefficient describes the entire overlap of the pair. The number of drug pairs with significantly large synergistic score is fairly big (4799 in Additional file 2 and 3163 on Additional file 3). For many of these pairs sharing one or more common pharmacological targets, experimental evidence of synergistic action is available. For instance the DNA antimetabolites Gemcitabine and Fludarabine share an inhibitory action on ribonuceotide reductase, hence we classify them as synergistic. They are known to act synergistically as anticancer agents in Acute Myeloid Leukemia [35]. For other pairs, there is no documented improvement by the simultaneous application of the two drugs. For instance, Panitumumab and Cetuximab are both monoclonal antibodies targeting the epidermal growth factor receptor, and used for treating EGFRexpressing metastatic colorectal cancers. Although they are said to have similar activity [36], they are also though to differ in their isotype and possibly in their mechanisms of action (Panitumumab is known to be inefficient for certain type of mutations [37]). In our metric the synergistic scores of the pair are both 1, just based on the common target. In other cases the improvement, as measured by the synergistic score, is positive but low because of the limited overlap among the pharmacological targets. One example among many is the pair Propofol and Sevoflurane. Both have the GABA_{A} receptor among their pharmacological targets. It is known [29] that the mechanism of action is “positive allosteric modulation” for both (here “ + ”), although their actions occur on different sites of the same target. In fact, their joint action is classified as additive in [29]. Given that Propofol has two more pharmacological targets and Sevoflurane has five more (plus, for both, several offtargets, none in common), their a _{ i }(j) coefficients are positive but low.
In oder to make a more systematic validation of the synergistic/compensatory effects described in this paper, we downloaded the DCDB database [38] (http://www.cls.zju.edu.cn/dcdb/), containing ∼1360 experimental combinations of two or more drugs. Around ∼1000 of these combinations are for human targets, and of these 118 overlap with the drug pairs sharing at least one target considered in this study (in the remaining experiments the drugs do not hit the same target, hence do not overlap with our subset of drug pairs). The drug combinations of DCDB are labeled as “efficacious” if experimental evidence is available of improved benefits (80 out of 118) “nonefficacious” for unsuccessful usages (28 out of 118), and “need further study” when unclear (10 out of 118). The majority of these 118 drug combinations are classified as (+, +) actions (31 cases) or (−, −) actions (72 cases). The vast majority of drug pairs sharing a common target for which experimental evidence is available confirms the synergistic effect we are suggesting in this paper. In particular this is true for the (+, +) pairs (24 out of 31), less for the (−, −) pairs (44 out of 72). The complete list is given in Additional file 10.
For what concerns the drug pairs influencing each other’s side effect, their number is remarkably low. This is probably due to the fact that often times when s _{ ij }>0 the synergistic score a _{ i }(j) is negative, i.e., the pair has a majority of compensatory effects on its common pharmacological targets. These pairs are not considered in our analysis (they belong to the incoherent tail of the histogram in Fig. 4). Nonetheless, if we look at pairs with positive synergistic score and acting nontrivially on each other’s side effect, a few interesting categories emerge. For what concerns the drug pairs having a high reciprocal benefit on the side effect, many of the cases listed in Additional file 6 deal with opioid receptors, where analgesics with agonist/antagonist action on different receptors are wellknown. In particular drug pairs having a μtype opioid receptor as coherent pharmacological target and some other opioid receptor (like δtype or κtype) as incoherent offtargets are quite frequent, see Additional file 6. Many different drug pairs (legal or less) show this kind of mutualistic benefit on the majority of their offtargets. Another family of examples found in Additional file 6, concerns drugs acting as agonists or partial agonists on β _{1} adrenergic receptors, but having opposite effects on β _{2} adrenergic receptors (here considered as offtarget of several drugs, according to DrugBank pharmacological action classification). An example of a pair of drugs having b _{ i }(j)=b _{ j }(i)=1 is given by Modafinil and Lisdexamfetamine. Both act as inhibitors on the sodiumdependent dopamine transporter, but the first drug is a partial agonist for the Alpha1B adrenergic receptor, while the second is an antagonist for the same target, hence these actions tend to cancel each other. Both drugs have no other known target. As an example of a pair for which b _{ i }(j)=1 and b _{ j }(i) positive but small (Additional file 7), we can mention Icosapent and Rosiglitazone. Both have an agonist action on Peroxisome proliferatoractivated receptor γ (pharmacological target), while they have incoherent actions on the offtarget Longchainfattyacid–CoA ligase 4 (one is an inducer, while the other is an inhibitor). While Rosiglitazone has no other target, Icosapent has 8 more targets, hence the two side effect score coefficients are 1 and 1/9. A pair for which the side effect worsen is given by Amitriptyline and Paliperidone (see Additional file 9). Both have a large number of targets (33 for Amitriptyline and 17 for Paliperidone), only a fraction of which are pharmacological (1 for Amitriptyline and 5 for Paliperidone). Both act as antagonist on the common pharmacological target 5hydroxytryptamine receptor 2A, but they also share 5 offtargets, on 4 of which (mostly adrenergic receptors) the action is coherent. Only on the Histamine H1 receptor they exert an incoherent action (one as antagonist, the other as agonist), hence overall both b _{ i }(j) and b _{ j }(i) are negative.
In DrugBank a set of “BioInteractions” (i.e., drugdrug interactions) is reported for each drug, based on common targets, but also on other mechanisms like assimilation and clearance (using information on enzymes and transporters affecting the drugs). These drugdrug interactions are limited to pharmacological targets, and essentially overlap with the pharmacological drug pairs used in our analysis and shown for instance in Additional file 1: Figure S4. Having lumped together many modes of action into positive and negative signs, however, allows us to perform a further step, namely to quantify the level of synergism of these drugdrug interactions and hence search for pairs with high synergistic benefit. Furthermore, DrugBank biointeractions do not take into account offtargets as we do here, hence no investigation of side effect is possible at all if we limit ourself to the information currently provided by DrugBank.
It is worth observing that the presence of positive and negative signs on the interactions and the limited knowledge of the drugtarget doseresponse curves (for single drugs and for drug pairs) mean that classical definitions of synergism, as given in the framework of epistasis analysis or in drug combination theory [39], cannot be applied in our context. The methodology adopted in this paper, counting the overlap of two drugs and splitting it into coherent and incoherent, seems to us the simplest possible way to generalize the notion of synergism to the present context.
Clearly, as also some of the examples mentioned above show, our definitions of synergism and side effect can be oversimplified in certain cases. More sophisticated variants of the approach discussed here consist for instance in replacing the main target/offtarget distinction that we use with some other criteria, like an ontological classification of the main targets of a drug (often known), or a colocalization on a specific pathway of interest for a disease. For instance many cases of synergistic therapeutic reinforcements happen because different drugs act on different targets, located for example on redundant branches of the same pathway or on distinct but complementary pathways [29]. As a rule of thumb, when the synergism is due to indirect mechanisms, it is more difficult to capture in largescale models. A classical example is augmentin, a combination of amoxicillin and clavulanic acid [7]. The calvulanic acid inhibits one of the degradation pathways of amoxicillin hence exerting an indirect, yet strong, synergistic effect. Another related shortcoming of the analysis performed here is that drugs are considered as on/off. Clearly dosage and timing of compound application are factors that could be taken into account in a model. For instance, the idea that two drugs having modes of action with the same sign lead to an enhanced effect on a common target is more likely to happen at low dosage than at saturation, where only the strongest drug is likely to bind to the target (and competitive, instead of cooperative, effects may appear). In some cases, the information on pharmacodynamics, assimilation and clearance mechanisms available in DrugBank could be used for this scope, but it is unlikely that a systematic analysis can be built on these bases and be feasible at network level. The caveat in this case is that very little can be said on how data on single drug response remain valid for mixing of drugs, unless extensive experimental evaluation is performed. The broad range of possibilities involving drug pairs, listed for instance in [29, 39] clearly shows that only direct experimental evidence can provide an irrefutable classification of combined modes of action. Nevertheless, we believe that our signbased classification can provide a number of clues to understand drugtarget networks at functional level, and help shed light into its (networkwide) functional organization.
Conclusions
In this paper we have introduced the idea that the mode of action of a drugtarget interaction can be used to construct a signed drugtarget network. Characterizing the edges of a drugtarget network with sign adds a connotation of functional nature to the network, and enables the (networkwide) investigation of a series of properties of importance in the context of drug combinatorics.
Methods
A drugtarget network is constructed using the DrugBank dataset (version 4.3, downloaded November 2015). It is a bipartite graph having two classes of nodes: drugs and targets. The edges represent known actions of a drug on a target. Both FDA approved and experimental drugs are considered. The data for this network are given in Table 1. About 68 % of the total of DrugBank drugtarget interactions are associated to human targets, see again Table 1. Selecting only human targets, a subnetwork can be extracted. It is on this subnetwork that the paper is focused.
Pharmacological actions and offtargets. The drugtarget actions (i.e., edges of the drugtarget graph) can be subdivided into “pharmacological actions” and “unknown” or “no pharmacological action”. Normally the pharmacological actions correspond to “main targets”, for which the drug has been designed and the remaining to “offtarget actions”, i.e., side effects of the drug, see Table 1 and Additional file 11.
Sign of the drugtarget actions. The vast majority of categories used by DrugBank to describe how a drug acts on a target can be classified into 2 modes of action: “positive” or “negative”. The specific categories falling into the two modes are listed in Table 2. Some categories, like “antibody”, “cofactor”, etc. cannot be classified as positive or negative edges, and hence are not associated to any sign, see again Table 2. See also the adjacency matrix in Additional file 11.
Skewness of the coherent/incoherent drug pairs distribution. Consider the signed human drugtarget network. For the drug pair (i, j), let c _{ ij } be the number of targets in common, split into \( c_{ij}^{++},\) \( c_{ij}^{} \) and \( c_{ij}^{+}\) according to the signs of the corresponding drugtarget edges. Denote \( \xi = c_{ij}^{++} + c_{ij}^{} \) the total number of coherent edge pairs for the drug pair, and
the probability of an edge pair to be coherent over the entire network. Then the probability of drawing at most ξ coherent edge pairs out of c _{ ij } edge pairs is the cumulative binomial
If α=0.05 is the threshold for statistical significance, when the pvalue 1−P[x≤ξ]<α the drug pair (i, j) is enriched for coherent actions. A similar calculation applies also for the significance of the incoherent actions of the pair (i, j).
Quantification of pharmacological synergism in drug pairs. If p _{ ij } is the number of signed pharmacological targets shared by the drugs i and j (split, according to the signs, into \( p_{ij}^{++}\), \( p_{ij}^{}\) and \( p_{ij}^{+}\), such that \( p_{ij}^{++}+ p_{ij}^{}+ p_{ij}^{+} = p_{ij}\)), then the synergistic score coefficient is defined as
Since p _{ ij }≤p _{ i }, it must be a _{ i }(j)∈[−1, 1], and, by construction, sign(a _{ i }(j))=sign(a _{ j }(i)). Targets that are classified as pharmacological for one drug of a pair but not for the other are not counted in p _{ ij }.
Side effect quantification in drug pairs. Let s _{ i } be the side effect of the ith drug, i.e., the number of human targets that are not classified in DrugBank as pharmacological targets. If s _{ ij } is the number of common targets of the drug pair (i, j) which are not pharmacological targets for at least one of the two drugs, and \( s_{ij}^{++}\), \( s_{ij}^{}\) and \( s_{ij}^{+}\) specify theirs signs, then the side effect score coefficient is defined as
Clearly b _{ i }(j)∈[−1, 1], and, by construction, sign(b _{ i }(j))=sign(b _{ j }(i)).
Abbreviations
FDA, Federal Drug Administration
References
 1
Frantz S. Drug discovery: Playing dirty. Nature. 2005; 437(7061):942–943.
 2
Imming P, Sinning C, Meyer A. Drugs, their targets and the nature and number of drug targets. Nat Rev Drug Discov. 2007; 6(2):821–34.
 3
Overington JP, AlLazikani B, Hopkins AL. How many drug targets are there?. Nat Rev Drug Discov. 2006; 5(12):993–996.
 4
Barabási AL, Gulbahce N, Loscalzo J. Network medicine: a networkbased approach to human disease. Nat Rev Genet. 2011; 12(1):56–68.
 5
Csermely P, Korcsmáros T, Kiss HJM, London G, Nussinov R. Structure and dynamics of molecular networks: A novel paradigm of drug discovery. Pharmacol Ther. 2013; 138(3):333–408.
 6
Hopkins AL. Network pharmacology. Nat Biotechnol. 2007; 25(10):1110–1111.
 7
Keith CT, Borisy AA, Stockwell BR. Innovation: Multicomponent therapeutics for networked systems. Nat Rev Drug Discov. 2005; 4(1):71–78.
 8
Tang J, Aittokallio T. Network pharmacology strategies toward multitarget anticancer therapies: From computational models to experimental design principles. Curr Pharm Des. 2014; 20(1):20–36.
 9
Sun X, Vilar S, Tatonetti NP. Highthroughput methods for combinatorial drug discovery. Sci Transl Med. 2013; 5(205):205rv1.
 10
Chen X, Yan CC, Zhang X, Zhang X, Dai F, Yin J, Zhang Y. Drugtarget interaction prediction: databases, web servers and computational models. Brief Bioinform. 2015;17(4).
 11
Lee S, Park K, Kim D. Building a drug–target network and its applications. Expert Opin Drug Discov. 2009; 4(11):1177–1189. PMID: 23480435.
 12
Yildirim MA, Goh KI, Cusick ME, Barabási AL, Vidal M. Drugtarget network. Nat Biotechnol. 2007; 25(10):1119–1126.
 13
Yamanishi Y, Araki M, Gutteridge A, Honda W, Kanehisa M. Prediction of drug–target interaction networks from the integration of chemical and genomic spaces. Bioinformatics. 2008; 24(13):i232–i240.
 14
Erlanson DA, Hansen SK. Making drugs on proteins: sitedirected ligand discovery for fragmentbased lead assembly. Curr Opin Chem Biol. 2004; 8(4):399–406.
 15
Surade S, Blundell TL. Structural biology and drug discovery of difficult targets: The limits of ligandability. Chem Biol. 2012; 19(1):42–50.
 16
Chen D, Zhang H, Lu P, Liu X, Cao H. Synergy evaluation by a pathwaypathway interaction network: a new way to predict drug combination. Mol BioSyst. 2016; 12:614–623.
 17
Keiser MJ, Setola V, Irwin JJ, Laggner C, Abbas AI, Hufeisen SJ, Jensen NH, Kuijer MB, Matos RC, Tran TB, Whaley R, Glennon RA, Hert J, Thomas KLH, Edwards DD, Shoichet BK, Roth BL. Predicting new molecular targets for known drugs. Nature. 2009; 462(7270):175–181.
 18
Young DW, Bender A, Hoyt J, McWhinnie E, Chirn GW, Tao CY, Tallarico JA, Labow M, Jenkins JL, Mitchison TJ, Feng Y. Integrating highcontent screening and ligandtarget prediction to identify mechanism of action. Nat Chem Biol. 2007; 4(1):59–68.
 19
Cheng F, Liu C, Jiang J, Lu W, Li W, Liu G, Zhou W, Huang J, Tang Y. Prediction of drugtarget interactions and drug repositioning via networkbased inference. PLoS Comput Biol. 2012; 8(5):e1002503.
 20
Chen X, Liu MX, Yan GY. Drugtarget interaction prediction by random walk on the heterogeneous network. Mol BioSyst. 2012; 8:1970–1978.
 21
Cichonska A, Rousu J, Aittokallio T. Identification of drug candidates and repurposing opportunities through compound–target interaction networks. Expert Opin Drug Discov. 2015; 10(12):1333–1345. PMID: 26429153.
 22
Facchetti G, Zampieri M, Altafini C. Predicting and characterizing selective multiple drug treatments for metabolic diseases and cancer. BMC Syst Biol. 2012; 6:115.
 23
Nguyen TTT, Chua JKK, Seah KS, Koo SH, Yee JY, Yang EG, Lim KK, Pang SYW, Yuen A, Zhang L, Ang WH, Dymock B, Lee EJD, Chen ES. Predicting chemotherapeutic drug combinations through gene network profiling. Sci Rep. 2016; 6:18658.
 24
Wishart DS. DrugBank: a comprehensive resource for in silico drug discovery and exploration. Nucleic Acids Res. 2006; 34(90001):D668–D672.
 25
Iacono G, Altafini C. Monotonicity, frustration, and ordered response: an analysis of the energy landscape of perturbed largescale biological networks. BMC Syst Biol. 2010; 4:83.
 26
Iacono G, Ramezani F, Soranzo N, Altafini C. Determining the distance to monotonicity of a biological network: a graphtheoretical approach. IET Syst Biol. 2010; 4(3):223–235.
 27
Milo R, ShenOrr S, Itzkovitz S, Kashtan N, Chklovskii D, Alon U. Network motifs: Simple building blocks of complex networks. Science. 2002; 298(5594):824–827.
 28
Sontag ED. Syst Synth Biol. 2007; 1(2):59–87.
 29
Jia J, Zhu F, Ma X, Cao ZW, Li YX, Chen YZ. Mechanisms of drug combinations: interaction and network perspectives. Nat Rev Drug Discov. 2009; 8(2):111–128.
 30
Campillos M, Kuhn M, Gavin AC, Jensen LJ, Bork P. Drug target identification using sideeffect similarity. Science. 2008; 321(5886):263–266.
 31
Lounkine E, Keiser MJ, Whitebread S, Mikhailov D, Hamon J, Jenkins JL, Lavan P, Weber E, Doak AK, Côté S, Shoichet BK, Urban L. Largescale prediction and testing of drug activity on sideeffect targets. Nature. 2012; 486:361–67.
 32
Mizutani S, Pauwels E, Stoven V, Goto S, Yamanishi Y. Relating drug–protein interaction network with drug side effects. Bioinformatics. 2012; 28(18):i522–i528.
 33
Takarabe M, Kotera M, Nishimura Y, Goto S, Yamanishi Y. Drug target prediction using adverse event report systems: a pharmacogenomic approach. Bioinformatics. 2012; 28(18):i611–i618.
 34
Facchetti G, Iacono G, Altafini C. Computing global structural balance in largescale signed social networks. Proc Nat Ac Sci. 2011; 108:20953–8.
 35
Shanks RH, Rizzieri DA, Flowers JL, Colvin OM, Adams DJ. Preclinical evaluation of gemcitabine combination regimens for application in acute myeloid leukemia. Clin Cancer Res. 2005; 11(11):4225–4233.
 36
Price TJ, Peeters M, Kim TW, Li J, Cascinu S, Ruff P, Suresh AS, Thomas A, Tjulandin S, Zhang K, Murugappan S, Sidhu R. Panitumumab versus cetuximab in patients with chemotherapyrefractory wildtype KRAS exon 2 metastatic colorectal cancer (ASPECCT): a randomised, multicentre, openlabel, noninferiority phase 3 study. Lancet Oncol. 2014; 15(6):569–579.
 37
Douillard JY, Oliner KS, Siena S, Tabernero J, Burkes R, Barugel M, Humblet Y, Bodoky G, Cunningham D, Jassem J, Rivera F, Kocákova I, Ruff P, BłasińskaMorawiec M, Šmakal M, Canon JL, Rother M, Williams R, Rong A, Wiezorek J, Sidhu R, Patterson SD. Panitumumab–FOLFOX4 treatment and RAS mutations in colorectal cancer. N Engl J Med. 2013; 369(11):1023–1034.
 38
Liu Y, Wei Q, Yu G, Gai W, Li Y, Chen X. DCDB 2.0: a major update of the drug combination database. Database. 2014; 2014(0):bau124–bau124.
 39
Chou TC. Theoretical basis, experimental design, and computerized simulation of synergism and antagonism in drug combination studies. Pharmacol Rev. 2006; 58(3):621–681.
Acknowledgments
Not applicable.
Funding
This work was supported by a grant from the Swedish Research Council (grant n. 201504390 to CA).
Availability of supporting data
All data are publicly available.
Authors’ contributions
Conceived the study: CA. Carried out research: NBT and CA. Contributed methods: NBT and CA. Wrote the manuscript: CA. All authors have read and approved the manuscript. All authors consent to participate to this paper, and consent to publish it. Both authors read and approved the final manuscript.
Competing interests
All authors declare to have no competing financial interest and no nonfinancial competing interests.
Consent for publication
Not applicable.
Ethics approval and consent to participate
Not applicable.
Author information
Affiliations
Corresponding author
Additional files
Additional file 1
Supplementary text and figures (pdf file). (PDF 125 kb)
Additional file 2
Text file containing all drug pairs having synergistic scores a _{ i }(j)≥0.5 and a _{ j }(i)≥0.5. (TXT 1903 kb)
Additional file 3
Text file containing all drug pairs having synergistic scores a _{ i }(j)≥0.5 but 0<a _{ j }(i)<0.5. (TXT 2664 kb)
Additional file 4
Text file containing all drug pairs having synergistic scores a _{ i }(j)=a _{ j }(i)=0. (TXT 40 kb)
Additional file 5
Text file containing all drug pairs having negative synergistic scores a _{ i }(j)<0 and a _{ j }(i)<0. (TXT 1082 kb)
Additional file 6
Text file containing all drug pairs having side effect scores b _{ i }(j)≥0.5 and b _{ j }(i)≥0.5. (TXT 19 kb)
Additional file 7
Text file containing all drug pairs having side effect scores b _{ i }(j)≥0.5 but 0<b _{ j }(i)<0.5. (TXT 76 kb)
Additional file 8
Text file containing all drug pairs having side effect scores b _{ i }(j)=b _{ j }(i)=0. (TXT 44 kb)
Additional file 9
Text file containing all drug pairs having negative side effect scores b _{ i }(j)<0 and b _{ j }(i)<0. (TXT 113 kb)
Additional file 10
Drugbank_Overlap_drugpairs_with_DCDB. Text file containing all drug pairs overlapping with the experimental drug combination of the database DCDB. (TXT 8 kb)
Additional file 11
Drug_Target_Signed_Network.mat. Matlab file containing the drugtarget adjacency matrix. (TXT 28 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Torres, N.B., Altafini, C. Drug combinatorics and side effect estimation on the signed human drugtarget network. BMC Syst Biol 10, 74 (2016). https://doi.org/10.1186/s1291801603268
Received:
Accepted:
Published:
Keywords
 Null Model
 Modafinil
 Paliperidone
 Common Target
 Pharmacological Target