Compartmental model ODEs
Compartmental models for PET data can be described in terms of a set of linear, first-order, constant-coefficient, ordinary differential equations (ODEs). Change of tracer concentration in one of the tissue compartments is a linear function of the concentrations in all other compartments:
From the linearity follows that kinetic measurements are the convolution of the tracer input function and the response function of the system.
For example, for two-tissue compartment model with tissue compartments C1 and C2, input function C0(t), and rate constants K1, k2, k3, and k4, can be described mathematically with a system of two ODEs:
The simulated tissue TACs are used to
- test and validate analysis methods and software
- estimate compartmental model rate constants (model fitting).
Analytic solutions using Laplace transformation
Analytic solutions for linear first-order differential equations can be derived using Laplace transformation. For this 2TCM example the solution is (Phelps et al., 1979):
, where ⊗ denotes the operation of convolution. For one-tissue compartmental model the solution is a lot simpler:
For simplified reference tissue model the solution is:
Alternative solutions for discrete-time data
Multistep predictor-corrector algorithms for ODE solution avoids the computationally slow convolution operation. Adams-Bashforth Predictor requires present and past values of x (xm and xm-1) to predict xm+1:
Δt is the sample time difference. Fourth-order Adams Predictor requires four previous values to predict xm+1:
As a fourth-order algorithm its errors are similar to the fourth-order Runge-Kutta algorithm (DiStefano, 2013). The implicit Second-order Adams-Moulton Corrector requires only one value of the past:
The Second-order Adams-Moulton method has been applied to PET data modelling (Kuwabara et al., 1993). Solutions for the most common compartmental models using this method are shown in TPC modelling reports.
As an example, the solution of the one-tissue compartmental model ODE,
starts by integrating the equation, assuming that at time zero all concentrations are zero:
Integrals are calculated using trapezoidal method. Definite integral of n th compartment is implicitly estimated based on second-order Adams-Moulton method with Euler integration:
In this example we have just one compartment (n=1). After substitution and rearrangement:
For simplified reference tissue model the solution is:
Using this method the TACs from compartmental models can be easily simulated even in spreadsheet software. For an example, you can download an Excel file where PET tissue data can be simulated using two-tissue compartmental model from measured plasma data and user definable rate constants. Excel files for simulating tissue data with three-tissue compartmental models are also available (compartments in parallel and in series).
Comparison of solution methods
The convolution method provides precise results only if the sampling time interval is very short (relative to the rate constants). The alternative method is computationally faster, and produces correct results (as far as it is possible from discrete data) without any interpolation; the convolution method only approaches the correct result with shorter interpolated sampling intervals (that is, smaller Δt is required in the convolution method). As an example, convolution and the alternative Adams-Moulton method are compared in a simulation of one-tissue compartmental model.
With data sampled at discrete times some errors in simulation will always be present. The longer the sampling time frame durations are, the more important it is to appropriately account for the frame durations in calculating AUC (integral).
- Fitting compartmental models
- Compartmental model analysis of regional TACs
- Reference tissue input compartmental models
- PET data
- Input function
- Model equations for compartmental model with dual plasma input: models in parallel and in series
Budinger TF, Huesman RH, Knittel B, Friedland RP, Derenzo SE (1985): Physiological modeling of dynamic measurements of metabolism using positron emission tomography. In: The Metabolism of the Human Brain Studied with Positron Emission Tomography. (Eds: Greitz T et al.) Raven Press, New York, 165-183.
DiStefano III J. Dynamic Systems Biology Modeling and Simulation. Academic Press, 2013. ISBN: 9780124104112.
Gunn RN, Gunn RS, Cunningham VJ. Positron emission tomography compartmental models. J Cereb Blood Flow Metab. 2001; 21: 635-652. doi: 00004647-200106000-00002.
Huang SC, Phelps ME (1986): Principles of tracer kinetic modeling in positron emission tomography and autoradiography. In: Positron Emission Tomography and Autoradiography: Principles and Applications for the Brain and Heart. (Eds: Phelps,M; Mazziotta,J; Schelbert,H) Raven Press, New York, 287-346.
Kuwabara H, Cumming P, Reith J, Léger G, Diksic M, Evans AC, Gjedde A. Human striatal L-DOPA decarboxylase activity estimated in vivo using 6-[18F]fluoro-DOPA and positron emission tomography: error analysis and application to normal subjects. J Cereb Blood Flow Metab. 1993; 13: 43-56. doi: 10.1038/jcbfm.1993.7.
Phelps ME, Huang S-C, Hoffman EJ, Selin C, Sokoloff L, Kuhl DE. Tomographic measurement of local cerebral glucose metabolic rate in humans with [F-18]2-fluoro-2-deoxy-D-glucose: validation of method. Ann Neurol. 1979; 6: 371-388. doi: 10.1002/ana.410060502.
Tags: Modeling, Compartmental model, Simulation, ODE, Convolution
Updated at: 2020-01-20
Created at: 2004-06-03
Written by: Vesa Oikonen