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

Pharmacokinetic Modeling and Simulation with Pharmacogenetic Insights Support the Relevance of Therapeutic Drug Monitoring for Myeloablative Busulfan Dosing in Adult HSCT

Therapeutic drug monitoring (TDM) of busulfan (Bu) is well-established in pediatric hematopoietic stem cell transplantation (HSCT), but its use in adults is limited due to a lack of clear recommendations and scarcity of evidence regarding its utility. GSTA1 promoter variants are reported to affect Bu clearance in both adults and pediatric patients. This study aimed to evaluate the value of preemptive genotyping GSTA1 and body composition (obesity) in individualizing Bu dosing in adults, through pharmacokinetic (PK) modeling and simulations. A population pharmacokinetic (PopPK) model was developed and validated with data from 60 adults who underwent HSCT. Simulations assessed different dosing scenarios based on body size metrics and GSTA1 genotypes. Due to the limited number of obese patients in the cohort, the effect of obesity on Bu pharmacokinetics (PK) was evaluated in silico using a physiologically-based pharmacokinetic (PBPK) model and relevant virtual populations from Simcyp software. Patients with at least 1 GSTA1*B haplotype had 17% lower clearance on average. PopPK simulations indicated that adjusting doses based on genotype increased the probability of achieving the target exposure (3.7 to 5.5 mg.h/L) from 53% to 60 % in GSTA1*A homozygous patients, and from 50% to 61% in *B carriers. Still, Approximately 40% of patients would not achieve this therapeutic window without TDM. A 2-sample optimal design was validated for routine model-based Bu first dose AUC0-∞ estimation, and the model was implemented in the Tucuxi user-friendly TDM software. PBPK simulations confirmed body surface area-based doses of 29 to 31 mg/m2/6h as the most appropriate, regardless of obesity status. This study emphasizes the importance of individualized Bu dosing strategies in adults to achieve therapeutic targets. Preemptive genotyping alone may not have a significant clinical impact, and routine TDM may be necessary for optimal transplantation outcomes.

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.

Population pharmacokinetic model for once‐daily intravenous busulfan in pediatric subjects describing time‐associated clearance

This study aimed to characterize the population pharmacokinetics (PK) of busulfan focusing on how busulfan clearance (CL) changes over time during once‐daily administration and assess different methods for measuring busulfan exposure and the ability to achieve target cumulative exposure under different dosing adjustment scenarios in pediatric stem cell transplantation recipients. Daily serial blood sampling was performed and concentration‐time data were analyzed using a nonlinear mixed‐effects approach. The developed PK model was used to assess achievement of target exposure under six dose‐adjustment scenarios based on simulations performed in RStudio (RxODE package)®. A total of 2491 busulfan plasma concentration–time measurements were collected from 95 patients characterizing 379 dosing days. A two‐compartment model with time‐associated CL best described the data with a typical CL of 14.5 L/h for an adult male with 62 kg normal fat mass (NFM; equivalent to 70 kg total body weight), typical volume of distribution central compartment (V1) of 40.6 L/59 kg NFM (equivalent to 70 kg total body weight), and typical volume of distribution peripheral compartment of 3.57 L/62 kg NFM. Model interindividual variability in CL and V1 was 14.7% and 34.9%, respectively, and interoccasional variability in CL was 6.6%. Patient size described by NFM, a maturation component, and time since start of treatment significantly influenced CL. Simulations demonstrated that using model‐based exposure estimates with each dose, and either a proportional dose‐adjustment calculation or model‐based calculated individual CL estimates to support dose adjustments, increased proportion of subjects attaining cumulative exposure within 5% of target compared with using noncompartmental analysis (100% vs. 0%). A time‐associated reduction in CL during once‐daily busulfan treatment was described.

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

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

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

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

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

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

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

Impact of busulfan pharmacokinetics on outcome in adult patients receiving an allogeneic hematopoietic cell transplantation

Busulfan (Bu) is widely used in conditioning regimens before allogeneic hematopoietic cell transplantation, with variable metabolism due to interindividual differences of pharmacokinetics (PK). The purpose of this study was to correlate pharmacokinetics and clinical outcomes. Lower-AUC, in range-AUC and higher-AUC were defined as ±25% of the targeted Bu-AUC. In 2019, we changed Bu dosing from 4×/day (Bu-4) to 1×/day (Bu-1) for ease of application. AUC-target range was reached in 46% of patients; 40% were in low-AUC and 14% in high-AUC. Among all toxicities, viral and fungal infections were significantly more frequent in high-AUC compared with low-AUC (20% vs. 8%; p = 0.01 and 37% vs. 17%; p = 0.03). Bu-1 showed lower PK values (66% vs. 36% of Bu-4 in low-AUC; p < 0.01) and higher incidence of mucositis (p = 0.02). Long-term outcomes at 2 years showed a higher non-relapse mortality (NRM) (p < 0.01) and higher relative risk of death in the high-AUC group compared to the other groups. Cumulative incidence of relapse and acute/chronic GvHD were not significantly different. The optimal cut-off in Bu-AUC associated with low NRM was 969 µmol/l*min (ROC AUC 0.67, sensitivity 0.86 and specificity 0.47) for Bu-4. In conclusion, low-AUC BU-PK seems of benefit regarding NRM and survival.

Effect of GSTA1 Variants on Busulfan-Based Conditioning Regimen Prior to Allogenic Hematopoietic Stem-Cell Transplantation in Pediatric Asians

