Population pharmacokinetics of busulfan in children: increased evidence for body surface area and allometric body weight dosing of busulfan in children

Identification, Isolation, and Characterization of a Novel Degradation Impurity of Busulfan Using Preparative Chromatography, NMR, and LC-MS

BACKGROUND

Busulfan is the most effective medication for treating chronic myelogenous or granulocytic leukemia because it has cytotoxic properties that harm or kill hematopoietic cells. It cannot absorb light in the Ultraviolet range due to its structure. Because of this, it is very challenging to quantify using traditional HPLC coupled with UV/Photodiode Array detectors. So, using sodium diethyldithiocarbamate, a derivatization method was developed to quantify related impurities. A significant unknown impurity was identified in derivatized samples of busulfan and a noticeably high percentage level was discovered during routine drug testing.

OBJECTIVE

We aimed to isolate, and characterize the unknown impurity observed in the samples and to identify its root cause.

METHODS

Preparative HPLC was used to isolate the unidentified, derivatized impurity, and 1H NMR, 13C NMR, and MS were used to decipher its structural components.

RESULTS

The spectral characterization data analysis showed that the unknown impurity was related to busulfan. Additionally, it was noted that the impurity developed as a result of the residual buffer used to prepare the derivatizing reagent.

CONCLUSION

The isolated impurity was found to be same as comparable to that found in busulfan drug substances, according to the results of the characterization tools. An alternative method of reagent preparation was optimized and deemed satisfactory because the buffer used in reagent preparation is the only factor contributing to the formation of impurities.

HIGHLIGHTS

Using cutting-edge analytical characterization tools, it was possible to explain the structural characteristics of an unknown impurity and discover that it was a novel impurity, which undoubtedly contributes to the comprehension of drug substance reaction properties.

Population Pharmacokinetic Model of Intravenous Busulfan in Hematopoietic Cell Transplantation: Systematic Review and Comparative Simulations

BACKGROUND

Busulfan is commonly used in the chemotherapy prior to hematopoietic cell transplantation (HCT). Busulfan has a narrow therapeutic window and a well-established exposure-response relationship with important clinical outcomes. Model-informed precision dosing (MIPD) based on population pharmacokinetic (popPK) models has been implemented in the clinical settings. We aimed to systematically review existing literature on popPK models of intravenous busulfan.

METHODS

We systematically searched Ovid MEDLINE, EMBASE, Cochrane Library, Scopus, and Web of Science databases from inception to December 2022 to identify original popPK models (nonlinear mixed-effect modeling) of intravenous busulfan in HCT population. Model-predicted busulfan clearance (CL) was compared using US population data.

RESULTS

Of the 44 eligible popPK studies published since 2002, 68% were developed predominantly in children, 20% in adults, and 11% in both children and adults. The majority of the models were described using first-order elimination or time-varying CL (69% and 26%, respectively). All but three included a body-size descriptor (e.g., body weight, body surface area). Other commonly included covariates were age (30%) and GSTA1 variant (15%). Median between-subject and between-occasion variabilities of CL were 20% and 11%, respectively. Between-model variabilities in predicted median CL were < 20% in all of the weight tiers (10-110 kg) in the simulation based on US population data.

CONCLUSION

Busulfan PK is commonly described using a first-order elimination or time-varying CL. A simple model with limited covariates were generally sufficient to attain relatively small unexplained variabilities. However, therapeutic drug monitoring may still be necessary to attain a narrow target exposure.

Individualizing busulfan dose in specific populations and evaluating the risk of pharmacokinetic drug-drug interactions

INTRODUCTION

Busulfan is an alkylating agent widely used in the conditioning of hematopoietic stem cell transplantation possessing a complex metabolism and a large interindividual and intra-individual variability, especially in children. Combined with the strong rationale of busulfan PK/PD relationships, factors altering its clearance (e.g., weight, age, and GST-A genetic polymorphism mainly) can also affect clinical outcomes.

AREAS COVERED

This review aims to provide an overview of the current knowledge on busulfan pharmacokinetics, its pharmacokinetics variabilities in pediatric populations, drug-drug interactions (DDI), and their consequences regarding dose individualization. This review was based on medical literature up until October 2021.

EXPERT OPINION

To ensure effective busulfan exposure in pediatrics, different weight-based nomograms have been established to determine busulfan dosage and provided improved results (65 - 80% of patients correctly exposed). In addition to nomograms, therapeutic drug monitoring (TDM) of busulfan measuring plasmatic concentrations to estimate busulfan pharmacokinetic parameters can be used. TDM is now widely carried out in routine practices and aims to ensure the targeting of the reported therapeutic windows by individualizing busulfan dosing based on the clearance estimations from a previous dose.

Pre-existing anti-factor VIII immunity alters therapeutic platelet-targeted factor VIII engraftment following busulfan conditioning through cytotoxic CD8 T cells

BACKGROUND

We previously demonstrated that busulfan preconditioning enabled sustained therapeutic platelet-derived factor VIII (FVIII) expression in naïve FVIIInull mice transplanted with 2bF8-transduced Sca-1+ cells. However, in mice with pre-existing inhibitors, platelet-FVIII expression was lost.

OBJECTIVE

In this study, we aimed to describe the mechanism of this platelet-FVIII loss.

METHODS

We monitored platelet-FVIII expression in FVIIInull mice that were immunized with rhFVIII to induce inhibitors and subsequently conditioned with busulfan before whole bone marrow transplantation or Sca-1+ hematopoietic stem cell transplantation (HSCT) from 2bF8 transgenic (2bF8Tg) mice. Busulfan with or without antithymocyte globulin or anti-CD8 antibody was employed before 2bF8Tg HSCT. Interferon gamma-ELISpot assay was used to assess which subset of cells was the target in platelet-FVIII loss. B-cell-deficient homozygous mutant mice were used to determine whether platelet-FVIII loss in FVIII-primed mice was mediated by antibody-dependent cellular cytotoxicity.

RESULTS

Platelet-FVIII expression was sustained in 2bF8Tg bone marrow transplantation but not in 2bF8Tg HSCT recipients. CD8 T-cell depletion in addition to busulfan preconditioning restored platelet-FVIII expression in 2bF8Tg-HSCT recipients. ELISpot analyses showed that FVIII-primed CD8 T cells were efficiently restimulated by 2bF8Tg-Sca-1+ cells and secreted interferon gamma, but were not stimulated by 2bF8Tg platelets/megakaryocytes, suggesting that 2bF8Tg-Sca-1+ cells are targets for FVIII-primed CD8 T cells. When 2bF8Tg-Sca-1+ cells were transplanted into FVIII-primed homozygous mutant mice preconditioned with busulfan, no FVIII expression was detected, suggesting that antibody-dependent cellular cytotoxicity was not the mechanism of platelet-FVIII loss in FVIII-primed mice.

CONCLUSION

Pre-existng immunity can alter the engraftment of 2bF8Tg-Sca-1+ cells through the cytotoxic CD8 T-cell-mediated pathway. Sufficient eradication of FVIII-primed CD8 T cells is critical for the success of platelet gene therapy in hemophilia A with inhibitors.

Individualized busulfan dosing improves outcomes compared to fixed-dose administration in pre-transplant minimal residual disease-positive acute myeloid leukemia patients with intermediate-risk undergoing allogeneic stem cell transplantation in CR

