## p2t_v3c - tpcclib 0.7.8 © 2022 by Turku PET Centre

```Deprecated, use sim_3tcm instead.

Simulation of PET tissue time-radioactivity concentration curve (TAC)
from arterial plasma (Ca) and blood (Cb) TACs, based on three-tissue
compartmental model, where the compartments are in series (default):

____    K1   ____    k3   ____    k5   ____
| Ca | ----> | C1 | ----> | C2 | ----> | C3 |
|____| <---- |____| <---- |____| <---- |____|
k2           k4           k6

dC1(t)/dt = K1*Ca(T) - (k2+k3)*C1(T) + k4*C2(T)
dC2(t)/dt = k3*C1(T) - (k4+k5)*C2(T) + k6*C3(T)
dC3(t)/dt = k5*C2(T) - k6*C3(T)
Ct(T) = C1(T) + C2(T) + C3(T)
Cvb(T) = Cab(T) - dCt(t)/dt / f
Cpet(T)= Vb*fA*Cab(T) + Vb*(1-fA)*Cvb(T) + (1-Vb)*Ct(T)

, or, optionally, three-tissue compartmental model, where the 2nd and 3rd
tissue compartments are parallel, often used to represent specific and
non-specific binding:

____    K1   ____    k3   ____
| Ca | ----> | C1 | ----> | C2 |
|____| <---- |____| <---- |____|
k2           k4
| ^
k5 | | k6
v |
____
| C3 |
|____|

dC1(t)/dt = K1*Ca(T) - (k2+k3+k5)*C1(T) + k4*C2(T) + k6*C3(T)
dC2(t)/dt = k3*C1(T) - k4*C2(T)
dC3(t)/dt = k5*C1(T) - k6*C3(T)
Ct(T) = C1(T) + C2(T) + C3(T)
Cvb(T) = Cab(T) - dCt(t)/dt / f
Cpet(T)= Vb*fA*Cab(T) + Vb*(1-fA)*Cvb(T) + (1-Vb)*Ct(T)

Usage: p2t_v3c [options] plasmafile bloodfile K1 k2 k3 k4 k5 k6 f Vb fA simfile

Options:
-paral[lel]
Model with parallel compartments C2 and C3 is applied.
-ser[ies]
Model with compartments C1, C2, and C3 in series is applied (default).
-sub | -nosub
TACs of sub-compartments (C1, C2 and C3) are written (-sub, default)
or not written (-nosub) into the simfile.
-h, --help
Display usage information on standard output and exit.
-v, --version
Display version and compile information on standard output and exit.
-d[n], --debug[=n], --verbose[=n]
Set the level (n) of debugging messages and listings.
-q, --quiet
Suppress displaying normal results on standard output.
-s, --silent
Suppress displaying anything except errors.

Enter blood flow (f) in units (ml/(min*ml)); set f=0 to assume that f>>K1
and Cvb=Cab. Enter the vascular volume fraction (Vb) and its arterial
fraction (fA) as percentages.
If the times in input files are in seconds, the units of rate constants
(k's) and blood flow (f) must also be specified as 1/sec.
For accurate results, plasma TAC should have very short sampling intervals.
To reduce the model, k5 or k3 can be set to 0.

Example:
p2t_v3c -nosub plasma.dat blood.dat 0.3 0.2 0.1 0.2 0 0 0.6 5 30 sim.dat

Simulated TACs are written in ASCII format with columns:
1) Sample time
2) Total tissue activity concentration (Cpet)
3) Activity concentration in 1st tissue compartment, (1-Vb)*C1
4) Activity concentration in 2nd tissue compartment, (1-Vb)*C2
5) Activity concentration in 3rd tissue compartment, (1-Vb)*C3
6) Arterial contribution to tissue activity, Vb*fA*Cab
7) Venous contribution to tissue activity, Vb*(1-fA)*Cvb

References:
1. TPCMOD0001.