Abstract
Background
Efflux is a widespread mechanism of reversible drug resistance in bacteria that can be triggered by environmental stressors, including many classes of drugs. While such chemicals when used alone are typically toxic to the cell, they can also induce the efflux of a broad range of agents and may therefore prove beneficial to cells in the presence of multiple stressors. The cellular response to a combination of such chemical stressors may be governed by a tradeoff between the fitness costs due to drug toxicity and benefits mediated by inducible systems. Unfortunately, disentangling the costbenefit interplay using measurements of bacterial growth in response to the competing effects of the drugs is not possible without the support of a theoretical framework.
Results
Here, we use the wellstudied multiple antibiotic resistance (MAR) system in E. coli to experimentally characterize the tradeoff between drug toxicity (“cost”) and druginduced resistance (“benefit”) mediated by efflux pumps. Specifically, we show that the combined effects of a MARinducing drug and an antibiotic are governed by a superposition of cost and benefit functions that govern these tradeoffs. We find that this superposition holds for all drug concentrations, and it therefore allows us to describe the full dose–response diagram for a drug pair using simpler cost and benefit functions. Moreover, this framework predicts the existence of optimal growth at a nontrivial concentration of inducer. We demonstrate that optimal growth does not coincide with maximum induction of the mar promoter, but instead results from the interplay between drug toxicity and mar induction. Finally, we derived and experimentally validated a general phase diagram highlighting the role of these opposing effects in shaping the interaction between two drugs.
Conclusions
Our analysis provides a quantitative description of the MAR system and highlights the tradeoff between inducible resistance and the toxicity of the inducing agent in a multicomponent environment. The results provide a predictive framework for the combined effects of drug toxicity and induction of the MAR system that are usually masked by bulk measurements of bacterial growth. The framework may also be useful for identifying optimal growth conditions in more general systems where combinations of environmental cues contribute to both transient resistance and toxicity.
Background
The resistance of bacteria to antibiotics has prompted intense scientific research in the last several decades because it directly underlies the clinical treatment of infections [1]. While a large number of studies have focused on mutationdriven resistance, recent attention has also shifted to transient, or “inducible”, drug resistance taking place on much shorter timescales [29]. This transient resistance does not rely on mutations, but can be induced by a large class of chemicals commonly found in drugs (e.g. antibiotics and painkillers) and food preservatives. These chemicals are typically toxic to the cell when used alone, but they can also induce resistance to a broad range of agents. Consequently, they may prove beneficial to cells in the presence of multiple stressors. The net effect of a combination of chemical stressors can therefore be counterintuitive, because it is governed by the interplay between inducible resistance and drug toxicity. Such situations may arise, for example, in the human digestive tract, where bacteria face a cocktail of diverse chemical cues. The combined effects of multiple stressors, in general, have been studied for many decades in hopes of optimizing the clinical efficacy of combinatorial therapies [1012]. More recently, the effects of drug interactions on the evolution of irreversible (mutationdriven) drug resistance have also been recognized [1316]. Drugs that interact synergistically to produce a strong toxic effect can accelerate the acquisition of mutations conferring drug resistance [14]. On the other hand, antagonistic drug pairs produce a weaker toxic effect but can slow the acquisition of resistance [13,15]. These results demonstrate an inherent tradeoff between the toxicity of the drug combination and its potential to facilitate drug resistance [17]. They also raise an interesting question: do tradeoffs between drug toxicity and resistance also play a role in inducible drug resistance? Furthermore, when cells are exposed a combination of toxic agents that potentially induce transient resistance, how are these tradeoffs related to synergy or antagonism between the given agents?
To address these questions, here we study inducible resistance mediated by the MAR (multiple antibiotic resistance) system. The MAR system, present in many bacterial species, consists of an operon that confers effluxmediated [1821] resistance to a broad range of antibiotics and can be activated by a host of chemical agents, including analgesics and food preservatives [25,79,2225]. For example, in E. coli, high concentrations (5 mM) of salicylate, a wellknown analgesic, lead to an increase in the minimum inhibitory concentration (MIC) of several antibiotics, and this effect has been partially attributed to the MAR system [4,26,27]. However, MIC alone masks the distinct effects of salicylate toxicity, antibiotic toxicity, and induction of the MAR system (the “benefit” of salicylate) on the growth, because it provides only an aggregate measure of these effects. Therefore, the underlying nature of the costbenefit interplay remains hidden. Here we combine experiments with a simple phenomenological model to fully characterize a range of effects resulting from the tradeoff between drug toxicity and the benefit mediated by efflux pumps.
Results
Experimental characterization of salicylateinduced antibiotic resistance
To study the tradeoff between drug toxicity and induced resistance, we first measured the effects of salicylate and two protein synthesis inhibitors, chloramphenicol and tetracycline, on cell growth. As expected, the effect of each drug alone is to slow cell growth as concentration increases. To quantitatively characterize this effect, we define growth cost as the reduction in growth rate of cells treated with one drug relative to the growth rate of untreated cells. In general, we find that cost functions are welldescribed by Hill functions with K_{i}, the concentration of drug i at which cost is half maximal, and n, the Hill coefficient (Figure 1a). This mathematical form is consistent with standard dose–response models [10]. In the presence of a high concentration of chloramphenicol or tetracycline, however, we find that adding salicylate can increase growth (Figure 1b, Additional file 1: Figure S1). This effect was previously reported for a single concentration of salicylate [26]; here we provide the entire nonmonotonic dependence of this suppressive effect on concentration, allowing us to tease apart contributions of salicylate’s cost from its benefit (below). We also performed detailed quantitative measurements of another suppressive interaction between an antibiotic, chloramphenicol, but with a different inducer, sodium benzoate (Additional file 1: Figure S1).
Additional file 1. Supplemental Figures, Supplemental Methods, and Supplemental Notes [[10,23,3642]].
Format: PDF Size: 2MB Download file
This file can be viewed with: Adobe Acrobat Reader
Figure 1. Single and Multiple Drug Growth Costs.a. Growth cost is defined as the reduction in growth rate of cells treated with a drug relative to the growth rate of untreated cells. Solid lines, best fits to Hill functions h(x) = x^{n}/(K_{i}^{n} + x^{n}), with K_{i} equal to the concentration of drug i (i = S for salicylate, i = A for both antibiotics) at which growth is inhibited by 50%. K_{A} = 1.80 (1.66, 1.93) μg/mL, n = 1.97 (1.67, 2.26) for chloramphenicol; K_{A} = 0.41 (0.37, 0.45) μg/mL, n = 1.90 (1.61, 2.19) for tetracycline; and K_{S} = 6.06 (5.50, 6.62) mM, n = 1 for salicylate. 95% confidence intervals from nonlinear least squares fitting in parentheses. Error bars are standard deviation of four replicates. b. In the presence of high concentrations of chloramphenicol (Cm), growth is maximal for a nonzero concentration of salicylate (upper insets: OD time series). Contours of constant growth in concentration space indicate suppression. Dashed line, contours of 45% and 55% maximum growth.
Phenomenological costbenefit model for MARinduced drug resistance
To quantitatively model the interplay between inducer cost and benefit in a twodrug environment that includes an inducer, we assume that the effects of the two drugs, in the absence of MAR induction, are independent in the sense that the relative growth rate of cells in the presence of both drugs (g_{SA}) equals the product of individual relative growth rates of cells in the presence of each drug alone (g_{SA} = g_{S} g_{A}). This assumption, known as Bliss independence [11], provides exact results for simple situations, such as a single enzymatic reaction disrupted by mutually non exclusive inhibitors [10], in which case the effects of the drugs are defined in terms of reaction fluxes. However, the assumption of Bliss independence often fails for more complex systems, and it is therefore primarily used as a phenomenological null model for the physiological effects of noninteracting drugs on cell proliferation or growth [10,14,2830].
Here, we extend the concept of Bliss independence to include drug interactions mediated by the induction of the MAR system. Specifically, we assume that the presence of one drug (S) provides a fitness benefit by reducing the effective concentration of the second drug (A), thereby coupling the effects of the drugs. We incorporate this rescaling and rewrite the model in terms of the combined growth cost, C_{SA} (by definition 1 g_{SA}) to arrive at
where f(A) and g(S) are the growth costs of drugs alone (Figure 1a). Equation 1 expresses Bliss independence in terms of growth costs, and also generalizes it to include an Sdependent reduction of A to A_{eff}. Therefore, the approximate additivity of drug costs implied by Bliss independence is modified to an approximate additivity of effective costs (Figure 2a). The model assumes that the presence of inducer reduces the effective concentration of drug A and adds an additional cost, but otherwise does not affect growth rate. This assumption considerably limits the spectrum of possible cellular responses to the drug pair, because drug S can only change the effective concentration of drug A, but will not change the shape of its cost function. The model can be readily extended to include two drugs which both act as inducers, in which case S is also changed to S_{eff}. However, because the antibiotics used are poor inducers of the MAR system (Additional file 1: Figure S7), we neglect their effects on the concentration S and limit this study to the directional model implied by Equation [1].
Figure 2. CostBenefit Model for Interactions between MAR inducers and Antibiotics.a. The combined growth cost of two drugs (S and A) is assumed to be Bliss independent, as long as the effective concentration of one drug is reduced by the presence of the second drug. Thus, the combined cost of both drugs is approximately equal to the sum of effective growth costs, f and g, for each individual drug. The reduction of each drug's growth cost depends on the concentration of the other drug, which in turn dictates the growth benefit provided. b. mar promoter activity, defined as the steady state timeaveraged Fluorescence/OD times relative growth rate, increases as a function of salicylate concentration. Error bars are +/− one estimated standard deviation. Inset: mar fluorescence concentration (YFP fluorescence/optical density) time series for concentrations of salicylate ranging from 0 mM (purple) to 6 mM (blue). Data points are means over two replicates. Promoter activity is estimated by averaging fluorescence concentration in the steady state and multiplying by relative growth rate.
A simple model, which assumes that the concentration of antibiotic A is reduced by the induction of efflux pumps (Additional file 1), suggests a form for A_{eff}:
We call the function β(S) the inducible benefit, and it contains all quantitative information about the resistance mechanisms induced by drug S. To estimate the inducible benefit, β(S), for the MAR system, we measured the activity of the mar promoter as a function of S (Figure 2b; Additional file 1). MAR promoter activity could dictate the functional form for β(S) if, for example, we assume it is proportional to the number of efflux pumps produced in response to salicylate (Additional file 1). We scale this promoter activity by an adjustable parameter β_{max}, so that , where K_{ind} = 0.8 is the concentration of inducer that yields half maximal promoter activity (Figure 2b). The adjustable parameter β_{max} links mar activity induced by drug S to a phenotypic response (resistance) to a second drug A. This parameter will be specific to drug A, and it provides a measure of the efficiency with which MAR induction eliminates A from the cell.
The costbenefit theory assumes that the effects of the inducing drug and the antibiotic are independent, up to a reduction in the concentration of the antibiotic. The aim of our model is not to achieve a microscopic theoretical description of the system, but rather to provide a minimal phenomenological model that quantitatively captures the measured behavior. The success or failure of the model must therefore be determined experimentally. Specifically, our model does not attempt to elucidate the microscopic variables governing the multidrug effects—which would require dozens, if not hundreds, of microscopic parameters—but rather posits that a simple relationship should exist between cell growth in the presence of one drug (which is a function of the cell’s internal state) and cell growth in the presence of two drugs (which is a function of an entirely different intracellular state). Specifically, the model requires equality between cell growth in the presence of A and cell growth in the presence of the combination S and A' (with A'> A), once we account for the costs of drug S. To directly verify this hypothesis, we can rearrange equation 1 to express this concept as
The left hand side represents the "adjusted" multidrug cost, once the toxic effects of S (cost) have been removed. The model implies that this adjusted cost of the drug combination, as a function of A_{eff}, is functionally equivalent to the cost function f(x) for a single drug. Hence, if one removes from the multidrug costs (C_{SA}) the effects of g(S) according to the lefthand side of Equation 3 and then properly chooses the concentration reduction A→A_{eff} (determined by the single parameter β_{max}), all twodrug data from a given drug combination should collapse to the single curve determined by f(x). Note that because a Hill function describes costs for many individual drugs as well as slices through many drug combination effect surfaces[10], collapsing the multidrug cost curves from a given drug pair (Figure 3a) to a common Hill form can be simply achieved using independent parameters for each salicylate concentration. By contrast, equation 3 suggests that a nontrivial collapse of all curves for a given drug pair requires only a single parameter β_{max}. Once this parameter is determined, the beneficial effects of different concentrations of salicylate are linked by the mar promoter activity curve. We verified this constraint (Figure 3b) by fitting the single free parameter, β_{max}, using twodrug cost curves for salicylate and chloramphenicol (β_{max} = 1.15 +/− 0.15) and salicylate and tetracycline (β_{max} = 6.04 +/− 0.16). In all cases, the model provides an excellent fit to the data (R^{2} = 0.98; see Additional file 1: Figure S2 for a direct comparison of growth rates from experiment and model). Interestingly, this analysis demonstrates that a superposition of cost and benefit functions quantitatively describes the combined effects of salicylate and an antibiotic. Surprisingly, we find that the benefit function depends only on the inducer concentration and is dictated by the mar induction curve.
Figure 3. CostBenefit Analysis Quantitatively Describes MultiDrug Growth Cost Functions.a. Multidrug growth cost curves for salicylate and tetracycline (top) and salicylate and chloramphenicol (bottom) are increasing functions of antibiotic concentration, but also depend on the concentration of salicylate. Each growth cost curve represents growth in increasing concentrations of antibiotic at a single salicylate concentration, ranging from 0 to 7 mM. Black lines, ([salicylate] = 0 mM), which correspond to the single drug growth cost functions. Data points, means of four replicates. Error bars, +/− one sample standard deviation. b. The costbenefit model requires rescaled versions of the multidrug cost functions to have the same mathematical form as single drug cost functions (Equation 3). To directly test this assumption, we rescaled the multidrug cost functions in panel a by removing the contribution from salicylate toxicity (cost) to each multidrug growth cost function. Then, we rescaled the concentration of each antibiotic (chloramphenicol or tetracycline) using the single parameter β_{max}, which accounts for the salicylateinduced reduction of intracellular antibiotic concentration (benefit). We trivially exploit the common Hill form of the single drug costs functions to show data from both drug pairs on a single plot, achieved by replacing drug concentration with effective toxicity, defined as ([A]/K_{A})^{n}, where K_{A} and n characterize the single drug growth cost function for the antibiotics (see Figure 1). We use the salicylate induction curve of the mar promoter (Figure 2b) as the benefit function β(S), for both drug interactions, with only the overall scale of the benefit function (β_{max}) specific to the two antibiotics. See also Additional file 1: Figure S2 for a direct comparison of growth rates from experiments and models.
Transition from salicylate as toxic to salicylate as beneficial
Our results uncover an interesting range of cellular behavior that emerges from the tradeoffs of salicylate toxicity and the simultaneous induction of multidrug resistance. First, for each concentration of antibiotic, we find that growth is maximal at a single salicylate concentration S*. In addition, we see an apparent transition between two regimes—one where the presence of the inducing drug is harmful, and another where it is beneficial—as A eclipses a threshold A_{crit} (Figure 4a). That is, S* becomes nonzero as A crosses A_{crit}. In other words, the cell benefits from the presence of salicylate, but only when the antibiotic environment is sufficiently toxic to offset the inherent toxicity of the inducer. Interestingly, we find that mar promoter activity saturates at S ≈ 4 mM (Figure 2b), which is greater than S* for all concentrations of antibiotics used. Therefore, maximum suppression at S* results from a nontrivial interplay between MAR expression and inducer cost, not merely a maximum in MAR expression.
Figure 4. Phase Diagram for SalicylateChloramphenicol Interactions.a. The concentration of salicylate, S* (in mM), at which growth is maximal depends on the external concentration of antibiotic. A transition occurs between the low antibiotic regime, where the presence of salicylate is harmful, to the high antibiotic where the presence of salicylate increases the growth of cells. Red, tetracylcline; Blue, chloramphenicol. Concentrations for both drugs measured in μg/mL. Curves, numerical calculation from model; circles, estimates of S* from experiments at single concentrations of A. Maxima were estimated using cubic spline interpolation between data points at different concentrations of salcilyate. Error bars, concentrations corresponding to +/− 1% of maximum. b. Apparent MIC, defined as the minimum concentration of antibiotic (red, tetracycline; blue, chloramphenicol; both in units of μg/mL) at which growth is reduced to δ = 0.5. While the quantitative results depend on the precise definition of MIC, the qualitative features do not depend on the choice of δ. c. Costbenefit theory predicts general properties of drug interactions as a function of the maximum inducible benefit β_{max} and the concentration of inducer, S. Solid line, phase boundary between synergistic and antagonistic drug interactions. Dashed line, phase boundary between antagonistic and suppressive interactions. In general, a higher value of the β_{max} increases the antagonism between drug pairs at a given concentration of drug S. Parameters characterizing individual drugs include K_{S} and K_{ind}, which characterize, respectively, the cost of drug S and the corresponding induction of resistance systems (e.g. efflux pumps), and K_{A} and n, which characterize the cost of drug A. By contrast, β_{max} couples the individual effects of two drugs. Insets, heat maps of twodimensional growth surfaces for 3 cell strains in the presence of Salicylate (Sal) and Chloramphenicol (Cm); top: β_{max} = 1.15 (WT cells; suppressive), β_{max} = 0.19 (mar mutant, antagonistic), and β_{max} = −0.15 (tolC mutant, synergistic). See also Additional file 1: Figure S3.
Contributions of inducer cost and benefit to apparent MIC’s of antibiotics
Given our characterization of salicylateinduced multidrug resistance, it is straightforward to calculate the apparent MIC of an antibiotic in the presence of any concentration of salicylate. Here, we define the apparent MIC to be the minimum concentration of drug A at which relative growth is reduced to a value of δ. For small concentrations of salicylate, the apparent MIC’s of both tetracycline and chloramphenicol increase dramatically (Figure 4b). However, at higher concentrations, the toxicity of salicylate begins to overwhelm its potential benefits, eventually leading to a decrease in apparent MIC. It is therefore clear that MIC reflects contributions of both the cost and benefits of the inducing drug. By contrast, the interdependence of inducible resistance and mar promoter activity can be described using only a single parameter, β_{max}, which provides a quantitative measure of inducible benefit that is not masked by the effects of inducer cost.
Phase diagram for MARmediated drug interactions
The interactions between salicylate and tetracycline, salicylate and chloramphenicol, and sodium benzoate and tetracycline are all suppressive, but this class of models describes a range of interactions ranging from synergistic to suppressive, based on the interplay of induction benefit and drug toxicity. Using Equations 1 and 2, we quantitatively determine a phase diagram (Figure 4, Additional file 1: Figure S3) that specifies how a drug interaction depends on the tradeoffs between inducible benefit (described by the parameters K_{ind} and β_{max}), the cost of inducer (K_{S}), and the steepness of the antibiotic doseresponse curve (n). Specifically, the cellular response to the first drug (S) alone includes both the growth cost of the drug (characterized by K_{S}) and the induction of MAR system (characterized by K_{ind}). The dynamics of the efflux pumps and their specificity for drug A determine the concentration reduction (A to A_{eff}), which is governed by β_{max}. The phase diagram demonstrates that properties of the drugs alone (n, K_{S}, K_{ind}) determine the level of drug coupling, contained in β_{max}, required to achieve antagonism or suppression. Drugs that strongly induce growth benefit and have low associated cost (K_{ind} <<K_{S}) are always suppressive. By contrast, high cost inducers (K_{ind} >> K_{S}) can never be suppressive because the cost of induction is too high (Figure 4c).
The inducible benefit parameter (β_{max} = 1.15) characterizing the salicylate and chloramphenicol combination in wild type cells is far above the suppressiveantagonistic boundary (β_{max} > > 2K_{ind}/(K_{1}n) = 0.45), and the interaction is clearly suppressive for all concentrations of salicylate. To explore different regions of the phase diagram experimentally, we measured the effects of two different mutations on the suppressive drug interaction between salicylate and chloramphenicol. The first mutant, ΔtolC [31], lacks the protein TolC required for efflux pumping [32]. While we expect that the ΔtolC deletion will partially suppress the benefit induced by salicylate, we cannot predict how the costs of the drugs and, therefore, the precise nature of the drug interaction, will be altered. Experimentally, this mutant showed decreased resistance to chloramphenicol (Figure 4a, 4b), corresponding to a rescaling of the single drug cost (i.e. K_{A} smaller than in wild type), but the cost of salicylate remains unchanged (K_{S} approximately same as wild type). β_{max} was measured to be −0.15 +/− 0.01. According to the phase diagram (Figure 4), the interaction should therefore be weakly synergistic, and in fact the contours of constant growth appear slightly concave (lower inset). The suppressive drug interaction has been eliminated because the benefit associated with salicylate has been decreased while its cost remains unchanged, resulting in significantly less antagonistic interaction.
Since the suppressive interaction between salicylate and chloramphenicol is partially associated with the synthesis of AcrABTolC efflux pumps, we also hypothesized that cells with constitutive mar promoter activity would disrupt the costbenefit interplay required for drug suppression. We selected such a mutant, herecalled a tet mutant, by growing cells in 1 μg/mL tetracycline for 48 hours. We observed that the MAR system is a common target for mutations that confer resistance to tetracycline (see also [23]) (Additional file 1: Figure S6). In terms of our costbenefit analysis, these cells synthesize resistance systems (benefit) without the associated cost of an inducing drug. To verify the molecular basis of resistance in the tet mutant, we sequenced the full genome and found only a single point mutation in the α5 region of the MAR repressor marR, which is linked to dimer formation and subsequent binding to the mar promoter [33]. In the salicylatechloramphenicol combination, the tet mutant showed increased resistance to chloramphenicol (K_{A} is larger than in wildtype) and a nearadditive drug interaction (β_{max} = 0.19 +/− 0.03) between salicylate and chloramphenicol (Figure 4). Because the cost of salicylate is unchanged, the elimination of suppression suggests that the mutation has blocked the inducible benefit of salicylate (see also Additional file 1: Figure S5 for similar results with a different drug pair). Biologically, the inactivation of the MarR repressor no longer requires the presence of salicylate, and therefore the cost of the drug exceeds its associated inducible benefit. Interestingly, we also found that the mar activity in the mutant is higher in the absence of salicylate than the maximum mar activity induced by even high concentrations of salicylate in the wild type strain. However, adding salicylate further increases mar activity in the mutant (Additional file 1: Figure S7A). Moreover, the functional dependence of the induction on salicylate is similar to that of the wildtype (Additional file 1: Figure S7A, inset). This result is not intuitive because the beneficial effects of salicylate on the growth of the mutant in the presence of chloramphenicol are quite small (β_{max} is approximately 17% of that of the wild type, and the antagonism between the drugs is markedly decreased).
Discussion
Efflux is a widespread cellular defense mechanism with clinical consequences for antimicrobial treatments. Costbenefit principles provide a convenient tool for analyzing these ubiquitous systems. We have shown that such a framework provides a quantitative description of drug interactions associated with the classic MAR system, which mediates an inducible resistance to many different toxic agents. The model allows us to quantify the relative importance of the induction mediating reversible resistance and the toxicity of the inducing agent. In particular, we observe a transition between two regimes, one where salicylate is harmful and one where it is beneficial, as the concentration of antibiotic crosses a threshold A_{crit}. By disentangling the contributions of inducer cost and benefit, the analysis also provides a quantitative measure (β_{max}) of the reversible (MARmediated) resistance to each antibiotic used. This measure complements previous MIC measurements [26], which provide information about the net effect of both inducer costs and benefits but that mask the relative contributions of each. We demonstrate that it is precisely this balance between cost and benefit that determines whether the antimicrobial properties of multiple agents are synergistic, antagonistic, or suppressive. Therefore, in addition to providing a new quantitative model of the MAR system, we also lay the groundwork for a link between the shortterm tradeoffs of inducible resistance, which determine the type of interaction between pairs of toxic agents, and the acquisition of irreversible resistance mediated by these drug interactions.
This framework may be applicable to other systems where environmental signals can be toxic but also induce transient resistance. For example, cancer cells can activate prosurvival stress responses upon treatment with chemotherapies or radiation[34], and subpopulations of tumor cells can access new transcriptional/differentiation states to become reversibly drug resistant[35]. Our analysis implements effluxmediated resistance using a concentration reduction A → A_{eff} for the effective concentration of drug A, and this form could be used to describe other common mechanisms of resistance, including enzymatic decay and target modification. As a consequence, the costbenefit framework developed here may prove useful for modeling and understanding a wide range of systems governed by the inherent tradeoffs of inducible resistance.
Conclusions
Our analysis provides a quantitative description of the MAR system and demonstrates that optimal cell growth in the presence of an inducer and an antibiotic does not coincide with maximum induction of the mar promoter, but instead results from the interplay between drug toxicity and mar induction. In addition, we show that this interplay determines whether the drugs are synergistic, antagonistic, or suppressive. The framework may also be extensible to more general systems where combinations of environmental cues contribute to both transient resistance and toxicity.
Methods
Strains
We used the E. coli Frag1B strain, which is wildtype for the AcrAB multiple drug efflux pump system, and also the wildtype E.coli BW25113 strain. The ΔtolC mutant is strain JW55031 from the Keio collection.
Drugs
Drug solutions were made from solid stocks (sodium salicylate, SigmaAldrich; chloramphenicol, MP Biomedicals; tetracycline hydrochloride, Acros Organics; doxycycline hyclate, Sigma Alderich; kanamycin, Fisher; ciprofloxacin, SigmaAlderich). All antibiotic stock solutions were stored in the dark at −20°C. All drugs were thawed (when necessary) and diluted in sterilized Lennox LB broth (Fisher) for experimental use.
Growth conditions and drug treatments
We inoculated LB media with a single colony and grew the cells overnight (12 h at 30°C with shaking at 200 rpm). Following overnight growth, stationary phase cells were diluted (600–1000 fold) in LB media containing 0.2% glucose and grown for an additional 3 h (30°C with shaking at 200 rpm). Following the initial dilution period, we transferred cells to 96 well microplates (round bottom, polystyrene, Corning) and set up a twodimensional matrix of drug concentrations in each 96well microplate (165 μl media per well). For the remainder of the experiment (4–8 h), cells were grown at 30°C with shaking at 1200 rpm on a Heidolph Titramax vibrating microplate shaker (Brinkmann). Plates were wrapped in aluminum foil to minimize evaporation. A600 (absorbance at 600 nm, proportional to optical density OD) and YFP fluorescence were measured at 15–25 min intervals for 4–8 h using a Wallac Victor2 1420 Multilabel Counter (PerkinElmer).
Growth rate estimation and MIC determination
Growth rates were determined by fitting multiple regions of growth curves (A_{600} vs. time, Additional file 1: Figure S4) in early exponential phase to an exponential function (MATLAB 7.6.0 curve fitting toolbox, The Mathworks) using a variable length (45–150 minutes) sliding window. Growth rate corresponded to the window position associated with best fit to an exponential function, determined by a minimum root mean squared error (RMSE). For each experimental condition, the final growth rate was the sample mean of 4 replicates (See also Additional file 1). Growth rate contours were estimated from experimental data using the Matlab csaps function to interpolate between data points (typically 96) spaced roughly evenly throughout the two dimensional space of drug concentrations. MIC was taken to be the drug concentration(s) at which growth was inhibited by 50%.
Drug resistant mutants
We diluted liquid cultures of wildtype cells (in stationary phase) 1000fold into 96 individual 150 μL cultures on a single microplate, with each well supplemented with 1.0 μg/mL tetracycline. After approximately 48 h of growth at 30°C with shaking, optical density and fluorescence were measured from each well. Samples from randomly selected cultures exhibiting high OD and fluorescence (Additional file 1: Figure S6) were frozen in 15% glycerol at −80°C. Frozen cultures were used to streak LB plates, and all experiments (including sequencing) were performed on cultures grown from a single colony isolated from these plates.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
KW and PC designed research, analyzed data, and wrote the paper. KW performed all experiments. Both authors read and approved the final manuscript.
Acknowledgements
We thank L. Bruneaux, E. Balleza and A. Subramaniam for technical guidance, C. Guet for advice and for providing the Frag1B pZS*2 MARvenus strain, and all members of the Cluzel lab for many helpful discussions. We also thank K. Dave for editorial advice and assistance and Kris Wood for comments on the manuscript. This work was supported in part by the NIH P50 award P50GM08189202 to the University of Chicago and the NSF Postdoctoral Fellowship 0805462 (to K.W.).
References

