Spatial effects in modeling pharmacokinetics of rapid action drugs.

We develop a model to describe the time course of plasma concentration of neuromuscular blocking agents used as anesthetics during surgery. This model, which overcomes the limitations of the classical compartment models routinely used in pharmacokinetics, incorporates spatial effects due to heterogeneity in the circulation: it takes the form of a dispersion equation on a circular domain, with a time-dependent leakage term. This term is fitted to the functional form desired once first-stage transients have died out. Comparisons are made with clinical data by adjusting three parameters in the model.