try ai
科普
编辑
分享
反馈
  • 基于生理的药代动力学(PBPK)建模

基于生理的药代动力学(PBPK)建模

SciencePedia玻尔百科
核心要点
  • PBPK模型通过将身体表示为一个由相互连接的器官组成的系统来模拟药物的体内过程,该系统受生理参数和质量守恒定律支配。
  • 模型的预测由关键的机理驱动,如非结合药物假说,该假说认为只有游离药物具有活性并能够穿过细胞膜。
  • 通过体外到体内的外推(IVIVE),PBPK建模将实验室测量数据转化为对药物在整个人体内行为的预测。
  • PBPK通过整合关于生理(如器官功能)、遗传(药物遗传学)乃至微生物组的个体数据来定制药物剂量,从而实现个性化医疗。
  • 该框架可以创建“虚拟人群”以进行计算机模拟临床试验,从而在人体测试前预测一系列反应并识别高风险个体。

引言

了解药物在人体内的作用是医学领域最关键的挑战之一。传统方法通常依赖于抽象的统计模型,但如果我们能够构建一个虚拟患者——一个尊重身体实际解剖学和生理学的数学模型呢?这正是基于生理的药代动力学(PBPK)建模的核心思想。它提供了一个机理性和预测性的框架,区别于那些生物学细节较少的方法,使科学家和临床医生能够以惊人的准确性模拟药物在体内的过程。PBPK填补了简单实验室实验与复杂人体反应之间的关键知识鸿沟,为开发更安全、更有效的药物提供了强大的工具。

本文将引导您进入PBPK建模的世界。在第一部分​​原理与机制​​中,我们将解构虚拟人体,探讨构成模型基础的质量平衡、器官房室、药物分配和清除率等基本概念。随后,在​​应用与跨学科联系​​部分,我们将看到这个强大的引擎如何实际运作,探索它如何用于将实验室数据外推至人体预测,为独特个体实现个性化医疗,并为孕妇等特殊人群提供关键见解,从而变革药物开发和临床实践。

原理与机制

想象一下,您是一位物流专家,任务是预测一百万个包裹在一座广阔而复杂的城市中的运送轨迹。您不会将这座城市视为一个单一、均质的盒子,而是会拿出一张地图。您会识别出主要的区域(器官)、连接它们的公路和街道(循环系统)、每条路线上的交通流量(血流量),以及每个地址的特定投递规则(细胞屏障和代谢酶)。然后,您将基于一条简单而不可违背的定律来追踪这些包裹:包裹守恒。进入一个区域的包裹,要么留在那里,要么在那里被处理,要么离开。

基于生理的药代动力学(PBPK)建模做的正是这件事,只不过这里的“包裹”是药物分子,而“城市”则是一个生命有机体——无论是一只小鼠还是一个人。它是解剖学、生理学和数学的美妙结合,创建了一个动态的、机理性的药物体内过程图谱。与更为抽象的方法不同,PBPK以真实、可测量的生物学特性为基础。其核心是关于运动和转化的故事,并由​​质量守恒​​这一基本原则所支配。

作为充分搅拌容器的器官

让我们从一个单一的区域或器官——比如肌肉——开始构建这张图谱。我们可以将其想象成一个简单的桶。一根代表动脉血液供应的软管以该器官的血流量 QmQ_mQm​ 决定的速率,持续向桶中注入药物溶液。一根排水管,即静脉回流,以相同的速率将溶液排出。桶中药物的量 AmA_mAm​ 根据一个简单的平衡关系发生变化:

变化率 = 输入率 - 输出率

输入率很简单:即血流量乘以动脉血中的药物浓度,Qm⋅CartQ_m \cdot C_{\text{art}}Qm​⋅Cart​。但输出率呢?要弄清楚这一点,我们需要做出第一个也是最常见的简化假设:​​充分搅拌模型​​(well-stirred model)。我们想象药物一进入肌肉,就会瞬间且完美地混合在整个组织体积 VmV_mVm​ 中。因此,肌肉内的药物浓度 Cm=Am/VmC_m = A_m / V_mCm​=Am​/Vm​ 在各处都是均匀的,包括即将通过静脉流出的血液中。