Levy SB, Marshall B: Antibacterial resistance worldwide: causes, challenges and responses.
Nat Med 2004, 10(12 Suppl):S122S129. PubMed Abstract  Publisher Full Text

Alekshun MN, Levy SB: The mar regulon: multiple resistance to antibiotics and other toxic chemicals.
Trends Microbiol 1999, 7(10):410413. PubMed Abstract  Publisher Full Text

Cohen SP, Hachler H, Levy SB: Genetic and functional analysis of the multiple antibiotic resistance (mar) locus in Escherichia coli.
J Bacteriol 1993, 175(5):14841492. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Cohen SP, Levy SB, Foulds J, Rosner JL: Salicylate induction of antibiotic resistance in Escherichia coli: activation of the mar operon and a mar independent pathway.
J Bacteriol 1993, 175(24):78567862. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Critzer FJ, Dsouza DH, Golden DA: Transcription analysis of stx1, marA, and eaeA genes in Escherichia coli O157:H7 treated with sodium benzoate.
J Food Prot 2008, 71(7):14691474. PubMed Abstract  Publisher Full Text

Herbert D, Paramasivan CN, Venkatesan P, Kubendiran G, Prabhakar R, Mitchison DA: Bactericidal action of ofloxacin, sulbactamampicillin, rifampin, and isoniazid on logarithmic and stationaryphase cultures of Mycobacterium tuberculosis.
Antimicrob Agents Chemother 1996, 40(10):22962299. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Potenski CJ, Gandhi M, Matthews KR: Exposure of Salmonella Enteritidis to chlorine or food preservatives decreases [corrected] susceptibility to antibiotics.
FEMS Microbiol Lett 2003, 220(2):181186. PubMed Abstract  Publisher Full Text

