Semi-Mechanistic Population Pharmacokinetic Model Incorporating Glutathione S-Transferase Activity for Personalized Busulfan Dosing in Pediatric Allogeneic Hematopoietic Cell Transplantation

Di Cao¹Xiaowen Qian¹Ping Wang¹Xinyi Zheng²Shan Huang³Zhonglin Wei⁴Wenjin Jiang¹Ling Yu¹Xin Jiang¹Ying Yu¹Junjun Mao^2*Xiaowen Zhai^1*

• ¹Children's Hospital, Fudan University, Shanghai, Shanghai Municipality, China
• ²Huashan Hospital, Fudan University, Shanghai, China
• ³Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
• ⁴Children's Hospital of Soochow University, Suzhou, Jiangsu Province, China

Background Busulfan is known for its high inter- and intra-individual pharmacokinetics/pharmacodynamics (PK/PD) variability, especially in children. Therefore, we aimed to identify factors affecting PK variability of busulfan in pediatric allogeneic hematopoietic cell transplantation (HCT) recipients and investigate the effect of glutathione S-transferase (GST) activity on busulfan metabolism using a semi-mechanistic population PK model. Methods Overall, 636 whole-blood busulfan concentrations from 65 pediatric HCT recipients were analyzed using nonlinear mixed-effects modeling. A semi-mechanistic population PK model was developed to describe busulfan metabolism in response to glutathione (GSH) depletion. The effects of potential covariates were selected based on previous study and physiologically-based theoretical mechanisms. Virtual clinical trials were conducted to compare different dosing strategies, and model-based optimal dosing regimen was recommended. Results A two-compartment model with first-order absorption was selected to describe busulfan PK. A GSH compartment was added to represent the relative amount of GSH available at any time. The estimated mean clearance of busulfan was 9.57 L h-1 (relative standard error: 10.8%). Busulfan disposition was best described by including normal fat mass (NFM) allometrically and GST enzyme activity on SGSH exponentially. The SGSH increased by 40.6% as GST enzyme activity increased from 0.9 nmol/min/mL to 20.7 nmol/min/mL. Patients with weights (WT) of 9-16 kg are at high risk of sinusoidal obstructive syndrome (SOS) when receiving WT-based dosing strategy. Conclusions NFM, age-dependent maturation function, and GST enzyme activity may contribute to busulfan PK variability. The WT-based dosing strategy showed a higher risk of SOS than the age-based dosing strategy in 9-16 kg patients.