药物动力学模型 数学建模

药物动力学模型

一般说来，一种药物要发挥其治疗疾病的作用，必须进入血液，随着血流到达作用部位。药物从给药部位进入血液循环的过程称为药物的吸收，而借助于血液循环往体内各脏器组织转运的过程称为药物的分布。

药物进入体内以后，有的以厡型发挥作用，并以厡型经肾脏排出体外；有的则发生化学结构的改变--称为药物的代谢。代谢产物可能具有药理活性，可能没有药理活性。不论是厡型药物或其代谢产物，最终都是经过一定的途径(如肾脏、胆道、呼吸器官、唾液腺、汗腺等)离开机体，这一过程称为药物的排泄。有时，把代谢和排泄统称为消除。

药物动力学(Pharmacokinetics)就是研究药物、毒物及其代谢物在体内的吸收、分布、代谢及排除过程的定量规律的科学。它是介于数学与药理学之间的一门新兴的边缘学科。自从20世纪30年代Teorell 为药物动力学奠定基础以来，由于药物分析技术的进步和电子计算机的使用，药物动力学在理论和应用两方面都获得迅速的发展。至今，药物动力学仍在不断地向深度和广度发展。药物动力学的研究方法一般有房室分析；矩分析；非线性药物动力学模型；生理药物动力学模型；药物药效学模型。下面我们仅就房室分析作一简单介绍。

为了揭示药物在体内吸收、分布、代谢及排泄过程的定量规律，通常从给药后的一系列时间(t) 采取血样，测定血(常为血浆，有时为血清或全血)中的药物浓度( C )；然后对血药浓度——时间数据数

据(C ——t 数据)进行分析。

一 一室模型

最简单的房室模型是一室模型。采用一室模型，意味着可以近似地把机体看成一个动力学单元，它适用于给药后，药物瞬间分布到血液、其它体液及各器官、组织中，并达成动态平衡的情况。下面的图

(一)表示几种常见的给药途径下的一室模型，其中C 代表在给药后时间t 的血药浓度，V 代表房室的容积，常称为药物的表观分布容积，K 代表药物的一级消除速率常数，故消除速率与体内药量成正比，D 代表所给刘剂量。

图(a)表示快速静脉注射一个剂量D ，由于是快速，且药物直接从静脉输入，故吸收过程可略而不计；图(b)表示以恒定的速率K ，静脉滴注一个剂量D ；若滴注所需时间为丅，则K=D/丅。图(c)表示口服或肌肉注射一个剂量D ，由于存在吸收过程，故图中分别用F 和0K 代表吸收分数和一级吸收速率常数。

1. 快速静脉注射

在图(a)中所示一室模型的情况下，设在时间t ，体内药物量为x ，则按一级消除的假设，体内药量减少速率与当时的药量成正比，故有下列方程：

dx Kt dt

=- (5.1) 快速静脉注射 恒速静脉滴注 口服或肌肉注射 0K F 0K

图(一)

初始条件为t=0，x=0，容易解得

-Kt x=De ……………………..(5.2)

注意到房室的容积为V ，故c=x/V ；记t=0时血药浓度为0C ，因此0C =D/V ，则有

0C=C Kt e -…………………….(5.3)

这就是快速静脉注射(简称静注)一个剂量D 时，符合一室模型的药物及其血药浓度随时间递减的方程。对方程3两边取对数得

0ln ln C C Kt =-

这表明在一室模型的情况下，将实测的C_t 数据在以t 为横轴，ln C 为纵轴的坐标系上作图，各个数据点应呈直线散布趋势。据此，用图测法或最小二乘法拟合一条直线，其斜率为K ，截距为0ln C ，于是K 和0C 便可求得。当然，如果数据点的散布明显地不是呈直线趋势，则可断言不宜采用一室模型来解释该药物在快速静脉注射时的体内动力学过程。在实际应用中，表征药物消除快慢常用的参数是生物半衰期，记为1t /2，它是指药物浓度降至原定值的一半所需的时间。在方程(3)中令t=1t /2，C =0C /2，可得

1ln 20.6922t K K

=?………………(5.4)

可见半衰期是常数，且与消除速率常数成反比。

例如，给一名志愿者一次静脉注射某药物100mg ，测得给药后一些时刻的血药浓度见下表，和在坐标系上作出各数据点，它们是呈直线散布趋势，故可采用一室模型。

一次静注100mg 所得数据

如用最小二乘法拟合如下的直线方程

ln C a bt =+……………..(5.5)

利用实测的C 一t 数据计算直线斜率和截距的公式为：

1112211111ln ln 11ln n n n i i i i i i i n n i i i i n n i i i i t C t C n b t t n a C b t n =======ì骣骣ï鼢ï珑-鼢镧?鼢珑ï桫桫ïï=ï骣ïï÷ç-÷íç÷çï桫ïïï骣ï÷ïç=-÷ïç÷ç秣桫ïî邋?邋邋……………..(5.6)

其中n 为C 一t 数据点的个数。

将上表中的有关数据代入 (6) 式得

b=-0.02744 a=1.7386

于是，拟合数据点的直线方程为

lnC=1.7386-0.02744

与方程 (4) 对照，便得0C 和K 的估计值为

()()10 5.689/,0.0274C g ml K h m -==

进而，可得该药物的生物半衰期1/2t 和表观分布容积V 为

1/2t =()0.69325.30.0274

h = ()010017.65.689D V l C === 2.恒速静脉滴注

在图 (b) 所示一室模型的情况不，体内药量x 随时间t 变化的微分方程如下：

0dx K Kx dt =- …………… (5.7)

在初始条件t=0，x=0之下，可得其解为

()01Kt K x e K

-=- ………….. (5.8) 其中0,t T ＃，这里T 为滴注持续的时间。利用x=VC ，由 (8) 式得

()01Kt K C e VK

-=- ………… (5.9) 这就是恒速静脉滴注期间，符合一室模型的药物浓度随时间递增的方程。

假如t=丅时，所给剂量D 滴注完毕，则此后的血药浓度便按静注射时的规律下降 (如图 二)，