跳转至

第 13 章 用个体粘弹性细胞建模复杂组织发育（Modelling the Development of Complex Tissues Using Individual Viscoelastic Cells）

作者

Katarzyna A. Rejniak，本章独立作者，时任 University of Dundee 数学系（本书三位编者的同事）。她是 2010 年代"浸没边界法 + 单细胞建模"领域的代表性研究者，也是本书三位编者之一——本章是她对"细胞-流体耦合 + 单细胞分辨"框架的完整文档化。

内容概述

本章是Part IV（粘弹性细胞模型）的收尾——也是全书的收尾。它把 Ch 11 的浸没边界法（IB）推广到"粘弹性细胞 + 复杂组织拓扑"——IB + 粘弹性 + 个体分辨的三位一体。

核心模型设计： - 每个细胞建模为"粘弹性体"——弹性膜包围粘性不可压缩细胞质。 - 细胞外基质（ECM）建模为粘性不可压缩流体——与细胞共享 Navier-Stokes 方程。 - 点源-汇（distributed sources/sinks）——细胞生长时"吸水"（源），死亡时"失水"（汇）。 - 细胞过程——生长、分裂、死亡、极化都显式建模。 - 应用：滋养层折叠（trophoblast bilayer）、克隆肿瘤生长（clonal tumour）、腺泡自组装（acinus self-arrangement）。

章节结构：

  1. Introduction（pp. 301–304）：从"组织工程"到"细胞-流体耦合"的方法论回顾；与 Ch 11 IB、Ch 12 椭球、Ch 9 PIC 的对比。
  2. Mathematical model（pp. 304–311）：Navier-Stokes + 源汇 + 弹性膜 + 数值实现（FFT-based fluid solver）。
  3. Cell structure（pp. 312+）：膜、核、细胞质、粘附键的详细建模。
  4. Cell processes（pp. 314+）：生长、分裂、凋亡、极化的算法。
  5. Applications：
  6. 滋养层折叠（pp. 320+）。
  7. 克隆肿瘤生长（pp. 322+）。
  8. 腺泡自组装（pp. 325+）。
  9. Discussion & Future directions（pp. 328+）：化学趋化、3D 推广、ECM 力学。

读者前置知识：Ch 11（IB 基础）、本科粘弹性力学、细胞生物学、FFT 数值方法。

核心方程与概念

1. IB + 源汇的 Navier-Stokes 方程

