Systematic discovery of drug interaction mechanisms
Drug combinations are increasingly important in disease treatments, for combating drug resistance, and for elucidating fundamental relationships in cell physiology. When drugs are combined, their individual effects on cells may be amplified or weakened. Such drug interactions are crucial for treatment efficacy, but their underlying mechanisms remain largely unknown. To uncover the causes of drug interactions, we developed a systematic approach based on precise quantification of the individual and joint effects of antibiotics on growth of genome-wide Escherichia coli gene deletion strains. We found that drug interactions between antibiotics representing the main modes of action are highly robust to genetic perturbation. This robustness is encapsulated in a general principle of bacterial growth, which enables the quantitative prediction of mutant growth rates under drug combinations. Rare violations of this principle exposed recurring cellular functions controlling drug interactions. In particular, we found that polysaccharide and ATP synthesis control multiple drug interactions with previously unexplained mechanisms, and small molecule adjuvants targeting these functions synthetically reshape drug interactions in predictable ways. These results provide a new conceptual framework for the design of multidrug combinations and suggest that there are universal mechanisms at the heart of most drug interactions.
A general principle of bacterial growth enables the prediction of mutant growth rates under drug combinations. Rare violations of this principle expose cellular functions that control drug interactions and can be targeted by small molecules to alter drug interactions in predictable ways.
Drugs play a crucial role in elucidating fundamental relationships in cell physiology (Falconer et al, 2011). When drugs are combined, interactions like synergism and antagonism can occur (Loewe, 1928; Keith et al, 2005) (Fig 1A); such drug interactions are often critical for the success of multidrug treatments (Pillai et al, 2005) and can slow or accelerate antibiotic resistance evolution (Chait et al, 2007; Hegreness et al, 2008). A case in point is the synergistic combination of trimethoprim and sulfa drugs which has been applied successfully for decades (Pillai et al, 2005) even though the mechanism of synergism has long remained elusive (Nichols et al, 2011). Synergism and antagonism occur frequently between antimicrobials and are largely determined by the primary cellular target of the drugs that are combined (Yeh et al, 2006; Ocampo et al, 2014). However, synergistic drug interactions are rarely explained by the genetic interactions between the corresponding drug target genes (Cokol et al, 2011). To design combinations exploiting the full potential of existing drugs, a deeper understanding of the underlying mechanisms of drug interactions is urgently needed. Due to the vast number of possible drug combinations, general principles that are valid across diverse drug pairs could greatly facilitate the identification of drug interaction mechanisms. Recent work in this direction revealed scaling laws describing the effects of resistance mutations on drug interactions (Chait et al, 2007; Wood et al, 2014); further, the effects of three or more drugs appear largely predictable from the two-drug effects (Wood et al, 2012). However, the underlying mechanisms of most two-drug interactions remain unknown.
Drug interactions could be caused by physicochemical effects, for example when one drug simply enhances the permeability of the cell envelope for another (Jawetz & Gunnison, 1953); alternatively, they may have more complex causes, specifically if one drug triggers a regulatory response, which affects the action of another. While many genes affect the cell's sensitivity to individual drugs (Hillenmeyer et al, 2008; Nichols et al, 2011) and recent work suggested that certain drug resistance mutations may affect drug interactions (Munck et al, 2014; Wood et al, 2014; Rodriguez de Evgrafov et al, 2015), it is unclear to what extent genetic perturbations can alter drug interactions. Likewise, the cellular functions that control these interactions are largely unknown. To pinpoint their underlying causes, we developed a systematic approach for identifying genes that reshape drug interactions: using precise growth rate measurements of a genome-wide set of E. coli gene deletion mutants (Baba et al, 2006), we show that drug interactions are robust to most, but not all, genetic changes. We present a general principle encapsulating this robustness, which enables the quantitative prediction of mutant growth under drug combinations. Rare mutants violating this principle expose cellular functions governing each drug interaction. We establish that diverse drug interactions are recurrently controlled by central cellular functions, in particular polysaccharide synthesis and ATP synthesis.
Results and Discussion
We began by quantifying the growth rates of ~4,000 nonessential E. coli gene deletion mutants (Baba et al, 2006) under six different antibiotic pairs and their constituent individual drugs (Materials and Methods). We selected antibiotics with diverse modes of action (Table 1) and drug pairs covering all interaction types (Fig 1A and B). Drug concentrations were adjusted to inhibit wild-type growth under individual drugs by ~30%; the same concentrations were used when drugs were combined, leading to different levels of growth inhibition (Materials and Methods). In total, we measured over 50,000 growth curves (optical density increase over time) using a dedicated robotic system (Fig 1C–F; Materials and Methods). These growth rate measurements were highly reproducible (Supplementary Fig S1), and consistent with established changes in mutant sensitivity to individual antibiotics (Tamae et al, 2008; Liu et al, 2010; Nichols et al, 2011): for example, DNA repair mutants were sensitive to ciprofloxacin. Drugs with related mode of action had similar effects on growth of genome-wide mutants (e.g., Pearson correlation ρ = 0.68 for chloramphenicol and tetracycline), and the effects of drug combinations were usually most similar to those of the constituent drugs (Supplementary Fig S2; Materials and Methods). These observations confirmed our expectation that perturbations of cell physiology caused by drug combinations are mostly an overlay of effects caused by the constituent drugs.
|Abbreviation||Drug||Mode of action (known target)||Concentration|
|CHL||Chloramphenicol||Protein synthesis (50S ribosome subunit)||1 µg/ml|
|CPR||Ciprofloxacin||DNA replication (gyrase)||4 ng/ml|
|MEC||Mecillinam||Cell wall (Penicillin Binding Protein)||38 ng/ml|
|NIT||Nitrofurantoin||Multiple mechanisms||2 µg/ml|
|TET||Tetracycline||Protein synthesis (30S ribosome subunit)||150 ng/ml|
|TMP||Trimethoprim||Folic acid synthesis (DHFR)||80 ng/ml|
We next identified a general principle characterizing drug interactions in mutants. We hypothesized that the shape of the two-drug growth response surface g(a,b), which defines the interaction (Fig 1A), does not change qualitatively in most mutants (here, g denotes growth rate and a, b the drug concentrations). To test this hypothesis, we measured 108 mutant response surfaces in two-dimensional drug concentration matrices covering different drug pairs; this selection included mutants with strongly altered sensitivity to the individual drugs (Materials and Methods). We found that the vast majority of mutant response surfaces gmut(a,b) were well approximated by a linearly rescaled wild-type surface: gmut(a,b) =γgWT(αa,βb) with scaling factors for maximum growth rate γ and for drug concentrations α, β (Fig 2A and B; Supplementary Figs S3 and S4). The sensitivity to one or both of the drugs often changed considerably, yet the response surface shape was generally preserved (Fig 2A); these observations held for all drug pairs and for mutants affecting diverse cellular functions (Fig 2B; Supplementary Fig S4), suggesting that most genetic perturbations do not affect drug interactions.
To test whether this conservation of drug interactions holds generally, we devised a strategy for predicting genome-wide mutant growth responses to drug combinations. For all ~4,000 mutants, we calculated the expected response to the drug combination by first rescaling the wild-type response surface according to the individual drug responses measured at fixed concentrations. At the rescaled drug concentrations, we then used the interaction coefficient of the wild-type, which quantifies the response to the drug combination relative to the Bliss additive expectation (Yeh et al, 2006), to calculate each mutant's expected growth rate under the drug combination (Fig 2C; Materials and Methods). The central assumption of this procedure is that the drug interaction is universally invariant, that is it is the same in mutants as in the wild-type upon rescaling of the drug concentrations. For all drug pairs and the vast majority of mutants, the growth rates measured at fixed concentration of the drug combination (Fig 1C–F) faithfully followed this prediction (Fig 2D; Supplementary Fig S5). These observations thus revealed a general principle of bacterial growth under drug combinations, which encapsulates the high robustness of drug interactions to genetic perturbations and enables the quantitative prediction of mutant growth rates under drug combinations.
The identification of this general principle empowered us to pinpoint ‘outlier’ mutants with unexpected growth response to drug combinations. These outliers are of key interest as they could have altered drug interaction; together with functional information on the mutated gene, they can thus point at the underlying drug interaction mechanism. Clear outliers for which the observed growth rate under the drug combination (Fig 1F) deviated significantly from the expected growth rate were rare (typically < 1% of mutants; Fig 2D; Supplementary Fig S5), facilitating this investigation. We measured the response surface of the strongest outliers for each drug combination in fine resolution 12 × 8 concentration matrices (Materials and Methods). For each drug pair, we thus identified several mutants with clearly reshaped drug interaction (Fig 3, Supplementary Fig S6). Drug interactions were often weakened or removed in these mutants; they were also amplified in certain mutants and, in some cases, entirely new ‘synthetic’ drug interactions appeared (Fig 3C). In a thiamin synthesis hypomorph, chloramphenicol–trimethoprim even became reciprocally suppressive; that is, addition of chloramphenicol on top of trimethoprim increased growth and vice versa (Supplementary Fig S7). We further observed clear biases in interaction changes: chloramphenicol–nitrofurantoin suppression was weakened or entirely removed in most mutants affecting this interaction (Fig 3A); in contrast, chloramphenicol–trimethoprim antagonism was often amplified to suppression (Fig 3B and C; Supplementary Fig S6L–Q), while other drug combinations showed more balanced interaction changes in both directions (Supplementary Fig S5). These data show that different drug combinations have varying potential for changing their drug interaction type, even if the wild-type interaction is similar. Interestingly, for the additive chloramphenicol–tetracycline combination, none of the outlier mutants showed any change in drug interaction: while the sensitivity to the constituent drugs often changed, additivity was generally preserved (Supplementary Fig S3). This observation is consistent with previous results that this interaction changed little in strains that had evolved resistance to chloramphenicol–tetracycline (Munck et al, 2014). Thus, the chloramphenicol–tetracycline drug interaction appears robust to genetic perturbations. Together, these data show that most drug interactions can be removed, amplified, and even qualitatively changed to a different interaction type by rare genetic perturbations, while other interactions are robust.
We next identified cellular functions controlling drug interactions; among these, the synthesis of secreted polysaccharides (capsular and lipopolysaccharides, LPS) and ATP synthesis stood out in that they affected virtually all drug interactions. Firstly, manual inspection of the functional annotation of the outlier mutants (Keseler et al, 2005) for each drug pair, and systematic gene ontology enrichment analysis congruently exposed specific functions controlling the respective interaction (Fig 3D, Supplementary Table S1, Materials and Methods): for example, perturbing tRNA processing consistently removed chloramphenicol–nitrofurantoin suppression (Fig 3A and C; Supplementary Fig S6X and Y); similarly, ribosome production and assembly altered ciprofloxacin–tetracycline suppression (dksA, rsgA, ksgA in Supplementary Fig S5C; Supplementary Table S1), confirming previous results (Bollenbach et al, 2009). Secondly, we noticed that certain functions recurrently control multiple drug interactions. Besides polysaccharide and ATP synthesis (discussed below), chaperone deletions (ppiD, dnaK) consistently caused amplified or synthetic suppression for distinct drug pairs and perturbing amino acid synthesis (metL) amplified chloramphenicol–nitrofurantoin suppression but removed trimethoprim–chloramphenicol antagonism (Fig 3C, Supplementary Fig S6). Thus, multiple drug interactions are controlled by few recurring cellular functions.
Polysaccharide synthesis affects the majority of drug interactions (Fig 3C and D; Supplementary Fig S6): perturbing this function removed chloramphenicol–nitrofurantoin suppression and trimethoprim–mecillinam synergy (gmhB, rcsD in Fig 3C; lpcA in Supplementary Fig S6); in contrast, it led to synthetic suppression between trimethoprim and chloramphenicol (lpxM in Fig 3C; Supplementary Fig S6M and N). Polysaccharide synthesis mutants have modified outer membrane composition, which affects the uptake of molecules dependent on their chemical properties (Nikaido & Vaara, 1985). Hence, a plausible cause of these drug interactions is that bacteria regulate polysaccharide synthesis in response to certain antibiotics, which then affects the uptake of other drugs. Consistent with this mechanism, antibiotics are known to affect polysaccharide synthesis (Rothfield & Pearlman-Kothencz, 1969) and LPS synthesis mutants have increased sensitivity to chloramphenicol (Fig 3C; Supplementary Fig S6V); thus, stimulation of LPS synthesis by nitrofurantoin could explain chloramphenicol–nitrofurantoin suppression. To directly test this scenario, we removed outer membrane LPS using ethylenediaminetetraacetic acid (EDTA) (Nikaido & Vaara, 1985). Strikingly, LPS removal increased sensitivity to chloramphenicol, abolished chloramphenicol–nitrofurantoin suppression, and rendered this drug interaction purely additive (Fig 4A and B). Together, these data support that regulated changes in cell envelope composition, which affect drug uptake, are a recurring mechanism underlying chloramphenicol–nitrofurantoin suppression and other antibiotic interactions.
ATP synthesis also controls multiple drug interactions (Fig 3C and D). Specifically, the ATP synthase mutant atpF was more sensitive to trimethoprim and this sensitivity was reduced by chloramphenicol or mecillinam, leading to suppression (Fig 3C); a similar effect occurred for ciprofloxacin–tetracycline (Supplementary Fig S6S). A thiamin synthesis hypomorph, which also has perturbed energy metabolism, behaved similarly (Supplementary Fig S7). Further, ATP synthase expression increased two-fold in response to trimethoprim (Supplementary Fig S8), suggesting that cells respond homeostatically to ATP deficiency. To test independently whether ATP synthesis affects drug interactions, we specifically blocked the proton pore of the ATP synthase FO subunit using drugs (oligomycin, venturicidin; Materials and Methods). Indeed, inhibiting ATP synthase led to suppression between trimethoprim and chloramphenicol (Fig 4C and D; Supplementary Fig S9), supporting that imbalances in energy metabolism cause this synthetic drug interaction. Together, these data suggest a mechanistic scenario in which impaired ATP synthase function leads to decreased intracellular ATP levels, which may become growth-limiting in the presence of trimethoprim or ciprofloxacin; concurrent translation inhibition by chloramphenicol would reduce global ATP turnover, replenish the intracellular ATP pool (Schneider et al, 2002), and thus lead to increased growth, which would explain the observed suppressive interaction. No suppression would occur in the wild-type where ATP is likely in excess and not growth-limiting. At a molecular level, decreased ATP concentrations might limit growth by aggravating DNA repair and synthesis (Waldstein et al, 1974), which is likely the growth-limiting process in the presence of drugs targeting DNA synthesis; however, other ATP-dependent processes could also contribute to suppression. Overall, our results show that perturbations of central cellular functions, unrelated to the common antibiotic targets, can reshape diverse drug interactions.
We established a general principle of bacterial growth, which enables the prediction of mutant growth rates under drug combinations from their growth under the individual drugs alone (Fig 2, Supplementary Fig S5). This principle may hold more generally and should be tested for combinations of other challenges such as osmotic, temperature, or pH stress in future work. While conceptually similar empirical laws are an integral part of physics, they are still scarce in biology (Scott & Hwa, 2011). Even without understanding their molecular origins, such principles are powerful since they enable the prediction of quantitative phenotypes. Here, such a principle was crucial for systematically revealing antibiotic interaction mechanisms.
The identification of cellular functions controlling drug interactions offers new strategies for the rational design of multidrug combinations. Specifically, we identified targets for potential adjuvants, which could reshape antibiotic interactions: thiamin synthesis inhibitors could render the chloramphenicol–trimethoprim combination reciprocally suppressive (Supplementary Fig S7); such reciprocal suppression may slow down resistance evolution but is extremely scarce among natural antibiotic interactions (Chait et al, 2007). LPS synthesis inhibitors could remove chloramphenicol–nitrofurantoin suppression (Fig 4A and B), thus preserving advantages of an untapped drug combination while increasing its potency. Drugs inhibiting cellular functions that control antibiotic interactions (LPS synthesis, ATP synthesis, and chaperoning; Fig 3C and D) are in development (Moreau et al, 2008; Evans et al, 2010; Du et al, 2011; Balemans et al, 2012). These inhibitors could reshape drug interactions even if they have poor antimicrobial activity alone since most mutants we identified have only mild growth defects. Finally, our approach revealed that certain drug combinations are robust to mutations (Supplementary Fig S3) or change primarily toward weakened antagonism (Fig 3A). The origins of such robustness and biases in drug interaction changes are unknown. Still, such insights can be used to avoid loss of synergism due to mutations occurring in treatments, which is a serious concern (Pena-Miller et al, 2013; Munck et al, 2014). It will be exciting to extend the systematic approach presented here to drug interactions in other systems including the most worrisome pathogenic microbes and cancer.
Materials and Methods
Strains, media, and drugs
Deletion strains are from the Keio collection of 3,985 nonessential gene deletions (Baba et al, 2006). Since the strains in this collection have a kanamycin resistance marker, we introduced kanamycin resistance on a low-copy-number plasmid (pUA66; Zaslaver et al, 2006) into the parent strain (BW25113, ‘WT’). All gene deletion mutants with clear effects on drug interactions (pgpA, gmhB, metL, ppiD, lpxM, atpF, rcsD, dnaK, atpC, rfaG, rfaP, rfaC, lpcA, rep, spr) were verified by sequencing; the correct gene deletion was confirmed in all cases. All experiments were performed in lysogeny broth (LB) medium. Drugs were obtained from Sigma-Aldrich (catalogue numbers: ciprofloxacin, 17850; chloramphenicol, C0378; mecillinam, 33447; nitrofurantoin, N7878; tetracycline, 268054; trimethoprim, 92131). Drug stocks were prepared in water (ciprofloxacin, mecillinam), ethanol (chloramphenicol, tetracycline, trimethoprim) or dimethylformamide (nitrofurantoin), passed through a 0.22-μm filter, and stored in the dark at −20°C. Drugs were used at fixed concentrations that inhibit wild-type growth by ~30% (Table 1); growth remained exponential at the concentrations used (Fig 1D and E). In drug combination experiments, the same concentration as in the single drug experiments was used except for mecillinam–trimethoprim for which no growth occurred at these concentrations due to its synergistic interaction; therefore, the concentration of both drugs was reduced in the combination experiment (20 ng/ml for mecillinam and 50 ng/ml for trimethoprim). The resulting growth inhibition of the wild-type was ~60% for chloramphenicol–tetracycline, ~50% for chloramphenicol–trimethoprim, ~20% for ciprofloxacin–tetracycline, ~55% for mecillinam–trimethoprim, ~35% for trimethoprim–tetracycline, and ~30% for nitrofurantoin–chloramphenicol. Thiamin pyrophosphate (Sigma-Aldrich catalogue number C8754) was dissolved in water and stored in the dark at −20°C. Oligomycin A and venturicidin A (Szabo Scandic catalogue numbers SACSC-201551A and SACSC-202380A) were dissolved in ethanol and stored in the dark at −20°C.
Growth rate measurements
Each strain was incubated for ~20 h in one well of a 96-well plate (nontreated transparent flat bottom, Nunc) containing 200 μl medium. Cultures were inoculated using a replicator (V&P Scientific) transferring ~0.2 μl from a (thawed) overnight culture kept at −80°C with 15% glycerol. Optical density (OD) at 600 nm was measured every ~30 min in a plate reader (Tecan Infinite F500, 5 flashes, 10 ms settle time; filter: D600/20×; Chroma). The plates were incubated in an automated incubator (Liconic Storex) kept at 30°C, > 95% humidity, and shaken at 720 rpm. In addition, directly before each measurement, plates were shaken on a magnetic shaker (Teleshake; Thermo Scientific) at 900 rpm for 20 s. A customized liquid handling robot (Tecan Freedom Evo 150) was used to automate these experiments and measure over 2,000 growth curves per day. To achieve nearly identical growth conditions for all strains in each condition, the growth curves of all ~4,000 deletion strains were measured over two consecutive days using the same freshly prepared drug solution. The growth rate in exponential phase was quantified from the OD increase over time by a linear fit of log(OD) in the range 0.022 < OD < 0.22 (magenta lines in Fig 1). Late growth occurring after 1,000 min was discarded because in rare cases, fast growing strains (likely resistant mutants) overtook the population. For mecillinam, only early growth (happening before 450 min) was considered because many instances of late fast growth occurred for this drug; this effect may be due to drug decay as mecillinam is relatively unstable (Baltzer, 1979). All growth rates were normalized to the growth rate of the parent strain in the absence of drug measured on the same day. These automated measurements led to highly reproducible growth rates: replicate measurements of the entire deletion collection under chloramphenicol on different days showed a Pearson correlation of 0.94 (Supplementary Fig S1A); replicates of growth rate measurements had a variation coefficient (standard deviation over mean) of typically < 5%. Media evaporation from plates and edge effects were virtually undetectable. Mutant sensitivities to antibiotics determined from these data were consistent with published data (Supplementary Fig S10).
Two-drug response surfaces
For each drug pair, response surfaces were measured for wild-type and 18 gene deletion mutants, covering outlier mutants with strong differences between observed and expected growth rate (Fig 3A and B; Supplementary Fig S5) and additional mutants that showed clear changes in sensitivity to at least one of the constituent drugs. Response surfaces were measured using 12 × 8 or 24 × 24 two-dimensional drug concentration matrices set up with a liquid handling robot across one (12 × 8) or six 96-well plates (24 × 24), respectively. The concentration profile for each drug was set up according to where cmax was the highest concentration used, x was linearly spaced from 0 to 1 with 8, 12, or 24 steps depending on the experiment, and a = 1/3. This concentration profile was chosen to adequately sample the relevant part of the two-drug response surface where growth rate changes significantly. The points in two-dimensional concentration space where growth was measured are shown by small gray dots in all Figures. For the representation of two-dimensional response surfaces, we used the optical density 12 h after inoculation instead of the growth rate because this quantity was slightly more reproducible and yielded smoother response surfaces; this representation does not affect any of the conclusions on drug interaction changes (Supplementary Fig S1B–D). Smooth surfaces and isoboles (contour lines) were calculated by linear interpolation (Matlab function interp2) of the experimental data. The IC50 line is the isobole of 50% growth inhibition, and in practice, we used the isobole of 90% growth inhibition as the MIC (Minimal Inhibitory Concentration) line. We measured the response surface of the wild-type and of all mutants that showed a clear change in drug interaction (Fig 3C; Supplementary Fig S7) at least in duplicate; replicate response surfaces measured on different days were generally highly reproducible (Supplementary Fig S1B–D) and, in particular, all drug interaction changes in mutants were confirmed.
Expected growth rate in drug combinations
This procedure yielded accurate predictions of growth rate under drug combinations (Fig 2D; Supplementary Fig S5). For mecillinam–trimethoprim, we had to slightly adjust this procedure since the concentrations used in the drug combination had to be reduced (see ‘Strains, media, and drugs’ above). We took this into account by multiplying and with a constant factor capturing the reduced drug concentrations. Mutants that deviated from this expectation were used to identify altered drug interactions (Fig 3, Supplementary Fig S5); some mutants deviated from this expectation for other reasons, in most cases because they had extremely long lag phase. To estimate the error of the expected growth rate, we added a 5% normal distributed relative error (empirically determined, see ‘Growth rate measurements’ above), and an estimated absolute error of 0.01 (capturing the limited reproducibility of extremely low growth rates) to each growth rate measurement. The standard deviation at each expected growth rate was then numerically calculated from 10,000 randomly sampled growth rates g(a) and g(b) (from a uniform distribution between 0 and 1) for each drug pair. The dashed lines in Figs 2D and 3A and B, and Supplementary Fig S5 show two standard deviations (corresponding to a 95% confidence interval) for both the expected (x-axis) and the measured growth rate in the combination (y-axis). Due to differences in response surface shape, the resultant error estimates depend strongly on the drugs used. The density scatterplots in these and other Figures were generated using the scatplot function available at http://www.mathworks.com/matlabcentral/fileexchange/8577-scatplot.
Gene ontology enrichment analysis
To identify mutants whose growth rate in the drug combination deviated strongly from the expectation (outliers in Fig 3A and B; Supplementary Fig S5), we performed an orthogonal regression using principal component analysis (Matlab function princomp) and used the orthogonal distance to the regression line to quantify deviations from the expectation. We performed gene ontology enrichment analysis (Fig 3D; Supplementary Table S1) on the 30 outliers with the strongest deviation from the expectation. We excluded outliers for which both the expected and the observed growth rates were below 0.15 because these extremely low growth rates are hard to quantify reliably. The gene ontology database used in our analysis was retrieved from geneontology.org (released 07/15/2014) and the gene association file linking gene names to GO numbers from ecocyc.org (GOC validation date: 06/26/2014) (Keseler et al, 2005). The P-values were obtained using a custom implementation of Sherlock and Weng's GO:Termfinder software (Tavazoie et al, 1999) and Bonferroni corrected for the number of GO terms tested.
Gene expression measurements
Transcriptional regulation of the atpI promoter was measured as a proxy for ATP synthase expression in concentration gradients of different antibiotics using a promoter–GFP reporter strain (Zaslaver et al, 2006) and quantified as described (Bollenbach et al, 2009).
We thank Rosalind Allen, Eric Brown, Remy Chait, Vincent Danos, Meriem El Karoui, Roy Kishony, Alex de Luna, Adam Palmer, Teuta Pilizota, Georg Rieckh, Nassos Typas, and all members of the Bollenbach laboratory for fruitful discussions, comments on the manuscript, and technical help. This work was supported in part by Austrian Science Fund (FWF) grant No. P 27201-B22, HFSP program grant No. RGP0042/2013, and Marie Curie Career Integration grant No. 303507.
Conflict of interest
The authors declare that they have no conflict of interest.
- Supplementary Figure S1 (application/PDF, 113.6 KB)
- Supplementary Figure S2 (application/PDF, 147.4 KB)
- Supplementary Figure S3 (application/PDF, 192.4 KB)
- Supplementary Figure S4 (application/PDF, 226.7 KB)
- Supplementary Figure S5 (application/PDF, 489.4 KB)
- Supplementary Figure S6 (application/PDF, 451 KB)
- Supplementary Figure S7 (application/PDF, 193.1 KB)
- Supplementary Figure S8 (application/PDF, 73.6 KB)
- Supplementary Figure S9 (application/PDF, 95.8 KB)
- Supplementary Figure S10 (application/PDF, 124.1 KB)
- Supplementary Table S1 (MS Excel, 29 KB)
- Review Process File (application/PDF, 264.7 KB)
- Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, Baba M, Datsenko KA, Tomita M, Wanner BL, Mori H (2006) Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol 2: 2006.0008PubMedWeb of Science®Google Scholar
- Balemans W, Vranckx L, Lounis N, Pop O, Guillemont J, Vergauwen K, Mol S, Gilissen R, Motte M, Lançois D, De Bolle M, Bonroy K, Lill H, Andries K, Bald D, Koul A (2012) Novel antibiotics targeting respiratory ATP synthesis in Gram-positive pathogenic bacteria. Antimicrob Agents Chemother 56: 4131–4139CASPubMedWeb of Science®Google Scholar
- Baltzer B (1979) Degradation of mecillinam in aqueous solution. J Pharm Sci 68: 1207–1215CASPubMedWeb of Science®Google Scholar
- Bollenbach T, Quan S, Chait R, Kishony R (2009) Nonoptimal microbial response to antibiotics underlies suppressive drug interactions. Cell 139: 707–718CASPubMedWeb of Science®Google Scholar
- Chait R, Craney A, Kishony R (2007) Antibiotic interactions that select against resistance. Nature 446: 668–671CASPubMedWeb of Science®Google Scholar
- Cokol M, Chua HN, Tasan M, Mutlu B, Weinstein ZB, Suzuki Y, Nergiz ME, Costanzo M, Baryshnikova A, Giaever G, Nislow C, Myers CL, Andrews BJ, Boone C, Roth FP (2011) Systematic exploration of synergistic drug pairs. Mol Syst Biol 7: 544CASPubMedWeb of Science®Google Scholar
- Du Q, Wang H, Xie J (2011) Thiamin (vitamin B1) biosynthesis and regulation: a rich source of antimicrobial drug targets? Int J Biol Sci 7: 41–52CASPubMedWeb of Science®Google Scholar
- Evans CG, Chang L, Gestwicki JE (2010) Heat shock protein 70 (hsp70) as an emerging drug target. J Med Chem 53: 4585–4602CASPubMedWeb of Science®Google Scholar
- Falconer SB, Czarny TL, Brown ED (2011) Antibiotics as probes of biological complexity. Nat Chem Biol 7: 415–423CASPubMedWeb of Science®Google Scholar
- Hegreness M, Shoresh N, Damian D, Hartl D, Kishony R (2008) Accelerated evolution of resistance in multidrug environments. Proc Natl Acad Sci USA 105: 13977–13981CASPubMedWeb of Science®Google Scholar
- Hillenmeyer ME, Fung E, Wildenhain J, Pierce SE, Hoon S, Lee W, Proctor M, St Onge RP, Tyers M, Koller D, Altman RB, Davis RW, Nislow C, Giaever G (2008) The chemical genomic portrait of yeast: uncovering a phenotype for all genes. Science 320: 362–365CASPubMedWeb of Science®Google Scholar
- Jawetz E, Gunnison JB (1953) Antibiotic synergism and antagonism; an assessment of the problem. Pharmacol Rev 5: 175–192CASPubMedWeb of Science®Google Scholar
- Keith CT, Borisy AA, Stockwell BR (2005) Multicomponent therapeutics for networked systems. Nat Rev Drug Discov 4: 1–8CASWeb of Science®Google Scholar
- Keseler IM, Collado-Vides J, Gama-Castro S, Ingraham J, Paley S, Paulsen IT, Peralta-Gil M, Karp PD (2005) EcoCyc: a comprehensive database resource for Escherichia coli. Nucleic Acids Res 33: D334–D337CASPubMedWeb of Science®Google Scholar
- Liu A, Tran L, Becket E, Lee K, Chinn L, Park E, Tran K, Miller JH (2010) Antibiotic sensitivity profiles determined with an Escherichia coli gene knockout collection: generating an antibiotic bar code. Antimicrob Agents Chemother 54: 1393–1403CASPubMedWeb of Science®Google Scholar
- Loewe S (1928) Die quantitativen Probleme der Pharmakologie. Ergebnisse der Physiol. 27: 47–187Google Scholar
- Moreau F, Desroy N, Genevard JM, Vongsouthi V, Gerusz V, Le Fralliec G, Oliveira C, Floquet S, Denis A, Escaich S, Wolf K, Busemann M, Aschenbrenner A (2008) Discovery of new Gram-negative antivirulence drugs: structure and properties of novel E. coli WaaC inhibitors. Bioorg Med Chem Lett 18: 4022–4026CASPubMedWeb of Science®Google Scholar
- Munck C, Gumpert HK, Wallin AIN, Wang HH, Sommer MO (2014) Prediction of resistance development against drug combinations by collateral responses to component drugs. Sci Transl Med 6: 262ra156CASPubMedWeb of Science®Google Scholar
- Nichols RJ, Sen S, Choo YJ, Beltrao P, Zietek M, Chaba R, Lee S, Kazmierczak KM, Lee KJ, Wong A, Shales M, Lovett S, Winkler ME, Krogan NJ, Typas A, Gross CA (2011) Phenotypic landscape of a bacterial cell. Cell 144: 143–156CASPubMedWeb of Science®Google Scholar
- Nikaido H, Vaara M (1985) Molecular basis of bacterial outer membrane permeability. Microbiol Rev 49: 1–32CASPubMedWeb of Science®Google Scholar
- Ocampo PS, Lázár V, Papp B, Arnoldini M, Abel Zur Wiesch P, Busa-Fekete R, Fekete G, Pál C, Ackermann M, Bonhoeffer S (2014) Antagonism is prevalent between bacteriostatic and bactericidal antibiotics. Antimicrob Agents Chemother 58: 4573–4582CASPubMedWeb of Science®Google Scholar
- Pena-Miller R, Laehnemann D, Jansen G, Fuentes-Hernandez A, Rosenstiel P, Schulenburg H, Beardmore R (2013) When the most potent combination of antibiotics selects for the greatest bacterial load: the smile-frown transition. PLoS Biol 11: e1001540CASPubMedWeb of Science®Google Scholar
- Pillai SK, Moellering RC, Eliopoulos GM (2005) Antimicrobial combinations. In Antibiotics in Laboratory Medicine, V Lorian (ed), pp 365–440. Philadelphia: Lippincott Williams and WilkinsGoogle Scholar
- Rodriguez de Evgrafov M, Gumpert H, Munck C, Thomsen TT, Sommer MO. (2015) Collateral resistance and sensitivity modulate evolution of high-level resistance to drug combination treatment in Staphylococcus aureus. Mol Biol Evol 32: 1175–1185CASPubMedWeb of Science®Google Scholar
- Rothfield L, Pearlman-Kothencz M (1969) Synthesis and assembly of bacterial membrane components: a lipopolysaccharide-phospholipid-protein complex excreted by living bacteria. J Mol Biol 44: 477–492CASPubMedWeb of Science®Google Scholar
- Schneider DA, Gaal T, Gourse RL (2002) NTP-sensing by rRNA promoters in Escherichia coli is direct. Proc Natl Acad Sci USA 99: 8602–8607CASPubMedWeb of Science®Google Scholar
- Scott M, Hwa T (2011) Bacterial growth laws and their applications. Curr Opin Biotechnol 22: 559–565CASPubMedWeb of Science®Google Scholar
- Tamae C, Liu A, Kim K, Sitz D, Hong J, Becket E, Bui A, Solaimani P, Tran KP, Yang H, Miller JH (2008) Determination of antibiotic hypersensitivity among 4,000 single-gene-knockout mutants of Escherichia coli. J Bacteriol 190: 5981–5988CASPubMedWeb of Science®Google Scholar
- Tavazoie S, Hughes JD, Campbell MJ, Cho RJ, Church GM (1999) Systematic determination of genetic network architecture. Nat Genet 22: 281–285CASPubMedWeb of Science®Google Scholar
- Waldstein E, Sharon R, Ben-Ishali R (1974) Role of ATP in excision repair of ultraviolet radiation damage in Escherichia coli. Proc Natl Acad Sci USA 71: 2651–2654CASPubMedWeb of Science®Google Scholar
- Wood K, Nishida S, Sontag ED, Cluzel P (2012) Mechanism-independent method for predicting response to multidrug combinations in bacteria. Proc Natl Acad Sci USA 109: 12254–12259CASPubMedWeb of Science®Google Scholar
- Wood KB, Wood KC, Nishida S, Cluzel P (2014) Uncovering scaling laws to infer multidrug response of resistant microbes and cancer cells. Cell Rep 6: 1073–1084CASPubMedWeb of Science®Google Scholar
- Yeh P, Tschumi AI, Kishony R (2006) Functional classification of drugs by properties of their pairwise interactions. Nat Genet 38: 489–494CASPubMedWeb of Science®Google Scholar
- Zaslaver A, Bren A, Ronen M, Itzkovitz S (2006) A comprehensive library of fluorescent transcriptional reporters for Escherichia coli. Nat Methods 3: 623–628CASPubMedWeb of Science®Google Scholar