一室药物动力学模拟器 返回
医工学

一室药物动力学模拟器 — 静脉推注和口服给药

假设体内为一个均匀房室的最基本药物动力学(PK)模型。改变给药量、分布容积、消失速度常数、吸收速度常数,血中浓度时间特征、半衰期、清除率、AUC、Cmax/Tmax即时计算,直观理解静脉推注和口服给药的区别。

参数设置
给药途径
静脉推注:瞬间给药。口服:从消化道吸收。
给药量 D
mg
分布容积 V_d
L
药物均匀分布时的见表容积
消失速度常数 k_e
1/hr
代谢、排泄的一阶消失速度常数
吸收速度常数 k_a
1/hr
仅对口服有效。消化道→血液的转移速度
口服生物利用度 F
%
口服给药的全身循环到达率(含首过效应)
计算结果
半衰期 t₁/₂ (h)
清除率 CL (L/h)
最高浓度 C_max (mg/L)
峰值时间 T_max (h)
AUC₀₋∞ (mg·h/L)
12h血中浓度 (mg/L)
体内动态示意图 — 吸收·分布·消失

从给药途径流入分布容积,按消失速度常数ke向体外排出。中央脉冲表示当前血中浓度。

血中浓度 C(t) — 0~24小时
AUC₀₋∞ vs 给药量
理论·主要公式

$$\text{IV: } C(t)=\frac{D}{V_d}e^{-k_e t},\quad \text{Oral: } C(t)=\frac{F D k_a}{V_d(k_a-k_e)}(e^{-k_e t}-e^{-k_a t})$$

血中浓度时间特征。D:给药量mg、Vd:分布容积L、F:生物利用度、ka:吸收速度常数1/h、ke:消失速度常数1/h。

$$t_{1/2}=\frac{\ln 2}{k_e},\qquad CL=k_e\cdot V_d,\qquad \text{AUC}_{0-\infty}=\frac{F\cdot D}{k_e\cdot V_d}$$

半衰期、清除率、曲线下面积。IV时F=1。AUC代表全身暴露量,是疗效和毒性评估的基本指标。

$$T_{\max}=\frac{\ln(k_a/k_e)}{k_a-k_e}\quad(\text{口服仅有})$$

口服给药时的峰值到达时间。当ka → ke时可用洛必达法则(实现中以微小ε回避)。

一室药物动力学简介

🙋
老师,吃的药进了体内后怎么运动的呢?能像X光一样"这里还有3g药"这样测出来吗?
🎓
基本靠采血测浓度。但人体有脑、肌肉、脂肪等不同器官,分布速度也不一样,真要逐个建模就太复杂了。所以干脆把"整个人体看成一个大烧杯",给进去的药均匀溶解在分布容积Vd里,然后以ke的速度漏出去,这就是一室模型。
🙋
可是人体哪有那么简单呢?这样怎么能准确决定用药量呢…
🎓
好问题!其实对于组织间转移快的药物,比如咖啡因、乙醇、多数抗菌药,一室模型已经够用了。血中浓度的半对数图纸上成直线就说明符合。但是脑子里堆积的麻醉药、脂肪里储存的地西泮,就需要2室或3室模型。不过实际工作中80%的药都能用一室搞定。就因为简单,所以"半衰期多少?清除率多少?"这些问题马上就能答。
🙋
我注意到左边从"给药途径"切换静脉推注到口服,图形完全变了!静脉是一开始就最高,然后直线下降;口服是个山形。
🎓
对对对,这就是药物动力学的精妙之处。静脉推注一针下去药直接进血管,t=0就是Cmax=D/Vd,然后乖乖按ke指数衰减。口服得从胃肠吸收,最开始吸收>消失,浓度上升;某个时刻Tmax两者平衡,达到峰值;之后消失占优,浓度下降。这是双指数曲线。比如ka=1.0、ke=0.1,Tmax=ln(10)/0.9≈2.6小时。
🙋
生物利用度F是什么?口服的时候才出现的参数呢。
🎓
就是"吃进去的药有百分之几到达全身血液"的比例。有的在胃里溶不了,有的肠子吸不了,还有的被肝脏首过代谢干掉,所以百分比达不到100%。比如心得安,F才30%左右,所以注射和吃的用量差3倍多。反过来F很高(90%以上)的药,注射和吃的剂量就差不多。决定用量时,漏掉F这个参数那就是大事故。
🙋
AUC(曲线下面积)在统计学里也听过,在用药里也这么重要吗?
🎓
超级重要!AUC是"药在体内停留的浓度·时间总和",也就是全身暴露量。公式特别简洁:AUC=D/CL,跟给药量正比,跟清除率反比。比如D=500mg、ke=0.1/hr、Vd=30L,算出来AUC≈166.7mg·h/L。对抗癌药这种治疗窗口窄的,靶向AUC值来调剂量(比如卡铂要AUC=5)。后发医药品和原研药的同等性评价,也是靠AUC和Cmax在±20%以内。

常见问题

