跳转至

第 6 章 细胞 Potts 模型在生物医学中的应用（The Cellular Potts Model in Biomedicine）

作者

Nicholas J. Savill（Edinburgh, 现 Centre for Inflammation Research）、Roeland M. H. Merks（Netherlands Institute for Systems Biology, Centrum Wiskunde & Informatica, Leiden）。本章是 CPM/GGH 在生物医学三大领域（表皮干细胞、肿瘤、血管生成）的方法论应用综述，体现了"CPM 作为机制发现工具"的范式。

内容概述

本章是 CPM/GGH 的"应用章"。它从三个具体生物医学问题出发，展示 CPM 如何提供"自下而上"的机制理解：

  1. 表皮干细胞集群的稳态调控（pp. 137–139）——Savill & Sherratt 用 CPM 模拟角质细胞集群的"自限大小"机制，发现只有"周边干细胞分化"+"干细胞内聚 > 干细胞对其他类型粘附"才能形成稳定集群。
  2. 肿瘤生长与侵袭（pp. 139–141）——Stott et al. 用 3D CPM 模拟良性无血管肿瘤；Jiang et al. 加入亚细胞布尔网络模拟多细胞肿瘤球体；Turner & Sherratt 研究 haptotaxis + ECM 降解的侵袭机制；Knewitz & Mombach 提出"恶性细胞-健康细胞粘附 > 恶性细胞内聚"是侵袭的充分条件。
  3. 血管生成与血管形成（pp. 141–149）——Merks, Glazier 等系统化地"反向工程"出四种可能机制，用 CPM 筛选出最可能解释体外 HUVEC 实验的机制——"主动伸长 + 接触抑制趋化"是核心。

读者前置知识：Ch 4–5 的 CPM/GGH 基础、肿瘤生物学（ECM、MMP、haptotaxis）、血管生物学（VEGF、VE-cadherin、angiogenesis vs. vasculogenesis）。

核心方程与概念

1. 表皮干细胞集群模型（Savill-Sherratt）

CPM 中三种角质细胞类型：干细胞（stem）、过渡扩增细胞（transit-amplifying）、有丝分裂后细胞（postmitotic）。分化规则： - 干细胞只在集群周边（接触非干细胞）才分化为过渡扩增细胞。 - 过渡扩增细胞分裂 3–5 代后变为 postmitotic。 - Postmitotic 细胞迁移离开基底层。

Delta-Notch 信号作为"分化诱导"——高 Delta 表达的细胞"使邻居分化、自己对分化信号不敏感"。

关键结论：只有"周边干细胞分化"+"干细胞内聚 > 对其他类型粘附"才形成稳定集群；集群大小由"周边面积 / 总体积"自动调控（本质是几何调控）。

2. 良性无血管肿瘤的 3D CPM（Stott et al.）

  • 营养（葡萄糖、氧）扩散方程 \(\partial_t n = D\nabla^2 n - c\)。
  • 细胞生长率 \(r = r_{\max} \cdot n/(n + K)\)（Michaelis-Menten）。
  • 3D 模拟显示球形结构：外层增殖 → 中层静止 → 内层坏死。

3. 多细胞肿瘤球体的多尺度模型（Jiang et al. 2005）

三层耦合： - 亚细胞层：Boolean network 模拟细胞周期 arrest 蛋白调控。 - 细胞层：CPM 模拟增殖、粘附、活力。 - 细胞外层：反应-扩散方程模拟营养、代谢物、生长促进/抑制因子。

拟合实验得到：生长促进因子 \(D \sim 10^{-6}\) cm²/h、抑制因子 \(D \sim 10^{-7}\) cm²/h。

4. 侵袭模型（Turner-Sherratt）

CPM + 蛋白酶 + ECM 方程： - 蛋白酶 \(m\) 降解 ECM，ECM \(f\) 形成梯度。 - 恶性细胞 haptotaxis 沿 ECM 梯度上移：\(\Delta H\) 中加入 \(-\mu_{\text{hap}} (f_{\text{neighbor}} - f_{\text{site}})\)。

关键观察：当 haptotaxis 强度 > 内聚强度时，侵袭前缘会从主体断裂；增殖反过来说抑制侵袭（更厚的"锚"将前缘与主体连接）。TGF-β 的双面性：低 haptotaxis 敏感度→良性；高敏感度→肿瘤更小但更侵袭性。

5. 血管生成的 Gamba-Serini 连续模型（背景）

