Drug combinations are a promising approach to achieve high efficacy at low doses and to overcome resistance. Drug combinations are especially useful when drugs cannot achieve effectiveness at tolerable doses, as occurs in cancer and tuberculosis (TB). However, discovery of effective drug combinations faces the challenge of combinatorial explosion, in which the number of possible combinations increases exponentially with the number of drugs and doses. A recent advance, called the dose model, uses a mathematical formula to overcome combinatorial explosion by reducing the problem to a feasible quadratic one: using data on drug pairs at a few doses, the dose model accurately predicts the effect of combinations of three and four drugs at all doses. The dose model has not yet been tested on higher-order combinations beyond four drugs. To address this, we measured the effect of combinations of up to ten antibiotics on E. coli growth, and of up to five tuberculosis (TB) drugs on the growth of M. tuberculosis. We find that the dose model accurately predicts the effect of these higher-order combinations, including cases of strong synergy and antagonism. This study supports the view that the interactions between drug pairs carries key information that largely determines higher-order interactions. Therefore, systematic study of pairwise drug interactions is a compelling strategy to prioritize drug regimens in high-dimensional spaces.
Drug combinations are important to increase efficacy and reduce resistance of treatment for infection and cancer. The major challenge is the vast number of experiments needed to scan the space of combination in order to find rare synergistic drugs and their optimal doses. In the past few years there has been an advance in the ability to predict the effects of drug cocktails, using a small number of experiments on drug pairs. These approaches have not been tested on combinations of more than a few drugs. Thus, it remains unclear whether there are useful combinations of 5–10 drugs that work at low doses. Here we show that a mathematical model can use data for drug pairs to predict ultra-high-order cocktails for E. coli and an important pathogen, M. tuberculosis. (1) We measured the effect of 124 new combinations of 3–10 antibiotic drugs, each at 13 doses. (2) We show that our model accurately predicts the effect of these cocktails using pair measurements alone. (3) We predict and verify new high-order combinations for TB that provide high efficacy at low doses, overcoming the well-known problem in TB that each drug alone has low efficacy.
Citation: Katzir I, Cokol M, Aldridge BB, Alon U (2019) Prediction of ultra-high-order antibiotic combinations based on pairwise interactions. PLoS Comput Biol 15(1):
Editor: Roland R. Regoes, ETH Zurich, SWITZERLAND
Received: July 28, 2018; Accepted: January 10, 2019; Published: January 30, 2019
Copyright: © 2019 Katzir et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This work was supported by an NIH Director’s New Innovator Award (1DP2LM011952-01) and an NIH grant (P50 GM107618-01A1) to B.B.A. This work was supported by European Research Council Executive Agency https://ec.europa.eu (grant number 693436), Minerva foundation http://www.minerva.mpg.de/weizmann/ (grant number 7125960101) and Israel Science foundation www.ISF.org.il (grant number 1349/15) to U.A. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Drugs in combination, compared to single drugs, are attractive therapies because of their potential for increased potency and slowed acquisition of drug resistance [1–14]. Combination therapy is important in cases where single drugs are not effective enough at tolerable doses, such as in the case of cancer and TB. In TB, the standard of care includes combinations of three or four drugs. For example, one current regimen for drug-sensitive TB uses a four-drug combination for two months, followed by a two-drug combination for at least four months . Drug-resistant TB therapy is more arduous, involving long-term, high-order combinations of antibiotics with more severe side effects .
It is therefore of interest to find additional combinations of drugs which are effective at low doses. However, it is difficult to experimentally test all possible combinations. This is because the number of experiments needed grows exponentially with the number of drugs and doses, a challenge known as the combinatorial explosion problem [17–20]. For example, if there are N = 10 drugs at D = 10 possible doses, there are DN ~1010 combinations. Hence, effective low-dose cocktails of multiple drugs might exist but cannot be found by exhaustive experimental search.
To overcome the combinatorial explosion problem, one approach is to measure a small number of combinations and use mathematical models to predict the rest. Recent studies explored the possibility of measuring only drug pairs to predict triplets and higher-order cocktails. Using pairs involves much fewer experiments than an exhaustive screen: the number of experiments scales quadratically with the number of drugs/doses rather than exponentially. For example, Wood et al  provided a formula based on single and pair measurements that gave excellent predictions for up to 4 antibiotics in E coli. This formula was less accurate in datasets with larger experimental noise . An alternative approach, called the pairs model, is more noise-resistant but less precise . These approaches can predict high-order cocktails only at the same doses at which the pairs and single drugs are measured. They therefore preclude a computational scan over all possible doses, which is important for finding effective cocktails.
A recent advance, the dose model , provides excellent predictions for all doses based on measuring singles and pairs at a small number of doses, even in the presence of sizable experimental noise. For example, for ten drugs, the dose model typically requires measuring the 45 possible pairs at about 10 dose combinations. This reduces the number of measurements needed from 1010 to about 500.
The dose model works by assuming that each drug i changes the effective concentration of the other drug j in a way that is described by an interaction parameter, aij. The model uses measurements of pairs at a few dose combinations to estimate the parameters aij and aji for each pair. It then uses these parameters to predict high-order combinations of three or more drugs (Methods).
The dose model was tested so far only on triplets and quadruplets of antibiotics and a triplets of cancer drugs . It is therefore unclear whether the dose model can predict interactions for higher-order combinations, which can potentially include more high-order interactions that are not captured at the level of pairs. Here we ask whether the dose model can predict combinations of up to ten drugs, using antibiotic effect on E. coli and the clinically relevant pathogen M. tuberculosis. We find that in both E. coli and M. tuberculosis, the dose model provides excellent predictions for high-order combinations. The dose model identifies several synergistic combinations of TB drugs.
The dose model predicts well the effect of combinations of 3–10 antibiotics on E. coli growth
To test the dose model, we chose 10 antibiotics with diverse mechanisms of action (Table 1). We define the effect of the drugs, g, as the reduction in growth of E. coli relative to growth with no drugs, evaluated by the optical density (OD 600nm) after 12h of growth in the presence of the drugs.
We first evaluated the dose-response curves of individual antibiotics at 13 doses (Fig 1, panels on the diagonal). Doses for each drug were spaced linearly such that the halfway effect D50 is approximately at the middle dose. These dose-response curves were well described by Hill curves (Methods) ,  (R2>0.95) with Hill coefficients in the range n = 0.95–3.25.
Fig 1. Effect of single and pairs of ten antibiotics on E coli growth was measured to determine the dose model parameters.
The panels on the diagonal shows the relative growth reduction of E. coli grown with individual antibiotics (Table 1) at 13 doses (dots). Doses for each drug were spaced linearly such that the halfway effect D50 is approximately at the middle dose. These dose-response curves are well described by Hill curves (solid line). The off-diagonal panels show the two-drug responses for all 45 pairs of the ten drugs, at 13 dose combinations (dots). For example, the panel located at the 4th row and 7th column shows the pair response of drug 4(CHL) and drug 7(MOX) (also shown in the inset figure). We mixed the drugs at the 13 doses used for the single-drug dose-response curves, thus diluting each drug by a factor of two. The pair response curves were used to determine the dose model parameters (aij). After fitting the parameters, the response curves of the pairs are well described by the dose model (solid lines). The pair interactions strength I = log(g12-g1g2+1), which indicates the mean deviation from the Bliss model (dashed line) is visualized by the color of each panel. The response curves presented here are the average of two (for combination 1–28) or three (for combination 29–115) biological replicates.
We next measured two-drug responses for all 45 pairs of the ten drugs, at 13 dose combinations. Here we used the method of Cokol et al , by using “diagonal” dose combinations: we mixed the drugs at the 13 doses used for the single-drug dose-response curves, thus diluting each drug by a factor of two. Thus, for each of the 45 drug pairs, we sampled the two-drug response surface at 13 dose combinations on the diagonal.
Fig 1 shows the 45 drug pair measurements, and compares them to the Bliss independence model ,  in which the effects (relative reduction in growth) of the two individual drugs is multiplied (Methods). Bliss independence was shown to be a reasonable model for high-order drug combinations, more than other baseline models such as the Loewe formula.
Deviations from Bliss indicate drug interaction, namely synergy and antagonism. Similarly to -, we calculated the interaction strength for all pairs by the logarithm of the deviation from the Bliss model, I = log2(g12-g1*g2+1) (visualized by color in Fig 1). Note that this formula is asymmetric, with higher |I| for synergy than antagonism, especially when interactions are strong. This asymmetry is small (<5%) for most (95%) of the interactions in this study.
We find that 17 of the drug pairs were antagonistic, 7 were synergistic and 21 showed no appreciable interaction within experimental error (zero interaction is within the 95% confidence interval, see S1 Fig).
We used these measurements to calibrate the dose model formula (Methods), by estimating interaction parameters for each pair, a12 and a21, which describe how drug 1 affects the dose of drug 2 and vise-versa:
The responses of the pairs are well described by the dose model (Fig 1, R2 = 0.9–0.98).
With the interaction parameters aij for each pair, we next used the dose model to predict the effect of all combinations of the ten drugs (Methods).
To test the predictions, we experimentally tested higher-order combinations. We chose 115 combinations of three to ten drugs (S1 Table) that span the range of predicted synergy and antagonism. The combinations tested included 27 triplets, 25 quadruplets, 16 quintuplets, 12 combinations of 6, 7 and 8 drugs, all ten 9-drug combinations and the full combination of the ten drugs.
For each combination, we measured the effect at 13 diagonal dose combinations (dots in Fig 2A). Thus, the dose combinations mixed the 13 doses used for the response curve, diluting each drug by a factor M for an M-drug combination.
Fig 2. The dose model predicts well the effect of combinations of 2–10 antibiotics on E. coli growth.
(A) The response curves for all 115 combinations (dots) compared to the dose model predictions (solid lines) and the bliss model predictions (dashed lines). The Interactions strength I = log(g12-g1g2+1), which indicates the mean deviation from the Bliss model (dashed line) is visualized by the color of each panel. The black and white squares on the top-right corner of each panel indicated which drugs where included in the cocktail. For example, in the first panel drugs AMP(1), MER(6) and MOX(7) were included in the cocktail. (B) The normalized growth compared to the dose model (x) and the Bliss model (o) predictions for 2–10 antibiotics. (B) R2 and RMSE values for the Dose model (x) and the Bliss model (o). The combination legend is also in S1 Table. The response curves presented here are the average of two (for combination 1–28) or three (for combination 29–115) biological replicates.
We find that the high-order combinations show sizable interactions when compared to the Bliss model, with R2 agreement in the range 0.5–0.85 (Fig 2B and 2C). At nine and ten drugs Bliss predictions were generally poor, with R2 = 0.6, indicating sizable interaction. We find, in agreement with Russ et al, that the Bliss multiplicative null model is better than the Loewe  null model for high-order combinations (S2 Fig and S3 Fig).
In contrast to the null models, the dose model shows excellent agreement for any number of drugs, with R2 in the range 0.88–0.96 and RMSE of 0.04–0.07. Remarkably, it makes excellent prediction for nine and ten drugs (R2>0.9).
We also compared the model residual distribution (the difference between model and measurements) to the experimental noise distribution (estimated form the difference between measurements in different biological repeats). We find that the dose model residual distribution is very similar to the experimental noise distribution, with the RMSE of the model residuals being slightly smaller (0.077 for the model and 0.078 for the experimental noise). In contrast the residuals of the Bliss or Lowe models are very different from the noise distribution (RMSE of 0.113 and 0.173), indicating model errors beyond experimental noise (S4 Fig).
We next asked whether there are specific combinations that significantly deviate from the dose model predictions. Such combinations can offer an example of high-order interactions that cannot be accounted for using the pair interactions in the dose model. We calculated the RMSE values of the model for each of the 115 combinations and compared it to the RMSE distribution evaluated from experimental noise (S5 Fig). We found that there is no significant difference between the two distributions in contrast to the RMSE distribution of the Bliss or Loewe models.
Additionally, we performed t-test for the mean difference between the model and the experiment for each of the 115 combinations (two or three repeats for each). None of the combinations showed a significant difference after correcting for multiple hypothesis testing errors by means of FDR  with 95% confidence. Although there is no combination for which the model deviate beyond the statistically expected deviation, one can test by repeated experiments the most deviating combinations to increase the statistical power and search for statistically significant combinations. In order to identify combination for further study we made a list of the RMSE values for each of the combinations, ordered by the difference between the model error and the experimental error (see S2 Table). One can also find the residual plots (S6 Fig) useful to identify potential combinations.
After testing the model’s predictive power, we sought to search for synergistic drug combinations. We used the model to numerically navigate the space of drug-dose combinations. We asked what fraction of combinations of 3–10 drugs from the current drug panel are predicted to show synergy or antagonism. We find that predicted synergy is rare in cocktails above four drugs, with antagonism being much more common (S7A Fig). This stems from the fact that most of the pairs in the present study are antagonistic. This result is general for the dose model and does not depend on the choice of the experimental design (see S7B Fig). Simulations of a hypothetical scenario in which we set most drug pairs to be synergistic result in mostly synergistic high-order combinations (S8 Fig).
The most synergistic cocktail relative to the Bliss model in the present drug panel is predicted to be the combination Ampicillin(1), Aztroenam(3), Meropenem(6) and Moxifloxacin(7). This combination was among those that we tested experimentally (Fig 2, combination 33). It was indeed the most synergistic cocktail among the 115 combinations tested. It provided reduction of 50% in overall dose needed to achieve a 50% reduction in growth, g = 0.5 relative to the Bliss model.it is possible that the synergy relative to Bliss is a result of the fact that three of the four drugs (AMP, AZT, MER) are beta-lactams, which probably bind the same target. Hence, a more appropriate null model for comparison in this case is Loewe (as suggested in ref ). The interaction of this combination using the Loewe model as a null model is additive.
The dose model predicts well high-order interactions in M. tuberculosis
We also tested the dose model for the effects of TB drugs on M. tuberculosis. Drug effect was defined by growth inhibition after 5 days of growth in the presence of the drugs normalized to growth with no drugs. We used nine TB drugs with various mechanisms of action (S2 Table)
The single-drug dose-response curves were well described by Hill curves ,  (R2>0.95) with Hill coefficients in the range n = 0.95–2.3. The responses of the drug pairs (experimental data originally published by Cokol et al ) were well described by the dose model (Fig 3, R2 = 0.85–0.98).
Fig 3. The effect of single and pair of 9 TB drugs on M. tuberculosis growth was measured to determine the dose model parameters.
The panels on the diagonal shows the individual antibiotics (Table 2) response at 13 doses (dots). Doses for each drug were spaced linearly such that the halfway effect D50 is approximately at the middle dose. These dose-response curves were well described by Hill curves (solid line). The off-diagonal panels show the two-drug responses for all 36 pairs of the 9 drugs, at 13 dose combinations (dots). For example, the panel located at the 4th row and 6th column shows the pair response of drug 4 and drug 6 (also shown in the inset figure). We mixed the drugs at the 13 doses used for the single-drug dose-response curves, thus diluting each drug by a factor of two. The pair response curves were used to determine the dose model parameters (aij). After fitting the parameters, the response curves of the pairs are well described by the dose model (solid line). The pairs interactions strength I = log(g12-g1g2+1), which indicates the mean deviation from the Bliss model (dashed line) is visualized by the color of each panel. The response curves presented here are the average of two biological replicates.
As above, we calculated the interaction strength I = log2(g12-g1g2+1) for all 36 pairs (visualized by color in Fig 3). We find that 16 of the pairs were antagonistic, and 6 were synergistic and 14 showed no appreciable interaction within experimental error (zero interaction is within the 95% confidence interval).
We used the interaction parameters aij to predict the effects of combinations of all combinations of drugs at all doses. Based on these predictions, we chose combinations predicted to be most synergistic and most antagonistic. We experimentally tested 9 TB drug combinations, including 4 triplets, 4 quadruplets and 1 quintuplet. Each combination was tested in the 13 diagonal dose combination, with two biological replicates.
We find that the combinations show synergy and antagonism—they differ from the Bliss model predictions (Fig 4). For example, the synergistic combination Bedaquiline(1), Clofazimine(2), Ethionamide(3) and Isoniazid(5) provided reduction of 40% in dose needed to achieve a 50% reduction in growth, g = 0.5 relative to the Bliss model (and also relative to the Loewe model that gives almost identical prediction to the Bliss in this case).
Fig 4. The dose model predicts well the effect of combinations of 3–5 TB drugs on M. tuberculosis growth.
(A) The response curves for all 9 combinations (dots) compared to the dose model predictions (solid lines) and the bliss model predictions (dashed lines). The Interactions strength I = log(g12-g1g2+1), which indicates the mean deviation from the Bliss model (dashed line) is visualized by the color of each panel. (B) The normalized growth compared to the Dose model (x) and the Bliss model (o) predictions for 3–5 TB drugs. (C) R2 and RMSE values for the Dose model (x) and the Bliss model (o). The full names of the drugs are in Table 2. The response curves presented here are the average of two biological replicates.
The dose model provides excellent prediction for all measured combinations, with R2 values exceeding 0.9 and RMSE of about 0.06 (Fig 4).
We find that the dose model of ref  provides excellent predictions for combinations of up to ten drugs in E coli and M tuberculosis, based on measurements of drug pairs. Previously this model was tested only up to combinations of four drugs.
Because the dose-model uses information from pairwise interactions, it provides great efficiency in the search for ultra-high-order drug synergies. The efficiency stems from the fact that there is a quadratic number of pairwise combinations n(n-1)/2, which is much less than the exponential number of all possible combinations among n drugs, 2^n. Therefore, conducting pairwise experiments and making predictions on high-order combinations provides an exponential increase in efficiency. For example, there are around 20 antibiotics specifically effective against M. tuberculosis. There are 190 pairwise interactions among these antibiotics, whereas there are 1 million different combinations involving any number of drugs. When one adds the number of possible doses, the efficiency of using pairs increases even more.
The use of the present approach may help to discover ultra-high-order antibiotic synergies against M. tuberculosis, with lower toxicity and higher efficacy than current treatment. Combinations may also slow down antibiotic resistance, because multiple bottlenecks are introduced by drugs with different targets. In this study, we identified one such four-antibiotic synergy against M. tuberculosis (Bdq+Clz+Eta+Inh).
The present results bear on the question of high-order interactions beyond pairs. High-order drug interactions can be considered as an amalgam of all possible interactions among the constituents –, . For example, the observed interaction of the combination A+B+C is the result of three pairwise interactions (A+B, A+C, B+C) and one emergent interaction among all three drugs (A+B+C). In this setting, the emergent interaction is a quantity that cannot be predicted using pairwise interactions. Thus, the emergent interactions depend on the null model used, such as the Bliss model, Lowe model or Isserlis-like formula used in . For example, using different null models, several previous studies found prevalent emergent interactions–, . The present finding that pairwise interactions can successfully predict high-order interactions using the dose-model is in agreement with the hypothesis that emergent interactions are rare.
The majority of drug interactions in high-order combinations studied here were antagonistic compared to Bliss and Loewe. This antagonism follows from the fact that most pairwise interactions are antagonistic in this study, as well as in previous studies . In cases where antagonism is common, prediction methods are especially important to discover rare synergies in the vast space of antagonistic ultra-high-order combinations.
Materials and methods
Strains and cultures: E. coli MG1655 strain and a pantothenate and leucine auxotrophic strain of M. tuberculosis H37Rv were used for experiments. E. coli cells were grown overnight to saturation in LB, M. tuberculosis cells were grown to mid-log phase in 7H9 medium supplemented with 0.05% Tween 80, 0.2% glycerol, 10% BBL Middlebrook OADC enrichment (VWR), 50μg/ml leucine, and 24μg/ml pantothenate.
Growth inhibition assays: Drugs were dissolved in DMSO and stored at -20°C. The top concentration used for each drug is given in Table 1. Assays were performed in clear, flat bottom 384-well microplates by dispensing nanoliter volumes of drugs using a digital drug dispenser (D300e Digital Dispenser, HP). Dispense locations were randomized within each plate in order to minimize plate position effect. E. coli and M. tuberculosis cells were diluted using fresh media to an OD of 0.01 and 0.05, respectively. 50 μL diluted cells were added to each well. Plates were sealed with aluminum plate seals and incubated without shaking at 37°C. OD600 for E. coli or M. tuberculosis plates were measured after 12 hours or 5 days, respectively (see example of growth curves in S9 Fig). Each measurement was done in two or three biological replicates.
Data analysis: Using the randomization map from the digital dispense software, plate measurement data were reconstituted and analyzed (MATLAB, Mathworks) as described above.
Single-drug dose-response curves: The single-drug dose-response curves were fitted by Hill curves (, ), characterized by a steepness parameter (Hill coefficient) ni and a halfway point D0i equal to the drug concentration of 50% effect:
Where gi,Di are the effect (growth inhibition) and the dose of drug i.
Where gi, is the single-drug effect of drug i and Di is the dose of drug i.
Dose Model: The dose model  is an extension of the Bliss model to include the effect of drugs interactions. It is based on the product of the effects of all drugs in the cocktail, but not at their true doses but rather at effective doses that differ from the true doses due to interactions with the other drugs in the cocktail. This interaction is modeled by introducing Michaelis-Menten-like interaction terms between drug pairs, where drug i changes the effective dose of drug j (and vice versa):
Where gi,gj are the single-drug dose responses (described by Hill curves); gij is the pair response of drug i and drug j; Di,Dj are the doses of drug i and drug j in the cocktail; Dieff,Djeff are the effective doses; D0i,D0j are the halfway doses of drug i and drug j; aij,aji are the interaction parameters.
Using the single and pair dose response data the interaction parameters aij,aji are fitted using MATLAB “fit” function. The effect of any drug combination is predicted by the dose model formula, by assuming the effects of the drugs on each other’s effective concentration are multiplicative:
Where the second equation is numerically solved using MATLAB “fmincon” function.
To minimize the effect of experimental noise, we used for predictions single-dose responses and combination-dose responses measured on the same plate.
Loewe additivity Model : The model assumes that the effect of drugs in a combination is additive in the dose, meaning that their combined effect is the same across all combinations that have the same total normalized dose:
Where Di,Dj are the doses of drug i and drug j which gives X% effect in the pair combination and DiX are the doses that gives X% effect when applied alone.
S1 Fig. The experimental error of the interaction index can be estimated from repeated measurements.
(A) The interaction index (I = log2(g12-g1*g2+1)) for pairs of replicate 1 versus replicate 2. The errorbars are two standard deviations (σI). (B) The distribution of the differences between the two replicates. The line is the pdf of normal disterbution with σI1−I2 = 0.066. This gives estimation of σI = 0.047 for one replicate, and for the mean value of the two replicates. The threshold for interaction was set to be to give confidence level of 95%.
S2 Fig. The Bliss multiplicative null model is better than the Loewe null model for high-order combinations against E.coli.
R2 (A) and RMSE values (B) for the Bliss (o), Loewe (*) and Dose(x) Models. As in  the Bliss model outperform the Loewe model when the number of drugs increase.
S3 Fig. The Bliss multiplicative null model is better than the Loewe null model for high-order combinations against TB.
R2 (A) and RMSE values (B) for the Bliss (o), Loewe (*) and Dose(x) Models. As in  the Bliss model outperforms the Loewe model when the number of drugs increase.
S4 Fig. The Dose model error is consistent with the experimental error.
(A) Comparison between the distribution of residuals of the dose model, to the distribution of residuals of the experimental repeats compared to each other. In the calculation of the residuals of the experimental repeats we corrected for experimental dose variations and other batch effects in the different repeats using the Bliss model. This means that we calculated the residuals of the deviation from Bliss (g-B) in the different biological repeats. This reduced the experimental RMSE from 0.13 to 0.078. Dose model RMSE is 0.077, very similar to the experimental value. This is in contrast to the RMSE of the bliss model (B) which is 0.113 and Loewe model which is 0.173. Kolmogorov–Smirnov test finds that the Dose model residual distribution is narrower than the experimental distribution(p = 0.97), and the Bliss and Loewe are wider than the experimental distribution (p<10−10)
S5 Fig. The distribution of errors for the different combinations is consistent with the experimental error.
(A) RMSE values of the dose model for each of the 115 combinations (two to three repeats) compared to the expected RMSE distribution evaluated from the repeated measurements (corrected for dose variation, see S4 Fig). We found that there is no significant difference between the two distributions (Kolmogorov–Smirnov test, p = 0.2) in contrast to the RMSE distribution of the Bliss (B) or Loewe models (C) (Kolmogorov–Smirnov test, p < 10−8).
S7 Fig. Predicted synergy is rare in cocktails above four drugs, with antagonism being much more common.
(A) Histogram of the predicted interaction parameter I (at g = 0.5) for all possible 3–8 antibiotic combinations in this study. (B) The same as (A) but with non-diagonal design. Here I is the average value of 10 random doses with g = 0.4–0.6.
S8 Fig. Predicted synergy is common in multi-drug cocktails when the pair interaction is mostly synergistic.
Histogram of the interaction parameter I(at g = 0.5) for all possible 3–8 drug combinations of simulated data, in which most of the pair interaction were set to be synergistic.
S9 Fig. We used growth curves of E. coli and Mtb to choose the end-time of the experiments.
(A) For E. coli experiments, we recorded the growth for 16 hours, with measurements every 5 minutes. (B) For Mtb experiments, we measured growth in four time points: day 3, 4, 5 and 7. 96 growth curves recorded in each of these experiments are shown with thin colored lines, and the median OD600 at each time point is shown with a thick black line. For E. coli and Mtb experiments, cultures reached saturation around 12 hours and 5 days, respectively. We chose these time points as the end-time for our drug interaction screen experiments as this end-point closely corresponds to growth rate.
S2 Table. List of RMSE values for the dose model and for the experimental repeats and their differences for all 115 three to ten antibiotic combinations tested against E. coli in this study.
The list is ordered by the maximal difference between the model and the experiment. See S1 Table for combination legend.
S1 Data. The data is OD at 600nm.
(1) allpairsof10Ecoli—the effect of all pair combinations of ten antibiotics on Ecoli growth.
(2) allpairsof10Mtb—the effect of all pair combinations of ten antibiotics on MTB growth.
(3) Ecoli030817—the effect of 20 three antibiotics combinations on Ecoli growth.
(4) Ecoli0609—the effect of up to five antibiotics combinations on Ecoli growth.
(5) Mtb0609results—the effect of up to five antibiotics combinations on MTB growth.
(6) 0822 –(a) 0822platessetup—the setup of up to 10 antibiotics combinations on Ecoli growth. (b)-(d) 0822rep1/2/3/—the results of three replicates of the experiment.
Chait R., Vetsigian K., and Kishony R., “What counters antibiotic resistance in nature?,” Nat. Chem. Biol., vol. 8, no. 1, pp. 2–5, Jan. 2012.
DeVita V. T., Young R. C., and Canellos G. P., “Combination versus single agent chemotherapy: a review of the basis for selection of drug treatment of cancer,” Cancer, vol. 35, no. 1, pp. 98–110, Jan. 1975. pmid:162854
Beppler C. et al., “Uncovering emergent interactions in three-way combinations of stressors.,” J. R. Soc. Interface, vol. 13, no. 125, p. 20160800, Dec. 2016. pmid:27974577
Tekin E., Beppler C., White C., Mao Z., Savage V. M., and Yeh P. J., “Enhanced identification of synergistic and antagonistic emergent interactions among three or more drugs.,” J. R. Soc. Interface, vol. 13, no. 119, p. 20160332, Jun. 2016. pmid:27278366
Beppler C. et al., “When more is less: Emergent suppressive interactions in three-drug combinations,” BMC Microbiol., vol. 17, no. 1, p. 107, Dec. 2017. pmid:28477626
Fitzgerald J. B., Schoeberl B., Nielsen U. B., and Sorger P. K., “Systems biology and combination therapy in the quest for clinical efficacy,” Nat. Chem. Biol., vol. 2, no. 9, pp. 458–466, Sep. 2006. pmid:16921358
Hopkins A. L., “Network pharmacology: the next paradigm in drug discovery,” Nat. Chem. Biol., vol. 4, no. 11, pp. 682–690, Nov. 2008. pmid:18936753
Lehár J. et al., “Synergistic drug combinations tend to improve therapeutically relevant selectivity,” Nat. Biotechnol., vol. 27, no. 7, pp. 659–666, Jul. 2009. pmid:19581876
Xu L. et al., “Combined EGFR/MET or EGFR/HSP90 inhibition is effective in the treatment of lung cancers codriven by mutant EGFR containing T790M and MET,” Cancer Res., vol. 72, no. 13, pp. 3302–3311, Jul. 2012. pmid:22552292
Bock C. and Lengauer T., “Managing drug resistance in cancer: lessons from HIV therapy,” Nat. Rev. Cancer, vol. 12, no. 7, pp. 494–501, Jul. 2012. pmid:22673150
Holohan C., Van Schaeybroeck S., Longley D. B., and Johnston P. G., “Cancer drug resistance: an evolving paradigm,” Nat. Rev. Cancer, vol. 13, no. 10, pp. 714–726, Oct. 2013. pmid:24060863
Hu C.-M. J. and Zhang L., “Nanoparticle-based combination therapy toward overcoming drug resistance in cancer,” Biochem. Pharmacol., vol. 83, no. 8, pp. 1104–1111, Apr. 2012. pmid:22285912
Yeh P. J., Hegreness M. J., Aiden A. P., and Kishony R., “Drug interactions and the evolution of antibiotic resistance,” Nat. Rev. Microbiol., vol. 7, no. 6, pp. 460–466, Jun. 2009. pmid:19444248
Tekin E. et al., “Prevalence and patterns of higher-order drug interactions in Escherichia coli,” npj Syst. Biol. Appl., vol. 4, no. 1, p. 31, Dec. 2018.
Horsburgh C. R., Barry C. E., and Lange C., “Treatment of Tuberculosis,” N. Engl. J. Med., vol. 373, no. 22, pp. 2149–2160, Nov. 2015. pmid:26605929
MA E., Kim S., PG S., and et al, “Standard short-course chemotherapy for drug-resistant tuberculosis: Treatment outcomes in 6 countries,” JAMA, vol. 283, no. 19, pp. 2537–2545, May 2000. pmid:10815117
Dancey J. E. and Chen H. X., “Strategies for optimizing combinations of molecularly targeted anticancer agents,” Nat. Rev. Drug Discov., vol. 5, no. 8, pp. 649–659, Aug. 2006. pmid:16883303
Lundstedt T. et al., “Experimental design and optimization,” Chemom. Intell. Lab. Syst., vol. 42, no. 1–2, pp. 3–40, Aug. 1998.
Calzolari D. et al., “Search Algorithms as a Framework for the Optimization of Drug Combinations,” PLoS Comput Biol, vol. 4, no. 12, p. e1000249, Dec. 2008. pmid:19112483
Feala J. D., Cortes J., Duxbury P. M., Piermarocchi C., McCulloch A. D., and Paternostro G., “Systems approaches and algorithms for discovery of combinatorial therapies,” Wiley Interdiscip. Rev. Syst. Biol. Med., vol. 2, no. 2, pp. 181–193, 2010. pmid:20836021
Wood K., Nishida S., Sontag E. D., and Cluzel P., “Mechanism-independent method for predicting response to multidrug combinations in bacteria,” Proc. Natl. Acad. Sci., vol. 109, no. 30, pp. 12254–12259, Jul. 2012. pmid:22773816
Zimmer A., Katzir I., Dekel E., Mayo A. E., and Alon U., “Prediction of multidimensional drug dose responses based on measurements of drug pairs.,” Proc. Natl. Acad. Sci. U. S. A., vol. 113, no. 37, pp. 10442–7, Sep. 2016. pmid:27562164
Zimmer A., Tendler A., Katzir I., Mayo A., and Alon U., “Prediction of drug cocktail effects when the number of measurements is limited,” PLOS Biol., vol. 15, no. 10, p. e2002518, Oct. 2017. pmid:29073201
Chou T. C. and Talalay P., “Quantitative analysis of dose-effect relationships: the combined effects of multiple drugs or enzyme inhibitors,” Adv. Enzyme Regul., vol. 22, pp. 27–55, 1984. pmid:6382953
Chou T.-C., “Theoretical Basis, Experimental Design, and Computerized Simulation of Synergism and Antagonism in Drug Combination Studies,” Pharmacol. Rev., vol. 58, no. 3, pp. 621–681, Sep. 2006. pmid:16968952
Cokol M., Kuru N., Bicak E., Larkins-Ford J., and Aldridge B. B., “Efficient measurement and factorization of high-order drug interactions in Mycobacterium tuberculosis,” Sci. Adv., vol. 3, no. 10, p. e1701881, Oct. 2017. pmid:29026882
Bliss C. I., “The Toxicity of Poisons Applied Jointly1,” Ann. Appl. Biol., vol. 26, no. 3, pp. 585–615, 1939.
BLISS C. I., “The calculation of microbial assays,” Bacteriol. Rev., vol. 20, no. 4, pp. 243–258, Dec. 1956. pmid:13403845
Russ D. and Kishony R., “Additivity of inhibitory effects in multidrug combinations,” Nat. Microbiol., p. 1, Oct. 2018.
Bollenbach T., Quan S., Chait R., and Kishony R., “Nonoptimal microbial response to antibiotics underlies suppressive drug interactions,” Cell, vol. 139, no. 4, pp. 707–718, Nov. 2009. pmid:19914165
Wood K. B., Wood K. C., Nishida S., and Cluzel P., “Uncovering Scaling Laws to Infer Multidrug Response of Resistant Microbes and Cancer Cells,” Cell Rep., vol. 6, no. 6, pp. 1073–1084, Mar. 2014. pmid:24613352
LOEWE S., “The problem of synergism and antagonism of combined drugs,” Arzneimittelforschung., vol. 3, no. 6, pp. 285–290, Jun. 1953. pmid:13081480
Benjamini Y. and Hochberg Y., “Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing,” J. R. Stat. Soc. Ser. B, vol. 57, no. 1, 1995.
Baeder D. Y., Yu G., Hozé N., Rolff J., and Regoes R. R., “Antimicrobial combinations: Bliss independence and Loewe additivity derived from mechanistic multi-hit models.,” Philos. Trans. R. Soc. Lond. B. Biol. Sci., vol. 371, no. 1695, p. 20150294, May 2016. pmid:27160596
Sampson S. L. et al., “Protection elicited by a double leucine and pantothenate auxotroph of Mycobacterium tuberculosis in guinea pigs.,” Infect. Immun., vol. 72, no. 5, pp. 3031–7, May 2004. pmid:15102816