METHODS: Data from heart transplant recipients (n = 87) administered the oral immediate-release formulation of tacrolimus (Prograf®) were collected. Routine drug monitoring data, principally trough concentrations, were used for model building (n = 1099). A published tacrolimus model was used to inform the estimation of Ka , V2 /F, Q/F and V3 /F. The effect of concomitant azole antifungal use on tacrolimus CL/F was quantified. Fat-free mass was implemented as a covariate on CL/F, V2 /F, V3 /F and Q/F on an allometry scale. Subsequently, stepwise covariate modelling was performed. Significant covariates influencing tacrolimus CL/F were included in the final model. Robustness of the final model was confirmed using prediction-corrected visual predictive check (pcVPC). The final model was externally evaluated for prediction of tacrolimus concentrations of the fourth dosing occasion (n = 87) from one to three prior dosing occasions.