Gamba, Serini 等从天体物理学的 Burgers 方程出发，把内皮细胞视为"分泌化学吸引物、沿其梯度迁移"： $$ \frac{\partial n}{\partial t} + \nabla \cdot (n \mathbf{v}) = 0 $$ $$ n \frac{\partial \mathbf{v}}{\partial t} + n (\mathbf{v} \cdot \nabla) \mathbf{v} = \nu \nabla^2 \mathbf{v} - \nabla p $$ $$ \partial_t c = D\nabla^2 c + a n - \gamma c $$ 问题：方程含有 \(n (\mathbf{v} \cdot \nabla) \mathbf{v}\) 项，意味着细胞加速度正比于梯度——但真实细胞没有惯性（黏度主导，Aristotelian regime）。Merks 等修改为 \(\mathbf{v} \propto \nabla c\)。

6. 血管形成的 CPM 替代（Merks-Glazier）

修改 Savill-Hogeweg 的趋化公式（Ch 1, Eq. 30）： - 化学吸引物方程：\(\partial_t c = 0.1 \cdot \text{boz}_{c,0}(1 - \text{boz}_{c,0})c + D\nabla^2 c\)（只在细胞内部释放，ECM 中衰变）。 - \(\mu = 500\) 的趋化强度驱动细胞沿 \(c\) 梯度迁移。 - 关键实验（Fig. 3）：仅靠趋化 → 圆形 cluster（Keller-Segel 集聚）；加主动伸长（target length = 50 sites）→ 网络形成；短距离梯度（慢扩散）→ 被动伸长 → 网络；强粘附→ 被动伸长 → 网络；接触抑制趋化（VE-cadherin 介导）→ 与形变无关的网络。

7. 接触抑制趋化（contact inhibition of chemotaxis）

VE-cadherin 在细胞-细胞界面聚集，抑制该处 filopodia 形成——即"细胞只能向自由表面伸出突起"。在 CPM 中实现为： - \(\mu = 0\) 在细胞-细胞界面；\(\mu = 500\) 在细胞-基底界面。 - 即使细胞被强制为圆形，也能形成血管网络。 - 这一机制独立于细胞形变——是"细胞-通讯"机制的直接体现。

8. 萌发（sprouting）不稳定性

在 cluster 初始条件下，化学吸引物形成 Gaussian 状剖面，中心浓度最高。无接触抑制时所有细胞都向中心移动，挤压内部细胞；有接触抑制时内部细胞对趋化信号"不敏感"，被外层细胞推开——结果是 cluster 表面起皱（buckling instability），形成"萌发"结构。

关键结论

  • 结论 1：表皮干细胞集群的稳态大小由几何调控（周边面积 / 总体积）自动控制，而非传统的"分裂率调控"——这一发现挑战了 20 世纪的细胞动力学共识。
  • 结论 2：3D CPM 能准确复现良性无血管肿瘤的"球形分层结构"（增殖外层 + 静止中层 + 坏死核心），多尺度模型可以定量匹配体外肿瘤球体的生长曲线。
  • 结论 3：肿瘤侵袭中细胞-ECM 粘附比细胞-细胞粘附更重要；蛋白酶分泌 + haptotaxis 的协同是侵袭的关键。TGF-β 的"低剂量抑癌、高剂量促癌"双面性是 GGH 模型给出的可证伪预测。
  • 结论 4：血管生成可以由四种机制分别解释——主动伸长、短距离梯度、强粘附、接触抑制趋化。只有"主动伸长 + 接触抑制趋化"能够同时复现 HUVEC 网络形成 + 模式粗化动力学——这一方法论"反向工程"为生物实验提供了清晰的优先顺序。
  • 结论 5：Gamba-Serini 连续模型中的"细胞惯性"项是生物学不合理的（细胞处于 Aristotelian regime）；Merks-Glazier 的修改 \(\mathbf{v} \propto \nabla c\) 是必要的、且与"细胞最大速度受限"的实验一致。
  • 结论 6：接触抑制趋化机制（VE-cadherin 介导）可以独立于细胞形变解释血管网络——这一发现支持 VE-cadherin knockout 小鼠表现出的"血管岛"表型。
  • 结论 7：血管萌发（sprouting angiogenesis）和血管形成（vasculogenesis）在 CPM 模拟中是同一细胞机制（接触抑制趋化）的两种表现形式——统一了 1980 年代以来被认为不同的两个生物学过程。

