跳转至

第 12 章 带细胞粘附和信号传导的三维变形椭球细胞模型（A 3-D Deformable Ellipsoidal Cell Model with Cell Adhesion and Signaling）

作者

Eirikur Palsson，本章独立作者。时任 University of Utah 数学系（与 Fogelson、Dallon 同一所大学——这是 Utah 数学系"生物建模派"的集体成果）。Palsson 与 Othmer 2000 年在 PNAS 发表的"Dd 运动"模型（Ch 1 已引用）是本章的雏形。资助：NSERC RGPIN249789。

内容概述

本章是"椭球 + 粘弹性 + 信号转导"模型的完整文档化版本。它把 Ch 8 的"中心基"模型扩展为"椭球 + 形变"——这是从"点"到"形变体"的关键步骤。

核心模型设计： - 细胞建模为可变椭球（3 个半轴 \(a, b, c\)，保持体积守恒）。 - 粘弹性通过 Kelvin-Voigt 单元（弹簧 + 阻尼器）实现。 - 粘附通过椭球表面"键"实现，键有 on/off 动力学。 - 信号转导通过 Tang-Othmer 7 维 ODE（与 Ch 9 共享）。 - 应用系统：Dictyostelium 发育（从分散阿米巴 → 螺旋波聚集 → 蛞蝓 → 子实体）。

章节结构：

  1. Introduction（pp. 271–275）：为什么用椭球（vs 球或点）——"细胞形状、刚度、信号传导的个体性"是关键。
  2. Mathematical model（pp. 275–288）：方程细节（运动、形变、粘附、信号转导）。
  3. Simulation results（pp. 288+）：Dd 发育的多阶段模拟。
  4. Discussion（pp. 297+）：与实验的对照、模型局限。

读者前置知识：本科粘弹性力学（Kelvin-Voigt 模型）、Langevin 方程、Tang-Othmer 信号转导 ODE（Ch 9）、细胞骨架生物物理。

核心方程与概念

1. 细胞建模：椭球 + 3 个半轴

每个细胞用中心 \(\mathbf{x}_i\)、3 个半轴长 \((a_i, b_i, c_i)\)、3 个单位主轴矢量 \((\mathbf{u}_i, \mathbf{v}_i, \mathbf{w}_i)\) 完全描述。约束：体积守恒 \(V_i = (4/3)\pi a_i b_i c_i = V_0\)（双轴或三轴独立变化即可）。

2. 粘弹性：Kelvin-Voigt 单元

每个半轴用"弹簧 + 阻尼器"并联（图 3）建模。形变动力学： $$ \eta_i \frac{da_i}{dt} = -k_i (a_i - a_{i,0}) + F_a^{(\text{adhesion})} + F_a^{(\text{chemotaxis})} $$ 其中 \(\eta_i\) 是阻尼系数，\(k_i\) 是弹簧刚度，\(a_{i,0}\) 是静止半轴长，\(F_a^{(\cdot)}\) 是来自粘附和趋化的力沿 \(a\) 方向的分量。这是理论方程 (T)。

3. 关键概念：椭球的"拉伸"与"压缩"

当两个椭球接触时，接触区域会"形变"——但作者选择保持椭球形状，只改变 3 个半轴长。这意味着"形变"被限制为 3 自由度（3 个半轴）——比 SEM（\(3N\) 个自由度）少很多。

优点：计算便宜（每个细胞 3 标量 + 3 矢量 = 12 个状态）。 缺点：无法描述"高度非椭球"形状（如丝状伪足、神经元）。

4. 细胞-细胞力：弹性 + 粘附

两椭球 \(i, j\) 接触时，接触区域是椭球表面的局部区域。当相对位移 \(\delta = d_{ij} - (R_i + R_j)\) < 0（重叠）时： - 弹性排斥：\(\mathbf{F}_{ij}^{\text{rep}} = k_{\text{rep}} \delta \mathbf{n}_{ij}\)，方向沿中心连线。 - 粘附吸引：通过 \(\alpha_{IIb}\beta_3\) 类似受体（这里用"细胞-细胞键"）形成——键的 on/off 动力学： $$ k_{\text{on}} = k_{\text{on}}^0 \exp(-E_b / k_B T), \quad k_{\text{off}} = k_{\text{off}}^0 \exp(E_b / k_B T) $$ 其中 \(E_b\) 是键能。键有 finite lifetime（典型毫秒-秒级）。