Pre-transplant minimal residual disease (MRD) impacts negatively on post-transplant relapse risk in acute myeloid leukemia (AML). Therapeutic drug monitoring by calculating area under the curve (AUC) was developed to optimize busulfan (Bu) exposure. Here, we compared post-transplant outcomes after individualized versus fixed busulfan dosage in intermediate-risk AML who achieved CR prior to allograft focusing on pre-transplant flow-MRD. Eighty-seven patients (median, 56 years) with intermediate-risk AML and pre-transplant flow-MRD ("different from normal") were included. Thirty-two patients received individualized busulfan; 54 fixed dosages. Individualized dosage was adjusted in 25/32 patients: increased, n = 18/25 (72%); decreased: n = 7/25 (28%). After median follow-up of 27 months, we observed lower 3-year relapses (6%, 2%-19% vs. 35%, 23%-49% p = 0.02), improved 3-year leukemia-free survival (LFS) (78%, 54%-91% vs. 55%, 40%-70% p = 0.009) and - overall survival (OS) (82%, 60%-93% vs. 69%, 54%-81% p = 0.05) after individualized compared to fixed Bu. Non-relapsed mortality (NRM) and acute graft versus host disease (GvHD) were not different. In multivariate analysis, fixed Bu showed unfavorable impact on OS (hazard ratio [HR] 4.6, p = 0.044), LFS (HR 3.6, p = 0.018) and relapses (HR 3.6, p = 0.033). Fixed Bu also had unfavorable impact on LFS (3.6, 1.1-12.6, p = 0.041) in pre-transplant MRD-positive patients. Individualized, AUC-based, busulfan is associated with lower relapses in intermediate-risk AML patients allografted in CR and may overcome pre-transplant MRD-positivity.

Implementation and Cross-Validation of a Pharmacokinetic Model for Precision Dosing of Busulfan in Hematopoietic Stem Cell Transplanted Children

Busulfan, a drug used in conditioning prior to hematopoietic stem cell transplantation (HSCT) in children, has a narrow therapeutic margin. The model-informed precision dosing (MIPD) of busulfan is desirable, but there is a lack of validated tools. The objective of this study was to implement and cross-validate a population pharmacokinetic (PK) model in the Tucuxi software for busulfan MIPD in HSCT children. A search of the literature was performed to identify candidate population PK models. The goodness of fit of three selected models was assessed in a dataset of 178 children by computing the mean error (ME) and root-mean-squared error of prediction (RMSE). The best model was implemented in Tucuxi. The individual predicted concentrations, the area under the concentration-time curve (AUC), and dosage requirements were compared between the Tucuxi model and a reference model available in the BestDose software in a subset of 61 children. The model from Paci et al. best fitted the data in the full dataset. In a subset of 61 patients, the predictive performance of Tucuxi and BestDose models was comparable with ME values of 6.4% and -2.5% and RMSE values of 11.4% and 13.6%, respectively. The agreement between the estimated AUC and the predicted dose was good, with 6.6% and 4.9% of the values being out of the 95% limits of agreement, respectively. To conclude, a PK model for busulfan MIPD was cross-validated and is now available in the Tucuxi software.

External Evaluation of Population Pharmacokinetic Models of Busulfan in Chinese Adult Hematopoietic Stem Cell Transplantation Recipients

Objective: Busulfan (BU) is a bi-functional DNA-alkylating agent used in patients undergoing hematopoietic stem cell transplantation (HSCT). Over the last decades, several population pharmacokinetic (pop PK) models of BU have been established, but external evaluation has not been performed for almost all models. The purpose of the study was to evaluate the predictive performance of published pop PK models of intravenous BU in adults using an independent dataset from Chinese HSCT patients, and to identify the best model to guide personalized dosing.

Methods: The external evaluation methods included prediction-based diagnostics, simulation-based diagnostics, and Bayesian forecasting. In prediction-based diagnostics, the relative prediction error (PE%) was calculated by comparing the population predicted concentration (PRED) with the observations. Simulation-based diagnostics included the prediction- and variability-corrected visual predictive check (pvcVPC) and the normalized prediction distribution error (NPDE). Bayesian forecasting was executed by giving prior one to four observations. The factors influencing the model predictability, including the impact of structural models, were assessed.

Results: A total of 440 concentrations (110 patients) were obtained for analysis. Based on prediction-based diagnostics and Bayesian forecasting, preferable predictive performance was observed in the model developed by Huang et al. The median PE% was -1.44% which was closest to 0, and the maximum F₂₀ of 57.27% and F₃₀ of 72.73% were achieved. Bayesian forecasting demonstrated that prior concentrations remarkably improved the prediction precision and accuracy of all models, even with only one prior concentration.

Conclusion: This is the first study to comprehensively evaluate published pop PK models of BU. The model built by Huang et al. had satisfactory predictive performance, which can be used to guide individualized dosage adjustment of BU in Chinese patients.

The Correlation Between Busulfan Exposure and Clinical Outcomes in Chinese Pediatric Patients: A Population Pharmacokinetic Study

Aims: The aims of the study were to 1) establish a population pharmacokinetic (Pop-PK) model for busulfan in Chinese pediatric patients undergoing hematopoietic stem cell transplantation (HSCT) and then estimate busulfan exposure and 2) explore the association between busulfan exposure and clinical outcomes.

Methods: A total of 128 patients with 467 busulfan concentrations were obtained for Pop-PK modeling using nonlinear mixed effect model (NONMEM) software. Sixty-three patients who received the 16-dose busulfan conditioning regimen were enrolled to explore the correlations between clinical outcomes and the busulfan area under the concentration–time curve (AUC) using the Cox proportional hazards regression model, Kaplan–Meier method and logistic regression.

Results: The typical values for clearance (CL) and distribution volume (V) of busulfan were 7.71 L h⁻¹ and 42.4 L, respectively. The allometric normal fat mass (NFM) and maturation function (Fmat) can be used to describe the variability in CL, and the fat-free mass (FFM) can be used to describe the variability in V. Patients with AUCs of 950–1,600 µM × min had 83.7% (95% CI: 73.3–95.5) event-free survival (EFS) compared with 55.0% (95% CI: 37.0–81.8) for patients with low or high exposure (p = 0.024). The logistic regression analysis results showed no association between transplant-related toxicities and the busulfan AUC (p > 0.05).

Conclusions: The variability in busulfan CL was related to the NFM and Fmat, while busulfan V was related to the FFM. Preliminary analysis results suggested that a busulfan AUC of 950–1,600 µM × min was associated with better EFS in children receiving the 16-dose busulfan regimen.

Population Pharmacokinetics of Busulfan and Its Metabolite Sulfolane in Patients with Myelofibrosis Undergoing Hematopoietic Stem Cell Transplantation

For patients with myelofibrosis, allogeneic hematopoietic stem cell transplantation (allo-HSCT) remains the only curative treatment to date. Busulfan-based conditioning regimens are commonly used, although high inter-individual variability (IIV) in busulfan drug exposure makes individual dose selection challenging. Since data regarding the IIV in patients with myelofibrosis are sparse, this study aimed to develop a population pharmacokinetic (PopPK) model of busulfan and its metabolite sulfolane in patients with myelofibrosis. The influence of patient-specific covariates on the pharmacokinetics of drug and metabolite was assessed using non-linear mixed effects modeling in NONMEM^®. We obtained 523 plasma concentrations of busulfan and its metabolite sulfolane from 37 patients with myelofibrosis. The final model showed a population clearance (CL) and volume of distribution (V_d) of 0.217 L/h/kg and 0.82 L/kg for busulfan and 0.021 L/h/kg and 0.65 L/kg for its metabolite. Total body weight (TBW) and a single-nucleotide polymorphism of glutathione-S-transferase A1 (GSTA1 SNP) displayed a significant impact on volume of distribution and metabolite clearance, respectively. This is the first PopPK-model developed to describe busulfan’s pharmacokinetics in patients with myelofibrosis. Incorporating its metabolite sulfolane into the model not only allowed the characterization of the covariate relationship between GSTA1 and the clearance of the metabolite but also improved the understanding of busulfan’s metabolic pathway.

Clinical outcomes of individualized busulfan-dosing in hematopoietic stem cell transplantation in Chinese children undergoing with therapeutic drug monitoring

To identify relationships between busulfan (Bu) exposure and outcomes of a cohort pediatric patients receiving hematopoietic stem cell transplantation (HSCT), along with a targeted busulfan-based conditioning regimen. We retrospectively evaluated targeted busulfan concentrations in 53 pediatric patients (age 0.4-16 years) who received busulfan 4 times daily according to recommended weight-based doses in a single-center analysis between 2018 and 2020. In this trial, individual busulfan pharmacokinetics were performed following dose 5 of the conditioning regimen. Twenty four of 53 patients (45.3%) studies did not require dose adjustments. Equal number of patients (24/53) required one dose adjustments while two-dose adjustment applied for 5 of 53 (9.4%). Twenty-one percent of the patients exhibited ll-lV aGVHD. The incidence of veno-occlusive disease (VOD) was in 3.8% of the 53 patients, while incidence of hemorrhagic cystitis (II-III) reached to 9.7%. Engraftment was successful in 98% of the 53 patients with relapse in 2% of cases. The probability of overall survival and disease-free survival at day 100 was 96% and 94%, respectively. In conclusion, therapeutic drug monitoring (TDM) and individualization of Bu dosage are essential to improve the efficacy and safety of busulfan-based regimen in Chinese pediatric HSCT recipients.

