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

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

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

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

PURPOSE

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

PATIENTS AND METHODS

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

RESULTS

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

CONCLUSION

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

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

INTRODUCTION

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

AREAS COVERED

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

EXPERT OPINION

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

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

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

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

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

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

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

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

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.

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.

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.

New dosing nomogram and population pharmacokinetic model for young and very young children receiving busulfan for hematopoietic stem cell transplantation conditioning

BACKGROUND

Busulfan (Bu) is the cornerstone of conditioning regimens prior to hematopoietic stem cell transplantation, widely used in both adults and children for the treatment of malignant and nonmalignant diseases. Despite an intravenous formulation, interindividual variability (IIV) remains high and optimal exposure difficult to achieve, especially in neonates and infants.

PROCEDURE

To ensure both efficacy and safety, we set up in 2005 an observational study designed for children not fully assessed during the drug registration procedure. From a large cohort of 540 patients, we developed a Bu population pharmacokinetic model based on body weight (BW) and maturation concepts to reduce IIV and optimize exposure. A new dosing nomogram was evaluated to better fit the population pharmacokinetic model.

RESULTS

Bu clearance IIV was significantly decreased from 61.3% (covariate-free model) to 28.6% when combining BW and maturation function. Median Bu area under the curve (AUC) was 1179 µmol/L × min compared to 1025 with the EMA dosing nomogram for children <9 kg. The target AUC was reached for each BW strata, significantly increasing the percentages of patients achieving reaching the targeted AUC as compared to FDA schedule.

CONCLUSION

This new model made it possible to propose a novel dosing nomogram that better considered children below 16 kg of BW and allowed better initial exposure as compared to existing dosing schedules. This nomogram, which would be easy to use to determine an optimal dosing schedule in daily practice, will need to be validated in clinical routine. Therapeutic drug monitoring remains strongly advisable for small children and those with specific diseases.

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.

Therapeutic Drug Monitoring of Busulfan for the Management of Pediatric Patients: Cross-Validation of Methods and Long-Term Performance

BACKGROUND

Busulfan (Bu) is an alkylating agent used as part of the conditioning regimen in pediatric patients before hematopoietic stem cell transplantation. Despite intravenous (IV) administration and dosing recommendations based on age and weight, reports have revealed interindividual variability in Bu pharmacokinetics and the outcomes of hematopoietic stem cell transplantation. In this context, adjusting doses to Bu's narrow therapeutic window is advised. We aimed to assess the utility of therapeutic drug monitoring (TDM) of Bu in children, the reliability of Bu quantification methods, and its stability in plasma when stored for up to 5 years.

METHODS

Eighteen patients from our TDM center (252 samples) were included. All of them received a 2-hour Bu IV infusion 4 times daily for a total of 16 doses. The first dose of Bu was age/weight-based, and the subsequent doses were adjusted from third or fifth dose onward based on the estimated first dose pharmacokinetic parameters to target steady-state concentrations (Css) of 600-900 ng/mL. The performance of our unit's high-performance liquid chromatography with tandem mass spectrometry method was assessed using a quality control (QC, 35 series) chart. International, multicenter, cross-validation test (n = 21) was conducted to validate different analytical methods. To assess Bu stability, regression analyses and Bland-Altman plots were performed on measurements at repeated time points on samples stored at -80°C for up to 5 years.

RESULTS

We observed a 4.2-fold interindividual variability in Bu Css after the first dose, with only 28% of children having a Css within the target range. During the 4 days of conditioning, 83% of children had their doses modified according to TDM recommendations. This achieved a Css within the target range in 75% of the children. Routine QC measurements were generally within the ±15% range around theoretical values, showing the optimal robustness of our center's analytical method. Two of the 21 Bu TDM centers returned inadequate results during cross-validation testing; both used a UV detection method. Storage at -80°C led to a fall in Bu content of 14.9% ± 13.4% at 2-4 years and of 20% ± 5% by 5 years (roverall = 0.92).

CONCLUSIONS

We conclude that TDM is an effective method of achieving targeted Bu levels in children. QC programs are crucial to monitoring and maintaining the quality of an analytical method.

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.