Logan plot

A Logan plot (or Logan graphical analysis) is a graphical analysis technique based on the compartment model that uses linear regression to analyze pharmacokinetics of tracers involving reversible uptake. It is mainly used for the evaluation of nuclear medicine imaging data after the injection of a labeled ligand that binds reversibly to specific receptor or enzyme.

In conventional compartmental analysis, an iterative method is used to fit the individual model parameters in the solution of a compartmental model of specific configuration to the measurements with a measured plasma time-activity curve that serves as an forcing (input) function, and the binding of the tracer can then be described. Graphical analysis is a simplified method that transforms the model equations into a linear equation evaluated at multiple time points and provides fewer parameters (i.e., slope and intercept). Although the slope and the intercept can be interpreted in terms of a combination of model parameters if a compartmental model configuration is assumed, the graphical methods are independent of any specific model configuration. In case of irreversible tracers, certain fraction of the radioactivity is trapped in the tissue or the binding site during the course of the experiment, whereas reversible tracers show uptake and loss from all compartments throughout the study. The theoretical foundation of graphical analysis for irreversible tracers (also called Patlak graphical analysis or Patlak plot) was laid by Clifford Patlak and his colleagues at NIH. Based on the original work of Patlak, Jean Logan and her colleagues from Brookhaven National Laboratory extended the method to tracers with reversible kinetics.

Description
The kinetics of radiolabeled compounds in a compartmental system can be described in terms of a set of first-order, constant-coefficient, ordinary differential equations. The time course of the activity in the multicompartmental system driven by a metabolite-corrected plasma input function $$C_p(t)$$ can be described by:


 * $$\frac{d\mathbf{A}}{dt} = \mathbf{KA} + \mathbf{Q}C_p(t)$$

where $$\mathbf{A}$$ is a column vector of activity concentration for each compartment at time $$t$$, $$\mathbf{K}$$ is the matrix of the transfer constants between compartments, and $$\mathbf{Q}$$ is the vector of plasma-to-tissue transfer constants. Patlak and Blasberg showed that the above equation can be written as:


 * $$\int_0^t A(\tau) \, d\tau = -\mathbf{U}_n^T \mathbf{K}^{-1} \mathbf{Q} \int_0^t C_p(\tau) \, d\tau + \mathbf{U}_n^T \mathbf{K}^{-1} \mathbf{A}$$

where $$\mathbf{U}_n^T$$ represents a row vector of 1s and $$A(t) = \mathbf{U}_n^T \mathbf{A}$$. The total activity in the region of interest, $$\mathrm{ROI}(t)$$, is a combination of radioactivities from all compartments plus a plasma volume fraction ($$V_p$$) and thus:


 * $$\int_0^t \mathrm{ROI}(\tau) \, d\tau = (-\mathbf{U}_n^T \mathbf{K}^{-1} \mathbf{Q} +V_p) \int_0^t C_p(\tau) \, d\tau + \mathbf{U}_n^T \mathbf{K}^{-1} \mathbf{A}$$

By dividing both sides by $$\mathrm{ROI}(t)$$, one obtains the following linear equation:


 * $${{\int_0^t \mathrm{ROI}(\tau) \, d\tau} \over \mathrm{ROI}(t)} = (-\mathbf{U}_n^T \mathbf{K}^{-1} \mathbf{Q} +V_p) {{\int_0^t C_p(\tau) \, d\tau} \over \mathrm{ROI}(t)} + {{\mathbf{U}_n^T \mathbf{K}^{-1} \mathbf{A}} \over {\mathbf{U}_n^T \mathbf{A} + V_p C_p}}$$

For $$t > t'$$, Patlak and his colleagues showed that $$\mathbf{A} = -\mathbf{K}^{-1} \mathbf{Q} C_p(t)$$, i.e., the steady-state condition. When this condition is satisfied, the intercept has reached its constant value so that after some time a plot of $${{\int_{0}^{t} C_p(\tau)d\tau} \over \mathrm{ROI}(t)}$$ versus $${{\int_{0}^{t} \mathrm{ROI}(\tau)d\tau} \over \mathrm{ROI}(t)}$$ becomes a straight line with slope $$(-\mathbf{U}_n^T \mathbf{K}^{-1} \mathbf{Q} + V_p)$$ and intercept $${{\mathbf{U}_n^T \mathbf{K}^{-1} \mathbf{A}} \over {\mathbf{U}_n^T \mathbf{A} + V_p C_p}}$$.

For a catenary two-tissue compartment model with transfer constants $$K_1$$ (forward transport from plasma to tissue), $$k_2$$ (reverse transport from tissue to plasma), $$k_3$$ (binding parameter proportional to $$B_\max k_\mathrm{on}$$), and $$k_4$$ (dissociation constant) to analyze enzyme or receptor system, the slope represents the total distribution volume ($$V_d$$) and is given by $$\frac{K_1}{k_2}(1+\frac{k_3}{k_4}) + V_p$$, where $$k_3 = B_\max k_\mathrm{on}$$, $$k_4 = k_\mathrm{off}$$, $$\frac{k_3}{k_4} = \frac{B_\max}{K_d}$$, and $$K_d = k_\mathrm{off}/k_\mathrm{on}$$, in which $$B_\max$$ is the concentration of ligand binding sites, $$K_d$$ is the equilibrium dissociation constant for the ligand-binding site complex, $$k_\mathrm{on}$$ is the ligand-binding association constant, $$k_\mathrm{off}$$ is the ligand-binding dissociation constant. For a one-tissue compartment model with transfer constants $$K_1$$ and $$k_2$$, the slope is $$\lambda + V_p$$, where $$\lambda$$ is the partition coefficient ($${K_1}/{k_2}$$) and the intercept is $$\frac{-1}{k_2(1+V_p/\lambda)}$$.