5. 趋化力

每个细胞感知局部 cAMP 浓度梯度 \(\nabla c\)。如果细胞被激活且有伪足沿某方向，趋化力为： $$ \mathbf{F}_{\text{chem}} = \chi \nabla c(\mathbf{x}_i, t) $$ 其中 \(\chi\) 是趋化灵敏度。关键：伪足只能"拉"（不能"推"）——这意味着细胞只能向有"附着"的方向移动。

6. 关键概念：力平衡

每个细胞的总力必须平衡——细胞不能"自我推进"（违反 Newton 第三定律）。在多细胞聚集中，如果一个细胞拉另一个细胞，第二个细胞会向后拉——如果它没有稳定的"锚"（如基底），就会被拉走。

这与 Ch 8 的"JKR 模型"形成方法论对照——JKR 中力是"对"的；本章显式做了力平衡。

7. 运动方程（过阻尼）

$$ \gamma_i \frac{d\mathbf{x}i}{dt} = \sum_j \mathbf{F}}^{\text{rep}} + \sum_j \mathbf{F{ij}^{\text{adh}} + \mathbf{F} $$ 其中 }} + \mathbf{F}_{\text{random}\(\gamma_i\) 是有效摩擦系数（来自细胞-基质或细胞-细胞粘性）。随机力是白噪声（Langevin 形式）。

8. 信号转导 ODE（Tang-Othmer 7 维）

每个细胞 \(i\) 有 7 个内部状态 \(v_1^i, ..., v_7^i\)（cAMP 信号转导网络的瞬态）。这些 ODE 与 Ch 9 Dallon 共享——即本章与 Ch 9 在信号转导部分有共同的代码基础。

cAMP 反应-扩散方程（3D 推广）： $$ \frac{\partial c}{\partial t} = D \nabla^2 c + \sum_i S_i(\mathbf{x} - \mathbf{x}_i) - \gamma c $$ 插值：cAMP 在 Eulerian 格子上解，细胞从中插值；细胞的输出按权重 \(w\) 分配到格子。

9. 伪足动力学

  • 伪足形成：每 ~3 分钟（生物观测）重新选择方向。
  • 方向选择偏好：原方向 + 趋化梯度方向（加权）。
  • 形成后：维持 ~ 数分钟（persistence time）。
  • 拉力：作用在伪足附着点，力被传递到基底（表面）或邻居细胞（聚集内）。

10. 关键概念："锚"（anchor）问题

当细胞在"内聚集体"内（如 Dd 蛞蝓）时，没有基底——所有力都通过"细胞-细胞粘附"传递。这意味着一个细胞的运动会牵动整个聚集体——这是力学耦合的"长程"效应。

这与 Ch 11（IB 方法）的方法论对照：Ch 11 把"细胞-基底"和"细胞-细胞"分别作为 Navier-Stokes 流体的"边界"——避免了椭球模型的"锚"问题。

关键结论

  • 结论 1：椭球是"球"和"点"之间的合理折中——3 个半轴足够描述大多数细胞形状变化（球 → 拉长 → 圆盘），又比 SEM（256+ 元素）计算便宜。
  • 结论 2：Kelvin-Voigt 单元是模拟细胞粘弹性的最简方式——弹簧（瞬时弹性）+ 阻尼器（粘性耗散）足以捕捉细胞骨架的粘弹性响应。
  • 结论 3：力平衡是"伪足只能拉"的结果——这意味着 Dd 蛞蝓中"孢子前体细胞沿柄细胞向上爬"的运动需要严格的力平衡。这一约束是"流体型"模型（无显式力）所欠缺的。
  • 结论 4：键的 on/off 动力学使粘附成为"动态"过程——这与 cadherin 等真实受体的特性一致。这比"静态 \(J\) 矩阵"的 GGH 模型更接近生物实在。
  • 结论 5：伪足方向重新选择 + 趋化梯度驱动 Dd 聚集的"流"模式——与 Ch 9 Dallon 的"100s 重新定向"规则数学上等价。
  • 结论 6：椭球模型能够复现 Dd 蛞蝓的形成（10⁵ 个细胞）——这需要力平衡、伪足动力学、信号转导 ODE 的精确耦合。
  • 结论 7：细胞刚度可以变化 2 个数量级（纤维细胞 0.6 mdyne/μm vs 红细胞 10⁻³ mdyne/μm）——椭球模型中可以独立设置每个细胞的刚度，这支持"细胞类型多样性"的建模。

挑战和开放性问题

  • 挑战 1：椭球形变自由度的限制。椭球只能"伸长/缩短 3 个轴"——无法描述"丝状伪足"、"分枝形态"（如神经元）、"高曲率突起"。对真核细胞类型的多样性，椭球是过度简化的。
  • 挑战 2：力平衡的数值实现。显式力平衡看似简单，但当多个细胞同时拉动时，可能出现"震荡"（高刚度、低阻尼）。时间步长受 CFL 类约束。
  • 挑战 3：信号转导 ODE 的参数标定。Tang-Othmer 7 维 ODE 仍是"平均细胞"模型——真实 Dd 细胞在聚集过程中表型分化（pre-spore vs pre-stalk）需要双 ODE 系统。
  • 挑战 4：与流体动力学的耦合。当前椭球模型不显式包含流体——这在"细胞在液态介质中运动"时是不完整的。未来应与 Ch 11 的 IB 框架整合。
  • 挑战 5：3D 模拟的规模。10⁵ 个椭球细胞模拟 3D 力学耦合——今天 GPU + 自适应邻域列表可以把规模推到 10⁶ 个细胞。
  • 挑战 6：与单细胞组学的接口。今天的 Dd 单细胞 RNA-seq 显示"pre-spore vs pre-stalk"分化有 100+ 差异基因——椭球模型应当能容纳这些分子差异。
  • 挑战 7：与 Dd 突变体的对照。streamer F、CF 缺失/过表达等突变体（Ch 9 讨论过）——椭球模型中如何复现这些表型？可能需要调节键动力学、伪足频率等。
  • 挑战 8：临床相关系统。肿瘤转移（Ch 1 HDC）、血管生成（Ch 6 GGH）、血液凝固（Ch 11 IB）——椭球模型可以应用于这些系统，但需要适配其特异的细胞类型。

个人反思与批判性分析

本章是"椭球 + 粘弹性 + 信号"建模的完整文档化。它最打动我的是对"力平衡"和"伪足只能拉"这两个约束的显式处理——这让模型有了真正的物理可解释性。

几点批判性观察：

  1. "椭球"是几何上的简化。真实细胞在迁移时是前端 lamellipodium + 后端 uropod 的非对称形状——椭球无法表达这种"头尾分化"。因此椭球模型在描述"细胞极化迁移"时是过度简化的。

  2. "伪足只能拉"是生物学真实的，但模型实现是简化的。真实细胞有 lamellipodium（前推）、收缩丝（后拉）、整合素介导的"机械转导"。椭球模型把这些简化为"沿伪足方向的拉力"——这忽略了"细胞如何在伪足形成时合成 actin 网络"的细节。

  3. "键的 on/off 动力学"是粘附的显式表达——这与 Ch 4-7 的 GGH（"软参数 \(J\) 矩阵"）形成鲜明对比。前者在可证伪性上更优（\(k_{\text{on}}, k_{\text{off}}\) 可以从单分子实验测量）——后者在计算效率上更优。

  4. "细胞刚度 2 个数量级"的事实让模型灵活——但也意味着参数空间大。多细胞系统对"刚度异质性"的敏感性是一个未充分研究的方向。

  5. Tang-Othmer 7 维 ODE 的"平均细胞"特征是 2007 年的局限。今天的 Dd 单细胞 RNA-seq 数据显示 pre-spore 和 pre-stalk 细胞的信号转导网络本身就不同——未来应当用"双 ODE 系统"分别描述两个亚群。

  6. "力平衡"在 Dd 蛞蝓形成中至关重要——这是 GGH 等"流体型"模型（无显式力）所欠缺的。这一方法论优势使椭球模型特别适合研究"力学传导"在发育中的作用。

  7. 与 Ch 8（Drasdo 球形）的对比。Ch 8 的球形是椭球的退化——椭球可以描述"形变"，球形只能描述"重叠"。但球形 + JKR 比椭球 + 简单弹簧在物理上更严格（JKR 有明确的热力学基础）。两种方法各有优势。

  8. 与 Ch 10（Newman SEM）的对比。椭球有 3+3=6 自由度（3 半轴 + 3 方向）；SEM 有 \(3N\) 自由度（\(N\) 个元素）。椭球对"非椭球"形状无能；SEM 可以表达任意形状——但 SEM 没有"显式力"概念。这是本书中"几何 vs 力学"的方法论张力。

  9. 与 Ch 11（IB）的对比。Ch 11 的"血小板聚集体"中没有显式细胞边界——血小板是 Navier-Stokes 流体的边界。这种"无边界"建模在数学上更干净，但在"细胞形状、粘附、信号转导"上不如椭球模型直观。两种方法各有适用场景。

  10. 本书的方法谱地位。本章是 Part IV（粘弹性）的第 2 章——与 Ch 11（IB）、Ch 13（粘弹性组织发育）形成"力学-几何"的方法论三角。椭球模型是其中"最显式的力学 + 最显式的几何"组合——这是它成为 Dd 发育建模主流方法的原因。

重要参考文献

[X1] Palsson E. A 3-D model used to explore how cell adhesion and stiffness can affect cell sorting and movement in multicellular systems. PhD thesis, University of Utah, 1998.

[X2] 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.10408

[X3] Dallon J.C., Othmer H.G. A discrete cell model with adaptive signalling for aggregation of Dictyostelium discoideum. Phil. Trans. R. Soc. Lond. B 1997;352(1351):391-417. DOI: 10.1098/rstb.1997.0029

[X4] Bonner J.T. The Cellular Slime Molds, 2nd ed. Princeton University Press, 1967.

[X5] Bryant P.J., Schmidt O. The genetic control of cell proliferation in Drosophila imaginal discs. J. Cell Sci. Suppl. 1990;13:169-189. DOI: 10.1242/jcs.1990.supplement_13.15

[X6] Davidenko J.M., Pertsov A.V., Salomonsz R., Baxter W.T., Jalife J. Stationary and drifting spiral waves of excitation in isolated cardiac muscle. Nature 1992;355(6358):349-351. DOI: 10.1038/355349a0

[X7] Winfree A.T. When Time Breaks Down: The Three-Dimensional Dynamics of Electrochemical Waves and Cardiac Arrhythmias. Princeton University Press, 1987.

[X8] Lechleiter J., Girard S., Peralta E., Clapham D. Spiral calcium wave propagation and annihilation in Xenopus laevis oocytes. Science 1991;252(5002):123-126. DOI: 10.1126/science.2011747

[X9] Sherratt J.A. Periodic travelling waves in a family of deterministic cellular automata. In: Mathematical and Computer Modelling, 1991. DOI: 10.1016/0895-7177(91)90021-G

[X10] Tomchik K.J., Devreotes P.N. Chemoattractant signaling in Dictyostelium discoideum. Annu. Rev. Cell Dev. Biol. 1999;15:469-517. DOI: 10.1146/annurev.cellbio.15.1.469

[X11] Raper K.B. The Dictyostelids. Princeton University Press, 1984.

[X12] Lee K.J., Cox E.C., Goldstein R.E. Competing patterns of signaling activity in the Dictyostelium signaling pathway. Phys. Rev. Lett. 1996;77(6):1083-1086. DOI: 10.1103/PhysRevLett.77.1083

[X13] Kessin R.H. Dictyostelium: Evolution, Cell Biology, and the Development of Multicellularity. Cambridge University Press, 2001. ISBN: 978-0521583640

[X14] Bonner J.T. The Social Amoebae: The Biology of Cellular Slime Molds. Princeton University Press, 2008.

[X15] Palsson E. A 3-D model used to explore how cell adhesion and stiffness can affect cell sorting and movement in multicellular systems. (cited in this chapter)

[X16] 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.

[X17] Bray D. Cell Movements. Garland Publishing, 1992.

[X18] Pollard T.D., Wu J.Q. Understanding the actin cytoskeleton. (cited in this chapter)

[X19] Fletcher D.A., Mullins R.D. Cell mechanics and the cytoskeleton. Nature 2010;463(7280):485-492. DOI: 10.1038/nature08908

[X20] Engler A.J., Sen S., Sweeney H.L., Discher D.E. Matrix elasticity directs stem cell lineage specification. Cell 2006;126(4):677-689. DOI: 10.1016/j.cell.2006.06.044

[X21] Lauffenburger D.A., Horwitz A.F. Cell migration: a physically integrated molecular process. Cell 1996;84(3):359-369. DOI: 10.1016/S0092-8674(00)81280-5

[X22] Ananthakrishnan R., Ehrlicher A. The forces behind cell movement. Biol. Cell 2007;99(11):547-555. DOI: 10.1042/BC20070066

[X23] Wang Y.L. Exchange of actin subunits at the leading edge of living fibroblasts: possible role of treadmilling. J. Cell Biol. 1985;101(2):597-602. DOI: 10.1083/jcb.101.2.597

[X24] Lin C.H., Espreafico E.M., Mooseker M.S., Forscher P. Myosin drives retrograde F-actin flow in neuronal growth cones. Neuron 1996;16(4):769-782. DOI: 10.1016/S0896-6273(00)80097-5

[X25] Pollard T.D., Borisy G.G. Cellular motility driven by assembly and disassembly of actin filaments. Cell 2003;112(4):453-465. DOI: 10.1016/S0092-8674(03)00120-X

[X26] Ridley A.J., Schwartz M.A., Burridge K., Firtel R.A., Ginsberg M.H., Borisy G., Parsons J.T., Horwitz A.R. Cell migration: integrating signals from front to back. Science 2003;302(5651):1704-1709. DOI: 10.1126/science.1092053

[X27] Parent C.A., Devreotes P.N. A cell's sense of direction. Science 1999;284(5415):765-770. DOI: 10.1126/science.284.5415.765

[X28] Chung C.Y., Firtel R.A. Signaling pathways at the leading edge of chemotaxing cells. J. Muscle Res. Cell Motil. 2002;23(7-8):773-779. DOI: 10.1023/A:1024494818682

[X29] Pollard T.D., Blanchoin L., Mullins R.D. Molecular mechanisms controlling actin filament dynamics in nonmuscle cells. Annu. Rev. Biophys. Biomol. Struct. 2000;29:545-576. DOI: 10.1146/annurev.biophys.29.1.545

[X30] Chen L., Iijima M., Tang M., Landree M.A., Huang Y.E., Xiong Y., Iglesias P.A., Devreotes P.N. Two phases of actin polymerization display different dependencies on PI(3,4,5)P3 accumulation and have unique roles during chemotaxis. Mol. Biol. Cell 2007;18(1):26-38. DOI: 10.1091/mbc.e06-05-0463

[X31] Li L., Cox E.C., Flyvbjerg H. 'Dicty dynamics': chemotaxis, signaling, and directed migration. (cited in this chapter)

[X32] van Haastert P.J.M. How cells use pseudopods for persistent movement and navigation. Sci. Signal. 2011;4(159):pe6. DOI: 10.1126/scisignal.2001108

[X33] Marel A.K., Zorn M., Klingner C., Saban S., Schilde C., Pleschka H., Pálmai Z., Rätz K., Prassler S., Remmler D., Csaba G., Esch T., Plattner H., S. (cited in this chapter regarding the tug-of-war problem in Dd aggregation)