Positive Feedback Regulation Results in Spatial Clustering and Fast Spreading of Active Signaling Molecules on a Cell Membrane

Positive feedback regulation is ubiquitous in cell signaling networks, often leading to binary outcomes in response to graded stimuli. However, the role of such feedbacks in clustering, and in spatial spreading of activated molecules, has come to be appreciated only recently. We focus on the latter, using a simple model developed in the context of Ras activation with competing negative and positive feedback mechanisms. We find that positive feedback, in the presence of slow diffusion, results in clustering of activated molecules on the plasma membrane, and rapid spatial spreading as the front of the cluster propagates with a constant velocity (dependent on the feedback strength). The advancing fronts of the clusters of the activated species are rough, with scaling consistent with the Kardar-Parisi-Zhang (KPZ) equation in one dimension. Our minimal model is general enough to describe signal transduction in a wide variety of biological networks where activity in the membrane-proximal region is subject to feedback regulation.

💡 Research Summary

The paper investigates how positive feedback combined with slow diffusion can generate spatial clustering and rapid propagation of active signaling molecules on a cell membrane. Using a minimal reaction scheme inspired by Ras activation—where Z represents Ras‑GDP, X represents Ras‑GTP, and Y represents the guanine‑exchange factor SOS—the authors model three processes: basal activation (Z + Y → X + Y), feedback‑enhanced activation (Z + X + Y → 2X + Y), and SOS removal (Y → ∅). The reaction rates (k₁, k₂, k₃) and diffusion constant D are chosen to match experimentally measured values for Ras and SOS, with diffusion on the membrane being orders of magnitude slower than in the cytosol.

A kinetic Monte Carlo algorithm on a lattice is employed in both one‑ and two‑dimensional geometries. In 1‑D, topological constraints cause Y particles to become trapped between X particles, leading to a frozen state; this scenario is biologically irrelevant for a membrane but illustrates the importance of dimensionality. In 2‑D, the system exhibits nucleation of X‑rich domains, growth of these domains, and eventual takeover of the whole lattice by X particles.

The authors define a critical nucleus size l_c by balancing diffusion‑driven loss (τ_D ≈ l²/D) against feedback‑driven gain (τ_fb ≈ 1/(k₂ ρ_Y ρ_Z)). Simple dimensional analysis predicts l_c ≈ √(D/(k₂ ρ_Y ρ_Z)). Simulations confirm that larger diffusion constants increase l_c, but quantitative deviations arise due to non‑linearity and stochastic fluctuations. The dynamic structure factor S(q,t) is used to extract the critical wave‑vector q_c = 2π/l_c; curves of S(q,t) at different times intersect at q_c after correcting for a uniform vertical shift caused by the continual increase of X particles.

The fraction of X particles, f_X(t) = N_X(t)/M, follows a sigmoidal rise and saturates at unity when k₃ = 0 (no SOS removal). By assuming a fixed number n of circular nuclei that grow with radius R(t), the authors derive df_X/dt = k n dA/dt (1‑f_X), where A = πR². Integrating yields f_X(t) = 1 − (1‑f_X(t₀)) exp