Tavio MM, Vila J, Perilli M, Casanas LT, Macia L, Amicosante G, de Jimenez Anta MT: Enhanced active efflux, repression of porin synthesis and development of Mar phenotype by diazepam in two enterobacteria strains.
J Med Microbiol 2004, 53(Pt 11):11191122. PubMed Abstract  Publisher Full Text

Viveiros M, Jesus A, Brito M, Leandro C, Martins M, Ordway D, Molnar AM, Molnar J, Amaral L: Inducement and reversal of tetracycline resistance in Escherichia coli K12 and expression of proton gradientdependent multidrug efflux pump genes.
Antimicrob Agents Chemother 2005, 49(8):35783582. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Greco WR, Bravo G, Parsons JC: The search for synergy: a critical review from a response surface perspective.
Pharmacol Rev 1995, 47(2):331385. PubMed Abstract  Publisher Full Text

Bliss CI: The calculation of microbial assays.
Bacteriol Rev 1956, 20(4):243258. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Loewe S: The problem of synergism and antagonism of combined drugs.
Arzneimittelforschung 1953, 3(6):285290. PubMed Abstract

Chait R, Craney A, Kishony R: Antibiotic interactions that select against resistance.
Nature 2007, 446(7136):668671. PubMed Abstract  Publisher Full Text

