Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, October 2007, p. 3714-3719, Vol. 51, No. 10
0066-4804/07/$08.00+0 doi:10.1128/AAC.00398-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Pediatric Oncology Branch, National Cancer Institute, National Institutes of Health, Bethesda, Maryland 20892,1 Emerging Pathogen Section, Ordway Research Institute, Albany, New York 12208,2 Center for Cancer and Blood Disorders, Children's National Medical Center, George Washington University School of Medicine and Public Health, Washington, D.C. 20010,3 Warren Albert Medical School of Brown University, Providence, Rhode Island 02903,4 Children's Hospital of Orange County, Orange, California 92868-3874,5 Department of Infectious Diseases, St. Jude Children's Research Hospital, Memphis, Tennessee 38105,6 Department of Pediatrics, Division of Pediatric Hematology/Oncology, Georgetown University Medical Center, Washington, D.C. 20007-2197,7 The Children's Hospital, Denver, Colorado 80218,8 Astellas Pharma US, Deerfield, Illinois9
Received 23 March 2007/ Returned for modification 4 May 2007/ Accepted 9 July 2007
|
|
|---|
|
|
|---|
The echinocandins represent a novel class of antifungal agents with in vitro, in vivo, and clinical activities against the medically important opportunistic fungal pathogens Candida spp. and Aspergillus spp. (10). The pharmacokinetics of echinocandins in pediatric patients have been studied in phase I trials (4, 13, 16, 17). However, there are no population pharmacokinetic models of the echinocandins for pediatric patients. The safety and descriptive pharmacokinetics of micafungin in neonates and children have recently been established in two phase I trials conducted in the United States (13, 16).
Important pharmacokinetic questions remain unresolved for the dosing of micafungin in pediatric patients. Specifically, what is the dosage of micafungin which produces an equivalent magnitude of drug exposure to that observed in adults, and what is the degree of expected pharmacokinetic variability? In order to address these questions, we first described the population pharmacokinetics of micafungin in children aged 2 to 17 years. We then compared three structural pharmacokinetic models, two of which explicitly incorporated weight as a covariate. These data provide further valuable information regarding appropriate pediatric dosing of micafungin and serve as a critical step in the optimization of antifungal therapy for invasive fungal diseases in this vulnerable population.
(Presented, in part, at the 46th Interscience Conference on Antimicrobial Agents and Chemotherapy, San Francisco, CA, 27 to 30 September 2006, abstr. M-299.)
|
|
|---|
Pharmacokinetic models and modeling.
All data were examined using a population methodology with the use of the Big and Little population pharmacokinetic nonparametric adaptive grid with adaptive
(NPAG) programs developed by Leary, Schumitzky, Jelliffe, and van Guilder (15). The inverse of the estimated assay variance was used as the weighting throughout the modeling process.
Three structural pharmacokinetic models were used in this study. The first represents a standard two-compartment pharmacokinetic model with time-delimited zero-order intravenous infusion and first-order elimination from the central compartment. The model is described by the differential equations 1a and 1b, below.
![]() | (1a) |
![]() | (1b) |
The second (equations 2a and 2b) and third (equations 3a and 3b) models, described below, were developed in order to further elucidate the effect of size (measured in terms of weight) on the pharmacokinetics of micafungin. The initial step in the development of these models involved an examination of the Bayesian estimates for clearance obtained from the standard model. These estimates were then plotted against weight, using both linear and logarithmic scales. Since both linear and logarithmic relationships appeared tenable (see Fig. 2, below), linear and power models which explicitly incorporated weight as a covariate were developed. The linear model took the following form:
![]() | (2a) |
![]() | (2b) |
![]() View larger version (11K): [in a new window] |
FIG. 2. Bayesian estimates for clearance obtained from the standard model. There is a weight-dependent increase in clearance, depicted on a linear scale in panel A and on a log scale in panel B. The slope of the regression in panel B approximates 0.75, which has been extensively used as an allometric power scaling parameter in a wide range of biological processes.
|
![]() | (3a) |
![]() | (3b) |
Model evaluation, comparison, and performance.
For all models, the Bayesian estimates for each patient were obtained using the "population of one" utility within NPAG and Big NPAG. The population mean, median, and modal values were evaluated in the maximum a posteriori Bayesian analysis. For each model, scatter plots of the observed-predicted relationships for each individual patient and the population as a whole were examined. Goodness of fit was evaluated on the basis of a visual inspection and coefficient of determination of the observed-predicted scatter plot after the Bayesian step, as well as the log-likelihood values for each model. Models were compared with reference to the respective log-likelihood values, and statistical comparisons were made using the likelihood ratio test, where twice the likelihood difference was evaluated against a
2 distribution with the appropriate number of degrees of freedom. Predictive performance was based upon the weighted mean error and the bias-adjusted weighted mean squared error.
Monte Carlo simulations. Monte Carlo simulations were performed using the allometric power model. The structural models were implemented within the simulation module of the pharmacokinetic program ADAPT II (9), and the full covariance matrix was inserted into subroutine PRIOR of ADAPT II. A subroutine within ADAPT II (courtesy of David D'Argenio, University of Southern California) enabled micafungin to be administered to each of the 9,999 simulated patients on a weight basis (in milligram per kg). For each simulated patient, the weight-based dose of drug in milligrams per kg was converted internally to an absolute dose of micafungin (in milligrams) by multiplying by the simulated weight. Thus, the simulation process mimicked drug administration as it occurred at the bedside in the original clinical trial, in which the dose of micafungin was planned on a mg per kg basis, but the absolute amount of drug administered to each patient was determined with reference to weight. To ensure consistency with the clinical trial, micafungin was infused over 1 h to simulated children. Both normal and log-normal parameter distributions were explored and discriminated on the basis of their ability to recapitulate the original parameter means and their dispersions. All simulations were performed at steady state between days 13 and 14 post-initiation of therapy. The area under the concentration-time curve from 0 to 24 h (AUC0-24) was determined by integration.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Micafungin population pharmacokinetic estimates for children aged 2 to 17 years from the NPAG and Big NPAG analysesa
|
![]() View larger version (12K): [in a new window] |
FIG. 1. Observed-predicted relationships for serum drug concentrations of micafungin after the Bayesian step for the standard model (A), the linear model (B), and the allometric power model (C).
|
|
View this table: [in a new window] |
TABLE 2. Predictive performance of micafungin population models
|
The estimates for the mean and median values for the parameters from the linear and allometric power models are summarized in Table 1. For both models, the fit to the data was excellent (r2 of 0.953 to 0.9600) (Fig. 1B and C), with measures of precision and bias which were comparable to the standard model (Table 2). While the differences were relatively small, the better (more positive) log-likelihood values of the linear and allometric power models compared with the standard model suggested that the incorporation of weight as a covariate added explanatory power to the models. While the predictive performances of each of the three models were comparable, the allometric power model had the largest log-likelihood value, suggesting that it best accounted for the data; for this reason this model was used in the Monte Carlo simulations. The full covariance matrix used for these analyses is shown in Table 3. The mean parameter values and their standard deviations could be recapitulated with a 9,999-patient Monte Carlo simulation in which a log-normal distribution was used for each model parameter (Table 4). The validity of the allometric power model was further confirmed when the values for the constants normalized to a 70-kg adult from the allometric power model (Vstd and SCLstd) fitted to the pediatric data set closely approximated estimates for the volume of distribution and clearance obtained from a separate adult data set (10.43 ± 5.60 liters and 1.17 ± 0.38 liters/h, respectively; T. Gumbo, submitted for publication).
|
View this table: [in a new window] |
TABLE 3. Covariance matrix for micafungin from the allometric power model
|
|
View this table: [in a new window] |
TABLE 4. Parameter values from the 9,999-patient Monte Carlo simulation with the allometric power modela
|
![]() View larger version (22K): [in a new window] |
FIG. 3. Monte Carlo simulations in 9,999 children aged 2 to 17 years receiving 2 mg/kg at steady state. (A) Mean micafungin concentrations along with values for 5% and 95% of simulated patients. (B) Range and distribution of the AUC0-24 that are predicted to develop following the administration of micafungin at 2 mg/kg.
|
![]() View larger version (9K): [in a new window] |
FIG. 4. Monte Carlo simulations in 9,999 simulated patients showing the micafungin dosages in children,weighing 10 to 80 kg, which are required to produce a mean AUC0-24 equivalent to the AUC0-24 observed in adults receiving 100 mg (A), 150 mg (B), and 200 mg (C) at steady state. To achieve an adult equivalent of 100, 150, and 200 mg/day, children aged 2 to 17 years should receive a dosage (in mg/kg) of 3.38 x weight–0.25, 5.07 x weight–0.25, and 6.77x weight–0.25, respectively. As weight declines, a progressively higher dosage is required to produce the same AUC; this phenomenon appears especially apparent in children weighing less than 20 kg. Children weighing 10 kg require nearly double the dosage of micafungin as children weighing 80 kg. A fixed dosage (e.g., 2 mg/kg) would result in lower levels of drug exposure in smaller children.
|
|
|
|---|
Population modeling of pediatric patients presents a number of challenges. Most important is the effect of size on drug disposition and elimination. The sizes of pediatric patients vary by more than an order of magnitude. To enable a more complete understanding of the importance of weight as a covariate, we modeled absolute (rather than weight-adjusted) dosages, volumes, and clearances. The improved log-likelihood values of the two models in which weight was incorporated as a covariate suggested that it accounted for a portion of the residual interpatient variability observed using the standard model (Table 2).
Allometry is the field of study which relates bodily function and morphology to body size (7). An understanding of allometric relationships is important, as adult pharmacokinetic models cannot necessarily be used to predict appropriate dosing regimens for children. An abundance of data suggests that organisms do not exhibit simple geometric scaling; this is a consequence of powerful biological constraints on structure and function that do not allow organisms to maintain the same geometric relationships as size changes (7). The 3/4 power law has been extensively used to model the relationship between physiological functions, such as metabolic rate and clearance, and size (14, 18). Body surface area represents a measure of size and has been used to guide pediatric dosing. Body surface area can be related to weight by using a 2/3 power scaling exponent (14). Predictions of clearance based upon body surface area and the 3/4 power model are similar and certainly superior to assuming that clearance is directly proportional to weight (14). In the current study, the decision to use a 3/4 power model, rather than a surface area model, was based upon two observations: first, many body functions scale predictably with an allometric 3/4 power model rather than with body surface area (7); second, the slope of the regression line in Fig. 2B better approximated 3/4 rather than 2/3. Furthermore, there is a theoretical basis for the allometric 3/4 power law which provides for a deeper understanding for the observation that clearance is proportionally higher in smaller children. There appear to be inherent limitations in the efficiency with which a progressively larger mass of tissue (in this case the liver) can be supplied with nutrients (or in the case of drug clearance, with the drug itself). Consequently, the rate of clearance in larger organisms is slower than predicted on the basis of tissue mass alone (conversely, the rate of clearance in smaller children is higher than predicted from data derived in adults). These concepts are presented in detail elsewhere (7). The nonlinear relationship between weight and clearance (resulting from inherent physiological constraints) means that as weight decreases, progressively higher dosages of micafungin (on a mg/kg basis) are required to achieve equivalent drug exposure; these increased dosages are higher than predicted on the basis of weight alone and are especially so in children weighing <10 to 15 kg.
A reasonable aim of pediatric dosing is to ensure levels of drug exposure which are comparable to those achievable in adults and which approximate those for which antifungal efficacy has been established. The current study demonstrates that smaller children require higher dosages of micafungin to ensure mean levels of drug exposure equivalent to those observed in larger children and adults. This can be achieved by increasing the dosage of micafungin at some (arbitrarily chosen) cutoff weight. In the current study, however, we were able to describe the continuous relationship between dose and weight and use this to precisely define antifungal dosing in pediatric patients. The Monte Carlo simulations also demonstrated that there is a degree of interpatient variability in serum drug concentrations and exposures following the administration of micafungin (Fig. 4); this observation may prompt an increase in dosage in the circumstance of a suboptimal therapeutic response despite the administration of a seemingly appropriate dose of micafungin.
In summary, the population pharmacokinetics of micafungin in children aged 2 to 17 years demonstrate the importance of considering and incorporating weight as a covariate in order to adequately describe drug behavior. The allometric power model developed in this study enables the identification of pediatric dosage regimens of micafungin that, based upon Monte Carlo simulations, result in drug exposure that is equivalent to that observed in adults, for whom antifungal efficacy has been established.
Published ahead of print on 16 July 2007. ![]()
|
|
|---|
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»