Ezetimibe is a potent cholesterol absorption inhibitor, with an erratic pharmacokinetic (PK) profile, attributed to an extensive enterohepatic recirculation (EHC). The aim of this study was to develop a population PK model able to adequately characterize the complex distribution processes of total ezetimibe. The analysis was performed on the individual concentration-time data obtained from 28 healthy subjects who participated in a bioequivalence study comparing two oral ezetimibe formulations. The population PK analysis was performed using nonlinear mixed effect modeling, where different EHC models were developed and evaluated for their performance. Total ezetimibe pharmacokinetics was best described by a four-compartment model featuring EHC through the inclusion of an additional gallbladder compartment, which was assumed to release drug at specific time-intervals consistent with food intake. The final PK model was able to adequately estimate the population pharmacokinetic parameters and to allow for a formal characterization of the pharmacokinetic profile and the secondary peaks due to enterohepatic recirculation.