房室模型的确定及参数计算.ppt

,第4章房室模型的确定及参数计算,第1节房室模型的确定,药代动力学系通过“速率类型”和“数学模型与隔室”这两大要素来分析药物体内动态规律的，这里十分重要的问题是要建立一个合适的房室模型，亦既房室数的确定问题，同样一组血药浓度资料，房室模型确定不当，将导致错误的结果。,一、最佳房室数确定原则,1.希望测定值能够均匀而随机地分布在拟合曲线的两侧。2.适当地使残差平方和（s）或加权残差平方和（sw）达到最小。,线性数学模型的血药浓度-时间曲线关系式一般为：,(4-1),i)2=,(4-2),(4-3),n为房室数，为房室序数，m为采血时间次数，xj、j为待定参数，ct为时刻血药浓度，ci为第i次取样时的血药浓度实测值，为第i次取样时的血药浓度的理论计算值，s为残差平方和，sw为加权残差平方和，wi为权重系数,c-t散点图判断法iv后血药浓度（c）对时间（t）在半对数坐标纸上绘出散点图，由散点图形估计房室数。如各数据点可用一条直线拟合，可初估为一室模型，拟合单指数方程（方程式4-4）:c=c0e-kt(4-4),二、确定房室数的具体方法,如图形在一处或两处出现转折，血药浓度呈现先快后慢的衰减曲线，可初估为二室或三室模型，拟合双指数或三指数方程（方程式4-5或4-6）,c=ae-t+be-t(4-5)c=ae-t+be-t+ge-t(4-6),为血药浓度，t为时间，、在二室模型中分别为分布速率常数和消除速率常数，在三室模型中分别为快分布相和慢分布相速率常数。为三室模型消除相速率常数。、b、g为、和相延伸线在纵轴的截距。e为自然对数的底(1+1/n)n的极限,血管外给药后的房室模型可根据c-t曲线中吸收后相的曲线形状加以估计，如吸收后相曲线为一直线，则可估计属于口服一室模型，拟合双指数方程（方程式4-7）。,如吸收后相的曲线形状表现为先快后慢的衰减曲线，则可估计属于口服二室模型，拟合三指数方程式（方程式4-8）。,（4-8）,（4-7）,应注意的是，一个静脉注射为二室模型的药物，如果分布速率并不十分快，其口服吸收曲线可能表现为口服一室模型，这是因为缓慢的分布相为吸收相所掩盖；如该药分布十分迅速，则可表现为口服二室模型。散点图法简单，但比较粗糙，不够准确，需采用以下方法进一步确证。,2.残差平方和或加权残差平方和判断法将c-t数据分别按一室、二室或三室模型拟合，求出相应的c-t方程式），然后按此方程式计算出不同时间的理论血药浓度，称之为计算值，实测值与计算值之差称为残差，求出s或sw和，s越小说明计算值与实验值契合程度就越高，因此，拟合的房室模型中，s或sw最小者即为所求的房室模型。,例如，以此法算得口服锂盐后按一室和二室拟合的残差平方和分别为0.75和0.00513，后者比前者小146倍，说明患者对锂盐的药代动力学具有二室模型特征。,在sw的计算中需要对数据进行权重，权重的含义在于：线性动力学中，药代动力学参数的计算常将血药浓度转换为对数浓度后对时间进行直线拟合(logct)，在计算中常用最小二乘法，但经转换后的数据残差平方和达最小，不等于原来数据的残差平方和达最小。在测得的一组血药浓度数值中，高低浓度相差较大（如相差两个数量级）时，若实验值与计算值残差不经权重，则在最小二乘法计算过程中，低浓度数据的作用将会被忽视，而实际上低浓度实验数据对曲线的拟合都是十分重要的。由于以上原因，在应用最小二乘法拟合药-时曲线时，需先对数据进行权重。,权重系数的确定有不同的方法,可采用:1/c1/c2多数文献采用1/c2（浓度倒数的平方）进行权重。,3.拟合优度r12值判断法,根据实验值与计算值按下式求得r12，在所拟合的房室模型中，r12值大的为最佳房室模型。,(4-9),4.f测验判断法,f值按下式计算：sw1及sw2分别为第一种和第二种模型的加权残差平方和df为自由度，即各自的实验数据点的数目减去参数的数目，iv一室模型的c-t方程式需确2个参数(c=c0e-kt)iv二室:4个参数(c=ae-t+be-t)iv三室:6个参数(c=ae-t+be-t+ge-t)如某实验测得12个实验数据点，则上述三种模型的df分别为10、8、6。如算得的f值比相应自由度的界值（显著水平）大，便可认为将参数的数目从i增至i+1是有意义的。,(4-10),5.aic判断法（akaikesinformationcriterion）,该法首次由日本统计学家赤池弘次（akaike）提出，该氏从信息理论出发，提出一种信息标准（informationcriterion），以便对信息量作出数字上的表达，并用统计学方法确定拟合于一组实验数据的数学方程的参数数目，故称aic法。akaike及tanabe根据随机误差遵从gaussion分布的假设，以下列方程式定义aic：为实验数据点的数目为拟合的房室模型的相应数学方程式中参数数目（p=2n，n为房室数）re为加权残差平方和（与方程式4-3中sw含义相同）,aic=nre+2p（4-11）,在拟合的房室模型中，aic小为好，aic最小的数学方程式被认为是对血药浓度时程的最佳表达，这种统计学方法谓之最低aic测定（minimumaicestimation，maice）。对于maice来说，不要求进行测验及显著水平测定。如按re及aic判断结果不一致，而re相差不大时，以aic为判定标准。,6.其他,多数药物属于二室模型其速率常数具有下列特征：且相差较大，k12、k21、ke（药物自中央室消除的一级速率常数）均为正值，k12+k2120ke。如或=，可视为一室模型加以处理。属于前者的药物，向组织的分布相当迅速，分布相比消除相快得多，以致在药动学处理中，分布相可以忽略。属于后者的药物，分布相和消除相速率几乎相等，故可认为药物在机体内瞬即达平衡，因而也可作一室处理。一般来说，当k12+k2120ke，二室模型可作一室处理。,三、举例,1.石吊兰素表4-1及图4-1分别为石吊兰素给大鼠静脉注射后的血药浓度及药-时曲线，由图可知，血药浓度初期下降很快，后期下降缓慢，曲线在30分钟有一转折，故可认为该药的药-时曲线可能为符合方程式4-5的双指数衰减曲线，其药-时资料可按二室开放模型处,表4-1大鼠尾静脉注射石吊兰素15mg.kg-1后的血浆药物浓度,图4-1大鼠静脉注射石吊兰素后血浆药物浓度的衰减,为验证房室数，曾用f测验值、加权残差平方和sw、拟合度r12值，f测验值及aic等方法对于一室及二室的契合程度进行了比较。具体步骤如下：由表4-1所列数据经最小二乘法分别按一室及二室模型拟合曲线，得相应的药-时曲线方程式：c=235.170e0.026t（一室）c=210.86e0.222t+164.36e0.020t（二室）按上述c-t方程式分别标出两种模型不同时间的血药浓度，谓之血药浓度的计算值，结果见表4-3中的第3栏、第5栏数值。按方程式4-2和4-3分别求出两种模型血药浓度的测定值ci与计算值的残差平方和s（表4-3中的第4栏，第6栏数值）以及加权残差平方和sw（见表4-4f测验值计算项下）。将以上数值代入有关公式计算r12值、f值及aic值，结果见表4-2。,由以上计算结果可知，单室与二室模型比较后sw值以二室为小r12值以二室为大，f计算值大于f界值，差异显著aic值以二室为小故四种方法均表明大鼠静脉注射石吊兰素后的药-时数据属二室模型。,一室,ci2=264588.15,二室,二室r12一室r12,一室,查f值表，相应自由度的f界值（机率5%）为6.96，f计算值f界值。,2.磺胺异噁唑,图4-2为人静注磺胺异噁唑后的血药浓度-时间曲线，该药-时曲线分别按一室、二室及三室拟合，其药时方程式及所计算的re，aic以及f测验值如图4-2所示，由图可知三种模型中，二室模型的aic最小，,一室与二室比较后的f测验值f12为488，大于f界值（f28（0.01）=8.65），具有显著性，而二室与三室比较后的f测验值f23为0.741，小于f界值（f26（0.05）=5.14），因此，由aic值及f测验值均表明该药-时曲线拟合二室模型，由图中数据还可以看到，二室模型的re值远小于一室模型但较三室模型略大，按maic法，在以aic值与re值进行判断时如有矛盾，以aic法为准。由本例还可以说明aic法与f测验值结果相符，因此在用aic法时可不必进行f测验。,应说明的是：房室模型的确定主要取决于该药本身的化学结构及药动学性质，然而，由于实验目的以及实验设计（如采血时间，给药途径，检测方法的灵敏度）等的不同，同一药物也可被判断为不同的房室模型，如头孢孟多(cefamandole)，以往文献报道为二室模型，aziz等缩短了取样间隔，延长了取样时间（6小时，取样15次）以及应用高灵敏度的hplc测定方法，结果证明还存在一个十分快速的分布相（分布半衰期为5.2min，持续约20min），由于这一发现，aziz认为该药属三室模型。另如核黄素的房室数就是单室、二室及多室模型等不同报道。saunders建议，在采集血样的安排方面，峰浓度前应有3个采样点，达峰时到二倍达峰时之间应不少于2点，此后到10倍达峰时应不少于4点。在具体确定某一药物的房室数时，有时用不同方法得出不同的房室数，这时应采用上述几种方法综合判断。如石吊兰素给猴静注石吊兰素后，分别拟合一室及二室模型，f测验值不显著，但r12值以二室为大，另以两种模型的计算值与实测值直观比较，以二室模型的符合程度较佳，故以二室开放模型来描述较为恰当。,f82(0.01)=8.65,f62(0.05)=5.14,第2节主要药物代谢动力学参数的求算,药代动力学参数的计算主要是根据给药后血药浓度测定数据以及尿药测定数据进行的，前者资料可靠、结果准确，但需抽血测定，后者准确性差，因尿排药速率受排尿量的影响很大。此外，尚可应用唾液药物浓度计算药动学参数。,一、以血药浓度数据求算药物代谢动力学参数,（一）一室模型药物快速静脉注射1.消除速率常数k与半衰期t1/2由方程式c=c0k/2.303t可以看出，从实验所得的血药浓的对数对时间作图，即“c-t”图应为一条直线，该直线可用最小二乘方回归拟合，其斜率为k/2.303，由此可求得k，进而应用公式t1/2=0.693/k便能很方便的求出t1/2。如将此直线延伸，其与纵轴交点所代表的血药浓度数值为c0。（图4-3）,图4-3一室模型药物iv后的药-时曲线,举例：见表-表-某药按300mg.kg-1剂量静脉快速注射后的血药浓度及有关计算,x=91x2=2213y=27.703xy=216.573（x）2=8281，xy=9127.703=2520.973,设回归方程为y=a+bx，a为截距，b为斜率。,k=0.139h-1（因y为浓度的自然对数，故c-t作图的斜率，即为消除速率常数，不必再乘以2.303）t1/2=0.693/0.139=4.98h截距a=ybx=5.765c0=5.765c0=318.9g.ml-1,2.表观分布容积vd根据静脉注射剂量a0及零时血药浓度c0求vd。,(4-11)c0用延伸法求得，既将ct数据按最小二乘法回归拟合的直线延伸，与纵轴交点所表示的血药浓度，在表4-5所举的例子中，可通过下列计算求得：因直线的截距a=ybx=5.765即c0=5.765c0=318.9g.ml-1因静脉剂量为300mg.kg-1，所以在上题中的表观分布容积vd=300103/318.9=941ml.kg-1,3.药-时曲线下面积auc,药代动力学研究中，特别是生物利用度测定时，常需计算auc。当一室模型药物静脉推注给药，可用公式很方便地计算出auc。根据定义有下列方程：,又由于c=c0e-kt，则,积分得：,在表4-5的例子中，c0为318.9g.ml-1，k为0.139h-1，故auc为2294.2g.h-1（318.9g.ml-10.139h-1）。,(4-12),（4-13）,4.清除率,清除率可分为总体清除率（cl）和肾清除率（clr）等。（1）总体清除率（cl）可定义为相对于血药浓度的药物消除速率，由此可得方程：,在一级动力学中，,又由于a为体内药量，a与c的比值等于表观分布容积，故上式可改写为：cl=kvd（4-14）由式4-14可知清除率为消除速率常数和表观分布容积的乘积。式4-14可转变为：,因c0/k=auc（方程式4-13），故,式4-15表明总体清除率为静脉注射剂量与药-时曲线下面积的比值。,(4-15),（2）肾清除率（clr）可定义为相对于血药浓度的尿药排泄速率。故,cu为尿中药物浓度（mg.ml-1），vu为每分钟尿量（ml.min-1），cpmid为集尿时间中点的血药浓度，上式中cuvu亦即集尿间隔期内尿药排泄速率（ml.min-1）。又clr=clfefe为尿中排出的原形药的总量占给药剂量的分数，fe可根据下列公式计算,u为尿中排出的原形药的总量（实际上尿收集时间应至少在5个半衰期以上）。,在应用式4-18时应注意，集尿期应尽可能长，通常应大于5个半衰期。,clr的另一求算方法是测定集尿期内尿中累计排出的原形药的总量ae再测定集尿期间药-时曲线下面积auc，然后按下列公式求算clr：,(4-16),(4-17),(4-18),（二）二室模型药物快速静脉注射,1.速率常数及半衰期通常用残差法（methodofresiduals）进行求算，该法是求解药动学参数时最常用的一种方法，所谓残差法是将一个多项指数式解析成其中的不同成分，即当时间t足够大时，其中速率常数较大的指数项趋近于零，从而略去以致剩下速率常数较小的单项指数，经对数转换，作线性回归，然后逐项剥脱，最后留下一个单项指数式，以残差浓度对时间作回归直线，最终各项参数均可一一解出。由于它能将一个多项指数式解析成其中的不同的指数成分，故又称之为羽状法（featheringmethod）或剥去法（peelingmethod）。二室模型药物快速静注时，血药浓度与时间关系可用下式描述：c=aet+bet（4-19）,在式（4-19）中，、分别为表观一级分布速率常数和表观一级消除速率常数。通常，故当t充分大时，aet0，而bet仍保持一定数值，此时方程式（4-19）可简化为：c=be-t（4-20）上式两边取常用对数得：(4-21),由式（4-21）可知，c对t作图应得一条直线（称线），其斜率为/2.303，故可求，其相应的半衰期可进一步根据=0.693/求得，将线延伸至与纵轴相交，从其截b可求得b。,式（4-19）和（4-20）相减得残差方程：cr=aet（4-22）cr为残差浓度，即将分布相中血药浓度实测值（c）减去线上个相应点的数值。式（4-22）两边取常用对数得：,(4-23),同理，cr对t作图应得一条直线（称线），由该线的斜率（/2.303）可求得，截距a可求a，以及通过公式=0.693/进一步求出