Guided by Signals: The Mathematics of Chemotaxis
by Hewan Shemtaga
Advisors: Dr. Wenxian Shen & Dr. Selim Sukhtaiev
Figure 1: Illustration of Chemotaxis in the presence of a Chemo-Attractant
Where u denotes the density of the cell/organism and v denotes the concentration of the chemical. In the above system of equations, (0.1) describes the dependence of the cell density variation over time on the migration of the population, the growth and death of the population, and the movement of the cell/organism towards or away from the chemical gradient. On the other hand, equation (0.2) represents the dependence of the chemical concentration variation over time on the diffusion of the chemical as well as the production and degradation of the chemical signal.
Our research concerns the above chemotaxis model on a one-dimensional graph (network) structure. Although the study of mathematical models on networks goes back to the 1980s, it has seen a recent surge in interest, primarily attributed to its numerous applications in the fields of biology and neurobiology. One notable application is in the domain of tissue engineering, concerning the motion of fibroblasts, the stem cells responsible for the regeneration of dermal tissue, on synthetic scaffolds during dermal wound healing. Simulating this process using a chemotaxis model on a network is suitable because the edges of the network symbolize the fibers of the scaffold and the underlying transport equations provide the evolution of fibroblast densities along each individual fiber . We began by examining the validity of the general chemotaxis model on a network structure. We then proved the existence of a unique solution whose behavior changes continuously with the initial conditions. Furthermore, under realistic biological conditions, we proved solutions stay within expected biological limits at all times. Subsequently, we carried out a mathematical analysis to investigate the behavior of cells and microorganisms in response to chemotaxis where cells/organisms follow a logistic growth. Our research aimed to answer a range of biologically relevant questions, including analyzing the dynamics of the system as the chemotaxis sensitivity crosses certain thresholds. Specifically, we aimed to understand whether there are distinct shifts in the system dynamics at certain chemotaxis sensitivity thresholds, commonly referred to as bifurcation points. Our findings revealed a shift from a constant stable solution when chemotaxis sensitivity is weak, to a bifurcation state characterized by existence of other nontrivial stationary solutions in the presence of strong chemotaxis sensitivity. This outcome implies, while weak chemotaxis sensitivity does not significantly alter long-term system dynamics, the dominance of chemotaxis sensitivity over logistic growth introduces a variety of intricate and noteworthy patterns.
Figure 2: (a) and (c) show initial population of cell/bacteria on a simple network structure. (b) With weak chemotaxis sensitivity, the population stabilizes at the carrying capacity uniformly in space. In contrast, with strong chemotaxis sensitivity, initial population may converge to nonconstant stationary solution (we eventually observe pattern formations in the network structure), as depicted in (d).
In conclusion, the mathematical study of chemotaxis models not only enhances our foundational knowledge about cells and organisms behavior in response to chemotaxis but also yields practical insights applicable in fields like medicine, bioengineering, and ecology. Looking ahead, our focus will be on how chemotaxis sensitivity influences the spreading speed of cells and organisms within network structures.
Acknowledgement: This material is based upon work supported by the National Science Foundation under Grant Number (NSF DMS-2243027), the Office of Vice-President for Research & Economic Development (OVPRED) through Research Support Program grant, and the Simons Collaboration Grant MP-TSM-00002897. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the funding organizations.