一室模型是最基本的药物动力学模型,假设体内全体是"均匀混合的单一容器(房室)"。给药的药物瞬间分布到分布容积Vd全体,以消失速度常数ke按一阶指数衰减。对于静脉推注,C(t)=D/Vd · e^(-ke·t),血中浓度的半对数图为直线。虽然实际人体有脑、脂肪、肌肉等分布速度不同的组织,但对于组织间转移快的药物(咖啡因、乙醇、多数抗菌药),一室模型近似充分。
半衰期定义为t₁/₂ = ln2 / ke,清除率定义为CL = ke · Vd。两者结合得t₁/₂ = ln2 · Vd / CL,说明半衰期与分布容积成正比,与清除率成反比。例如Vd=30L、ke=0.1/hr,则t₁/₂≈6.9hr、CL=3.0L/hr。肾功能低下患者CL下降,同样Vd下半衰期延长,需调整给药间隔或减量。清除率可理解为"单位时间内药物被清除的血液体积",代表肝代谢加肾排泄的总和。
口服给药中药物以吸收速度常数ka从消化道进入血流,以ke消失,形成"双指数模型"。浓度表达式为C(t) = F·D·ka / [Vd·(ka - ke)] · (e^(-ke·t) - e^(-ka·t))。Cmax到达时刻Tmax由dC/dt=0求解得Tmax = ln(ka/ke) / (ka - ke)。例如ka=1.0、ke=0.1时,Tmax=ln(10)/0.9≈2.56hr。F是口服生物利用度,包括消化道吸收率和首过效应损失。F低的药物(心得安F≈30%)需大幅减量相比注射剂。
AUC₀₋∞(血中浓度时间曲线从0到无穷的积分)代表"全身暴露量",是预测药物疗效和毒性的最重要指标。一室模型中AUC = D/(ke·Vd) = D/CL,AUC与给药量正比例,与清除率反比例。例如D=500mg、ke=0.1/hr、Vd=30L,则AUC≈166.7mg·h/L。口服给药中AUC = F·D/CL,已知F即可直接比较静脉和口服的AUC来评估生物利用度。抗癌药等治疗指数窄的药物,需调整给药量使AUC达到治疗目标值(如卡铂AUC=5)。

实际应用

抗菌药给药设计(TDM:治疗药物监测):万古霉素、氨基糖苷类抗菌药治疗窗口窄,需边测边调。用一室模型估计患者个体的Vd和ke,调整下次给药量和间隔,使AUC/MIC(如万古霉素400~600)在目标范围。肾功能用Cockcroft-Gault式或贝叶斯估计推算ke。

抗癌药AUC靶向给药(Calvert式):卡铂靶向AUC=5~7 mg·min/mL。Calvert式Dose = AUC × (GFR + 25)是把一室模型的AUC=D/CL用肾功能改写。按肾功能而非体表面积决定剂量,能更好控制血液毒性。至今还是铂类的标准做法。

后发药(学名药)的生物等效性试验:上市学名药需证明与原研药的AUC和Cmax在90% CI范围80~125%。约24~36人吃一遍、抽15左右次血,用一室或非房室法解析AUC。F差大了就有效果差异风险,不批上市。

新药研发I期试验(首次人体):第一次给人用新药时,从低剂量逐级递增,测血中浓度推估人的Vd、ke、CL、F。与动物数据对比验证,为后续试验定基础剂量。不符一室的进2/3室或群体PK(NONMEM)。

常见误解和注意

最大误解是把分布容积Vd当成真实体液量。Vd只是"从血浆浓度反算出的表观体积",和真实体液无关。成人总体液约42L,但Vd因药而异可以从5L(华法林,强结合白蛋白)到7000L(氯喹,强结合组织)差百倍。Vd超过体重的药物不是"血浆稀释"而是"积累在组织",解毒透析都没用,临床意义完全不同。本工具Vd范围1~200L是常用药的经验值。

次大误解是简单地说"半衰期长=一日一次就行"。虽然t₁/₂≈24h的定常态维持一天一次可以,但取决于"治疗窗有多窄、允许Cmax/Cmin变化多大"。治疗指数窄的药即使t₁/₂长也要分2~3次;反过来t₁/₂短的也能做缓释制一天一次。半衰期只是"达稳定所需时间(约4~5倍t₁/₂)"和"停药后残留时间"的参考,给药间隔要综合治疗窗、依从性、剂型来判断。

最后要注意一阶速论不是所有浓度都成立。本模型假设ke不随浓度变,但乙醇、苯妥英那样酶饱和的药就不行,要用Michaelis-Menten的零/混合阶。乙醇血浓无关约7g/h消失;苯妥英中毒稍加量浓度就爆升"饱和困局"。本工具是一阶消失的教学模拟器,饱和区的临床判断用不了。

使用指南

  1. 选择给药方式。在静脉推注(IV Bolus)和口服(Oral)间切换
  2. 设定给药量(Dose)、分布容积(Vd)、消失速度常数(ke)。口服时还要输入吸收速度常数(ka)和生物利用度(F)
  3. 按模拟运行,血中浓度曲线绘制,t₁/₂、CL、Cmax、Tmax、AUC₀₋∞、12h血中浓度自动计算

具体计算例子

茶碱(治哮喘的药):给药量500 mg、Vd=50 L、ke=0.11 h⁻¹、口服F=90%、ka=1.2 h⁻¹,算得t₁/₂=6.3h、CL=5.5 L/h、Cmax=6.8 mg/L(Tmax=2.1h)、AUC₀₋∞=82.6 mg·h/L。换静脉推注ka项消失,Cmax升到10.0 mg/L

实务注意项

  1. 华法林等肝代谢主体的药在肝功能低下患者,ke下降半衰期延长,需把ke改成0.05 h⁻¹级别重算
  2. 口服的Cmax受进食影响。高脂饮食下ka下降(如1.2→0.7 h⁻¹)、Cmax下降、Tmax延长,要反映这些变化验证
  3. 利用AUC₀₋∞=给药量×F/CL关系,可对不同生物利用度的剂型(如地高辛片F=70% vs 溶液F=90%)进行用量调整