Supporting data for characterization of the busulfan metabolite EdAG and the Glutaredoxins that it adducts

The presence of busulfan metabolites and pharmacometabolomics in plasma drawn immediately before allograft infusion in hematopoietic cell transplant recipients

Busulfan is hepatically metabolized through glutathione (GSH) conjugation; in vitro, this process depletes hepatocyte GSH stores and generates the cytotoxic metabolite γ-glutamyldehydroalanylglycine, which is too unstable to be quantitated in vivo. We sought to evaluate if pre-graft (i.e., immediately before allograft infusion) concentrations of busulfan metabolites' and of endogenous metabolomic compounds (EMCs) representing the glutathione pathway were associated with clinical outcomes in hematopoietic cell transplant (HCT) recipients receiving busulfan. The clinical outcomes evaluated were relapse, acute graft versus host disease (GVHD), chronic GVHD, non-relapse mortality, and neutrophil nadir. In pre-graft samples obtained from patients immediately before allograft infusion, our objectives were to evaluate for: (1) the presence of busulfan and its metabolites tetrahydrothiophenium ion (THT+), tetrahydrothiophene 1-oxide, sulfolane, and 3-hydroxysulfolane (N = 124); (2) EMCs using a global metabolomics assay (N = 77); and (3) the association of the busulfan metabolites and the EMCs with clinical outcomes. In the pre-graft samples, busulfan and THT+ could not be detected. THT 1-oxide, sulfolane, and 3-hydroxysulfolane were quantitated in 9.6%, 26%, and 58% of pre-graft samples; their concentrations were not associated with clinical outcomes. Four pre-graft EMCs were statistically significantly associated with the neutrophil nadir. The pre-graft EMCs were not associated with the other clinical outcomes. In conclusion, busulfan's metabolites are present in patients' plasma immediately before allograft infusion; the neutrophil nadir is associated with pre-graft EMCs. Future research should investigate the association of clinical outcomes with the concentrations of busulfan's metabolites and EMCs in the pre-graft plasma from allogeneic HCT recipients.

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.

Characterization and quantitation of busulfan DNA adducts in the blood of patients receiving busulfan therapy

DNA alkylating drugs have been used as cancer chemotherapy with variable outcomes. The establishment of predictive biomarkers to identify patients who will effectively respond to treatment would allow for the development of personalized therapies. As the degree of interaction of alkylating drug with DNA plays a key role in their mechanism of action, our hypothesis is that the measurement of the DNA adducts formed by alkylating drugs could be used to inform patient stratification. Beginning with busulfan, we took advantage of our DNA adductomic approach to characterize DNA adducts formed by reacting busulfan with calf-thymus DNA. Samples collected from six patients undergoing busulfan-based chemotherapy prior to allogeneic hematopoietic cell transplantation were analyzed for the presence of busulfan-derived DNA adducts. Among the 15 adducts detected in vitro, 12 were observed in the patient blood confirming the presence of a large profile of DNA adducts in vivo. Two of the detected adducts were structurally confirmed by comparison with synthetic standards and quantified in patients. These data confirm our ability to comprehensively characterize busulfan-derived DNA damage and set the stage for the development of methods to support personalized chemotherapy.

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.

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.

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.

The Myeloablative Drug Busulfan Converts Cysteine to Dehydroalanine and Lanthionine in Redoxins

The myeloablative agent busulfan (1,4-butanediol dimethanesulfonate) is an old drug that is used routinely to eliminate cancerous bone marrow prior to hematopoietic stem cell transplant. The myeloablative activity and systemic toxicity of busulfan have been ascribed to its ability to cross-link DNA. In contrast, here we demonstrate that incubation of busulfan with the thiol redox proteins glutaredoxin or thioredoxin at pH 7.4 and 37 °C results in the formation of putative S-tetrahydrothiophenium adducts at their catalytic Cys residues, followed by β-elimination to yield dehydroalanine. Both proteins contain a second Cys, in their catalytic C-X-X-C motif, which reacts with the dehydroalanine, the initial Cys adduct with busulfan, or the S-tetrahydrothiophenium, to form novel intramolecular cross-links. The reactivity of the dehydroalanine (DHA) formed is further demonstrated by adduction with glutathione to yield a lanthionine and by a novel reaction with the reducing agent tris(2-carboxyethyl)phosphine (TCEP), which yields a phosphine adduct via Michael addition to the DHA. Formation of a second quaternary organophosphonium salt via nucleophilic substitution with TCEP on the initial busulfan-protein adduct or on the THT(+)-Redoxin species is also observed. These results reveal a rich potential for reactions of busulfan with proteins in vitro, and likely in vivo. It is striking that several of the chemically altered protein products retain none of the atoms of busulfan, in contrast to typical drug-protein adducts or traditional protein modification reagents. In particular, the ability of a clinically used drug to convert Cys to dehydrolanine in intact proteins, and its subsequent reaction with biological thiols, is unprecedented.

The busulfan metabolite EdAG irreversibly glutathionylates glutaredoxins

The DNA alkylating agent busulfan is used to 'precondition' patients with leukemia, lymphomas and other hematological disorders prior to hematopoietic stem cell transplants. Busulfan is metabolized via conjugation with glutathione (GSH) followed by intramolecular rearrangement to the GSH analog γ-glutamyl-dehydroalanyl -glycine (EdAG). EdAG contains the electrophilic dehydroalanine, which is expected to react with protein nucleophiles, particularly proteins with GSH binding sites such as glutaredoxins (Grx's). Incubation of EdAG with human Grx-1 or Grx-2 results in facile adduction of cys-23 and cys-77, respectively, as determined by ESI-MS/MS. The resulting modified proteins are catalytically inactive. In contrast, the glutathione transferase A1-1 includes a GSH binding site with a potentially reactive tyrosinate (Tyr-9) but it does not react with EdAG. Similarly, Cys-112 of GSTA1-1, which lies outside the active site and is known to form disulfides with GSH, does not react with EdAG. The results provide the first demonstration of the reactivity of any busulfan metabolites with intact proteins, and they suggest that GSH-binding sites containing thiolates are most susceptible. The adduction of Grx's by EdAG suggests the possible alteration of proteins that are normally regulated via Grx-dependent reversible glutathionylation or deglutathionylation. Dysregulation of Grx-dependent processes could contribute to cellular toxicity of busulfan.