A mathematical model for neutrophil gradient sensing and polarization

Chemotaxis is a key process in immune response and wound healing, enabling cells to move toward chemical signals. While the general mechanism of cell migration is understood, the precise way a single chemoattractant signal establishes both front and back of the cell remains unclear. Earlier models have proposed the need for global inhibitors to manage polarization, but these have not been confirmed experimentally. Recent biochemical findings suggest localized signaling pathways may be sufficient. This gap motivated the development of a new mathematical model to test whether known biochemical mechanisms can explain gradient sensing and polarization. Prior research has shown that actin and myosin dynamics are essential for movement. However, the integration of these dynamics with chemoattractant signals is not well understood. This paper addresses the need for a model that incorporates experimentally verified biochemical components. The unresolved question of how a cell distinguishes front from back in response to a gradient is central to this work.

The study aims to develop a mathematical model that explains how neutrophils sense chemical gradients and establish cell polarity. The model is based on biochemical mechanisms that have been experimentally characterized. The researchers propose to test whether known dynamics of Rho GTPase and PI3K are sufficient for gradient sensing and polarization. This approach avoids the need for unverified global inhibitors used in previous models. The study also seeks to determine whether the model can function in actin-inhibited cells. The motivation is to provide a framework that aligns with experimental evidence. The model is designed to be robust across a range of parameter values. The ultimate goal is to propose testable experiments that can validate or refine the model.

The researchers constructed a mathematical model using known biochemical interactions in neutrophils. The model integrates Rho GTPase and PI3K activation dynamics. These components were selected based on their established roles in cell polarization. The model simulates how these mechanisms respond to uniform and gradient concentrations of chemoattractants. The model also tests whether polarization occurs in cells with inhibited actin polymerization. The researchers used computational simulations to assess the model's behavior. The model's predictions were evaluated for robustness across parameter variations. The study proposes specific experiments to test the model's predictions in real cells.

The model shows that Rho GTPase and PI3K dynamics are sufficient for gradient sensing and polarization. These mechanisms correctly localize front and rear pathways in response to chemoattractant gradients. The model works even when actin polymerization is inhibited. This suggests that actin is not essential for polarization in this context. The model's predictions remain stable across a wide range of parameter values. A key finding is the proposed coincidence circuit involving PI3K and Ras. This circuit eliminates the need for global inhibitors used in earlier models. The model provides a testable framework for future experimental validation.

The model demonstrates that known biochemical mechanisms can explain gradient sensing and polarization in neutrophils. The researchers propose that global inhibitors are not necessary for this process. The model's predictions are consistent with experimental data on Rho GTPase and PI3K. The model also functions in actin-inhibited cells, suggesting actin is not essential for polarization. The coincidence circuit involving PI3K and Ras is a novel contribution of the model. The model's robustness across parameter variations supports its reliability. The study concludes that further experiments are needed to test the model's predictions. The proposed experiments aim to validate the model and advance understanding of neutrophil chemotaxis.