Hegreness M, Shoresh N, Damian D, Hartl D, Kishony R: Accelerated evolution of resistance in multidrug environments.
Proc Natl Acad Sci USA 2008, 105(37):1397713981. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Michel JB, Yeh PJ, Chait R, Moellering RC, Kishony R: Drug interactions modulate the potential for evolution of resistance.
Proc Natl Acad Sci USA 2008, 105(39):1491814923. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Chait R, Shrestha S, Shah AK, Michel JB, Kishony R: A differential drug screen for compounds that select against antibiotic resistance.
PLoS One 2010, 5(12):e15179. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Torella JP, Chait R, Kishony R: Optimal drug synergy in antimicrobial treatments.
PLoS Comput Biol 2010, 6(6):e1000796. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Levy SB: Active efflux mechanisms for antimicrobial resistance.
Antimicrob Agents Chemother 1992, 36(4):695703. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Levy SB: Active efflux, a common mechanism for biocide and antibiotic resistance.
J Appl Microbiol 2002, 92(Suppl):65S71S. PubMed Abstract  Publisher Full Text

Okusu H, Ma D, Nikaido H: AcrAB efflux pump plays a major role in the antibiotic resistance phenotype of Escherichia coli multipleantibioticresistance (Mar) mutants.
J Bacteriol 1996, 178(1):306308. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Wall ME, Markowitz DA, Rosner JL, Martin RG: Model of transcriptional activation by MarA in Escherichia coli.
PLoS Comput Biol 2009, 5(12):e1000614. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Cohen SP, Hooper DC, Wolfson JS, Souza KS, McMurry LM, Levy SB: Endogenous active efflux of norfloxacin in susceptible Escherichia coli.
Antimicrob Agents Chemother 1988, 32(8):11871191. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Cohen SP, McMurry LM, Hooper DC, Wolfson JS, Levy SB: Crossresistance to fluoroquinolones in multipleantibioticresistant (Mar) Escherichia coli selected by tetracycline or chloramphenicol: decreased drug accumulation associated with membrane changes in addition to OmpF reduction.
Antimicrob Agents Chemother 1989, 33(8):13181325. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Cohen SP, McMurry LM, Levy SB: marA locus causes decreased expression of OmpF porin in multipleantibioticresistant (Mar) mutants of Escherichia coli.
J Bacteriol 1988, 170(12):54165422. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

