药物动力学模型

1、药物动力学模型一般说来，一种药物要发挥其治疗疾病的作用，必须进入血液，随着血流到达作用 部位。药物从给药部位进入血液循环的过程称为药物的吸收，而借助于血液循环往体内 各脏器组织转运的过程称为药物的分布。药物进入体内以后，有的以原型发挥作用，并以原型经肾脏排出体外；有的则发生 化学结构的改变-称为药物的代谢。代谢产物可能具有药理活性，可能没有药理活性。 不论是原型药物或其代谢产物，最终都是经过一定的途径(如肾脏、胆道、呼吸器官、 唾液腺、汗腺等)离开机体，这一过程称为药物的排泄。有时，把代谢和排泄统称为消 除。药物动力学(Pharmacokinetics)就是研究药物、毒物及其代谢物在体内的吸收

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

3、时间数据数据(Ct数据)进行分析。一 一室模型最简单的房室模型是一室模型。采用一室模型，意味着可以近似地把机体看成一个 动力学单元，它适用于给药后，药物瞬间分布到血液、其它体液及各器官、组织中，并 达成动态平衡的情况。下面的图(一)表示几种常见的给药途径下的一室模型，其中C代 表在给药后时间t的血药浓度，V代表房室的容积，常称为药物的表观分布容积，K代 表药物的一级消除速率常数，故消除速率与体内药量成正比，D代表所给刘剂量。图(a)表示快速静脉注射一个剂量D,由于是快速，且药物直接从静脉输入，故吸收 过程可略而不计；图(b)表示以恒定的速率K,静脉滴注一个剂量D；若滴注所需时间为 丁，则K=D

4、/T。图(c)表示口服或肌肉注射一个剂量D,由于存在吸收过程，故图中分 别用F和K。代表吸收分数和一级吸收速率常数。快速静脉注射在图(a)中所示一室模型的情况下，设在时间t,体内药物量为x,则按一级消除的 假设，体内药量减少速率与当时的药量成正比，故有下列方程：dx TOC o 1-5 h z = -Kt(5.1)dt初始条件为t=0, x=0,容易解x=De -Kt(5.2)注意到房室的容积为V,故c=x/V；记t=0时血药浓度为C0，因此C0=D/V，则有C=C eK(5.3)这就是快速静脉注射(简称静注)一个剂量D时，符合一室模型的药物及其血药浓度随时 间递减的方程。对方程3两边取对数得

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

6、一次静脉注射某药物100mg，测得给药后一些时刻的血药浓 度见下表，和在坐标系上作出各数据点，它们是呈直线散布趋势，故可采用一室模型。一次静注100mg所得数据t (h)C以g / ml)lnCtlnC120.55.521.70840.85420.2525.421.69013.3802435.321.67155.0144964.801.56869.411736124.101.411016.9318144242.941.078425.881857647.59.128061.4741769.25如用最小二乘法拟合如下的直线方程ln C U -bt(5.5)利用实测的C 一 t数据计算直线斜率和截距

7、的公式为:将上表中的有关数据代入(6)式得b=-0.02744a=1.7386于是，拟合数据点的直线方程为lnC=1.7386-0.02744与方程(4)对照，便得C0和K的估计值为C0 =5.689 Ng / ml , K= 0.0274 h 1进而，可得该药物的生物半衰期11/2和表观分布容积V为/2 =-0-69 = 25.3 h V =日= 17.6 l10.0274C 05.689恒速静脉滴注在图(b)所示一室模型的情况不，体内药量x随时间t变化的微分方程如下:0K (5.7)在初始条件t=0, x=0之下，可得其解为1 = % 1-et(5.8)其中。 tT,，这里T在初始条件t=

8、0, x=0之下，可得其解为1 = % 1-et(5.8)其中。 tT,，这里T为滴注持续的时间。利用x=VC,由（8）式得K- 0 1VK(5.9)e Kt这就是恒速静脉滴注期间，符合一室模型的药物浓度随时间递增的方程。假如t=T时，所给剂量D滴注完毕，则此后的血药浓度便按静注射时的规律下降（如图二）,不过此时初始浓度为KeKT /VK，故滴注停止后的C0一 t方程（为区别起见，特记为c t）如下:(5.10)由此可见，我们可以从滴注停止后测得c,_t,数据，求得K和V的估计值（K0和T 皆已知）假如滴注总是持续进行，则由（10）式可知，血药浓度将趋于一个极限，记作c =lim Kri-e=

9、Kt(5.11)(5.12)这个血药浓度称为稳态浓度，又称坪水平。记在时刻t的血药浓度达到坪水平的分(5.12)f =C = 1-eK= 1-e 饥 s css可见达到稳态的快慢取决于消除速率常数K或半衰期，与滴注速率K无关。例如，当滴注持续时间等于5倍半衰期时，由（12）式算得=0969，此时血药浓度约为坪水平fss彻 97%。口服或肌肉注射4(5.13)在图（C）所示一室模型的情况下，设在时刻t，体内药量为x，吸收部位的药量为a， 则可建立如下的微分方程组(5.13)匕=K4 一 K dx ” 了 - K 4在初始条件t=0，工 =FD，x=0之下，可解得 a4 = a e Kt e Ka

