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

A mass conserved reaction-diffusion system captures properties of cell polarity

Various molecules exclusively accumulate at the front or back of migrating eukaryotic cells in response to a shallow gradient of extracellular signals. Directional sensing and signal amplification highlight the essential properties in the migrating cells, known as cell polarity. In addition to these, such properties of cell polarity involve unique determination of migrating direction (uniqueness of axis) and localized gradient sensing at the front edge (localization of sensitivity), both of which may be required for smooth migration. Here we provide the mass conservation system based on the reaction-diffusion system with two components, where the mass of the two components is always conserved. Using two models belonging to this mass conservation system, we demonstrate through both numerical simulation and analytical approximations that the spatial pattern with a single peak (uniqueness of axis) can be generally observed and that the existent peak senses a gradient of parameters at the peak position, which guides the movement of the peak. We extended this system with multiple components, and we developed a multiple-component model in which cross-talk between members of the Rho family of small GTPases is involved. This model also exhibits the essential properties of the two models with two components. Thus, the mass conservation system shows properties similar to those of cell polarity, such as uniqueness of axis and localization of sensitivity, in addition to directional sensing and signal amplification.Comment: PDF onl

A Comparison of Mathematical Models for Polarization of Single Eukaryotic Cells in Response to Guided Cues

Polarization, a primary step in the response of an individual eukaryotic cell to a spatial stimulus, has attracted numerous theoretical treatments complementing experimental studies in a variety of cell types. While the phenomenon itself is universal, details differ across cell types, and across classes of models that have been proposed. Most models address how symmetry breaking leads to polarization, some in abstract settings, others based on specific biochemistry. Here, we compare polarization in response to a stimulus (e.g., a chemoattractant) in cells typically used in experiments (yeast, amoebae, leukocytes, keratocytes, fibroblasts, and neurons), and, in parallel, responses of several prototypical models to typical stimulation protocols. We find that the diversity of cell behaviors is reflected by a diversity of models, and that some, but not all models, can account for amplification of stimulus, maintenance of polarity, adaptation, sensitivity to new signals, and robustness

Chemotaxis: a feedback-based computational model robustly predicts multiple aspects of real cell behaviour

The mechanism of eukaryotic chemotaxis remains unclear despite intensive study. The most frequently described mechanism acts through attractants causing actin polymerization, in turn leading to pseudopod formation and cell movement. We recently proposed an alternative mechanism, supported by several lines of data, in which pseudopods are made by a self-generated cycle. If chemoattractants are present, they modulate the cycle rather than directly causing actin polymerization. The aim of this work is to test the explanatory and predictive powers of such pseudopod-based models to predict the complex behaviour of cells in chemotaxis. We have now tested the effectiveness of this mechanism using a computational model of cell movement and chemotaxis based on pseudopod autocatalysis. The model reproduces a surprisingly wide range of existing data about cell movement and chemotaxis. It simulates cell polarization and persistence without stimuli and selection of accurate pseudopods when chemoattractant gradients are present. It predicts both bias of pseudopod position in low chemoattractant gradients and-unexpectedly-lateral pseudopod initiation in high gradients. To test the predictive ability of the model, we looked for untested and novel predictions. One prediction from the model is that the angle between successive pseudopods at the front of the cell will increase in proportion to the difference between the cell's direction and the direction of the gradient. We measured the angles between pseudopods in chemotaxing Dictyostelium cells under different conditions and found the results agreed with the model extremely well. Our model and data together suggest that in rapidly moving cells like Dictyostelium and neutrophils an intrinsic pseudopod cycle lies at the heart of cell motility. This implies that the mechanism behind chemotaxis relies on modification of intrinsic pseudopod behaviour, more than generation of new pseudopods or actin polymerization by chemoattractant

Synthetic spatially graded Rac activation drives directed cell polarization and locomotion

Migrating cells possess intracellular gradients of Rho GTPases, but it is unknown whether these shallow gradients themselves can induce motility. Here we describe a new method to present cells with induced linear gradients of active, endogenous Rac without receptor activation. Gradients as low as 15% were sufficient to not only trigger cell migration up the synthetic gradient, but also to induce both cell polarization and repolarization. Response kinetics were inversely proportional to Rac gradient values, in agreement with a new mathematical model, suggesting a role for natural input gradient amplification upstream of Rac. Increases in Rac levels beyond a well-defined threshold dramatically augmented polarization and decreased sensitivity to the gradient value. The threshold was governed by initial cell polarity and PI3K activity, supporting a role for both in defining responsiveness to natural or synthetic Rac activation. Our methodology suggests a general way to investigate processes regulated by intracellular signaling gradients

Adaptive-Control Model for Neutrophil Orientation in the Direction of Chemical Gradients

AbstractNeutrophils have a remarkable ability to detect the direction of chemoattractant gradients and move directionally in response to bacterial infections and tissue injuries. For their role in health and disease, neutrophils have been extensively studied, and many of the molecules involved in the signaling mechanisms of gradient detection and chemotaxis have been identified. However, the cellular-scale mechanisms of gradient sensing and directional neutrophil migration have been more elusive, and existent models provide only limited insight into these processes. Here, we propose a what we believe is a novel adaptive-control model for the initiation of cell polarization in response to gradients. In this model, the neutrophils first sample the environment by extending protrusions in random directions and subsequently adapt their sensitivity depending on localized, temporal changes in stimulation levels. Our results suggest that microtubules may play a critical role in integrating all the sensing events from the cellular periphery through their redistribution inside the neutrophils, and may also be involved in modulating local signaling. An unexpected finding was that model neutrophils exhibit significant randomness in timing and directionality of activation, comparable to our experimental observations in microfluidic devices. Moreover, their responses are robust against alterations of the rate and amplitude of the signaling reactions, and for a broad range in chemoattractant concentrations and spatial gradients

Effects of 3D Geometries on Cellular Gradient Sensing and Polarization

During cell migration, cells become polarized, change their shape, and move in response to various internal and external cues. Cell polarization is defined through the spatio-temporal organization of molecules such as PI3K or small GTPases, and is determined by intracellular signaling networks. It results in directional forces through actin polymerization and myosin contractions. Many existing mathematical models of cell polarization are formulated in terms of reaction-diffusion systems of interacting molecules, and are often defined in one or two spatial dimensions. In this paper, we introduce a 3D reaction-diffusion model of interacting molecules in a single cell, and find that cell geometry has an important role affecting the capability of a cell to polarize, or change polarization when an external signal changes direction. Our results suggest a geometrical argument why more roundish cells can repolarize more effectively than cells which are elongated along the direction of the original stimulus, and thus enable roundish cells to turn faster, as has been observed in experiments. On the other hand, elongated cells preferentially polarize along their main axis even when a gradient stimulus appears from another direction. Furthermore, our 3D model can accurately capture the effect of binding and unbinding of important regulators of cell polarization to and from the cell membrane. This spatial separation of membrane and cytosol, not possible to capture in 1D or 2D models, leads to marked differences of our model from comparable lower-dimensional models.Comment: 31 pages, 7 figure