McMurry LM, George AM, Levy SB: Active efflux of chloramphenicol in susceptible Escherichia coli strains and in multipleantibioticresistant (Mar) mutants.
Antimicrob Agents Chemother 1994, 38(3):542546. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Rosner JL: Nonheritable resistance to chloramphenicol and other antibiotics induced by salicylates and other chemotactic repellents in Escherichia coli K12.
Proc Natl Acad Sci USA 1985, 82(24):87718774. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Price CT, Lee IR, Gustafson JE: The effects of salicylate on bacteria.
Int J Biochem Cell Biol 2000, 32(10):10291043. PubMed Abstract  Publisher Full Text

Fitzgerald JB, Schoeberl B, Nielsen UB, Sorger PK: Systems biology and combination therapy in the quest for clinical efficacy.
Nat Chem Biol 2006, 2(9):458466. PubMed Abstract  Publisher Full Text

Yeh P, Tschumi AI, Kishony R: Functional classification of drugs by properties of their pairwise interactions.
Nat Genet 2006, 38(4):489494. PubMed Abstract  Publisher Full Text

Lehar J, Zimmermann GR, Krueger AS, Molnar RA, Ledell JT, Heilbut AM, Short GF, Giusti LC, Nolan GP, Magid OA, et al.: Chemical combination effects predict connectivity in biological systems.
Mol Syst Biol 2007, 3:80. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, Baba M, Datsenko KA, Tomita M, Wanner BL, Mori H: Construction of Escherichia coli K12 inframe, singlegene knockout mutants: the Keio collection.
Mol Syst Biol 2006, 2:2006 0008. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Fralick JA: Evidence that TolC is required for functioning of the Mar/AcrAB efflux pump of Escherichia coli.
J Bacteriol 1996, 178(19):58035805. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Alekshun MN, Levy SB, Mealy TR, Seaton BA, Head JF: The crystal structure of MarR, a regulator of multiple antibiotic resistance, at 2.3 A resolution.
Nat Struct Biol 2001, 8(8):710714. PubMed Abstract  Publisher Full Text

