Glutathione S-transferase polymorphisms are not associated with population pharmacokinetic parameters of busulfan in pediatric patients

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.

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.

The Impact of Pharmacogenetics on Pharmacokinetics and Pharmacodynamics in Neonates and Infants: A Systematic Review

In neonates, pharmacogenetics has an additional layer of complexity. This is because in addition to genetic variability in genes that code for proteins relevant to clinical pharmacology, there are rapidly maturational changes in these proteins. Consequently, pharmacotherapy in neonates has unique challenges. To provide a contemporary overview on pharmacogenetics in neonates, we conducted a systematic review to identify, describe and quantify the impact of pharmacogenetics on pharmacokinetics and -dynamics in neonates and infants (PROSPERO, CRD42022302029). The search was performed in Medline, Embase, Web of Science and Cochrane, and was extended by a PubMed search on the ‘top 100 Medicines’ (medicine + newborn/infant + pharmacogen*) prescribed to neonates. Following study selection (including data in infants, PGx related) and quality assessment (Newcastle–Ottawa scale, Joanna Briggs Institute tool), 55/789 records were retained. Retained records relate to metabolizing enzymes involved in phase I [cytochrome P450 (CYP1A2, CYP2A6, CYP2B6, CYP2C8/C9/C18, CYP2C19, CYP2D6, CYP3A5, CYP2E1)], phase II [glutathione-S-transferases, N-acetyl transferases, UDP-glucuronosyl-transferase], transporters [ATP-binding cassette transporters, organic cation transporters], or receptor/post-receptor mechanisms [opioid related receptor and post-receptor mechanisms, tumor necrosis factor, mitogen-activated protein kinase 8, vitamin binding protein diplotypes, corticotrophin-releasing hormone receptor-1, nuclear receptor subfamily-1, vitamin K epoxide reductase complex-1, and angiotensin converting enzyme variants]. Based on the available overview, we conclude that the majority of reported pharmacogenetic studies explore and extrapolate observations already described in older populations. Researchers commonly try to quantify the impact of these polymorphisms in small datasets of neonates or infants. In a next step, pharmacogenetic studies in neonatal life should go beyond confirmation of these associations and explore the impact of pharmacogenetics as a covariate limited to maturation of neonatal life (ie, fetal malformations, breastfeeding or clinical syndromes). The challenge is to identify the specific factors, genetic and non-genetic, that contribute to the best benefit/risk balance.

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.

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.

Clinical pharmacology of cytotoxic drugs in neonates and infants: Providing evidence-based dosing guidance

Cancer in neonates and infants is a rare but challenging entity. Treatment is complicated by marked physiological changes during the first year of life, excess rates of toxicity, mortality, and late effects. Dose optimisation of chemotherapeutics may be an important step to improving outcomes. Body size–based dosing is used for most anticancer drugs used in infants. However, dose regimens are generally not evidence based, and dosing strategies are frequently inconsistent between tumour types and treatment protocols. In this review, we collate available pharmacological evidence supporting dosing regimens in infants for a wide range of cytotoxic drugs. A systematic review was conducted, and available data ranked by a level of evidence (1–5) and a grade of recommendation (A–D) provided on a consensus basis, with recommended dosing approaches indicated as appropriate. For 9 of 29 drugs (busulfan, carboplatin, cyclophosphamide, daunorubicin, etoposide, fludarabine, isotretinoin, melphalan and vincristine), grade A was scored, indicating sufficient pharmacological evidence to recommend a dosing algorithm for infants. For busulfan and carboplatin, sufficient data were available to recommend therapeutic drug monitoring in infants. For eight drugs (actinomycin D, blinatumomab, dinutuximab, doxorubicin, mercaptopurine, pegaspargase, thioguanine and topotecan), some pharmacological evidence was available to guide dosing (graded as B). For the remaining drugs, including commonly used agents such as cisplatin, cytarabine, ifosfamide, and methotrexate, pharmacological evidence for dosing in infants was limited or non-existent: grades C and D were scored for 10 and 2 drugs, respectively. The review provides clinically relevant evidence-based dosing guidance for cytotoxic drugs in neonates and infants.

•

Treating cancer in neonates and infants is challenging.

•

Dose optimisation of cytotoxic drugs is an important step to improving outcomes.

•

Clinical pharmacological evidence supporting dosing regimens in infants was collated.

•

All available pharmacological evidence was ranked by a level of evidence.

•

A grade of recommendation was derived and a recommended dose per agent provided.

Influence of GSTM1, GSTT1, and GSTP1 genetic polymorphisms on disorders in transplant patients: a systematic review

