A limited sampling schedule to estimate individual pharmacokinetic parameters of fludarabine in hematopoietic cell transplant patients

Prospective Validation and Refinement of a Population Pharmacokinetic Model of Fludarabine in Children and Young Adults Undergoing Hematopoietic Cell Transplantation

Fludarabine is a nucleoside analog with antileukemic and immunosuppressive activity commonly used in allogeneic hematopoietic cell transplantation (HCT). Several fludarabine population pharmacokinetic (popPK) and pharmacodynamic models have been published enabling the movement towards precision dosing of fludarabine in pediatric HCT; however, developed models have not been validated in a prospective cohort of patients. In this multicenter pharmacokinetic study, fludarabine plasma concentrations were collected via a sparse-sampling strategy. A fludarabine popPK model was evaluated and refined using standard nonlinear mixed effects modelling techniques. The previously described fludarabine popPK model well-predicted the prospective fludarabine plasma concentrations. Individuals who received model-based dosing (MBD) of fludarabine achieved significantly more precise overall exposure of fludarabine. The fludarabine popPK model was further improved by both the inclusion of fat-free mass instead of total body weight and a maturation function on fludarabine clearance. The refined popPK model is expected to improve dosing recommendations for children younger than 2 years and patients with higher body mass index. Given the consistency of fludarabine clearance and exposure across its multiple days of administration, therapeutic drug monitoring is not likely to improve targeted exposure attainment.

Population Pharmacokinetics of Fludarabine in Children and Adults during Conditioning Prior to Allogeneic Hematopoietic Cell Transplantation

Background

Fludarabine is often used as an important drug in reduced toxicity conditioning regimens prior to hematopoietic cell transplantation (HCT). As no definitive pharmacokinetic (PK) basis for HCT dosing for the wide age and weight range in HCT is available, linear body surface area (BSA)-based dosing is still used.

Objective

We sought to describe the population PK of fludarabine in HCT recipients of all ages.

Methods

From 258 HCT recipients aged 0.3–74 years, 2605 samples were acquired on days 1 (42%), 2 (17%), 3 (4%) and 4 (37%) of conditioning. Herein, the circulating metabolite of fludarabine was quantified, and derived concentration-time data were used to build a population PK model using non-linear mixed-effects modelling.

Results

Variability was extensive where area under the curve ranged from 10 to 66 mg h/L. A three-compartment model with first-order kinetics best described the data. Actual body weight (BW) with standard allometric scaling was found to be the best body-size descriptor for all PK parameters. Estimated glomerular filtration rate (eGFR) was included as a descriptor of renal function. Thus, clearance was differentiated into a non-renal (3.24 ± 20% L/h/70 kg) and renal (eGFR × 0.782 ± 11% L/h/70 kg) component. The typical volumes of distribution of the central (V1), peripheral (V2), and second peripheral (V3) compartments were 39 ± 8%, 20 ± 11%, and 50 ± 9% L/70 kg respectively. Intercompartmental clearances between V1 and V2, and V1 and V3, were 8.6 ± 8% and 3.8 ± 13% L/h/70 kg, respectively.

Conclusion

BW and eGFR are important predictors of fludarabine PK. Therefore, current linear BSA-based dosing leads to highly variable exposure, which may lead to variable treatment outcome.

Pharmacokinetics and Pharmacodynamics of Treosulfan in Patients With Thalassemia Major Undergoing Allogeneic Hematopoietic Stem Cell Transplantation

A treosulfan (Treo)-based conditioning regimen prior to hematopoietic stem cell transplantation (HSCT) has been successfully used in treating hematological malignant and nonmalignant diseases. We report Treo pharmacokinetics (PK) in patients with thalassemia major undergoing HSCT (n = 87), receiving Treo at a dose of 14 g/m² /day. Median Treo AUC and clearance (CL) was 1,326 mg*h/L and 10.8 L/h/m² , respectively. There was wide interindividual variability in Treo AUC and CL (64 and 68%) which was not explained by any of the variables tested. None of the Treo PK parameters were significantly associated with graft rejection or toxicity; however, Treo CL <7.97 L/h/m² was significantly associated with poor overall (hazard ratio (HR) 2.7, confidence interval (CI) (1.09-6.76), P = 0.032) and event-free survival (HR 2.4, CI (0.98-5.73), P = 0.055). Further studies in a larger cohort are warranted to identify the factors explaining the variation in Treo PK as well as to establish a therapeutic range of Treo for targeted dose adjustment to improve HSCT outcome.

