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

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.

Pharmacogenetic Aspects of Drug Metabolizing Enzymes and Transporters in Pediatric Medicine: Study Progress, Clinical Practice and Future Perspectives

As the activity of certain drug metabolizing enzymes or transporter proteins can vary with age, the effect of ontogenetic and genetic variation on the activity of these enzymes is critical for the accurate prediction of treatment outcomes and toxicity in children. This makes pharmacogenetic research in pediatrics particularly important and urgently needed, but also challenging. This review summarizes pharmacogenetic studies on the effects of genetic polymorphisms on pharmacokinetic parameters and clinical outcomes in pediatric populations for certain drugs, which are commonly prescribed by clinicians across multiple therapeutic areas in a general hospital, organized from those with the most to the least pediatric evidence among each drug category. We also further discuss the research status of the gene-guided dosing regimens and clinical implementation of pediatric pharmacogenetics. More and more drug-gene interactions are demonstrated to have clinical validity for children, and pharmacogenomics in pediatrics have shown evidence-based benefits to enhance the efficacy and precision of existing drug dosing regimens in several therapeutic areas. However, the most important limitation to the implementation is the lack of high-quality, rigorous pediatric prospective clinical studies, so adequately powered interventional clinical trials that support incorporation of pharmacogenetics into the care of children are still needed.

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.

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.

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.