但等一下——静脉血中的浓度就等于组织中的浓度吗?不完全是。药物有自己的偏好。一种油性的、亲脂性的药物可能会觉得组织的脂肪环境远比血液的水性血浆更“舒适”。它会积聚或​​分配​​(partition)到组织中。我们用​​组织-血浆分配系数 KpK_pKp​​​ 来量化这种偏好。KpK_pKp​ 为10意味着在平衡状态下,组织中的浓度将是血浆中浓度的十倍。

因此,在充分搅拌的假设下,离开组织的静脉血与组织本身处于平衡状态。其浓度 CvenousC_{\text{venous}}Cvenous​ 与组织浓度的关系为 Cvenous=Cm/KpC_{\text{venous}} = C_m / K_pCvenous​=Cm​/Kp​。现在,我们用于肌肉的质量平衡方程就变成了一个强大的常微分方程:

VmdCmdt=Qm(Cart−CmKp)V_m \frac{dC_m}{dt} = Q_m \left( C_{\text{art}} - \frac{C_m}{K_p} \right)Vm​dtdCm​​=Qm​(Cart​−Kp​Cm​​)

这个优雅的方程是PBPK建模的基石。它告诉我们,基于四个关键参数:器官体积(VmV_mVm​)、血流量(QmQ_mQm​)、药物的动脉浓度(CartC_{\text{art}}Cart​)以及药物对该器官的亲和力(KpK_pKp​),器官中的药物浓度如何随时间变化。通过为每个主要器官编写类似的方程,并通过循环系统将它们全部连接起来(其中动脉浓度是离开所有其他器官的血液混合后的结果),我们就创建了一个完整的、动态的身体图谱。

当桶模型不够用时:灌注与渗透

充分搅拌模型非常简洁,但自然界钟爱复杂性。该模型的核心假设是,药物从血管进入组织细胞毫无困难——该过程仅受限于血液输送药物的速度。这被称为​​灌注限制​​(perfusion-limited)分布。对于许多小的、亲脂性的药物来说,这是一个完全合理的近似。

但如果药物分子很大,或带电荷,并发现细胞膜是一个难以逾越的屏障呢?在这种情况下,瓶颈就不是药物的递送,而是穿过毛细血管壁的缓慢过程。这被称为​​渗透限制​​(permeability-limited)分布。此时,关键在于由血流量 QQQ 决定的对流递送速率与由渗透性-表面积乘积 PSPSPS 表征的跨膜扩散转运能力之间的较量。

  • 如果 PS≫QPS \gg QPS≫Q,细胞膜就像一扇敞开的门。转运是灌注限制的。桶模型成立。
  • 如果 PS≪QPS \ll QPS≪Q,细胞膜就像繁忙时段的旋转栅门。进入速率受限于栅门,而非人群。药物在血管中的浓度可能很高,但在组织中仍然很低。

在渗透限制的情况下,我们的单桶模型就不再足够。我们必须完善我们的图谱,将器官表示为至少两个房室:毛细血管内的血管空间和血管外的组织空间。这两个房室现在通过一个通量项 PS⋅(浓度差)PS \cdot (\text{浓度差})PS⋅(浓度差) 连接起来,该项明确描述了跨越屏障的缓慢转运过程。PBPK框架的精妙之处就在于,它能够根据药物和组织的基础物理和生物学特性,选择合适的细节层次。

非结合药物的秘密

我们已经讨论过分配系数 KpK_pKp​ 作为药物对组织偏好的度量。但什么决定了这种偏好呢?答案在于药理学最基本的原则之一:​​非结合药物假说​​(unbound drug hypothesis)。

想象一下血液中的药物分子。它们中的许多并非“自由”的;它们与白蛋白等大分子血浆蛋白结合,就像游轮上的乘客。同样,在组织内,药物可以与细胞蛋白和脂肪结合。非结合药物假说指出,只有“自由”的或​​非结合​​的药物分子——那些独自游动的分子——能够穿过细胞膜,与靶点相互作用,并被清除。结合的药物只是随波逐流。

这个简单的想法具有深远的影响。在平衡状态下,跨越血浆-组织屏障达到平衡的是非结合药物的浓度。让我们将血浆中的非结合药物分数表示为 fu,pf_{u,p}fu,p​,组织中的表示为 fu,tf_{u,t}fu,t​。然后我们可以写出:

