Abstract
This paper proposes an innovative method, named b-ntPET, for solving a competition model in PET. The model is built upon the state-of-the-art method called lp-ntPET. It consists in identifying the parameters of the PET kinetic model relative to a reference region that rule the steady state exchanges, together with the identification of four additional parameters defining a displacement curve caused by an endogenous neurotransmitter discharge, or by a competing injected drug targeting the same receptors as the PET tracer. The resolution process of lp-ntPET is however suboptimal due to the use of discretized basis functions, and is very sensitive to noise, limiting its sensitivity and accuracy. Contrary to the original method, our proposed resolution approach first estimates the probability distribution of the unknown parameters using Markov-Chain Monte-Carlo sampling, distributions from which the estimates are then inferred. In addition, and for increased robustness, the noise level is jointly estimated with the parameters of the model. Finally, the resolution is formulated in a Bayesian framework, allowing the introduction of prior knowledge on the parameters to guide the estimation process toward realistic solutions. The performance of our method was first assessed and compared head-to-head with the reference method lp-ntPET using well-controlled realistic simulated data. The results showed that the b-ntPET method is substantially more robust to noise and much more sensitive and accurate than lp-ntPET. We then applied the model to experimental animal data acquired in pharmacological challenge studies and human data with endogenous releases induced by transcranial direct current stimulation. In the drug challenge experiment on cats using [18F]MPPF, a serotoninergic 1A antagonist radioligand, b-ntPET measured a dose response associated with the amount of the challenged injected concurrent 5-HT1A agonist, where lp-ntPET failed. In human [11C]raclopride experiment, contrary to lp-ntPET, b-ntPET successfully detected significant endogenous dopamine releases induced by the stimulation. In conclusion, our results showed that the proposed method b-ntPET has similar performance to lp-ntPET for detecting displacements, but with higher resistance to noise and better robustness to various experimental contexts. These improvements lead to the possibility of detecting and characterizing dynamic drug occupancy from a single PET scan more efficiently.
Highlights
Positron emission tomography (PET) is a functional 3D in vivo imaging technique that allows to visualize and quantify with a very high sensitivity the local concentration of an injected radiotracer molecule
We introduce a novel resolution method for ntPET models, called b-ntPET, whereby the parameter estimation relies on a Markov-Chain Monte-Carlo (MCMC) sampling in a Bayesian framework
We introduced a novel resolution method for ntPET models, called b-ntPET, that addresses some shortcomings and limitations of the current lp-ntPET model, and whereby the parameter estimation relies on a Markov-Chain Monte-Carlo (MCMC) sampling in a Bayesian framework, allowing the joint estimation of the noise level and the model parameters as well as the use of prior knowledge to guide the estimation process toward realistic solutions
Summary
Positron emission tomography (PET) is a functional 3D in vivo imaging technique that allows to visualize and quantify with a very high sensitivity the local concentration of an injected radiotracer molecule. In 1995, both Fisher et al (1995) and Morris et al (1995) advanced the possibility to use PET to detect dynamic changes in receptor binding or receptor occupancy occurring during activation studies. Their theory relied on the hypothesis that a cognitive task increases the firing rate of the involved neurons, leading to a release of endogenous neurotransmitter in the synaptic level with a measurable effect on the PET kinetics. One of the challenges for PET neuroimaging experiments became the design of robust and reliable kinetic analysis approaches with an integrated competition model to account for transient changes in kinetic binding and receptor occupancy in both low and high target density regions
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.