第九章 伤口愈合 (Wound Healing)
学习笔记
9.1 伤口愈合简史 (Brief History of Wound Healing)
Murray的《数学生物学》第九章深入探讨了伤口愈合的数学建模方法。本章从古代医学文献开始，引用了约公元前1300年的Edwin Smith外科纸草文献，这是已知最早的外科医学文献之一，描述了伤口处理的基本原则。
中世纪外科医生Henri de Mondeville（1260-1320年）的工作具有重要影响。他的著作《Chirurgie》于1306年开始撰写，直到1316年仍未完成。该著作对当时的权威Galen学说提出质疑，展现了科学思维在医学中的早期萌芽。Mondeville主张伤口愈合应采用新鲜的肉类敷料，随后使用油脂和蜂蜜每日处理，这一方法在当时具有创新意义。
历史意义：伤口愈合研究是最古老的科学学科之一，其历史可追溯至古代文明。这为后续理解伤口愈合的生物学机制奠定了基础。
9.2 生物学背景：表皮伤口 (Biological Background: Epidermal Wounds)
2.1 伤口愈合的基本阶段
伤口愈合过程传统上分为三个阶段：
- 炎症期（Inflammation）：血凝块形成，白细胞涌入
- 伤口闭合期（Wound Closure）：表皮细胞迁移覆盖伤口
- 细胞外基质重塑期（Extracellular Matrix Remodelling）：瘢痕组织形成
2.2 表皮迁移机制
正常表皮细胞是非运动性的。然而，在伤口附近，细胞经历显著的表型改变，通过指状伪足（lamellipodia）获得运动能力。主要的细胞运动调控因素包括：
- 接触抑制（Contact Inhibition）：细胞运动的主要控制因素
- 趋化作用（Chemotaxis）：通过生长因子谱系进行调控
- 接触导向（Contact Guidance）
- 有丝分裂效应（Mitogenetic Effects）
- 自体抑制剂（Autoinhibitors）
2.3 两种愈合机制的比较
| 特征 | 滚动机制（Rolling） | 滑动机制（Sliding） |
|---|---|---|
| 主导细胞 | 边缘细胞依次植入 | 内部细胞被动响应 |
| 运动方式 | 边缘牵引 | 内部跟随 |
| 证据 | 哺乳动物表皮伤口符合 | 两栖动物表皮伤口 |
2.4 有丝分裂活性
损伤后，表皮有丝分裂活性在伤口边缘附近的1毫米带宽内增加。有丝分裂活性在伤口边缘处最高，可达正常表皮的15倍。刺激有丝分裂增加的因素包括：
- 接触抑制的缺失（适用于有丝分裂和细胞运动）
- 细胞形状改变：细胞扩展变平会增加分裂速率
生长调节因子： - 抑制剂：已有两种确定的表皮抑制剂，在细胞周期的不同点发挥作用 - 激活剂：角化细胞生长因子（KGF）是表皮伤口愈合中的重要调节因子
9.3 表皮伤口愈合模型 (Model for Epidermal Wound Healing)
Sherratt和Murray（1990）提出了两个简单模型，通过与实验比较得出更符合生物学实际的模型。该模型包含两个守恒方程：
3.1 模型方程的一般形式
其中 \(n\) 是单位面积上皮细胞密度，\(c\) 是调节有丝分裂的化学物质浓度。
3.2 化学物质产生项
激活剂情况： $$ f(n) = \lambda c_0 \cdot \frac{n}{n_0} \cdot \frac{n_0^2 + \alpha^2}{n^2 + \alpha^2} \quad \text{对于激活剂} \quad (9.3) $$
抑制剂情况： $$ f(n) = \frac{\lambda c_0}{n_0} \cdot n \quad \text{对于抑制剂} \quad (9.4) $$
3.3 有丝分裂的化学控制项 \(s(c)\)
激活剂情况（包含在大c时s(c)下降至s(0)的实验发现）： $$ s(c) = k \cdot \left[\frac{2c_m(h - \beta)c}{c_m^2 + c^2} + \beta\right] \quad \text{其中} \quad \beta = \frac{c_0^2 + c_m^2 - 2hc_0c_m}{(c_0 - c_m)^2} \quad (9.5) $$
抑制剂情况： $$ s(c) = \frac{(h-1)c + hc_0}{2(h-1)c + c_0} \cdot k \quad (9.6) $$
3.4 最终模型方程
3.5 初始条件和边界条件
初始条件（在伤口区域内）： $$ n = c = 0 \quad \text{在} \quad t = 0 \quad (9.9) $$
边界条件（在伤口边界上）： $$ n = n_0, \quad c = c_0 \quad \text{在伤口边界上，对所有} t \quad (9.10) $$
9.4 无量纲化、线性稳定性和参数值 (Nondimensional Form, Linear Stability and Parameter Values)
4.1 无量纲化
使用长度尺度 \(L\)（伤口的典型线性维度）和时间尺度 \(1/k\)（细胞周期时间），引入无量纲量：
4.2 无量纲化后的模型方程
激活剂情况： $$ g(n) = \frac{n(1 + \alpha^2)}{n^2 + \alpha^2}, \quad s(c) = \frac{2c_m(h - \beta)c}{c_m^2 + c^2} + \beta \quad (9.16) $$
抑制剂情况： $$ g(n) = n, \quad s(c) = \frac{(h-1)c + h}{2(h-1)c + 1} \quad (9.17) $$
4.3 线性稳定性条件
要求未受伤状态稳定（小扰动下），而受伤状态不稳定。通过线性分析得：
4.4 参数估计
| 参数 | 激活剂 | 抑制剂 |
|---|---|---|
| 化学衰减半衰期 | ~12小时（\(\lambda \approx 0.3\) h\(^{-1}\)） | ~12小时（\(\lambda \approx 0.05\) h\(^{-1}\)） |
| 细胞周期时间 | 100小时（\(k = 0.01\) h\(^{-1}\)） | 相同 |
| 细胞扩散系数 \(D\) | \(3.5 \times 10^{-10}\) cm\(^2\) s\(^{-1}\) | \(6.9 \times 10^{-11}\) cm\(^2\) s\(^{-1}\) |
| 化学扩散系数 \(D_c\) | \(3.1 \times 10^{-7}\) cm\(^2\) s\(^{-1}\) | \(5.9 \times 10^{-6}\) cm\(^2\) s\(^{-1}\) |
9.5 表皮伤口修复模型的数值解 (Numerical Solution for the Epidermal Wound Repair Model)
5.1 数值方法
Sherratt和Murray（1990, 1991）在径向对称几何条件下数值求解了模型方程，并将结果与体内外实验数据进行了比较。
5.2 愈合阶段特征
数值解呈现出两个特征阶段：
- 滞后期（Lag Phase）：初始阶段
- 线性期（Linear Phase）：伤口半径随时间线性减小
5.3 愈合速度
对于半径0.5 cm的伤口，量纲化波速为： - 激活剂：\(2.6 \times 10^{-3}\) mm h\(^{-1}\) - 抑制剂：\(1.2 \times 10^{-3}\) mm h\(^{-1}\)
与Van den Brenk（1956）研究中的\(8.6 \times 10^{-3}\) mm h\(^{-1}\)进行比较。
5.4 关键发现
- 激活剂和抑制剂机制与实验数据吻合良好
- 两种机制在愈合时间上存在差异
- 伤口几何形状对愈合时间有影响
9.6 表皮模型的行波解 (Travelling Wave Solutions for the Epidermal Model)
6.1 行波形式
在线性期，解的形式为具有恒定形状和恒定速度的行波。设： $$ n(x,t) = N(z), \quad c(x,t) = C(z), \quad z = x + at $$
代入模型方程得： $$ aN' = DN'' + s(C)N(2-N) - N \quad (9.19) $$ $$ aC' = D_cC'' + \lambda g(N) - \lambda C \quad (9.20) $$
6.2 边界条件
6.3 两种近似
由于ODE系统是四阶的，考虑两个近似来降低系统阶数：
近似1：\(\lambda = \infty\)（化学平衡）
此时方程简化为二阶ODE： $$ aN' = \frac{1}{D}N - \frac{1}{D}\psi(N) \quad \text{其中} \quad \psi(N) = sNg(N)(2-N) - N \quad (9.22) $$
波速下界：\(a_{min} = 2\sqrt{D(2s(0)-1)}\)
对于所用参数： - 激活剂：\(a_{min} = 0.01\)（数值解为0.05） - 抑制剂：\(a_{min} = 0.09\)（数值解为0.03）
近似2：\(D = 0\)（无细胞扩散）
此时系统降为三阶（方程9.24-9.25），通过数值求解验证了与完整PDE系统的良好一致性。
6.4 激活剂的扰动分析
对于激活剂模型，将解展开为： $$ N(z;\epsilon) = N_0(z) + \epsilon N_1(z) + \epsilon^2 N_2(z) + \cdots \quad (9.28) $$
其中 \(N_0 = \xi/(1+\xi)\)，\(\xi = e^{z/a}\)
化学浓度解涉及复杂积分表达式（方程9.29）。
9.7 表皮伤口模型的临床意义 (Clinical Implications of the Epidermal Wound Model)
7.1 有丝分裂调节因子的局部应用
单次添加：无显著效果，化学物质快速扩散和指数衰减
渐进释放（使用浸有化学物质的敷料）：效果显著
修改后的模型在c方程右边添加常数项 \(c_{dress}\)： $$ \frac{\partial c}{\partial t} = D_c\nabla^2 c + \lambda g(n) - \lambda c + c_{dress} \quad \text{（修改后）} $$
抑制剂情况下的释放速率影响更大。
7.2 伤口几何形状的影响
使用参数化伤口形状家族研究愈合时间变化：
矩形伤口：愈合时间随长宽比增加而减少
参数化形状函数（方程9.30）： $$ y = f_{shape}(x;\alpha) = \frac{1}{2}\sqrt{1+\frac{1}{\alpha^2}} - \text{sign}(\alpha)\sqrt{\frac{1}{4}\left(1+\frac{1}{\alpha^2}\right) - \left(x-\frac{1}{2}\right)^2} + \frac{1}{2\alpha} \quad (9.30) $$
当 \(\alpha\) 从-1增加到+1时，伤口形状从尖点形经菱形变为卵形。
7.3 关键预测
- 尖头伤口：在愈合过程中变圆（激活剂和抑制剂均如此）
- 卵形伤口：在愈合过程中变平，抑制剂机制下变平更显著
- 更显著的变平导致更大的伤口周长，从而加速愈合
7.4 简单模型（忽略生化效应）
其中 \(\epsilon\) 是有丝分裂相对迁移作用重要性的无量纲参数。
9.8 胚胎表皮修复机制 (Mechanisms of Epidermal Repair in Embryos)
8.1 胚胎与成体愈合的关键差异
| 特征 | 胚胎愈合 | 成体愈合 |
|---|---|---|
| 愈合速度 | 快（约20小时） | 较慢 |
| 运动机制 | 肌动蛋白绳（Actin Cable） | 伪足爬行 |
| 瘢痕形成 | 无 | 有 |
| 生长因子调控 | 水平较低 | 水平较高 |
8.2 肌动蛋白绳的发现
Martin和Lewis（1991, 1992）的实验研究表明： - 在伤口边缘观察到粗大的肌动蛋白绳 - 绳位于伤口边缘的基底层细胞内 - 伤口后1小时内形成，持续至伤口闭合 - 通过荧光标记的鬼笔环肽（phalloidin）染色观察到
8.3 "钱包绳"机制
伤口边缘的肌动蛋白绳产生环向张力，像拉绳一样将边缘向内拉。这种机制可以解释： - 表皮岛的收缩 - 伤口的快速闭合
8.4 临床意义
胚胎伤口愈合无瘢痕，这一理解对临床伤口管理具有深远意义。胎儿手术是一个高风险领域，已用于治疗脊柱裂（自1998年起），并有望用于面部畸形等缺陷的修复。
9.9 胚胎伤口中的肌动蛋白排列：力学模型 (Actin Alignment in Embryonic Wounds: Mechanical Model)
9.1 模型背景
修改Murray-Oster（1984）的机械化学模型，忽略粘性效应，专注于微丝各向异性的影响。
9.2 平衡应力张量
其中： - \(G(\mathbf{r})\)：细胞内肌动蛋白丝密度 - \(\mathbf{u}\)：位移向量 - \(\boldsymbol{\epsilon}\)：应变张量 - \(E\)、\(\nu\)：弹性参数 - \(\tau_0\)：主动牵引力参数 - \(\beta\)：肌动蛋白丝密度影响牵引力的关键参数
9.3 力平衡方程
边界条件： $$ \mathbf{u}(\infty) = 0, \quad \boldsymbol{\sigma} \cdot \mathbf{n} = 0 \quad \text{在自由边界上} \quad (9.34) $$
9.4 肌动蛋白密度关系
假设肌动蛋白总量守恒： $$ G(\mathbf{r})(1 + \nabla \cdot \mathbf{u}) = \kappa \quad (9.35) $$
9.5 一维模型分析
对于一维slash伤口，方程简化为： $$ \frac{d}{dx}\left[\frac{1}{1+u_x}(E+\nu)u_x + \frac{1}{1+\beta u_x}\right] - \frac{\lambda u}{1+u_x} = 0 \quad (9.40) $$
边界条件： $$ (E+\nu)u_x + \frac{1}{1+\beta u_x} = 0 \quad \text{在} x=0, \quad u=0 \quad \text{在} x=\infty \quad (9.41) $$
9.6 参数约束条件
存在唯一单调递减解的条件（物理约束 \(u_x > -1\)）： $$ \beta < 1 - \frac{1}{E+\nu}, \quad 1 \leq E+\nu \leq 2 \quad (9.47) $$
或 $$ \beta \leq \frac{1}{4}(E+\nu), \quad E+\nu > 2 $$
9.7 关键发现
- 参数必须接近临界值才能捕捉实验观察到的行为
- 对于 \(\beta > 1/2\)，无法同时捕捉收缩和肌动蛋白聚集
- 肌动蛋白绳形成于伤口边缘几个细胞长度范围内
- 回缩约为14%，与实验观察一致
9.8 径向对称解
对于圆形伤口，在柱坐标系下分析。应力张量的主值为： $$ \sigma_{rr} = \frac{Eu + (u+u/r) + [1+\beta(u+u/r)]^{-1}}{1+u+u/r} \quad (9.49) $$
9.9 参数的生物学解释
- \(E\)、\(\nu\)：与表皮 sheet 的弹性特性直接相关
- \(\lambda\)：表皮与基底层的附着强度
- \(\beta\)：肌动蛋白丝密度影响牵引力的程度，通过肌球蛋白横桥形成的协同效应体现
9.10 二维应力排列的肌动蛋白丝机械模型 (Mechanical Model with Stress Alignment of the Actin Filaments in Two Dimensions)
10.1 改进的收缩模型
采用更精确的收缩形式： $$ \vartheta = -(p_1 + p_2 + p_1p_2) = -(u_r + u/r + uu_r/r) \quad (9.54) $$
10.2 主应力分量
10.3 肌动蛋白丝密度函数
其中 \(F\) 满足： 1. \(F(\phi; 0) = \delta(\phi)\)：单向应力时，所有丝沿应力方向排列 2. \(F(\phi; 1)\) 常数：各向同性应力时，微丝网络随机取向 3. 归一化条件 4. 对称性条件
近似形式： $$ \int_0^{\pi/2} F(\phi; \sigma_1/\sigma_2)\cos\phi \, d\phi \approx \frac{\sqrt{\pi\sigma_1}}{2\sigma_1 + (\pi-2)\sigma_2} \quad (9.61) $$
10.4 最终方程
得到高度非线性的ODE（方程9.64），必须数值求解。
10.5 模型预测与实验比较
- 肌动蛋白绳位于伤口边缘
- 预测的细胞回缩约为30-40 μm，与实验观察吻合
- 应力诱导的微丝网络优先排列得以验证
10.6 肌动蛋白绳形成的意义
模型表明，胚胎表皮伤口中肌动蛋白绳的初始形成可能是表皮细胞片层创后机械平衡的副产品。应力诱导的肌动蛋白丝排列是绳索形成的核心机制。
公式汇总表 (Formula Summary Table)
| 编号 | 公式名称 | 公式 | 用途 |
|---|---|---|---|
| (9.1) | 细胞密度变化 | \(\frac{\partial n}{\partial t} = \text{迁移} + \text{生成} - \text{损失}\) | 表皮细胞守恒 |
| (9.2) | 化学浓度变化 | \(\frac{\partial c}{\partial t} = \text{扩散} + \text{产生} - \text{衰减}\) | 有丝分裂调节因子守恒 |
| (9.3) | 激活剂产生项 | \(f(n) = \lambda c_0 \cdot \frac{n(1+\alpha^2)}{n^2+\alpha^2}\) | 激活情况 |
| (9.4) | 抑制剂产生项 | \(f(n) = \frac{\lambda c_0}{n_0} \cdot n\) | 抑制情况 |
| (9.5) | 激活剂 \(s(c)\) | \(s(c) = k[\frac{2c_m(h-\beta)c}{c_m^2+c^2}+\beta]\) | 有丝分裂化学控制 |
| (9.6) | 抑制剂 \(s(c)\) | \(s(c) = \frac{(h-1)c+hc_0}{2(h-1)c+c_0} \cdot k\) | 有丝分裂化学控制 |
| (9.7) | n的模型方程 | \(\frac{\partial n}{\partial t} = D\nabla^2 n + s(c)n(2-\frac{n}{n_0})-kn\) | 完整表皮模型 |
| (9.8) | c的模型方程 | \(\frac{\partial c}{\partial t} = D_c\nabla^2 c + f(n)-\lambda c\) | 完整化学模型 |
| (9.12) | 无量纲n方程 | \(\frac{\partial n}{\partial t} = D\nabla^2 n + s(c)n(2-n)-n\) | 无量纲化模型 |
| (9.13) | 无量纲c方程 | \(\frac{\partial c}{\partial t} = D_c\nabla^2 c + \lambda g(n)-\lambda c\) | 无量纲化模型 |
| (9.18) | 稳定性条件 | \(s(0) > 1/2\) 给出参数约束 | 线性稳定性要求 |
| (9.19) | 行波-N方程 | \(aN' = DN'' + s(C)N(2-N)-N\) | 行波分析 |
| (9.20) | 行波-C方程 | \(aC' = D_cC'' + \lambda g(N)-\lambda C\) | 行波分析 |
| (9.32) | 应力张量 | \(\boldsymbol{\sigma} = G[E\boldsymbol{\epsilon}+(\nabla\cdot\mathbf{u})\mathbf{I}]+\frac{\tau_0}{1+\beta\nabla\cdot\mathbf{u}}\mathbf{I}\) | 力学模型 |
| (9.33) | 力平衡方程 | \(\nabla\cdot\boldsymbol{\sigma}-\lambda G\mathbf{u}=0\) | 胚胎伤口平衡 |
| (9.35) | 肌动蛋白守恒 | \(G(\mathbf{r})(1+\nabla\cdot\mathbf{u})=\kappa\) | 肌动蛋白密度关系 |
| (9.40) | 一维力平衡 | \(\frac{d}{dx}[\frac{(E+\nu)u_x}{1+u_x}+\frac{1}{1+\beta u_x}]-\frac{\lambda u}{1+u_x}=0\) | 一维模型 |
| (9.47) | 参数约束 | \(\beta < 1-\frac{1}{E+\nu}\) 或 \(\beta \leq \frac{1}{4}(E+\nu)\) | 物理解存在条件 |
| (9.54) | 收缩精确形式 | \(\vartheta = -(u_r+u/r+uu_r/r)\) | 二维改进模型 |
| (9.59) | 主应力分量 | \(\sigma_1=G_1[E_1+\nabla\cdot\mathbf{u}]+G_1\tau\) | 应力诱导排列 |
| (9.60) | 肌动蛋白密度 | \(G_1=G_0\int_0^{\pi/2}F(\phi;\sigma_1/\sigma_2)\cos\phi\,d\phi\) | 各向异性密度 |
核心概念总结
成人与胚胎伤口愈合对比
| 方面 | 成人表皮 | 胚胎表皮 |
|---|---|---|
| 运动机制 | 伪足爬行 | 肌动蛋白绳"钱包绳" |
| 调控方式 | 生化生长因子 | 力学响应 |
| 愈合速度 | 较慢 | 快（约20小时） |
| 瘢痕形成 | 有 | 无 |
| 数学模型 | 反应扩散模型 | 力学平衡模型 |
主要建模贡献者
- Sherratt, J.A.
- Murray, J.D.
- Maini, P.K.
- Cook, J.
- Tracqui, P.
- Martin, P.
- Lewis, J.
参考文献要点
- Sherratt & Murray (1990, 1991, 1992): 表皮伤口愈合的反应扩散模型
- Martin & Lewis (1991, 1992): 胚胎伤口的肌动蛋白绳机制
- Sherratt (1991, 1993): 胚胎伤口愈合的力学模型
- Oster & Odell (1984): 细胞牵引力的力学基础
笔记完成日期：2026年5月10日
章节总页数：约100页（原文）
核心公式数量：25+个主要方程
学习要点： 1. 掌握表皮伤口愈合的反应扩散建模方法 2. 理解行波解的物理意义和数学推导 3. 认识胚胎与成体愈合的本质差异 4. 学会应用力学模型解释肌动蛋白绳的形成