Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, April 2005, p. 1483-1494, Vol. 49, No. 4
0066-4804/05/$08.00+0 doi:10.1128/AAC.49.4.1483-1494.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Department of Biology, Emory University, Atlanta, Georgia,1 Department of Microbiology, Ramón y Cajal University Hospital, Madrid, Spain2
Received 16 July 2004/ Returned for modification 7 October 2004/ Accepted 28 December 2004
|
|
|---|
|
|
|---|
With the aid of mathematical models, we develop these three hypotheses to account for the decrease in bacterial mortality: (i) antibiotic decay (since the efficacy of antibiotics increases with their concentration, a decay in the effective concentration of the antibiotics leads to a decrease in mortality), (ii) inherited resistance (an ascent of genetically resistant mutants decreases the mortality of the total bacterial population), and (iii) phenotypic tolerance (the bacterial population, though genetically homogeneous, is physiologically heterogeneous with respect to its susceptibility, which during antibiotic exposure leads to an enrichment of the fraction of phenotypically tolerant bacteria and thus a decrease in the overall bacterial mortality).
We investigate this phenomenon in vitro, using Escherichia coli and five classes of antibiotics (ciprofloxacin, ampicillin, rifampin, streptomycin, and tetracycline).
We present evidence against the first two of these hypotheses and in support of the third hypothesis (phenotypic tolerance). Using a mathematical model of antibiotic treatment, we demonstrate that phenotypic tolerance can impair the efficacy of treatment.
|
|
|---|
We also used variants of E. coli CAB1 resistant to antibiotics. We used spontaneous mutants and transconjugants. Spontaneous mutants that were resistant to nalidixic acid, rifampin, and streptomycin were obtained by plating approximately 2 x 109 E. coli CAB1 from cultures grown overnight onto LB agar with either 64 µg of nalidixic acid per ml, 25 µg of rifampin per ml, or 40 µg of streptomycin per ml. Transconjugants were generated by conjugative transfer of plasmid R1 (Cm, Km, Am, Sm/Sp) and/or plasmid RK2 (Am, Tc) into E. coli CAB1. To test whether the time-kill results observed in this study were restricted to one strain (E. coli CAB1), we repeated these experiments with E. coli K-12 MG1655.
Culture and sampling media. Bacteria were grown at 37°C with aeration and shaking (200 rpm) in Luria-Bertani (LB) broth (Difco) with 10 ml of broth in 50-ml Erlenmeyer flasks or 50 to 100 ml of broth in 250-ml flasks. Total cell densities and the densities of antibiotic-resistant mutants and transconjugants were estimated from CFU data by diluting the bacteria (in 0.85% saline) and plating on LB agar and LB agar with the selecting antibiotics, respectively.
Antibiotics. Ciprofloxacin was from Mediatech, Inc. (Herndon, Va.), and ampicillin, rifampin, streptomycin, and tetracycline were from Sigma-Aldrich (St. Louis, Mo.). The stock solution of rifampin (10 mg/ml) was dissolved in methanol, and the stock solution of tetracycline (25 mg/ml) was dissolved in 50% ethanol. The stock solutions of streptomycin, ampicillin, and ciprofloxacin were dissolved in sterile distilled water (at 40, 25, and 10 mg/ml, respectively). The stock solutions of rifampin, ciprofloxacin, and tetracycline were maintained at 20°C, and those of ampicillin and streptomycin were maintained at 4°C.
Antibiotic time-kill curves. Cultures of E. coli CAB1 grown overnight in LB broth were diluted in 10 ml of fresh warm (37°C) LB broth and incubated for 2 h to initiate exponential growth. These cultures were grown to a final density of approximately 2 x 106 cells per ml before antibiotics were added (dissolved in LB broth). The cultures were incubated with shaking at 200 rpm at 37°C and sampled to estimate the number of CFU per milliliter every 15 min for the first 2 h, every 30 min for the next 2 h, and at 4, 6, and 24 h. Time-kill curves were obtained in each case for two antibiotic concentrations, selected to illustrate the possibility of different concentration-dependent killing dynamics.
Effective concentrations of antibiotics. Bioassays were used to monitor the changes in the effective concentration of antibiotics. At specific times after the start of the time-kill experiments, the cultures were passed through 0.45-µm-pore-size membrane filters (Tyffryn; Pall Corporation) to remove the bacteria. Aliquots of 0.1 ml of exponentially growing monocultures of E. coli CAB1 or mixed cultures of exponentially growing E. coli CAB1 containing low frequencies of E. coli CAB1 or C600 resistant to the antibiotic in the medium were added to the supernatant and incubated for 2 h at 37°C. Samples were then taken and plated on LB agar with and without the selecting antibiotic to estimate the density of viable cells and the relative frequency of the resistant strain. In the experiments with monocultures, the decline in viable cell density at 2 h was measured in fresh antibiotic-containing medium and in the filtrates of used medium with antibiotics (with new bacteria added). Changes in bactericidal activity were statistically evaluated by paired Wilcoxon signed-rank test using the R program (a P value of less than 0.05 was considered statistically significant). In the experiments with mixed cultures (sensitive and resistant strains), the decline in the viable density of the total population and the increase in the frequency of resistance in the mixed cultures were used as measures of the activity of antibiotics, a resistance competition assay (6).
Tests for resistant mutants. (i) Test for inherited resistance. Single colonies were tested for mutation to resistance by streaking on LB agar plates containing resistance breakpoint concentrations of streptomycin (40 µg/ml), rifampin (25 µg/ml), and nalidixic acid (30 µg/ml) (i.e., concentrations allowing growth of only resistant mutants).
(ii) MIC estimation by Etest. Aliquots of 0.2 ml of single-colony cultures grown overnight were spread on LB agar onto which Etest strips (AB Biodisk, Solna, Sweden) were placed. Estimates of the MICs for these cultures were obtained the next day using the Etest criteria.
Reexposure experiments. At different times in the course of time-kill experiments, 10-ml samples of the cultures were passed through 0.45-µm-pore-size membrane filters and washed with fresh LB broth, and the bacteria on the filters were suspended in LB broth containing different concentrations of either the same or different antibiotics. Samples of the antibiotic-exposed bacteria were obtained during the exponential decline phase (at 30 min) and later at 1.5 and 3 h. For controls, 10-ml samples from exponentially growing cultures with no history of antibiotic exposure were filtered and resuspended in fresh medium with antibiotic. The viable density of cells (in reexposed and control cultures) was determined at time zero when the drug was added and at the end of each reexposure period, at 2 or 3 h. In the experiments in which the cells were reexposed to other antibiotics (the generality experiments) or higher concentrations of the same antibiotic (the robustness experiments), the samples were removed and reexposed at 4 h.
|
|
|---|
![]() View larger version (28K): [in a new window] |
FIG. 1. Time-kill curves. Changes in estimated densities of viable E. coli CAB1 (in CFU per milliliter) exposed to ciprofloxacin, ampicillin, rifampin, streptomycin, and tetracycline. The concentrations of the antibiotics (in micrograms per milliliter) are shown in the symbol key boxes. The MICs of these antibiotics for E. coli CAB1 are listed in Table 3.
|
|
View this table: [in a new window] |
TABLE 1. Densities of viable bacteria at three different times
|
To ascertain whether these results were restricted to E. coli CAB1, we performed similar time-kill experiments with E. coli K-12 (MG1655). We obtained analogous results with a decline in the rate of mortality with persistence of a fraction of viable bacteria (data not shown).
Three hypotheses and their predictions by computer simulations. We consider three hypotheses to account for the decline in the rate of mortality with persistence of a fraction of viable bacteria in the time-kill experiments presented in Fig. 1: (i) decay in the concentration of the antibiotic (as time proceeds, the antibiotic becomes decreasingly bactericidal), (ii) inherited resistance (the population initially includes a minority of genetically resistant cells which increases in number and frequency as the sensitive population dies), and (iii) phenotypic tolerance (there is variation among the bacteria in their susceptibility to the different antibiotics, and as time proceeds, the relative frequency of the more tolerant members of the population increases). To illustrate how these different hypotheses can account for the observations and the predictions they make, we solved differential equation models of each of these processes numerically, the details of which are included in the Appendix. These simulations were programmed by using the Berkeley Madonna program and can be obtained from the www.eclf.net website. For the analysis of the properties of these models, we use parameters that provide initial mortality rates similar to those observed for E. coli CAB1 in vitro (Fig. 1). The parameters are listed in the legend to Fig. 7.
![]() View larger version (29K): [in a new window] |
FIG. 7. Computer simulation of antibiotic treatment with and without tolerance. Changes in the density of viable bacteria with different levels of tolerance and treatment regimens are shown. For pharmacokinetics, the initial concentration of antibiotic is a0 and decays exponentially with a rate parameter, , of 0.5 h1. At either 8- or 12-h intervals 10 µg of that antibiotic per ml is added. For pharmacodynamics, we assume a Hill function for the pharmacodynamics of the sensitive subpopulation µS(a) with Smax = 1.0 h1, Smin = 10.0 h1, zMIC = 4 µg/ml, and the Hill coefficient = 1 (see Appendix). The maximum replication rate of the tolerant subpopulation is Tmax = 0.5 h1, and the tolerant population is unaffected by the antibiotic. With these parameters, in the absence of antibiotic decay, the time-kill curve for the tolerant population is identical to that in Fig. 2c. (a) Effects of different values of the tolerance parameter, f, and antibiotic dose (10 µg/ml) every 8 h. (b) Effects of different values of the tolerance parameter, f, and antibiotic dose (10 µg/ml) every 12 h.
|
(zMIC) = 0], and the bacterial population begins to grow. With the parameters used to generate the simulations in Fig. 2a, the density of viable bacteria at 24 h exceeds 2 x 109 cells/ml. One testable prediction of the antibiotic decay hypothesis is that when the rate of mortality levels off or the population begins to increase, the antibiotic would no longer be bactericidal.
![]() View larger version (23K): [in a new window] |
FIG. 2. Computer simulations. Predicted results of time-kill experiments anticipated under the three hypotheses are shown. (A) Antibiotic decay. Rate of antibiotic decay = 0.4 h1. Hill function parameters for the pharmacodynamics follow: maximum and minimum growth rates, max = 1.0 h1; min = 10.0 h1; zMIC = 4 µg/ml; Hill coefficient = 1. The right x axis is the antibiotic concentration. (B) Inherited resistance. Growth rate of sensitive cells rS = 1; death rate of sensitive cells mS = 6; growth rate of resistant cells rR = 1. The initial density of resistant cells is 1 cell per ml. (C) Phenotypic tolerance. Growth rate of sensitive cells rS = 1; growth rate of tolerant cells rT = 0.5; mortality of sensitive cells mS = 10; death rate of tolerant cells mT = 0; tolerance parameter f = 0.001. The descriptions of these models and the equations used in these simulations are presented in the Appendix.
|
In Fig. 2c, we present the results of simulations of the phenotypic tolerance model. In this model, we assume that the bacteria are of two states, one sensitive to killing by the antibiotic and one tolerant. In contrast to genetic resistance, these states are not inherited; when a new bacterial cell arises, the probability that it will be tolerant (f) is the same for sensitive and tolerant mother cells. We refer to f as the tolerance parameter.
In the parameters used in the simulations presented in Fig. 2c, we assume that the rate at which tolerant cells are produced, f, is substantially less than the rate of production of sensitive cells (1 f). After addition of antibiotics, the sensitive cells decline exponentially, tolerant cells become enriched, and the rate of mortality of the total population, N, declines. The number of viable cells present when the time-kill curve levels off depends on the tolerance parameter, f. The larger the value of f is, the greater the density of phenotypic tolerant cells will be when the total population reaches a low density. Another prediction of the phenotypic tolerance model is that samples taken during the exponential decline phase would be dominated by sensitive cells while those taken after the rate of mortality has declined would be dominated by tolerant cells.
Tests of the hypotheses by in vitro experiments. While we believe these three classes of hypotheses are exhaustive, they are not mutually exclusive and could be operating simultaneously.
(i) Recovery of viable cells at 24 h.
In the majority of cases, when cultures containing approximately 106 E. coli CAB1 per ml are exposed to the five antibiotics, the number of viable cells at 24 h is less than 105 CFU/ml (Table 1). This relatively low density of cells at 24 h is inconsistent with the predictions of both the antibiotic decay hypothesis and the inherited resistance hypothesis. Both of these hypotheses predict bacterial densities at 24 h approaching those of stationary-phase cultures (
109), either because of a decline in bactericidal activity followed by growth of sensitive bacteria or because of growth of genetically resistant cells (assuming these cells do not have too large of a growth disadvantage). Finally, in accord with the inherited resistance hypothesis, it should be possible to recover cells with inherited resistance soon after the total cell density has stopped decreasing. While rifampin-, nalidixic acid-, and streptomycin-resistant mutants were recovered from some exposed cultures, this was a relatively rare outcome in the time-kill experiments and could not account for the majority of these results.
(ii) Effective antibiotic concentration. We compared the rate of killing of exponentially growing E. coli CAB1 in fresh LB broth with antibiotics (fresh medium) with the rate of killing in 5-h-old cell-free filtrates of previously exposed cultures (filtrates). The results of these experiments are presented in Table 2. For ciprofloxacin, ampicillin, rifampin, and streptomycin, there was either no decline or a small decline in the bactericidal activity, while for tetracycline, the bactericidal activity appeared to be greater in the filtrates than in the fresh medium. None of these differences were statistically significant. We interpret these results to suggest that decay in the effective concentration of the antibiotic could account for some of the decline in the rate of mortality of bacteria exposed to ciprofloxacin and streptomycin but could not account for the decline in mortality of bacteria exposed to ampicillin, rifampin, and tetracycline. After 5 h, all antibiotics tested remained bactericidal.
|
View this table: [in a new window] |
TABLE 2. Changes in the bactericidal activity of the antibioticsa
|
(iii) Time-kill curves of cultures derived from exposed cells. If inherited resistance were responsible for the decline in the rate of mortality and the survival of bacteria in the presence of antibiotics, the shape of time-kill curves of cultures derived from the survivors of time-kill experiments should be different from those of cultures derived from bacteria that have not been exposed to antibiotics. More specifically, these cultures should be dominated by resistant cells, and there should be little if any decline in the density of the total bacterial population in the presence of antibiotics. However, this was not the case. The trajectories of the time-kill curves produced from cultures of bacteria isolated from bacteria with prior exposure to antibiotics were effectively the same as those derived from cells that had no prior exposure to the antibiotic (Fig. 3).
![]() View larger version (28K): [in a new window] |
FIG. 3. Test for inherited resistance in whole cultures. (A) Exposure of whole cultures of naïve E. coli CAB1. (B) Reexposure of cultures surviving 3 h of exposure. The surviving cells from panel A at 3 h were caught on filters and subcultured overnight in the absence of antibiotics before they were reexposed to the same drug. The antibiotic concentrations are shown in micrograms per milliliter. Antibiotic abbreviations: cipro, ciprofloxacin; amp, ampicillin, rif, rifampin; strep, streptomycin; tet, tetracycline.
|
(v) Test for inherited resistance by Etest. To ascertain whether there may be inherited resistance to low concentrations of the antibiotics, we used Etest to estimate the MICs of antibiotics for unexposed bacteria and bacteria obtained after 4 h of antibiotic exposure from four independent time-kill experiments. The MICs from these Etests are listed in Table 3. The results of these experiments were also inconsistent with the inherited resistance hypothesis. The MICs of antibiotics for cultures derived from single colonies of bacteria that survived exposure to antibiotics were the same as those derived from unexposed bacteria. Most colonies tested had either no increase or small increases in MIC, but all remained well below the resistance breakpoints.
|
View this table: [in a new window] |
TABLE 3. MICs of antibiotics (by Etest) for bacteria that survived exposure to antibioticsa
|
![]() ![]() View larger version (29K): [in a new window] |
FIG. 4. Test for phenotypic tolerance. Reexposure of cells previously exposed to antibiotics for different periods of time (30 min, 1.5 h, and 3 h). Controls were exponentially growing cells with no history of antibiotic exposure. The antibiotic concentrations were as follows: 0.0625 µg/ml for ciprofloxacin and 24 µg/ml for ampicillin, rifampin, streptomycin, and tetracycline. Survival after the reexposure is depicted here. The values represent the log10 ratios of the density of surviving cells after 2 h of exposure and their densities prior to reexposure (surv dens./init dens.). Each value represents the average ratio ± standard error (error bar) of three experiments.
|
![]() ![]() View larger version (55K): [in a new window] |
FIG. 5. Test for robustness. E. coli CAB1 was exposed to a low concentration of one antibiotic and then reexposed to increasing concentrations of the same antibiotic. Survival after the reexposure is depicted here. Antibiotic concentrations at the first exposure were 0.0625 µg/ml for ciprofloxacin and 24 µg/ml for ampicillin, rifampin, streptomycin, and tetracycline. The bars represent the ratios of the density of surviving cells (after the second exposure) relative to the initial density at reexposure (surv dens/init dens.). The time of reexposure is shown on the x axis. The boxed control values are the ratios of the density of surviving naïve cells after 2 h of exposure relative to the initial density. The concentrations of the antibiotics (in micrograms per milliliter) at the second exposure are noted in the symbol key boxes.
|
![]() ![]() View larger version (31K): [in a new window] |
FIG. 6. Cross-tolerance to other antibiotics. After 4 h of exposure to the antibiotics indicated at the x axis, the bacteria were reexposed to five different antibiotics, indicated at the top of each graphs. Two repetitions were performed, both of which are shown. The bars represent the ratios of the density of surviving cells (after the second exposure) relative to the initial density at reexposure. The control values are the ratios of the density of surviving naïve cells after 2 h of exposure (to the antibiotic at the top of the graph) relative to the initial density. The concentrations of the antibiotics at the first and second exposure were as follows: 0.0625 µg/ml for ciprofloxacin, and 24 µg/ml for ampicillin, rifampin, streptomycin, and tetracycline. The controls in this experiment are the same as those in Fig. 5. Abbreviations: CIP, ciprofloxacin; AMP, ampicillin; RIF, rifampin; STR, streptomycin; TET, tetracycline.
|
Potential clinical implications of phenotypic tolerance. To explore the possible clinical implications of phenotypic tolerance, we developed a model of antibiotic treatment similar to that employed in references 25 and 30. For the pharmacokinetics of this model, we assume that the effective concentration of the antibiotic declines exponentially between doses and administered a fixed dose of antibiotic every 8 or 12 h. For the pharmacodynamics, we combine a Hill function (30) with the phenotypic tolerance model. We assume the relationship between the concentration of the antibiotic and the rate of growth (or death) of the sensitive subpopulation is defined by a Hill function, while the rate of growth (or death) of the tolerant subpopulation is not affected by the antibiotic concentration. As a consequence of antibiotic treatment, the tolerant subpopulation is enriched. However, since the concentration of the antibiotic is continually changing, the intensity of enrichment for the tolerant subpopulation also varies with time. A more detailed description of this phenotypic tolerance model is in the Appendix.
To illustrate the effects of tolerance, we monitor the changes in the density of bacteria over the course of 5 days with 8-h antibiotic dosing regimens in the absence of tolerance and with two values of the tolerance parameter f (Fig. 7a). In Fig. 7b, we consider the effects of a 12-h dosing regimen in the absence of tolerance and with two values of f. The results of this analysis clearly illustrate that the microbiological course of treatment can be profoundly affected by phenotypic tolerance. A treatment regimen that rapidly clears an infection in the absence of tolerance can fail to do so if there is phenotypic tolerance.
|
|
|---|
With the aid of mathematical models, we developed and explored the properties of three hypotheses to account for this phenomenon: (i) antibiotic decay, (ii) inherited resistance, and (iii) phenotypic tolerance. Although these hypotheses are not mutually exclusive, our experiments show that only the phenotypic tolerance hypothesis provides a sufficient explanation for the decreasing rate of mortality of bacteria exposed to all five antibiotics and for the survival of substantial fractions of their populations.
The antibiotic decay hypothesis cannot explain the decrease in the rate of mortality of E. coli CAB1 exposed to ampicillin, rifampin, and tetracycline, because the effective concentrations of these antibiotics did not decline with time. A small reduction in bactericidal activity was observed for ciprofloxacin and streptomycin during the period of exposure; therefore, we cannot rule out a contribution of antibiotic decay to the decrease in bacterial mortality for these two drugs. However, even after 5 h, these two antibiotics remain strongly bactericidal, and the decay in their effective concentration is not sufficient to account for the observed decrease in mortality or the survival of bacteria in medium with these antibiotics. While mutants resistant to streptomycin and rifampin arose in some of our time-kill experiments, these phenomena cannot be explained by inherited resistance. A similar rate of decrease in the rate of mortality and survival of the exposed bacteria was also observed in the absence of inherited resistance or partial resistance to streptomycin and rifampin as well as the other three antibiotics. Moreover, when we repeated these time-kill experiments with bacteria surviving exposure and regrown in antibiotic-free medium, the dynamics were the same as those observed with naïve cells. This would not be anticipated if there was evolved resistance in the previously exposed populations.
To explore the predictions of the phenotypic tolerance hypothesis, we constructed a mathematical model of this process in which we assumed that the bacterial population consists of a phenotypically sensitive subpopulation and a phenotypically tolerant subpopulation. We assume that during cell division there is a probability that a daughter cell will be tolerant and one minus that probability that it will be sensitive. This probability is the same for dividing tolerant and sensitive bacteria. In this model, the decrease in mortality of the total bacterial population is caused by an enrichment of phenotypically tolerant bacteria in the presence of antibiotics. The results of our experiments with all five antibiotics were fully consistent with the two predictions of this model. (i) Bacteria isolated after 30 min of exposure were more susceptible to killing by antibiotics upon reexposure than bacteria isolated after 1.5 to 3 h of exposure. (ii) As noted above, when the bacteria surviving exposure were cultured in the absence of antibiotics, their susceptibility was the same as that of a naïve population.
Our experiments also suggest that the tolerant state is robust and general and can extend to higher concentrations of the antibiotic than to the concentration to which they were exposed and to other antibiotics. Except for bacteria grown in higher concentrations of tetracycline, bacteria present later in time-kill experiments, during the tolerant phase, are substantially less susceptible to higher concentrations of the antibiotic to which they were exposed than naïve bacteria. Also, depending on the antibiotic the bacteria were exposed to, bacteria isolated during the tolerant phase are less susceptible to antibiotics of other classes than naïve bacteria are.
We believe that the phenotypic tolerance observed in these studies with exponentially growing bacteria is different from the antibiotic tolerance previously described for drugs that inhibit cell wall synthesis (32-36). This phenomenon appears to be due to the exposed bacteria being at stationary phase and for that reason refractory to antibiotics. On the other hand, the phenotypic tolerance observed here with E. coli CAB1 and these five different classes of antibiotics may be the same as what has been referred to as adaptive resistance in the studies of Pseudomonas aeruginosa exposed to aminoglycosides (2, 3, 10, 11, 20, 39). In those investigations, a lower sensitivity to the bactericidal effects of these antibiotics was observed for antibiotic-free cultures derived from the survivors of exposed cells relative to that for unexposed cells. The level of this sensitivity increased and eventually became similar to that of naïve cells with an increase in the amount of time these preexposed bacteria were maintained in antibiotic-free medium. It is also possible that the small-colony variants isolated among the survivors of Salmonella enterica serovar Typhimurium and Staphylococcus aureus exposed to aminoglycosides may be a manifestation of phenotypic tolerance as considered here (31, 37). When recultured in the absence of antibiotics, these small-colony variants give rise to normal colonies which are as sensitive to the antibiotics as naïve cells.
In this study, we focused on the population dynamics and the potential clinical implications of phenotypic tolerance and have not explored the cellular, physiological, or molecular processes responsible. The mathematical model of phenotypic tolerance we developed here is not mechanistic, as it does not specify the biological mechanisms responsible for either physiological variation in sensitivity to antibiotics within an isogenic population or for the transition between the susceptible and tolerant states. Moreover, this model is only a minimalist caricature of the phenotypic tolerance phenomenon in that it assumes that the bacterial population consists of only two distinct homogeneous subpopulations, susceptible and tolerant. In reality, there would almost certainly be a continuous distribution in the extent of susceptibility. It should be noted, however, that the same qualitative outcome would occur with a continuous distribution; exposure to antibiotics would enrich the less-susceptible cells and thereby reduce the overall susceptibility to the lethal effects of these antibiotics. This variation in susceptibility could be a reflection of the variation among the cells with respect to their stage in the replication cycle. It could also be a consequence of stochasticity (noise) in the cellular, physical, and biochemical processes affecting gene expression (14). Even for a single bacterium, temporal variation in cellular processes can occur in the absence of external stimuli (23). Also contributing to this variation in susceptibility or possibly reflecting it are the physiological (metabolic) and morphological changes observed in bacteria exposed to antibiotics (13, 24, 26, 38).
Two recent studies have investigated the phenomenon of phenotypic tolerance. Miller et al. (28) showed that induction of the SOS response by ampicillin can protect E. coli against the bactericidal effects of ampicillin by disturbing cell wall synthesis and growth. Balaban et al. (1) showed that there is variation in growth rate in a genetically homogeneous population and that the cells with lower growth rates (before exposure) survived exposure to ampicillin better than normally growing cells. While the killing kinetics of different classes of antibiotics vary due to due to the different mechanisms of action (30), exposure of bacteria to any one of these classes of antibiotics results in survival of a minority of the population. However, the mechanisms of tolerance (persistence) demonstrated for ampicillin (1, 28) might not account for phenotypic tolerance for all classes of antibiotics. Our observation of cross-tolerance in some combinations of antibiotic exposures and reexposures suggests that tolerance mechanisms could be shared by some classes of antibiotics but probably not by all. Massive changes in gene expression leading to changes in the syntheses of proteins of metabolic and stress response pathways and cell division during exposure of E. coli to ampicillin and ofloxacin have recently been observed (19). A number of these alterations in the gene expression levels were shared between bacteria exposed to ampicillin and ofloxacin. Keren et al. (21, 22) suggest that random fluctuations in gene expression are responsible for the formation of specialized persister cells. These cells overexpress toxins which, at high concentrations, can shut down the function of multiple drug targets (such as protein synthesis and DNA replication) and thereby prevent the antibiotic from corrupting the function of the target molecules. The HipA and RelE toxins have been shown to be directly involved in multidrug tolerance in E. coli (22).
Does phenotypic tolerance matter clinically? Although we did not investigate the consequences of phenotypic tolerance in a clinical setting, we assessed its potential microbiological consequences by using a mathematical model of antibiotic treatment. In this model, we combine the phenotypic tolerance model with the pharmacodynamics and pharmacokinetics of antibiotics. Our simulations of a periodic dosing antibiotic treatment regimen (Fig. 7) illustrates that phenotypic tolerance could have a substantial effect on the rate of clearance of a bacterial population and could actually prevent clearance in a situation where treatment would be effective in the absence of tolerance. We believe that these in vitro experimental results and this a priori consideration of the microbiological consequences of phenotypic tolerance in patients treated with antibiotics are sufficiently compelling to warrant the study of this process in more detail and explore its implications for the rational design of antibiotic treatment regimens.
|
|
|---|
![]() | (1) |
(a) is the pharmacodynamic function (30) which describes the relationship between the bacterial growth coefficient,
, (which can take on positive or negative values) and the antibiotic concentration, a:
![]() | (2) |
max is the maximum growth rate coefficient in the absence of antibiotics,
min is the minimum net growth rate coefficient at the highest possible antibiotic concentrations, zMIC is the pharmacodynamic MIC of the antibiotic [i.e., the antibiotic concentration that causes
(zMIC) = 0], and
is the Hill coefficient, a measure of how fast the net bacterial growth rate coefficient,
, changes around the zMIC.
We assume that the antibiotic concentration, a, decays exponentially over time at a rate
per hour. The rates of change in the concentration of the antibiotic are given by the differential equation
![]() | (3) |
![]() | (4) |
ln(a0/zMIC), and starts to increase after that. The rate of growth of N converges to
max (Fig. 2a). (ii) Inherited resistance. In this model, we divide the bacterial population into two subpopulations, one that is genetically sensitive and one that is genetically resistant to the antibiotic with densities of S and R CFU per milliliter, respectively. We assume that, in the absence of antibiotics, sensitive and resistant bacteria grow at rates rS and rR, respectively. In the presence of antibiotics, sensitive bacteria die at a rate mS, while there is no mortality of resistant bacteria.
The rates of change in the densities (CFU per milliliter) of sensitive and resistant cells in this population are given by the differential equations
![]() | (5) |
![]() | (6) |
(iii) Phenotypic tolerance.
In this model, we assume that bacterial cells can exist in two physiological states, sensitive and resistant, with densities S and T CFU per milliliter, respectively. In contrast to the inherited resistance model, these states are not inherited; when a new bacterial cell arises, the probability, f, that it will be tolerant is independent of the state of the cell from when it came. We assume that these phenotypically sensitive and tolerant bacteria replicate at rates rS and rT, respectively, and are killed at rates mS and mT, respectively. The rates of change in the densities of these two phenotypes are given by
![]() | (7) |
![]() | (8) |
Mathematical model to assess the clinical implications of phenotypic tolerance. To illustrate the impact of phenotypic tolerance on treatment efficacy, we use a model that combines the phenotypic tolerance model with models of pharmacodynamics and pharmacokinetics of the antibiotic.
In this model, we assume bacteria of two physiological states, sensitive and tolerant, with densities S and T CFU per milliliter. For the pharmacokinetics, we assume that a dose of a0 is administered every 8 or 12 h and that the antibiotic concentration declines exponentially at a rate,
.
In the above phenotypic tolerance model (equations 7 and 8), the replication rate of sensitive and tolerant subpopulations are not independent. Hence, to combine the pharmacodynamics with the pharmacokinetics, we need to separate the pharmacodynamic function of these two populations into replication and mortality rates. To describe the pharmacodynamic relation between the mortality rate constant, µS, of the phenotypically sensitive bacteria and the antibiotic concentration, a, we use a Hill function:
![]() | (9) |
max
(a). As in the above phenotypic tolerance model, we assume that the tolerant bacteria are not killed by the antibiotic. With these definitions and assumptions, the rates of change in the concentration of the antibiotic between doses and densities of the sensitive and tolerant populations are given by
![]() | (10) |
![]() | (11) |
![]() | (12) |
This research was funded in part by grants from the U.S. National Institutes of Health (GM33782 and AI40662), The Wellcome Trust (IPRAVE project) (to B.R.L.), and the Spanish Pneumococcal Infection Study Network (G03/103) (to F.B.).
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»