透皮吸收量计算模型 - 中析研究所生物检测中心

透皮吸收量计算模型：原理、方法与考量

透皮给药系统因其非侵入性、避免首过效应、提供可控释放等优势，在药物递送领域占据重要地位。准确预测药物透过皮肤屏障进入体循环的量（透皮吸收量）是系统设计与评价的核心。透皮吸收量计算模型为这一过程提供了重要的理论指导和定量工具。

一、透皮吸收的生理学与物理学基础

皮肤是人体最大的器官，也是药物透皮吸收的主要屏障，其结构直接影响药物的渗透路径与速率：

皮肤结构屏障：
- 角质层 (Stratum Corneum, SC)： 最外层，由致密排列的死细胞（角质细胞）和细胞间脂质（主要为神经酰胺、胆固醇、游离脂肪酸）构成。这是药物透皮吸收最主要的限速屏障，呈亲脂性。
- 活性表皮 (Viable Epidermis)： 位于角质层之下，包含活细胞，水性环境，药物在此扩散阻力相对较小。
- 真皮 (Dermis)： 主要由结缔组织构成，富含血管和淋巴管。药物扩散至此通常即被认为被吸收进入体循环。
- 皮下组织 (Hypodermis)： 主要为脂肪组织，对大多数药物的吸收贡献相对较小，但可作为脂溶性药物的储库。
药物渗透途径：
- 细胞间途径 (Intercellular route)： 药物分子穿过角质层细胞间的脂质基质。这是大多数药物，尤其是小分子亲脂性药物最主要的渗透路径。
- 跨细胞途径 (Transcellular route)： 药物分子连续穿过角质细胞和细胞间脂质。对亲水性药物可能更重要，但阻力较大。
- 附属器途径 (Appendageal route)： 通过毛囊、汗腺和皮脂腺等皮肤附属器。这些通道在皮肤表面所占比例很小（通常<0.1%），但在渗透初始阶段或特定化合物（如大分子、离子）的渗透中可能起一定作用，但稳态下贡献通常较小且常被忽略。
渗透驱动力： 药物分子主要通过被动扩散透过皮肤，其驱动力是药物在皮肤给药部位的浓度与其在皮肤深层（如真皮层毛细血管）的浓度差（浓度梯度）。

二、核心数学模型：菲克扩散定律

药物经皮肤的被动扩散过程通常基于菲克第一扩散定律 (Fick's First Law of Diffusion) 进行建模。在假设皮肤为均质膜、扩散为稳态的条件下，定律表述为：

$$ J = -D \frac{dC}{dx} $$

其中：

$$ J $$\)： 通量 (Flux)，单位时间单位面积内通过的药物量（如 μg/cm²/h）。这是计算透皮吸收量的核心参数。
$$ D $$\)： 扩散系数 (Diffusion Coefficient)，反映药物分子在皮肤屏障（主要是角质层）中迁移能力的参数（如 cm²/h）。
$$ \frac{dC}{dx} $$\)： 浓度梯度 (Concentration Gradient)，药物浓度 $$ C $$\) 沿渗透路径 $$ x $$\) 方向的变化率（如 μg/cm³/cm）。

对于透皮吸收，在稳态 (Steady State) 条件下（即扩散速率恒定，浓度梯度不随时间变化），并假设皮肤屏障厚度为 $$ h $$\)（主要是角质层有效厚度），药物在皮肤外表面（供体侧）浓度为 $$ C_v $$\)，在皮肤内表面（受体侧）浓度为 $$ C_r $$\)（通常假设受体侧浓度远低于供体侧，$$ C_r \approx 0 $$\)），则菲克定律可简化为：

$$ J_{ss} = \frac{D}{h} \times C_v = K_p \times C_v $$

其中：

$$ J_{ss} $$\)： 稳态透皮通量 (Steady-State Flux)。
$$ \frac{D}{h} $$\)： 渗透系数 (Permeability Coefficient, $$ K_p $$\))（如 cm/h）。$$ K_p $$\) 是描述特定药物-皮肤组合渗透特性的关键参数，综合了药物在皮肤屏障中的扩散能力和分配特性。

