Active waves from non-reciprocity and cytoplasmic exchange

Basic Information

• Paper ID: 2505.09740
• Title: Active waves from non-reciprocity and cytoplasmic exchange
• Authors: Jason R. Picardo (IIT Bombay), V. Jemseena (ICTS-TIFR), K. Vijay Kumar (ICTS-TIFR)
• Classification: physics.bio-ph nlin.PS physics.flu-dyn
• Publication Date: May 14, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2505.09740

Abstract

This study investigates pattern formation mechanisms in active biological matter, with particular focus on the actomyosin cortex—a prototypical active thin film system. The authors establish a minimal model containing a single stress-regulating chemical species that exchanges material with the cytoplasm through linear turnover reactions. By deriving low-dimensional amplitude-phase models applicable to one-dimensional periodic domains and spherical surfaces, the work reveals that nonlinear wave formation originates from a pitchfork bifurcation of the static pattern, which arises from non-reciprocal interactions between spatial modes in the concentration field. Numerical analysis validates theoretical predictions and identifies similar pulsating patterns on impermeable domains.

Research Background and Motivation

Problem Definition

Pattern formation in active biological matter typically originates from feedback mechanisms between chemical concentration fields and mechanical stress. The actomyosin cortex, as a prototypical example of an active membrane, exhibits complex spatiotemporal patterns that are crucial for biological processes such as cell division and morphogenesis.

Research Significance

1. Biological Relevance: Pulsating patterns in the actomyosin cortex form the basis for many cellular processes, including establishment of cell polarity and embryonic development
2. Physical Significance: Provides a novel theoretical framework for understanding non-equilibrium pattern formation in active matter
3. Applied Value: Offers potential mechanisms for controlling and regulating cellular behavior

Limitations of Existing Approaches

1. Existing models typically require multiple active species or complex interactions to generate oscillatory patterns
2. Most theoretical predictions attribute pulsating patterns to Hopf bifurcation of the uniform state, a mechanism difficult to realize in single-component systems
3. Lack of deep understanding regarding the role of non-reciprocal interactions in pattern formation

Research Motivation

This work aims to demonstrate that pulsating patterns can be generated through cytoplasmic exchange mechanisms even in the simplest single-component active fluid systems, and to elucidate the underlying non-reciprocal interaction mechanism.

Core Contributions

1. Theoretical Innovation: First demonstration that a single active species can drive pulsating patterns through linear turnover reactions
2. Mechanism Clarification: Reveals that the pitchfork bifurcation induced by non-reciprocal interactions is the fundamental mechanism for generating traveling waves
3. Model Construction: Establishes low-dimensional amplitude-phase models applicable to different geometries (one-dimensional periodic domains, spherical surfaces)
4. Control Strategy: Discovers that pulsating pattern switching can be controlled by modulating activity strength and turnover rate
5. Numerical Verification: Validates theoretical predictions through full model numerical simulations

Methodology Details

Task Definition

Investigate how spatiotemporal pulsating patterns are generated through cytoplasmic exchange in active fluid thin films containing a single stress-regulating chemical species, and analyze their formation mechanisms and control methods.

Model Architecture

Fundamental Equations

The system is described by two coupled partial differential equations:

Concentration Evolution Equation:

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

Force Balance Equation:

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

Where:

• c(x,t): surface concentration field
• v(x,t): fluid velocity
• D: diffusion coefficient
• κ: turnover rate
• c₀: baseline concentration
• η: viscosity coefficient
• γ: friction coefficient
• σᵃ(c) = ζΔμf(c): active stress, f(c) = c/(cₛ + c)

Dimensionless Parameters

The model is controlled by four dimensionless parameters:

1. Péclet Number P ≡ ζΔμ/(γD): ratio of active advection to diffusion
2. Damköhler Number R ≡ κL²/D: ratio of turnover rate to diffusive timescale
3. Viscosity Parameter μ ≡ √(η/(γL²))
4. Saturation Parameter S ≡ cₛ/c₀

Reduced Low-Dimensional Model

Using Galerkin truncation, the concentration field is expanded as:

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

Yielding three coupled ordinary differential equations:

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

Where φ = 2(φ₂ - φ₁) is the phase difference.

Technical Innovations

1. Non-Reciprocal Interaction Identification: Discovers that the quadratic nonlinear term in equation (5) is ρ₁ρ₂cosφ, while in equation (6) it is ρ₁²cosφ, demonstrating non-reciprocal coupling between modes
2. Parity-Breaking Mechanism: When φ ≠ 0, the concentration distribution loses central symmetry, forming asymmetric traveling wave solutions
3. Secondary Bifurcation Theory: First identification that pitchfork bifurcation of static patterns can directly lead to traveling waves, rather than through traditional Hopf bifurcation
4. Geometric Universality: Demonstrates the universality of this mechanism across different geometries (linear, spherical)

Experimental Setup

Numerical Methods

1. One-Dimensional Case: Finite difference method for solving the PDE system
2. Spherical Case: Spherical harmonic expansion using the SHTns library for numerical computation
3. Parameter Scanning: Systematic scanning of the P-R parameter space to construct phase diagrams

Evaluation Metrics

1. Stability Analysis: Linear stability analysis to determine bifurcation points
2. Wave Speed Measurement: Calculate traveling wave velocity u by tracking concentration peak positions
3. Phase Difference: Measure phase relationships φ between different modes

Comparative Verification

