跳转至

第 8 章 中心基单细胞模型：基于与胶体粒子概念类比的多细胞组织方法（Center-based Single-cell Models: An Approach to Multi-cellular Organization Based on a Conceptual Analogy to Colloidal Particles）

作者

Dirk Drasdo，本章独立作者。时任德国 University of Leipzig（莱比锡大学）生物信息学教授，后转至 INRIA Paris（法国国家信息与自动化研究院）。本章是他对"中心基模型"框架的系统化论述——也是 Part III（Off-lattice Cell Models）的开篇。

内容概述

本章提出一个根本性的方法论选择：为什么把细胞近似为"各向同性弹性粘性球"是合理的？ 关键论据：(a) 体外悬浮细胞在无刺激下做随机游走，扩散系数 \(D \sim 10^{-7}\) cm²/s——与胶体粒子 Fokker-Planck 方程形式相同；(b) JKR（Johnson-Kendall-Roberts）模型对 ~180 种细胞的粘附力拟合良好（cytoskeleton 未破裂时）；(c) 细胞趋化行为可以通过 Langevin 方程描述，其集体极限给出 Keller-Segel 方程。

章节结构：

  1. Introduction（pp. 171–172）：细胞与胶体粒子的概念类比——包括"集体行为可类比"的论据和"细胞不是被动粒子"的差异。
  2. Model framework（pp. 173–180）：四大类假设（细胞形状、细胞-细胞相互作用、细胞迁移、细胞分裂/死亡）+ 6 个具体变体。
  3. Interaction energy（pp. 176–179）：JKR、Hertz、harmonic-like 三种势能的具体形式。
  4. Mathematical concepts（pp. 180–183）：三种 Metropolis 算法（Scheme I, II, III） + Langevin 方程。
  5. Selected simulation results（pp. 184+）：单层培养、多细胞球体、肠隐窝、胚胎发育。
  6. Discussion（pp. 196+）：与本书其他方法（HDC、CPM）的对比。

读者前置知识：本科弹性力学（Young modulus、Poisson ratio）、JKR 接触理论、Metropolis 算法、连续介质力学。

核心方程与概念

1. JKR 模型（Johnson-Kendall-Roberts）

两个各向同性弹性粘附球的相互作用。给定相对压缩 \(\delta\)（两球形心距离的减少），力 \(F_{ij}\) 由隐式方程 $$ \delta = \frac{a^2}{R} - \left[\frac{16\pi \gamma a^3}{3 E^* R} + \delta^2\right]^{1/2} $$ 其中 \(a\) 是接触圆半径，\(R = R_i R_j / (R_i + R_j)\) 是有效半径，\(E^* = (1 - \sigma_i^2)/E_i + (1 - \sigma_j^2)/E_j\)，\(\gamma\) 是粘附能密度。

关键特征：当 \(d_{ij} = R_i + R_j\)（细胞刚接触）时，立即形成有限接触面积（不像 Hertz 模型）；当拉开时，有 hysteresis 行为——接触在 \(d_{ij} > R_i + R_j\) 时才断裂。这是理论方程 (T)。

2. Hertz 模型（无 hysteresis）

对于小形变，两个弹性球的弹性相互作用能 $$ V_{ij}^{\text{rep}} = \frac{2 E^*}{5} (R_i R_j)^{1/2} (R_i + R_j - d_{ij})^{5/2} $$ 加上 JKR 类似的粘附能 \(-N_{\text{mol}} W_s = -\sigma W_s A_{ij}\)，其中 \(A_{ij}\) 是接触面积。

Hertz 模型没有 hysteresis 行为——接触一旦断裂就立即分开，不如 JKR 物理真实。

3. Harmonic-like 势能

把两个粘附细胞近似为带"不可变形核"的立方体，用线性弹簧连接： $$ V_{ij} = \begin{cases} (2.5 V^{\text{attr}} - V^{\text{attr}}) \cdot 2(1 - \delta_{ij}) & \text{if } -\delta < d_{ij} < 0 \ V^{\text{attr}} & \text{if } 0 \leq d_{ij} < \delta \ \infty & \text{if } d_{ij} < -\delta \ 0 & \text{otherwise} \end{cases} $$ 其中 \(\delta\) 是形变范围，\(V^{\text{attr}}\) 是粘附能。这个模型最简单、计算最快，但对球形细胞不严格。

4. 三种 Metropolis 算法（Scheme I, II, III）