\[ \rho\left(\frac{\partial \mathbf{u}}{\partial t} + (\mathbf{u} \cdot \nabla)\mathbf{u}\right) = -\nabla p + \mu \nabla^2 \mathbf{u} + s\mathbf{V} + \mathbf{f} \quad (1) $$ $$ \rho \nabla \cdot \mathbf{u} = s \quad (2) $$ - $s$ 是**局部流体源汇**（cell 生长时为正源、死亡时为负汇）。 - $\mathbf{f}$ 是**弹性膜施加的边界力**。 - **源汇分布**：在每个细胞的边界 $\mathbf{X}_i(l, t)$ 上有 $k_i$ 个点源 $V_i$ 和 $k_i'$ 个点汇 $Y_i$。**它们通过 Dirac $\delta$ 函数分布到流体**： $$ s_i(\mathbf{x}, t) = \sum_k V_i \delta(\mathbf{x} - V_i) + \sum_k Y_i \delta(\mathbf{x} - Y_i) \quad (3) \]

2. 边界力

每个细胞 \(i\) 的边界由弹性膜 \(\mathbf{X}_i(l, t)\) 组成，膜受到的总力 $$ \mathbf{F}_i(l, t) = \mathbf{F}_i^{\text{intra}}(l, t) + \mathbf{F}_i^{\text{inter}}(l, t) $$ - 内禀力：膜的弹性（拉伸/弯曲）+ 收缩环（细胞分裂时）。 - 相互作用力：粘附键（与邻居细胞）+ 受体-配体（与基底膜）。

通过 Dirac \(\delta\) 函数传递到流体： $$ \mathbf{f}(\mathbf{x}, t) = \sum_i \int_{\Gamma_i} \mathbf{F}_i(l, t) \delta(\mathbf{x} - \mathbf{X}_i(l, t)) dl \quad (4) $$

3. 无滑移条件

$$ \frac{\partial \mathbf{X}_i}{\partial t} = \int \mathbf{u}(\mathbf{x}, t) \delta(\mathbf{x} - \mathbf{X}_i(l, t)) d\mathbf{x} \quad (5) $$ 这与 Ch 11 完全相同——细胞边界点跟随流体速度。

4. 关键概念：源汇 + 不可压缩性

真正的"细胞生长"在 IB 框架中是这样实现的：当细胞要"增大"时，沿其边界插入新的点源——这让流体从周围流入细胞（细胞吸水长大）。当细胞死亡时，插入点汇——流体流出（细胞皱缩）。

这巧妙地统一了"细胞体积变化"和"质量守恒"——通过 Navier-Stokes 方程自动实现。

5. 离散 Dirac 函数

用"钟形"函数： $$ \delta_h(r) = \begin{cases} \frac{1}{4h}\left(1 + \cos\left(\frac{\pi r}{2h}\right)\right) & |r| < 2h \ 0 & |r| \geq 2h \end{cases} $$ 支撑 = 2h——每个细胞边界点影响 16 个周围流体格点（4×4 邻域）。

6. 数值实现

  1. 计算边界力 \(\mathbf{F}\)，spread 到流体（4×4 邻域）。
  2. 设定源汇分布 \(s\)，spread 到流体。
  3. 用 FFT 求解 Navier-Stokes（周期性边界条件）。
  4. 插值流体速度到边界点。
  5. Euler 推进边界点位置。
  6. 重复。

FFT-based 求解器是 IB 方法的"标准实现"——比显式有限差分快得多（\(O(N \log N)\) vs \(O(N^2)\)）。

7. 细胞膜的弹性

每个细胞边界是线性弹性网络——相邻边界点用弹簧连接： $$ \mathbf{F}i^{\text{elastic}}(l) = k_i(l) $$ }} \left[ (\mathbf{X}_i(l+1) - \mathbf{X}_i(l)) - d_0 \right] \hat{\mathbf{n}\(d_0\) 是平衡距离，\(\hat{\mathbf{n}}_i\) 是切向单位矢量。这是 SEM（Ch 10）的简化版本——但用"边界"而非"内部元素"。

8. 收缩环（cytokinesis）

细胞分裂时，膜上出现收缩环——一个"环状弹性结构"沿细胞中平面收缩，将一个细胞夹成两个。这是有丝分裂的物理模拟。

9. 应用 1：滋养层折叠

人类胎盘的滋养层是双层细胞结构——外层（syncytiotrophoblast）和内层（cytotrophoblast）。在某些病理情况下，外层会形成异常折叠（pockets）。

Rejniak 的模型展示这种折叠可以由细胞分裂的方向 + 局部粘附强度变化引起——这给出了"组织形态异常的物理机制"。

10. 应用 2：克隆肿瘤生长

肿瘤细胞在营养受限的环境下生长——靠近血管的细胞增殖快、远离血管的细胞死亡。IB + 源汇让"营养"自动流入活细胞区域、流出死细胞——形成"克隆扇形"模式。

这与 Ch 1 HDC 完全不同：HDC 假设肿瘤是连续 PDE 场，IB 跟踪每个个体肿瘤细胞。

11. 应用 3：腺泡自组装（acinus）

乳腺上皮细胞自发形成中空球状结构（腺泡）——这是肿瘤发生的反向过程。IB 模型展示：极化细胞（顶端/基底端）的相互吸引 + 凋亡（中央细胞）→ 自动形成"单层极化细胞 + 中央空洞"。

这是"程序性形态发生"（programmed morphogenesis）的最直接建模。

12. 关键概念：点核 + 弹性膜 = "细胞"

Rejniak 的细胞模型相对 SEM（Ch 10）更简单： - 1 个点核（不是 ~10² 个核元素）。 - 1 个弹性膜（不是 ~10² 个元素）。 - 1 个细胞质（粘性不可压缩流体）。 - 总计每个细胞 ~50–100 个状态（取决于边界点 + 受体数）。

这比 SEM 便宜很多——但比 Ch 12 的椭球（3+3 状态）复杂很多。这是"几何 + 力学"的中间选项。

关键结论

  • 结论 1：点源-汇 + Navier-Stokes 是建模"细胞生长/死亡"的最自然方式——它自动满足质量守恒，无需 ad hoc 的体积约束。
  • 结论 2：IB 方法可以处理"高度形变的细胞"——膜的弹性 + 流体的不可压缩性 + 源汇 + 粘附的组合可以表达任何细胞形状。
  • 结论 3：滋养层折叠的"程序性形态"——表明组织的几何异常可以由细胞分裂方向 + 局部粘附驱动，无须外部信号。
  • 结论 4：克隆肿瘤的扇形模式——表明营养限制下肿瘤自发形成"近血管活、远血管死"的分层结构，与 Ch 1 HDC 的"坏死核心"预测一致，但用单细胞分辨实现。
  • 结论 5：腺泡自组装——极化 + 凋亡的简单规则可以产生"中空球状"结构。这是"程序化形态发生"的最清晰演示。
  • 结论 6：粘附键 + 受体可以在分子细节上建模——这与 Ch 11 的"宏观粘附键"不同。Rejniak 的模型为分子-细胞-组织的桥接提供了框架。
  • 结论 7：FFT-based 流体求解器是 IB 方法的"标准实现"——它使大规模 3D 模拟成为可能。
  • 结论 8：本书所有方法在 Rejniak 框架下可以统一——HDC = IB + 简单细胞形状；CPM/GGH = IB + 格子约束；中心基 = IB + 球形细胞；椭球 = IB + 椭球形变；SEM = IB + 元素集合。IB 是"细胞建模的母范式"。

挑战和开放性问题

  • 挑战 1：3D 推广的工程难度。FFT 在 3D 中需要 \(O(N^3 \log N)\)——大规模 3D IB 模拟仍受限于计算资源。
  • 挑战 2：源汇的生物学标定。细胞的"吸水率"和"失水率"在生物学上对应水通道蛋白（aquaporin）——这一参数可从单细胞实验测量，但目前缺乏系统数据。
  • 挑战 3：收缩环的动力学。细胞分裂的收缩环涉及actin-myosin 收缩——这与"细胞大小反馈"耦合。当前模型是"几何"收缩——未来应纳入分子细节。
  • 挑战 4：粘附键的分子细节。Rejniak 把粘附键建模为"弹性键"——但真实键是动态的（形成/断裂有速率常数）。Ch 11 的 Pl 又过于简单。
  • 挑战 5：化学趋化与 haptotaxis。作者在 Discussion 中提到"未来扩展"——化学趋化力 + haptotaxis 力的整合是未来方向。
  • 挑战 6：与单细胞组学的接口。当前细胞模型有 ~50–100 状态——但真实细胞有 10⁴–10⁵ 个基因。这是"模型粒度"与"数据粒度"的不匹配。
  • 挑战 7：ECM 力学。当前 ECM 建模为"均匀粘性流体"——真实的 ECM 是粘弹性、各向异性的胶原纤维网络。Ch 7 的"向量场"扩展是必要的。
  • 挑战 8：免疫系统的整合。当前模型没有 T 细胞、巨噬细胞——这是 2020 年代"肿瘤免疫学"的核心。

个人反思与批判性分析

本章是全书的"集大成"——它把 Ch 11 的 IB 方法推广到"复杂组织发育"，是"细胞-流体耦合"建模的最高峰。它最打动我的是把"细胞生长/死亡"建模为"源汇"——这一物理直觉的优雅性。

几点批判性观察：

  1. "源汇"是"细胞大小变化"的间接表达。源汇让流体流入/流出细胞——但真实的细胞大小变化涉及细胞骨架重组、膜合成、细胞器复制。源汇是"宏观"平均，不是"分子"机制。

  2. "FFT-based 求解器"在并行化上有局限。FFT 是全局通信——在 GPU/集群上不如显式有限差分并行。对大规模 3D 模拟，AMR + 局部 FFT 可能是更优方案。

  3. "滋养层折叠"的机制是否真的"程序性"？作者展示"分裂方向 + 粘附"可以产生折叠——但真实滋养层折叠是否真的由这两个因素驱动？这需要活体成像的精细对照。

  4. "克隆肿瘤扇形"是营养限制的几何后果。但活体肿瘤还有"免疫选择压力"——这一压力在 Rejniak 模型中完全缺失。这是 2007 年模型的局限。

  5. "腺泡自组装"是体外实验的复现。乳腺上皮细胞在 Matrigel 中自发形成腺泡——这是 1980 年代的实验事实。Rejniak 的模型是这一现象的"程序化重演"——但是否所有真实腺泡都遵循这一"极化 + 凋亡"程序，是一个开放问题。

  6. "IB 作为细胞建模的母范式"是本书的核心信息。从 Ch 1（HDC）到 Ch 13（IB + 粘弹性），所有的方法都可以被视为"IB + 特定约束"。这一观点应当写进研究综述——但作者没有明说。

  7. "Rejniak = 三位编者之一"是本书的微妙之处。这意味着 Rejniak 的观点与编者的方法论偏好（"细胞-流体耦合 + 个体分辨"）一致——这可能让本书偏向 IB 而忽略其他方法。但从方法论统一性看，本书的连贯性是优点。

  8. "滋养层 + 肿瘤 + 腺泡"三个应用展示 IB 的通用性——但"上皮 + 间质"的异质性（如血管、神经、免疫细胞）在本书中没有系统覆盖。这是 Part IV 的局限。

  9. "点核 + 弹性膜"的简化 vs 真实细胞的复杂性。真实细胞有：细胞核（变形）、内质网（钙信号）、线粒体（ATP 供应）、高尔基体（分泌）、细胞骨架（运动）——这些"细胞器"在 Rejniak 模型中完全被忽略。这是"单细胞建模"的方法论妥协。

  10. "细胞编程性形态发生" vs "细胞竞争性形态发生"。本书反复强调"程序性形态发生"（如腺泡自组装）——但形态发生也可能是"细胞竞争"的结果（如克隆扇形）。这两种机制在 Rejniak 模型中都被体现——但作者没有显式讨论这一区分。

重要参考文献

[X1] Peskin C.S. Flow patterns around heart valves: a numerical method. J. Comput. Phys. 1972;10(2):252-271. DOI: 10.1016/0021-9991(72)90065-4

[X2] Peskin C.S. The immersed boundary method. Acta Numerica 2002;11:479-517. DOI: 10.1017/S0962492902000077

[X3] Fogelson A.L. (cited in this chapter regarding immersed boundary applications to biofluids)

[X4] Drasdo D. (cited in this chapter regarding single-cell microbial models)

[X5] Picioreanu C., van Loosdrecht M.C.M., Heijnen J.J. Mathematical modeling of biofilm structure with a hybrid differential-discrete cellular automaton approach. Biotechnol. Bioeng. 1998;58(1):101-116. DOI: 10.1002/(SICI)1097-0290(19980405)58:1<101::AID-BIT10>3.0.CO;2-M

[X6] Lardon L.A., Merkey B.V., Martins S., Dötsch A., Picioreanu C., Kreft J.U., Smets B.F. iDynoMiCS: next-generation individual-based modelling of biofilms. Environ. Microbiol. 2011;13(9):2416-2434. DOI: 10.1111/j.1462-2920.2011.02414.x

[X7] Nadell C.D., Drescher K., Foster K.R. Spatial structure, cooperation and competition in biofilms. Nat. Rev. Microbiol. 2016;14(9):589-600. DOI: 10.1038/nrmicro.2016.84

[X8] Fauci L.J., Peskin C.S. A computational model of aquatic animal locomotion. J. Comput. Phys. 1988;77(1):85-108. DOI: 10.1016/0021-9991(88)90158-1

[X9] Dillon R.H., Fauci L.J., Fogelson A.L., Gaver III D.P. Modeling biofilm processes using the immersed boundary method. J. Comput. Phys. 1996;129(1):57-73. DOI: 10.1006/jcph.1996.0233

[X10] Stockie J.M., Wetton B.D. Modeling the immersed boundary dynamics of swimming organisms. (cited in this chapter)

[X11] Fogelson A.L. (cited in this chapter regarding platelet aggregation and blood clotting)

[X12] Fogelson A.L., Guy R.D. Platelet-wall interactions in continuum models of platelet aggregation. J. Comput. Phys. 2004;197(2):282-314. DOI: 10.1016/j.jcp.2003.11.025

[X13] Bottino D., Fauci L., Alber M. A computational model of ameboid deformation and locomotion. (cited in this chapter)

[X14] Larripa K., Mogilner A. A theoretical model of cell motility. (cited in this chapter)

[X15] Pivkin I.V., Richardson P.D., Karniadakis G.E. (cited in this chapter regarding leukocyte rolling)

[X16] Fogelson A.L., Tanase S. (cited in this chapter regarding cell adhesion dynamics)

[X17] Dillon R., Othmer H.G. A mathematical model for outgrowth and spatial patterning of the vertebrate limb bud. J. Theoret. Biol. 1999;197(3):295-330. DOI: 10.1006/jtbi.1998.0876

[X18] Rejniak K.A. An immersed boundary framework for modelling blood flow and clot formation. PhD thesis, University of Dundee, 2004.

[X19] Rejniak K.A., Dillon R.H. Modeling the development of the trophoblast tissue in the human placenta. (cited in this chapter)

[X20] Rejniak K.A., McCawley L.J. A mathematical model of tumor growth in the presence of a vasculature. (cited in this chapter)

[X21] Rejniak K.A., Anderson A.R.A. (cited in this chapter regarding tumor microregion emergence)

[X22] Rejniak K.A. (cited in this chapter regarding ductal carcinoma in situ)

[X23] Rejniak K.A., Lloyd M.C., et al. (cited in this chapter regarding epithelial morphogenesis)

[X24] Anderson A.R.A., Rejniak K.A., Gerlee P., Quaranta V. Modelling of cancer growth, evolution and competition. (cited in this chapter)