Pharmacokinetics and Model-Based Dosing to Optimize Fludarabine Therapy in Pediatric Hematopoietic Cell Transplant Recipients

A prospective multicenter study was conducted to characterize the pharmacokinetics (PK) and pharmacodynamics (PD) of fludarabine plasma (f-ara-a) and intracellular triphosphate (f-ara-ATP) in children undergoing hematopoietic cell transplantation (HCT) and receiving fludarabine with conditioning. Plasma and peripheral blood mononuclear cells (PBMCs) were collected over the course of therapy for quantitation of f-ara-a and f-ara-ATP. Nonlinear mixed-effects modeling was used to develop the PK model, including identification of covariates impacting drug disposition. Data from a total of 133 children (median age, 5 years; range, .2 to 17.9) undergoing HCT for a variety of malignant and nonmalignant disorders were available for PK-PD modeling. The implementation of allometric scaling of PK parameters alone was insufficient to describe drug clearance, particularly in very young children. Renal impairment was predicted to increase drug exposure across all ages. The rate of f-ara-a entry into PBMCs (expressed in pmoles per million cells) decreased over the course of therapy, resulting in 78% lower f-ara-ATP after the fourth dose (1.7 pmoles/million cells [range, .2 to 7.2]) compared with first dose (7.9 pmoles/million cells [range, .7 to 18.2]). The overall incidence of treatment-related mortality (TRM) was low at 3% and 8% at days 60 and 360, respectively, and no association with f-ara-a exposure and TRM was found. In the setting of malignancy, disease-free survival was highest at 1 year after HCT in subjects achieving a systemic f-ara-a cumulative area under the curve (cAUC) greater than 15 mg*hour/L compared to patients with a cAUC less than 15 mg*hour/L (82.6% versus 52.8% P = .04). These results suggest that individualized model-based dosing of fludarabine in infants and young children may reduce morbidity and mortality through improved rates of disease-free survival and limiting drug-related toxicity. ClinicalTrials.gov Identifier: NCT01316549.

Neonatal Maturation of Paracetamol (Acetaminophen) Glucuronidation, Sulfation, and Oxidation Based on a Parent-Metabolite Population Pharmacokinetic Model

OBJECTIVES

This study aimed to model the population pharmacokinetics of intravenous paracetamol and its major metabolites in neonates and to identify influential patient characteristics, especially those affecting the formation clearance (CL_formation) of oxidative pathway metabolites.

METHODS

Neonates with a clinical indication for intravenous analgesia received five 15-mg/kg doses of paracetamol at 12-h intervals (<28 weeks' gestation) or seven 15-mg/kg doses at 8-h intervals (≥28 weeks' gestation). Plasma and urine were sampled throughout the 72-h study period. Concentration-time data for paracetamol, paracetamol-glucuronide, paracetamol-sulfate, and the combined oxidative pathway metabolites (paracetamol-cysteine and paracetamol-N-acetylcysteine) were simultaneously modeled in NONMEM 7.2.

RESULTS

The model incorporated 259 plasma and 350 urine samples from 35 neonates with a mean gestational age of 33.6 weeks (standard deviation 6.6). CL_formation for all metabolites increased with weight; CL_formation for glucuronidation and oxidation also increased with postnatal age. At the mean weight (2.3 kg) and postnatal age (7.5 days), CL_formation estimates (bootstrap 95% confidence interval; between-subject variability) were 0.049 L/h (0.038-0.062; 62 %) for glucuronidation, 0.21 L/h (0.17-0.24; 33 %) for sulfation, and 0.058 L/h (0.044-0.078; 72 %) for oxidation. Expression of individual oxidation CL_formation as a fraction of total individual paracetamol clearance showed that, on average, fractional oxidation CL_formation increased <15 % when plotted against weight or postnatal age.

CONCLUSIONS

The parent-metabolite model successfully characterized the pharmacokinetics of intravenous paracetamol and its metabolites in neonates. Maturational changes in the fraction of paracetamol undergoing oxidation were small relative to between-subject variability.