总血浆浓度,Cp,total=非结合血浆浓度, Cp,unboundfu,pC_{p, \text{total}} = \frac{\text{非结合血浆浓度, } C_{p, \text{unbound}}}{f_{u,p}}Cp,total​=fu,p​非结合血浆浓度, Cp,unbound​​

总组织浓度,Ct,total=非结合组织浓度, Ct,unboundfu,tC_{t, \text{total}} = \frac{\text{非结合组织浓度, } C_{t, \text{unbound}}}{f_{u,t}}Ct,total​=fu,t​非结合组织浓度, Ct,unbound​​

在平衡状态下,Cp,unbound=Ct,unboundC_{p, \text{unbound}} = C_{t, \text{unbound}}Cp,unbound​=Ct,unbound​。如果我们现在回顾分配系数的定义,Kp=Ct,total/Cp,totalK_p = C_{t, \text{total}} / C_{p, \text{total}}Kp​=Ct,total​/Cp,total​,一个异常简洁的关系便应运而生:

Kp=fu,pfu,tK_p = \frac{f_{u,p}}{f_{u,t}}Kp​=fu,t​fu,p​​

药物在组织和血浆之间的整体分配,竟然只是两个房室中非结合分数之比! 这个方程将蛋白质结合的微观世界与药物分布的宏观行为联系了起来。

但是什么决定了这些非结合分数呢?这是药物特性和组织环境共同作用的结果。亲脂性(通过 log⁡P\log PlogP 测量)发挥了作用,因为油性药物喜欢溶解在细胞脂肪中。但对于许多弱酸或弱碱性药物而言,它们的电荷状态至关重要。细胞内的pH值(约7.0)比血浆中(7.4)略酸。对于弱碱性药物,这种较低的pH值使其更容易质子化(即带正电荷)。这种带电形式会被“困”在细胞内,难以再次穿过细胞膜出去。它还可以与细胞中带负电的成分(如酸性磷脂)强力结合。这种现象被称为​​离子捕获​​(ion trapping),可导致碱性药物在组织中大量积聚,这是仅靠亲脂性无法预测的。机理PBPK模型,例如使用Rodgers-Rowland方法的模型,明确地考虑了这些pH梯度和静电相互作用,从而可以从第一性原理预测 KpK_pKp​。

身体的清除策略:清除率

药物的体内过程直到被清除才算完成。肝脏是身体主要的代谢工厂,通过化学修饰药物来促进其排泄。在我们的PBPK图谱中,我们必须为肝脏房室配备一个清除功能。

肝脏固有的代谢能力由​​内源性清除率 CLintCL_{int}CLint​​​ 来描述。它代表了肝酶在不考虑药物递送速度的情况下,处理非结合药物的最大速率。然而,实际观测到的​​肝清除率 CLhCL_hCLh​​​ 是这种内源性能力与将药物递送至酶的血流量 QhQ_hQh​ 相结合的产物。

正如我们对组织房室结构的讨论一样,药理学家们长期以来一直在争论模拟肝脏的最佳方式。它是一个单一的“充分搅拌”的容器,其中酶所接触到的浓度是离开肝脏的血液中的低浓度吗?还是一个“平行管”系统,其中药物浓度在流经一排排酶时逐渐下降? 每种模型都给出了关于 CLhCL_hCLh​ 如何依赖于 QhQ_hQh​ 和 fu⋅CLintf_u \cdot CL_{int}fu​⋅CLint​ 的不同数学表达式,而它们之间的选择取决于具体的药物和生理环境。

有了这些组件,我们就可以回答一些关键问题。对于以恒定速率 RinR_{in}Rin​ 输注的药物,其长期浓度会是多少?在稳态下,药物进入身体的速率必须等于其被清除的速率。总清除速率是全身总清除率(CLtotCL_{\text{tot}}CLtot​,即肝脏、肾脏等清除率的总和)乘以稳态血浆浓度 Cp,ssC_{p,ss}Cp,ss​。因此,在一个异常简洁的平衡关系中:

Cp,ss=RinCLtot=RinCLh+CLr+…C_{p,ss} = \frac{R_{in}}{CL_{\text{tot}}} = \frac{R_{in}}{CL_h + CL_r + \dots}Cp,ss​=CLtot​Rin​​=CLh​+CLr​+…Rin​​

PBPK建模使我们能够从头构建 CLhCL_hCLh​ 和 CLrCL_rCLr​ 的表达式,将生理学与这一关键临床指标直接联系起来。

超越小分子:一个适用于所有药物的框架

PBPK框架的力量和精妙之处在于其适应性。它不仅限于小而简单的分子。以一种大的治疗性单克隆抗体(mAb)为例。这种150 kDa的蛋白质的行为与阿司匹林截然不同。它太大,难以轻易地通过扩散穿过细胞膜;相反,它是通过​​对流​​被拖入组织的。它也太大,无法通过扩散回到血液中;相反,它通过缓慢、平稳的​​淋巴系统​​回流返回循环。它的清除也不是通过简单的代谢。它受到一种名为​​FcRn​​的特殊再循环受体的保护,免于降解。我们可以将所有这些独特的规则整合到我们的PBPK图谱中。房室和连接保持不变,但转运定律和清除函数被更新,以反映新分子的独有生物学特性。质量平衡的基本原则对所有情况都适用。

从个体到群体:虚拟人群与不确定性

到目前为止,我们已经为单个、平均的个体构建了一张图谱。但医学的最终目标是治疗真实、多样化的人群。这正是PBPK建模真正大放异彩的地方。它使我们能够区分​​变异性​​(variability)——个体之间真实存在的内在差异——和​​不确定性​​(uncertainty)——我们自身对模型参数知识的不完善。

  • ​​变异性​​是偶然的(aleatory,“像掷骰子”)。正是这种生物学上的随机性使您与邻居有所不同。PBPK模型通过将参数与年龄、体重、性别和器官功能等个体​​协变量​​(covariates)联系起来,来捕捉这种变异性。
  • ​​不确定性​​是认知的(epistemic,“与知识相关”)。这是我们参数估计中的模糊性,可以通过收集更多数据来减少。

通过描述目标人群中协变量的分布(例如,根据人口普查数据),我们可以使用PBPK模型生成一个​​虚拟人群​​——一个由数千个独特的、生理上一致的虚拟人组成的计算机模拟(in silico)队列。然后,我们可以对这个人群进行“虚拟临床试验”,不仅可以预测药物的平均反应,还可以预测完整的反应范围,从而在真实试验开始之前就识别出可能有毒性或疗效不佳风险的个体。

最后,面对如此复杂的模型,我们如何知道其众多参数中哪些最重要呢?我们通过​​敏感性分析​​来探究模型。我们可以“微调”参数,观察输出如何变化。例如,​​全局敏感性分析​​可以将预测中的总不确定性(如患者的药物暴露量)进行分解,并将其归因于每个输入参数的不确定性。这告诉我们哪些生理因素是药物行为的真正驱动因素,从而指导未来的研究,帮助我们建立更稳健、更可靠的模型。

从一个简单的桶到一个虚拟人,PBPK建模的历程证明了运用第一性原理来建立对生物学的量化和预测性理解的强大力量。这是一张不断完善的图谱,但它已经指导我们为全人类开发更安全、更有效的药物。

应用与跨学科联系

在我们之前的讨论中,我们煞费苦心地构建了基于生理的药代动力学(PBPK)模型的蓝图。我们看到,通过尊重身体自身的结构——其器官、血流和生化机制——我们如何能够构建一个生命体的数学画像。您可能会说,这是一个优美的理论练习,但它有用吗?我们能用这个“虚拟患者”来做些什么呢?

事实证明,答案的广度惊人。本章将带领我们踏上PBPK建模应用的旅程,从化学家的实验台到患者的床边,从遗传密码到我们肠道内的微生物宇宙。我们将看到,我们的PBPK模型不仅仅是一个描述性工具,它更是一个预测引擎,一个用于药物和化学品的“飞行模拟器”,让我们能够探索、预测和理解分子与人体之间错综复杂的相互作用。

从试管到人体:IVIVE的魔力