The glutathione-S-transferase (GST) enzymes are phase II isoenzymes responsible for protection against free radicals and xenobiotics. Since these proteins are described as polymorphic, polymorphisms in genes that encode them may alter enzymatic function and contribute to oxidative stress. In this context, such polymorphisms were already associated with several diseases and multiple therapeutic outcomes. A systematic review was performed to evaluate studies regarding the association between polymorphisms in three genes encoding enzymes of the GST family - GSTM1, GSTT1, and GSTP1 - and disorders in transplant patients. A total of 125 articles on which inclusion and exclusion criteria were applied were identified at PubMed database. Thirty-two studies met the target criteria and were included in the review. The mechanisms by which GST genotypes influence the development of disorders in transplant patients differ by disorder: they may participate in it by decreasing metabolism of drugs administered to patients undergoing transplantation, then exposing them to greater toxicity; by decreasing the repair ability against oxidative stress; or by encoding proteins that may be recognized as foreign, setting of an alloimmune reaction. Although some results are better established - such as GSTM1 null genotype's role in the development of toxicity events in transplant patients - others require further evidences, as GST influence on the development of pulmonary decline and posttransplant diabetes mellitus (PTDM). The importance of investigating these associations lies in a personalized medicine, in which the high-risk genotype patient has its treatment individualized and its care for prophylaxis and surveillance increased, potentially reducing this population's morbimortality.

Characterization of Drug-Drug Interactions on the Pharmacokinetic Disposition of Busulfan in Pediatric Patients during Hematopoietic Stem Cell Transplantation Conditioning

AIM

The study objective was to develop a population pharmacokinetic model for busulfan to comprehensively examine drug-drug interactions (DDIs) in pediatric patients undergoing hematopoietic stem cell transplantation. Currently, there is limited evidence to substantiate potential DDIs with busulfan.

METHODS

This retrospective study population was comprised of 250 patients receiving, on average, 0.8 mg/kg intravenous busulfan as pretreatment. All model analyses were conducted using non-linear mixed effects modeling in Pumas v2.0. The metabolic pathways of primary interest were glutathione conjugation and cytochrome P450 activity. Concomitant medications were categorized as CYP inhibitors, inducers, or GST depleters and included in the model as conditional covariates. A bootstrap simulation and visual predictive check were conducted to qualify the final model.

RESULTS

The final one-compartment model incorporates covariates of weight and age in relation to their effects on both total body clearance and volume of distribution. The estimated typical values of clearance and volume were 1.138 L/hr (CI: 1.095 - 1.179 L/hr) and 3.527 L (CI: 3.418 - 3.621 L), respectively. No significant changes in clearance were observed when medications that alter proposed hepatic and metabolic pathways of busulfan were coadministered.

CONCLUSION

To the best of our knowledge, this is the largest single center study of busulfan in children and the first to quantify the maturation effect of both clearance and volume. This study could not demonstrate a difference in busulfan clearance when comparing patients who received medications that alter the GST, CYP3A4, or CYP2C9 pathway to those who did not.

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.

Is Busulfan Clearance Different in Patients With Sickle Cell Disease? Let's Clear Up That Case With Some Controls

In busulfan-based conditioning regimen for hematopoietic stem cell transplantation in children, accurate a priori determination of the first dose is important because of its narrow therapeutic window. Sickle cell disease (SCD) influences pharmacokinetics of the commonly used drugs by affecting organs responsible for drug metabolism and elimination. This pharmacokinetics study assesses the influence of SCD on the metabolic pathway of busulfan that is mainly metabolized in the liver. In this retrospective cross-sectional case-control study, 16 patients with SCD were matched to 50 patients without SCD on known busulfan clearance's covariates (glutathione-S-transferase alpha1 polymorphisms, age, weight). Clearance of the first dose of busulfan was not significantly different independently of genetic or anthropometric factors in patients with or without SCD.

Genetic Predictors for Sinusoidal Obstruction Syndrome-A Systematic Review

Sinusoidal obstruction syndrome (SOS) is a potentially life-threatening complication after hematopoietic stem cell transplantation (HSCT) or antineoplastic treatment without HSCT. Genetic variants were investigated for their association with SOS, but the evidence is inconclusive. We performed a systematic literature review to identify genes, gene variants, and methods of association analyses of genetic markers with SOS. We identified 23 studies after HSCT and 4 studies after antineoplastic treatment without HSCT. One study (4%) performed whole-exome sequencing (WES) and replicated the analysis in an independent cohort, 26 used a candidate-gene approach. Three studies included >200 participants (11%), and six were of high quality (22%). Variants in 34 genes were tested in candidate gene studies after HSCT. Variants in GSTA1 were associated with SOS in three studies, MTHFR in two, and CPS1, CTH, CYP2B6, GSTM1, GSTP1, HFE, and HPSE in one study each. UGT2B10 and LNPK variants were identified in a WES analysis. After exposure to antineoplastic agents without HSCT, variants in six genes were tested and only GSTM1 was associated with SOS. There was a substantial heterogeneity of populations within and between studies. Future research should be based on sufficiently large homogenous samples, adjust for covariates, and replicate findings in independent cohorts.

The analysis of GSTA1 promoter genetic and functional diversity of human populations