Population pharmacokinetics of fludarabine in patients with aplastic anemia and Fanconi anemia undergoing allogeneic hematopoietic stem cell transplantation

Although hematopoietic stem cell transplantation (HSCT) with a conditioning regimen consisting of fludarabine (F-araA) and cyclophosphamide (Cy) is associated with improved outcome in young patients with aplastic anemia (AA) and Fanconi anemia (FA), several factors limit the success of the procedure. We evaluated the population pharmacokinetics (POPPK) of F-araA and its influence on HSCT outcome in patients (n=53) with AA and FA undergoing HSCT. Patients carrying a 5′-UTR polymorphism in NT5E gene (rs2295890 G>C) exhibited significantly lower plasma F-araA clearance compared to those with wild-type genotype (7.12 vs 5.03 L/h/m² (29%) P<0.05). F-araA clearance was significantly higher in patients with AA compared to FA (2.46 ×, P<1e−6). Of all the outcome parameters evaluated (engraftment, rejection/graft failure, GvHD, TRM, OS), high F-araA AUC (>29.4 μm*h) was the only significant factor associated with the development of aGvHD by both univariate and multivariate analysis (P=0.02). The influence of plasma F-araA levels need to be evaluated in a larger cohort of patients to propose the need for therapeutic drug monitoring.

Simultaneous quantification of busulfan, clofarabine and F-ARA-A using isotope labelled standards and standard addition in plasma by LC-MS/MS for exposure monitoring in hematopoietic cell transplantation conditioning

In allogeneic hematopoietic cell transplantation (HCT) it has been shown that over- or underexposure to conditioning agents have an impact on patient outcomes. Conditioning regimens combining busulfan (Bu) and fludarabine (Flu) with or without clofarabine (Clo) are gaining interest worldwide in HCT. To evaluate and possibly adjust full conditioning exposure a simultaneous analysis of Bu, F-ARA-A (active metabolite of Flu) and Clo in one analytical run would be of great interest. However, this is a chromatographical challenge due to the large structural differences of Bu compared to F-ARA-A and Clo. Furthermore, for the bioanalysis of drugs it is common to use stable isotope labelled standards (SILS). However, when SILS are unavailable (in case of Clo and F-ARA-A) or very expensive, standard addition may serve as an alternative to correct for recovery and matrix effects. This study describes a fast analytical method for the simultaneous analysing of Bu, Clo and F-ARA-A with liquid chromatography-tandem mass spectrometry (LC-MS/MS) including standard addition methodology using 604 spiked samples. First, the analytical method was validated in accordance with European Medicines Agency guidelines. The lower limits of quantification (LLOQ) were for Bu 10μg/L and for Clo and F-ARA-A 1μg/L, respectively. Variation coefficients of LLOQ were within 20% and for low medium and high controls were all within 15%. Comparison of Bu, Clo and F-ARA-A standard addition results correspond with those obtained with calibration standards in calf serum. In addition for Bu, results obtained by this study were compared with historical data analysed within TDM. In conclusion, an efficient method for the simultaneous quantification of Bu, Clo and F-ARA-A in plasma was developed. In addition, a robust and cost-effective method to correct for matrix interference by standard addition was established.

Personalized fludarabine dosing to reduce nonrelapse mortality in hematopoietic stem-cell transplant recipients receiving reduced intensity conditioning

Patients undergoing hematopoietic cell transplantation (HCT) with reduced intensity conditioning (RIC) commonly receive fludarabine. Higher exposure of F-ara-A, the active component of fludarabine, has been associated with a greater risk of nonrelapse mortality (NRM). We sought to develop a model for fludarabine dosing in adult HCT recipients that would allow for precise dose targeting and predict adverse clinical outcomes. We developed a pharmacokinetic model from 87 adults undergoing allogeneic RIC HCT that predicts F-ara-A population clearance (Clpop) accounting for ideal body weight and renal function. We then applied the developed model to an independent cohort of 240 patients to identify whether model predictions were associated with NRM and acute graft versus host disease (GVHD). Renal mechanisms accounted for 35.6% of total F-ara-A Clpop. In the independent cohort, the hazard ratio of NRM at day 100 was significantly higher in patients with predicted F-ara-A clearance (Clpred) <8.50 L/h (P < 0.01) and area under the curve (AUCpred) >6.00 μg × h/mL (P = 0.01). A lower Clpred was also associated with more NRM at month 6 (P = 0.01) and trended toward significance at 12 months (P = 0.05). In multivariate analysis, low fludarabine clearance trended toward higher risk of acute GVHD (P = 0.05). We developed a model that predicts an individual's systemic F-ara-A exposure accounting for kidney function and weight. This model may provide guidance in dosing especially in overweight individuals and those with altered kidney function.

