The effect of glutathione S-transferase variants on pediatric busulfan metabolism was investigated by noncompartmental and population pharmacokinetic modeling. Twenty-nine children who underwent related or unrelated bone marrow or umbilical cord blood hematopoietic cell transplant were retrospectively studied. GSTA1, GSTP1, and GSTM1 variants were explored for their effects on busulfan exposures. Noncompartmental pharmacokinetic analyses showed that carriers of GSTA1*B had a 2.6-fold higher busulfan area under the curve and concentration at steady state compared with noncarriers (P <or= .01). Population pharmacokinetic modeling demonstrated that carriers of GSTA1*B reduced busulfan clearance by 30%. Monte Carlo simulations were then performed to assess busulfan dosing regimens based on GSTA1 genotypes. Simulations determined that dosing based on GSTA1 genotype, weight, and age resulted in fewer children exceeding the upper therapeutic limit compared with dosing using age and weight only. Larger, prospective studies are needed to confirm these findings.