Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, November 2008, p. 4043-4049, Vol. 52, No. 11
0066-4804/08/$08.00+0 doi:10.1128/AAC.00569-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Department of Pediatrics, Division of Neonatology, Children's Hospital of Philadelphia, Philadelphia, Pennsylvania,1 Department of Pediatrics, Division of Clinical Pharmacology and Therapeutics, Children's Hospital of Philadelphia, Philadelphia, Pennsylvania,2 Department of Pediatrics, University of Virginia, Charlottesville, Virginia,3 Department of Pediatrics, University of Utah, Salt Lake City, Utah,4 Department of Pediatrics and Duke Clinical Research Institute, Duke University, Durham, North Carolina,5 Department of Pediatrics, University of Louisville, Louisville, Kentucky,6 KAI Research, Inc,7 Hospital for Children and Adolescents and Department of Clinical Pharmacology, University of Helsinki, Helsinki, Finland,8 Metrum Institute, Tariffville, Connecticut9
Received 1 May 2008/ Returned for modification 25 August 2008/ Accepted 9 September 2008
|
|
|---|
|
|
|---|
Fluconazole is a potent antifungal azole drug used for both prevention and treatment of candidiasis (10, 14, 17, and 19; Diflucan package insert, Pfizer, New York, NY). Fluconazole administered either orally or intravenously has an excellent safety profile in adults and children. It is minimally metabolized, is eliminated as active drug in the urine, and effectively penetrates tissues and cerebrospinal fluid. Although it is frequently used in the neonatal population, remarkably little fluconazole PK data have been published for neonates, with studies being limited to a small number of preterm 24- to 29-week-gestation neonates (9, 21, 23). Existing data are thus insufficient for the development of a generalized dosing guidance. We therefore performed a population PK study in young infants who were receiving fluconazole as routine clinical care for prevention or treatment of systemic candidiasis in the neonatal intensive care unit.
Portions of this work were presented at the Society for Pediatric Research Annual Meeting (Honolulu, HI, 2 to 6 May 2008).
|
|
|---|
34 weeks) and postnatal age (PNA) (<14 days and 14 to 119 days). The second study (study 2) was an open label PK study of a panel of antimicrobial drugs being conducted at Duke University. For both studies, fluconazole dosing was determined by the routine clinical practice in each unit and no exclusion criteria were used. These studies were approved by the institutional review boards at each respective institution, and informed consent was obtained. For covariate analysis, we collected the following information: indication for fluconazole, BGA, dates of positive Candida cultures, and daily assessments of PNA, postmenstrual age (PMA), weight, urine output (ml/24 h), and respiratory support. Covariates that exhibited time-dependent changes (e.g., weight, age, and creatinine) were permitted to change with time, and the actual value in the data set reflects the observations made at each patient visit. Missing weights were imputed with the last recorded value carried forward for up to 7 days.
Serum creatinine (SCRT) was recorded when obtained for clinical care. Infants are more likely to have SCRT values measured when they have clinical signs of renal insufficiency. Infants with no measured creatinine values were assumed to have a SCRT value of
1.0 mg/dl. To avoid the impact of missing SCRT data, we created a dichotomized variable, CR, and assigned a value of 0 for time points when the SCRT level was never measured or when it was
1.0 mg/dl. Likewise, a CR value of 1 was designated for infants with a measured SCRT of >1.0 mg/dl, and a SCRT value of > 1.0 was imputed for up to 7 days. We used the exponentiated expression SCRT(CR x
) to switch between the renal function subtypes.
PK sample collection. Infants in study 1 were randomly assigned to one of two PK sampling schedules (schedule A included preinfusion, end of infusion, and 1, 6 to 8, 24, and 48 h postinfusion; schedule B included preinfusion, end of infusion, and 3, 10 to 12, 24, and 48 h postinfusion). Sampling was initiated as soon as informed consent was obtained, typically corresponding to the first through fifth dose. If an infant remained on fluconazole, trough plasma samples were obtained on days 7, 14, and 21 (±2 days). Each sample was 300 µl of blood in EDTA Microtainers. Plasma was separated and removed within 30 min and stored at –20°C. To supplement these PK samples, plasma from residual discarded clinical specimens timed per clinical routine practice (EDTA Microtainers) was collected up to 72 h after the blood was collected. These scavenged plasma samples were frozen at –20°C. Infants in study 2 had samples scavenged from discarded blood in the clinical laboratory. Samples from all sites were shipped on dry ice to Children's Hospital of Philadelphia, where they were stored at –70°C for analysis. To evaluate fluconazole stability in scavenged samples, we determined the stability of fluconazole (0.03 to 9 µg/ml) when stored in whole blood in EDTA Microtainers. Fluconazole plasma concentrations were measured after whole blood triplicate samples were stored for 0.5, 24, 48, and 72 h at room temperature.
Liquid chromatography-tandem mass spectrometry assay for detection of fluconazole in neonates. We developed and validated an analytical method for fluconazole detection in human plasma suitable for the small plasma volumes obtainable in neonates. Liquid chromatography-tandem mass spectrometry analysis was conducted on an API 4000 Q TRAP (Sciex, Toronto, Canada) coupled with a Shimadzu high-performance liquid chromatography (HPLC) system (Kyoto, Japan), using an electrospray ionization source in the positive mode and under the following conditions: curtain gas, 19; gas 1 (nebulizer gas), 35; gas 2 (heater gas), 80; CAD gas, medium; TurboIonSpray voltage, 4,250 V; entrance potential, 11 V; collision energy, 25 V; source temperature, 550°C; and dwell time, 200 ms. The optimized declustering potential and collision cell exit potential were set at 51 and 5 V, respectively. The collision energy was optimized based on the individual fragmentation selected to obtain the most intense precursor to product ion transitions. Fluconazole-D4, with a molecular weight of 311.20, was obtained from SynFine (Ontario, Canada). HPLC separation was performed using a gradient mobile phase of acetonitrile and water in 0.1% formic acid on a Waters XTerra C18 HPLC column (100 by 2.1 mm; 3.5 µm) (Milford, MA) with a Waters XTerra guard column (10 by 2.1 mm; 3.5 µm) (Milford, MA) at a flow rate of 0.2 ml/min. Multiple reaction monitoring was used to detect fluconazole and fluconazole-D4 (internal standard) at 307.20/219.80 and 311.20/222.80 m/z, respectively. Analytical data were acquired by Analyst software (version 1.4.1). The lower limit of quantitation of fluconazole in plasma was 0.01 µg/ml. Intraday and interday coefficients of variation are <8.1% at concentrations ranging from 0.01 to 10 µg/ml.
Population PK analysis. Pharmacokinetic data were analyzed with a nonlinear mixed effect modeling (NONMEM) approach using the computer program NONMEM (version 5) (3). The first-order conditional estimation method was used for all model runs. One- and two-compartment structural models were evaluated. Interindividual random effects were evaluated on clearance (CL) and volume of distribution (V). Covariance was described by a block Omega matrix. We used an exponential model for interindividual variance for CL and V. A combined additive and proportional error model was deemed appropriate to describe residual variability. The potential impact of physiologically plausible, clinical covariates on PK parameters was explored in a forward stepwise manner as follows: weight (kg), allometric scaling, BGA (weeks), PNA (week of life; defined as day of life/7), PMA (defined as BGA plus PNA in weeks), and SCRT. Covariates were retained in the final model if there was an improved goodness of fit and a decreased minimum objective function. A drop in objective function of >10.83, considered significant at P < 0.001, was used to discriminate among alternative nested models. Continuous covariates were scaled to their median values. Empirical Bayesian estimates of individual infant PK parameters were generated from the final model using the POSTHOC subroutine.
Model evaluation. Models were evaluated based on the following six criteria: (i) successful minimization, (ii) goodness of fit as assessed by the Akaike information criterion (18), (iii) diagnostic plots, (iv) precision of parameter estimates, (v) an internal quantitative predictive check, and (vi) an external visual predictive check. We assessed precision of the final population PK model parameter estimates using stratified nonparametric bootstrapping (1,000 replicates) to generate the 95% confidence intervals (CIs) for parameter estimates. For the internal predictive check, the quantity of interest was the average observed fluconazole concentration in each infant. Summary metrics of the quantity of interest were calculated for the observed values in this study cohort and compared to values obtained from 1,000 Monte Carlo simulation replicates of the original data set using the final population PK model. For the external predictive check, the final model was used to generate 100 Monte Carlo simulation replicates of fluconazole exposure in the fluconazole PK preterm neonatal cohort described by Saxén et al. (21), and simulated results were compared with those observed in the study. For prediction into the Saxén data set, we excluded SCRT from the model, since only a single SCRT value was available per infant and a different creatinine assay was utilized for babies in the Saxén study (21).
Assessment of dose-exposure relationship.
Monte Carlo simulations using the final population PK model were used to explore the impact of PNA and BGA on dose-exposure relationships. For target exposure, we chose an area under the concentration-time curve from 0 to 24 h (AUC0-24) of 800 mg · h/liter. This AUC is consistent with exposures in adult intensive care units or immunocompromised patients treated with fluconazole at 800 mg/day and ensures that exposure exceeds the pharmacodynamic target of an AUC/MIC value of >50 for Candida species with a MIC of 8 µg/ml at the CLSI sensitivity breakpoint (2, 12, 13, 19). We determined the interval AUC24 to be the AUC for each 24-h interval after a given dose calculated within the NONMEM control file and derived from the differential equations defining the model. These simulations performed with NONMEM included infants with demographic characteristics of the study population and infants exposed to 12 mg/kg/day fluconazole for 14 days. We excluded exposure predictions when an infant's SCRT was
1.3 mg/dl because these infants would have a dose adjustment for renal insufficiency. The dose required to achieve a target steady-state AUC24 was calculated by multiplying the individual model-predicted CL (liter/kg · h) by the target AUC24 (800 mg · h/liter).
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Baseline characteristics of infants
|
Population PK model construction. A one-compartment model was the appropriate structural model for this data set based on diagnostic plots (see Fig. 2), and a large coefficient of variation around estimates of intercompartmental clearance. We began with a base model of CL and V with weight (1) (Table 2). In the base model, scatter plots revealed correlations between interindividual variances on CL (ETA1) and V (ETA2); therefore, a covariance term was added to the model. The residual unexplained variability of CL (ETA1) showed a correlation with indices of age and maturity (Fig. 1). Steps in subsequent model building are shown in Table 2. CL estimated by the model with BGA and PNA together was superior to CL estimated with PMA. No additional maturity affects were added to the model for V since no additional decline in partitioning of variance was seen.
![]() View larger version (13K): [in a new window] |
FIG. 2. Final model diagnostic plots of observed versus predicted fluconazole concentration and weighted residuals. For population prediction (A) and individual prediction (B), the individual data points are labeled by individual ID numbers with a dotted line connecting individual data points. Line of identity is included as a reference. For weighted residuals (C), individual data points are labeled by individual ID numbers with a dotted line as a smoothing spline trend line through the data. A solid line at y = 0 is included as a reference.
|
|
View this table: [in a new window] |
TABLE 2. Construction of population model
|
![]() View larger version (22K): [in a new window] |
FIG. 1. Scatter plot of variance term for CL (ETA1) using the base model. Correlation between variances on CL (ETA1) and the following variables: variance on V (ETA2) (A), BGA (B), PNA (C), and PMA (D).
|
We evaluated the impact of scavenged specimen acquisition with delayed separation of plasma from whole blood. Observed fluconazole concentrations from scavenged samples were indistinguishable from concentrations in PK samples on plots of observed versus predicted concentrations. We also confirmed fluconazole stability in whole blood stored at room temperature in the laboratory (data not shown). In the random effects error model, we incorporated a term (SCAV) to be equal to 1 for samples collected from scavenged discarded blood (Table 2). Scavenging introduced minimal bias and appeared to slightly underestimate fluconazole concentrations by 4% (95% CI, –11% to 2%) (Table 3).
|
View this table: [in a new window] |
TABLE 3. Final model population pharmacokinetic parameter estimates
|
![]() View larger version (20K): [in a new window] |
FIG. 3. Predictive check for plasma fluconazole concentration by day of life. Observed fluconazole concentrations (open circles) in a previously published Saxén trial (21) of twelve 26- to 29-week-gestation infants receiving 6 mg/kg fluconazole on day of life 0 (birthday), 3, 6, 9, and 12 are superimposed with the median (dotted line) and 90% population prediction interval (shaded area) for fluconazole concentrations from 100 Monte Carlo simulated trials, given the final model and parameters. Bold arrows indicate fluconazole doses that preceded plasma PK samples. Open lines indicate fluconazole doses that were not followed by plasma PK samples.
|
Dose-exposure relationship.
Monte Carlo simulation was used to evaluate the impact of prematurity, weight, and PNA on the dose-exposure relationship (Fig. 4). We predicted the concentration of fluconazole achieved after a 12-mg/kg daily dose regimen in infants with the age, weight, and prematurity represented in this cohort. For these simulations, we excluded exposure estimates when creatinine was
1.3 mg/dl since these infants would require a dose adjustment for renal insufficiency. Steady-state AUC is not achieved during the first week (Fig. 4A and B). The average half-life is 30 and 50 h for 23- to 29-week- and 30- to 40-week-gestation infants, respectively. Infants who were less than 30 weeks of gestation achieved a median AUC closer to the target AUC of 800 mg · h/liter (Fig. 4A); infants who were more than 30 weeks of gestation achieved a lower median AUC of 400 mg · h/day (Fig. 4B). In order to achieve the 800 mg · h/liter AUC target, dosing will need to be adjusted for BGA and PNA (Fig. 4C).
![]() View larger version (15K): [in a new window] |
FIG. 4. Evaluation of dose exposure relationship. The median (dark line) and population predicted interval from the 10th percentile to the 90th percentile (shaded area) for fluconazole AUC from 100 Monte Carlo simulated trials, given the final model and parameters. Simulated 24-h interval AUC for each day of therapy among 23- to 29-week-gestation infants (A) or 30- to 40-week-gestation infants (B) receiving 12 mg/kg/day fluconazole. (C) Median box plot of predicted dose required to achieve steady-state AUC target of 800 mg · h/liter in infants stratified by BGA and PNA.
|
|
|
|---|
We developed a population PK model for fluconazole disposition in preterm and term infants less than 90 days of age who were receiving fluconazole as part of clinically indicated treatment. Population PK analysis allowed the use of limited sampling schemes and the varied dosing regimens used in routine care. Stratified enrollment ensured broad distribution of both BGAs and PNAs. Randomization to two sampling time schedules allowed for more dense sampling times. Because fluconazole is a stable drug, we were able to increase the number of samples available for analysis by scavenging discarded plasma from clinical specimens. Our robust, sensitive quantitative assay allowed us to accurately measure fluconazole in microvolumes of plasma.
In this model, fluconazole CL increased with allometrically scaled weight (1), BGA, and PNA. PMA, however, did not perform as well as the combination of BGA and PNA. This is not surprising because PMA cannot distinguish a 3-week-old, 24-week-gestation infant from a 1-day-old, 27-week-gestation infant. The model was rigorously evaluated, showed good precision and minimal bias, and was able to estimate into a previously described PK data set in preterm infants (21).
Our model suggests that dose adjustment will be needed to account for the significant changes in fluconazole CL that occurs as a function of BGA and PNA. Although pharmacodynamic efficacy exposure targets for infants have not been described, critically ill or immunocompromised adults with invasive candidemia are typically treated with 800 mg/day of fluconazole to target a steady-state AUC of 800 mg · h/liter. This is double the standard dose required to achieve a pharmacodynamic target AUC24/MIC24 of >50 for Candida, with the CLSI sensitivity breakpoint MIC being
8 (2, 12, 13, 19). The higher plasma exposure is a reasonable target to ensure adequate tissue exposure in patients with widespread, invasive disease (11, 19, 20). Plasma concentrations may not always reveal tissue concentration, and tissue concentrations were not available.
For preterm infants with invasive candidiasis, we chose this target AUC of 800 mg · h/liter because these infants have immature immune systems, have a high projected mortality rate, and often have disseminated disease, including meningoencephalitis (7). A loading dose would be necessary to reach the desired steady-state concentration rapidly. Our data suggest that the lower fluconazole dosages frequently used in the neonatal population result in underexposure and may partially explain the prolonged periods of candidemia (5) and episodes of breakthrough candidemia (23).
Although the neonatologists caring for these infants reported no adverse events, and none of the infants developed a breakthrough fungal infection, a limitation of this study is that comprehensive efficacy and safety data (e.g., liver function tests results, electrocardiograms, etc.) were not collected and thus could not be linked to exposure. Furthermore, we were unable to evaluate drug-drug interactions.
Longitudinal prediction of CL changes with PNA beyond 14 days is difficult. Infants who enrolled when they were more than 14 days old were often receiving fluconazole during periods of clinical deterioration secondary to sepsis or necrotizing enterocolitis. This model may underpredict CL beyond 14 days of age in infants who are relatively well and receiving fluconazole prophylaxis from birth. We plan to perform sparse PK sampling in infants who will be enrolled into a multicenter trial of fluconazole prophylaxis to improve our longitudinal prediction of changes in CL with PNA in these relatively well infants.
Conclusions. Pharmacokinetic neonatal study designs that leverage routine clinical care can provide reliable PK information. The fluconazole population PK model from this nontraditional design was able to accurately predict the observed fluconazole concentrations in a historic, more traditional PK study. Incorporating both BGA and PNA improved the model fit beyond that achieved with PMA. BGA and PNA effects on CL are significant, and dose adjustments will likely be needed. Our future studies will use Monte Carlo simulation to fully explore dose-exposure relationships, evaluate the effects of a loading dose, and guide further development of pharmacologically rational neonatal dosing recommendations.
This project was supported by grant no. 1U10-HD037255-06 (K. C. Wade, J. S. Barrett, and P. C. Adamson) and 1U10-HD45962-05 (D. K. Benjamin, Jr.) from the National Institute of Child Health and Development Pediatric Pharmacology Research Unit. Additional support at Children's Hospital of Philadelphia was provided by the University of Pennsylvania Clinical and Translational Research Center (CTRC) UL1-RR-024134 from the National Center for Research Resources.
The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Center for Research Resources, the National Institute of Child Health and Development, or the National Institutes of Health. No potential conflicts of interest relevant to this article were reported.
The following investigators, research team members, and hospitals participated in the fluconazole PK study in the Pediatric Pharmacology Research Network Unit of the NICHD. Study sites are listed according to the number of children they evaluated. David A. Kaufman, Audrey Meyers, and Marci Williams, University of Virginia Hospital, Charlottesville, VA; Robert Ward, Stacy Morgan, and Jeanne Francis, Primary Children's Hospital, Salt Lake City, UT; Kelly Wade, Tonia Morrison, and Dustin J. Paul, Children's Hospital of Philadelphia, Philadelphia, PA; Dan L. Stewart, Janice Sullivan, and Elizabeth McDowell, Kosair Children's Hospital, Louisville, KY; Kristine Palmer, Laura James, and Angela Riggs, Arkansas Children's Hospital, Little Rock, AR; Jacob Aranda, Ginger Steinhilber, and Deanna Sypula, Michigan Children's Hospital, Detroit, MI; Mohan P. Venkatesh, Ann R. Stark, and Karen Jones, Texas Children's Hospital, Houston, TX; John van den Anker, Louis Scavo, and Elaine Williams, Children's National Medical Center, Washington, DC; Pablo J. Sanchex, George H. McCracken, Jr., and Luz Muniz, University of Texas Southwestern Medical Center, Dallas, TX. The fluconazole samples from the PPRU Antimicrobial Study were provided by Daniel K. Benjamin, Jr., P. Brian Smith, and Kim Fisher at Duke University Medical Center. Rene Kozloff, Nadia Ramey, and Celeste Crouse provided data management support at KAI Research, Inc.
Published ahead of print on 22 September 2008. ![]()
Members of the PPRU fluconazole study group are listed in Acknowledgements. ![]()
|
|
|---|
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»