Association of fludarabine pharmacokinetic/dynamic biomarkers with donor chimerism in nonmyeloablative HCT recipients

PURPOSE

Fludarabine monophosphate (fludarabine) is an integral component of many reduced-intensity conditioning regimens for hematopoietic cell transplantation (HCT). Fludarabine's metabolite, 9-β-D-arabinofuranosyl-2-fluoroadenine (F-ara-A), undergoes cellular uptake and activation to form the active cytotoxic metabolite fludarabine triphosphate (F-ara-ATP), which inhibits cellular DNA synthesis in CD4(+) and CD8(+) cells. In this study, we evaluated whether fludarabine-based pharmacologic biomarkers were associated with clinical outcomes in HCT recipients.

METHODS

Participants with hematologic diseases were conditioned with fludarabine and low-dose total body irradiation (TBI) followed by allogeneic HCT and post-grafting immunosuppression. After fludarabine administration, we evaluated pharmacological biomarkers for fludarabine-F-ara-A area under the curve (AUC) and the ratio of circulating CD4(+) and CD8(+) cells (CD4(+)/CD8(+) ratio) after fludarabine administration-in 102 patients; F-ara-ATP accumulation rate in enriched CD4(+) and CD8(+) cells was evaluated in 36 and 34 patients, respectively.

RESULTS

Interpatient variability in the pharmacological biomarkers was high, ranging from 3.7-fold (F-ara-A AUC) to 39-fold (F-ara-ATP in CD8(+) cells). Circulating CD8(+) cells were more sensitive to fludarabine administration. A population pharmacokinetic-based sampling schedule successfully allowed for estimation of F-ara-A AUC in this outpatient population. There was a poor correlation between the F-ara-AUC and the F-ara-ATP accumulation rate in CD4(+) (R (2) = 0.01) and CD8(+) cells (R (2) = 0.00). No associations were seen between the four biomarkers and clinical outcomes (day +28 donor T cell chimerism, acute graft-versus-host disease (GVHD), neutrophil nadirs, cytomegalovirus reactivation, chronic GVHD, relapse, non-relapse mortality, or overall mortality).

CONCLUSIONS

Considerable interpatient variability exists in pharmacokinetic and fludarabine-based biomarkers, but these biomarkers are not associated with clinical outcomes in fludarabine/TBI-conditioned patients.

Population pharmacokinetic/dynamic model of lymphosuppression after fludarabine administration

PURPOSE

Quantitative relationships between 9-β-D-arabinofuranosyl-2-fluoroadenine (F-ara-A) concentrations and lymphosuppression have not been reported, but would be useful for regimen design. A population pharmacokinetic/pharmacodynamic model was constructed in this study using data from 41 hematopoietic cell transplant (HCT) recipients conditioned with busulfan in combination with fludarabine (total dose 120 mg/m², Protocol 1519) or with fludarabine (total dose 250 mg/m²) with rabbit antithymocyte globulin (rATG, Protocol 2041).

METHODS

Individual pharmacokinetic parameters were fixed to post hoc Bayesian estimates, and circulating absolute lymphocyte counts (ALC) were obtained during the 3 weeks prior to graft infusion. A semi-physiological cell-kill model with three lymphocyte transit compartments was applied and aptly characterized the time course of suppression of circulating ALC by fludarabine administration. Drug- and system-specific parameters were estimated using a maximum likelihood expectation maximization algorithm, and the final model was qualified using an internal visual predictive check.

RESULTS

The final model successfully characterized the time course and variability in ALC. Pharmacodynamic parameters exhibited considerable between subject variability (38.9-211 %). The HCT protocol was the only covariate associated with the pharmacodynamic parameters, specifically the lymphocyte kill rate, the transit rate between lymphocyte compartments, and the baseline ALC.