R-praziquantel integrated population pharmacokinetics in preschool- and school-aged African children infected with Schistosoma mansoni and S. haematobium and Lao adults infected with Opisthorchis viverrini

Racemic praziquantel (PZQ) is the standard treatment for schistosomiasis and liver fluke infections (opisthorchiasis and clonorchiasis). The development of an optimal pediatric formulation and dose selection would benefit from a population pharmacokinetic (popPK) model. A popPK model was developed for R-PZQ, the active enantiomer of PZQ, in 664 subjects, 493 African children (2-15 years) infected with Schistosoma mansoni and S. haematobium, and 171 Lao adults (15-78 years) infected with Opisthorchis viverrini. Racemate tablets were administered as single doses of 20, 40 and 60 mg/kg in children and 30, 40 and 50 mg/kg in 129 adults, and as 3 × 25 mg/kg apart in 42 adults. Samples collected by the dried-blood-spot technique were assayed by LC-MS/MS. A two-compartment disposition model, with allometric scaling and dual first-order and transit absorption, was developed using Phoenix™ software. Inversely parallel functions of age described the apparent oral bioavailability (BA) and clearance maturation in children and ageing in adults. BA decreased slightly in children with dose increase, and by 35% in adults with multiple dosing. Crushing tablets for preschool-aged children increased the first-order absorption rate by 64%. The mean transit absorption time was 70% higher in children. A popPK model for R-PZQ integrated African children over 2 years of age with schistosomiasis and Lao adults with opisthorchiasis, and should be useful to support dose optimization in children. In vitro hepatic and intestinal metabolism data would help refining and validating the model in younger children as well as in target ethnic pediatric and adult groups.

Busulfan dose recommendation in inherited metabolic disorders: Population pharmacokinetic analysis

Busulfan is a commonly used alkylating agent in the conditioning regimens of hematopoietic cell transplantation (HCT). Population pharmacokinetic (popPK) models enable the description of busulfan PK and optimization of exposure, which leads to improvement of event-free survival after HCT. Prior busulfan popPK analysis has been limited by small numbers in patients with inherited metabolic disorders (IMD). The primary objective was to characterize population PK of busulfan in a large cohort of children and young adults undergoing HCT for IMD. PopPK analysis of busulfan drug concentrations was performed using data from 78 patients with IMD who received intravenous busulfan (every 24 hours, 4 doses) as part of pretransplantation combination chemotherapy. The final model for busulfan drug clearance was then used to estimate individual doses aimed to achieve a target cumulative area under the curve (cAUC) of 80 to 100 mg · h/L. We then compared the probability of cAUC within the range of 80 to 100 mg · h/L by the developed dosing regimen versus conventional regimen. A 1-compartment, linear elimination model best described the PK of busulfan. Significant covariates demonstrated to affect busulfan clearance included total body weight and the time (in days) from busulfan infusion start. The probability of target cumulative AUC attainment by the developed dosing versus the conventional dosing were 47% versus 43% for body weight <12 kg, and 48% versus 36% for body weight ≥12 kg. We described population PK of intravenous busulfan in a large IMD cohort. The proposed dosing regimen based on the developed model can improve the target cAUC attainment of busulfan for IMD.

Precision dosing of intravenous busulfan in pediatric hematopoietic stem cell transplantation: Results from a multicenter population pharmacokinetic study

Busulfan (Bu) is a common component of conditioning regimens before hematopoietic stem cell transplantation (HSCT) and is known for high interpatient pharmacokinetic (PK) variability. This study aimed to develop and externally validate a multicentric, population PK (PopPK) model for intravenous Bu in pediatric patients before HSCT to first study the influence of glutathione‐s‐transferase A1 (GSTA1) polymorphisms on Bu's PK in a large multicentric pediatric population while accounting for fludarabine (Flu) coadministration and, second, to establish an individualized, model‐based, first‐dose recommendation for intravenous Bu that can be widely used in pediatric patients. The model was built using data from 302 patients from five transplantation centers who received a Bu‐based conditioning regimen. External model validation used data from 100 patients. The relationship between body weight and Bu clearance (CL) was best described by an age‐dependent allometric scaling of a body weight model. A stepwise covariate analysis identified Day 1 of Bu conditioning, GSTA1 metabolic groups based on GSTA1 polymorphisms, and Flu coadministration as significant covariates influencing Bu CL. The final model adequately predicted Bu first‐dose CL in the external cohort, with 81% of predicted area under the curves within the therapeutic window. The final model showed minimal bias (mean prediction error, −0.5%; 95% confidence interval [CI], −3.1% to 2.0%) and acceptable precision (mean absolute prediction error percentage, 18.7%; 95% CI, 17.0%–20.5%) in Bu CL prediction for dosing. This multicentric PopPK study confirmed the influence of GSTA1 polymorphisms and Flu coadministration on Bu CL. The developed model accurately predicted Bu CL and first doses in an external cohort of pediatric patients.

Neuroprotective Efficacy of Prophylactic Antiepileptic Therapies in Busulfan Conditioning Regimen

Busulfan (BU) is commonly used in the anticancer treatment and bone marrow/hematopoietic stem cell transplantation (BM/HSCT) conditioning regimen. This study aimed to evaluate the neuroprotective efficacy of antiepileptic therapies in the BU conditioning regimen. Patients who underwent BM/HSCT between January 1, 2009, and January 1, 2014, were evaluated retrospectively. Diagnosis, demographic characteristics, neurological side effects, electroencephalography (EEG) abnormalities, and antiepileptic drugs were determined. Levetiracetam was given to 49% and phenytoin to 23% of patients as antiepileptic prophylaxis. After BU treatment, seizures developed in three patients, posterior reversible leukoencephalopathy in two patients, and peripheral neuropathy in two patients. Forty-five patients who underwent EEG before and after BU treatment were identified. The effects of antiepileptic prophylaxis over the electroencephalographic abnormalities after the BU conditioning regimen were found to be not different (p = 0.491). Neurological side effects and electroencephalographic abnormalities following BU conditioning regimen under antiepileptic prophylaxis are comparable.

Evaluation of two software using Bayesian methods for monitoring exposure and dosing once-daily intravenous busulfan in paediatric patients receiving haematopoietic stem cell transplantation

AIM

To assess the ability of model-based personalised dosing tools to estimate busulfan exposure (i) in comparison to clinically used intensive sampling exposure estimation procedure, (ii) using limited sampling strategies and (iii) to predict changes in busulfan clearance during busulfan treatment.

METHODS

Data on intravenous busulfan dosing for patients with 4 consecutive days were entered into Bayesian forecasting software, InsightRX and NextDose. Prediction of busulfan cumulative exposure was compared to current clinical practice estimation, aiming for pre-defined individualised target of cumulative exposure. Estimation performance was tested given several limited sampling strategies.

RESULTS

Thirty-two paediatric patients (0.2-16.5 years) provided a total of 103 daily exposure measurements estimated using 7 samples taken per day (full sampling), with 19 patients having sampling following all doses administered. Both software tools utilising Bayesian methods provided acceptable relative bias and precision of cumulative exposure estimations under the tested sampling scenarios. Relative bias ranged from median RE of 0.1-14.6% using InsightRX and from 3.4-7.8% using NextDose. Precision ranged from median RMSE of 0.19-0.32 mg·h·L^-1 for InsightRX and 0.08-0.1 mg·h·L^-1 for NextDose. A median reduction in busulfan clearance from day 1 to day 4 was observed in the clinical data (-10.9%), when using InsightRX (-18.6%) and with NextDose (-14.7%).

CONCLUSION

