《PLOS Computational Biology》：Membrane tension and cell–nucleus mechanical coupling regulate confined cell migration
编辑推荐：
细胞迁移是胚胎发育至癌症转移等众多生理过程的基础，主要由质膜（plasma membrane, PM）与细胞核的耦联力学以及细胞骨架动力学共同调控。作为最大且最坚硬的细胞器，细胞核被认为是受限迁移中的限速因素。然而，细胞核的定位因细胞类型和迁移模式而异，且协调
细胞迁移是胚胎发育至癌症转移等众多生理过程的基础，主要由质膜（plasma membrane, PM）与细胞核的耦联力学以及细胞骨架动力学共同调控。作为最大且最坚硬的细胞器，细胞核被认为是受限迁移中的限速因素。然而，细胞核的定位因细胞类型和迁移模式而异，且协调其与质膜力学耦联的力尚未被完全理解。近几十年来，数学模型在支持细胞迁移研究方面日益重要，但模型面临在机制细节与计算效率之间取得平衡的根本挑战。本研究开发了一个基于几何表面偏微分方程（geometric surface partial differential equations, GS-PDEs）的无量纲、生物物理基础模型用于细胞迁移研究。该方法得益于降维处理——质膜和核膜（nuclear envelope, NE）被建模为能量几何膜——并采用演化表面有限元方法（evolving surface finite element method, ESFEM）求解，在机制细节与计算可行性之间取得平衡。该模型再现了实验观察到的迁移模式：脂质储库耗竭和核纤层展开分别限制质膜和核膜的表面积，从而阻断细胞通过狭窄通道；细胞核速度在孔隙出口处达到峰值；细胞核运动在牵拉与推挤模式之间交替。此外，该建模方法提出了关键预测：细胞核穿越孔隙出口时的速度峰值由储存于核膜后部膜张力中的曲率弹性能释放所驱动；质膜—核膜力学耦联可能阻碍细胞移位，因为部分运动能量被投入于细胞核定位。该模型量化了实验上不易获取的力，包括膜张力、弯曲能、渗透压、质膜—核膜弹性及接触力，为检验和提出关于受限细胞迁移的机制假说提供了一个可解释的平台。
受限细胞迁移的力学机制：膜张力与细胞—核耦联的几何表面偏微分方程建模研究解读
一、研究背景与问题
细胞迁移是贯穿胚胎发育到癌症转移的 fundamental 生理过程。在体内，细胞必须穿越组织中的狭窄孔隙和通道，这一过程涉及质膜（plasma membrane, PM）与细胞核的协同变形。细胞核作为细胞中最大且最坚硬的细胞器，被广泛认为是受限迁移中的限速因素。然而，细胞核在不同细胞类型和迁移模式下的定位差异显著，协调其与质膜之间力学耦联的力尚未被完全阐明。传统实验手段难以直接测量膜张力、弯曲能、渗透压等关键力学参数，这限制了研究者对受限迁移机制的理解。
在计算建模方面，现有方法面临根本性挑战：基于个体的模型（agent-based models）虽能捕捉大规模细胞群体动力学，但无法描述细胞形态变化；相场模型（phase-field models）和体积尖锐界面模型虽能有效描述细胞形状演化，但需要完整域描述，计算成本高昂。因此，亟需一种既能保留关键生物物理机制细节又具备计算可行性的建模框架。
二、研究内容与核心贡献
研究人员开发了一个基于几何表面偏微分方程（geometric surface partial differential equations, GS-PDEs）的无量纲建模框架，将质膜和核膜（nuclear envelope, NE）建模为能量几何膜，通过演化表面有限元方法（evolving surface finite element method, ESFEM）数值求解。该框架整合了多种生物物理力：渗透压调控封闭面积、周长控制力模拟脂质储库耗竭、主动力模拟肌动蛋白丝驱动的迁移、质膜—核膜弹性耦联力模拟LINC复合体功能、外部压力以及排斥接触力。模型以无量纲形式表达，各力均以膜张力归一化，使迁移模式可通过一组可解释的无量纲群识别。
该模型成功再现了三个实验观察到的迁移模式：脂质储库耗竭和核纤层展开分别限制质膜和核膜表面积导致细胞被阻断于狭窄通道；细胞核速度在孔隙出口处达到峰值；细胞核运动在牵拉与推挤模式之间交替。更重要的是，模型提出了两个关键机制预测：第一，细胞核穿越孔隙出口时的速度峰值由储存于核膜后部高曲率区域的曲率弹性能释放所驱动；第二，质膜—核膜力学耦联可能阻碍细胞移位，因为牵拉细胞核向前的力会对细胞前端产生反向拉力，增加其延伸阻力，使部分运动能量消耗于细胞核定位。
三、主要技术方法
研究人员采用的核心方法是基于GS-PDEs的建模框架，将质膜和核膜视为二维闭合曲线，其能量由Helfrich泛函描述，包括张力能和弯曲能。通过L2梯度流推导形状演化方程，并结合多种外部力项。数值求解采用演化表面有限元方法，使用二次Lagrange有限元空间离散，通过FEniCSx开源库实现，采用Newton法和MUMPS直接求解器。网格质量通过等分布算法（equidistribution algorithm）维护。研究人员进行了系统的数值收敛性验证，包括曲线缩短流的解析解对比、沙漏形松弛、孔隙穿越和通道阻断四种构型。代码开源发布于GitHub并归档于Zenodo。
四、研究结果
2.1 几何表面PDE框架：建立了质膜和核膜作为能量几何膜的理论框架，推导了Helfrich能量泛函及其一阶变分，给出了形状演化的L2梯度流方程，并系统定义了六类外部力：尺寸控制力（渗透压和周长控制）、主动力、质膜—核膜弹性力、外部压力、排斥接触力（质膜—核膜间、外部障碍物、自接触）。
2.2 无量纲公式化：将模型转化为无量纲形式，引入关键无量纲参数，包括弯曲—张力比、粘度比、张力比和弯曲模量比，使各力均以膜张力归一化，便于识别不同迁移模式的控制参数。
2.3 演化表面有限元离散：采用二次Lagrange有限元空间离散闭合曲线，建立了半离散和全离散弱形式，通过Newton法求解非线性系统，并应用等分布算法维护网格质量。
2.4 数值收敛性验证：通过四种构型验证离散精度——曲线缩短流与解析解对比、沙漏形松弛、孔隙穿越和通道阻断，确认了空间和时间离散的收敛性。
五、讨论与结论
研究人员通过模型量化了实验上不易获取的力学参数，包括膜张力、弯曲能、渗透压、质膜—核膜弹性及接触力，为检验受限细胞迁移的机制假说提供了可解释平台。两个关键预测具有重要生物学意义：核膜后部储存的曲率弹性能释放驱动出口速度峰值，揭示了膜张力在迁移动力学中的主动作用；质膜—核膜力学耦联可能阻碍而非促进迁移，提示细胞需通过锚定于环境基质来将耦联反应力传递给细胞外基质，从而将耦联转化为有效的牵拉机制。该框架具有通用性，适用于任意数量的耦联演化表面，为未来整合更多细胞组分和更复杂微环境提供了扩展基础。