CONCLUSIONS

This model can be used to simulate the degree of lymphosuppression for design of future fludarabine-based conditioning regimens.

Optimizing drug therapy in pediatric SCT: focus on pharmacokinetics

Given age-related differences in drug metabolism and indications for hematopoietic SCT (HSCT), personalized drug dosing of the conditioning regimen and post-transplant immunosuppression may reduce graft rejection, relapse rates and toxicity in pediatric HSCT recipients. This manuscript summarizes the pharmacokinetic/dynamic data of HSCT conditioning and post-grafting immunosuppression, presented at the First Annual Pediatric Bone Marrow Transplant Consortium (PBMTC) meeting in April 2013. Personalized dosing of BU to a target plasma exposure reduces graft rejection in children and improves relapse/toxicity rates in adults. Current weight-based dosing achieves the target BU exposure in only a minority (24.3%) of children. The initial BU dose should be based on the European Medicines Agency nomogram or population pharmacokinetic models to improve the numbers of children achieving the target exposure. There are limited pharmacokinetic data for treosulfan, CY, fludarabine and alemtuzumab as HSCT conditioning in children. For post-grafting immunosuppression, mycophenolic acid (MPA) clearance may be increased in younger children (<12 years). The preferred MPA pharmacokinetic monitoring parameters and target range are still evolving in HSCT recipients. Multi-institutional trials incorporating properly powered pharmacokinetic/dynamic studies are needed to assess the effect of variability in the plasma exposure of drugs/metabolites on clinical outcomes in pediatric HSCT recipients.

Recipient pretransplant inosine monophosphate dehydrogenase activity in nonmyeloablative hematopoietic cell transplantation

Mycophenolic acid, the active metabolite of mycophenolate mofetil (MMF), inhibits inosine monophosphate dehydrogenase (IMPDH) activity. IMPDH is the rate-limiting enzyme involved in de novo synthesis of guanosine nucleotides and catalyzes the oxidation of inosine 5'-monophosphate to xanthosine 5'-monophosphate (XMP). We developed a highly sensitive liquid chromatography-mass spectrometry method to quantitate XMP concentrations in peripheral blood mononuclear cells (PMNCs) isolated from the recipient pretransplant and used this method to determine IMPDH activity in 86 nonmyeloablative allogeneic hematopoietic cell transplantation (HCT) patients. The incubation procedure and analytical method yielded acceptable within-sample and within-individual variability. Considerable between-individual variability was observed (12.2-fold). Low recipient pretransplant IMPDH activity was associated with increased day +28 donor T cell chimerism, more acute graft-versus-host disease (GVHD), lower neutrophil nadirs, and more cytomegalovirus reactivation but not with chronic GVHD, relapse, nonrelapse mortality, or overall mortality. We conclude that quantitation of the recipient's pretransplant IMPDH activity in PMNC lysate could provide a useful biomarker to evaluate a recipient's sensitivity to MMF. Further trials should be conducted to confirm our findings and to optimize postgrafting immunosuppression in nonmyeloablative HCT recipients.

Pharmacokinetic/pharmacodynamic modelling approaches in paediatric infectious diseases and immunology

Pharmacokinetic/pharmacodynamic (PKPD) modelling is used to describe and quantify dose-concentration-effect relationships. Within paediatric studies in infectious diseases and immunology these methods are often applied to developing guidance on appropriate dosing. In this paper, an introduction to the field of PKPD modelling is given, followed by a review of the PKPD studies that have been undertaken in paediatric infectious diseases and immunology. The main focus is on identifying the methodological approaches used to define the PKPD relationship in these studies. The major findings were that most studies of infectious diseases have developed a PK model and then used simulations to define a dose recommendation based on a pre-defined PD target, which may have been defined in adults or in vitro. For immunological studies much of the modelling has focused on either PK or PD, and since multiple drugs are usually used, delineating the relative contributions of each is challenging. The use of dynamical modelling of in vitro antibacterial studies, and paediatric HIV mechanistic PD models linked with the PK of all drugs, are emerging methods that should enhance PKPD-based recommendations in the future.

A pilot pharmacologic biomarker study in HLA-haploidentical hematopoietic cell transplant recipients

PURPOSE

