We explore the impact of saturable distribution over the central and the peripheral compartment in pharmacokinetic models, whilst assuming that back flow into the central compartiment is linear. Using simulations and analytical methods we demonstrate characteristic tell-tale differences in plasma concentration profiles of saturable versus linear distribution models, which can serve as a guide to their practical applicability. For two extreme cases, relating to (i) the size of the peripheral compartment with respect to the central compartment and (ii) the magnitude of the back flow as related to direct elimination from the central compartment, we derive explicit approximations which make it possible to give quantitative estimates of parameters. In three appendices we give detailed explanations of how these estimates are derived. They demonstrate how singular perturbation methods can be successfully employed to gain insight in the dynamics of multi-compartment pharmacokinetic models. These appendices are also intended to serve as an introductory tutorial to these ideas.