Pattern formation in a coupled membrane-bulk reaction-diffusion model for intracellular polarization and oscillations

Understanding how cells establish and maintain spatial organization remains a central challenge in cell biology. While reaction-diffusion models have long been used to explain pattern formation, their application to intracellular processes is still evolving. Prior research has shown that mass-conserved systems can generate stable spatial patterns through Turing instabilities. However, the role of spatial dimensionality and diffusion rates in shaping these patterns is not fully understood. Existing models often simplify cytosolic and membrane-bound dynamics separately. This gap motivated researchers to integrate bulk and surface domains into a single framework. No prior work had resolved how geometry affects oscillatory and stationary patterns in mass-conserved systems. The need for a unified approach to study both membrane-bound and cytosolic components is clear. This paper addresses that need by analyzing a coupled bulk-surface model. The study aims to clarify how spatial constraints influence intracellular polarization.

The study aims to explore how spatial geometry and dimensionality influence pattern formation in intracellular systems. Specifically, it focuses on the behavior of Cdc42, a key polarity protein. The researchers wanted to understand how mass conservation affects the redistribution of molecules between cytosolic and membrane-bound states. They sought to determine whether oscillations or stationary patterns dominate in different geometries. The motivation stemmed from the need to unify bulk and surface dynamics in a single model. The authors aimed to test how dimensionality affects the emergence of anti-phase and in-phase oscillations. They also wanted to assess the role of nonlinear reaction kinetics in pattern formation. By comparing 1-D and 2-D models, they aimed to reveal how spatial constraints shape cell polarization.

The researchers developed a membrane-bulk reaction-diffusion model to study Cdc42 dynamics. They first analyzed a 1-D model representing fission yeast cells. This model included two diffusion equations in the bulk domain and nonlinear ODEs at each end for binding kinetics. They extended the analysis to a 2-D model with circular geometry. In this model, species could diffuse in the cytosol or bind to the membrane. They also considered a nonlocal PDE system to approximate fast bulk diffusion. The models incorporated mass conservation constraints to simulate realistic redistribution. The researchers used bifurcation analysis to identify steady states and oscillatory behaviors. They examined how spatial dimensionality influences symmetry-breaking instabilities. The approach combined analytical and numerical methods to explore pattern formation.

The analysis revealed symmetric and asymmetric steady states in the 1-D model. Anti-phase oscillations were observed between the two ends of the cell. In-phase oscillations were excluded due to mass conservation constraints. The 2-D model showed no radially symmetric oscillations. Instead, symmetry-breaking instabilities led to stationary Turing patterns. Oscillatory Turing instabilities produced traveling and standing waves. Codimension-two Bogdanov-Takens bifurcations occurred when stationary and oscillatory instabilities coincided. These bifurcations caused traveling waves to slow down and become stationary patterns. The results highlight how spatial dimensionality shapes pattern formation. The findings suggest that geometry strongly influences the type of patterns observed.

The study clarifies how geometry and dimensionality affect intracellular pattern formation. The authors propose that mass conservation selects specific spatial modes for redistribution. They suggest that in 1-D, only anti-phase oscillations are possible. In higher dimensions, symmetry-breaking instabilities dominate. The findings indicate that radially symmetric oscillations are excluded in 2-D models. The researchers propose that Turing instabilities lead to stationary or oscillatory patterns. They suggest that codimension-two bifurcations cause transitions between wave and stationary states. The conclusions emphasize the role of spatial constraints in shaping cell polarization. The results provide insights into how geometry influences intracellular dynamics.