⚛️ biophysics

A mechanical model of multicellular remodelling in epithelial monolayers

本文在离格点细胞中心模型中引入了一个动态参考系（DRF），用于模拟上皮单层中的亚细胞重塑，从而能够重现多尺度力学响应，并分析变形历史如何影响组织恢复。

原作者： Osborne, J. M., Van Ammers, R., Davit, Y., Gavaghan, D. J., Byrne, H.

发布于 2026-09-14

📖 1 分钟阅读☕ 轻松阅读

原作者： Osborne, J. M., Van Ammers, R., Davit, Y., Gavaghan, D. J., Byrne, H.

原始论文采用 CC BY 4.0 许可（https://creativecommons.org/licenses/by/4.0/）。 ⚕️ 这是一篇未经同行评审的预印本的AI生成解释。这不是医疗建议。请勿根据此内容做出健康决定。阅读完整免责声明

我们身体的外层以及内脏器官的内衬是由上皮组织构成的，这些细胞片层充当着身体的保护屏障和守门人。这些组织并非静止不变的；随着器官的生长、愈合或对环境做出反应，它们不断地重塑自身、拉伸和放松。科学家早已知晓，当这些薄片被拉伸或挤压时，它们并不会简单地弹回或以单一、均匀的动作拉长。相反，它们会表现出两个截然不同的阶段：首先是几秒钟内发生的快速弹性弹跳，随后是可能持续数分钟甚至数小时的缓慢蠕变调整。这种缓慢的变化是一种重塑过程，即细胞内部的机制——特别是支撑和连接它们的蛋白质支架与粘性键——在重新组织以适应新的形状。理解这种双速行为至关重要，因为它揭示了组织如何在压力下维持其完整性，以及它们可能在疾病中如何失效；然而，现有的计算机模型一直难以在不变得异常复杂的情况下捕捉到这种复杂性。

现在，一个研究小组开发出了一种模拟这种行为的新方法，创建了一个既能捕捉快速弹跳又能捕捉缓慢蠕变的数字模型，用以模拟上皮片层。研究人员并没有试图模拟每个细胞内每一个微小的蛋白质（这需要巨大的计算能力），而是引入了一个被称为“动态参考框架”的概念。想象一下，组织是由通过弹簧连接的细胞组成的集合体。在旧模型中，这些弹簧的长度是固结定的；如果你拉伸组织，弹簧会拉伸，然后永远保持那个新的长度，除非细胞彼此分离。这个新模型增加了第二层“隐形”的弹簧，代表了细胞对其形状的内部记忆。当组织被拉伸时，可见的细胞会迅速移动，但这个隐形层会缓慢地改变自身的靶向长度，以匹配新的现实。这种转变代表了细胞重新组织其内部结构以适应其感受到的压力的生物学过程。

研究人员通过运行计算机模拟蠕变和应力松弛实验来测试这一想法，这些是测试材料在力作用下如何表现的标准方法。在蠕变实验中，研究人员对数字组织片的边缘施加恒定的拉力，并观察它随时间变化的拉伸程度。在应力松弛实验中，组织被拉伸到固定宽度并保持在那里，同时测量维持该拉伸状态所需的力。当模型仅包含快速的弹性响应时，组织的表现就像一根简单的橡皮筋：它会迅速拉伸然后停止，或者力会迅速下降然后保持稳定。然而，当研究人员开启动态参考框架以允许亚细胞重塑时，模拟发生了显著变化。组织最初仍会迅速反应，但随后会随时间缓慢持续拉伸，或者力会持续下降，完美地模仿了在真实生物实验中观察到的两阶段行为。

