A mathematical model for neutrophil gradient sensing and polarization

Directed cell migration in response to chemical cues, also known as chemotaxis, is an important physiological process involved in wound healing, foraging, and the immune response. Cell migration requires the simultaneous formation of actin polymers at the leading edge and actomyosin complexes at the sides and back of the cell. An unresolved question in eukaryotic chemotaxis is how the same chemoattractant signal determines both the cell's front and back. Recent experimental studies have begun to reveal the biochemical mechanisms necessary for this polarized cellular response. We propose a mathematical model of neutrophil gradient sensing and polarization based on experimentally characterized biochemical mechanisms. The model demonstrates that the known dynamics for Rho GTPase and phosphatidylinositol-3-kinase (PI3K) activation are sufficient for both gradient sensing and polarization. In particular, the model demonstrates that these mechanisms can correctly localize the front'' and rear'' pathways in response to both uniform concentrations and gradients of chemical attractants, including in actin-inhibited cells. Furthermore, the model predictions are robust to the values of many parameters. A key result of the model is the proposed coincidence circuit involving PI3K and Ras that obviates the need for the global inhibitors'' proposed, though never experimentally verified, in many previous mathematical models of eukaryotic chemotaxis. Finally, experiments are proposed to (in)validate this model and further our understanding of neutrophil chemotaxis.