医学和毒理学面临的最大挑战之一,是弥合实验室与生命有机体之间的鸿沟。我们可以轻易地在细胞培养皿中测试新药的效力,或在试管中测量肝酶分解化学物质的速度。但人不是一个巨大的培养皿。人有泵血的心脏、过滤的肾脏和储存物质的脂肪。我们如何可能仅凭一个简单的体外(in vitro,字面意思为“在玻璃中”)实验,就预测一种化学物质在整个人体内的命运呢?

这正是体外到体内外推(IVIVE)的魔力所在,它是PBPK建模的一个基石应用。PBPK模型扮演着终极“放大器”的角色。我们首先在实验室中测量一个基本参数,例如内源性清除率(CLintCL_{int}CLint​),它衡量了一小部分肝酶样本代谢药物的效率。然后,我们的PBPK模型,凭借其对标准人体生理学的知识——平均肝脏大小、肝细胞中的蛋白质含量以及流向肝脏的血流量——将这一微观测量值进行放大。它计算出整个器官的清除率会是多少,从而在我们对人体用药之前,就给出了关于人体将如何处理该药物的有力预测。

同样的原理也让我们能够解决安全和风险问题。想象一下,在基于细胞的检测中,发现某种环境化学品在特定浓度 AC50AC_{50}AC50​ 下会干扰激素信号传导。为了评估对人类的风险,我们需要知道现实世界中多大的暴露水平会在体内产生相同的活性浓度。从外部剂量到内部浓度的过程对化学品而言充满障碍:它必须被吸收,在穿过肠壁和肝脏时存活下来,并应对与血浆蛋白的结合。PBPK模型模拟了这整个过程,考虑了实验室检测与人体血液中蛋白质结合的差异,并使用“充分搅拌”肝脏模型等原理来模拟药物的清除。通过反向运行模拟,PBPK帮助我们确定“人体等效剂量”——即可能导致靶点处产生潜在有害浓度的外部暴露量。这对于监管机构为从杀虫剂到塑料等各种物质设定安全限值来说,是一个不可或缺的工具。

个性化医疗的艺术

或许PBPK最激动人心的前沿领域,是其在开创真正个性化医疗中的作用。“平均人”是一个有用的统计学构想,但在临床上,我们治疗的是个体,每个个体都有其独特的生理、基因构成和生活方式。PBPK模型因其本质,就是为了容纳这种个体性而构建的。

根据生理状况定制剂量

考虑一种能被肝脏高效清除的药物——即所谓的“高提取”药物。对于这类药物,限制其从体内清除的主要因素不是肝脏的代谢能力,而仅仅是血液将药物输送到肝脏的速率。其清除变为血流限制(flow-limited)。现在,想象一位患有中度心力衰竭的患者,其心脏泵血效率较低,导致心输出量下降。对于PBPK模型而言,这并非一个不可预见的并发症,而只是一个输入参数的简单变化。模型会立即预测,流向肝脏的血流量减少将降低药物的清除率。因此,为维持理想的治疗浓度,这位患者需要的输注速率要低于心脏健康的人。这不是猜测,而是基于系统基本机理的量化预测。

同样,模型也可以解释身体组成的差异。肥胖个体的脂肪组织体积远大于消瘦个体。对于一种亲脂性(即“喜好脂肪”)的药物,这部分额外的脂肪组织就像一个巨大的储库,它会储存药物,并将其缓慢释放回循环中。PBPK模型包含了独立的脂肪和肌肉房室,并有各自的体积和血流量,因此可以精确模拟这一现象。它能预测肥胖将如何改变药物的分布,可能会延长其药效,并为制定不同的给药策略提供信息。

从基因到药物反应:药物遗传学

药物反应的个体差异比解剖结构更深层——它写在我们的DNA中。许多负责代谢药物的酶,如著名的细胞色素P450家族,是由在人群中存在变异的基因编码的。一些人拥有导致酶活性极高的基因变体(“超快代谢者”),而另一些人的基因变体则产生缓慢或无功能的酶(“慢代谢者”)。