通过调整这个新模型的设置，团队发现这种缓慢重塑的速度取决于细胞对压力的敏感程度以及它们更新内部“记忆”的速度之快。他们发现，如果细胞适应得慢，组织表现得主要像一种简单的弹性材料。但如果细胞很敏感并且会随时间更新其内部结构，组织就会表现出自然界中观察到的复杂的、多速度的响应。至关重要的是，模型还显示了压力的历史非常重要。如果组织在张力下维持很长时间，其内部结构就会发生显著适应。当张力最终释放时，组织并不会立即弹回；相反，它会缓慢恢复，且这种恢复的速率取决于它被拉伸了多久。这是一个关键发现，因为许多旧的计算机模型预测，无论组织被拉伸了多久，恢复情况都是相同的，这一结果与现实世界的观察相矛盾。

这项研究证实，组织形状的缓慢蠕变变化是由细胞内部组件的持续重组驱动的，而不是由细胞分裂或移动到新位置驱动的。研究人员证明，无需模拟细胞内的每一个分子，就可以捕捉到这种行为，这使得研究包含数千个细胞的大型组织片成为可能。他们的工作表明，组织记住并适应过去压力的能力是其机械设计的一个基本属性。通过提供一个将细胞的微观重组与组织的宏观行为联系起来的框架，这个新模型为理解器官如何发育、如何愈合以及如何应对疾病带来的机械力提供了更清晰的路径。用于生成这些发现的代码和数据现已向其他科学家开放，允许学术界构建更详细、更准确的活体组织模型。

技术摘要：上皮单层中多细胞重塑的机械模型

问题陈述
上皮单层构成了许多哺乳动物器官的基础，表现出对外部力的复杂机械响应。实验证据（例如 Harris 等人 [6]）表明，这些组织在多个时间尺度上对变形做出反应：快速的弹性响应（秒级）以及随后的缓慢松弛或蠕变（分钟至小时级）。现有的数学模型难以在不引入细胞增殖或重排的情况下重现这种行为，而这些过程在如此短的时间尺度内是可以忽略不计的。

连续介质模型（例如标准线性固体模型）通常假设存在单一特征时间尺度。
复杂的粘弹性模型（例如多个 Maxwell 元）可以产生多个时间尺度，但通常是现象学的，其参数难以与特定的亚细胞机制联系起来。
基于细胞的模型（例如顶点动力学、重叠球体模型）通常依赖于固定的相互作用法则。当忽略增殖和重排时，它们也往往表现出单时间尺度的响应。虽然亚细胞元件（SE）模型可以通过内部重组来捕捉多时间尺度行为，但其高昂的计算成本使其仅限于极小规模的系统（<10 个细胞），这使得在大规模组织层级的模拟中并不实用。

方法论
作者提出了一种对离格（off-lattice）细胞中心模型（特别是利用 Voronoi 镶嵌/细分框架的 VT 模型）的新颖扩展。其核心创新在于引入了动态参考框架（Dynamic Reference Frame, DRF）。

模型架构： 组织由三组细胞中心表示：
1. 物理细胞 ( $r_i$ )： 实际发生变形的细胞位置。
2. 参考细胞 ( $\rho_i$ )： 一个随时间演化的动态构型，代表组织的“重塑”状态。
3. 静态细胞 ( $s_i$ )： 一个作为基准的固定初始构型。
控制方程：
- 物理运动： 物理细胞根据粘性阻力、外力和连接物理细胞与参考细胞的弹簧力之间的力平衡进行运动（方程 2）。弹簧力取决于物理距离 ( $||r_{ij}||$ ) 与参考距离 ( $||\rho_{ij}||$ ) 之间的差异。
- 参考演化： 参考细胞根据独立的运动方程进行演化（方程 3）。它们受两种竞争力的驱动：
  1. 将其拉向静态构型 ( $s_i$ ) 的恢复力。
  2. 将其拉向物理构型随时间平均值 ( $\bar{r}_{ij}$ ) 的重塑力，其中 $\bar{r}_{ij}$ 是基于记忆窗口 $T_M$ 计算的。