Bayesian methods were shown to have relatively low bias and precisely estimate busulfan exposure using intensive sampling and several limited sampling strategies, which provides evidence for prospective studies to evaluate these tools in clinical practice. A trend to overestimation of exposure using Bayesian methods was observed compared to clinical practice. Reduction of busulfan clearance from day 1 to 4 of once daily dosing was confirmed and should be considered when adjusting doses.

Optimization of Busulfan Dosing Regimen in Pediatric Patients Using a Population Pharmacokinetic Model Incorporating GST Mutations

Purpose

The aim of this study was to develop a novel busulfan dosing regimen, based on a population pharmacokinetic (PPK) model in Chinese children, and to achieve better area under the concentration-time curve (AUC) targeting.

Patients and Methods

We collected busulfan concentration-time samples from 69 children who received intravenous busulfan prior to allogeneic hematopoietic stem cell transplantation (allo-HSCT). A population pharmacokinetic model for busulfan was developed by nonlinear mixed effect modelling and was validated by an external dataset (n=14). A novel busulfan dosing regimen was developed through simulated patients, and has been verified on real patients. Limited sampling strategy (LSS) was established by Bayesian forecasting. Mean absolute prediction error (MAPE) and relative root mean Squared error (rRMSE) were calculated to evaluate predictive accuracy.

Results

