Population Pharmacokinetics of Intravenous Busulfan in Japanese Pediatric Patients With Primary Immunodeficiency Diseases

Population Pharmacokinetic Modeling for Twice-Daily Intravenous Busulfan in a Large Cohort of Pediatric Patients Undergoing Hematopoietic Stem Cell Transplantation-A 10-Year Single-Center Experience

Reaching target exposure of busulfan-based conditioning prior to hematopoietic stem cell transplantation is vital for favorable therapy outcomes. Yet, a wide inter-patient and inter-occasion variability in busulfan exposure has been reported, especially in children. We aimed to identify factors associated with the variability of busulfan pharmacokinetics in 124 consecutive patients transplanted at the University Children's Hospital Zurich between October 2010 and February 2020. Clinical data and busulfan plasma levels after twice-daily intravenous administration were analyzed retrospectively by population pharmacokinetic modeling. The volume of distribution correlated with total body water. The elimination rate constant followed an age-dependent maturation function, as previously suggested, and correlated with the levels of serum albumin. Acute lymphoblastic leukemia reduced busulfan clearance by 20%. Clearance significantly decreased by 17% on average from the start to the third day of busulfan administration, in agreement with other studies. An average reduction of 31% was found in patients with hemophagocytic lymphohistiocytosis and X-linked lymphoproliferative disease. In conclusion, we demonstrate that in addition to known factors, underlying disease and serum albumin significantly impact busulfan pharmacokinetics in pediatric patients; yet, substantial unexplained variability in some patients remained. Thus, we consider repeated pharmacokinetic assessment essential to achieve the desired target exposure in twice-daily busulfan administration.

Predicting busulfan exposure in patients undergoing hematopoietic stem cell transplantation using machine learning techniques

PURPOSE

This study aimed to establish an optimal model to predict the busulfan (BU) area under the curve at steady state (AUCss) by using machine learning (ML).

PATIENTS AND METHODS

Seventy-nine adult patients (age ≥18 years) who received BU intravenously and underwent therapeutic drug monitoring from 2013 to 2021 at Fujian Medical University Union Hospital were enrolled in this retrospective study. The whole dataset was divided into a training group and test group at the ratio of 8:2. BU AUCss were considered as the target variable. Nine different ML algorithms and one population pharmacokinetic (pop PK) model were developed and validated, and their predictive performance was compared.

RESULTS

All ML models were superior to the pop PK model (R2 = 0.751, MSE = 0.722, 14 and RMSE = 0.830) in model fitting and had better predictive accuracy. The ML model of BU AUCss established through support vector regression (SVR) and gradient boosted regression trees (GBRT) had the best predictive ability (R2 = 0.953 and 0.953, MSE = 0.323 and 0.326, and RMSE = 0.423 and 0.425).

CONCLUSION

All the ML models can potentially be used to estimate BU AUCss with the aim of facilitating rational use of BU on the individualized level, especially models built by SVR and GBRT algorithms.

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.

Conditioning regimens for inborn errors of immunity: current perspectives and future strategies

Inborn errors of immunity (IEI) are caused by germline genetic mutations, resulting in defects of innate or acquired immunity. Hematopoietic cell transplantation (HCT) is indicated for curative therapy especially in patients with IEI who develop fatal opportunistic infections or severe manifestations of immune dysregulation. The first successful HCT for severe combined immunodeficiency (SCID) was reported in 1968. Since then, the indications for HCT have expanded from SCID to various non-SCID IEI. In general, HCT for IEI differs from that for other hematological malignancies in that the goal is not to eradicate certain immune cells but to achieve immune reconstitution. European Society for Blood and Marrow Transplantation/European Society for Immunodeficiencies guidelines recommend reduced-intensity conditioning to avoid treatment-related toxicity, and the optimal conditioning regimen should be considered for each IEI. We review conditioning regimens for some representative IEI disorders in Japanese and worldwide cohort studies, and future strategies for treating IEI.

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.

Can Published Population Pharmacokinetic Models of Busulfan Be Used for Individualized Dosing in Chinese Pediatric Patients Undergoing Hematopoietic Stem Cell Transplantation? An External Evaluation

Busulfan is a bifunctional alkylating agent that is widely used before hematopoietic stem cell transplantation (HSCT), in combination with other chemotherapeutic drugs. As of 2020, there is no population pharmacokinetic (popPK) model for busulfan in Chinese pediatric patients. A systemic external evaluation of 11 published popPK models was conducted in Chinese pediatric patients undergoing HSCT. Forty pediatric patients were enrolled in this study, with a total of 183 blood concentrations. The relative prediction error (PE%), median PE%, median absolute PE%, and percentage of PE% within ±20% and ±30% were calculated in prediction-based diagnostics. Simulation-based diagnostics were conducted through a prediction- and variability-corrected visual predictive check and the normalized prediction distribution error. The relative individual prediction error was calculated using Bayesian forecasting with 1 to 3 concentration points. The 1-compartment open linear popPK model, which was built by Su-jin Rhee et al (model H), incorporating the patient's body surface area, age, dosing day, and aspartate aminotransferase as significant covariates had preferable predictability than other popPK models. In prediction-based diagnostics, the median PE%, percentage of PE% within ±20%, and percentage of PE% within ±30% of model H were 8.48%, 45.35%, and 59.56%, respectively. The normalized prediction distribution error of model H showed that it followed the normal distribution. Based on Bayesian forecasting, model H showed good predictive performance. Thus, model H was the most appropriate model that can be used clinically for individualized dosage adjustments in Chinese pediatric HSCT patients.

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.

Hematopoietic Cell Transplantation with Reduced Intensity Conditioning Using Fludarabine/Busulfan or Fludarabine/Melphalan for Primary Immunodeficiency Diseases