1. Simplified Model vs. Full Model: Compare results from Galerkin truncation model with complete PDE system
2. Different Geometries: Verify mechanism consistency across different domains
3. Parameter Dependence: Analyze behavior under different S and μ values

Experimental Results

Main Results

Phase Diagram Analysis

Four main regions identified in the P-R parameter space:

1. Uniform State: Stable uniform distribution at low P values
2. Static Pattern: Stable spatial patterns at moderate P values
3. Traveling Waves: Moving wave dynamics for specific P-R combinations
4. Complex Oscillations: Complex spatiotemporal dynamics at high P values

Traveling Wave Characteristics

1. Velocity Scaling: Wave speed u increases with P as u ∼ √P
2. Phase Saturation: Phase difference φ saturates as P increases
3. R Dependence: Both u and φ exhibit non-monotonic variation with R, reaching maximum at moderate R values

Geometric Effects

• Waves on Spheres: Manifest as asymmetric spots moving along the equator
• Impermeable Boundaries: Similar pulsating patterns observed on finite domains

Ablation Studies

1. R = 0 Case: No traveling waves exist when R = 0 in the full model, demonstrating necessity of cytoplasmic exchange
2. Single-Mode Truncation: N = 1 Galerkin truncation can only predict static patterns
3. Parameter Robustness: Results verified for robustness to S and μ parameters

Case Analysis

Figure 1 illustrates typical spatiotemporal evolution:

• Static Pattern: Concentration distribution maintains fixed spatial structure
• Traveling Wave: Asymmetric concentration peak propagates at constant velocity
• Complex Oscillations: Complex dynamics including spot nucleation and coalescence

Experimental Findings

1. Counterintuitive Phenomenon: Increasing turnover rate R (typically considered a stabilizing factor) actually promotes wave dynamics under certain conditions
2. Mode Competition: Traveling wave emergence requires comparable amplitudes of two spatial modes (ρ₂/ρ₁ ≈ 1)
3. Symmetry Breaking: φ ≠ 0 is a necessary condition for traveling wave existence

Related Work

Active Matter Theory

This research builds upon the theoretical foundation of active matter fluid mechanics, particularly the work of Marchetti et al. and Jülicher et al., but is the first to realize pulsating patterns in single-component systems.

Pattern Formation Theory

Unlike traditional Turing patterns, the patterns discovered here do not require an activator-inhibitor mechanism but are realized through mechanical-chemical coupling.

Cell Biology Applications

Related experimental observations include:

• PAR protein polarization waves observed by Goehring et al.
• Actin contraction waves discovered by Martin et al.
• Recent observations by Krishna et al. in reconstituted actomyosin networks

Conclusions and Discussion

Main Conclusions

1. Power of Minimal Models: A simple model containing only one active species and linear turnover reactions can generate rich spatiotemporal dynamics
2. Critical Role of Non-Reciprocity: Non-reciprocal interactions between spatial modes are the fundamental mechanism for generating traveling waves
3. Biological Relevance: This mechanism may explain various wave phenomena observed in cells

Limitations

1. Model Simplification: Neglects cell membrane elasticity, three-dimensional effects, and other factors
2. Parameter Range: Certain biologically relevant parameter regions may exceed the model's applicability
3. Experimental Verification: More quantitative experiments are needed to validate theoretical predictions

Future Directions

1. Model Extension: Consider multi-component systems, membrane elasticity, cytoplasmic flow, and other factors
2. Experimental Validation: Design experiments to measure the effects of key parameters such as turnover rate on wave dynamics
3. Application Exploration: Apply the mechanism to other biological systems such as tissue morphogenesis

In-Depth Evaluation

Strengths

1. Theoretical Innovation: First elucidation of the pattern formation mechanism in single-component active systems
2. Mathematical Rigor: Derivation from complete PDE system to simplified model is rigorous, with complete linear stability analysis
3. Comprehensive Numerical Verification: High consistency between simplified and full models enhances conclusion credibility
4. Clear Physical Picture: The physical images of non-reciprocal interactions and parity breaking are clearly articulated

Weaknesses

1. Insufficient Biological Verification: Lacks direct comparison with actual cellular experiments
2. Parameter Estimation: Certain model parameters (such as friction coefficient γ) are difficult to measure accurately in real systems
3. Higher-Order Effects: In strongly nonlinear regions, two-mode truncation may be insufficient

Impact

1. Theoretical Contribution: Provides a new theoretical framework and mechanism for active matter theory
2. Practical Value: Offers guidance for controlling cellular behavior and designing artificial active systems
3. Interdisciplinary Impact: Connects biophysics, fluid mechanics, and nonlinear dynamics

Applicable Scenarios

1. Cell Biology: Explains phenomena such as actin waves and establishment of cell polarity
2. Developmental Biology: Understands spatiotemporal pattern formation in embryonic development
3. Artificial Active Systems: Design active materials with controllable spatiotemporal dynamics

References

This paper cites 53 important references spanning multiple fields including active matter theory, pattern formation, and cell biology. Particularly noteworthy works include:

• Marchetti et al. (2013) - Comprehensive review of active soft matter fluid mechanics
• Bois et al. (2011) - Pattern formation in active fluids
• Kumar et al. (2014) - Pulsating patterns in active fluids
• Coullet & Iooss (1990) - Instability theory of one-dimensional cellular patterns

This paper makes important contributions to theoretical biophysics, not only revealing a new pattern formation mechanism in active matter but also providing new perspectives for understanding and controlling cellular behavior. Its rigorous mathematical analysis and clear physical picture make it an important reference in the field.