Eleven patients diagnosed with various hematologic malignancies receiving an HLA-haploidentical hematopoietic cell transplant (HCT) participated in an ancillary biomarker trial. The goal of the trial was to evaluate potential pharmacologic biomarkers pertinent to the conditioning regimen [fludarabine monophosphate (fludarabine) and cyclophosphamide (CY)] or postgrafting immunosuppression [CY and mycophenolate mofetil (MMF)] in these patients.

METHODS

We characterized the interpatient variability of nine pharmacologic biomarkers. The biomarkers evaluated were relevant to fludarabine (i.e., area under the curve (AUC) of 2-fluoro-ara-A or F-ara-A), CY (i.e., AUCs of CY and four of its metabolites), and MMF (i.e., total mycophenolic acid (MPA) AUC, unbound MPA AUC, and inosine monophosphate dehydrogenase (IMPDH) activity).

RESULTS

Interpatient variability in the pharmacologic biomarkers was high. Among those related to HCT conditioning, the interpatient variability ranged from 1.5-fold (CY AUC) to 4.0-fold (AUC of carboxyethylphosphoramide mustard, a metabolite of CY). Among biomarkers evaluated as part of postgrafting immunosuppression, the interpatient variability ranged from 1.7-fold (CY AUC) to 4.9-fold (IMPDH area under the effect curve). There was a moderate correlation (R (2) = 0.441) of within-patient 4-hydroxycyclophosphamide formation clearance.

CONCLUSIONS

Considerable interpatient variability exists in the pharmacokinetic and drug-specific biomarkers potentially relevant to clinical outcomes in HLA-haploidentical HCT recipients. Pharmacodynamic studies are warranted to optimize the conditioning regimen and postgrafting immunosuppression administered to HLA-haploidentical HCT recipients.

Improved quantitative method for fludarabine in human plasma by liquid chromatography and tandem mass spectrometry

An improved quantitative assay was developed and validated for fludarabine in human plasma. Fludarabine and its internal standard, cladribine, were separated on a C18 analytical column after sample purification by strong anion-exchange solid-phase extraction. Quantitation was performed by electrospray triple-quadrupole mass spectrometry in positive ionization mode using multiple-reaction monitoring. This assay had excellent inter- and intra-assay precisions within 8%, and accuracies ranging from 100 to 116%. The method was linear within the concentration range of 0.2-250ng/mL using 100μL of plasma with mean R(2)=0.9999. The extraction recoveries were 85% for fludarabine and 95% for the internal standard, which represent a significant improvement over the previously published methods. We utilized this method for pharmacokinetic (PK) investigations in 215 patients. Interference peaks constantly observed in each blank plasma sample were well resolved from fludarabine using our optimized LC-MS/MS conditions, demonstrating the reliability of this improved assay. The validated method will be further applied to PK studies within our bone marrow transplant program, which will allow for optimal dose and scheduling of fludarabine in these patients.

Population pharmacokinetics and dose optimization of mycophenolic acid in HCT recipients receiving oral mycophenolate mofetil

We sought to create a population pharmacokinetic model for total mycophenolic acid (MPA), to study the effects of different covariates on MPA pharmacokinetics, to create a limited sampling schedule (LSS) to characterize MPA exposure (i.e., area under the curve or AUC) with maximum a posteriori Bayesian estimation, and to simulate an optimized dosing scheme for allogeneic hematopoietic cell transplantation (HCT) recipients. Four thousand four hundred ninety-six MPA concentration-time points from 408 HCT recipients were analyzed retrospectively using a nonlinear mixed effects modeling approach. MPA pharmacokinetics was characterized with a two-compartment model with first-order elimination and a time-lagged first-order absorption process. Concomitant cyclosporine and serum albumin were significant covariates. The median MPA clearance (CL) and volume of the central compartment were 24.2 L/hour and 36.4 L, respectively, for a 70 kg patient receiving tacrolimus with a serum albumin of 3.4 g/dL. Dosing simulations indicated that higher oral MMF doses are needed with concomitant cyclosporine, which increases MPA CL by 33.8%. The optimal LSS was immediately before and at 0.25 hours, 1.25 hours, 2 hours, and 4 hours after oral mycophenolate mofetil administration. MPA AUC in an individual HCT recipient can be accurately estimated using a five-sample LSS and maximum a posteriori Bayesian estimation.

