Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, March 2006, p. 935-942, Vol. 50, No. 3
0066-4804/06/$08.00+0 doi:10.1128/AAC.50.3.935-942.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Peter J. Shaw,2
Christa E. Nath,3
Satya P. Yadav,2
Katherine R. Stephen,2
John W. Earl,3 and
Andrew J. McLachlan1*
Faculty of Pharmacy, University of Sydney,1 Department of Oncology,2 Department of Biochemistry, the Children's Hospital at Westmead, Sydney, NSW, Australia3
Received 21 June 2005/ Returned for modification 30 September 2005/ Accepted 27 December 2005
|
|
|---|
|
|
|---|
Children with malignant diseases who experience prolonged periods of myelosuppression due to cytotoxic chemotherapy or hematopoietic stem cell transplant are highly susceptible to invasive fungal infections (11). The uses of indwelling vascular catheters and broad-spectrum antibiotics further increase the risk of immunosuppressed children developing invasive fungal infections. The broad-spectrum antifungal activity and improved therapeutic properties of L-AmB may provide an important therapeutic advance in the treatment and prophylaxis of invasive fungal infections in pediatric oncology patients. However, there remains debate about the appropriate dosing regimens in children. No single study has prospectively investigated the pharmacokinetics of L-AmB that might guide rational design of dosing regimens in this special population.
The clinical pharmacokinetics of L-AmB have been investigated in healthy adult volunteers after a single intravenous dose (5) and in adult patients after multiple doses (26). All of the studies utilized a traditional pharmacokinetic data analysis approach. However, in the clinical setting where pediatric patients or those with severe illness are studied, only limited numbers of samples can be obtained due to ethical, logistical, and medical considerations. Therefore, the population approach was employed in the present study as it provides not only a solution for combining concentration information from both intensive and sparse sampling but also the valuable knowledge of the population distribution of pharmacokinetic parameters and factors influencing the variability in these parameters.
The aim of this study was to explore the pharmacokinetics of L-AmB in pediatric oncology patients by using a nonlinear mixed effects modeling approach in order to estimate the typical population pharmacokinetic parameters, the variability between patients, and the variability between occasions and within patients and to identify the covariates that are significant predictors of variability in L-AmB pharmacokinetics. As the developed population pharmacokinetic model will be used for the dosing strategy optimization, the stability and predictive performance of the model were evaluated using a bootstrap sampling technique.
|
|
|---|
Pharmacokinetic sampling. L-AmB (AmBisome; Gilead Sciences, VIC, Australia) was administered once daily as a 1-h infusion through one of the lumens of a central line catheter for each patient. The first day of L-AmB administration was defined as day 1. Heparinized whole-blood samples (1 to 2 ml) were collected from the opposite lumen to which L-AmB was delivered and after 5 ml blood was withdrawn and discarded to avoid contamination. The times of dosing and sample collection were accurately recorded on a sample collection sheet. Multiple blood samples, to characterize the full pharmacokinetic profile, were taken in 22 out of 48 courses at specific times in one dosing interval, including immediately before the infusion and 1, 2, 4, 7, 12, and 24 h after the infusion. Trough and peak samples were obtained from 25 patients on day 7 of therapy or, if those samples were not available, on day 6 or day 8. The remaining pharmacokinetic samples were collected with those for routine hospital care throughout the period of L-AmB therapy in order to minimize the risk of infection due to blood sampling from the central line catheter.
Plasma fractions were separated by centrifugation at 1,500 x g for 10 min at 4°C (Beckman CS-15R; Beckman Instruments, California) and were stored at 80°C until analysis.
Drug assay. High-performance liquid chromatography was employed for the quantification of AmB in plasma (5). Briefly, the internal standard p-nitroaniline at a concentration of 27 µg/ml (30 µl) and methanol (200 µl) were added to an aliquot of plasma sample (100 µl) in a microcentrifuge tube. The tube was vortex mixed for 1 min and stood in the dark for 10 min to allow protein precipitation. After centrifugation at 12,100 x g for 10 min (Eppendorf MiniSpin; Biolab, Australia), an aliquot of 20 µl of methanol extract was injected directly into the high-performance liquid chromatography system for analysis. During the extraction procedure, the liposomal bilayer was disrupted so that the total AmB concentration could be measured.
The separation was performed isocratically on a Phenomenex Spherex C18 (150 by 4.6 mm, 5 µm) analytical column, which was coupled by an Alltech Alltima C18 (7.5 by 4.6 mm, 5 µm) guard column. The mobile phase, consisting of 10 mM sodium acetate buffer, including 10 mM EDTA (pH 3.6) and acetonitrile (650:350, vol/vol), was delivered at 1.0 ml/min. Peaks were detected at a wavelength of 405 nm. The assay was linear up to 80 µg/ml, with a limit of quantitation of 0.1 µg/ml. The overall deviation from the true concentration was less than 11.0% for within run and 6.0% for between runs, and the precision was less than 4.5% for within run and 5.0% for between runs. The average recovery was 88.4% at the concentrations of quality control samples with a standard deviation of 6.2%.
Population pharmacokinetic analysis. A population pharmacokinetic model was developed and fitted to AmB concentration-time data by using NONMEM (version 5.1.1; NONMEM project group, University of California, San Francisco, CA). A Compaq visual Fortran compiler (version 6.6) was used. The first-order conditional estimation method was used throughout the model-building procedure. The modeling approach was implemented in a series of steps which are outlined below.
(i) Development of base population pharmacokinetic model.
The base pharmacokinetic model was defined as the model that best described the L-AmB concentration-time data without consideration of covariate effects. One-compartment and two-compartment open pharmacokinetic models with zero-order input were evaluated using subroutines from the NONMEM library (6). The exponential error model was employed to model the interindividual variability (IIV) for the pharmacokinetic parameters defined as follows:
![]() |
i represents the pharmacokinetic parameter for the ith individual,
is the typical value of pharmacokinetic parameter
in the population (e.g., population mean),
i quantifies the deviation of
i from
with a distribution of (0,
2). The difference between the jth observed concentration (Y) in the ith individual and its respective prediction (
) was modeled with a combined additive and proportional error model:
![]() |
1 (proportional component) and
2 (additive component) are random effects quantifying the errors between Y and
with a distribution of (0,
2). Between-occasion variability (BOV) for the parameter
was evaluated as an additional level of random effect and expressed as follows:
![]() |
1...n are random variables assumed to be normally distributed with a mean of 0, and variance is denoted by
2 (15). In the present modeling analysis, each "occasion" was defined as a period of time in which patients were administered L-AmB on a daily basis and was therefore also referred to as a "course" of therapy. The time interval between two occasions was at least two weeks. BOV was tested on clearance (CL) and volume of distribution in the central compartment (V1) and was estimated using NONMEM's BLOCK (2) SAME option, as it was assumed that BOV was the same for all occasions.
(ii) Development of full population pharmacokinetic model.
Factors, including body weight (WT) and age, were investigated by using a covariate model. Taking the covariate WT as an example, the typical model would be
![]() |
is the model predicted value for
given the covariate value for WT;
1 represents the population central tendency for the pharmacokinetic parameter
and corresponds to
evaluated at the value of WTmedian for WT;
2 represents the effect of WT on
; WTmedian is the median weight of the study population. The possible influence of various covariates on the pharmacokinetic parameter of L-AmB was modeled in a multiplicative fashion (14). Parameters were centered on the typical value (median) of each covariate in the study population. The relationships between pharmacokinetic parameters and covariates were evaluated using univariate analysis in NONMEM. Discrimination between hierarchical models was based on the objective function value (OFV) by using the log-likelihood ratio test (22). A P value of 0.001, corresponding to the decrease in OFV by more than 10.83 (df = 1), was chosen to decide the full model. (iii) Full model evaluation. The accuracy and robustness of the full model were evaluated by using a bootstrap procedure, which was performed in an automated fashion by using the bootstrap option in the software package Wings for NONMEM (WFN, version 408; N. Holford, University of Auckland) (20). The results from 1,000 successful runs were obtained. The mean and 95% predictive interval (PI) were reported, and the difference between the bootstrap mean and the estimate obtained from original data set for all the population parameters was calculated.
Pharmacodynamic analysis of L-AmB in antifungal therapy. Two PK-pharmacodynamic (PD) parameters, peak concentration at steady state in relation to the MIC (Cmax,ss/MIC) and area under the concentration-time curve at steady state in relation to the MIC (AUCss/MIC), were employed to evaluate their relationship with clinical antifungal efficacy of L-AmB (1). Cmax,ss was predicted using the final population pharmacokinetic model performed in NONMEM by assuming that the steady state of multiple daily dosing of L-AmB is achieved by day 7 (24). AUCss was estimated as dose/CL, where CL is the individual posterior Bayesian estimate of L-AmB clearance and dose is the dose of L-AmB on the 7th day of the course. The determination of MICs of fungal pathogen was performed in the Pacific Laboratory Medicine Services (Royal North Shore Hospital, NSW, Australia) by testing fungal species, which were isolated from individual patients who had proven invasive fungal infections, against D-AmB. The assessment of proven invasive fungal infection and antifungal efficacy were based upon the criteria proposed by the European Organization for Research and Treatment of Cancer/Invasive Fungal Infections Cooperative Group (EORTC/IFICG) (4). SPSS (version 11.5; SPSS, Inc., Illinois) was used to carry out the statistical analyses. A P value of less than 0.05 was considered to be statistically significant.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Patient characteristics in population pharmacokinetic study of L-AmBa
|
Population pharmacokinetic modeling.
The population pharmacokinetic model-building process is presented in Table 2. A two-compartment pharmacokinetic model with zero-order input and first-order elimination was found to describe L-AmB-observed concentration-time data better than a one-compartment model as it markedly decreased the OFV (
OFV, 673.8), leading to its use for the following model-building procedures.
|
View this table: [in a new window] |
TABLE 2. Pharmacokinetic model-building summary for L-AmB
|
The individual posterior Bayesian estimates of pharmacokinetic parameters that were generated from the base model were plotted against weight, height, age, and sex. The Mann-Whitney U test indicated that sex did not significantly affect the pharmacokinetics of L-AmB. Since weight and height were highly correlated (r2 = 0.93) and weight alone yielded a better correlation with CL and V1 than did height, weight and age were evaluated as potential covariates.
In the first step of the model-building procedure, each single combination of a covariate with one of the primary pharmacokinetic parameters (CL and V1) was included sequentially into the base model. The covariate effects on Q (intercompartmental clearance) and V2 were not evaluated because of a lack of a physiological rationale and clinical meaning. At the significance level of P < 0.05, weight had significant effect on both CL and V1 and age was found to significantly influence V1. The base model with the addition of weight as a covariate for CL was used as the reference model for the next step of the forward addition procedure.
In the second step of the model selection procedure, further inclusion of weight as a covariate for V1 produced a significant decrease in the OFV by 13.7 points compared to that produced by the reference model (P < 0.001). When this model was used as the reference model, the inclusion of the covariate age on V1 increased the OFV. Therefore, the inclusion of patient weight rather than age in the analysis provided the largest effect on the goodness-of-fit of the model. The structure of the full model is presented in Table 2 as model 9.
Table 3 summarizes the population parameter estimates obtained from the full population pharmacokinetic model. Figure 1 displays the diagnostic plots for the full model. All of the pharmacokinetic parameters were reliably estimated as the percentage of relative standard errors (RSE) was less than 50%. No trend was observed in the weighted residual plots, and the population-predicted concentrations were symmetrically distributed around the line of identity, indicating that the model adequately describes the pharmacokinetic profile of L-AmB in all the doses studied. The representative concentration-time profile in one patient treated with L-AmB at a dose of 1.8 is shown in Fig. 2, in which the model fit appears to be adequate.
|
View this table: [in a new window] |
TABLE 3. L-AmB population parameter estimates from full model
|
![]() View larger version (19K): [in a new window] |
FIG. 1. Diagnostic plots of full population pharmacokinetic model. Upper panel, observed AmB concentration (conc) versus population-predicted concentration (mg/liter) (r2 = 0.7159; the dash lined is the line of identity). Lower panel, weight residuals (WRES) versus time after first dose (h). Time is presented with reference to the first L-AmB dose in each study course.
|
![]() View larger version (12K): [in a new window] |
FIG. 2. Representative concentration (conc)-time profile of a patient treated with a 1-h intravenous infusion of L-AmB at a dose of 1.8 mg/kg over a dosing interval. Closed circles represent the individual observed concentrations, and the solid line represents the population-predicted concentrations obtained from the full pharmacokinetic model.
|
![]() |
![]() |
For the population investigated, the typical value of CL and V1 was 0.44 liters/h and 3.12 liters for a patient weighing 21 kg. CL and V1 increased with every 10 kg of weight by 16 and 27%, respectively. Compared with the base model, the IIV of CL was decreased from 24 to 10% but the IIV for V1 failed to be well estimated, as its
2 output was very low. This should not be interpreted as an absence of variability in V1 but instead was a consequence of the fact that NONMEM had difficulty in characterizing the differences in mean parameter values for individuals when BOV was much higher than IIV (15). The estimate of the proportional term for the residual variability, expressed as the percentage coefficient of variation, was 27% and the estimate of the additive component of residual variability was 0.034 mg/liter.
Full model evaluation. From the original data set, 1,000 replicate data sets were generated and used for the evaluation of the stability of the full model and accuracy of parameter estimates. Minimization successful and minimization terminated due to rounding errors with numbers of significance above two were included in the bootstrap calculations. Table 4 lists the results of the bootstrap procedure, presented as mean and 95% PI. The mean values of the bootstrap procedure were comparable to the parameter estimates from the original data set, indicating that the accuracy and robustness of the full model was acceptable. The secondary pharmacokinetic parameters of L-AmB, including t1/2ß (59.4 ± 36.5 h) and Vss (1.14 ± 0.68 liters/kg), were generated using the individual posterior Bayesian estimates of parameters obtained from the full model.
|
View this table: [in a new window] |
TABLE 4. Evaluation of the full model using the bootstrap procedure
|
|
View this table: [in a new window] |
TABLE 5. Microbiological, pharmacokinetic, and pharmacodynamic information for patients with proven fungal infection
|
![]() View larger version (13K): [in a new window] |
FIG. 3. The ratio of peak concentration at steady state to the MIC of the patients with partial response (left box, 40.2 ± 13.3, n = 6 cases) versus that of the patients with complete clinical response (right box, 67.9 ± 17.5, n = 4 cases).
|
|
|
|---|
The inclusion of the potential covariates in the model (e.g., patient demographic characteristics) was performed according to a stepwise forward addition procedure. Patient clinical characteristics, including the status of fungal infections and the clinical responses to L-AmB treatment, were not included into the covariate model-building process as these characteristics putatively have no relevance to the pharmacokinetics of L-AmB. Since renal elimination plays a minor role in the excretion of L-AmB (5), markers of renal function were not considered as potential covariates. It was found that body weight was an influential covariate, as it explained part of the variability in L-AmB pharmacokinetics with statistical significance. The incorporation of weight into the pharmacokinetic model indicated that weight-related changes in pharmacokinetic parameters would result in the alteration of drug exposure in patients receiving a given dose. A series of simulations were performed using the population pharmacokinetic parameters and variability to explore the Cmax,ss of L-AmB in a range of doses from 1 to 12.5 mg/kg for patients with weight ranges from 10 to 70 kg. It was found that achieving a comparable Cmax,ss in younger patients whose weight was less than 20 kg would necessitate the administration of a higher dose, suggesting that the current dosing regimens of L-AmB are likely to be suboptimal in this patient cohort.
The clearance and volume of distribution of L-AmB in the central compartment exhibited significant BOV with estimated values of 46 and 56%, respectively. The reason for this observation is unclear but is likely to reflect physiological changes that are occurring in this population over the period of the study. Two mechanisms may contribute to the observed BOV. Firstly, lipid-based formulations influence the association of drugs with lipoproteins in blood, which in turn may be significantly influenced by disease (27). Therefore, it is possible that changes in the lipoprotein profile of patients between occasions could lead to alterations in the pharmacokinetics of AmB. Secondly, dose-dependent hepatic uptake of L-AmB, which we have observed in the isolated perfused rat liver (13), suggests the total amount of drug a patient receives on occasion may impact on the pharmacokinetics of L-AmB. Karlsson and Sheiner (15) indicated that the relative sizes of IIV and BOV have implications for therapeutic drug monitoring (TDM) and concluded that the BOV of pharmacokinetic parameters needed to be lower than those of IIV for carrying out meaningful TDM. The observation that a drug in a particular patient cohort has a larger BOV than IIV suggests that the concentrations measured in the first occasion are not likely to inform the dose individualization in the subsequent occasions. A study which (i) has sufficient patients contributing data for more than one occasion and (ii) has the data on each occasion being informative about the parameter of interest would improve the estimate of IIV when BOV is higher than IIV. Further work is needed to clarify the possible impact of this observation on the TDM of L-AmB in children.
The encapsulation of AmB into liposomes significantly alters the pharmacokinetics of AmB. The clearance and volume of distribution in the central compartment decreased compared with those reported by our research group in pediatric patients who were treated with D-AmB (CL, 0.44 versus 0.88 liters/h and V1, 3.12 versus 9.97 liters) (19), suggesting a more rapid distribution of D-AmB into the tissues than L-AmB. Consistent with these observations, L-AmB exhibits higher plasma peak concentrations and prolonged retention in plasma compared with those that were attained when clinically relevant doses of D-AmB are administered. Taken together, the unique pharmacokinetics of L-AmB represent a potential clinical advantage for antifungal therapy with L-AmB as these pharmacokinetic properties can facilitate the localization of the drug to the sites of infection (9).
Three pharmacokinetic-pharmacodynamic parameters, including the percentage of time that drug concentration exceeds the MIC (%T > MIC), the Cmax/MIC, and the AUC/MIC, have been used to describe the association between drug exposure and antifungal effect (1, 3). For antifungal drugs exhibiting concentration-dependent killing and long postantifungal effects, such as polyenes, the PK-PD parameters that characterize this type of dosing strategy include Cmax/MIC and AUC/MIC. However, in this study, because the patients with proven fungal infection received a range of dosing regimens on each occasion, it was not possible to provide one meaningful estimate of AUC and AUC/MIC. In this relatively small patient cohort, Cmax/MIC, rather than AUC/MIC, was chosen as a potential PK-PD parameter which is predictive of clinical outcomes. A previous PK-PD study in neutropenic disseminated-candidiasis mice treated with D-AmB has demonstrated that Cmax/MIC was the PK-PD parameter best predictive of antifungal efficacy (as measured by a log reduction in CFU) relative to AUC/MIC and %T > MIC, and the breakpoint of Cmax/MIC equal to 10 was associated with the achievement of maximum antifungal efficacy (2). Furthermore, Groll et al. investigated the relationship between drug concentrations and the antifungal efficacy of each of four formulations of AmB, including D-AmB, AMB colloidal dispersion, AMB lipid complex, and L-AmB in a rabbit model of central nervous system candidiasis. A strong inverse correlation between Cmax/MIC (62.0 ± 1.3) and residual fungus burden in brain tissue was observed in L-AmB-treated group (r = 0.792; P < 0.0001) (12). Consistent with the animal studies, the present study demonstrated an impact of Cmax,ss/MIC on the clinical antifungal efficacy of L-AmB in patients with proven fungal infections with a higher Cmax,ss/MIC associated with the desired clinical response.
Using the population pharmacokinetic data and the preliminary pharmacodynamic data that were derived from this study, it is possible to evaluate the suitability of selected L-AmB dose regimens in typical patients by using a simulation approach. Based on the predicted Cmax,ss, the PK-PD parameter (Cmax,ss/MIC) was generated for typical patients over the weight range of 10 to 70 kg who were treated with L-AmB at doses ranging from 1 mg/kg to 12.5 mg/kg by assuming that MIC for Candida spp. and Aspergillus spp. was 0.5 and 1.0 mg/liter, respectively. Using the a priori criterion of achieving no less than 85% of patients attaining the proposed pharmacodynamic target (Cmax,ss/MIC > 50) following L-AmB antifungal therapy, it was found that patients weighing greater than 30 kg attained the target Cmax,ss/MIC against Candida and Aspergillus when receiving an L-AmB daily dose of 4 mg/kg and 7.5 mg/kg, respectively, whereas lighter patients should be given even higher doses. However, further clinical studies that are sufficiently powered and appropriately designed are warranted to define the PD target of L-AmB antifungal efficacy, which can be used to inform the dosing strategy optimization in conjunction with the validated population pharmacokinetic model. It should be noted also that these metrics have utilized plasma concentrations of L-AmB which combine unbound AmB, plasma protein-bound AmB, and liposomal encapsulated AmB. The implications of using the unbound concentration of AmB in these calculations have yet to be elucidated.
In conclusion, a population pharmacokinetic model was developed and evaluated for L-AmB in pediatric oncology patients. The clearance and volume of distribution were influenced by body weight and showed significant between-occasion variability. The application of the population pharmacokinetic model in designing the optimal dosing strategy is worthy of further investigation in an attempt to improve the antifungal management in pediatric patients with malignant diseases.
Present address: Department of Pharmaceutical Sciences, University at Buffalo, SUNY. ![]()
|
|
|---|
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»