Scheme I：随机选细胞做迁移、定向变化、生长、分裂、死亡试验，Metropolis 接受/拒绝，时间步 \(\Delta t / N\)（\(N\) 是细胞数）。

Scheme II：先建立"待处理"列表，按随机顺序处理，每处理一个从列表中移除，每步时间 \(\Delta t\)。

Scheme III：把不同事件（迁移、定向、生长、分裂、死亡）放入独立列表，按时间间隔 \(\Delta t_m, \Delta t_o, \Delta t_g, \Delta t_d\) 各自调度，保证每种事件的时间间隔是固定的。

Scheme III 最严格——它保证不同事件的相对时间尺度精确可控，适合做参数扫描和模型对比。

5. Metropolis 接受概率

$$ P_{\text{accept}} = \min\left{1, \exp\left(-\frac{\sum_{i<j}(V_{ij}(t + \Delta t) - V_{ij}(t))}{F_r}\right) / T\right} $$ \(F_r\) 是"参考能量"（类比热运动的 \(k_B T\)），但\(F_r\) 不是真正的热能——它反映细胞主动探索环境的"努力"。这里 detailed balance 仍被保留（与 GGH 不同），但这种 detailed balance 在生物系统中是否真的成立是一个开放问题。

6. Langevin 方程

细胞位置 \(\mathbf{r}_i\) 的演化 $$ \gamma \frac{d\mathbf{r}i}{dt} = -\nabla_i U(t) $$ 其中 }} + \mathbf{F}_{\text{noise}\(\gamma\) 是有效摩擦系数，\(\mathbf{F}_{\text{noise}}\) 是白噪声（满足 \(\langle F_{\text{noise}, \alpha}(t) F_{\text{noise}, \beta}(t') \rangle = 2\gamma k_B T \delta_{\alpha\beta} \delta(t - t')\)）。比 Metropolis 更"动力学"——直接求解运动方程。

7. 接触抑制的生物力学形式

定义"形变阈值" \(\epsilon = 0.75\)——若 \(d_{ij} > (l_i + l_j)/2\)（即两个细胞相距超过一半直径），则未形变；否则被压缩、细胞周期延长甚至停滞。

这是"接触抑制"的几何形式——它把"细胞不能分裂"的现象归结为"细胞被挤得太厉害"——避开了显式的"信号通路"建模。

8. 关键概念：参考能量 \(F_r\)

\(F_r\) 是细胞主动迁移的"驱动能量"——类比热运动 \(k_B T\) 但物理含义不同。Beysens et al. (2000) 给出经验值 \(F_r \sim 10^{-16}\) J（约 25 \(k_B T\)）。这远大于真实热能——但细胞的主动运动（细胞骨架重组）确实能产生这么大的力。

9. 细胞分裂的三种算法

  • 算法 (i)：细胞在分裂期间从球形形变为哑铃形，体积不断增长。
  • 算法 (ii)：先以球形增长到目标体积，然后形变为哑铃，体积不变。
  • 算法 (iii)：算法 (ii) 的变种，哑铃长度超过阈值立即分裂为两个独立球——可能产生重叠但很快消失。

所有算法都用小步长 \(d\alpha \ll l\)——保证形变"准静态"。

10. 与 PDE 极限的连接

中心基模型的群体极限在某些条件下可以导出经典 PDE： - 无相互作用 + 随机游走：Fokker-Planck。 - 加入繁殖 - 死亡：Fisher-Kolmogorov-Petrovskii-Piskounov (FKPP) 方程。 - 加入趋化：Keller-Segel 方程。

但作者警告：当细胞-细胞物理相互作用不可忽略时，FKPP 极限不正确——中心基模型的"挤压反馈"无法用线性扩散项描述。这是个体模型不可被 PDE 替代的核心论据。

关键结论

  • 结论 1：细胞在未受信号刺激、悬浮或单层培养状态下，可以近似为"各向同性弹性粘性球"——JKR 模型对 ~180 种细胞的粘附力拟合支持这一假设。
  • 结论 2：Metropolis 算法中，\(F_r \sim 10^{-16}\) J 是参考能量——它类比热能但物理含义不同（细胞主动迁移的"努力"），详细平衡在生物系统中不一定成立。
  • 结论 3：接触抑制的"生物力学形式"——细胞被压缩时（\(d_{ij} < (l_i + l_j)/2\)）停止分裂——是不需要显式信号通路的极简模型，在体外培养中能复现正常细胞在 confluence 时停止分裂的现象。
  • 结论 4：JKR 模型的 hysteresis 行为是细胞粘附的真实特征——细胞一旦接触，立即形成强粘附键；拉开时需要更大距离才能断裂。这种 hysteresis 是组织完整性的物理基础。
  • 结论 5：Scheme III 的"按事件类型分列表 + 固定时间间隔"是最严格的 Metropolis 变种——它保证每种事件类型的速率是精确控制的，适合做参数扫描。
  • 结论 6：中心基模型在体外单层培养、多细胞球体、肠隐窝、胚胎发育 4 个系统中都得到实验验证——这说明"胶体粒子"类比对广泛的生物学问题是适用的。
  • 结论 7：FKPP 极限在"细胞物理相互作用不可忽略"时不成立——这意味着连续 PDE 在描述致密组织时不可靠，必须使用个体基模型。
  • 结论 8：本模型可以扩展到早期胚胎——细胞质量不双倍增长（daughter 比 mother 小），这是"卵裂"型发育的特征。

挑战和开放性问题

  • 挑战 1：3D 形变的精确建模。当前模型把细胞近似为"球"或"哑铃"——真实细胞在 3D 中有复杂的形状变化（纺锤体、丝状伪足、极化）需要更精细的形变模型（如 Ch 12 的椭球模型）。
  • 挑战 2：粘附异质性。当前模型假设所有细胞对的 \(W_s\) 相同——但真实细胞表面 cadherin 表达是空间异质且时间可变的。需要把粘附参数 \(W_s\) 改为动态变量。
  • 挑战 3：细胞骨架动力学。当前模型用"形变阈值 \(\epsilon\)"描述接触抑制——但真实的接触抑制由 p27/Kip1 等细胞周期蛋白调控。需要一个从"分子"到"形变"的桥接。
  • 挑战 4：流体的显式建模。当前模型把"培养基"视为无形阻力的"摩擦"——但在血流、微血管、淋巴管等场景下，流体动力学与细胞运动耦合。这是 Ch 12 / Ch 13 的研究方向。
  • 挑战 5：与单细胞组学的接口。当前模型把"细胞类型"作为离散标签——单细胞 RNA-seq 给出"连续状态空间"。需要一个从"基因表达谱"到"细胞参数"的映射。
  • 挑战 6：参数标定的可复现性。\(F_r\)、\(W_s\)、\(E\)、\(\sigma\) 等参数在不同细胞类型、不同实验条件下都不同——这导致模型的可比性差。
  • 挑战 7：计算规模。2007 年的单 CPU 可以模拟 10⁴–4×10⁵ 个细胞。今天 GPU + 自适应邻域列表可以把规模推到 10⁶。
  • 挑战 8：与非球形细胞类型的整合。神经元、星形胶质细胞、上皮细胞都不是球形——需要为每种细胞类型建立专门的形变模型。

个人反思与批判性分析

本章是"中心基模型"的方法论宣言。它最打动我的是"细胞是胶体粒子"这一概念类比——它用一个具体的、可测量的物理对象（弹性粘附球）替代了抽象的"细胞"——这让建模的参数空间被工业界和物理学界已知的常数所约束。

几点批判性观察：

  1. "细胞是胶体粒子"的概念类比有明确的边界。作者在 §1 末段诚实地承认"细胞不是 Brownian 粒子"——它们有主动运动、主动生长分裂、主动调控。这意味着 JKR 模型只能描述"细胞行为的某个简化子集"——尤其在多细胞通讯（信号通路、力学传导）的场景下，"胶体粒子"类比不适用。

  2. JKR 模型的"hysteresis"对结果的影响需要更系统地研究。作者展示了 JKR 的 hysteresis 行为，但没给出"如果用 Hertz 模型（非 hysteresis）结果会怎样"的系统对比。这是一个被简化的方法论问题。

  3. Scheme III 的"按事件类型分列表"虽严格但与生物学时间尺度不完全对应。例如"细胞分裂"和"细胞迁移"在真实生物中是同时进行的事件——Scheme III 把它们分到不同时刻处理会引入系统偏差。但作者声称"结果对算法选择不敏感"——这一论断需要更广泛的数值实验支持。

  4. "形变阈值 \(\epsilon\)"作为接触抑制模型。这一极简假设避开了所有信号通路——这是它的优势（少参数）也是它的局限（无法解释为什么某些细胞在没有形变时也不分裂，如 p53 介导的细胞周期 arrest）。在 2007 年这是一个合理的简化，今天可能不再足够。

  5. FKPP 极限"在细胞物理相互作用不可忽略时不正确"。这一论断非常重要——它是个体基模型相对于连续 PDE 的核心方法论优势。但作者没有给出"什么时候 FKPP 极限可以"和"什么时候不可以"的清晰判据——这需要更系统的多尺度分析。

  6. "参考能量 \(F_r\)"的生物学含义。作者清楚 \(F_r\) 不是真正的热能，但没有给出 \(F_r\) 与"细胞骨架活性"的具体定量关系。这是所有"中心基"模型的方法论痛点——把"主动"参数化为"被动"的随机数。

  7. 本书的方法谱地位。本章是 Part III 的开篇——与 Ch 1（HDC, 连续场+离散细胞）、Ch 2-3（LGCA/CPM, 格子基）、Ch 11-13（粘弹性连续模型）形成对比。中心基模型的"几何明确"是它的最大优势——细胞之间的接触、挤压、形变都用欧几里得坐标描述，没有"格子伪影"。

  8. 与 Ch 12（Zaman 椭球模型）的关系。本章的球形假设是 Ch 12 椭球模型的退化——椭球模型可以描述"主动形变"，而球形模型不能。这是一个明确的扩展方向：把球形升级为椭球是"最小的可逆"步骤。

  9. "Metropolis vs. Langevin"的方法论选择。作者展示了两种动力学实现，但没给出"何时用 Metropolis、何时用 Langevin"的指导原则。Metropolis 适合做"准确定性"演化（热平衡类比），Langevin 适合做"时间精确"演化（惯性 + 噪声）。这种选择应当与生物学时间尺度（细胞周期、信号传导时间）匹配。

  10. "细胞-细胞相互作用" vs "细胞-基质相互作用"。本章主要关注前者，但体外培养中细胞-基质（collagen、Matrigel）相互作用同样重要。作者在 §1 提到 "substrate is treated as a cell with infinite radius \(R \to \infty\)"——这是一种简化，真实的基质是连续介质，其力学性质（弹性、塑性、粘弹性）需要显式建模。

重要参考文献

[X1] Beysens D.A., Forgacs G., Glazier J.A. Cell sorting is analogous to phase ordering in fluids. Proc. Natl. Acad. Sci. USA 2000;97(17):9467-9471. DOI: 10.1073/pnas.97.17.9467

[X2] Drasdo D., Forgacs G. Modeling the interplay between cell-cell adhesion, cell-cell signaling, and gene regulation in cell sorting. J. Theoret. Biol. 2000;204(4):467-494.

[X3] Johnson K.L., Kendall K., Roberts A.D. Surface energy and the contact of elastic solids. Proc. R. Soc. Lond. A 1971;324(1558):301-313. DOI: 10.1098/rspa.1971.0141

[X4] Hertz H. Über die Berührung fester elastischer Körper. Journal für die reine und angewandte Mathematik 1882;92:156-171.

[X5] Drasdo D., Höhme S. A single-cell-based model of tumor growth in vitro: monolayers and spheroids. Phys. Biol. 2005;2(3):133-147. DOI: 10.1088/1478-3975/2/3/001

[X6] Drasdo D., Höhme S., Block M. Macroscopic limits and mesoscopic models of cell motion in confluent monolayers. In: Collective Dynamics of Interacting Cell Populations. Wiley-VCH, 2011.

[X7] Palsson E., Othmer H.G. A model for individual and collective cell movement in Dictyostelium discoideum. Proc. Natl. Acad. Sci. USA 2000;97(19):10448-10453. DOI: 10.1073/pnas.97.19.10448

[X8] Keller E.F., Segel L.A. Model for chemotaxis. J. Theoret. Biol. 1971;30(2):225-234. DOI: 10.1016/0022-5193(71)90050-6

[X9] Galle J., Sittig D., Hanisch I., Wobus M., Anderer U., Schober A., Ehninger G. Single cell-based mathematical model of in vitro cell migration, cell proliferation and death. Biosystems 2006;83(2-3):107-120.

[X10] Brú A., Albertos S., López J.M., Brú I. Tumor cell proliferation under cyclic hypoxia. (cited in this chapter)

[X11] Schaller G., Meyer-Hermann M. Multicellular tumor spheroid in an off-lattice Voronoi-Delaunay cell model. Phys. Rev. E 2005;71(5):051910. DOI: 10.1103/PhysRevE.71.051910

[X12] Mendoza-Blanco A., Drasdo D. Physical model of cell migration in a 3D microenvironment. (cited in this chapter)

[X13] 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