Fractional kinetics in multi-compartmental systems

Fractional calculus, the branch of calculus dealing with derivatives of non-integer order (e.g., the half-derivative) allows the formulation of fractional differential equations (FDEs), which have recently been applied to pharmacokinetics (PK) for one-compartment models. In this work we extend that theory to multi-compartmental models. Unlike systems defined by a single ordinary differential equation (ODE), considering fractional multi-compartmental models is not as simple as changing the order of the ordinary derivatives of the left-hand side of the ODEs to fractional orders. The latter may produce inconsistent systems which violate mass balance. We present a rationale for fractionalization of ODEs, which produces consistent systems and allows processes of different fractional orders in the same system. We also apply a method of solving such systems based on a numerical inverse Laplace transform algorithm, which we demonstrate that is consistent with analytical solutions when these are available. As examples of our approach, we consider two cases of a basic two-compartment PK model with a single IV dose and multiple oral dosing, where the transfer from the peripheral to the central compartment is of fractional order alpha < 1, accounting for anomalous kinetics and deep tissue trapping, while all other processes are of the usual order 1. Simulations with the studied systems are performed using the numerical inverse Laplace transform method. It is shown that the presence of a transfer rate of fractional order produces a non-exponential terminal phase, while multiple dose and constant infusion systems never reach steady state and drug accumulation carries on indefinitely. The IV fractional system is also fitted to PK data and parameter values are estimated. In conclusion, our approach allows the formulation of systems of FDEs, mixing different fractional orders, in a consistent manner and also provides a method for the numerical solution of these systems.