关于用图像估算周期性血管外给药药物动力学参数重点

1、关于用图像估算周期性血管外给药药物动力学参数Throw away the apple because of the core.【摘要】 目的：直观、快捷地估算药物动力学参数的数值。方法：对周期 性血管外给药二室模型的药时曲线与参数的几何关系做了探讨，在MATLA的图形界面上，估算出了周期性血管外给药药物动力学参数的数值。结果：分别得 到了坪浓度、维持剂量、周期性给药的次数、达坪分数、积累系数等参数的数 值。结论：为药物动力学参数的估算提供了新的途径，该思路可拓广到同类问 题的估算上 () 。None but a wise man can employ leisure well.【关键词】 周期

2、性血管外给药 ; 药物动力学参数 ; 图像; 估算 利用数学模型确定药物动力学参数，需要做繁杂的数值计算，这对于基层医药 科研人员来说是一个不小的障碍。本文在分析模型参数几何意义的基础上，借 助MATLAB勺作图功能，探求周期性血管外途径给药二室模型药物动力学参数的 估算。Quality is better than quantity.1 周期性血管外给药二室模型与药物动力学参数的几何关系One cannot be in two place at once.1.1 周期性血管外给药二室模型及其参数周期性血管外 n 次给药后，体内血药浓度 C(n)1 与时间 t 的二室模型 为：C(n)1(t)

3、=A11-e- n a t 1-e- a t e- a t +A21-e-n B t 1-e-B t e- B t -(A1+A2)1-e- n K a t 1-e-K a t e-K a t (1)模型(1)中的t是从第n次给药后开始计算的时间(0 <t < t )，其它各系数的含义见 文献1 。若让模型 (1) 中的 n 趋向无穷，得到周期性血管外给药二室模型的 坪浓度：Css(t)=A111-e- a t e- a t +A211-e- B t e- B t -(A1+A2)11- e-K a t e-K a t由模型(1)、(2)在一个周期t内的平均值的比值得到 达坪分数：

4、 fss(n)=1 tJF(Dt 0C(n)1(t)dt JF)1t JF(Dt OCss(t)dt JF)=A1a (1-e- n a t )+A2 B (1-e- n B t )-+A2Ka(1-e- n Kat )A1 a +A2B -A1+A2Ka(3)最低坪浓度Cmin与第一次给药末端浓度 C1(t )的比是积累系数：R=A11-e-a t +A21-e- B t -A1+A21-e-K a tA1e- a 丫 +A2e-B t -(A1+A2)e-K a t 通常在已知周期性血管外给药二 室模型的情况下，根据上述各式和牛顿迭代法可以计算出坪浓度、维持剂量、 达坪状态时周期性给药的次

5、数等药物动力学参数。1.1 模型的曲线与药物动力学参数间的几何关联由于模型的函数为初等函数，故模型对应的曲线是连续、渐进变化的，不 会出现间断或剧烈波动的情况，从而为从图形中估算参数的数值提供了条件。 下面的几个参数与模型的曲线有如下关联：坪浓度Css是(1)的曲线稳定变化时一个周期内血药浓度与时间的对应关系。当给药次数n较大时，(1)的曲线在一个周期内所对应部分近似为坪浓度 (2) 的曲线(医药学/临床医学论文 ) 。 最低坪水平 Cmin 是(1) 的曲线在达坪状态一个周期内的最低点的纵坐 标，即达坪状态下一个周期内的最小值。最高坪水平Cmax是 的曲线在达坪状态一个周期内的最高点的纵坐标

6、，即达坪状态下一个周期内的最大 值。 峰时TmO是(1)的曲线在达坪状态下一个周期内最高点与最低点的横 坐标的差(0 当周期性给药剂量D0 一定时，改变给药周期 t,通过(1) 的曲线得到周期变化后的最高或最低坪浓度。反之，若给定安全浓度、有效坪 浓度，由 (2) 的曲线可以估算出符合条件的给药周期 。 由(3) 的曲线上点 的纵横坐标可以获得达坪分数 fss(n) 与周期性血管外给药次数 n 的对应关 系，从而估算出期望血药浓度与周期性给药次数的数量对应。积累系数R可由(1) 的曲线在第一个周期末端的纵坐标与 (1) 的曲线在达坪状态一个周期内 的最低点的纵坐标 (即最低坪水平 Cmin )

7、的比值计算出来。论文教育信息网 2 周期性血管外给药药物动力学参数的图形估算实例口服四环素溶液500mg者进行血药浓度测定，10人平均血药浓度如下1 ：取样时间 t(hour)1234681224 浓度 C(mg/ml).05504.78464.90764.46293.72492.50351.0096 若给药周期为每日 4 次，由残数法得到n次给药后药时模型为：C(n )1(t)=1.93051I-e-0.21849 n t 1-e- 0.21849 t e-0.21849 t +6.20837 l-e-0.07568 n t 1-e-0.07 568 - t e- 0.07568 t -8.

8、13888 l-e-0.668 n t 1-e- 0.668 t e-0.668 t 在 MATLA下画出和对应达坪分数 fss(n) 的图像，如图 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%。 累

9、积系数R=CminYC=11.354.31炭2.63036 。 当每日给药次数改变为 3次或2次， 即给药周期为8h或12h时，从图2中的F, G H L点的纵坐标能够得周期变 化后的最高或最低坪浓度。图 3是在峰时 T0=1.97(h) ，周期性给药剂量 D0 不变的状况下，安全浓度 Cmax 和有效浓度 Cmin 与变量 t 的曲线，其函数 表达式为：Cmax=1.93051l-e- 0.21849 T01-e- 0.21849 t +6.20837 l-e-0.07568 T01-e- 0.07568 t -8.13888 l-e- 0.668 T01-e-0.668 t (5)Cmin

10、=1.93051 11-e- 0.21849 t +6.20837 11-e- 0.07568 t - 8.13888 11-e- 0.668 t 从图3可以看出，若要求安全浓度Cmax=13.11(mg/ml)，则给药周期 t >xM=6.559(h);若要求有效浓度 Cmin=6(mg/ml)，则给药周期 t >xN=9.678(h)。图2 8小时和12小时周 期性给药的药时曲线图3安全浓度Cma和有效浓度Cmin与变动周期t 的曲线A faithful friend is hard to find.3 结语A thing is bigger for being shared.本研究探讨了借助图像估算周期性血管外给药状态下的药物动力学参数， 该方法有以下特点：由于MATLA作图快捷、方便，使得药物动力学参数的估算简洁，高效，回避了繁难的数学计算。便于基层临床科研人员和医药院校 师生借鉴与使用。 鉴于MATLABS形界面中，按钮Data Cursor仅能显示 到点的坐标的百分位，故此，估算参数会有不超过百分之一的误差。不过通过 图形界面的图形放大功能，能得到万分之一精度的估算值。本问题的处理方式可拓广到其它给药方式药物动力学参数的