对于一个超快代谢者,标准剂量的药物可能会被迅速清除,以至于永远达不到治疗水平。而对于一个慢代谢者,相同剂量则可能累积到毒性水平。PBPK建模为整合这些遗传信息提供了完美的框架。通过根据患者已知的基因型调整肝脏房室中的 CLintCL_{int}CLint​ 参数,模型可以预测由此产生的总药物暴露量(通常以曲线下面积,即AUC衡量)的差异。这使得临床医生能够超越“一刀切”的方法,根据患者的基因图谱预先调整剂量,从而在最大化疗效的同时最小化风险。

体内宇宙:人类微生物组

当我们考虑到我们并非独自存在时,个体性的故事又展开了引人入胜的新篇章。我们的身体是数万亿微生物的家园,尤其是在我们的肠道中,它们形成了一个被称为微生物组的复杂生态系统。这些细菌拥有庞大的酶库,其中一些可以代谢我们人类细胞无法代谢的药物。

对于一种口服药物,它从肠道吸收后面临的第一个挑战,是在到达体循环之前,必须经受住肠壁和肝脏“首过代谢”的考验。肠道细菌在这一过程中的贡献长期以来一直是个谜。随着宏基因组学的出现,我们现在可以对一个人微生物组的集体DNA进行测序,并量化特定细菌基因的丰度。PBPK建模可以将这些信息直接与药代动力学联系起来。通过将肠壁代谢率设定为特定细菌基因丰度的函数,模型可以预测患者独特的微生物组将如何影响药物的口服生物利用度。这为未来通过一份粪便样本就能帮助确定适合患者的药物和剂量打开了大门。

窥探特殊世界

当我们需要了解药物在那些难以、不道德或不可能直接进行实验的情况下的行为时,PBPK的力量就愈发凸显。

怀孕:同时保护两位患者

治疗孕妇意味着要治疗两位患者:母亲和发育中的胎儿。预测药物的安全性和有效性需要了解它如何穿过胎盘屏障以及在胎儿体内达到何种浓度。这是一个临床试验受到严格限制的领域。PBPK建模提供了一种强大且符合伦理的替代方案。母-胎PBPK模型通过增加胎盘和胎儿房室来扩展标准框架,其转运速率反映了怀孕期间独特的生理状况。这些模型可以针对不同的孕期进行调整,考虑血容量、心输出量和胎盘血流量的变化。它们可用于预测胎儿暴露量,甚至估算药物在脆弱的胎儿大脑中的分配情况,为指导孕期治疗提供关键见解。

从血药浓度到药物效应:PK/PD的联系

到目前为止,我们一直专注于药代动力学(PK)——即身体对药物的作用。但我们最终关心的是药效动力学(PD)——即药物对身体的作用。PBPK模型为建立这种联系提供了完美的基础。

以异丙酚(propofol)这样的麻醉剂为例。对于麻醉而言,重要的不仅是血液中的浓度,还有“效应部位”——大脑中 GABAAGABA_{\text{A}}GABAA​ 受体——的浓度。PBPK模型可以包含一个详细的大脑房室,模拟血流量以及药物从血液分配到脑组织的速率。通过这样做,它可以预测药物在作用部位的浓度时程。这使我们能够计算出具有临床意义的参数,例如达到稳态效应部位浓度90%所需的时间,从而为麻醉师提供关于药物起效时间的机理理解。

新前沿:诊断学与药物发现

PBPK建模的应用甚至超出了治疗学领域。在核医学中,放射性示踪剂被用于器官成像和疾病诊断。这些图像的安全性和质量取决于示踪剂的生物分布和剂量学——即每个器官接受了多少辐射。PBPK模型被用来预测这种分布。通常,初步研究在小鼠等小动物身上进行。一种名为异速标度法(allometric scaling)的技术,将代谢率等生理参数与体重相关联,然后可以在PBPK框架内使用,将小鼠数据转化为对人体器官滞留时间和辐射剂量的预测。这对于开发新的诊断显像剂至关重要。

在这次巡礼中,我们看到了PBPK模型的多种面貌:它是工程师放大实验室数据的工具包,是临床医生实现个性化医疗的水晶球,是毒理学家评估风险的盾牌,也是生物学家窥探复杂系统相互作用的显微镜。它证明了一个简单思想的力量:通过理解各个部分以及连接它们的规则,我们就能开始理解整体。 “虚拟人”不再是科幻小说;它是一个实用的工具,正在重塑我们对待健康、疾病以及构成我们个体本质的方式。