Busulfan is widely used as a chemotherapy treatment before hematopoietic stem-cell transplantation (HSCT). However, the response of busulfan is highly variable and unpredictable, whereby the pharmacogenetic interference of glutathione S-transferase (GST) has strong evidence in Caucasians and some adult Asians but not in pediatric Asian patients. This study was aimed at investigating the associations of GST genetic polymorphisms with variations in the pharmacokinetic (PK) properties of busulfan in pediatric Asian patients. This retrospective cohort study recruited 92 pediatric patients. The polymorphism of GSTA1 was genotyped by Sanger sequencing, and GSTM1 and GSTP1 were genotyped by real-time PCR. Drug concentration and PK estimation were identified using an LC-MS/MS method and a noncompartmental model. Statistical analysis was performed by R software. Out of 92 patients, 48 (53%) were males, the mean age was 8.4 ± 5.12 years old, and the average weight was 26.52 ± 14.75 kg. The allele frequencies of GSTA1*B and of GSTM1 and GSTP1* deletions were 16.9%, 68.5%, and 21.2%, respectively. Patients with GSTA1*B had a statistically significant impact on the PK of busulfan, whereas those with GSTM1 and GSTP1 did not (p > 0.05). The carriers of GSTA1*B showed a significant difference compared to noncarriers in terms of t1/2 (for first dose: 161.9 vs. 134.3 min, p = 0.0016; for second dose: 156.1 vs. 129.8, p = 0.012), CL (88.74 vs. 124.23 mL/min, p = 0.0089), Cmax (4232.6 vs. 3675.5 ng/mL, p = 0.0021), and AUC (5310.6 vs. 4177.1 µM/min, p = 0.00033). The augmentation of AUC was around 27.1% in patients carrying the GSTA1*B variant. The GSTA1 polymorphism was significantly associated with variations of the pharmacokinetic properties of busulfan treatment in pediatric Asian patients.

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.

Evaluation of the Robustness of Therapeutic Drug Monitoring Coupled with Bayesian Forecasting of Busulfan with Regard to Inaccurate Documentation

Background

Inaccurate documentation of sampling and infusion times is a potential source of error in personalizing busulfan doses using therapeutic drug monitoring (TDM). Planned times rather than the actual times for sampling and infusion time are often documented.

Therefore, this study aimed to evaluate the robustness of a limited sampling TDM of busulfan with regard to inaccurate documentation.

Methods

A pharmacometric analysis was conducted in NONMEM® 7.4.3 and “R” by performing stochastic simulation and estimation with four, two and one sample(s) per patient on the basis of a one-compartment- (1CMT) and two-compartment (2CMT) population pharmacokinetic model. The dosing regimens consisted of i.v. busulfan (0.8 mg/kg) every 6 h (Q6H) or 3.2 mg/kg every 24 h (Q24H) with a 2 h- and 3 h infusion time, respectively. The relative prediction error (rPE) and relative root-mean-square error (rRmse) were calculated in order to determine the accuracy and precision of the individual AUC estimation.

Results

A noticeable impact on the estimated AUC based on a 1CMT-model was only observed if uncertain documentation reached ± 30 min (1.60% for Q24H and 2.19% for Q6H). Calculated rPEs and rRmse for Q6H indicate a slightly lower level of accuracy and precision when compared to Q24H. Spread of rPE’s and rRmse for the 2CMT-model were wider and higher compared to estimations based on a 1CMT-model.

Conclusions

The estimated AUC was not affected substantially by inaccurate documentation of sampling and infusion time. The calculated rPEs and rRmses of estimated AUC indicate robustness and reliability for TDM of busulfan, even in presence of erroneous records.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11095-021-03115-8.

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.

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

AIM

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

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

Purpose

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

Patients and Methods

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

Results

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

Conclusion

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

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

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

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.

Pharmacogenomics as a Tool to Limit Acute and Long-Term Adverse Effects of Chemotherapeutics: An Update in Pediatric Oncology

In the past decades, new cancer treatments have been introduced in pediatric oncology leading to improvement in clinical outcomes and survival rates. However, due to inter-individual differences, some children experience severe chemotherapy-induced toxicities or a poor clinical outcome. An explanation for the diversity in response to chemotherapy is genetic variation, leading to differences in expression and activity of metabolizing and transport enzymes as well as drug targets. Pharmacogenetic testing has emerged as a promising tool to predict and limit acute and long-term adverse effects in patients. However, in pediatric oncology, limited number of patients and a considerable diversity in study results complicate the interpretation of test results and its clinical relevance. With this review, we provide an overview of new developments over the past four years regarding relevant polymorphisms related to toxicity in pediatric oncology. The following chemotherapeutics and associated toxicities are discussed: alkylating agents, anthracyclines, asparaginase, methotrexate, platinum compounds, steroids, thiopurines, topoisomerase inhibitors, and vinca alkaloids. Our review identifies several questions regarding the role of genetic variants in chemotherapy-induced toxicities. Ambiguities in the literature stem from small population sizes, differences in (statistical) interpretation and variations in sequencing technologies as well as different clinical outcome definitions. Standardization of clinical outcome data and toxicity definitions within electronic health records combined with the increased availability of genomic sequence techniques in clinical practice will help to validate these models in upcoming years.

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

BACKGROUND

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

OBJECTIVES

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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.

Pharmacogenetics: Applications to Pediatric Patients

Individual genomic differences may affect drug disposition and effects of many drugs, and identification of biomarkers are crucial to personalize dosage and optimize response. In children, developmental changes associated with growth and maturation translate into different relationships between genotype and phenotype and different responses to treatment compared to adults. This review aims to summarize some developmental aspects of pharmacogenetics, based on practical examples.

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.