GSTA1 encodes a member of a family of enzymes that function to add glutathione to target electrophilic compounds, including carcinogens, therapeutic drugs, environmental toxins, and products of oxidative stress. GSTA1 has several functional SNPs within its promoter region that are responsible for a change in its expression by altering promoter function. This study aims to investigate distributions of GSTA1 promoter haplotypes across different human populations and to assess their impact on the expression of GSTA1. PHASE 2.1.1 was used to infer haplotypes and diplotypes of six GSTA1 promoter SNPs on 2501 individuals from 26 populations classified by the 1000 Genomes Project into five super-populations that included Africa (N = 660), America (N = 347), East Asia (N = 504), Europe (N = 502), and South Asia (N = 488). We used pairwise FST analysis to compare sub-populations and luciferase reporter assay (LRA) to evaluate the impact of each SNP on activation of transcription and interaction with other SNPs. The distributions of GSTA1 promoter haplotypes and diplotypes were significantly different among the different human populations. Three new promoter haplotypes were found in the African super-population. LRA demonstrated that SNPs at -52 and -69 has the most impact on GSTA1 expression, however other SNPs have a significant impact on transcriptional activity. Based on LRA, a new model of cis-elements interaction is presented. Due to the significant differences in GSTA1 diplotype population frequencies, future pharmacogenomics or disease-related studies would benefit from the inclusion of the complete GSTA1 promoter haplotype based on the newly proposed metabolic grouping derived from the LRA results.

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.

Busulfan Pharmacokinetics in Adenosine Deaminase-Deficient Severe Combined Immunodeficiency Gene Therapy

The pharmacokinetics of low-dose busulfan (BU) were investigated as a nonmyeloablative conditioning regimen for autologous gene therapy (GT) in pediatric subjects with adenosine deaminase-deficient severe combined immunodeficiency disease (ADA SCID). In 3 successive clinical trials, which included either γ-retroviral (γ-RV) or lentiviral (LV) vectors, subjects were conditioned with BU using different dosing nomograms. The first cohort received BU doses based on body surface area (BSA), the second cohort received doses based on actual body weight (ABW), and in the third cohort, therapeutic drug monitoring (TDM) was used to target a specific area under the concentration-time curve (AUC). Neither BSA-based nor ABW-based dosing achieved a consistent cumulative BU AUC; in contrast, TDM-based dosing led to more consistent AUC. BU clearance increased as subject age increased from birth to 18 months. However, weight and age alone were insufficient to accurately predict the dose that would consistently achieve a target AUC. Furthermore, various clinical, laboratory, and genetic factors (eg, genotypes for glutathione-S-transferase isozymes known to participate in BU metabolism) were analyzed, but no single finding predicted subjects with rapid versus slow clearance. Analysis of BU AUC and the postengraftment vector copy number (VCN) in granulocytes, a surrogate marker of the level of engrafted gene-modified hematopoietic stem and progenitor cells (HSPCs), demonstrated gene marking at levels sufficient for therapeutic benefit in the subjects who had achieved the target BU AUC. Although many factors determine the ultimate engraftment following GT, this work demonstrates that the BU AUC correlated with the eventual level of engrafted gene-modified HSPCs within a vector group (γ-RV versus LV), with significantly higher levels of granulocyte VCN in the recipients of LV-modified grafts compared to recipients of γ-RV-transduced grafts. Taken together, these findings provide insight into low-dose BU pharmacokinetics in the unique setting of autologous GT for ADA SCID, and these dosing principles may be applied to future GT trials using low-dose BU to open the bone marrow niche.

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.

A semi-mechanistic model based on glutathione depletion to describe intra-individual reduction in busulfan clearance

Aim

To develop a semi‐mechanistic model, based on glutathione depletion and predict a previously identified intra‐individual reduction in busulfan clearance to aid in more precise dosing.

Methods

Busulfan concentration data, measured as part of regular care for allogeneic hematopoietic cell transplantation (HCT) patients, were used to develop a semi‐mechanistic model and compare it to a previously developed empirical model. The latter included an empirically estimated time effect, where the semi‐mechanistic model included theoretical glutathione depletion. As older age has been related to lower glutathione levels, this was tested as a covariate in the semi‐mechanistic model. Lastly, a therapeutic drug monitoring (TDM) simulation was performed comparing the two models in target attainment.

Results

In both models, a similar clearance decrease of 7% (range −82% to 44%), with a proportionality to busulfan metabolism, was found. After 40 years of age, the time effect increased with 4% per year of age (0.6–8%, P = 0.009), causing the effect to increase more than a 2‐fold over the observed age‐range (0–73 years). Compared to the empirical model, the final semi‐mechanistic model increased target attainment from 74% to 76%, mainly through better predictions for adult patients.

Conclusion

These results suggest that the time‐dependent decrease in busulfan clearance may be related to gluthathione depletion. This effect increased with older age (>40 years) and was proportional to busulfan metabolism. The newly constructed semi‐mechanistic model could be used to further improve TDM‐guided exposure target attainment of busulfan in patients undergoing HCT.

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).

Population pharmacokinetic analysis of intravenous busulfan: GSTA1 genotype is not a predictive factor of initial dose in Chinese adult patients undergoing hematopoietic stem cell transplantation

PURPOSE

This study aimed to develop a population pharmacokinetic (PPK) model to investigate the impact of GSTA1, GSTP1, and GSTM1 genotypes on busulfan pharmacokinetic (PK) variability in Chinese adult patients.

METHODS

Forty-three and 19 adult patients who underwent hematopoietic stem cell transplantation (HSCT) were enrolled for modeling group and validation group, respectively. All patients received twice-daily intravenous busulfan as part of conditioning regimen before HSCT. The PPK model was developed by nonlinear mixed-effect modeling. Covariates investigated were age, sex, actual body weight, body surface area, diagnoses, hepatic function markers, GST genotypes and conditioning regimen.