挑战和开放性问题

  • 挑战 1：CPM 与 ODE/PDE 的接口稳定性。Merks 的血管模型在每个 MCS 演化 PDE 30 s——这种"sub-cycling"是 ad hoc 的，缺乏严格的稳定性分析。收敛性证明是一个开放的理论问题。
  • 挑战 2：3D 中的生物医学应用。本章的肿瘤和血管模拟都是 2D 简化。3D 模拟在 2007 年还非常昂贵，今天 GPU 可以做到，但3D 模拟的"可解释性"远低于 2D——可视化与定量分析的接口尚未标准化。
  • 挑战 3：与组学数据对接。Savill & Sherratt 的干细胞模型只用了 Delta-Notch 一种信号通路；今天单细胞 RNA-seq 可以同时给出 100+ 信号通路的活性。CPM 如何整合这种"高通量"数据是当前研究热点。
  • 挑战 4：与药代动力学整合。Turner-Sherratt 的 TGF-β 模型是"信号-表型"耦合，但没有整合 PK——药效建模是临床应用的关键环节。
  • 挑战 5：肿瘤异质性。CPM 把"细胞类型"作为离散标签——但真实肿瘤是连续演化的（基因组不稳定、表型可塑性）。今天"单细胞分辨率的肿瘤演化"建模需要新的框架。
  • 挑战 6：血管生成的物理化学细节。VEGF-A 的浓度梯度、ECM 的硬度梯度、流体剪切力——这些都在 Merks 模型中简化。与体外 HUVEC 培养的精确匹配需要更精细的多物理场耦合。
  • 挑战 7：CPM 的"可证伪性"。本章展示的多个机制都能复现血管网络——但哪个是"真正"的机制？这要求实验上对单一基因/分子进行扰动并观察模拟预测。这种"扰动-预测-验证"循环是 2007 年后该领域的方法论进步。

个人反思与批判性分析

本章是 CPM/GGH 应用的"成熟范式"——它不是方法论创新，而是"把已有方法用于重要生物医学问题"的展示。这种"应用驱动"的工作对于把 CPM 从"学术圈玩具"推广到"临床相关工具"至关重要。

几点批判性观察：

  1. "反向工程"是 CPM 在血管研究中的核心方法论贡献。作者不只展示"模型能复现现象"，而是枚举 4 种可能机制 → 用同一 CPM 框架分别测试 → 排除不能复现的 → 留下"最小充分机制"。这是机制发现而非"现象演示"——是 CPM 方法论的成熟表现。

  2. VE-cadherin 介导的接触抑制趋化是 2000 年代的重要发现。它的发现过程体现了 CPM 作为"假设生成器"的角色——而非"假设验证器"。但这一发现至今未在临床上转化为靶向药物——这反映了"理论发现 → 临床应用"的漫长 gap。

  3. "percolation morphogenesis"（渗流形态发生）这一术语。Gamba-Serini 用天体物理学的 Burgers 方程类比血管形成——这种跨学科类比是 21 世纪初生物物理的"风格特征"。今天机器学习 + 单细胞测序的组合正在取代这种"风格"——但"用其他学科的方程理解生命"的精神应当保留。

  4. "haptotaxis + ECM degradation"是肿瘤侵袭的最小机制。本章的 Turner-Sherratt 模型显示这两个机制的协同即可产生指状侵袭——这比连续 PDE 模型（如 Ch 1 Anderson 的 HDC）给出更"机制"性的解释。但 HDC 的优势是"可以耦合到氧气/营养场"——这两种方法应当互补。

  5. 与单细胞实验的定量对照。本章多个模拟给出了"对实验的预测"（如 6 种肿瘤系的 \(D\)–\(\beta\) 线性相关），但这种预测是"事后拟合"还是"事前预测"在文献中常被混淆。真正的事前预测需要模型在数据公开之前完成——这是该领域"可重复性危机"的一个侧面。

  6. 临床相关性。本章的肿瘤模型基于"haptotaxis + ECM 降解"——但 2010 年代的临床肿瘤学发现，免疫微环境（T 细胞、巨噬细胞、PD-L1）才是预后的关键因素。CPM 在 2007 年几乎没有整合免疫——这是该章的一个时代局限。今天"免疫肿瘤学的 CPM"是一个新兴方向。

  7. 本书的方法谱地位。本章是 CPM/GGH 在 Part II 中的"应用章"，与 Ch 1（HDC 肿瘤）、Ch 2（LGCA 粘菌/肿瘤）、Ch 3（CPM-Fokker-Planck 趋化）、Ch 5（CPM 生物物理）、Ch 7（GGH 未来）形成完整的方法论对话。本书的最大价值正是让这些方法在同一框架下被并列讨论——这种"横向对比"在分散的期刊文章中很难出现。

  8. CPM vs ABM（agent-based modeling）的相对优势。今天 NetLogo、Repast、Mesa 等 ABM 框架可以做 CPM 的所有事（甚至更简单），且接口更友好。CPM 的相对优势在于：(a) 严格的统计力学血统（Ch 4）；(b) 显式的"形变-力"映射（Ch 5）；(c) 已经有 30 年的生物医学应用积累。但对新的研究者来说，ABM 可能是更实用的入门选择。

  9. "4 种机制都能产生血管网络"这一发现的解释。这反映了复杂生物系统中的"简并性"（degeneracy）——多个不同的机制可以产生相同的表型。这是生物学而非物理学的特征——物理学中"唯一性定理"是常态。CPM/GGH 揭示了生物系统的这一"软"特征。

  10. "可证伪性"与"可证真性"的张力。本章的工作流是"枚举假设 → 模拟 → 留下能解释现象的"——这在科学哲学上更接近"证真"（verification）而非"证伪"（falsification）。Popper 框架下，CPM 的预测需要"在新的实验条件下失败"才能算严格证伪。这是该领域方法论哲学的开放问题。

