Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, June 2005, p. 2189-2199, Vol. 49, No. 6
0066-4804/05/$08.00+0 doi:10.1128/AAC.49.6.2189-2199.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Division of Infectious Diseases, Department of Internal Medicine, College of Medicine,1 Division of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, Lexington, Kentucky2
Received 13 September 2004/ Returned for modification 20 January 2005/ Accepted 5 February 2005
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
We fit three pharmacodynamic models to the data generated from our experiments and found that a two-population model with unique growth (g) and killing (k) rate constants for the ciprofloxacin-susceptible (gS and kS, respectively) and -resistant (gR and kR, respectively) subpopulations best described the initial killing and subsequent regrowth patterns observed in the in vitro system (5). The model correctly described the enrichment of subpopulations with low-level resistance in the parent cultures. It confirmed the experimental observations that there was no clear relationship between the TMSW and the enrichment of resistant subpopulations. The model indicated that resistance depended upon the selection of low-level-resistant minor subpopulations in bacterial cultures. It predicted that resistant subpopulations would not emerge when a low-density culture with a low probability of mutants was exposed to a simulated clinical dosing regimen (400 mg every 12 h) or when a high-density culture with a higher probability of mutants was exposed to the high concentrations of an experimental regimen designed to rapidly eradicate grlA mutants, followed by the lower concentrations of the clinical regimen. The validity of these predictions was confirmed with in vitro system experiments.
In this study we extend our observations to levofloxacin, which has been shown to select resistant S. aureus variants less frequently than ciprofloxacin (10, 14, 24, 34). This was done by modeling the effect of simulated clinical and experimental levofloxacin regimens on two S. aureus strains and their grlA mutants in the in vitro system. The experiments provided additional information about the succession of mutations that occur in resistance loci as bacteria evolve and allowed us to test the robustness of our previously developed pharmacodynamic model.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Antimicrobial agent. Analytical-grade levofloxacin (R. W. Johnson Pharmaceutical Research Institute, Raritan, N.J.) was used to prepare stock solutions according to established guidelines (30). The stock solutions were frozen (80°C) in aliquots and used within 30 days.
Susceptibility tests. Levofloxacin MICs and minimal bactericidal concentrations (MBCs) were determined by the broth microdilution method (29, 30) in triplicate for each organism before exposure to levofloxacin and for organisms recovered at 0, 24, 48, and 96 h during in vitro system experiments. MICs were also measured in the presence of 20 µg/ml reserpine (Sigma, St. Louis, Mo.), a competitive inhibitor of NorA, to screen for potential efflux of levofloxacin (20).
Mutant prevention concentration determinations.
Levofloxacin MPCs for MRSA 8043 and MRSA 8282 were determined using a previously described agar dilution method (4). Concentrated bacterial suspensions (
1010 CFU/ml) were prepared, and 200 µl samples were applied to Mueller-Hinton agar (Difco Laboratories, Detroit, Mich.) containing levofloxacin concentrations of 0, 0.25, 0.5, 0.75, 1, 1.25, 1.5, 1.75, 2, 2.25, and 2.5 µg/ml. The MPC for each strain was defined as the lowest concentration at which no bacteria were detected following 48 h of incubation at 37°C. Nucleotide sequence analysis of the quinolone resistance-determining regions of grlA/B and gyrA/B (see below) was performed using DNA extracted from the most-resistant variants that appeared.
In vitro system experiments.
A two-compartment hollow-fiber in vitro system that allowed the study strains to be exposed to fluctuating concentrations of levofloxacin was used. The main features and general operation of the system and the method for preparing standardized inocula for the experiments have been described previously (4). Log-phase cultures (
107 CFU/ml) were exposed to a series of monoexponential levofloxacin pharmacokinetic profiles for 96 h. Three simulations were designed to reproduce the concentration-time profiles of clinical levofloxacin doses of 750 mg, 500 mg, and 250 mg administered as intermittent 1-hour infusions every 24 h. The simulated central and peripheral compartment concentrations were intended to mimic the total levofloxacin concentrations observed in plasma and skin blister fluid, respectively, of adult patients with normal renal function (6-9). An additional experimental regimen (125 mg every 24 h) was simulated to produce levofloxacin concentrations below the MPCs of the study strains for the entire duration of the experiment. Growth control experiments for each strain were conducted over 36 h. Every in vitro system experiment was performed in duplicate.
Pharmacokinetic analysis. Samples collected from the central and peripheral compartments of the in vitro system during each dosing regimen experiment were analyzed for levofloxacin concentrations using a previously validated ion pair high-performance liquid chromatography method (4). Standard curves were linear (R2 > 0.999) over the concentration range tested (0.1 to 50 µg/ml). The within-run relative standard deviations (SD) (n = 5) for levofloxacin quality control concentrations of 25, 2.5, and 0.25 µg/ml were 1.51%, 1.38%, and 0.97%. Between-run relative SD (n = 16) for the same quality control concentrations were 1.88%, 1.82%, and 1.23%. The reliable lower limit of quantification (defined as percent relative SD and percent deviation from nominal concentration that was <15%) was 0.06 µg/ml.
Central and peripheral compartment concentration-time data for levofloxacin were analyzed by compartmental and noncompartmental methods (5, 13). A one-compartment pharmacokinetic model, with first-order elimination from the central compartment, was fitted to levofloxacin concentrations in the central compartment using a nonlinear regression program (WinNonlin version 4.0; Pharsight Corporation, Mountain View, Calif.). All pharmacokinetic data were weighted by the inverse of squared predicted concentrations based on the precision of the high-performance liquid chromatography assay.
Pharmacodynamic analysis. Viable counts were determined for samples collected from the peripheral compartment of the in vitro system using previously described drop and membrane filter count methods (4). Experiments performed to assess the effect of levofloxacin carryover indicated that dilutions of at least 100-fold were necessary to prevent a significant (P < 0.05, paired t test) lowering of colony counts in samples containing levofloxacin concentrations of 0.5 to 16 µg/ml compared to those of drug-free controls. When bacterial counts were unreliable because of levofloxacin carryover, the membrane filter count method was used. Previous experiments established reliable lower limits of detection (coefficients of variation < 20%) of 3.0 x 102, 1.0 x 102, and 7.0 x 101 CFU/ml for filtered samples of 100, 250, and 500 µl, respectively (4).
Resistant subpopulation analysis.
Viable cells were counted to quantify resistant bacterial subpopulations in samples removed from the peripheral compartment of the in vitro system at 0, 24, and 36 h (growth control experiments) and at 0, 24, 48, and 96 h (levofloxacin experiments). The frequency of resistant bacteria in the starting cultures at each levofloxacin concentration was calculated by dividing the number of colonies that appeared by the inoculum applied to each plate. For the purposes of this paper, we defined levofloxacin-susceptible bacteria as those recovered only on drug-free agar. Bacteria with low-level levofloxacin resistance were defined as those recovered only on agar containing a levofloxacin concentration of 0.5 µg/ml. Bacteria with high-level resistance were defined as those recovered on agar containing levofloxacin concentrations of
1 µg/ml.
Nucleotide sequence analysis of grlA/B and gyrA/B. Colonies appearing on levofloxacin-containing Mueller-Hinton agar in the MPC studies and the starting cultures during in vitro system experiments were analyzed for nucleotide changes within the grlA/B and gyrA/B genes as described previously (4). Colonies recovered on drug-free agar at the end of the in vitro system experiments were replica plated onto agar containing increasing levofloxacin concentrations, and colonies that appeared on plates with the highest concentrations were selected for sequence analysis. All colonies were examined when three or fewer colonies were present on a plate. Three colonies were randomly chosen for sequence analysis when more than three colonies were present. Genomic DNA was isolated using lysostaphin (Sigma) and the DNeasy tissue kit (QIAGEN, Valencia, Calif.). Primers, PCR conditions, the sequencing procedure, and interpretation of results were identical to those reported previously (4).
Pharmacodynamic modeling. The candidate two-population pharmacodynamic models and their equations were presented in an earlier paper (5). The model equations were simultaneously fit to viable-cell count data obtained from the in vitro system experiments with the two levofloxacin-susceptible parent strains (MRSA 8043 and MRSA 8282) and from additional in vitro system studies that characterized the growth and killing of the most prevalent grlA mutants detected in the parent cultures (MRSA 8043L0-1 and MRSA 8282L0-1). Two-population pharmacodynamic models with equivalent (model 1) and different (model 2) net growth (g) and levofloxacin killing (k) rate constants for levofloxacin-susceptible (S) and -resistant (R) subpopulations were evaluated. The initial equation conditions and the approach for deriving initial parameter estimates were the same as previously described (5).
Pharmacodynamic model equations were fitted to equally weighted log-transformed viable counts as a function of time using a nonlinear regression program (WinNonlin). The goodness of fit was evaluated by determining correlation coefficients (r), performing residual analysis, and assessing the precision of parameter estimates (12). The runs test was used to determine whether model-fitted curves deviated systematically from observed viable-cell count data. The pharmacodynamic model of best fit was determined by the second-order Akaike information criterion (AICc) (3). Mathematical simulations were performed, using mean parameter estimates from the model of best fit, to explore the relationship between the levofloxacin TMSW and the emergence of resistant subpopulations.
| RESULTS |
|---|
|
|
|---|
|
Sequencing of DNA from 30 MRSA 8043 and 30 MRSA 8282 colonies recovered on drug-free agar showed silent mutations or wild-type grlA/B and gyrA/B in all instances. Sequencing of 25 MRSA 8043 and 24 MRSA 8282 colonies recovered on agar containing 0.5 µg/ml of levofloxacin revealed no mutations, silent mutations, or point mutations in grlA and grlB (Table 1). Point mutations in grlA were found in 48% of MRSA 8043 and 83% of MRSA 8282 variants with low-level resistance. The grlB mutations were detected in all three low-level-resistant MRSA 8282 colonies sequenced from the starting population during experiments with the regimen of 125 mg every 24 h. No bacteria with point mutations in gyrA/B were found in the starting cultures of either strain during any experiment.
|
|
Bacterial population dynamics.
Bacteria were in log-phase growth at the start of each experiment. In the absence of levofloxacin, the bacteria continued to grow at an exponential rate until the carrying capacity of the in vitro system was reached (
109 CFU/ml for MRSA 8043 and
1010 CFU/ml for MRSA 8282), and counts remained near this plateau for the remainder of each experiment (Fig. 3).
|
1.5 x 108 to 5.2 x 108 CFU in the 20-ml peripheral compartment for both bacteria. Exposure of cultures to the simulated regimens of 750 and 500 mg every 24 h resulted in a decline in viable counts to below the reliable limit of detection (70 CFU/ml). Small numbers of viable bacteria persisted below this limit for the remainder of the experiments. When the bacteria were exposed to the regimen of 250 mg every 24 h, viable counts decreased to the limit of detection and then increased to slightly above this limit toward the end of the experiments. When the cultures were exposed to a regimen of 125 mg every 24 h, viable counts declined initially but then increased and approached the carrying capacity of the system. The slopes of the killing and regrowth portions of the viable-cell count curves were similar for replicate dosage regimen simulations with a given strain (Fig. 3). Changes in susceptibilities of bacterial populations. No change in levofloxacin MICs (Fig. 1) was observed during growth control experiments with MRSA 8043 and MRSA 8282 or when the cultures were exposed to simulated levofloxacin regimens of 750 mg, 500 mg, and 250 mg every 24 h. In contrast, the MICs increased 2- to 16-fold when the bacteria were exposed to the levofloxacin regimen of 125 mg every 24 h. Reserpine did not decrease the MIC for resistant MRSA 8282 bacteria recovered from the latter experiments but lowered the MIC for resistant MRSA 8043 in one of the replicate experiments by twofold.
Changes in resistant subpopulations and grlA/B and gyrA/B genotypes. During growth control experiments, the number of low-level-resistant variants with or without mutations in grlA increased (Fig. 1 and Table 1). Resistant subpopulations did not emerge when the cultures were exposed to the 750 and 500 mg every 24 h levofloxacin regimens. Bacteria that persisted at the ends of these experiments did not have mutations in grlA/B or gyrA/B. The grlA mutants that were present in some of the starting cultures were still detected at 96 h after exposure to the regimen of 250 mg every 24 h, but no additional mutations were found. In contrast, when the cultures were exposed to the regimen of 125 mg every 24 h, subpopulations with high-level resistance and mutations in grlA, grlB, and gyrA emerged. Different patterns of nucleotide sequence changes in grlA/B and gyrA/B were observed during the latter experiments. In some cases, no mutants were detected in the starting cultures, and resistant bacteria were recovered at 96 h with mutations in grlB and gyrA. In other cases, grlA mutants were present in the starting populations, and resistant bacteria were recovered at 96 h with no additional mutations. In yet other experiments, grlA and grlB mutants were present in the starting populations, and resistant bacteria were recovered at 96 h with additional gyrA mutations. The pattern of nucleotide sequence changes varied with different dosage regimen simulations and between replicate experiments of the same simulation (Table 1).
Pharmacokinetic-pharmacodynamic values and the emergence of resistance. The simulated levofloxacin dosage regimens of 500 mg and 750 mg every 24 h, which provided Cmax/MIC ratios between 10.6 and 18.9 and AUC24/MIC ratios between 113 h and 188 h (Table 2), did not result in the emergence of resistance to levofloxacin. The TMSW for a 24-h period ranged from 19.6 to 45.9% during these experiments. The regimen of 250 mg every 24 h, which produced Cmax/MIC ratios between 5.33 and 6.17 and AUC24/MIC ratios between 56.6 h and 64.1 h, failed to eradicate grlA mutants that were detected in the starting cultures (Table 1) but did not result in a change in levofloxacin MICs at 96 h (Fig. 1). The TMSW was between 55.7 and 62.2% in these experiments. The regimen of 125 mg every 24 h, which produced the lowest Cmax/MIC (2.64 to 3.22) and AUC24/MIC (28.1 h to 32.9 h) ratios resulted in the emergence of subpopulations with mutations in gyrA and high-level resistance (Fig. 1 and Table 1). The TMSW ranged between 46.1 and 57.5% in these experiments.
|
Pharmacodynamic modeling.
The AICc and other model selection statistics indicated that model 2 best described the trends in the viable-cell count data (Table 3). Model 1 was associated with significant bias (P < 0.05, runs test) in model-predicted versus observed viable counts, with the model underpredicting the growth rate of MRSA 8043 and MRSA 8282 and overpredicting the killing rate of MRSA 8043L0-1 and MRSA 8282LC0-1 during the levofloxacin simulations with 1,000 mg every 24 h. There was also a systematic deviation (P < 0.05) in observed viable counts above or below the fitted regrowth curves for replicate experiments with regimens of 250 mg and 125 mg every 24 h. Model 2 better characterized the trends in bacterial killing and growth observed during in vitro system experiments (Fig. 3) and had less bias and more-precise parameter estimates than model 1. The coefficient of variation (CV) of all model parameter estimates, except R0 (the initial size of the R subpopulation), for both strains was
13.8% (Table 4). Estimates of R0 were less precise than the other estimates (CVs, 66.4% for MRSA 8043/MRSA 8043L0-1 and 46.1% for MRSA 8282/MRSA 8282L0-1). There was a strong positive correlation between gS and kS values for both MRSA 8043 (r = 0.990) and MRSA 8282 (r = 0.984). The high interdependence of these parameter estimates was anticipated since the model was fit to viable count data that reflected the net result of growth and killing. The univariate 95% confidence intervals for g and k estimates for S and R subpopulations did not coincide for MRSA 8282/MRSA 8282L0-1, but the confidence intervals of gS and gR estimates overlapped for MRSA 8043/MRSA 8043L0-1. The univariate and planar 95% confidence intervals of the EC50, the concentration that results in 50% of the maximal k, for susceptible and resistant isolates did not coincide.
|
|
7 to 8 log10 lower than those for the S subpopulation. In contrast, the model predicted that R populations grown in the absence of S populations would reach the carrying capacity of the in vitro system. Model predictions for the growth of parent strains and their corresponding grlA mutants were in concordance with total and resistant subpopulation profiles observed during growth control experiments (Fig. 1 and 3). The model predicted that the S subpopulation would be eradicated (i.e., decline to < 1 CFU) during all levofloxacin dosage simulations, with the time required for extinction decreasing as the intensity of levofloxacin exposure increased. For instance, 82 h (MRSA 8043) to 98 h (MRSA 8282) would be required to extinguish S populations at the lowest levofloxacin exposure (125 mg every 24 h) but only 36 h (MRSA 8043) to 39 h (MRSA 8282) would be necessary to eradicate S populations at the highest exposure (750 mg every 24 h).
The model predicted that the levofloxacin regimens of 125 mg and 250 mg every 24 h would result in enrichment of the R subpopulation, with the R populations constituting 0.62% and 0.55% of the total MRSA 8043 population at 24 h for each of the regimens, respectively. The observed percentage of the total MRSA 8043 population resistant to levofloxacin at
0.5 µg/ml for the two regimens was 3.3% and <8.1% (Fig. 1). The model predicted that R would comprise 0.21 and 0.26% of the total MRSA 8282 population by 24 h for the same levofloxacin exposures. The observed percentages were 0.4% and <1.6%. The model predicted that the R population would constitute the entire population by 96 h with both regimens. The observed percentage of bacteria resistant to levofloxacin concentrations of
0.5 µg/ml at 96 h with the regimen of 125 mg every 24 h ranged from 96 to 101%, but the percentage for the regimen of 250 mg every 24 h ranged from <1.6 to 39%.
Mathematical simulations with the pharmacodynamic model identified no clear relationship between TMSW and the enrichment of the R subpopulation. Constant and fluctuating pharmacokinetic profiles that failed to eradicate the R subpopulation and that produced concentrations entirely outside the mutant selection window (TMSW, 0%) or within the window for variable periods (0% < TMSW
100%) were predicted to result in the selection of R populations.
| DISCUSSION |
|---|
|
|
|---|
The starting cultures of all in vitro system experiments consisted mainly of levofloxacin-susceptible bacteria, although minor subpopulations with low-level resistance (MIC, 1 µg/ml) were reproducibly detected at low frequencies (5.7 x 108 to 1.5 x 106). Resistance may have been due to efflux, although our assays with reserpine suggested that this played a minor role. Genotyping of the levofloxacin-resistant variants in the starting cultures revealed a variety of mutations in grlA/B. The S80Y, S80F, E84K, and A116P grlA mutations observed have been associated with phenotypic levofloxacin resistance in S. aureus (15, 22, 35), but the P144S grlA and the E422D grlB mutations do not appear to contribute to resistance (16, 26, 33). Other resistant subpopulations did not have mutations in the sequenced regions of grlA/B or gyrA/B. Resistance may have been due to mutations in other regions of these genes. The presence of minor subpopulations of low-level-resistant variants in the starting cultures and the variety of mutations among these bacteria are consistent with the concept that mutations arise randomly from wild-type cells during nonselective growth, and these findings are similar to our experience with ciprofloxacin (4).
Different patterns of killing and regrowth were observed in the in vitro system during the levofloxacin simulations. The initial rate of bacterial killing approached a maximum, and the rate of regrowth decreased with increasing Cavg ss. The emergence of levofloxacin resistance also varied in the different pharmacokinetic environments. The greatest increases in MICs were observed with the simulated experimental dosing regimen of 125 mg every 24 h, providing Cavg ss that were just above the MIC but within the mutant selection windows of both strains. In these experiments, low-level-resistant variants present in the starting cultures became a greater proportion of the population within 24 h. During the remainder of the experiment, their numbers increased and bacteria with a variety of mutations in grlA/B or gyrA and high-level resistance appeared. With the simulated clinical regimen of 250 mg given every 24 h, which also produced a Cavg ss within the MSW but closer to the MPC, bacterial numbers initially fell to the limit of detection but then slowly increased. Bacteria recovered at 96 h did not have an increased MIC measured by broth microdilution, although a slight increase (less than twofold) was detected when the MIC was measured by the Etest (AB Biodisk, Piscataway, N.J.). No additional mutations in grlA/B or gyrA/B were detected. Resistance did not develop when the Cavg ss exceeded the MPC. With these clinical regimens (500 and 750 mg every 24 h), viable counts declined to the limit of detection within 24 h, and the minor subpopulations of grlA mutants that were present in the starting populations appeared to be eliminated. These findings are consistent with our ciprofloxacin experience (4) and the suggestion of others that the key to successful antimicrobial therapy is eradication of minor low-level-resistant subpopulations present in the inoculum (1, 2).
A variety of different grlA/B and gyrA/B genotypes that emerged with exposure to levofloxacin were seen in resistant bacteria, and there was no apparent pattern to the presence or type of mutations that appeared. Genotypic differences occurred in several cases, even though the bacteria had been exposed to the same simulated levofloxacin regimen, and bacteria with the same MIC often had different mutations. As with ciprofloxacin (4), these data suggest that bacteria follow heterogeneous pathways of evolution and that the sequence of resistance determinants is not predetermined.
Our work supports the concept that resistant subpopulations of S. aureus are selectively enriched when levofloxacin concentrations fall within the MSW, although we could not find a relationship between the TMSW and the degree of resistance as reported by others (11). In our in vitro system, high-level-resistant variants did not appear with the regimen of 250 mg every 24 h, producing a TMSW that was slightly higher than that of the regimen of 125 mg every 24 h, with which resistant variants appeared. It is also of note that the TMSW for the regimens of 125 mg every 24 h and 500 mg every 24 h were very similar, and yet resistance did not appear with the latter regimen. Drug concentrations outside the MSW were below the MIC in the former regimen but above the MPC in the latter regimen. These findings are consistent with those of our previous work with ciprofloxacin (4) and again suggest that pharmacokinetic parameters other than TMSW are important for the selection of resistant variants.
Our data showed that resistance did not occur with the Cmax/MIC ratios above 5.33 and AUC24/MIC ratios above 56.6 h afforded by the dosing simulations of 250, 500, and 750 mg every 24 h but occurred with Cmax/MICs up to 3.22 and AUC24/MICs up to 32.9 h provided by the regimen of 125 mg every 24 h. These findings are consistent with there being a relationship between antimicrobial exposure and the selection of resistant variants. Additional experiments using a wider range of concentrations and incremental changes in Cmax/MIC or AUC24/MIC would be useful to fully evaluate these pharmacodynamic concepts and determine optimal values for levofloxacin in the in vitro system. The Cmax/MIC or AUC24/MICs obtained in the in vitro system represent free (unbound) drug exposures. We acknowledge that our experimental design may have overestimated the killing activity for a given levofloxacin regimen because the pharmacokinetic profiles simulated may not reflect the free-drug exposures in normal volunteers or infected patients if protein binding (
30%) is considered.
Our findings are consistent with an earlier study where no resistance occurred when S. aureus was exposed to levofloxacin daily doses of 400 to 800 mg in an in vitro glass infection model (21). Our results are similar to those of another study where a levofloxacin-susceptible (MIC, 0.5 µg/ml) S. aureus strain was exposed to a simulated levofloxacin regimen of 500 mg every 24 h in a hollow-fiber in vitro system (25). Viable counts for this susceptible strain declined to below the lower limit of detection when cells were exposed to this dosage regimen simulation as in our experiments. In the same study, when strains with reduced susceptibility to levofloxacin (MICs of 2 to 4 µg/ml) were exposed to the identical regimen, viable counts decreased but then increased. These results are similar to ours with the levofloxacin simulation of 125 mg every 12 h, perhaps because the fluoroquinolone exposures relative to the MIC were similar. A different relationship between levofloxacin exposure and increases in MIC was reported by others (11). Increases in MICs of up to 2.5-fold occurred with simulated regimens producing Cmax/MIC and AUC24/MIC ratios closely approximating those of our 125-mg, 250-mg, and 500-mg levofloxacin simulations. The emergence of resistance within 72 h with the 250-mg and the 500-mg simulations in this study may have been due to the larger inocula used (10- to 40-fold higher than ours). This would increase the likelihood that more spontaneously occurring low-level-resistant variants would be present at the outset of the experiments. Our pharmacodynamic model predicted that if R0 for MRSA 8043 and MRSA 8282 was increased by 40-fold, resistant subpopulations would emerge within 72 h with a levofloxacin regimen of 500 mg every 24 h.
The two-population pharmacodynamic models of this study were developed using an approach similar to that used to model the evolution of ciprofloxacin resistance (5). Other investigators have published similar analyses of multipopulation responses to antimicrobial pressure both in animals (17) and in an in vitro system (18). Low-level-resistant grlA mutants of each strain were incorporated into our models as the resistant subpopulation because they represented the most prevalent resistant phenotype and genotype detected in the starting cultures. This approach differed from that of other investigations where model parameter estimates for the resistant subpopulation were derived solely from the regrowth phase of the total viable count-time curves and not from additional independent experiments using resistant bacteria (23, 28, 36). Although both of our candidate pharmacodynamic models described the patterns of bacterial growth and killing observed in the in vitro system, the simpler model (model 1), with equivalent growth and levofloxacin killing rate constants, had more biased parameter estimates. It underpredicted the growth of the susceptible population and overpredicted the killing rate of the resistant subpopulation by levofloxacin. The alternate model, with unique growth and killing rate constants for susceptible and resistant subpopulations (model 2), resulted in parameter estimates with less bias and greater precision. We did not consider more-complex pharmacodynamic models (e.g., ones that allowed for adaptive resistance during levofloxacin exposure or that had more than two subpopulations with different susceptibility phenotypes or genotypes) because our previous experience with ciprofloxacin suggested that it would be difficult to obtain precise estimates for the additional parameters from the viable count data (5).
The net growth and levofloxacin killing rate constants of the resistant subpopulation were lower, and the EC50 estimates were higher than the respective parameters of the susceptible subpopulation. This was similar to our findings with ciprofloxacin (5). The planar confidence intervals for EC50s of the susceptible and resistant subpopulations did not overlap, unlike those for g and k, suggesting that the lower rate of killing of the resistant bacteria at comparable levofloxacin exposures was mainly due to differences in EC50s. Estimates of R0s were less precise than other parameter estimates, as was the case when we applied the model to ciprofloxacin pharmacodynamics. The R0 estimates were consistent with the observed resistance frequencies in our starting cultures and with measured resistance frequencies reported by others (10, 14, 24, 34).
The parameter estimates developed with the model using grlA mutants may not be applicable to bacteria with different resistance mechanisms (e.g., efflux or small-colony variants). As with ciprofloxacin (5), however, the model correctly characterized the emergence of subpopulations with low-level resistance even though grlA mutants were not detected in all of the starting cultures. This suggests that knowing the susceptibility phenotypes of the resistant subpopulations may be sufficient for designing regimens to prevent the emergence of resistance.
Predictions of the two-population model were often correct, but limitations were evident. For instance, the model correctly predicted that resistance would not occur with the regimens of 750 mg and 500 mg every 24 h, but it failed to predict the persistence of low numbers of susceptible bacteria. The model correctly predicted the elimination of susceptible bacteria and the enrichment of low-level-resistant variants during exposure to 125 mg every 24 h exposures but did not explain the origin or appearance of additional subpopulations with high-level resistance and mutations in gyrA. The model also predicted that low-level-resistant variants would be enriched with the regimen of 250 mg every 24 h. Bacterial numbers did increase toward the end of the 96-h experiment (Fig. 3), but the number of grlA mutants and other resistant variants present in the starting cultures did not increase above the reliable limit of detection (Fig. 1). It is possible that the bacteria with slightly increased MICs (detected by the Etest but not broth microdilution) detected at the end of the experiment may have become the predominant subpopulation if the experiments had been continued beyond 96 h.
The ability of the model to describe the population dynamics of bacteria exposed to both levofloxacin and ciprofloxacin suggests that the model is robust and may have general applicability to describing and predicting the relationship between fluoroquinolone exposure and the evolution of resistance in S. aureus. The model may even be applicable to other antimicrobials, such as rifampin or linezolid, where resistance occurs by point mutations (27, 31, 32).
Resistance did not emerge with simulated levofloxacin pharmacokinetic profiles that provided lower Cmax/MIC and AUC24/MIC ratios than those of the ciprofloxacin profiles that resulted in the selection of resistant variants (4). This finding suggests that the pharmacokinetic-pharmacodynamic ratios required to prevent the emergence of resistance in S. aureus vary for different fluoroquinolones. Our robust pharmacodynamic model provides a theoretical framework for understanding the bacterial strain and antimicrobial characteristics that influence the emergence of resistance in S. aureus populations. Mathematical simulations with the model are warranted to examine strain and fluoroquinolone factors that influence the pharmacokinetic-pharmacodynamic indices that are required to prevent the evolution of resistance in S. aureus.
| ACKNOWLEDGMENTS |
|---|
We thank Glenn A. Kaatz (Wayne State University, School of Medicine) for kindly providing the S. aureus control strains, SA1199 and SA1199B, used in the fluoroquinolone efflux screening studies.
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Clin. Vaccine Immunol. | Clin. Microbiol. Rev. |
|---|---|
| J. Clin. Microbiol. | ALL ASM JOURNALS |