PURPOSE

The purpose of our study was to compare the safety and efficacy of hematopoietic cell transplantation (HCT) using fludarabine (Flu)-based reduced intensity conditioning (RIC) with busulfan (BU) or melphalan (Mel) for primary immunodeficiency diseases (PID).

METHODS

We retrospectively analyzed transplant outcome, including engraftment, chimerism, immune reconstitution, and complications in 15 patients with severe combined immunodeficiency (SCID) and 27 patients with non-SCID PID. The patients underwent Flu-based RIC-HCT with BU (FluBU: 7 SCID, 16 non-SCID) or Mel (FluMel: 8 SCID, 11 non-SCID). The targeted low-dose BU with therapeutic drug monitoring was set to 30 mg hour/L for SCID.

RESULTS

The 2-year overall survival of all patients was 79.6% and that of patients with SCID in the FluBU and FluMel groups was 100% and 62.5%, respectively. In the FluBU group, all seven patients achieved engraftment, good immune reconstitution, and long-term survival. All five patients receiving umbilical cord blood transplantation achieved complete or high-level mixed chimerism and sufficient specific IgG production. In the FluMel group, six of eight patients achieved complete or high-level mixed chimerism. Viral reactivation or new viral infection occurred in one FluBU group patient and four FluMel group patients. In the non-SCID group, 10 of 11 patients (91%) who received FluMel achieved complete or high-level mixed chimerism but had variable outcomes. Patients with WAS (2/2 patients), NEMO deficiency (2/2 patients), and X-linked hyper IgM syndrome (2/3 patients) who received FluBU achieved complete or high-level mixed chimerism and long-term survival.

CONCLUSIONS

RIC-HCT with FluBU is a safe and effective strategy for obtaining high-level donor chimerism, immune reconstitution including B cell function, and long-term survival in patients with SCID. In patients with non-SCID PID, the results varied according to the subtype of the disease. Further prospective studies are required to optimize the conditioning regimen for non-SCID PID.

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.

Review of the Pharmacokinetics and Pharmacodynamics of Intravenous Busulfan in Paediatric Patients

We aimed to review the pharmacokinetics (PK) of intravenous busulfan in paediatric patients, identify covariate factors influencing exposure, investigate evidence of changes in PK behaviour over time, and correlate exposure with efficacy and toxicity outcomes. A literature review was undertaken of original research published between 2007 and 2019, investigating the PK and pharmacodynamics (PD) of intravenous busulfan in patients ≤ 18 years of age. The review identified 41 publications characterising the PK, and 45 publications describing the PD, of busulfan. Median typical clearance (CL) was 0.22 L/h/kg and median typical volume of distribution was 0.69 L/kg. Patient weight, age, glutathione-S-transferase A1 (GSTA1) genotype and busulfan dosing day/time were the most commonly identified factors affecting CL. Of nine studies investigating changes in CL, seven reported reduced CL over the 4-day course of treatment. Exposure monitoring methods and therapeutic targets were heterogeneous across studies. Relationships between busulfan exposure and patient outcomes were observed in five studies. One study observed a cumulative area under the concentration-time curve over all days of treatment of between 78 and 101 mg/L·h, and two studies observed an average concentration at first dose of < 600 ng/mL improved overall survival, transplant-related mortality, or relapse. One study observed increased sinusoidal obstructive syndrome with maximum busulfan concentration > 1.88 ng/mL. Patient weight, age and GSTA1 genotype are important covariates to consider when individualising busulfan therapy. Reduced busulfan CL over time may need to be accounted for, particularly in patients not receiving phenytoin co-therapy. Standardised monitoring of busulfan exposure over the entire course of treatment and further investigation of the role of busulfan metabolites and pharmacogenomics is warranted.

Population pharmacokinetics of busulfan in Saudi pediatric patients undergoing hematopoietic stem cell transplantation

Background Busulfan is an antineoplastic drug that is used widely as part of a conditioning regimen in pediatric patients undergoing hematopoietic stem cell transplantation. It has a narrow therapeutic index and highly variable pharmacokinetics; therefore therapeutic drug monitoring is recommended to optimize busulfan dosing. Objective To study the population pharmacokinetics of busulfan in Saudi pediatric patients to optimize its dosing. Settings King Abdullah Specialist Children's Hospital in Riyadh, Saudi Arabia. Methods This pharmacokinetic observational study was conducted between January 2016 and December 2018. All pediatric patients receiving IV busulfan and undergoing routine therapeutic drug monitoring were included. Population pharmacokinetics modeling was conducted using Monolix2019R1. Pharmacokinetic data of busulfan in children. Results The study included 59 patients and 513 samples. The mean ± SD age was 6.10 ± 3.17 years, and the dose administered was 0.994 ± 0.15 mg/kg. The mean ± SD Cmax and area under the curve (AUC) were 900.60 ± 402.8 ng/mL and 1031.14 ± 300.75 µM min, respectively. Based on our simulations, the European Medicines Agency recommended dose were adequate for most patient's groups to achieve the conventional target of an AUC_0-tau of 900-1350 µM min. For patients in the lower weight group < 9 kg, higher doses were need at 1.2 mg/kg. With regards to the newly proposed target of AUC 78-101 mg h/mL, all of the doses we tested had low probability of achieving it. Conclusions Most of our patients had less than a proportional increase in busulfan concentration suggesting autoinduction. The high interindividual variability and autoinduction make dose adjustments challenging and AUC at steady state difficult to predict from the first dose. One approach to improve dose predictions is to use Bayesian dosing software. Based on our simulations, the European Medicines Agency recommended doses were adequate for most patient groups, except those in the lower (< 9 kg) and higher weight groups (> 34 kg).