10、K -Ka从而血药浓度随时间变化的方程为K FDC ae Kt K FDC ae Kt e 膈V K-Ka(5.15)令M= K FD /V K _K，则上式可写为=M ef 厂喝=M ef 厂喝(5.16)在通常情况下，吸收比消除快的多，即，故对于足够大的血药浓度实际KaK上是时间的单项指数函数，为区别起见，记为C*= Me- Kt (5.17)或 ln C *= In M - Kt (5-18)据此可得K和M的估计值，然后计算足够大的t之前各个实测浓度与按(5.17)式推算的C *与C值之差称为“剩余浓度” C : rc = c* c mKa (5-19)或 ln C = ln M - K

11、 t (5-20)据此可得K的估计值。上述这种估计消除和吸收速率常数的方法称为剩余法。(二)二室型二室模型是从动力学角度把机体设想为两部分，分别称为中央室和周边室。中央室 一般包括血液及血流丰富的组织(如心、肝、肾等)，周边室一般指血液供应少，药物不 易进入的组织(如肌肉、皮肤、某些脂肪组织等)。在快速静注的情况下常见的二室模型 如图4-2所示。图中V代表中央室的容积，k 10代表药物从中央室消除的一级速率常数，k12和k21分别 代表药物从中央室到周边室和反方向的一级转运速率常数，其余符号同前。设在时刻t，中央室和周边室中的药物量分别为尤和尤，则可写出下列微分方程 TOC o 1-5 h z

12、 12组：(5.14)(5.15)= k x 一( k + k )x(5.14)(5.15)dt 21 212101dx丁 = k x 一 k xdt 12 121 2在初始条件n n 。之下，可解得 t = 0, x = D , x = 0 x1x2D(a-k )D(k x1x2a 13a 13Dk (=12 e-P t ea ta-p其中a和B由下列关系式决定:(5.16)a + P= k + k + k aP = k k(5.16)通常规定aB。D(a-kD(a-k )C =7-2+ e-at 4U (a-B)i(5.17)D(k -B) na1e-Bt 1(5.17)令 A=D(a-

13、k2i)/(a-B), B=D(k -B)/Vx (a-p)则有C = Ae-t + Be-Bt 18)根据(4.18)式，利用实测Ct数据，用剩余法或电子计算机作曲线拟合，可 得a、B、及A、B的值，而后按下列公式计算模型参数：匕=TTBk = A B+ B a (4.19)21 A + Bk =也10k 21k 12 = a + p - k 1o - k 21这组公式不难从(4.17)、(4.18)式及A、B的定义导出。(三)多次给药在临床药物治疗中绝大多数药物都需要多次给药，以使血药浓度在足够长的一段 时间内处于安全，有效的治疗范围。因此，认识多次给药下血药浓度的变化规律是拟订 合理的给

14、药方案的基础。这里，我们只讨论一室模型多次重复静活的情况。T假定某药在快速静注下，符合一室模型的动力学规律，那末，每隔一段时间静注一个剂量D时，血药浓度C随时间t将如何变化呢？静注第一剂后，Ct关系为C = C o e - kt其中 C = D / V ,0 t 其中 C = D / V ,0 t t(C 1 )max = C 0不难理解，静注第二剂后，(C )= C + C e-Kt(C2 )min =(C2 )max 静注n剂后，就有显然，最高浓度为C ，最低浓度为C(C)= C e k t则有=C (1 + e - k t )=Co (1 + e - kt )e - kt(c )n ma

15、x=C|_1+ e - kt + e -(n-1)kt(CnLn = C=C011 一 e - ktr1-e-nKt、e-反11-e-KT J(5.22)(5.21)由此可知，重复静注n剂后，血药浓度随时间的变化规律为(5.23)假如n充分大，使血药浓没达到稳态，那么，对(5.22)式取n8的极限，使得稳态浓 度的变化规律为(C(C 、C 日 TTTKT Je-K t(5.24)(5.25)最高和最低稳态浓度分别为(5.25)DV (1 - e - K t )(C )8 minCe-Kt =eKT -1 V (C )8 minCe-Kt =eKT -1 V Ckt 一 1(5.26)在一个给药间隔时间内，平均稳在浓度为(5.27)1 e Kdt = 4.1 - e - k T JVK T图4-4表示每隔6小时重复静注一个剂量D产生 最后，我们举一个实例。卡那霉素的治疗血药浓度范围通常为10-25%/ml。假定该药在 其个病人的生物半衰期为3小时，表观分布容积为15l，试问多次重复静注方案应该怎 样？ss首先 有注意到最高和最低稳态浓度依赖于给药方案(D和T)，两者之比为eKT，故ln( C )max(C