透皮吸收量的计算通常是基于稳态通量 $$ J_{ss} $$\)：
$$ Amount_{absorbed} = J_{ss} \times A \times t $$
其中：

$$ A $$\)： 给药面积 (Area of Application)（cm²）。
$$ t $$\)： 时间 (Time)（h）。该公式适用于稳态持续时间内的吸收量计算。

三、关键模型参数及其获取

模型的准确性高度依赖于参数的可靠性：

渗透系数 ($$ K_p $$\)):
- 定义： 单位浓度梯度下的稳态通量 ($$ K_p = J_{ss} / C_v $$\)）。
- 获取方法： 主要通过体外扩散实验测得。
  - Franz扩散池 (Franz Diffusion Cell)： 最标准装置。将离体皮肤（人或动物）夹在供给室（含药物溶液或制剂）和接收室（含接收介质）之间，维持恒温（通常32℃）。定时采集接收液测定药物浓度，计算累积透过量 $$ Q $$\) 对时间 $$ t $$\) 作图。
  - 滞后时间 ($$ Lag Time, t_L $$\))： 累积透过量曲线线性部分外推至时间轴所得截距，反映药物在皮肤中建立稳态扩散所需时间 ($$ t_L \approx \frac{h^2}{6D} $$\))。
  - 稳态通量 ($$ J_{ss} $$\))： 累积透过量-时间曲线的线性区斜率 ($$ J_{ss} = \frac{dQ}{dt} $$\))。
  - $$ K_p $$\) 计算： $$ K_p = \frac{J_{ss}}{C_v} $$\)（$$ C_v $$\) 为供给室药物初始浓度）。
- 估算方法 (Potts-Guy等)： 当实验数据缺乏时，可根据药物的理化性质（分子量 $$ MW $$\)、辛醇/水分配系数 $$ \log P $$\)、熔点 $$ MP $$\) 等）通过经验方程估算 $$ K_p $$\)：
  $$ \log K_p = a \cdot \log P + b \cdot \log MW + c \cdot MP + d $$
  (系数a, b, c, d 由大量实验数据拟合得到)。这类模型精度有限，但对早期筛选有参考价值。
扩散系数 ($$ D $$\)) 与屏障厚度 ($$ h $$\)):
- 通常不是直接输入参数，而是整合在 $$ K_p $$\) ($$ K_p = \frac{D \cdot K_m}{h} $$\), $$ K_m $$\)为药物在皮肤/供给相间的分配系数）或通过滞后时间 $$ t_L $$\) 间接反映 ($$ D \approx \frac{h^2}{6 t_L} $$\))。角质层有效厚度 $$ h $$\) 常取经验值（如10-20μm）。
给药浓度 ($$ C_v $$\)):
- 对于溶液型给药系统，$$ C_v $$\) 通常取药物在溶剂中的饱和溶解度 ($$ C_s $$\))（假设为饱和溶液），因为最大通量 $$ J_{max} = K_p \times C_s $$\)。
- 对于复杂的半固体或贴剂制剂，$$ C_v $$\) 指药物在制剂-皮肤界面的有效浓度，可能与制剂中药物的热力学活度 ($$ \gamma \cdot C $$\)) 相关，估算更为复杂。
时间 ($$ t $$\)):
- 稳态模型假设给药时间远大于滞后时间 ($$ t >> t_L $$\))。
- 对于非稳态（给药初期或短时间接触）吸收计算，需要解菲克第二定律的非稳态扩散方程（通常需要数值解法）。

四、超越简单模型：复杂性与改进

基础稳态模型虽然应用广泛，但有其固有局限。更精细的模型考虑了更多复杂性：

