1、用图像估算周期性血管外给药药物动力学参数【摘要】 目的：直观、快捷地估算药物动力学参数的数值。方法：对周期性血管外给药二室模型的药 时曲线与参数的几何关系做了探讨，在 MATLAB 的图形界面上，估算出了周期性血管外给药药物动力学参数的数值。结果：分别得到了坪浓度、维持剂量、周期性给药的次数、达坪分数、积累系数等参数的数值。结论：为药物动力学参数的估算提供了新的途径，该思路可拓广到同类问题的估算上。 【关键词】 周期性血管外给药; 药物动力学参数; 图像; 估算利用数学模型确定药物动力学参数，需要做繁杂的数值计算，这对于基层医药科研人员来说是一个不小的障碍。本文在分析模型参数几何意义的基础上，
2、借助 MATLAB 的作图功能，探求周期性血管外途径给药二室模型药物动力学参数的估算。1 周期性血管外给药二室模型与药物动力学参数的几何关系1.1 周期性血管外给药二室模型及其参数周期性血管外 n 次给药后，体内血药浓度 C(n)1 与时间 t 的二室模型为：C(n)1(t)=A11-e-n1-e-e-+A21-e-n1-e-e-(A1+A2) 1-e-nK1-e-Ke-K(1)模型(1)中的 t 是从第 n 次给药后开始计算的时间(0t)，其它各系数的含义见文献1。若让模型(1)中的 n 趋向无穷，得到周期性血管外给药二室模型的坪浓度：Css(t)=A111-e-e-+A211-e-e-(A
3、1+A2)11-e-Ke-K(2)由模型(1)、(2)在一个周期 内的平均值的比值得到达坪分数：fss(n)=1JF(D0C(n)1(t)dtJF)1JF(D0Css(t)dtJF)=A1(1-e-n)+A2(1-e-n)-+A2Ka(1-e-nKa)A1+A2-A1+A2Ka(3)最低坪浓度 Cmin 与第一次给药末端浓度 C1() 的比是积累系数：R=A11-e-+A21-e-A1+A21-e-K A1e-+A2e-(A1+A2)e-K(4)通常在已知周期性血管外给药二室模型的情况下，根据上述各式和牛顿迭代法可以计算出坪浓度、维持剂量、达坪状态时周期性给药的次数等药物动力学参数。1.1 模
4、型的曲线与药物动力学参数间的几何关联由于模型的函数为初等函数，故模型对应的曲线是连续、渐进变化的，不会出现间断或剧烈波动的情况，从而为从图形中估算参数的数值提供了条件。下面的几个参数与模型的曲线有如下关联： 坪浓度Css 是(1)的曲线稳定变化时一个周期内血药浓度与时间的对应关系。当给药次数 n 较大时，(1)的曲线在一个周期内所对应部分近似为坪浓度(2)的曲线。 最低坪水平 Cmin 是(1)的曲线在达坪状态一个周期内的最低点的纵坐标，即达坪状态下一个周期内的最小值。 最高坪水平 Cmax 是(1)的曲线在达坪状态一个周期内的最高点的纵坐标，即达坪状态下一个周期内的最大值。 峰时 Tm0 是
5、(1)的曲线在达坪状态下一个周期内最高点与最低点的横坐标的差(0 当周期性给药剂量 D0 一定时，改变给药周期 ，通过(1)的曲线得到周期变化后的最高或最低坪浓度。反之，若给定安全浓度、有效坪浓度，由(2)的曲线可以估算出符合条件的给药周期 。 由(3)的曲线上点的纵横坐标可以获得达坪分数 fss(n) 与周期性血管外给药次数 n 的对应关系，从而估算出期望血药浓度与周期性给药次数的数量对应。 积累系数 R 可由(1)的曲线在第一个周期末端的纵坐标与(1)的曲线在达坪状态一个周期内的最低点的纵坐标(即最低坪水平 Cmin )的比值计算出来。2 周期性血管外给药药物动力学参数的图形估算实例口服四
6、环素溶液 500mg 者进行血药浓度测定，10 人平均血药浓度如下1：取样时间 t(hour)1234681224 浓度 C(mg/ml).05504.78464.90764.46293.72492.50351.0096 若给药周期为每日 4 次，由残数法得到 n 次给药后药 时模型为：C(n)1(t)=1.930511-e-0.21849n1-e-0.21849e-0.21849+6.208371-e-0.07568n1-e-0.07568e-0.07568-8.138881-e-0.668n1-e-0.668e-0.668(4)在 MATLAB 下画出(4)和对应达坪分数 fss(n) 的
7、图像，如图 1 所示。图 1 药 时曲线，达坪分数周期性给药次数曲线由图 1 可见： 第 8 次周期性给药后血药浓度近似处于达坪状态，之后每个给药周期内的图像即为坪浓度 Css 与给药时间的对应关系。 最高坪水平 Cmax=YA=1414(mg/ml) ，最低坪水平Cmin=YB=1135(mg/ml) ，峰时 T0=xA-xB=1.97(h) 。 当给药次数为 4或 10 时，达坪分数分别为给药 fss(4)=YD=82.99% ，fss(10)=YE=98.89% 。 累积系数 R=CminYC=11.354.3152.63036 。 当每日给药次数改变为 3 次或 2 次，即给药周期为
8、8h 或 12h 时，从图 2 中的 F，G，H，L点的纵坐标能够得周期变化后的最高或最低坪浓度。图 3 是在峰时T0=1.97(h) ，周期性给药剂量 D0 不变的状况下，安全浓度 Cmax 和有效浓度 Cmin 与变量 的曲线，其函数表达式为：Cmax=1.930511-e-0.21849T01-e-0.21849+6.208371-e-0.07568T01-e-0.07568 -8.138881-e-0.668T01-e-0.668(5)Cmin=1.9305111-e-0.21849+6.2083711-e-0.07568-8.1388811-e-0.668(6)从图 3 可以看出，若
9、要求安全浓度 Cmax=13.11(mg/ml) ，则给药周期 xM=6.559(h) ;若要求有效浓度 Cmin=6(mg/ml) ，则给药周期 xN=9.678(h) 。图 2 8 小时和12 小时周期性给药的药 时曲线图 3 安全浓度 Cmax 和有效浓度Cmin 与变动周期 的曲线3 结语本研究探讨了借助图像估算周期性血管外给药状态下的药物动力学参数，该方法有以下特点： 由于 MATLAB 作图快捷、方便，使得药物动力学参数的估算简洁，高效，回避了繁难的数学计算。便于基层临床科研人员和医药院校师生借鉴与使用。 鉴于 MATLAB 图形界面中，按钮 Data Cursor 仅能显示到点的坐标的百分位，故此，估算参数会有不超过百分之一的误差。不过通过图形界面的图形放大功能，能得到万分之一精度的估算值。 本问题的处理方式可拓广到其它给药方式药物动力学参数的估算上。【参考文献】1 刘昌孝,刘定运.药物动力学概论.北京：中国学术出版社，1984，157169.2 梁文权,方晓玲，等.生物药剂学与药物动力学.北京：人民卫生出版社，2004，214219.3 M.吉伯尔迪, D.佩里尔. 药物动力学.北京:科学出版社，1981，127131.