A one-compartment model with first-order elimination best described the data. GSTA1 genotypes, body surface area (BSA) and aspartate aminotransferase (AST) were found to be significant covariates of Bu clearance, and BSA had significant impact of the volume of distribution. Moreover, two equations were obtained for recommended dose regimens: dose (mg)=34.14×BSA (m²)+3.75 (for GSTA1 *A/*A), Dose (mg)=30.99×BSA (m²)+3.21 (for GSTA1 *A/*B). We also presented a piecewise dosage based on BSA categories for each GSTA1 mutation. A two-point LSS, two hours and four hours after dosing, behaved well with acceptable prediction precision (rRMSE=1.026%, MAPE=6.55%).

Conclusion

We recommend a GSTA1-BSA and BSA-based dosing (Q6 h) based on a PPK model for personalizing busulfan therapy in pediatric population. Additionally, an optimal LSS (C_2h and C_4h) provides convenience for therapeutic drug monitoring (TDM) in the future.

Total Body Irradiation Forever? Optimising Chemotherapeutic Options for Irradiation-Free Conditioning for Paediatric Acute Lymphoblastic Leukaemia

Total-body irradiation (TBI) based conditioning prior to allogeneic hematopoietic stem cell transplantation (HSCT) is generally regarded as the gold-standard for children >4 years of age with acute lymphoblastic leukaemia (ALL). Retrospective studies in the 1990's suggested better survival with irradiation, confirmed in a small randomised, prospective study in the early 2000's. Most recently, this was reconfirmed by the early results of the large, randomised, international, phase III FORUM study published in 2020. But we know survivors will suffer a multitude of long-term sequelae after TBI, including second malignancies, neurocognitive, endocrine and cardiometabolic effects. The drive to avoid TBI directs us to continue optimising irradiation-free, myeloablative conditioning. In chemotherapy-based conditioning, the dominant myeloablative effect is provided by the alkylating agents, most commonly busulfan or treosulfan. Busulfan with cyclophosphamide is a long-established alternative to TBI-based conditioning in ALL patients. Substituting fludarabine for cyclophosphamide reduces toxicity, but may not be as effective, prompting the addition of a third agent, such as thiotepa, melphalan, and now clofarabine. For busulfan, it's wide pharmacokinetic (PK) variability and narrow therapeutic window is well-known, with widespread use of therapeutic drug monitoring (TDM) to individualise dosing and control the cumulative busulfan exposure. The development of first-dose selection algorithms has helped achieve early, accurate busulfan levels within the targeted therapeutic window. In the future, predictive genetic variants, associated with differing busulfan exposures and toxicities, could be employed to further tailor individualised busulfan-based conditioning for ALL patients. Treosulfan-based conditioning leads to comparable outcomes to busulfan-based conditioning in paediatric ALL, without the need for TDM to date. Future PK evaluation and modelling may optimise therapy and improve outcome. More recently, the addition of clofarabine to busulfan/fludarabine has shown encouraging results when compared to TBI-based regimens. The combination shows activity in ALL as well as AML and deserves further evaluation. Like busulfan, optimization of chemotherapy conditioning may be enhanced by understanding not just the PK of clofarabine, fludarabine, treosulfan and other agents, but also the pharmacodynamics and pharmacogenetics, ideally in the context of a single disease such as ALL.

Assessment of a Model-Informed Precision Dosing Platform Use in Routine Clinical Care for Personalized Busulfan Therapy in the Pediatric Hematopoietic Cell Transplantation (HCT) Population

Introduction

Population pharmacokinetic (PK) studies demonstrate model-based dosing for busulfan that incorporates body size and age improve clinical target attainment as compared to weight-based regimens. Recently, for clinical dosing of busulfan and TDM, our institution transitioned to a cloud-based clinical decision support tool (www.insight-rx.com). The goal of this study was to assess the dose decision tool for the achievement of target exposure of busulfan in children undergoing hematopoietic cell transplantation (HCT).

Patients and Methods

Patients (N = 188) were grouped into cohorts A, B, or C based on the method for initial dose calculation and estimation of AUC: Cohort A: Initial doses were based on the conventional dosing algorithm (as outlined in the manufacturers' package insert) and non-compartmental analysis (NCA) estimation using the trapezoidal rule for estimation of AUC following TDM. Cohort B: Initial doses for busulfan were estimated by a first-generation PK model and NCA estimation of AUC following TDM. Cohort C: Initial doses were calculated by an updated, second-generation PK model available in the dose decision tool with an estimation of AUC following TDM.

Results

The percent of individuals achieving the exposure target at the time of first PK collection was higher in subjects receiving initial doses provided by the model-informed precision dosing platform (cohort C, 75%) versus subjects receiving initial doses based on either of the two other approaches (conventional guidelines/cohort A, 25%; previous population PK model and NCA parameter estimation, cohort B, 50%). Similarly, the percent of subjects achieving the targeted cumulative busulfan exposure (cAUC) in cohort C was 100% vs. 66% and 88% for cohort A and B, respectively. For cAUC, the variability in the spread of target attainment (%CV) was low at 4.1% for cohort C as compared to cohort A (14.8%) and cohort B (17.1%).

Conclusion

Achievement of goal exposure early on in treatment was improved with the updated model for busulfan and the Bayesian platform. Model-informed dosing and TDM utilizing a Bayesian-based platform provides a significant advantage over conventional guidelines for the achievement of goal cAUC exposure.

Pharmacokinetic analysis for model-supported therapeutic drug monitoring of busulfan in Japanese pediatric hematopoietic stem cell transplantation recipients

This prospective observational study analyzed the pharmacokinetics of busulfan in Japanese children and evaluated the predicting accuracy of previous pediatric PPK models of busulfan. This study enrolled five patients (aged 2-12 years, BW 14-48 kg) receiving a busulfan-based conditioning regimen for hematopoietic stem cell transplantation at our hospital between January 2017 and December 2018. All patients received a 2-hour intravenous busulfan infusion four times daily for a total of 16 doses. After the infusions, 51 plasma samples were collected with the plasma busulfan concentration measured by liquid chromatography-tandem mass spectrometry. PPK model fitting was analyzed using the (%MPE) and the (%MAPE). Limited sampling strategies for estimating busulfan AUC were also evaluated. High interpatient variability was observed in the PK parameters. The most suitable PPK model that reflected our data was McCune's two-compartment model (%MPE -8.7, %MAPE 19.3). A combination sampling method using the busulfan concentration at 2 and 6 hours after the start of the first busulfan dose was found to be able to estimate AUC_4 day . These results provide useful information on busulfan therapeutic drug monitoring in the Japanese pediatric population.

Intra-individual Pharmacokinetic Variability of Intravenous Busulfan in Hematopoietic Stem Cell-Transplanted Children

BACKGROUND

Busulfan therapeutic drug monitoring (TDM) is necessary to better achieve the target exposure in children before hematopoietic stem cell transplantation (HSCT). However, TDM-based dosing may be challenging if intra-individual pharmacokinetic variability (also denoted inter-occasion variability [IOV]) occurs during therapy.

OBJECTIVES

The objectives of this study were to describe and quantify busulfan IOV in children, and to investigate its potential determinants.

METHODS

We performed a new analysis of published data from children who received intravenous busulfan over 4 days before HSCT. We calculated individual pharmacokinetic parameters on each day of therapy using a published population pharmacokinetic model of busulfan and analyzed their changes. Population estimation of IOV was also performed with non-linear mixed effects (NLME) modeling. Potential predictors of significant decrease in busulfan clearance (CL) were assessed by using machine learning approaches.

RESULTS

IOV could be assessed in 136 children. Between day (D) 1 and D2, most patients (80%) experienced a decrease in busulfan CL, with a median change of - 7.9%. However, both large decreases (minimum, -  48.5%) and increases in CL (maximum, + 44%) were observed. Over D1-D3 of therapy, mean CL significantly decreased (-  15%), with a decrease of ≥ 20% in 22% of patients. Some patients also showed unstable CL from day to day. NLME modeling of IOV provided a coefficient of variation of 10.6% and 13.1% for volume of distribution (V_d) and CL, respectively. Some determinants of significant decreases in busulfan CL were identified, but predictive performance of the models was limited.

CONCLUSIONS

Significant busulfan intra-individual variability may occur in children who receive a HSCT and is hardly predictable. The main risk is busulfan overexposure. Performing TDM repeatedly over therapy appears to be the best way to accurately estimate busulfan exposure and perform precision dosing.

Modelling of neutrophil dynamics in children receiving busulfan or treosulfan for haematopoietic stem cell transplant conditioning

AIMS

Busulfan and treosulfan are cytotoxic agents used in the conditioning regime prior to paediatric haematopoietic stem cell transplantation (HSCT). These agents cause suppression of myeloid cells leaving patients severely immunocompromised in the early post-HSCT period. The main objectives were: (i) to establish a mechanistic pharmacokinetic-pharmacodynamic (PKPD) model for the treatment and engraftment effects on neutrophil counts comparing busulfan and treosulfan-based conditioning, and (ii) to explore current dosing schedules with respect to time to HSCT.

METHODS

Data on 126 patients, 72 receiving busulfan (7 months-18 years, 5.1-47.0 kg) and 54 treosulfan (4 months-17 years, 3.8-35.8 kg), were collected. In total, 8935 neutrophil count observations were recorded during the study period in addition to drug concentrations to develop a mechanistic PKPD model. Absolute neutrophil count profiles were modelled semimechanistically, accounting for transplant effects and differing set points pre- and post-transplant.

RESULTS

PK were best described by 2-compartment models for both drugs. The Friberg semimechanistic neutropenia model was applied with a linear model for busulfan and a maximum efficacy model for treosulfan describing drug effects at various stages of neutrophil maturation. System parameters were consistent across both drugs. The HSCT was represented by an amount of progenitor cells enhancing the neutrophils' proliferation and maturation compartments. Alemtuzumab was found to enhance the proliferative rate under which the absolute neutrophil count begin to grow after HSCT.

CONCLUSION

A semimechanistic PKPD model linking exposure to either busulfan or treosulfan to the neutrophil reconstitution dynamics was successfully built. Alemtuzumab coadministration enhanced the neutrophil proliferative rate after HSCT. Treosulfan administration was suggested to be delayed with respect to time to HSCT, leaving less time between the end of the administration and stem cell infusion.

Population pharmacokinetic modelling of busulfan and the influence of body composition in paediatric Fanconi anaemia patients

Aims

Fanconi anaemia (FA) is a rare disorder characterized by progressive bone marrow failure that requires haematopoietic cell transplantation (HCT). Busulfan is used in conditioning regimens prior to HCT. Doses used in non‐FA patients cause life‐threatening toxicities in FA patients and data on busulfan pharmacokinetics (PK) in this population are limited. This study characterized busulfan PK in paediatric FA patients using population PK modelling and evaluated the effect of body composition on steady‐state concentrations (C_ss).

Methods

A total of 200 busulfan plasma concentrations in 29 FA patients from a recent study (http://Clinicaltrials.gov; NCT01082133) were available for population PK modelling. The effect of different body size‐scaled doses and body compositions on C_ss was investigated using population PK modelling.

Results

Fat free mass (FFM) was identified as the best size descriptor in a two‐compartment busulfan PK model in FA patients. Conventional dosing, based on an amount of busulfan per kilogram of total body mass, resulted in higher C_ss in FA patients with higher body mass index (BMI). A newly proposed FFM‐based dosing strategy would eliminate the observed trend of higher concentrations in high BMI patients, and achieve consistent C_ss across a wide BMI spectrum.

Conclusions

This is the first study to describe the population PK of busulfan in paediatric FA patients. The proposed model will facilitate PK model‐informed precision dosing. FFM‐based dosing is expected to improve the probability of achieving target C_ss, particularly in obese patients, while minimizing the risk of overdosing.

Pharmacotherapy in Children and Adolescents: Oncology

Pharmacotherapy in paediatric oncology is a difficult task. It is challenging to determine the optimal dose in children of different age groups. In addition, anticancer drugs display severe side effects reducing the quality of life. Late effects like secondary tumours and cardiotoxicity can be apparent years after treatment and must be taken into account when planning treatment schedules. Classical cytoreducing agents are still of great importance in treating children with leukaemia and solid tumours. In addition, drugs developed by rational drug design (targeted drugs) are a very important part of many treatment protocols, and newer drugs are emerging in several types of cancer. Unfortunately, there is only limited experience with newer drugs in children, because new drugs are mostly developed for adults. Complicated therapy regimens require a solid knowledge of the pharmacology of the drugs applied. This chapter attempts to introduce some pharmacological knowledge for the most important anticancer drugs in children with a focus on side effects and age-specific considerations.

Population Pharmacokinetics and Safety of Solithromycin following Intravenous and Oral Administration in Infants, Children, and Adolescents

Solithromycin is a novel fluoroketolide antibiotic which was under investigation for the treatment of community-acquired bacterial pneumonia (CABP). A phase 1 study was performed to characterize the pharmacokinetics (PK) and safety of solithromycin in children. Eighty-four subjects (median age, 6 years [age range, 4 days to 17 years]) were administered intravenous (i.v.) or oral (capsules or suspension) solithromycin (i.v., 6 to 8 mg/kg of body weight; capsules/suspension, 14 to 16 mg/kg on days 1 and 7 to 15 mg/kg on days 2 to 5). PK samples were collected after the first and multidose administration. Data from 83 subjects (662 samples) were combined with previously collected adolescent PK data (n = 13; median age, 16 years [age range, 12 to 17 years]) following capsule administration to perform a population PK analysis. A 2-compartment PK model characterized the data well, and postmenstrual age was the only significant covariate after accounting for body size differences. Dosing simulations suggested that 8 mg/kg i.v. daily and oral dosing of 20 mg/kg on day 1 (800-mg adult maximum) followed by 10 mg/kg on days 2 to 5 (400-mg adult maximum) would achieve a pediatric solithromycin exposure consistent with the exposures observed in adults. Seventy-six treatment-emergent adverse events (TEAEs) were reported in 40 subjects. Diarrhea (6 subjects) and infusion site pain or phlebitis (3 subjects) were the most frequently reported adverse events related to treatment. Two subjects experienced TEAEs of increased hepatic enzymes that were deemed not to be related to the study treatment. (The phase 1 pediatric studies discussed in this paper have been registered at ClinicalTrials.gov under identifiers NCT01966055 and NCT02268279.).

Incorporation of GSTA1 genetic variations into a population pharmacokinetic model for IV busulfan in paediatric hematopoietic stem cell transplantation

AIMS

The aim of this study is to develop a population pharmacokinetic (PopPK) model for intravenous busulfan in children that incorporates variants of GSTA1, gene coding for the main enzyme in busulfan metabolism.

METHODS

Busulfan concentration-time data was collected from 112 children and adolescents (median 5.4 years old, range: 0.1-20) who received intravenous busulfan during the conditioning regimen prior to stem cell transplantation. Weight, sex, baseline disease (malignant vs. non-malignant), age, conditioning regimen and GSTA1 diplotypes were evaluated as covariates of pharmacokinetic parameters by using nonlinear mixed effects analysis. The ability to achieve the target AUC_24h (3600-6000 μM min^-1 ) was assessed by estimating the first dose based on the present PopPK model and by comparing the results with other available models in children.

RESULTS

A one-compartment model with first-order elimination best described the data. Allometric scaling of weight and a factor of busulfan metabolism maturation were included in the base model. GSTA1 diplotypes were found to be a significant covariate of busulfan clearance, which was 7% faster in rapid metabolizers and 12% slower in poor metabolizers, in comparison with normal ones. Busulfan doses calculated using the parameters of the proposed PopPK model were estimated to achieve the target AUC in 85.2% of the cases (95% CI 78.7-91.7%).

CONCLUSION

This is the first PopPK for busulfan that successfully incorporated GSTA1 genotype in a paediatric population. Its use may contribute to better prediction of busulfan exposure in children and adolescents since the first dose, by tailoring the dose according to the individual metabolic capacity.

Population pharmacokinetics analysis of intravenous busulfan in Chinese patients undergoing hematopoietic stem cell transplantation

There are several reports describing population pharmacokinetic (popPK) models of busulfan (BU). However, limited information is available in Chinese hematopoietic stem cell transplantation (HSCT) patients. The present study aimed to establish a popPK model of intravenous BU in Chinese HSCT patients for individualized drug therapy. The popPK model of BU was developed from a total of 284 concentration-time points from 53 patients. The effects of demographic and biochemical covariates were investigated by nonlinear mixed effect model (NONMEM) software. Plots, visual predictive check (VPC), bootstrap and normalized prediction distribution error (NPDE) were performed to determine the stability and the reliability of the final model. A one-compartment model with first-order elimination process was confirmed as the final structural model for BU. For a typical patient whose body surface area (BSA) is 1.7 m² , the population typical values of CL and Vd were 11.86 L/h, and 48.2 L, respectively. The result suggested BSA showed significant influence on CL and Vd (P<.001). Plots revealed the final model was performing a goodness fit. The steady rate verified by bootstrap was 100%, relative deviation was less than 4.00%, estimated value of final model was in the 95% confidence interval (CI). The VPC results showed the observed values were almost all positioned within the 5th and 95th CIs. The mean and variance of the NPDE were 0.0363 (Wilcoxon signed-rank test, 0.298) and 0.877 (Fisher variance test, 0.134; SW test of normality, 0.108), respectively. The global adjusted P value was 0.305, which indicated that the prediction of the BU popPK model was adequate. A physician-friendly Microsoft Excel-base tool was implemented using the final popPK model for designing individualized dosing regimens.

Population Pharmacokinetics/Pharmacodynamics of 3,4-Diaminopyridine Free Base in Patients With Lambert-Eaton Myasthenia

Lambert-Eaton myasthenia (LEM) is a rare autoimmune disorder associated with debilitating muscle weakness. There are limited treatment options and 3,4-diaminopyridine (3,4-DAP) free base is an investigational orphan drug used to treat LEM-related weakness. We performed a population pharmacokinetic/pharmacodynamic (PK/PD) analysis using 3,4-DAP and metabolite concentrations collected from a phase II study in patients with LEM. The Triple Timed Up & Go (3TUG) assessment, which measures lower extremity weakness, was the primary outcome measure. A total of 1,270 PK samples (49 patients) and 1,091 3TUG data points (32 randomized patients) were included in the PK/PD analysis. A two-compartment and one-compartment model for parent and metabolite, respectively, described the PK data well. Body weight and serum creatinine partially explained the variability in clearance for the final PK model. A fractional inhibitory maximum effect (Emax ) model characterized the exposure-response relationship well. The PK/PD model was applied to identify a suggested dosing approach for 3,4-DAP free base.

GSTA1 Genetic Variants and Conditioning Regimen: Missing Key Factors in Dosing Guidelines of Busulfan in Pediatric Hematopoietic Stem Cell Transplantation

Busulfan (Bu) is a key component of conditioning regimens used before hematopoietic stem cell transplantation (SCT) in children. Different predictive methods have been used to calculate the first dose of Bu. To evaluate the necessity of further improvements, we retrospectively analyzed the currently available weight- and age-based guidelines to calculate the first doses in 101 children who underwent allogenic SCT in CHU Sainte-Justine, Montreal, after an intravenous Bu-containing conditioning regimen according to genetic and clinical factors. The measured areas under the curve (AUCs) were within target (900 to 1500 µM/min) in 38.7% of patients after the administration of the first dose calculated based on age and weight, as locally recommended. GSTA1 diplotypes linked to poor Bu metabolism (G3) and fludarabine-containing regimens were the only factors associated with AUC within target (OR, 4.7 [95% CI, 1.1 to 19.8, P = .04]; and OR, 9.9 [95% CI, 1.6 to 61.7, P = .01], respectively). From the 11 methods selected for dose calculation, the percentage of AUCs within the target varied between 16% and 74%. In some models G3 was associated with AUCs within the therapeutic and the toxic range, whereas rapid metabolizers (G1) were correlated with subtherapeutic AUCs when different methods were used. These associations were confirmed by clearance-prediction analysis, in which GSTA1 diplotypes consistently influenced the prediction errors of the methods. These findings suggest that these factors should be considered in Bu dose prediction in addition to the anthropometric data from patients. Furthermore, our data indicated that GSTA1 diplotypes was a factor that should be included in future population pharmacokinetic models, including similar conditioning regiments, to improve the prediction of Bu exposure after its initial dose.

Pediatric patients undergoing hematopoietic stem cell transplantation can greatly benefit from a novel once-daily intravenous busulfan dosing nomogram

Busulfan, a bifunctional alkylating agent, has been used as a conditioning regimen prior to allogeneic hematopoietic stem cell transplantation (HSCT). The aim of this study was to derive a novel once-daily intravenous (IV) busulfan dosing nomogram for pediatric patients undergoing HSCT using a population pharmacokinetic (PK) model. A population PK analysis was performed using 2183 busulfan concentrations in 137 pediatric patients (age: 0.6-22.2 years), who received IV busulfan once-daily for 4 days before undergoing HSCT. Based on the final population PK model, an optimal once-daily IV busulfan dosing nomogram was derived. The percentage of simulated patients achieving the daily target area under the concentration-time curve (AUC) by the new nomogram was compared with that by other busulfan dosing regimens including the FDA regimen, the EMA regimen, and the empirical once-daily regimen without therapeutic drug monitoring (TDM). A one-compartment open linear PK model incorporating patient's body surface area, age, dosing day, and aspartate aminotransferase as a significant covariate adequately described the concentration-time profiles of busulfan. An optimal dosing nomogram based on the PK model performed significantly better than the other dosing regimens, resulting in >60% of patients achieving the target AUC while the percentage of patients exceeding the toxic AUC level was kept <25% during the entire treatment period. A novel once-daily busulfan dosing nomogram for pediatric patients undergoing HSCT is useful for clinicians, particularly in a setting where TDM service is not readily available or to optimize the dose on day 1.

Paediatric Drug Development and Formulation Design-a European Perspective

The availability of licensed paediatric drugs is lagging behind those for adults, and there is a lack of safe formulations in suitable doses that children are able and willing to take. As a consequence, children are commonly treated with off-label or unlicensed drugs. As off-label and unlicensed drug use are associated with a greater risk for harm than on-label drug use, a range of global initiatives have been developed to realize "better" medicines for children. This review describes the challenges and achievements of the European Union to realize this goal, with a focus on paediatric drug development and formulation design. In 2007, a European Paediatric Regulation was installed enforcing companies to consider children in the early development of drugs with a new drug substance, for a new indication or with a new route of administration. The Regulation, e.g. requires companies to develop a paediatric investigation plan discussing the proposed clinical trials in children of different ages and the formulations for future marketing. Since 2013, the pharmaceutical design of any newly marketed paediatric drug should comply with the "Guideline on the Pharmaceutical Development of Medicines for Paediatric Use." Companies should, e.g. justify the route of administration, dosage form, formulation characteristics, safety of excipients, dosing frequency, container closure system, administration device, patient acceptability and user information. In this review, the guideline's key aspects are discussed with a focus on novel formulations such as mini-tablets and orodispersible films, excipients with a potential risk for harm such as azo dyes and adequate user instructions.

Pharmacokinetics, Pharmacodynamics, and Pharmacogenomics of Immunosuppressants in Allogeneic Hematopoietic Cell Transplantation: Part II

Part I of this article included a pertinent review of allogeneic hematopoietic cell transplantation (alloHCT), the role of postgraft immunosuppression in alloHCT, and the pharmacokinetics, pharmacodynamics, and pharmacogenomics of the calcineurin inhibitors and methotrexate. In this article (Part II), we review the pharmacokinetics, pharmacodynamics, and pharmacogenomics of mycophenolic acid (MPA), sirolimus, and the antithymocyte globulins (ATG). We then discuss target concentration intervention (TCI) of these postgraft immunosuppressants in alloHCT patients, with a focus on current evidence for TCI and on how TCI may improve clinical management in these patients. Currently, TCI using trough concentrations is conducted for sirolimus in alloHCT patients. Several studies demonstrate that MPA plasma exposure is associated with clinical outcomes, with an increasing number of alloHCT patients needing TCI of MPA. Compared with MPA, there are fewer pharmacokinetic/dynamic studies of rabbit ATG and horse ATG in alloHCT patients. Future pharmacokinetic/dynamic research of postgraft immunosuppressants should include '-omics'-based tools: pharmacogenomics may be used to gain an improved understanding of the covariates influencing pharmacokinetics as well as proteomics and metabolomics as novel methods to elucidate pharmacodynamic responses.

Population pharmacokinetics of busulfan in pediatric and young adult patients undergoing hematopoietic cell transplant: a model-based dosing algorithm for personalized therapy and implementation into routine clinical use

BACKGROUND

Population pharmacokinetic (PK) studies of busulfan in children have shown that individualized model-based algorithms provide improved targeted busulfan therapy when compared with conventional dose guidelines. The adoption of population PK models into routine clinical practice has been hampered by the tendency of pharmacologists to develop complex models too impractical for clinicians to use. The authors aimed to develop a population PK model for busulfan in children that can reliably achieve therapeutic exposure (concentration at steady state) and implement a simple model-based tool for the initial dosing of busulfan in children undergoing hematopoietic cell transplantation.

PATIENTS AND METHODS

Model development was conducted using retrospective data available in 90 pediatric and young adult patients who had undergone hematopoietic cell transplantation with busulfan conditioning. Busulfan drug levels and potential covariates influencing drug exposure were analyzed using the nonlinear mixed effects modeling software, NONMEM. The final population PK model was implemented into a clinician-friendly Microsoft Excel-based tool and used to recommend initial doses of busulfan in a group of 21 pediatric patients prospectively dosed based on the population PK model.

RESULTS

Modeling of busulfan time-concentration data indicates that busulfan clearance displays nonlinearity in children, decreasing up to approximately 20% between the concentrations of 250-2000 ng/mL. Important patient-specific covariates found to significantly impact busulfan clearance were actual body weight and age. The percentage of individuals achieving a therapeutic concentration at steady state was significantly higher in subjects receiving initial doses based on the population PK model (81%) than in historical controls dosed on conventional guidelines (52%) (P = 0.02).

CONCLUSIONS

When compared with the conventional dosing guidelines, the model-based algorithm demonstrates significant improvement for providing targeted busulfan therapy in children and young adults.

Busulfan dosing algorithm and sampling strategy in stem cell transplantation patients

AIM

The aim of this investigation was to develop a model-based dosing algorithm for busulfan and identify an optimal sampling scheme for use in routine clinical practice.

METHODS

Clinical data from an ongoing study (n = 29) in stem cell transplantation patients were used for the purposes our analysis. A one compartment model was selected as basis for sampling optimization and subsequent evaluation of a suitable dosing algorithm. Internal and external model validation procedures were performed prior to the optimization steps using ED-optimality criteria. Using systemic exposure as parameter of interest, dosing algorithms were considered for individual patients with the scope of minimizing the deviation from target range as determined by AUC(0,6 h).

RESULTS

Busulfan exposure after oral administration was best predicted after the inclusion of adjusted ideal body weight and alanine transferase as covariates on clearance. Population parameter estimates were 3.98 h(-1), 48.8 l and 12.3 l h(-1) for the absorption rate constant, volume of distribution and oral clearance, respectively. Inter-occasion variability was used to describe the differences between test dose and treatment. Based on simulation scenarios, a dosing algorithm was identified, which ensures target exposure values are attained after a test dose. Moreover, our findings show that a sparse sampling scheme with five samples per patient is sufficient to characterize the pharmacokinetics of busulfan in individual patients.

CONCLUSION

The use of the proposed dosing algorithm in conjunction with a sparse sampling scheme may contribute to considerable improvement in the safety and efficacy profile of patients undergoing treatment for stem cell transplantation.

Effect of age on the pharmacokinetics of busulfan in patients undergoing hematopoietic cell transplantation; an alliance study (CALGB 10503, 19808, and 100103)

PURPOSE

Older patients with acute myeloid leukemia (AML) and myelodysplastic syndrome have often been excluded from myeloablative-conditioning regimens containing busulfan because of non-disease-related morbidity and mortality. We hypothesized that busulfan clearance (BuCL) in older patients (>60 years) would be reduced compared to that in younger patients, potentially explaining observed differences in busulfan tolerability.

METHODS

AML patients in three CALGB hematopoietic cell transplantation studies were treated with a conditioning regimen using IV busulfan, dosed at 0.8 mg/kg. Plasma busulfan concentrations were determined by LC-MS and analyzed by non-compartmental methods. BuCL was normalized to actual (ABW), ideal (IBW), or corrected (CBW) body weight (kg). Differences in BuCL between age groups were examined using the Wilcoxon rank sum test.

RESULTS

One hundred and eighty-five patients were accrued; 174 provided useable pharmacokinetic data. Twenty-nine patients ≥ 60 years old (median 66; range 60-74) had a significantly higher BuCL versus those <60 years old (median 50; range 18-60): BuCL 236 versus 168 mL/min, p = 0.0002; BuCL/ABW 3.0 versus 2.1 mL/min/kg, p = 0.0001; BuCL/IBW 3.8 versus 2.6 mL/min/kg, p = 0.0035; BuCL/CBW 3.4 versus 2.6 mL/min/kg, p = 0.0005. Inter-patient variability in clearance (CV %) was up to 48 % in both age groups. Phenytoin administration, a potential confounder, did not affect BuCL, regardless of weight normalization (p > 0.34).

CONCLUSIONS

Contrary to our hypothesis, BuCL was significantly higher in older patients compared to younger patients in these studies and does not explain the previously reported increase in busulfan toxicity observed in older patients.

Effect of genetic variants GSTA1 and CYP39A1 and age on busulfan clearance in pediatric patients undergoing hematopoietic stem cell transplantation

BACKGROUND

Busulfan is used in preparative regimens prior to stem cell transplantation in pediatric patients. There is significant interpatient variability in busulfan pharmacokinetics (PK) and exposure is related to outcome. To date, only polymorphisms in genes encoding for glutathione-S-transferases were studied, but could only explain a small portion of the variability in PK.

AIM

To investigate the effect of seven genetic markers on busulfan clearance and the effect of ontogenesis on these genetic variants in a pediatric population.

MATERIALS & METHODS

In an earlier study of our group seven genetic markers in GSTA1, CYP2C19, CYP39A1, ABCB4, SLC22A4 and SLC7A8 were associated with busulfan clearance in adult patients. Eighty four pediatric patients were genotyped for these markers and genotype was associated with busulfan clearance.

RESULTS & CONCLUSION

GSTA1 and CYP39A1 were found to be associated with busulfan clearance. When combined, the two haplotypes explained 17% of the variability in busulfan clearance. Furthermore, the effect of GSTA1 haplotype on clearance was dependent on age.

Optimizing drug development of anti-cancer drugs in children using modelling and simulation

Modelling and simulation (M&S)-based approaches have been proposed to support paediatric drug development in order to design and analyze clinical studies efficiently. Development of anti-cancer drugs in the paediatric population is particularly challenging due to ethical and practical constraints. We aimed to review the application of M&S in the development of anti-cancer drugs in the paediatric population, and to identify where M&S-based approaches could provide additional support in paediatric drug development of anti-cancer drugs. A structured literature search on PubMed was performed. The majority of identified M&S-based studies aimed to use population PK modelling approaches to identify determinants of inter-individual variability, in order to optimize dosing regimens and to develop therapeutic drug monitoring strategies. Prospective applications of M&S approaches for PK-bridging studies have scarcely been reported for paediatric oncology. Based on recent developments of M&S in drug development there are several opportunities where M&S could support more informative bridging between children and adults, and increase efficiency of the design and analysis of paediatric clinical trials, which should ultimately lead to further optimization of drug treatment strategies in this population.

Evaluation of effects of busulfan and DMA on SOS in pediatric stem cell recipients

BACKGROUND

Busulfan (Bu) is a DNA-alkylating agent used for myeloablative conditioning in stem cell transplantation in children and adults. While the use of intravenous rather than oral administration of Bu has reduced inter-individual variability in plasma levels, toxicity still occurs frequently after hematopoietic stem cell transplantation (HSCT). Toxicity (especially hepatotoxic effects) of intravenous (IV) Bu may be related to both Bu and/or N,N-dimethylacetamide (DMA), the solvent of Bu. In this study, we assessed the relation between the exposure of Bu and DMA with regards to the clinical outcome in children from two cohorts.

METHODS

In a two-centre study Bu and DMA AUC (area under the curve) were correlated in pediatric stem cell recipients to the risk of developing SOS and to the clinical outcome.

RESULT

In patients receiving Bu four times per day Bu levels >1,500 µmol/L minute correlate to an increased risk of developing a SOS. In the collective cohort, summarizing data of all 53 patients of this study, neither high area under the curve (AUC) of Bu nor high AUC of DMA appears to be an independent risk factor for the development of SOS in children.

CONCLUSION

In this study neither Bu nor DMA was observed as an independent risk factor for the development of SOS. To identify subgroups (e.g., infants), in which Bu or DMA might be risk factors for the induction of SOS, larger cohorts have to be evaluated.

Busulfan in infant to adult hematopoietic cell transplant recipients: a population pharmacokinetic model for initial and Bayesian dose personalization

PURPOSE

Personalizing intravenous busulfan doses to a target plasma concentration at steady state (Css) is an essential component of hematopoietic cell transplantation (HCT). We sought to develop a population pharmacokinetic model to predict i.v. busulfan doses over a wide age spectrum (0.1-66 years) that accounts for differences in age and body size.

EXPERIMENTAL DESIGN

A population pharmacokinetic model based on normal fat mass and maturation based on postmenstrual age was built from 12,380 busulfan concentration time points obtained after i.v. busulfan administration in 1,610 HCT recipients. Subsequently, simulation results of the initial dose necessary to achieve a target Css with this model were compared with pediatric-only models.

RESULTS

A two-compartment model with first-order elimination best fit the data. The population busulfan clearance was 12.4 L/h for an adult male with 62 kg normal fat mass (equivalent to 70 kg total body weight). Busulfan clearance, scaled to body size-specifically normal fat mass, is predicted to be 95% of the adult clearance at 2.5 years postnatal age. With a target Css of 770 ng/mL, a higher proportion of initial doses achieved the therapeutic window with this age- and size--dependent model (72%) compared with dosing recommended by the U.S. Food and Drug Administration (57%) or the European Medicines Agency (70%).

CONCLUSION

This is the first population pharmacokinetic model developed to predict initial i.v. busulfan doses and personalize to a target Css over a wide age spectrum, ranging from infants to adults.

Effect of weight and maturation on busulfan clearance in infants and small children undergoing hematopoietic cell transplantation

Little information is currently available regarding the pharmacokinetics (PK) of busulfan in infants and small children to help guide decisions for safe and efficacious drug therapy. The objective of this study was to develop an algorithm for individualized dosing of i.v. busulfan in infants and children weighing ≤12 kg, that would achieve targeted exposure with the first dose of busulfan. Population PK modeling was conducted using intensive time-concentration data collected through the routine therapeutic drug monitoring of busulfan in 149 patients from 8 centers. Busulfan PK was well described by a 1-compartment base model with linear elimination. The important clinical covariates affecting busulfan PK were actual body weight and age. Based on our model, the predicted clearance of busulfan increases approximately 1.7-fold between 6 weeks to 2 years of life. For infants age <5 months, the model-predicted doses (mg/kg) required to achieve a therapeutic concentration at steady state of 600-900 ng/mL (area under the curve range, 900-1350 μM·min) were much lower compared with standard busulfan doses of 1.1 mg/kg. These results could help guide clinicians and inform better dosing decisions for busulfan in young infants and small children undergoing hematopoietic cell transplantation.

Variation in prescribing patterns and therapeutic drug monitoring of intravenous busulfan in pediatric hematopoietic cell transplant recipients

Personalizing intravenous (IV) busulfan doses in children using therapeutic drug monitoring (TDM) is an integral component of hematopoietic cell transplant. The authors sought to characterize initial dosing and TDM of IV busulfan, along with factors associated with busulfan clearance, in 729 children who underwent busulfan TDM from December 2005 to December 2008. The initial IV busulfan dose in children weighing ≤12 kg ranged 4.8-fold, with only 19% prescribed the package insert dose of 1.1 mg/kg. In those children weighing >12 kg, the initial dose ranged 5.4-fold, and 79% were prescribed the package insert dose. The initial busulfan dose achieved the target exposure in only 24.3% of children. A wide range of busulfan exposures were targeted for children with the same disease (eg, 39 target busulfan exposures for the 264 children diagnosed with acute myeloid leukemia). Considerable heterogeneity exists regarding when TDM is conducted and the number of pharmacokinetic samples obtained. Busulfan clearance varied by age and dosing frequency but not by underlying disease. The authors- group is currently evaluating how using population pharmacokinetics to optimize initial busulfan dose and TDM (eg, limited sampling schedule in conjunction with maximum a posteriori Bayesian estimation) may affect clinical outcomes in children.

Body weight-dependent pharmacokinetics of busulfan in paediatric haematopoietic stem cell transplantation patients: towards individualized dosing

BACKGROUND AND OBJECTIVES

The wide variability in pharmacokinetics of busulfan in children is one factor influencing outcomes such as toxicity and event-free survival. A meta-analysis was conducted to describe the pharmacokinetics of busulfan in patients from 0.1 to 26 years of age, elucidate patient characteristics that explain the variability in exposure between patients and optimize dosing accordingly.

PATIENTS AND METHODS

Data were collected from 245 consecutive patients (from 3 to 100 kg) who underwent haematopoietic stem cell transplantation (HSCT) in four participating centres. The inter-patient, inter-occasion and residual variability in the pharmacokinetics of busulfan were estimated with a population analysis using the nonlinear mixed-effects modelling software NONMEM VI. Covariates were selected on the basis of their known or theoretical relationships with busulfan pharmacokinetics and were plotted independently against the individual pharmacokinetic parameters and the weighted residuals of the model without covariates to visualize relations. Potential covariates were formally tested in the model.

RESULTS

In a two-compartment model, body weight was the most predictive covariate for clearance, volume of distribution and inter-compartmental clearance and explained 65%, 75% and 40% of the observed variability, respectively. The relationship between body weight and clearance was characterized best using an allometric equation with a scaling exponent that changed with body weight from 1.2 in neonates to 0.55 in young adults. This implies that an increase in body weight in neonates results in a larger increase in busulfan clearance than an increase in body weight in older children or adults. Clearance on the first day was 12% higher than that of subsequent days (p < 0.001). Inter-occasion variability on clearance was 15% between the 4 days. Based on the final pharmacokinetic-model, an individualized dosing nomogram was developed.

CONCLUSIONS

The model-based individual dosing nomogram is expected to result in predictive busulfan exposures in patients ranging between 3 and 65 kg and thereby to a safer and more effective conditioning regimen for HSCT in children.