A limited sampling schedule to estimate mycophenolic Acid area under the concentration-time curve in hematopoietic cell transplantation recipients

Mycophenolate mofetil (MMF) is a key component of postgrafting immunosuppression in hematopoietic cell transplant (HCT) recipients. The plasma area under the curve (AUC) of its active metabolite, mycophenolic acid (MPA), is associated with MMF efficacy and toxicity. This study developed a population pharmacokinetic model of MPA in HCT recipients and created limited sampling schedules (LSSs) to enable individualized pharmacotherapy. A retrospective evaluation of MPA concentration-time data following a 2-hour MMF intravenous (IV) infusion was conducted in 77 HCT recipients. The final model consisted of 1 and 2 compartments for MMF and MPA pharmacokinetics, respectively. The mean estimated values (coefficient of variation, %) for total systemic clearance, distributional clearance, and central and peripheral compartment volumes of MPA were 36.9 L/h (34.5%), 15.3 L/h (80.4%), 11.9 L (71.7%), and 182 L (127%), respectively. No covariates significantly explained variability among individuals. Optimal LSSs were derived using a simulation approach based on the scaled mean squared error. A 5-sample schedule of 2, 2.5, 3, 5, and 6 hours from the start of the infusion precisely estimated MPA AUC(0-12 h) for Q12-hour IV MMF. A comparable schedule (2, 2.5, 3, 4, and 6 hours) similarly estimated MPA AUC(0-8) (h) for Q8-hour dosing.

Optimal sampling strategy development methodology using maximum a posteriori Bayesian estimation

Maximum a posteriori Bayesian (MAPB) pharmacokinetic parameter estimation is an accurate and flexible method of estimating individual pharmacokinetic parameters using individual blood concentrations and prior information. In the past decade, many studies have developed optimal sampling strategies to estimate pharmacokinetic parameters as accurately as possible using either multiple regression analysis or MAPB estimation. This has been done for many drugs, especially immunosuppressants and anticancer agents. Methods of development for optimal sampling strategies (OSS) are diverse and heterogeneous. This review provides a comprehensive overview of OSS development methodology using MAPB pharmacokinetic parameter estimation, determines the transferability of published OSSs, and compares sampling strategies determined by MAPB estimation and multiple regression analysis. OSS development has the following components: 1) prior distributions; 2) reference value determination; 3) optimal sampling time identification; and 4) validation of the OSS. Published OSSs often lack all data necessary for the OSS to be clinically transferable. MAPB estimation is similar to multiple regression analysis in terms of predictive performance but superior in flexibility.

A pilot pharmacologic biomarker study of busulfan and fludarabine in hematopoietic cell transplant recipients

PURPOSE

Sixteen patients diagnosed with various hematologic malignancies participated in a phase II study evaluating the addition of rabbit antithymocyte globulin (rATG, Thymoglobulin(®)) to the hematopoietic cell transplant (HCT) conditioning regimen of IV fludarabine monophosphate (fludarabine) and targeted intravenous (IV) busulfan (fludarabine/(T)busulfan). Our goal was to evaluate pharmacologic biomarkers pertinent to both medications in these patients.

METHODS

We characterized the interpatient variability of pharmacologic biomarkers relevant to busulfan, specifically busulfan concentration at steady state, and fludarabine, specifically F-ara-A area under the curve (AUC) and fludarabine triphosphate (F-ara-ATP) intracellular accumulation and concentration in separate CD4(+) and CD8(+) T-lymphocyte populations.

RESULTS

Acute and chronic graft versus host disease (GvHD) occurred in 11 patients and one patient, respectively. Four patients died before day +100 of non-relapse causes, which met the protocol stopping guidelines. The cumulative incidence of relapse was 25% at 3 year post-HCT. Interpatient variability in the busulfan- and fludarabine-relevant pharmacologic biomarkers was 2.1- to 2.5-fold. F-ara-A AUC and accumulated F-ara-ATP in CD8(+) cells had the highest hazard ratio for non-relapse mortality and overall survival, respectively. However, neither achieved statistical significance.

CONCLUSIONS

The low rates of GvHD, particularly in its chronic form, were encouraging, and further biomarker studies are warranted to optimize the fludarabine/(T)busulfan/rATG conditioning regimen.