无量纲化： 系统经过无量纲化处理以定义关键参数：
- $\alpha$ （重塑比率）：重塑力相对于恢复力的相对强度。
- $\eta$ （重塑速率）：参考框架相对于物理框架的演化速率。
- $T_M$ （重塑记忆）：用于对物理变形进行平均的时间窗口。
模拟设置： 该模型模拟了在两种实验方案下的二维上皮单层（105 个细胞）：
1. 蠕变松弛： 施加恒定力；记录应变。
2. 应力松弛： 将应变增加至最大值并保持；记录应力。
分析： 时间序列数据通过指数函数之和进行拟合。通过判断是否可以使用单一指数提供简约拟合（误差 $E_1 < 0.05$ ）来确定时间尺度的数量。

核心贡献

动态参考框架 (DRF)： 本文引入了一种机制，将细胞的瞬时变形与其机械参考构型解耦。这使得模型能够模拟亚细胞重塑（细胞骨架重组、肌球蛋白周转、黏附键重塑），而无需显式建模亚细胞组件，从而避免了 SE 模型的高昂计算成本。
多时间尺度的涌现： 作者证明了 DRF 框架能自然地产生多时间尺度响应（快速弹性 + 缓慢松弛/蠕变），这是以往的标准 VT 或细胞中心模型在没有增殖的情况下所缺失的特征。
参数敏感性： 研究量化了重塑参数 ( $\alpha, \eta, T_M$ ) 如何控制单时间尺度与多时间尺度行为之间的转变。
历史依赖性： 该模型成功重现了实验观察到的现象，即释放后的组织恢复取决于施加应变的持续时间，而标准的弹性模型无法捕捉到这一现象。

结果

蠕变实验：
- 在没有重塑的情况下 ( $\alpha=0$ )，组织表现出单时间尺度响应（快速弹性变形后进入平台期）。
- 存在重塑时 ( $\alpha > 0$ )，随着参考框架的演化，会出现第二个较慢的时间尺度，从而允许进一步的变形。
- 多时间尺度的涌现对于施加力的量级是鲁棒的，但对 $\alpha$ （重塑强度）和 $\eta$ （重塑速率）高度敏感。增加记忆窗口 ( $T_M$ ) 也会促进多时间尺度行为。
应力松弛实验：
- 与蠕变类似，重塑引入了二次松弛阶段。
- 升压时间 ( $T_R$ ) 作为控制参数；中间强度的升压时间会暴露不同的后升压松弛时间尺度，而极短或极长的升压时间可能会掩盖这些特征。
- 增加 $T_M$ 会扩大观察到多时间尺度响应的参数空间。
应变恢复：
- 在没有重塑的模型中，释放保持应变后的恢复速率与保持时长 ( $T_H$ ) 无关。
- 在 DRF 模型中，恢复速率与 $T_H$ 成反比。模型预测，较长的保持时间允许参考框架发生更深的变形（代表累积的亚细胞重塑），导致释放时的应力较低且随后的恢复较慢。这种行为符合 Harris 等人 [6] 的普遍实验发现，即恢复取决于变形持续时间，而特定的 $1/T_H$ 关系是该提议模型的定量结果。

意义与主张
本文声称提供了一个在组织尺度研究亚细胞重塑的高效计算框架。通过使用带有 DRF 的扩展标准细胞中心模型，作者可以在不诉诸于极其昂贵的亚细胞模拟或具有不可测量参数的现象学连续介质模型的情况下，重现上皮单层中观察到的复杂多时间尺度机械行为。

作者指出，该框架允许“对生物组织的发展和功能进行更详细的描述和分析”。他们强调，该模型是模块化的：DRF 概念可以应用于其他离格模型（如顶点动力学），并可以通过添加中间参考状态来扩展到维或包含细胞生长和分化等过程。这项工作表明，由亚细胞周转介导的变形历史是组织力学的一个关键因素，现在可以使用这种方法对其进行系统性的研究。

类似论文