Li F, Sethi G: Targeting transcription factor NFkappaB to overcome chemoresistance and radioresistance in cancer therapy.
Biochim Biophys Acta 2010, 1805(2):167180. PubMed Abstract  Publisher Full Text

Sharma SV, Lee DY, Li B, Quinlan MP, Takahashi F, Maheswaran S, McDermott U, Azizian N, Zou L, Fischbach MA, et al.: A chromatinmediated reversible drugtolerant state in cancer cell subpopulations.
Cell 2010, 141(1):6980. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

George AM, Levy SB: Amplifiable resistance to tetracycline, chloramphenicol, and other antibiotics in Escherichia coli: involvement of a nonplasmiddetermined efflux of tetracycline.
J Bacteriol 1983, 155:531540. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Alekshun MN, Levy SB: Regulation of chromosomally mediated multiple antibiotic resistance: the mar regulon.
Antimicrob Agents Chemother 1997, 41:20672075. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Maneewannakul K, Levy S: Identification of mar mutants among quinoloneresistant clinical isolates of Escherichia coli.
Antimicrob Agents Chemother 1996, 40:16951698. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Le TT, Emonet T, Harlepp S, Guet CC, Cluzel P: Dynamical determinants of druginducible gene expression in a single bacterium.
Biophys J 2006, 90:33153321. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Lutz R, Bujard H: Independent and tight regulation of transcriptional units in Escherichia coli via the LacR/O, the TetR/O and AraC/I1I2 regulatory elements.
Nucleic Acids Res 1997, 25:12031210. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Elf J, Nilsson K, Tenson T, Ehrenberg M: Bistable bacterial growth rate in response to antibiotics with low membrane permeability.
Phys Rev Lett 2006, 97:258104. PubMed Abstract  Publisher Full Text