Measurement of BuChE activity with
[11C]MP4B PET
Model
1-[11C]methyl-4-piperidinyl n-butyrate ([11C]MP4B or [11C]nBMP or [11C]BMP) was suggested as an optimal tracer for measuring butyrylcholinesterase (BuChE) activity in the brain (Snyder et al., 2001) and was further validated by Roivainen et al. (2004).
For the quantification of BuChE activity, a three-compartment (two-tissue compartment) model is applied. The rate constant k3 represents the rate of hydrolysis of [11C]MP4B by BuChE. The diffusion of radioactive metabolite of [11C]MP4B through the blood-brain-barrier is assumed negligible during the PET study, thus k4=0.
Kuhl et al. (2006) estimated K1, k3 and VB for [11C]MP4B pixel-by-pixel, constraining K1/k2 to a value determined as the mean estimate across cortical regions of all subjects from the unconstrained fits (Koeppe et al. 1999; Kuhl et al., 1999).
Analysis method in TPC
Pre-processing plasma input
Make sure that you have all the necessary data files:
- On-line sampler data file
- Dynamic PET image file from this study
This file is used to correct for possible start time mismatch between PET and blood sampling. Normally, both were started simultaneously, and in that case you do not need the dynamic image file yet. - Count-rate curve
- Plasma curve from manual sampling
- Plasma metabolite data file
Then, follow the instructions in http://www.turkupetcentre.net/analysis/doc/tracer_input.html. Note that for [C-11]MP4B:
- Conversion from blood to plasma or from plasma to blood is not necessary, because blood and plasma TACs are similar
- Hill-type function should be used to fit the fractions of parent tracer in plasma
- In the time delay correction, do not use the whole data that was collected but only the first 10 minutes after injection. Time delay must be fitted between PET count-rate curve and metabolite corrected plasma TAC, but remember to correct simulatenously also the total plasma TAC.
Image processing
PET images are summed over frames, coregistered with MRI, and regions of interest are defined.
For regional analysis, TACs are calculated from the dynamic PET image.
Before proceeding, make sure that both the plasma and tissue data are in the same calibration units (preferably kBq/ml) and that the time unit is min. Image data from HR+ and PET-CT may originally be in units Bq/ml.
Regional BuChE activity (k3)
After all the previous steps have been done succesfully, the enzyme activity k3 can be calculated using programs lhsol or fitk3. Program fitk3 allows constraining K1/k2 to a predetermined value and fitting the vascular volume fraction. If K1/k2 is not constrained, you should consider reporting (K1/k2)*k3 as an index of enzyme activity, instead of k3.
Examples:
lhsol -k3 ua2826ap_comb_pure.delay.kbq ua2826dy1.dft 999 ua2826k3.res ua2826k3fit.dat
fitk3 -i=constraints.set ua2826ap_comb_pure.delay.kbq ua2826ab_comb.delay.kbq ua2826dy1.dft 999 ua2826k3.res ua2826k3fit.dat
Note that you must use the delay fitted (filename usually contains *.delay.*) and metabolite corrected plasma curves.
If necessary, there are tools for processing the result file further.
Butyrylcholinesterase activity maps
To be added later.