重要参考文献

[X1] Savill N.J., Sherratt J.A. Control of epidermal stem cell clusters by Notch-mediated lateral induction. Dev. Biol. 2003;258(1):141-153. DOI: 10.1016/S0012-1606(03)00128-2

[X2] Lowell S., Jones P., Le Roux I., Dunne J., Watt F.M. Stimulation of human epidermal differentiation by Delta-Notch signalling at the boundaries of stem-cell clusters. Curr. Biol. 2000;10(9):491-500. DOI: 10.1016/S0960-9822(00)00451-6

[X3] Stott E.L., Britton N.F., Glazier J.A., Zajac M. Stochastic simulation of benign avascular tumour growth using the Potts model. Math. Comput. Model. 1999;30(5-6):183-198. DOI: 10.1016/S0895-7177(99)00156-8

[X4] Jiang Y., Pjesivac-Grbovic J., Cantrell C., Freyer J.P. A multiscale model for avascular tumor growth. Biophys. J. 2005;89(6):3884-3894. DOI: 10.1529/biophysj.105.060129

[X5] Turner S., Sherratt J.A. Intercellular adhesion and tumour invasion: a discrete model. IMA J. Math. Appl. Med. Biol. 2002;19(3):175-204. DOI: 10.1093/imammb/19.3.175

[X6] Turner S., Sherratt J.A., Cameron D., Murray J.D. Tamoxifen stimulates the invasion of endometriotic cells in a three-dimensional co-culture. J. Theoret. Med. 2001;3(3):179-193. DOI: 10.1080/10273660108833074

[X7] Knewitz M.A., Mombach J.C.M. Computer simulation of the influence of cellular adhesion on the morphology of the interface between tissues of proliferating and quiescent cells. Comp. Biol. Med. 2006;36(1):59-69. DOI: 10.1016/j.compbiomed.2004.08.004

[X8] Merks R.M.H., Glazier J.A. A cell-centered approach to developmental biology. Physica A 2005;352(1):113-130. DOI: 10.1016/j.physa.2004.12.028

[X9] Merks R.M.H., Brodsky S.V., Goligorksy M.S., Newman S.A., Glazier J.A. Cell elongation is a key to in silico blood vessel growth. Dev. Biol. 2006;290(2):371-383. DOI: 10.1016/j.ydbio.2005.11.042

[X10] Merks R.M.H., Perryn E.D., Shirinifard A., Glazier J.A. Contact-inhibited chemotaxis in de novo and sprouting blood-vessel growth. PLoS Comput. Biol. 2008;4(9):e1000163. DOI: 10.1371/journal.pcbi.1000163

[X11] Gamba A., Ambrosi D., Coniglio A., De Candia A., Di Talia S., Giraudo E., Serini G., Preziosi L., Bussolino F. Percolation morphogenesis and Burgers dynamics in blood vessel formation. Phys. Rev. Lett. 2003;90(11):118101. DOI: 10.1103/PhysRevLett.90.118101

[X12] Keller E.F., Segel L.A. Initiation of slime mold aggregation viewed as an instability. J. Theoret. Biol. 1970;26(3):399-415. DOI: 10.1016/0022-5193(70)90092-5