| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Cancer Therapy: Clinical |
Authors' Affiliations: Departments of 1 Oncology and 2 Pediatrics, Johns Hopkins University and Department of Pharmacy, The Johns Hopkins Hospital, Baltimore, Maryland; 3 The Cancer Therapy and Research Center, Institute for Drug Development, San Antonio, Texas; 4 The University of Texas M.D. Anderson Cancer Center, Houston, Texas; 5 Princess Margaret Hospital, Toronto, Ontario, Canada; 6 Pediatric Pharmacology Research Unit, University of California at San Diego, La Jolla, California; and 7 BioChem Pharma, Inc., Laval, Quebec, Canada
Requests for reprints: Sharyn D. Baker, The Sidney Kimmel Comprehensive Cancer Center at Johns Hopkins, 1650 Orleans Street, Room 1M86, Baltimore, MD 21231. Phone: 410-502-7149; Fax: 410-614-9006; E-mail: sbaker7{at}jhmi.edu.
| Abstract |
|---|
|
|
|---|
Experimental Design: Plasma samples from 111 cancer patients receiving troxacitabine (0.12-12.5 mg/m2) as a 30-minute infusion in phase I trials were used to develop the model with NONMEM. Clinical covariates evaluated included creatinine clearance, body surface area, age, and sex. From the model, a troxacitabine dosage rate of 2.0 to 3.0 mg/m2/d was expected to achieve a target concentration of 0.1 µmol/L; plasma samples were obtained during the infusion from eight patients receiving troxacitabine as a 3-day infusion.
Results: Troxacitabine pharmacokinetics were characterized by a three-compartment linear model. The mean value for systemic clearance [interindividual variability (CV%)] from the covariate-free model was 9.1 L/h (28%). Creatinine clearance and body surface area accounted for 36% of intersubject variation in clearance. Troxacitabine 2.0 mg/m2/d (n = 3) and 3.0 mg/m2/d (n = 5) for 3 days produced mean ± SD end of infusion concentrations of 0.12 ± 0.03 and 0.15 ± 0.03 µmol/L, respectively.
Conclusions: Renal function and body surface area were identified as sources of troxacitabine pharmacokinetic variability. The population pharmacokinetic model modelderived dosage rates for continuous infusion administration successfully achieved predetermined target plasma concentrations. The present model may be used to optimize treatment with troxacitabine by developing a dosing strategy based on both renal function and body size.
|
In phase I studies, the pharmacokinetic behavior of troxacitabine was characterized by a mean plasma clearance of 10 L/h, with
3-fold variability between patients, and renal excretion resulting in the elimination of 61% to 77% of the given dose as unchanged drug (1517). Plasma exposure toxicity assessments in patients with refractory leukemia receiving 30-minute infusions showed that elevated troxacitabine plasma area under the concentration-time curve was related to the severity of mucositis and hand-foot syndrome (13).
Recent preclinical data suggest that troxacitabine given as a continuous infusion may be more efficacious than intermittent infusions. In a HT-29 human xenograft model, prolonged exposure to low micromolar concentrations (
0.1 µmol/L) of troxacitabine for 3 to 5 days leads to superior growth inhibition compared with single high-dose bolus infusions (18). Based on this data, troxacitabine given by continuous infusion is being evaluated in patients with refractory acute myeloid leukemia (19). The starting dosage rate was 10.1 mg/m2/d, which was determined using the present population pharmacokinetic model to achieve a minimum target plasma concentration of 0.4 µmol/L. The initial infusion duration was 3 days with an increase in infusion to 5 to 6 days. In the first 13 patients treated, two complete responses were observed. With these promising results, accelerated development of continuous infusion troxacitabine is planned. A parallel phase I study of continuous infusion troxacitabine is being done in patients with advanced solid malignancies. In this study, the starting dosage rate was 2.0 to 3.0 mg/m2/d for 3 days, which was expected to reach a lower target plasma concentration of 0.1 µmol/L (20).
Population pharmacokinetic analysis is a useful tool for identification of sources of pharmacokinetic variability during anticancer drug development and can aid in the design of alternative dosing regimens to enhance their efficacy and safety (21). The objective of this study was to develop and validate a population pharmacokinetic model for troxacitabine to characterize clinical covariates that influence troxacitabine pharmacokinetics and to design a continuous dosing regimen to achieve a target plasma concentration of 0.1 µmol/L for phase I evaluation in patients with advanced solid malignancies.
| Materials and Methods |
|---|
|
|
|---|
Four covariates were considered in the analysis, including continuous [age, body surface area (BSA), and creatinine clearance] and categorical (sex) variables. BSA was calculated using Mosteller's formula: BSA = [height (cm) x weight (kg) / 3,600]0.5 (22). Creatinine clearance was estimated using the Crockcroft and Gault method and expressed in units of mL/min; entrance criteria to the clinical protocols included an estimated creatinine clearance of
45 mL/min (23). Plasma sampling was extensive; a mean of 12.6 samples (range, 9-16 samples) per patient were obtained between 4 minutes to 168.5 hours after the start of infusion. Troxacitabine was quantitated in plasma using a method involving high-performance liquid chromatography with tandem mass spectrometry detection and calibration curves ranging from 0.6 to 99.9 ng/mL, as previously described (12, 15, 16), and all concentrations were measured in the same laboratory.
Pharmacokinetic model
The population pharmacokinetic model for troxacitabine was developed in two steps: (a) covariate-free (base) model development and (b) covariate model development. Pharmacokinetic models were fitted to the data using the computer software programs NONMEM version V, level 1.1 (GloboMax LLC, Hanover, MD) and PDx POP version 1.1 (Globomax LLC). Xpose 3.1/S-PLUS 6.0 (Insightful Corporation, Seattle, WA) was used for graphical diagnostics and covariate screen.
Step 1: covariate-free (base) model development. The troxacitabine population pharmacokinetic model was built using two drug input methods: (a) actual dose (mg/h; model A) and (b) dose normalized to BSA (mg/m2/h; model B). Model B was explored because the dose of troxacitabine given was normalized to a patient's individual BSA. This approach would provide an independent assessment of the contribution of BSA to interindividual variability on troxacitabine clearance when compared with the covariate-free model for model A. The best structural pharmacokinetic model for troxacitabine was defined by fitting the data to predefined pharmacokinetic models provided in the NONMEM library. Linear compartmental pharmacokinetic models were evaluated because linear pharmacokinetic characteristics were noted in previous phase I trials (12, 15, 16). The pharmacokinetic model incorporated unknown intersubject and intrasubject variability and defined the number and form of
, a measure of intersubject variability, and
ij, a measure of intrasubject or residual variability. Intersubject variability for each pharmacokinetic variable was modeled with the exponential function. Residual error models of the additive, proportional, exponential, and combination methods were evaluated for the best structural pharmacokinetic model. The computational, intense, first-order conditional estimation method with or without interaction (FOCE + INTER) was used (24, 25).
For model discrimination, the NONMEM objective function (2 times the log-likelihood function or 2 log-likelihood) was used to calculate the Aikake Information Criterion = (2 log-likelihood) + 2p, for each nonhierarchical model, where 2p is twice the number of pharmacokinetic variables appearing in the model. Diagnostic plots for observed versus predicted plasma concentrations, residuals versus predicted concentrations, and weighed residuals versus predicted concentration were also considered during model discrimination.
Step 2: covariate model development. Covariate models were developed separately for models A and B as described above. Each candidate covariate identified from the screen was evaluated for inclusion to the selected structural pharmacokinetic model in a stepwise fashion by adding covariates in the order of highest to lowest level of significant correlation. General additive models and tree models were used in this analysis as implemented in the Xpose software (26).
Continuous variables were evaluated with the following models on pharmacokinetic variables as described in the manual for the software program PDx Pop: linear additive effect, normalized by median weight, centered on median weight, and power or allometric. The order of selecting pharmacokinetic variables for covariate model building was systemic clearance, volumes of distribution, and then intercompartmental clearance constants. Covariates were included into the model provided that the run resulted in successful minimization and covariate run status, a minimum reduction in objection function value of 3.875 (
2 distribution for 1 degree of freedom with
= 0.05), and a reduction in the intersubject variance of the associated pharmacokinetic variable, and the %RSE (relative SE of estimation) associated with each variable estimate was <50%.
A stepwise backward deletion procedure was done for both full covariate models. Covariates remained in the final population pharmacokinetic model when the removal of the covariate resulted in a
10.828 increase in the objective function (
2 distribution for 1 degree of freedom with
= 0.001). The final covariate models were run with the FOCE + INTER method.
Statistical assessment of final population pharmacokinetic models
To determine the influence of covariates on model prediction, a predictive performance evaluation was done for pharmacokinetic variables derived from both group I and II final population models as described previously (27, 28). Briefly, prediction errors (bias) were determined for pharmacokinetic variables estimated by the covariate-free (PE-N) and the final covariate (PE) models. For all models, percentage prediction error (PEj or PE-Nj) was calculated as (TVPKj PKj) / TVPKj x 100, where TVPKj is the typical (population) value for the pharmacokinetic variable for patient j, and PKj is the true value for the pharmacokinetic variable for patient j. The individual calculated PEj and PE-Nj of the pharmacokinetic variable was compared to see if the final population pharmacokinetic model had improved on predicting the respective pharmacokinetic variable. An improvement was defined as having the absolute value of PEj to be less than the absolute value of PE-Nj.
Validation of final population pharmacokinetic models
Validation of the final population pharmacokinetic models was undertaken using internal cross-validation techniques as recommended by the guidance document from the U.S. Department of Health and Human Services Food and Drug Administration Center for Drug Evaluation and Research Center for Biologics Evaluation and Research, Guidance for industry: population pharmacokinetics, February 1999 (http://www.fda.gov/cder/guidance/1852fnl.pdf) and described previously by Vozeh et al. (29).
Internal cross-validation. Briefly, 20 different data sets were generated by randomly removing the observed plasma concentrations of five to six subjects in each data set from the entire data set of 111 subjects, where data from each subject was removed only once. Each of the 20 data sets were then fitted with the final population pharmacokinetic models to provide individual predicted plasma concentrations from which an individual's standardized mean prediction error was calculated as previously described (29). The mean and 95% confidence interval for standardized mean prediction error were assessed; the model was considered to provide unbiased estimates for predicting troxacitabine plasma concentrations if the 95% confidence interval included 0.
Application of population pharmacokinetic model for the design of a continuous infusion regimen. One objective of the present study was to develop a population pharmacokinetic model to design a continuous dosing regimen to sustain a troxacitabine plasma concentration of 0.1 µmol/L for 3 days. Based on the final population pharmacokinetic models, a continuous infusion of troxacitabine 2.0 to 3.0 mg/m2/d was estimated to achieve a minimum target plasma concentration of 0.1 µmol/L, which is being evaluated in patients with advanced solid tumors (20). In this study, plasma samples were obtained at 4 hours after the start of infusion then daily at 24-hour intervals, with a final sample immediately before termination of the infusion, to determine if the target concentration of 0.1 µmol/L was achieved. Plasma troxacitabine concentrations for this study were measured using the same assay and laboratory used for the phase I studies (12, 1517).
| Results |
|---|
|
|
|---|
|
|
|
25% of intersubject variance. The further addition of BSA as a covariate on clearance decreased the CV% from 21% to 18%, indicating that BSA accounted for 14% of interindividual variability on troxacitabine clearance. These results are consistent with those obtained from the base models, where interindividual variability on clearance was decreased by 14% when drug in put was based on drug normalized to BSA (model B) compared with drug input based on the actual dose (model A). The objection function value continued to decrease with the successive addition of BSA on V1, V2, and Q3, but interindividual variability for clearance and other variables were not reduced (Fig. 3A
). The stepwise backward deletion process resulted in the removal of no covariates from the covariate model.
|
|
During the covariate screen for model B, where drug input was dose normalized to BSA in units of mg/m2/h, creatinine clearance was identified as a covariate for troxacitabine clearance. The addition of creatinine clearance on troxacitabine clearance decreased interindividual variability (CV%) for this variable from 24% to 20% (Fig. 2B). Subsequent models included the addition of AGE on V1, creatinine clearance on Q2, and BSA on Q3, with a significant reduction in the objective function, but provided no further explanation of interindividual variability on troxacitabine clearance (Fig. 2B). In addition, when AGE was added as a covariate on V1, the %RSE for V1 increased from 10% to 55% and increasingly worsened with addition of covariates on other variables. Likewise, with the addition of creatinine clearance on Q2, the %RSE for Q2 increased from 13% to 42%, and the %RSE for interindividual variability on this variable increased from 20% to 250%. Similarly, the addition of BSA to Q3 resulted in worsened %RSE (>100%) for this variable.
Table 3 presents troxacitabine typical population pharmacokinetic variables for the base model and final model for model B. From base to final model, the value for objection function value decreased by 32 (from 8,104 to 8,072), and the addition of creatinine clearance as a covariate accounted for 18% of unexplained variation on clearance (CV% decreased from 24% to 20%). The goodness of fit between observed and final model-predicted troxacitabine concentrations are shown in Fig. 3B (R2 = 0.8823).
Statistical assessment of final population pharmacokinetic models
Table 4
summarizes the prediction errors associated with pharmacokinetic variables estimated from the covariate-free (PE-N) and two final covariate (PE) population models, and the percentage of patients whose pharmacokinetic variable value predictions were improved by the covariate models. Addition of covariates reduced the prediction error and the SD for all pharmacokinetic variables and similar prediction errors were noted for both covariate models. Predictions were improved for clearance in >50% of patients for both models A and B.
|
Application of population pharmacokinetic model for the design of a continuous infusion regimen. Eight patients received troxacitabine 2.0 mg/m2/d (n = 3) for 3 days or 3.0 mg/m2/d for 3 days (n = 5). At 2.0 mg/m2/d, observed mean ± SD concentrations at 4, 24, 48, and 72 hours (end of infusion) were 0.06 ± 0.002, 0.10 ± 0.006, 0.10 ± 0.01, and 0.12 ± 0.03 µmol/L, respectively; at 3.0 mg/m2/d, concentrations were 0.08 ± 0.01, 0.13 ± 0.02, 0.15 ± 0.04, and 0.15 ± 0.03 µmol/L, respectively. Assessment of end of infusion concentrations show that the minimum target concentration of 0.1 µmol/L was achieved in all patients with the starting dose of 2.0 to 3.0 mg/m2/d, which was derived from the described population pharmacokinetic model.
| Discussion |
|---|
|
|
|---|
Troxacitabine is eliminated principally via renal excretion. During model development, it was anticipated that creatinine clearance would be a covariate explaining a significant percentage of interindividual variability in drug clearance. The addition of BSA as an independent, significant covariate on clearance is generally consistent with the correlation of body size with glomerular filtration rate (30) and supports the calculation of troxacitabine dose based on BSA. During covariate screening, age was not selected as a significant covariate on clearance. Because age was included in the formula for calculated creatinine clearance, inclusion of this variable may have accounted for age-related decline in renal function. Perhaps the inclusion of younger pediatric and older geriatric patients into the population may have resulted in retaining age into the final model. Sex was not included into the model most likely due to the predominance of male subjects (80%).
Both final models, one where drug input was the actual dose (model A) and one where drug input was the dose normalized to BSA (model B), equally accounted for a significant percentage of interindividual variability for troxacitabine clearance. Final estimates for interindividual variability were 18% and 20% for models A and B, respectively. SEs (%RSE) on variable estimates and predictive performance were similar for models A and B.
Recent preclinical data suggest that prolonged troxacitabine exposure at low micromolar concentrations (
0.1 µmol/L) was superior to bolus infusion administration, and that a 3- to 5-day infusion achieved optimal efficacy (18). Troxacitabine penetrates into cells via passive diffusion, which may explain, in part, why troxacitabine showed enhanced cytotoxicity with prolonged versus bolus administration. To optimize the therapeutic potential of troxacitabine, several phase I trials were designed to evaluate prolonged infusion schedules (19, 20). The population pharmacokinetic model was used to derive a starting dosage rate to achieve a target troxacitabine plasma concentration of 0.1 µmol/L in patients with advanced solid malignancies, and the target plasma concentration was achieved in all patients receiving 2.0 to 3.0 mg/m2/d, albeit with accompanying intersubject pharmacokinetic variability. The equations for troxacitabine clearance as described in Tables 2 and 3 (final models), could be used to individualize the dosage rate for continuous infusion administration. For example, an individualized dose could be determined using values for creatinine clearance and BSA and the clearance equation in Table 2 according to the following equation: dose (µg/h) = [(4.8 x creatinine clearance0.431) + (BSA 1.95) x 2.84] x target Css (µg/L). If the target Css is 0.1 µmol/L (21 µg/L), the individualized dosage rate for an individual with a BSA of 1.80 m2 and creatinine clearance of 45 mL/min (2.7 L/h) would be 146 µg/h (3.5 mg/d); the individualized dosage rate for an individual with the same BSA but a creatinine clearance of 100 mL/min (6.0 L/h) would be 209 µg/h (5.0 mg/d). If both patients received the same dose normalized to their BSA (e.g., 2.0 mg/m2/d or 3.6 mg/d), the patient with the higher creatinine clearance may have been undertreated. The clinical applicability and safety of the proposed dosing strategy would require prospective evaluation compared with the currently used BSA-based dosing.
In conclusion, a population pharmacokinetic model has been developed that accounted for 36% of unexplained intersubject variance for troxacitabine clearance using measures of renal function (calculated creatinine clearance) and body size (BSA). Further refinement of dosing strategies to decrease pharmacokinetic variability will likely include more accurate assessment of glomerular filtration rate and renal function. In addition, the population pharmacokinetic model was used to determine a dosage rate for continuous infusion administration that achieved predetermined target plasma concentrations. The present analysis shows that population pharmacokinetic modeling can be incorporated in the drug development process to identify sources of pharmacokinetic variability, aid in the design of alternative dosing regimens, and optimize therapy.
| Footnotes |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Conflict of interest statement: S.D. Baker and M. Hidalgo received research funding from Shire Pharmaceuticals.
Received 10/17/05; revised 12/15/05; accepted 1/17/06.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. Jimeno, W. A. Messersmith, C. K. Lee, W. W. Ma, D. Laheru, R. C. Donehower, S. D. Baker, and M. Hidalgo Phase I study of troxacitabine administered by continuous infusion in subjects with advanced solid malignancies Ann. Onc., February 1, 2008; 19(2): 374 - 379. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Sparreboom, A. C. Wolff, R. H.J. Mathijssen, E. Chatelut, E. K. Rowinsky, J. Verweij, and S. D. Baker Evaluation of Alternate Size Descriptors for Dose Calculation of Anticancer Drugs in the Obese J. Clin. Oncol., October 20, 2007; 25(30): 4707 - 4713. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |