基于核医学影像的药代动力学分析

王洪凯wang.hongkai@

基于核医学影像的药代动力学分析生物医学工程原理I1整理ppt自我介绍99-03本科,北京航空航天大学 电子工程与信息技术03-09博士,清华大学

生物医学工程，医学影像处理09-14研究员,美国加州大学洛杉矶分校 分子影像，医学成像系统研究14-今副教授，大连理工大学医学影像分析，核医学影像2整理ppt报告内容核医学影像药代动力学简介房室模型原理介绍其他药代动力学模型经典案例分析3整理ppt报告内容核医学影像药代动力学简介房室模型原理介绍其他药代动力学模型经典案例分析4整理ppt1.核医学影像药代动力学简介医学影像模式(MedicalImageModality)X-RayCTMRB-ModeUltrasoundDopplerUltrasoundPET5整理ppt1.核医学影像药代动力学简介药代动力学（Pharmacokinetics）是研究药物在动物体内的含量随时间变化规律的科学，是药理学的一种。研究药物在机体内的吸收、分布、代谢及排泄的过程。应用在药物治疗、临床药理、分子药理、生物化学、生物药剂、分析化学、药剂、药理及毒理等多种科学领域中。药物动力学对指导新药设计，优化给药方案，改进剂型，提供高效、速效（或缓释）、低毒（或低副作用）的药物制剂，已经发挥了重大作用。6整理ppt1.核医学影像药代动力学简介传统药代动力学的测量方法血液：可以反映药物即刻的变化与药物的疗效、不良反应紧密相关。尿液：反映药物经过肾脏排泄的量粪便：反映药物经过肠道排泄的量其它：组织分布的浓度和量7整理ppt1.核医学影像药代动力学简介定量化的核医学影像为我们提供了研究药代动力学的无创方法可在图像中测量每个器官中的药物浓度随时间变化曲线（timeactivitycurve,TAC）15s25s45s1.5min6min15min30min60min[18F]FDG在体内分布的动态PET序列影像YunfengC.et.al,IJIG,2011.8整理ppt1.核医学影像药代动力学简介基于核医学影像的药代动力学分析流程注射示踪药物动态PET成像影像数据分析感兴趣区域划分血液时间活动曲线(BTAC)&组织时间活度曲线(TTAC)房室模型药代动力学参数曲线拟合参数计算9整理ppt1.核医学影像药代动力学简介核医学影像药代动力学应用神经影像分析心脏影像分析新药试制小动物影像分析肿瘤的评估与治疗分子影像、生物信息学10整理ppt报告内容核医学影像药代动力学简介房室模型原理介绍其他药代动力学模型经典案例分析11整理ppt2.房室模型原理介绍什么是房室模型注：dC(t)/dt表示浓度随时间的变化速度K1：药物由血浆流入组织的速率系数k2：药物由组织流入血浆的速率系数血浆中药物浓度组织中药物浓度这些比率代表了药物与组织结合和分离的能力，是我们要求解的动力学参数。通过测得C0(t)和C1(t)来求解K1和k2，这是典型的求解微分方程组12整理ppt2.房室模型原理介绍基础的四房室模型（1个血浆房室+3个组织房室）

3-tissuecompartmentmodel(3TCM)Tissue=Freetracer+Specificallyboundedtracer+non-specificallyboundedtracerFree:游离在组织中的没有跟受体结合的示踪剂No-specific:虽然与组织结合了，但是没有跟目标受体结合的示踪剂Bound:与目标受体结合的示踪剂13整理ppt2.房室模型原理介绍简化的三房室模型（1个血浆房室+2个组织房室）2-tissuecompartmentmodel(2TCM)

Free与non-specificbinding被合并成一个房室合并的前提是free与non-specificbinding之间的交换速率足够快，二者如同一体房室C1也被叫做non-displaceablecompartment（没有被置换药物的房室）

14整理ppt2.房室模型原理介绍简化的二房室模型（1个血浆房室+1个组织房室）1-tissuecompartmentmodel(1TCM)

Displaceablebinding与non-displaceablebinding被合并成一个房室合并的前提是二者之间的交换速率足够快，如同一体

15整理ppt2.房室模型原理介绍房室模型中各参数的生理含义

K1k2k3,k4,k5,k6VT(Distributionvolume,DV)

VND(Non-DisplaceableUptakeVolume)BPND(BindingPotential)

药物从血浆流入组织的能力药物从组织流回血浆的能力药物在不同房室之间流入流出的能力药物流动达到平衡后，组织与血浆中药物浓度的比值对于1TCM，VT=K1/k2；对于2TCM，VT=K1/k2(1+k3/k4)。

针对2TCM而言VND=K1/k2，表示药物达到平衡后组织中未与受体结合的药物与血浆中药物浓度比值

BPND=k3/k4，表示药物达到平衡后,组织中与受体结合的药物与未与受体结合的药物的浓度比值

16整理ppt2.房室模型原理介绍如何利用动态核医学影像来求解房室模型

要求解K1,k2,k3,k4,需要先测得CP(t),C1(t),C2(t)，即每个房室的时间活度曲线CP(t)可以通过采动脉血的方式（临床上比较难实现，常用于科研）C1(t),C2(t)咋测？C1(t),C2(t)用动态PET影像测。以2TCM为例17整理ppt2.房室模型原理介绍如何利用动态核医学影像来求解房室模型

血液的分离与测量先在分离机中分离出血浆通过过滤或HPLC方法进一步得出血浆中未被代谢的药物比例从血浆中未被代谢的药物比例中，进一步测量没有和血浆蛋白结合的药物比例，通常比较难测，默认假设这个比例是1血浆中没有被代谢掉且没有与血浆蛋白结合的药物浓度才是CP(t)18整理ppt2.房室模型原理介绍如何利用动态核医学影像来求解房室模型

胰腺用户在动态PET图像中勾勒各个器官的感兴趣区域（ROI，对于三维区域也叫VOI）软件会计算每个ROI中像素值随时间变化的曲线，转化为器官中药物浓度随时间变化的曲线CT(t)CT(t)=C1(t)+C2(t),因此在图像中没法分别测出C1(t)和C2(t)，只能测出二者之和有了二者之和，再通过数学的方法，求解K1,k2,k3,k4基于动态图像组织时间活度曲线的提取肝脏脾脏19整理ppt2.房室模型原理介绍如何利用动态核医学影像来求解房室模型

具体计算过程线性时不变系统初始条件：CE(0)=0，CM(0)=0拉普拉斯变换法求解方程：卷积运算20整理ppt报告内容核医学影像药代动力学简介房室模型原理介绍其他药代动力学模型经典案例分析21整理ppt3.其它药代动力学模型图模型（GraphicalModel，GraphicalPlot）像素级药代动力学参数图(Pixel-wisekineticmodeling)22整理ppt3.其它药代动力学模型图模型（GraphicalModel，GraphicalPlot）

通过直线拟合的方式来求解动力学参数，计算过程中会画出一幅直线图(graph)，因此得名Graphicalmodel,也叫Graphicalplot。Graphicalplot可分为Loganplot和Patlakplot两种模型。如果药物与受体是不可逆的结合(irreversiblebinding)，则用Patlakplot。如果药物与受体是可逆的结合(reversiblebinding)，则用Loganplot。23整理ppt3.其它药代动力学模型图模型之

PatlakPlot

Patlakplot用于irreversiblebinding因为是不可逆的结合，因此组织中的药物浓度与血浆中药物浓度的累积量（即对时间的积分）成正比测得CP(t)和CT(t)后，K和V可通过直线拟合算法求出K反映药物被组织代谢的速率(metabolicflux)，V反映了分布容积(distributionvolume)24整理ppt3.其它药代动力学模型图模型之

LoganPlot

Loganplot用于reversiblebinding因为是可逆的结合，因此组织中的药物浓度的时间积分与血浆中药物浓度的时间积分成正比测得CP(t)和CT(t)后，K和b可通过直线拟合算法求出K近似等于分布容积VT25整理ppt3.其它药代动力学模型像素级药代动力学模型之前的动力学参数都是以组织为单位的能否针对每一个像素计算动力学参数，从而形成一幅“参数图”？答案yes有助于在更高的空间分辨率上了解动力学参数分布每个像素就是一个ROI便于和标准参数模板进行对比PMOD软件得到的全脑像素级Vt参数图26整理ppt4.PMOD软件的药代动力学分析功能与特点像素级药代动力学模型计算每个像素的多种药代动力学参数比较两幅动力学参数图像27整理ppt报告内容核医学影像药代动力学简介房室模型原理介绍其他药代动力学模型经典案例分析28整理ppt4.经典案例分析脑部影像分析对患有多系统萎缩的老人和健康的老人分别进行123I-SPECT显像比较仿射几何法和弹性模型法这两种方法在核磁共振的脑部配准的应用同时研究了两种方法对胆碱能神经递质分布定量分析的影响。图像分析过程中进行了图像配准、图谱配准，各脑区活度曲线的测量，计算RefenceTissueModel,使用PXMOD模块，求得每个像素的BPND参数29整理ppt4.经典案例分析脑部影像分析脑部MRI图像的自动三维分割脑部标准图谱的配准MR中ROI到PET图像中的映射每个脑区的活度测量30整理ppt4.经典案例分析阿尔兹海默病影像分析比较了三种动态PET数据分析技术，即PMOD的PALZ模块、AD相关低代谢收敛指数（HCI）和meta分析感兴趣区平均代谢法（metaROI）在鉴别中度或轻度阿尔茨海默病患者和轻度认知障碍患者应用中的能力。三种分析技术的敏感度和准确性的鉴别能力相当PMOD阿尔茨海默病分析工具没有附加需求，自动化程度高，支持更多的图像数据格式。PMOD不仅可以实现全自动的PALZ指数计算，而且还有批处理模式，一次性自动处理大量的数据。31整理ppt4.经典案例分析阿尔兹海默病影像分析病人图像与正常人的自动配准与对比给出每个脑区的异常分数32整理ppt4.经典案例分析药物机理研究Nature子刊文章关于氯胺酮的抗抑郁症作用机制研究计算了的多种药代动力学参数，包括BindingPotential,LoganPlot,k2’等生成了像素级的Bindingpotential参数图33整理ppt4.经典案例分析生物分布和放射剂量研究关于新型转位蛋白显像剂F18-PBR06辐射剂量吸收的研究，以及评估F18