Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, April 2007, p. 1246-1252, Vol. 51, No. 4
0066-4804/07/$08.00+0 doi:10.1128/AAC.00685-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Cognigen Corporation, Buffalo, New York,1 UCSD Pediatric ACTG Pharmacology Laboratory, San Diego, California,2 Division of Infectious Diseases, Children's Hospital, San Diego, California,3 Division of Pediatric Pharmacology and Critical Care, Rainbow Babies and Children's Hospital, Cleveland, Ohio,4 Department of Pediatrics and Pharmacology, School of Medicine, Case Western Reserve University, Cleveland, Ohio,5 Departments of Pediatrics and Pharmacology, University of MissouriKansas City, and Pediatric Pharmacology Research Unit, Children's Mercy Hospitals and Clinics, Kansas City, Missouri,6 Department of Pediatrics, School of Medicine, Case Western Reserve University, Cleveland, Ohio,7 Pediatric Infectious Disease, University of Arkansas for Medical Sciences, and Arkansas Children's Hospital, Little Rock, Arkansas,8 Department of Clinical Pharmacology/Experimental Medicine, Bristol-Myers Squibb Pharmaceutical Research Institute, Princeton, New Jersey9
Received 2 June 2006/ Returned for modification 31 October 2006/ Accepted 28 December 2006
|
|
|---|
|
|
|---|
In adults, orally administered gatifloxacin is well absorbed in a dose-independent fashion, with an absolute bioavailability of approximately 96% (23). The rate and extent of absorption are not significantly affected by food (22). Gatifloxacin undergoes limited biotransformation and is eliminated primarily unchanged in the urine (70 to 80%), with an elimination half-life of approximately 8 to 10 h (14, 28). Gatifloxacin exhibits linear pharmacokinetics (PK), with clearance being independent of dose (14). The mean area under the concentration-time curve from 0 to 24 h (AUC0-24) ± standard deviation for adults given multiple daily doses of 400 mg is 34.4 ± 4.1 mg·h/liter (28).
The pharmacokinetics of gatifloxacin in pediatric patients were investigated in a single-dose safety/pharmacokinetic study utilizing intensive blood sampling and noncompartmental pharmacokinetic analysis methods (9). While that study did examine the potential impact of gatifloxacin disposition on pediatric-dose selection, it did not fully consider all of the potentially relevant factors that can contribute to variability in pharmacokinetics, which are afforded by population-based techniques. In general, population pharmacokinetic models have the additional advantage of providing precise estimates of the variance-covariance between pharmacokinetic parameters, useful for subsequent simulation experiments, and are applicable to sparse sampling strategies often employed in phase II/III pediatric studies (5, 31, 32). Thus, the data collected in the above-described study were used to develop a population pharmacokinetic model describing the patient factors responsible for the intersubject variability in gatifloxacin pharmacokinetics. This model is suitable for use, in combination with sparse-concentration sampling strategies often employed in pediatric studies, to estimate exposure in future studies, potentially leading to the identification of PK and pharmacodynamic (PD) targets for efficacy against pediatric respiratory tract infections. Furthermore, this model should also be suitable for the conduct of robust simulation experiments critical to defining appropriate dosing regimens as well as to evaluating the comparative utility of anti-infective agents (2, 7).
|
|
|---|
2.0 years, >2.0 to
6.0 years, >6.0 to
12.0 years, and >12.0 to
16.0 years. Within each age stratum, patients were administered one dose of 5, 10, or 15 mg/kg of body weight of gatifloxacin as an oral suspension, up to a maximum dose of 600 mg. All patients meeting the eligibility criteria as defined in the protocol (82 of the 84 patients) were included in the analysis. Eleven blood samples were collected from each patient at the following times after dosing: 0.5, 0.75, 1, 1.5, 2, 3, 4, 6, 8, 12, and 24 h. Additional information was collected regarding patient demographics and covariate factors. This information included weight, height, age, gender, race, and serum creatinine concentration. Creatinine clearance was calculated from the measured serum creatinine concentration using the Schwartz formula (30). Body surface area (BSA) was calculated using the method of Gehan and George (12). Based on the distribution of patient race, each patient was categorized as white, black, or other.
Drug analysis.
Plasma specimens were analyzed for gatifloxacin using a previously described, validated high-performance liquid chromatography fluorescence method (23). The lower limit of quantitation of gatifloxacin in plasma was 0.01 µg/ml; standard curves were linear (r2
0.983) over the concentration range of 0.01 to 5 µg/ml. The between-run and within-run levels of precision for quality control samples were no greater than coefficients of variation of 9.5 and 10.1% respectively, with deviations from nominal concentrations being within ±5.6%.
Population pharmacokinetic model development. The gatifloxacin concentration-time data were fit to a one-compartment model using the NONMEM computer program (version 5) (6). Due to the large number of samples per patient, the first-order conditional estimation method with interaction was used (4). When a one-compartment model with first-order absorption and elimination was evaluated (subroutines ADVAN2 and TRANS2), the apparent oral clearance (CL/F, where F is oral bioavailability), apparent volume of distribution (V/F), and first-order absorption rate constant (ka) parameters were estimated. Interindividual variability in CL/F, V/F, and ka values was modeled using an exponential-error model. The two error models evaluated to estimate residual variability were a proportional-error model and a combined model of additive errors plus proportional errors.
Patient covariates were added to the basic model using a forward-selection procedure. The covariates evaluated included age, weight, body surface area, serum creatinine, creatinine clearance, gender, and race. Initially, Bayesian estimates of clearance generated from the basic model were plotted against each covariate to determine potential functional forms for the relationship (i.e., linear, piece-wise linear, and nonlinear). Each covariate of interest was then added to the model for clearance (or elimination rate) and volume of distribution separately, and the results of these runs were compared. All covariate effects whose addition resulted in a significant reduction in the value of the objective function (a decrease of at least 3.84 [
= 0.05, 1 degree of freedom]) were included in the full multivariable model. The resulting multivariable model was used to examine the appropriateness of the interindividual-variability and residual-error models. A recursive backward-elimination procedure was then performed to further refine the model and eliminate correlated predictors. During each step of backward elimination, the covariate whose removal resulted in the smallest nonsignificant increase in the objective function (a change of no more than 10.83 [
= 0.001, 1 degree of freedom]) was removed from the model. This process was repeated until no further covariate effects could be removed from the model without a statistically significant difference in the minimum value of the objective function. The final model was then assessed for the appropriateness of the interindividual-variability and residual-error models.
Statistical analysis.
Statistical significance was assessed by the change in the log likelihood value obtained for various models. For each analysis, NONMEM computed the minimum value of the objective function, a statistic that is proportional to minus twice the log likelihood of the data. In the case of hierarchical models, the change in the minimum value of the objective function produced by the inclusion of a parameter is asymptotically distributed as
2, with the number of degrees of freedom equal to the number of parameters added to or deleted from the model. When alternative models could not be cast as hierarchical, the change in the objective function was used only as a qualitative measure of statistical significance.
The goodness-of-fit of each NONMEM analysis was assessed by the examination of the following: scatter plots of predicted concentrations versus measured concentrations and predicted concentrations versus weighted residuals, the precision of the parameter estimates as measured by the percent standard error of the mean (standard error/parameter estimate x 100%), and changes in the estimates of interindividual variability and residual variability for the specified model.
Model validation.
Due to the relatively small sample size, a cross-validation approach was taken. In this approach, the data were split into 10 groups of similar demographic characteristics in order to provide 10 distinct validation data sets ("test data sets") of eight patients each. For the first test set, the final population model was fit to data from the remaining 74 patients (the "training set") to obtain mean population pharmacokinetic parameter estimates. The model, with parameter values fixed to those determined from the training data set, was then applied to the test data set in order to predict the measured concentrations. Because the total sample size was 82 subjects, 2 subjects were never included in a test set but were included in the "training data sets." This process was then repeated for the remaining nine test data sets. The bias and precision of the predicted concentrations were assessed by means of the percent prediction error (PEP) and the absolute PEP (APEP) calculated using equations 1 and 2 below:
![]() | (1) |
![]() | (2) |
Simulation of exposure using alternative dosing regimens. The pediatric dosing regimen for gatifloxacin oral suspension previously under consideration was as follows: 10 mg/kg/dose given once daily up to a maximum dose of 400 mg once daily. As the maximum dose allowed in this study was 600 mg once daily, it was necessary to conduct simulations in order to evaluate the suitability of the regimen under consideration. These simulations were conducted by creating a data set of 1,134 hypothetical pediatric subjects with age-appropriate weight and body surface area characteristics. The age appropriateness of the patient characteristics was achieved by first simulating age using a uniform distribution and then simulating the patient's height and weight using the CDC reference charts.
Using the fixed- and random-effect models from the final population PK model, two AUC0-24 values were generated for each simulated patient, one using a 10-mg/kg regimen with no upper dose limit and one using the proposed regimen of 10 mg/kg up to 400 mg. An estimate for the free-drug AUC-to-MIC ratio (fraction unbound AUC [fuAUC]/MIC90) was obtained by dividing the free-drug AUC (0.8·AUC0-24) estimated for the proposed regimen (10 mg/kg up to 400 mg) by the MIC90 of gatifloxacin to Streptococcus pneumoniae (0.5 µg/ml) (19).
|
|
|---|
2.0, 18 being between 2.0 and 6.0, 22 being between 6.1 and 12.0, and 20 being >12.0 years of age. Scatter plots of gatifloxacin plasma concentrations versus time since the last dose are illustrated in Fig. 1, stratified by dose group. |
View this table: [in a new window] |
TABLE 1. Demographic summary of the patients included in the population pharmacokinetic model development
|
![]() View larger version (13K): [in a new window] |
FIG. 1. Scatter plots (by dose) of observed gatifloxacin concentrations versus time since the last dose.
|
Covariate and interindividual-variability models. Table 2 details the results of covariate model building. After a screening of various size functions on CL/F and V/F, a linear relationship between body surface area and CL/F combined with a linear relationship between weight and V/F best described the data. Forward-selection and backward-elimination procedures failed to yield any covariates that significantly improved the model and confirmed that both of these relationships improved the model.
|
View this table: [in a new window] |
TABLE 2. Summary of covariate-model building for the clearance-based modela
|
![]() | (3) |
![]() | (4) |
![]() | (5) |
|
View this table: [in a new window] |
TABLE 3. Parameter and variance estimates for final population pharmacokinetic modela
|
![]() View larger version (25K): [in a new window] |
FIG. 2. Goodness-of-fit plots for the final population pharmacokinetic model. The solid line represents the line of unity for the two upper panels and the 0 reference line for the bottom two panels.
|
Simulated exposure using an alternative dosing regimen. Results of the simulation suggest that, if the dose is not truncated at 400 mg, older (i.e., heavier) pediatric patients are likely to experience gatifloxacin exposure well above the range experienced by younger (i.e., lighter) patients. For those patients with a body weight above 40 kg, the median (minimum to maximum) predicted AUC0-24 using a uniform dose of 10 mg/kg would be 39.5 mg·h/liter (20.3 to 94.2) compared to the median (minimum to maximum) of 29.8 mg·h/liter (12.4 to 84.3) for patients weighing less than 40 kg. The median (minimum to maximum) predicted AUC0-24 using a hypothetical fixed dose of 400 mg for the patients weighing more than 40 kg was 29.8 mg·h/liter (13.9 to 63.2) (see Fig. 3).
![]() View larger version (20K): [in a new window] |
FIG. 3. Scatter plot of the individual AUC estimates (AUC0-24, µg·h/ml) versus patient age (years) for the simulated population. The dose used for this simulation was 10 mg/kg, with an upper limit of 400 mg. The solid line represents the best-fit line of AUC versus age.
|
|
|
|---|
A primary benefit of population analysis methods over traditional methods is the ability to quantify the impact of patient factors on the pharmacokinetic parameters of the drug (5, 31, 32, 33). In this analysis, relationships between body size and gatifloxacin pharmacokinetic parameters were of the most importance. Body surface area best explained the intersubject variability in gatifloxacin clearance, while weight best described the variability in volume. The linear relationship between BSA and apparent clearance is consistent with the predominant renal elimination of gatifloxacin and the flat linear relationship between BSA and the glomerular filtration rate in older infants and children (30). Beyond these size-related relationships, no other covariate, including age, serum creatinine, and estimated creatinine clearance, improved the model. In this case, body size is a likely surrogate for developmental variability in the clearance organ, and no further influence of age could then be identified. While one might expect estimated creatinine clearance to be a significant independent covariate in the model, the patients in this study all had normal renal function, reducing our ability to evaluate the impact of reduced renal function. In addition, all patients had very low serum creatinine values (the maximum value was 1.0 g/dl) reported to one decimal place, reducing the precision of creatinine clearance estimates. Finally, both BSA and creatinine clearance are estimated from subject height; thus, the high degree of correlation seen between these two covariates was expected.
Although BSA was found to be a better predictor of clearance than body weight, the two measures are highly correlated, and body weight is simpler to use for dosing pediatric patients. As the nonlinear relationship between weight and BSA is most pronounced in the youngest patients, there is a potential for underdosing of the smallest children or overdosing of the largest children when using a weight-based dosing scheme in situations where clearance is best correlated with BSA. However, given the inherent interindividual variability in clearance and the relatively wide therapeutic index of gatifloxacin, any added precision gained from BSA-based dosing is unlikely to be clinically significant. A more clinically relevant issue is that of the maximum dose to be used in pediatric patients. In children weighing less than 40 kg, a dose of 10 mg/kg provides exposures similar to those seen in adults from studies utilizing 400 mg. However, children weighing more than 40 kg have higher gatifloxacin exposure (in terms of AUC0-24) than the younger age groups. As can be seen in Fig. 3, limiting the maximum dose to 400 mg once daily brings the gatifloxacin exposure in the older (i.e., heavier) children in line with the exposure seen in the younger (i.e., lighter) children. As the established PK-PD target for efficacy in Streptococcus pneumoniae infections is a fuAUC/MIC of 30 (1, 10, 11, 27), the dose of 10 mg/kg appears to be appropriate from an effectiveness standpoint as well, given the present MIC90 for gatifloxacin in large cohorts of Streptococcus pneumoniae isolates (19).
The results of this analysis demonstrate the applicability of a population-based pharmacokinetic analysis to gatifloxacin data obtained from a pediatric population. Validation of this approach is supported by both internal assessment (i.e., cross-validation) and external comparison to later studies (13, 25). The value of the population-based approach for gatifloxacin in pediatrics, as well as for pediatric drug development in general, resides with the potential to make use of sparse data sets obtainable from future studies to evaluate exposure-response relationships for safety and/or efficacy in pediatric subpopulations defined by disease state and/or development (e.g., patients with renal compromise, neonates, and infants). Such analyses can be extremely useful in providing dose justification and may help to streamline the pediatric drug development process, allowing children earlier access to potentially useful therapeutic agents.
Published ahead of print on 12 January 2007. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»