A computational model of chemotaxis-based cell aggregation

Dynamic cluster field modeling of collective chemotaxis

Collective migration of eukaryotic cells is often guided by chemotaxis, and is critical in several biological processes, such as cancer metastasis, wound healing, and embryogenesis. Understanding collective chemotaxis has challenged experimental, theoretical and computational scientists because cells can sense very small chemoattractant gradients that are tightly controlled by cell-cell interactions and the regulation of the chemoattractant distribution by the cells. Computational models of collective cell migration that offer a high-fidelity resolution of the cell motion and chemoattractant dynamics in the extracellular space have been limited to a small number of cells. Here, we present Dynamic cluster field modeling (DCF), a novel computational method that enables simulations of collective chemotaxis of cellular systems with O ( 1000 ) cells and high-resolution transport dynamics of the chemoattractant in the time-evolving extracellular space. We illustrate the efficiency and predictive capabilities of our approach by comparing our numerical simulations with experiments in multiple scenarios that involve chemoattractant secretion and uptake by the migrating cells, cell migration in confined spaces, regulation of the attractant distribution by cell motion, and interactions of the chemoattractant with an enzyme. The proposed algorithm opens new opportunities to address outstanding problems that involve collective cell migration in the central nervous system, immune response and cancer metastasis.

A mathematical model of cell movement and clustering due to chemotaxis

This paper presents a numerical method for modelling cell migration and aggregation due to chemotaxis where the cell is attracted towards the direction in which the concentration of a chemical signal is increasing. In the model presented here, each cell is represented by a system of springs connected together at node points on the cell's membrane and on the boundary of the cell's nucleus. The nodes located on a cell's membrane are subject to a force which is proportional to the gradient of the concentration of the chemical signal which mimics the behaviour of the chemical receptors in the cell's membrane. In particular, the model developed here will consider what happens when two (or more) cells collide and how their membranes connect to each other to form clusters of cells. The methods described in this paper will be illustrated with a number of typical examples simulating cells moving in response to a chemical signal and how they combine to form clusters.

Modelling the motion of clusters of cells in a viscous fluid using the boundary integral method

In experiments clusters of cells are often observed to move in response to a chemical signal which is present in the fluid surrounding the cells. This process is known as chemotaxis. This paper presents a method for modelling the motion of clusters of cells moving through a viscous fluid in response to a known chemical signal using a boundary integral formulation of the governing equations rather than the more usual differential equation formulation. The numerical results presented in this paper show that the boundary integral method can be used to simulate the motion of cell clusters through the fluid. The results of the simulations are compared to some experimental observations of cell and cluster motion.

A simple mathematical model of cell clustering by chemotaxis

Chemotaxis is the process by which cells and clusters of cells follow chemical signals in order to combine and form larger clusters. The spreading of the chemical signal from any given cell can be modeled using the linear diffusion equation, and the standard equations of motion can be used to determine how a cell, or cluster of cells, moves in response to the chemical signal. The resulting differential equations for the cell locations are integrated through time using the fourth-order Runge-Kutta method. The effect which changing the initial concentration magnitude, diffusion constant and velocity damping parameter has on the shape of the final clusters of cells is investigated and discussed.

Computational modeling of epithelial-mesenchymal transformations

An epithelial-mesenchymal transformation (EMT) involves alterations in cell-cell and cell-matrix adhesion, the detachment of epithelial cells from their neighbors, the degradation of the basal lamina and acquisition of mesenchymal phenotype. Here we present Monte Carlo simulations for a specific EMT in early heart development: the formation of cardiac cushions. Cell rearrangements are described in accordance with Steinberg's differential adhesion hypothesis, which states that cells possess a type-dependent adhesion apparatus and are sufficiently motile to give rise to the tissue conformation with the largest number of strong bonds. We also implement epithelial and mesenchymal cell proliferation, cell type change and extracellular matrix production by mesenchymal cells. Our results show that an EMT is promoted more efficiently by an increase in cell-substrate adhesion than by a decrease in cell-cell adhesion. In addition to cushion tissue formation, the model also accounts for the phenomena of matrix invasion and mesenchymal condensation. We conclude that in order to maintain epithelial integrity during EMT the number of epithelial cells must increase at a controlled rate. Our model predictions are in qualitative agreement with available experimental data.