非稳态模型：
- 菲克第二定律： $$ \frac{\partial C}{\partial t} = D \frac{\partial^2 C}{\partial x^2} $$\)
- 用于描述药物在皮肤内浓度分布随时间变化的过程，计算任意时刻的吸收速率和累积吸收量。需结合初始条件和边界条件求解（常需数值方法如有限差分法、有限元法）。
多层皮肤模型：
- 将皮肤视为由角质层、活性表皮、真皮等组成的多层结构，各层具有不同的厚度 $$ h_i $$\)、扩散系数 $$ D_i $$\)、分配系数 $$ K_i $$\)。
- 药物在层间界面需满足浓度连续 ($$ C_{i} = K_{i/i+1} C_{i+1} $$\)) 和通量连续 ($$ J_i = J_{i+1} $$\)) 条件。
- 更贴近皮肤生理结构，能更精确模拟药物渗透动力学，但参数获取困难。
角质层脂质路径模型：
- 将角质层视为由亲水性角质细胞（“砖”）嵌入亲脂性细胞间脂质基质（“灰浆”）的复合介质。
- 药物扩散主要发生在曲折的脂质路径中，模型需考虑路径的曲折度 ($$ \tau $$\)) 和脂质相的体积分数 ($$ \epsilon $$\))。
- 通量表达式更复杂：$$ J \propto \frac{\epsilon}{\tau} D_{lipid} \frac{\Delta C}{h} $$\)。有助于理解渗透促进剂作用机制。
结合清除动力学：
- 真皮层的药物会被毛细血管和淋巴管清除进入体循环。精细模型会将真皮层的扩散与一级清除动力学相结合：
  $$ \frac{\partial C_d}{\partial t} = D_d \frac{\partial^2 C_d}{\partial x^2} - k_{clear} C_d $$
  ($$ C_d $$\)为真皮层药物浓度，$$ D_d $$\)为真皮层扩散系数，$$ k_{clear} $$\)为清除速率常数)。
- 适用于需要精确模拟系统暴露量（如药代动力学预测）的情况。
数值模拟与计算机模型：
- 利用计算机软件（如COMSOL Multiphysics, MATLAB PDE Toolbox, 专用药动学/药效学软件）求解描述复杂几何结构、异质介质、时变边界条件等的偏微分方程组。
- 可整合药物的理化性质、皮肤生理参数、制剂特性、甚至血流动力学，构建高度个性化的虚拟渗透模型，用于制剂优化和体内吸收预测。

五、模型应用与局限性

应用：

制剂设计与优化： 预测不同处方（渗透促进剂、基质、药物浓度等）对透皮速率的影响。
体外-体内相关性 (IVIVC)： 建立体外释放/渗透数据与体内吸收结果的关联，减少临床前研究。
生物等效性评价支持： 理解关键质量属性（CQAs）对渗透的影响。
暴露评估与安全性： 评估外源性化学物质（化妆品、环境污染物）透过皮肤的可能吸收量及潜在风险。
机理研究： 探索药物渗透途径、屏障功能、渗透促进剂作用机制。

局限性与挑战：

皮肤异质性： 个体差异（年龄、性别、种族、解剖部位、健康状况）、离体皮肤的处理与储存均显著影响渗透性参数。体外实验常用动物皮肤（如猪耳、鼠皮）或有限来源的人体皮肤，与实际人体应用存在差异。
参数获取困难： 准确的 $$ D $$\)、$$ h $$\)、$$ K_m $$\)、多层/脂质路径参数、清除速率等难以精确测定，常依赖估算或假设。
模型简化： 基础模型忽略了许多复杂因素（如皮肤代谢、渗透促进剂的动力学作用、制剂在皮肤表面的变化、皮肤屏障的动态变化如水合作用）。
非稳态复杂性： 短时间接触或脉冲式给药预测需要求解复杂方程，不易获得解析解。
体内环境模拟： 体外实验难以完全模拟体内血流、皮肤温度梯度、机械摩擦等动态生理条件。
制剂复杂性： 半固体、贴剂等制剂中药物的释放和皮肤渗透相互耦合，建模难度大。

六、结论

透皮吸收量计算模型是基于物理学定律构建的定量工具，从简单的稳态菲克扩散模型到复杂的多层数值模型，为理解、预测和优化药物及化学物质经皮渗透过程提供了强大的框架。稳态模型因其简洁实用，在制剂筛选和初步评估中广泛应用。然而，皮肤屏障的复杂结构和生理、个体差异性以及制剂-皮肤相互作用的动态性，使得精确预测体内透皮吸收量仍面临挑战。未来模型的发展方向在于更深入地整合皮肤生理学、药剂学、计算科学和实验数据，提高个性化预测能力，并更好地表征非稳态和复杂制剂行为，以服务于更安全、有效的透皮药物递送系统和更精准的化学物质暴露风险评估。