Abstract Busulfan is an alkylating drug routinely used in conditioning regimens for allogeneic hematopoietic cell transplantation (allo-HCT). A myeloablative conditioning regimen, including busulfan, is commonly used in patients undergoing T-cell depletion (TCD) and allo-HCT, but data on optimal busulfan pharmacokinetic (PK) exposure in this setting are limited. Between 2012 and 2019, busulfan PK was performed to target an area under the curve exposure between 55 and 66 mg × h/L over 3 days using a noncompartmental analysis model. We retrospectively re-estimated busulfan exposure following the published population PK (popPK) model (2021) and correlated it with outcomes. To define optimal exposure, univariable models were performed with P splines, wherein hazard ratio (HR) plots were drawn, and thresholds were found graphically as the points at which the confidence interval crossed 1. Cox proportional hazard and competing risk models were used for analyses. 176 patients were included, with a median age of 59 years (range, 2-71). Using the popPK model, the median cumulative busulfan exposure was 63.4 mg × h/L (range, 46.3-90.7). The optimal threshold was at the upper limit of the lowest quartile (59.5 mg × h/L). 5-year overall survival (OS) with busulfan exposure ≥59.5 vs <59.5 mg × h/L was 67% (95% CI, 59-76) vs 40% (95% CI, 53-68), respectively (P = .02), and this association remained in a multivariate analyses (HR, 0.5; 95% CI, 0.29; 0.88; P = .02). In patients undergoing TCD allo-HCT, busulfan exposure is significantly associated with OS. The use of a published popPK model to optimize exposure may significantly improve the OS.