《药物动力学模型数学建模》由会员分享，可在线阅读，更多相关《药物动力学模型数学建模（13页珍藏版）》请在金锄头文库上搜索。
1、1 药物动力学模型一般说来，一种药物要发挥其治疗疾病的作用，必须进入血液，随着血流到达作用部位。 药物从给药部位进入血液循环的过程称为药物的吸收，而借助于血液循环往体内各脏器组织转运的过程称为药物的分布。药物进入体内以后， 有的以厡型发挥作用， 并以厡型经肾脏排出体外；有的则发生化学结构的改变-称为药物的代谢。代谢产物可能具有药理活性，可能没有药理活性。不论是厡型药物或其代谢产物，最终都是经过一定的途径(如肾脏、胆道、呼吸器官、唾液腺、汗腺等)离开机体，这一过程称为药物的排泄。有时，把代谢和排泄统称为消除。药物动力学 (Pharmacokinetics) 就是研究药物、毒物及其代谢物在体内的吸
2、收、分布、代谢及排除过程的定量规律的科学。它是介于数学与药理学之间的一门新兴的边缘学科。自从 20 世纪 30 年代 Teorell为药物动力学奠定基础以来， 由于药物分析技术的进步和电子计算机的使用，药物动力学在理论和应用两方面都获得迅速的发展。至今，药物动力学仍在不断地向深度和广度发展。药物动力学的研究方法一般有房室分析；矩分析；非线性药物动力学模型；生理药物动力学模型；药物药效学模型。下面我们仅就房室分析作一简单介绍。为了揭示药物在体内吸收、分布、代谢及排泄过程的定量规律，通常从给药后的一系列时间(t) 采取血样，测定血 (常为血浆，有时为血清或全血 )中的药物浓度 ( C )；然后对血
3、药浓度时间数据数2 据(Ct 数据)进行分析。一一室模型最简单的房室模型是一室模型。采用一室模型， 意味着可以近似地把机体看成一个动力学单元，它适用于给药后， 药物瞬间分布到血液、其它体液及各器官、组织中，并达成动态平衡的情况。下面的图(一)表示几种常见的给药途径下的一室模型，其中 C 代表在给药后时间 t 的血药浓度， V 代表房室的容积，常称为药物的表观分布容积，K 代表药物的一级消除速率常数，故消除速率与体内药量成正比，D代表所给刘剂量。图(a)表示快速静脉注射一个剂量D，由于是快速， 且药物直接从静脉输入，故吸收过程可略而不计； 图(b)表示以恒定的速率K，静脉滴注一个剂量 D；若滴注
4、所需时间为丅，则K=D/丅。图 (c)表示口服或肌肉注射一个剂量D，由于存在吸收过程， 故图中分别用 F 和0K代表吸收分数和一级吸收速率常数。1. 快速静脉注射在图(a)中所示一室模型的情况下，设在时间t，体内药物量为 x，则按一级消除的假设， 体内药量减少速率与当时的药量成正比，故有下列方程：dxKtdt= -(5.1) 3 快速静脉注射恒速静脉滴注口服或肌肉注射0KF 0K图(一) 初始条件为 t=0，x=0，容易解得-Ktx=De(5.2)注意到房室的容积为V，故 c=x/V；记 t=0 时血药浓度为0C，因此0C=D/V，则有0C=CKte- .(5.3) 这就是快速静脉注射 (简称
5、静注 )一个剂量 D 时，符合一室模型的药物及其血药浓度随时间递减的方程。对方程3 两边取对数得0lnlnCCKt=-这表明在一室模型的情况下， 将实测的 C_t 数据在以 t 为横轴，ln C为纵轴的坐标系上作图，各个数据点应呈直线散布趋势。据此，用图测法或最小二乘法拟合一条直线，其斜率为K，截距为0ln C，于是K 和0C便可求得。当然，如果数据点的散布明显地不是呈直线趋势，则可断言不宜采用一室模型来解释该药物在快速静脉注射时的体内动力学过程。 在实际应用中， 表征药物消除快慢常用的参数是生物半（a）(b) （c）K K K V,C V,C V,C 4 衰期，记为1t/2，它是指药物浓度降
6、至原定值的一半所需的时间。在方程(3)中令 t=1t/2，C =0C/2，可得1ln 20.6922tKK=?(5.4) 可见半衰期是常数，且与消除速率常数成反比。例如，给一名志愿者一次静脉注射某药物100mg，测得给药后一些时刻的血药浓度见下表， 和在坐标系上作出各数据点， 它们是呈直线散布趋势，故可采用一室模型。一次静注 100mg所得数据t (h) C(/g mlm) lnC tlnC 2t0.5 2 3 6 12 24 5.52 5.42 5.32 4.80 4.10 2.94 1.7084 1.6901 1.6715 1.5686 1.4110 1.0784 0.8542 3.380
7、2 5.0144 9.4117 16.9318 25.8818 0.25 4 9 36 144 576 47.5 9.1280 61.4741 769.25 如用最小二乘法拟合如下的直线方程ln Cabt=+ (5.5) 利用实测的 C 一 t 数据计算直线斜率和截距的公式为：111 2211111lnln11lnnnniiii iiinnii iinnii iitCtCnbttnaCbtn=骣骣?鼢?珑-鼢镧?鼢珑?桫桫? =?骣?-?桫?骣?=-?秣桫?邋?邋邋(5.6) 5 其中 n 为 C 一 t 数据点的个数。将上表中的有关数据代入(6) 式得b=-0.02744 a=1.7386
8、于是，拟合数据点的直线方程为lnC=1.7386-0.02744 与方程 (4) 对照，便得0C和 K 的估计值为()()1 05.689/,0.0274Cg mlKhm-=进而，可得该药物的生物半衰期1/ 2t和表观分布容积 V 为1/ 2t=( )0.69325.3 0.0274h=( )010017.65.689DVlC=2.恒速静脉滴注在图 (b) 所示一室模型的情况不， 体内药量 x 随时间 t 变化的微分方程如下：0dxKKxdt=- (5.7) 在初始条件 t=0，x=0 之下，可得其解为()01KtKxeK-=- (5.8) 其中0,tT， 这里 T 为滴注持续的时间。利用 x
9、=VC， 由 (8) 式得()01KtKCeVK-=-(5.9) 这就是恒速静脉滴注期间， 符合一室模型的药物浓度随时间递增的方程。6 假如 t=丅时，所给剂量D 滴注完毕，则此后的血药浓度便按静注射时的规律下降(如图 二)，不过此时初始浓度为()01/KTKeVK-， 故滴注停止后的 C 一t 方程(为区别起见，特记为Ct-) 如下：()01KTKtKe CeVK- = (5.10) 由此可见，我们可以从滴注停止后测得Ct-数据，求得 K 和V 的估计值 (0K和丅皆已知 ) 假如滴注总是持续进行，则由(10)式可知，血药浓度将趋于一个极限，记作()001limKt ss tKKCe VKV
10、K-=-= (5.11) 这个血药浓度称为稳态浓度，又称坪水平。记在时刻t 的血药浓度达到坪水平的分数为ssf，则有1/ 20.693 11ttKt ssssCfeeC-=-=- (5.12) 7 可见达到稳态的快慢取决于消除速率常数K 或半衰期，与滴注速率 K 无关。例如，当滴注持续时间等于5 倍半衰期时，由 (12)式算得0.969ssf=，此时血药浓度约为坪水平彻97。3. 口服或肌肉注射在图（ c）所示一室模型的情况下，设在时刻t，体内药量为 x，吸收部位的药量为ax，则可建立如下的微分方程组aaa aadxK xKxdt dxK xdt?=-?-=? ?(5.13) 在初始条件 t=
11、0，ax=FD，x=0 之下，可解得()aK tKtaaK FDxeeKK-=-,(5.14) 从而血药浓度随时间变化的方程为()()aK tKtaaK FDCeeV KK-=-,(5.15) 令 M=()/aaK FD V KK-，则上式可写为()aK tKtCMee-=-,(5.16) 在通常情况下，吸收比消除快的多，即aKK，故对于足够大的 t，血药浓度实际上是时间的单项指数函数，为区别起见，记为KtCMe,(5.17) 或lnlnCMKt,(5.18) 8 据此可得 K 和 M 的估计值，然后计算足够大的t 之前各个实测浓度与按(5.17) 式推算的C*与 C 值之差称为“剩余浓度”r
12、C: aK t rCCCMe-*=-=,(5.19) 或l nl nraCMK t,(5.20) 据此可得 K 的估计值。上述这种估计消除和吸收速率常数的方法称为剩余法。(二) 二室型二室模型是从动力学角度把机体设想为两部分，分别称为中央室和周边室。中央室一般包括血液及血流丰富的组织(如心、肝、肾等)，周边室一般指血液供应少，药物不易进入的组织(如肌肉、皮肤、某些脂肪组织等 )。 在快速静注的情况下常见的二室模型如图4-2 所示。图中1V代表中央室的容积，10k代表药物从中央室消除的一级速率常数，12k和21k分别代表药物从中央室到周边室和反方向的一9 级转运速率常数，其余符号同前。设在时刻
13、t，中央室和周边室中的药物量分别为1x和2x，则可写出下列微分方程组：1 212121012 121212dxk xkkxdt dxk xk xdt,(5.14) 在初始条件120,0txD x之下，可解得2112 112 2ttttDkD kxeeDkxee,(5.15) 其中和由下列关系式决定：1012211021kkkk k,(5.16) 通常规定 。由于11xV C，故描述血药浓度随时间变化的方程为212111ttDkD kCeeVV,(5.17) 令211211/,/ADkVBD kV则有attCAeBe,(4.18) 10 根据(4.18)式，利用实测 Ct 数据，用剩余法或电子计
14、算机作曲线拟合，可得 、及 A、B 的值，而后按下列公式计算模型参数：1211021121021DVAB ABkABkkkkk,(4.19) 这组公式不难从 (4.17)、(4.18)式及 A、B 的定义导出。(三) 多次给药在临床药物治疗中绝大多数药物都需要多次给药，以使血药浓度在足够长的一段时间内处于安全，有效的治疗范围。因此，认识多次给药下血药浓度的变化规律是拟订合理的给药方案的基础。这里，我们只讨论一室模型多次重复静活的情况。假定某药在快速静注下，符合一室模型的动力学规律，那末，每隔一段时间，静注一个剂量D 时，血药浓度C 随时间 t 将如何变化呢？静注第一剂后， Ct 关系为0ktC
15、C e其中0/,0CD Vt，显然，最高浓度为0C，最低浓度为0kC e，记为11 10maxCC10minkCC e不难理解，静注第二剂后，则有2000max220minmax11KKKKKCCC eCeCCeCee静注 n 剂后，就有1 00max111nK nKK nKeCCeeCe,(5.21)0minmax1 1nK KK nnKeCCeCee,(5.22) 由此可知，重复静注n 剂后，血药浓度随时间的变化规律为01 1nK K nKeCCee0t,(5.23) 假如 n 充分大，使血药浓没达到稳态，那么，对(5.22) 式取 n的极限，使得稳态浓度的变化规律为0 1K KCCee,(5.24) 最高和最低稳态浓度分别为0 m a x11KKCDCeVe,(5.25) 00 min111K KKKCCDCeeeV e,(5.26) 在一个给药间隔时间内，平均稳在浓度为12 000111TT Kt ssKCDCC dtedteVK,(5.27) 图 4-4 表示每隔 6 小时重复静注一个剂量D 产生的 Ct 曲线最后， 我们举一个实例。卡那霉素的治疗血药浓度范围通常为10-25/ml。假定该药在其个病人的生物半衰期为3 小时，表观分布容积为15l，试问多次重复静注方案应该怎样？首先，注意到最高和最低稳态浓度依赖于给药方案(D 和)，两者之比为Ke，故有maxmin
