2025-11-24T20:16:17.055816

Active waves from non-reciprocity and cytoplasmic exchange

Picardo, Jemseena, Kumar

Pattern formation in active biological matter typically arises from the feedback between chemical concentration fields and mechanical stresses. The actomyosin cortex of cells is an archetypal example of an active thin film that displays such patterns. Here, we show how pulsatory patterns emerge in a minimal model of the actomyosin cortex with a single stress-regulating chemical species that exchanges material with the cytoplasm via a linear turnover reaction. Deriving a low-dimensional amplitude-phase model, valid for a one-dimensional periodic domain and a spherical surface, we show that nonlinear waves arise from a secondary parity-breaking bifurcation that originates from the nonreciprocal interaction between spatial modes of the concentration field. Numerical analysis confirms these analytical predictions, and also reveals analogous pulsatory patterns on impermeable domains. Our study provides a generic route to the emergence of nonreciprocity-driven pulsatory patterns that can be controlled by both the strength of activity and the turnover rate.

academic

Active waves from non-reciprocity and cytoplasmic exchange

基本信息

论文ID: 2505.09740
标题: Active waves from non-reciprocity and cytoplasmic exchange
作者: Jason R. Picardo (IIT Bombay), V. Jemseena (ICTS-TIFR), K. Vijay Kumar (ICTS-TIFR)
分类: physics.bio-ph nlin.PS physics.flu-dyn
发表时间: 2025年5月14日 (arXiv预印本)
论文链接: https://arxiv.org/abs/2505.09740

生物学意义：肌动球蛋白皮层的脉动模式是许多细胞过程的基础，如细胞极性建立、胚胎发育等
物理学意义：为理解活性物质中的非平衡模式形成提供了新的理论框架
应用价值：为控制和调节细胞行为提供了潜在的机制

现有方法局限性

现有模型通常需要多个活性物质或复杂的相互作用才能产生振荡模式
大多数理论预测脉动模式源于均匀态的Hopf分岔，但这种机制在单组分系统中难以实现
缺乏对非互惠相互作用在模式形成中作用的深入理解

研究动机

本文旨在证明即使在最简单的单组分活性流体系统中，通过胞质交换机制也能产生脉动模式，并揭示其背后的非互惠相互作用机制。

核心贡献

理论创新：首次证明单一活性物质通过线性周转反应可以驱动脉动模式
机制阐明：揭示了非互惠相互作用导致的奇偶性破缺分岔是产生行波的根本机制
模型构建：建立了适用于不同几何形状（一维周期域、球面）的低维幅度-相位模型
控制策略：发现可以通过调节活性强度和周转率来控制脉动模式的开关
数值验证：通过全模型数值模拟验证了理论预测的准确性

∂c/∂t = -∂(vc)/∂x + D∂²c/∂x² - κ(c - c₀)

力平衡方程：

η∂²v/∂x² - γv = -∂σᵃ(c)/∂x

其中：

c(x,t)：表面浓度场
v(x,t)：流体速度
D：扩散系数
κ：周转率
c₀：基础浓度
η：粘性系数
γ：摩擦系数
σᵃ(c) = ζΔμf(c)：活性应力，f(c) = c/(cₛ + c)

无量纲参数

模型由四个无量纲参数控制：

Péclet数 P ≡ ζΔμ/(γD)：活性对流与扩散的比值
Damköhler数 R ≡ κL²/D：周转率与扩散时间尺度的比值
粘性参数 μ ≡ √(η/(γL²))
饱和参数 S ≡ cₛ/c₀

低维化简模型

采用Galerkin截断方法，将浓度场展开为：

c(x,t) = 1 + ρ₁(t)cos[x - φ₁(t)] + ρ₂(t)cos[2(x - φ₂(t))]

得到三个耦合的常微分方程：

ρ̇₁ = -(R + m₁)ρ₁ + (P/S²)[ρ₁(A₁ - B₁ρ₁² + C₁ρ₂²) - D₁ρ₁ρ₂cosφ]
ρ̇₂ = -(R + m₂)ρ₂ + (P/S²)[ρ₂(A₂ - B₂ρ₁² - C₂ρ₂²) + D₂ρ₁²cosφ]
φ̇ = -(P/S²)(Aφρ₁² - Bφρ₂²)sinφ/ρ₂

其中φ = 2(φ₂ - φ₁)是相位差。