RESULTS

A total of 488 busulfan concentrations from 43 patients were obtained for the PPK model. The PK of intravenous busulfan was described by one-compartment model with first-order elimination with estimated clearance (CL) of 14.2 L/h and volume of distribution of 64.1 L. Inclusion of GSTA1 genotype as a covariate accounted for 1.1% of the inter-individual variability of busulfan CL (from 17.8% in the basic model to 16.7% in the final model). The accuracy and applicability of the final model were externally validated in the independent group. The difference of busulfan PK between Chinese patients and Caucasian patients existed because of the rarity of haplotype *B in Chinese population.

CONCLUSIONS

Although the GSTA1 genotype-based PPK model of intravenous busulfan was successfully developed and externally validated, the GSTA1 genotype was not considered to be clinically relevant to busulfan CL. We did not suggest the guidance of GSTA1 genotype on initial busulfan dose in Chinese adult patients.

Impact of GSTA1 Polymorphisms on Busulfan Oral Clearance in Adult Patients Undergoing Hematopoietic Stem Cell Transplantation

Background: Busulfan pharmacokinetics exhibit large inter-subject variability. Our objective was to evaluate the influence of glutathione S-transferase A1 (GSTA1) gene variants on busulfan oral clearance (CLo) in a population of patients undergoing hematopoietic stem cell transplantation. Methods: This is a quasi-experimental retrospective study in adult patients (n = 87 included in the final analyses) receiving oral busulfan. Pharmacokinetics data (area under the plasma concentration-time curve (AUC) determined from 10 blood samples) were retrieved from patients’ files and GSTA1 *A and *B allele polymorphisms determined from banked DNA samples. Three different limited sampling methods (LSM) using four blood samples were also compared. Results: Carriers of GSTA1*B exhibited lower busulfan CLo than patients with an *A/*A genotype (p < 0.002): Busulfan CLo was 166 ± 31, 187 ± 37 vs. 207 ± 47 mL/min for GSTA1*B/*B,*A/*B and *A/*A genotypes, respectively. Similar results were obtained with the tested LSMs. Using the standard AUC method, distribution of patients above the therapeutic range after the first dose was 29% for GSTA1*A/*A, 50% for *A/*B, and 65% for *B/*B. The LSMs correctly identified ≥91% of patients with an AUC above the therapeutic range. The misclassified patients had a mean difference less than 5% in their AUCs. Conclusion: Patients carrying GSTA1 loss of function *B allele were at increased risk of overdosing on their initial busulfan oral dose. Genetic polymorphisms associated with GSTA1 explain a significant part of busulfan CLo variability which could be captured by LSM strategies.

Personalized pharmacokinetic targeting with busulfan in allogeneic hematopoietic stem cell transplantation in infants with acute lymphoblastic leukemia

Individual busulfan (BU) dosing based on pharmacokinetic (PK) data is preferable for hematopoietic stem cell transplantation (HSCT) conditioning, but information on BU PK in infants is scarce. We report BU PK data on HSCT conditioning for infants with KMT2A-gene-rearrangement-positive acute lymphoblastic leukemia (MLL-r ALL). Infants showed wide variation in BU PK indices, such as clearance (CL) and volume of distribution (V_d) value, which are distributed more widely among those who received oral, rather than intravenous (IV), BU. Because the steady state concentration (C_ss) fluctuates readily in infants, dose re-adjustment based on PK at the initial administration was important even if the initial dose was determined by a PK test. HSCT can be performed safely within the C_ss range of 600-900 ng/mL per dose, although it was difficult to fit within the therapeutic index of BU. The clinical outcome of engraftment, graft-versus-host disease, adverse events, including sinusoidal obstruction syndrome, and survival did not correlate with the BU PK data, which paradoxically suggests that remaining within this C_ss range helped minimize transplant-related toxicities, while securing engraftment in infants with MLL-r ALL.

Effect of glutathione S-transferase genetic polymorphisms on busulfan pharmacokinetics and veno-occlusive disease in hematopoietic stem cell transplantation: A meta-analysis

This meta-analysis was conducted to derive an integrated conclusion about the influence of glutathione S-transferase (GST) genetic polymorphisms on busulfan pharmacokinetic (PK) parameters and veno-occlusive disease (VOD). Studies which analysed the effect of GST genetic polymorphisms on area under the curve (AUC), clearance (CL) or VOD were searched for and selected. A pooled analysis was conducted using Comprehensive Meta-Analysis programme. Nineteen studies were included in this meta-analysis. GSTA1*B and GSTM1 null genotypes significantly decreased CL_IV of busulfan (standardized difference in means (SDM) = -1.103; P = 0.019 and SDM = -0.418; P = 0.002, respectively). GSTA1*B significantly increased AUC_IV of busulfan (SDM = 0.832; P = 0.046), whereas GSTM1 did not (SDM = 0.155; P = 0.478). The PK parameters of oral busulfan did not differ according to GST genotype. GSTA1, GSTM1 and GSTP1 were not significantly associated with VOD occurrence. GSTA1 and GSTM1 genotypes affected CL_IV of busulfan, but only GSTA1 affected AUC_IV . There was no significant difference in the PK parameters of oral busulfan (CL_PO and AUC_PO ) and VOD when only GST genotypes were considered.

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.

GSTA1 diplotypes affect busulfan clearance and toxicity in children undergoing allogeneic hematopoietic stem cell transplantation: a multicenter study

Busulfan (BU) dose adjustment following therapeutic drug monitoring contributes to better outcome of hematopoietic stem cell transplantation (HSCT). Further improvement could be achieved through genotype-guided BU dose adjustments. To investigate this aspect, polymorphism within glutathione S transferase genes were assessed. Particularly, promoter haplotypes of the glutathione S transferase A1 (GSTA1) were evaluated in vitro, with reporter gene assays and clinically, in a pediatric multi-center study (N =138) through association with BU pharmacokinetics (PK) and clinical outcomes. Promoter activity significantly differed between the GSTA1 haplotypes (p<0.001) supporting their importance in capturing PK variability. Four GSTA1 diplotype groups that significantly correlated with clearance (p=0.009) were distinguished. Diplotypes underlying fast and slow metabolizing capacity showed higher and lower BU clearance (ml/min/kg), respectively. GSTA1 diplotypes with slow metabolizing capacity were associated with higher incidence of sinusoidal obstruction syndrome, acute graft versus host disease and combined treatment-related toxicity (p<0.0005). Among other GST genes investigated, GSTP1 313GG correlated with acute graft versus host disease grade 1-4 (p=0.01) and GSTM1 non-null genotype was associated with hemorrhagic cystitis (p=0.003). This study further strengthens the hypothesis that GST diplotypes/genotypes could be incorporated into already existing population pharmacokinetic models for improving first BU dose prediction and HSCT outcomes. (N^o Clinicaltrials.gov identifier: NCT01257854. Registered 8 December 2010, retrospectively registered).

Clarifying busulfan metabolism and drug interactions to support new therapeutic drug monitoring strategies: a comprehensive review

INTRODUCTION

Busulfan (Bu) is an alkylating agent with a limited therapeutic margin and exhibits inter-patient variability in pharmacokinetics (PK). Despite decades of use, mechanisms of Bu PK-based drug-drug interactions (DDIs), as well as the negative downstream effects of these DDIs, have not been fully characterized. Areas covered: This article provides an overview of Bu PK, with a primary focus on how known and potentially unknown drug metabolism pathways influence Bu-associated DDIs. In addition, pharmacogenomics of Bu chemotherapy and Bu-related DDIs observed in the stem cell transplant clinic (SCT) are summarized. Finally the increasing importance of Bu therapeutic drug monitoring is highlighted. Expert opinion: Mechanistic studies of Bu metabolism have shown that in addition to GST isoenzymes, other oxidative enzymes (CYP, FMO) and ABC/MDR drug transporters likely contribute to the overall clearance of Bu. Despite many insights, results from clinical studies, especially in polypharmacy settings and between pediatric and adult patients, remain conflicting. Further basic science and clinical investigative efforts are required to fully understand the key factors determining Bu PK characteristics and its effects on complications after SCT. Improved TDM strategies are promising components to further investigate, for instance DDI mechanisms and patient outcomes, in the highly complex SCT treatment setting.

Pharmacometabolomics for predicting variable busulfan exposure in paediatric haematopoietic stem cell transplantation patients

Owing to its narrow therapeutic range and high pharmacokinetic variability, optimal dosing for busulfan is important to minimise overexposure-related systemic toxicity and underexposure-related graft failure. Using global metabolomics, we investigated biomarkers for predicting busulfan exposure. We analysed urine samples obtained before busulfan administration from 59 paediatric patients divided into 3 groups classified by area under the busulfan concentration-time curve (AUC), i.e., low-, medium-, and high-AUC groups. In the high-AUC group, deferoxamine metabolites were detected. Phenylacetylglutamine and two acylcarnitines were significantly lower in the high-AUC group than in the low-AUC group. Deferoxamine, an iron-chelating agent that lowers serum ferritin levels, was detected in the high-AUC group, indicating that those patients had high ferritin levels. Therefore, in a retrospective study of 130 paediatric patients, we confirmed our hypothesis that busulfan clearance (dose/AUC) and serum ferritin level has a negative correlation (r = −0.205, P = 0.019). Ferritin, acylcarnitine, and phenylacetylglutamine are associated with liver damage, including free radical formation, deregulation of hepatic mitochondrial β-oxidation, and hyperammonaemia. Our findings reveal potential biomarkers predictive of busulfan exposure and suggest that liver function may affect busulfan exposure.

Energy landscape of a GSTP1 polymorph linked with cytological function decay in response to chemical stressors

Gene polymorphisms lead to varied structure and functional properties. A single nucleotide polymorphism (SNP) i.e. Ile105Val (rs1695) in glutathione S-transferase P1 (GSTP1) gene influences cytological toxicity and modulates the risk to occupational diseases. Apart from this, cancer, neuropathy, NOx, SOx and ozone mediated respiratory function decline including lung inflammation, asthma, allergy etc., have been reported in people with this missense mutation. Here, the functional properties of rs1695 polymorph are revisited through a computational approach. Changes incurred by GSTP1 antioxidant protein as a result of alteration in its sequence, have been studied through docking followed by Poisson-Boltzmann electrostatic equation interpretation, grid and coulombic energy profile mapping for protein polymorphs with DelPhi. Molecular docking simulation of variant and wild type (WT) protein was carried out with eight FDA approved compounds that target GSTP1 for treatment of various diseases. This was to observe binding pattern variation upon mutation induction. Grid, reaction field and coulombic energy calculation of WT and mutated polymorph, complexed with and without these moieties was then attempted. Alteration in conformation and energy was observed in apo- and holo- form of GSTP1 and their ligand-bound complexes as a result of this mutation. This study is a demo of appraising gene-environment interaction based deleteriousness through molecular docking and dynamics simulation approach.

Prediction of Intravenous Busulfan Clearance by Endogenous Plasma Biomarkers Using Global Pharmacometabolomics

INTRODUCTION

High-dose busulfan (busulfan) is an integral part of the majority of hematopoietic cell transplantation conditioning regimens. Intravenous (IV) busulfan doses are personalized using pharmacokinetics (PK)-based dosing where the patient's IV busulfan clearance is calculated after the first dose and is used to personalize subsequent doses to a target plasma exposure. PK-guided dosing has improved patient outcomes and is clinically accepted but highly resource intensive.

OBJECTIVE

We sought to discover endogenous plasma biomarkers predictive of IV busulfan clearance using a global pharmacometabolomics-based approach.

METHODS

Using LC-QTOF, we analyzed 59 (discovery) and 88 (validation) plasma samples obtained before IV busulfan administration.

RESULTS

In the discovery dataset, we evaluated the association of the relative abundance of 1885 ions with IV busulfan clearance and found 21 ions that were associated with IV busulfan clearance tertiles (r² ≥ 0.3). Identified compounds were deoxycholic acid and/or chenodeoxycholic acid, and linoleic acid. We used these 21 ions to develop a parsimonious seven-ion linear predictive model that accurately predicted IV busulfan clearance in 93% (discovery) and 78% (validation) of samples.

CONCLUSION

IV busulfan clearance was significantly correlated with the relative abundance of 21 ions, seven of which were included in a predictive model that accurately predicted IV busulfan clearance in the majority of the validation samples. These results reinforce the potential of pharmacometabolomics as a critical tool in personalized medicine, with the potential to improve the personalized dosing of drugs with a narrow therapeutic index such as busulfan.

Pharmacometabonomic Prediction of Busulfan Clearance in Hematopoetic Cell Transplant Recipients

Intravenous (IV) busulfan doses are often personalized to a concentration at steady state (Css) using the patient's clearance, which is estimated with therapeutic drug monitoring. We sought to identify biomarkers of IV busulfan clearance using a targeted pharmacometabonomics approach. A total of 200 metabolites were quantitated in 106 plasma samples, each obtained before IV busulfan administration in hematopoietic cell transplant (HCT) recipients. Both univariate linear regression with false discovery rate (FDR) and pathway enrichment analyses using the Global test were performed. In the univariate analysis, glycine, N-acetylglycine, 2-hydroxyisovaleric acid, creatine, serine, and tyrosine were statistically significantly associated with IV busulfan clearance at P < 0.05, with the first three satisfying the FDR of q < 0.1. Using pathway enrichment analysis, the glycine, serine, and threonine metabolism pathway was statistically significantly associated with IV busulfan clearance at P < 0.05 and q < 0.1, and a pathway impact >0.1. Glycine is a component of glutathione, which is conjugated with busulfan via glutathione transferase enzymes. These results demonstrate the potential utility of pharmacometabonomics to inform IV busulfan dosing. Future studies are required to validate these findings.

Glutathione Transferase Gene Variants Influence Busulfan Pharmacokinetics and Outcome After Myeloablative Conditioning

Supplemental Digital Content is Available in the Text.

Background:

Busulfan (Bu) and cyclophosphamide (Cy) are frequently included in conditioning regimens before hematopoietic stem cell transplantation (HSCT). Both drugs are detoxified by glutathione transferases (GST), and GST gene variants may explain some of the interindividual variability in pharmacokinetics and drug toxicity.

Methods:

The study investigated adult patients (n = 114) receiving oral Bu pre-HSCT. Bu doses were adjusted to obtain an average steady-state concentration (C_ss) of 900 mcg/L.

Results:

Median first dose Bu C_ss was 1000 mcg/L (600–1780 mcg/L). Patients carrying 1 and 2 GSTA1*B (rs3957357) alleles demonstrated median 12% and 16% higher Bu C_ss (P ≤ 0.05). Bu exposure (average C_ss; odds ratio = 1.009, 95% confidence interval = 1.002–1.017, P = 0.013) and GSTM1 gene copy number (odds ratio = 17.1, 95% confidence interval = 1.46–201, P = 0.024) were significant predictors of mortality ≤30 days. The mortality was 25% versus 2% among carriers of 2 versus no GSTM1 copies (P = 0.021). Mortality ≤3 months was associated with higher first dose Bu exposure (1090 versus 980 mcg/L, P = 0.021). GSTM1 expression and high Bu exposure may increase Cy toxicity by reducing intracellular glutathione.

Conclusions:

GST genotyping before HSCT may allow better prediction of Bu pharmacokinetics and drug toxicity, and thereby improve outcome after BuCy conditioning.

Influence of glutathione S-transferase gene polymorphisms on busulfan pharmacokinetics and outcome of hematopoietic stem-cell transplantation in thalassemia pediatric patients

Hematopoietic stem-cell transplantation (HSCT) is currently the only curative therapeutic option for the treatment of thalassemia. In spite of the high cure rate, HSCT can lead to life-threatening adverse events in some patients. Busulfan (Bu) is a key component of the conditioning regimen prior to HSCT. Inter-individual differences in Bu pharmacokinetics (PK) are hypothesized to influence Bu efficacy and toxicity. Since Bu is mainly metabolized by glutathione S-transferase (GST), we investigated the relationship of GSTA1 and GSTM1 genotypes with first-dose PK and HSCT outcomes in 44 children with thalassemia intermedia and thalassemia major. All children received a myeloablative conditioning regimen with IV Bu. Association analysis revealed a relationship between GSTA169C>T (or haplotype *A/*B) and first Bu dose PK that was dependent on sex and Pesaro risk classification (PRC). Among female patients and patients with PRC I-II, homozygous individuals for the GSTA1T-69 allele defining haplotype *B, had higher Bu exposure and lower clearance (P⩽0.01). Association with HSCT outcomes showed that patients with the GSTM1 null genotypes had higher occurrence of regimen-related toxicity (P=0.01). These results suggest that GST genotypes could be useful to tailor the first Bu dose accordingly to improve HSCT outcome.

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.

Gonadotoxic effects of busulfan in two strains of mice

Busulfan is a chemotherapy drug that has side effects on spermatogonial stem cells (SSC). The effects of bulsufan treatment on male germ cells and fertility vary significantly between individuals. In this study, we have used molecular, cellular and histopathology approaches to investigate the effects of a single intraperitoneal dose of busulfan (40mgkg(-1)) in two mice strains, Balb/C and Swiss, at two different periods after treatment, 30 and 90 days. Testicular degeneration was observed in both Balb/C and Swiss mice after busulfan injection. Interestingly, testicular functions and fertility recovered spontaneously post busulfan treatment in Swiss mice, but not in Balb/C mice. Abnormal fertility induced by busulfan in Balb/C mice was associated with altered seminiferous tubules, sperm morphology and transcript levels of Nanos2, Nanos3, Gdnf and Plzf genes. These findings revealed that SSC of Balb/C mice are more sensitive to the toxic effects of busulfan then those of Swiss mice.

Population pharmacokinetics and pharmacodynamics of busulfan with GSTA1 polymorphisms in patients undergoing allogeneic hematopoietic stem cell transplantation

AIM

A population pharmacokinetic (PPK) analysis was conducted to describe the influence of GSTA1 polymorphisms on intravenous busulfan in adults undergoing allogeneic hematopoietic stem cell transplantation.

PATIENTS & METHODS

A PPK model was developed from 36 patients by a one-compartment model with first-order elimination.

RESULTS

The typical value of clearance and volume of distribution were 11.0 l/h and 42.4 l, respectively. Clearance decreased by 15% and area under the concentration-time curves (AUCs) increased with GSTA1 variants compared with wild-type (both p < 0.05). Subtherapeutic AUCs were seen only in wild-type patients.

CONCLUSION

To our knowledge, this is the first PPK study to suggest that GSTA1 polymorphisms in adults are associated with busulfan PK.

Role of Pharmacogenetics in Hematopoietic Stem Cell Transplantation Outcome in Children

Hematopoietic stem cell transplantation (HSCT) is an established therapeutic procedure for several congenital and acquired disorders, both malignant and nonmalignant. Despite the great improvements in HSCT clinical practices over the last few decades, complications, such as graft vs. host disease (GVHD) and sinusoidal obstructive syndrome (SOS), are still largely unpredictable and remain the major causes of morbidity and mortality. Both donor and patient genetic background might influence the success of bone marrow transplantation and could at least partially explain the inter-individual variability in HSCT outcome. This review summarizes some of the recent studies on candidate gene polymorphisms in HSCT, with particular reference to pediatric cohorts. The interest is especially focused on pharmacogenetic variants affecting myeloablative and immunosuppressive drugs, although genetic traits involved in SOS susceptibility and transplant-related mortality are also reviewed.

New insights into the pharmacokinetics of intravenous busulfan in children with sickle cell anemia undergoing bone marrow transplantation

BACKGROUND

Busulfan (Bu) is an integral part of conditioning regimens for patients with sickle cell anemia (SCA) undergoing transplantation. Patients with SCA might predispose to transplant-related neurological and pulmonary toxicities due to pre-existing disease-related cerebrovascular and lung injury. Bu therapy appears to be an important contributing factor in this context.

PROCEDURE

We studied the pharmacokinetics of intravenous Bu and clinical outcomes of 36 children with SCA undergoing bone marrow transplantation. Most patients had pre-existing organ system damage. Busulfan was administered every 6 hr for 4 days with pharmacokinetic-guided dose adjustment to target a conservative area under the concentration versus time curve (AUC) range of 900-1,350 µMol*min.

RESULTS

We found that the first-dose Bu clearance was significantly higher (P < 0.0005) than the subsequent daily clearance, which remained unchanged during the following days. After the first-dose, 69% of patients achieved the target range. We adapted a new dose-adjustment strategy targeting exposures to the lower end (900 µMol*min) of the AUC range after the first dose of Bu to avoid unnecessary dose increases on subsequent days due to differences in clearance. This strategy enabled most patients to maintain the AUC within therapeutic range following dose adjustments.

CONCLUSIONS

Differences in Bu clearance after the first-dose and subsequent daily doses in patients with SCA should be considered for pharmacokinetic-guided dose adjustment. Conservative AUC range and targeting exposures to the lower end of the range after the first dose was associated with negligible toxicity, and high engraftment and sickle cell-free survival rates.

An analysis of the therapeutic benefits of genotyping in pediatric hematopoietic stem cell transplantation

Hematopoietic stem cell transplantation is a high-risk procedure that is offered, with curative intent, to patients with malignant and nonmalignant disease. The clinical benefits of personalization of therapy by genotyping have been demonstrated by the reduction in transplant related mortality from donor-recipient HLA matching. However, defining the relationship between genotype and transplant conditioning agents is yet to be translated into clinical practice. A number of the therapeutic agents used in stem cell transplant preparative regimens have pharmacokinetic parameters that predict benefit of incorporating pharmacogenomic data into dosing strategies. Busulfan, cyclophosphamide, thio-TEPA and etoposide have well-described drug metabolism pathways, however candidate gene studies have identified there is a gap in the identification of pharmacogenomic data that can be used to improve transplant outcomes. Incorporating pharmacogenomics into pharmacokinetic modeling may demonstrate the therapeutic benefits of genotyping in transplant preparative regimen agents.

Predictive performance of a physiologically based pharmacokinetic model of busulfan in children

A physiologically based pharmacokinetic (PBPK) model of the DNA-alkylating agent busulfan was slightly modified and scaled from adults to children in order to predict the systemic busulfan drug exposure in children. Capitalizing on the recent major software release of PK-Sim®, we refined our PBPK model by implementing glutathione S transferase (GST) in 11 organs using the software integrated enzyme expression database. In addition, two irreversible binding processes (i.e., DNA and plasma protein binding) were applied by using Koff and KD values. The model was scaled from adults to children. Simulations were computed and compared to concentration-time data after intravenous (i.v.) busulfan administration to 36 children. Based on the results, an age-dependent enzyme activity and maturation ratio was tailored and evaluated with an external dataset consisting of 23 children. Initial adult to children scaling indicated lower clearance values for children in comparison to adults. Subsequent age-dependent maturation ratio resulted in three different age groups: Activity of busulfan-glutathione conjugate formation was 80%, 61%, and 89% in comparison to adults for children with an age of up to 2 years, > 2-6 years, and > 6-18 years, respectively. Patients of the evaluation dataset were simulated with a mean percentage error (MPE) for all patients of 3.9% with 3/23 children demonstrating a MPE of > ±30%. The PBPK model parameterization sufficiently described the observed concentration-time data of the validation dataset while showing an adequate predictive performance. This PBPK model could be helpful to determine the first dose of busulfan in children.

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.

The association of cytochrome P450 genetic polymorphisms with sulfolane formation and the efficacy of a busulfan-based conditioning regimen in pediatric patients undergoing hematopoietic stem cell transplantation

Cytochrome P450 enzymes (CYPs) and flavin-containing monooxygenases (FMOs) likely have a role in the oxidation of intermediate metabolites of busulfan (Bu). In vitro studies to investigate the involvement of these enzymes are cumbersome because of the volatile nature of the intermediate metabolite tetrahydrothiophene (THT) and the lack of sensitive quantitation methods. This study explored the association between the CYP2C9, CYP2C19, CYP2B6 and FMO3 genotypes and sulfolane (Su, a water soluble metabolite of Bu) plasma levels in children undergoing hematopoietic stem cell transplantation (HSCT). The relationship between these genotypes and the effectiveness of myeloablative conditioning was also analyzed. Sixty-six children receiving an intravenous Bu-based myeloablative conditioning regimen were genotyped for common functional variant alleles in CYP2C9 (*2 and *3), CYP2C19 (*2 and *17), FMO3 (rs2266780, rs2266782 and rs1736557) and CYP2B6 (*5 and *9). The plasma levels of Bu and its metabolite Su were measured after the ninth Bu dose in a subset of 44 patients for whom plasma samples were available. The ratio of Bu to Su was considered the metabolic ratio (MR) and was compared across the genotype groups. Higher MRs were observed in CYP2C9*2 and *3 allele carriers (mean±s.d.: 7.8±3.6 in carriers vs 4.4±2.2 in non-carriers; P=0.003). An increased incidence of graft failure was observed among patients with an MR>5 compared with those with MR values <5 (20% vs 0%; P=0.02). In contrast, a significantly higher incidence of relapse and graft failure (evaluated as event-free survival) was observed in patients with malignant disease who carried CYP2B6 alleles with reduced function on both chromosomes compared with carriers of at least one normal allele (100% vs 40%; P=0.0001). These results suggest that CYP2C9 has a role in the oxidation reactions of THT and indicate that it may be possible to predict the efficacy of Bu-based myeloablative conditioning before HSCT on the basis of CYP genotypes and Bu MRs.