Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, November 2008, p. 3987-3993, Vol. 52, No. 11
0066-4804/08/$08.00+0 doi:10.1128/AAC.01468-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Department of Clinical Sciences and Administration, College of Pharmacy,1 Department of Chemical and Biomolecular Engineering, College of Engineering, University of Houston, Houston, Texas,3 Pharmacy Department, Singapore General Hospital, Republic of Singapore2
Received 13 November 2007/ Returned for modification 9 March 2008/ Accepted 9 July 2008
|
|
|---|
|
|
|---|
When evaluating various dosing regimens of antimicrobial agents, the total daily dose, the dosing frequency, the dose given at each dosing interval, the length of (intravenous) administration, and the duration of therapy may have a significant influence on the killing activity and propensity to suppress resistance emergence, depending on the pharmacodynamic properties of the agents and clinically achievable concentrations (associated with acceptable toxicity). The numerous combinations of these variables involved in designing dosing regimens are prohibitive for comprehensive evaluation of all the different scenarios. For example, to evaluate six daily doses (e.g., 0.5, 1, 2, 4, 6, and 8 g), four dosing frequencies (e.g., every 6 hours [q6h], q8h, q12h, and q24h), five intravenous dosing administrations (e.g., intermittent infusion of 0.5, 1, 2, and 4 h and continuous infusion over 24 h), and three durations of treatment (e.g., 5, 10, and 14 days) would result in 360 (6 x 4 x 5 x 3) regimens to be investigated. In view of the labor-intensiveness of each investigation, a few dosing regimens are often empirically chosen to be studied, and the potential of new agents may not be thoroughly realized. The results could be dramatic. The development of daptomycin was unsuccessful in the 1970s due to an empirical dosing selection (8-h dosing interval). However, with an improved understanding of pharmacokinetics/pharmacodynamics as the foundation for a dosing strategy, the same agent was redeveloped and approved by the FDA for clinical use in 2003 by switching to a once-daily and weight-based dosing regimen. Therefore, a robust method to guide the selection of the most effective dosing regimen(s) would be valuable.
We previously developed a mathematical model to capture the dynamic relationship between a heterogeneous microbial population and drug concentrations (12, 15). The model was further refined to efficiently predict the microbial response to multiple antimicrobial agent dosing regimens (16). This modeling approach could be used as a decision support tool for dosing regimen design, and it may be used at different stages of drug development. One such application could be at the interface between late discovery and early clinical development (i.e., lead compound optimization). The model could provide a prediction of the new agent's potential clinical utility based on simple time-kill studies and allometrically scaled pharmacokinetic data. Another possibility is the comparison of several related drug candidates with conflicting in vitro potencies and pharmacokinetics (e.g., a more potent agent with a shorter elimination half-life versus a less potent agent with a longer elimination half-life). However, before such mathematical modeling can be used routinely to guide highly targeted investigation of dosing regimens in preclinical studies and clinical trials, its utility should be investigated in more than one drug-pathogen combination. Therefore, the objective of this study was to examine the flexibility of our mathematical-modeling approach using additional drug-pathogen combinations. Two wild-type gram-negative bacteria deemed to be "particularly problematic" and in need of new drug treatment were chosen for our investigations (11).
(This study was presented in part at the 16th European Conference of Clinical Microbiology and Infectious Diseases [ECCMID], Munich, Germany, 31 March to 3 April 2007, and the 17th ECCMID, Barcelona, Spain, 19 to 22 April 2008.)
|
|
|---|
Microorganisms. Pseudomonas aeruginosa ATCC 27853 and Acinetobacter baumannii ATCC BAA 747 (American Type Culture Collection, Rockville, MD) were used in the study. The bacteria were stored at –70°C in Protect (Key Scientific Products, Round Rock, TX) storage vials. Fresh isolates were subcultured twice on 5% blood agar plates (Hardy Diagnostics, Santa Maria, CA) for 24 h at 35°C prior to each experiment.
Susceptibility studies.
MICs/minimum bactericidal concentrations (MBCs) were determined in Ca-MHB using a modified broth macrodilution method as described by the CLSI (1). The final concentration of bacteria in each broth macrodilution tube was approximately 5 x 105 CFU/ml of Ca-MHB. Serial twofold dilutions of the aminoglycosides were used. The MIC was defined as the lowest concentration of drug that resulted in no visible growth after 24 h of incubation at 35°C in ambient air. Samples (50 µl) from clear tubes and the cloudy tube with the highest drug concentration were plated on Mueller-Hinton agar (MHA) plates (Hardy Diagnostics, Santa Maria, CA). The MBC was defined as the lowest concentration of drug that resulted in a
99.9% kill of the initial inoculum. The drug carryover effect was assessed by visual inspection of the distribution of colonies on medium plates. The studies were conducted in duplicate and repeated at least once on a separate day.
Time-kill studies and mathematical modeling. Time-kill study data of P. aeruginosa over 24 h have been reported previously (12); a clinically relevant concentration range of gentamicin (from 0 to 32 mg/liter) was used. Similarly, the experiment was repeated with A. baumannii using a clinically relevant concentration range of amikacin (from 0 to 128 mg/liter). The mathematical structure of the growth dynamics model is shown in Fig. 1. Briefly, the rate of change of bacteria over time was expressed as the difference between the intrinsic bacterial growth rate and the (sigmoidal) kill rate provided by the antimicrobial agent. A decline in the kill rate over time and regrowth were attributed to adaptation, which was explicitly modeled as an increase in the concentration necessary to achieve a 50% maximal kill rate (C50k), using a saturable function of antimicrobial agent selective pressure (both the aminoglycoside concentration [C] and time). The time-kill study data were used as inputs to derive the best-fit model parameter estimates, as previously described (12, 15). The modeling estimation process involved two steps. For each bacterium, the intrinsic bacterial growth rate (Kg) and maximal bacterial population size (to account for contact inhibition) were first determined from placebo (control) experiments, using the ADAPT II program (2). Using these parameter estimates, the parameter values in the kill function were subsequently determined using data from all active treatment experiments simultaneously. The baseline inoculum (107 CFU/ml) was deemed to be dense enough to constitute a heterogeneous bacterial population.
![]() View larger version (35K): [in a new window] |
FIG. 1. Bacterial growth dynamics model and various model parameters (adapted with permission from reference 15).
|
Experimental validation. Regardless of the prediction method used, the computer simulations were compared to experimental data from an in vitro hollow-fiber infection model with similar antimicrobial agent exposures. The experimental setup has been described elsewhere (16). A human-like elimination half-life (approximately 2.5 h) for both gentamicin and amikacin was simulated in the infection models. Serial samples were obtained from the infection models over time to ascertain the simulated pharmacokinetic exposures. The aminoglycoside concentrations in these samples were assayed using validated methods, as detailed below. A one-compartment linear model was fitted to the observed time-concentration profiles using the ADAPT II program (2).
In addition, serial samples were obtained at baseline, 4, and 8 h and daily (predose) in duplicate from each hollow-fiber system for quantitative culture to define the effects of various drug exposures on the bacterial population. Prior to culturing the bacteria quantitatively, the bacterial samples were centrifuged at 10,000 x g for 15 min and reconstituted with sterile normal saline in order to minimize the drug carryover effect. Total bacterial populations were quantified by spiral plating (Spiral Biotech, Bethesda, MD) 10x serial dilutions of the samples (50 µl) onto drug-free MHA plates. Subpopulations with reduced susceptibility (resistant) were quantified by culturing them on cation-adjusted MHA plates supplemented with the exposed agent (gentamicin or amikacin) at a concentration of 3x MIC. Since susceptibility testing is performed in twofold dilutions and one tube (2x concentration) difference is commonly accepted as reasonable interday variation, quantitative cultures on drug-supplemented medium plates (at 3x MIC) would allow reliable detection of bacterial subpopulations with reduced susceptibility. The medium plates were incubated at 35°C for up to 24 (total population) and 72 (subpopulations with reduced susceptibility) h, and the bacterial density of each sample was enumerated visually. The theoretical lower limit of detection was 400 CFU/ml. Determination of the susceptibilities of the resistant isolates (recovered from the drug-supplemented medium plates at the end of the experiments) to the exposed agent was repeated to confirm the emergence of resistance.
Pharmacokinetic profiles investigated. All dosing regimens investigated were guided by computer model predictions. In view of the pharmacodynamic property of gentamicin, preliminary simulations revealed that suppression of resistance would be unlikely using any clinically achievable exposures. Therefore, as a proof of concept, several supraphysiologic dosing regimens were investigated, in addition to two clinically relevant dosing regimens (using conventional nomenclature, a Cmax/MIC of 4 q8h and a Cmax/MIC of 12 q24h). On the other hand, four clinically achievable dosing regimens of amikacin were examined, corresponding to a Cmax /MIC of 5 q12h, a Cmax/MIC of 6 q8h, a Cmax/MIC of 13 q12h, and a Cmax /MIC of 20 q24h.
Bioassays.
Gentamicin concentrations were determined by a microbioassay utilizing Klebsiella pneumoniae ATCC 13883 as the reference organism. The bacteria were incorporated into 30 ml of molten cation-adjusted MHA (at 50°C) to achieve a final concentration of approximately 1 x 105 CFU/ml. The agar was allowed to solidify in 150-mm medium plates. A size 3 cork bore was used to create nine wells in the agar per plate. Standards and samples were tested in duplicate with 40 µl of the appropriate solution in each well. The gentamicin standard solutions ranged from 1 to 32 mg/liter in Ca-MHB. The medium plates were incubated at 35°C for 24 h, and the zones of inhibition were measured. The assay was linear (correlation coefficient,
0.99), using the zone diameter versus the log of the standard drug concentration. The intraday and interday coefficients of variation for all standards were <4% and <6%, respectively. Similarly, amikacin concentrations were determined by a microbioassay utilizing Escherichia coli ATCC 25922 as the reference organism. The assay was linear (correlation coefficient,
0.98), using the zone diameter versus the log of the amikacin standard concentrations from 4 to 256 mg/liter. The intraday and interday coefficients of variation for all standards were <7% and <12%, respectively.
|
|
|---|
Time-kill studies. Data from the time-kill studies and model fits to the data are shown in Fig. 2 and 3, respectively. The estimates of the best-fit model parameters are shown in Table 1. Taken as a whole, the observations in bacterial burdens over time (under constant antimicrobial agent concentrations) were reasonably described by the model.
![]() View larger version (16K): [in a new window] |
FIG. 2. Time-kill studies of gentamicin against P. aeruginosa ATCC 27853 (A) and of amikacin against A. baumannii ATCC BAA 747 (B). The data are shown as means ± standard deviations. Complete bacterial eradication was observed with amikacin concentrations of >16 mg/liter after 2 h of drug exposure.
|
![]() View larger version (9K): [in a new window] |
FIG. 3. Model fits to the experimental data in time kill studies. Shown are gentamicin against P. aeruginosa ATCC 27853 (A) and amikacin against A. baumannii ATCC BAA 747 (B).
|
|
View this table: [in a new window] |
TABLE 1. Susceptibilities of isolates and final estimates of best-fit model parameters in time kill studiesa
|
![]() View larger version (15K): [in a new window] |
FIG. 4. Comparison of computer-simulated and experimental bacterial (P. aeruginosa) responses to various gentamicin exposures. The dosing frequencies are in parentheses. The data are shown as means ± standard deviations. The horizontal dotted line [in Cmax/MIC = 30 (q12h)] depicts the reliable lower limit of detection.
|
![]() View larger version (37K): [in a new window] |
FIG. 5. Likelihood of emergence of bacterial (A. baumannii) resistance to various dosing regimens of amikacin as predicted by a response surface analysis. Two intersecting planes are shown: a translucent mesh surface (representing different dosing regimens) and an opaque surface (where D/Kg = 1). The three-dimensional mesh surface is made up of a collection of data points; each datum point is characterized by a value on the x, y, and z axes corresponding to the dose (x), dosing interval (y), and D/Kg (z). The average kill rate of various dosing regimens against the most resistant bacterial subpopulation was quantified by D, an index of dosing intensity (D Kk [the maximal kill rate], and D is dependent on the Cmax and dosing frequency [an intrinsic property of a dosing regimen]). For a dosing regimen to suppress resistance amplification, it is imperative that the average kill rate (D) be more than the intrinsic growth rate (Kg) of the bacterial population. Both D and Kg are expressed in h–1, so D/Kg is a dimensionless ratio. To identify promising dosing regimens (combinations of dose and dosing interval) to prevent resistance development, the corresponding values of D/Kg should be >1 (the region where the translucent mesh surface is above the opaque plane), as indicated by the arrow (e.g., 2,000 mg [Cmax/MIC = 20] given every 24 h or 1,500 mg [Cmax/MIC = 15] every 12 h). Note: Cmax/MIC = dose/(volume of distribution x MIC) = dose/(70 x 0.35 x 4).
|
![]() View larger version (7K): [in a new window] |
FIG. 6. Validation of microbial responses to various amikacin dosing regimens. The dosing frequencies are in parentheses. The data are shown as means ± standard deviations. The horizontal dotted lines depict the reliable lower limits of detection.
|
|
|
|---|
Rational dosing regimen design involves multiple variables. Since comprehensive evaluation of all combinations is prohibitive in view of the labor-intensiveness of each investigation, the initial choice of the dosing regimens to be tested preclinically is often empirical (mostly trial and error). Given that studies with these infection models are costly and time-consuming to perform, poorly guided exploration studies to examine the potentials of various dosing regimens (e.g., dose escalation or dose fractionation studies) may not be very efficient and cost-effective. Of interest, conventional pharmacodynamic indices predicting outcomes may also be subjected to the range of drug concentrations examined (see the supplemental material labeled "Appendix 2").
In contrast, our proposed modeling approach does not require the use of surrogate pharmacodynamic indices to make useful predictions of microbial response to antimicrobial exposures (the Cmax/MIC was empirically chosen to represent different dosing regimens in Fig. 4 and 6 due to convention; other indices could also be used). As such, it offers a method to improve dosing regimen selection, which could suppress the development of resistance during therapy. With this method, standard time-kill study data over 24 h are used as model inputs. The utility of a large number of dosing regimens can be effectively screened in a comprehensive fashion, but only promising ones would be investigated in preclinical studies and clinical trials. In addition, because the dosing regimens investigated are designed to prevent resistance emergence, the clinical-utility life span of new agents would also likely be prolonged.
In all cases, bacterial regrowth after an initial decline in the bacterial burden was attributed to selective amplification of a resistant subpopulation(s). We were able to confirm the presence of phenotypic resistance by repeated susceptibility testing. However, we did not ascertain the specific mechanism of aminoglycoside resistance in this study, as it was not the primary focus of the study. Since the acquisition of additional genetic material for A. baumannii was highly unlikely in our in vitro infection model, initial screening attempts did not suggest a common mechanism of resistance. Nonetheless, we believe that the experimental data shown are adequate to demonstrate the utility of our modeling approach.
As with previous studies, the proposed modeling approach was expected to be a general evaluative tool. The utility of the model was illustrated in this study by experimental data for P. aeruginosa and A. baumannii. However, the proposed model-based approach is not confined to a specific antimicrobial agent-pathogen combination. An attractive feature of our model is that it does not rely on detailed knowledge of the mechanism of resistance. The different pharmacologic and microbiologic characteristics of various drug-pathogen combinations can be represented by different model parameter values (e.g., the growth rate constant of the pathogen and the dispersion of the maximal kill rate constant of the antimicrobial agent under investigation) derived from actual (short-term) experiments. While various antimicrobial agents may have different mechanisms of action or killing profiles (e.g., bactericidal or bacteriostatic), and various pathogens may have different biological characteristics, the same mathematical model structure can still be used (e.g., different extents of susceptibility reduction due to different mechanisms of resistance can be reflected in the value of
, and the rate of resistance selection can be reflected in the value of
). As the MIC increases, the C50k for the surviving population also increases in a nonlinear fashion until "maximal adaptation" has been achieved. In the past, the same (a two-subpopulation) mathematical-model structure has been used to examine the effects of quinolone exposures on P. aeruginosa (14), Mycobacterium tuberculosis (4), and Staphylococcus aureus (16). These studies suggest that an appropriate mathematical-model structure capturing essential features can be developed, which would be flexible enough to characterize the dynamic interaction of more than one antimicrobial agent-pathogen combination. Consequently, the proposed model would likely be extrapolated to other antimicrobial agents (e.g., antibacterials, antifungals, and antivirals) with different mechanisms of action, as well as to other pathogens (e.g., human immunodeficiency virus, tuberculosis, anthrax, and avian influenza) with different biological characteristics.
In conclusion, using limited data from time-kill studies over 24 h, our mathematical model was reasonable in qualitatively predicting extended microbial responses to various concentration-time profiles of both gentamicin and amikacin over 3 days. This approach appears promising as a decision support tool to guide highly targeted investigation of dosing regimens in preclinical investigations and clinical studies. The in vivo relevance and sensitivity of the mathematical-model predictions are currently under investigation.
We thank George H. Miller for advice on the phenotypic screening and the mechanism of aminoglycoside resistance.
Published ahead of print on 25 August 2008. ![]()
Supplemental material for this article may be found at http://aac.asm.org/. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»