Travelling Waves in Hyperbolic Chemotaxis Equations

A Mathematical Modelling Study of Chemotactic Dynamics in Cell Cultures: The Impact of Spatio-temporal Heterogeneity

As motivated by studies of cellular motility driven by spatiotemporal chemotactic gradients in microdevices, we develop a framework for constructing approximate analytical solutions for the location, speed and cellular densities for cell chemotaxis waves in heterogeneous fields of chemoattractant from the underlying partial differential equation models. In particular, such chemotactic waves are not in general translationally invariant travelling waves, but possess a spatial variation that evolves in time, and may even oscillate back and forth in time, according to the details of the chemotactic gradients. The analytical framework exploits the observation that unbiased cellular diffusive flux is typically small compared to chemotactic fluxes and is first developed and validated for a range of exemplar scenarios. The framework is subsequently applied to more complex models considering the chemoattractant dynamics under more general settings, potentially including those of relevance for representing pathophysiology scenarios in microdevice studies. In particular, even though solutions cannot be constructed in all cases, a wide variety of scenarios can be considered analytically, firstly providing global insight into the important mechanisms and features of cell motility in complex spatiotemporal fields of chemoattractant. Such analytical solutions also provide a means of rapid evaluation of model predictions, with the prospect of application in computationally demanding investigations relating theoretical models and experimental observation, such as Bayesian parameter estimation.

Mathematical modeling of chemotaxis guided amoeboid cell swimming

Cells and microorganisms adopt various strategies to migrate in response to different environmental stimuli. To date, many modeling research has focused on the crawling-basedDictyostelium discoideum(Dd) cells migration induced by chemotaxis, yet recent experimental results reveal that even without adhesion or contact to a substrate, Dd cells can still swim to follow chemoattractant signals. In this paper, we develop a modeling framework to investigate the chemotaxis induced amoeboid cell swimming dynamics. A minimal swimming system consists of one deformable Dd amoeboid cell and a dilute suspension of bacteria, and the bacteria produce chemoattractant signals that attract the Dd cell. We use themathematical amoeba modelto generate Dd cell deformation and solve the resulting low Reynolds number flows, and use a moving mesh based finite volume method to solve the reaction-diffusion-convection equation. Using the computational model, we show that chemotaxis guides a swimming Dd cell to follow and catch bacteria, while on the other hand, bacterial rheotaxis may help the bacteria to escape from the predator Dd cell.

Structured models of cell migration incorporating molecular binding processes

The dynamic interplay between collective cell movement and the various molecules involved in the accompanying cell signalling mechanisms plays a crucial role in many biological processes including normal tissue development and pathological scenarios such as wound healing and cancer. Information about the various structures embedded within these processes allows a detailed exploration of the binding of molecular species to cell-surface receptors within the evolving cell population. In this paper we establish a general spatio-temporal-structural framework that enables the description of molecular binding to cell membranes coupled with the cell population dynamics. We first provide a general theoretical description for this approach and then illustrate it with three examples arising from cancer invasion.

Emergence of spatial patterns in a mathematical model for the co-culture dynamics of epithelial-like and mesenchymal-like cells

Accumulating evidence indicates that the interaction between epithelial and mesenchymal cells plays a pivotal role in cancer development and metastasis formation. Here we propose an integro-differential model for the co-culture dynamics of epithelial-like and mesenchymal-like cells. Our model takes into account the effects of chemotaxis, adhesive interactions between epithelial-like cells, proliferation and competition for nutrients. We present a sample of numerical results which display the emergence of spots, stripes and honeycomb patterns, depending on parameters and initial data. These simulations also suggest that epithelial-like and mesenchymal-like cells can segregate when there is little competition for nutrients. Furthermore, our computational results provide a possible explanation for how the concerted action between epithelial-cell adhesion and mesenchymal-cell spreading can precipitate the formation of ring-like structures, which resemble the fibrous capsules frequently observed in hepatic tumours.

A Discrete Velocity Kinetic Model with Food Metric: Chemotaxis Traveling Waves

We introduce a mesoscopic scale chemotaxis model for traveling wave phenomena which is induced by food metric. The organisms of this simplified kinetic model have two discrete velocity modes, [Formula: see text] and a constant tumbling rate. The main feature of the model is that the speed of organisms is constant [Formula: see text] with respect to the food metric, not the Euclidean metric. The uniqueness and the existence of the traveling wave solution of the model are obtained. Unlike the classical logarithmic model case there exist traveling waves under super-linear consumption rates and infinite population pulse-type traveling waves are obtained. Numerical simulations are also provided.

Travelling wave analysis in chemotaxis: case of starvation

In this paper we investigate the existence of travelling wave solutions for a chemotaxis model under the scenarios of zero growth and constant growth rate. We use Lie symmetry analysis to generate generalized travelling wave solutions, a wider class of solutions than that obtained from the standard ansatz. Unlike previous approaches, we allow for diffusivity and signal degradation. We study the influence of cell growth, diffusivity and signal degradation on the behaviour of the system. We apply realistic boundary conditions to explicitly provide biologically relevant solutions. Our results generalize known results.

Bacterial chemotaxis without gradient-sensing

Chemotaxis models are based on spatial or temporal gradient measurements by individual organisms. The key contribution of Keller and Segel (J Theor Biol 30:225-234, 1971a; J Theor Biol 30:235-248, 1971b) is showing that erratic measurements of individuals may result in an accurate chemotaxis phenomenon as a group. In this paper we provide another option to understand chemotactic behavior when individuals do not sense the gradient of chemical concentration by any means. We show that, if individuals increase their dispersal rate to find food when there is not enough food, an accurate chemotactic behavior may be obtained without sensing the gradient. Such a dispersal has been suggested by Cho and Kim (Bull Math Biol 75:845-870, 2013) and was called starvation driven diffusion. This model is surprisingly similar to the original Keller-Segel model. A comprehensive picture of traveling bands and fronts is provided.

Travelling waves in hybrid chemotaxis models

Hybrid models of chemotaxis combine agent-based models of cells with partial differential equation models of extracellular chemical signals. In this paper, travelling wave properties of hybrid models of bacterial chemotaxis are investigated. Bacteria are modelled using an agent-based (individual-based) approach with internal dynamics describing signal transduction. In addition to the chemotactic behaviour of the bacteria, the individual-based model also includes cell proliferation and death. Cells consume the extracellular nutrient field (chemoattractant), which is modelled using a partial differential equation. Mesoscopic and macroscopic equations representing the behaviour of the hybrid model are derived and the existence of travelling wave solutions for these models is established. It is shown that cell proliferation is necessary for the existence of non-transient (stationary) travelling waves in hybrid models. Additionally, a numerical comparison between the wave speeds of the continuum models and the hybrid models shows good agreement in the case of weak chemotaxis and qualitative agreement for the strong chemotaxis case. In the case of slow cell adaptation, we detect oscillating behaviour of the wave, which cannot be explained by mean-field approximations.

Travelling waves for a velocity-jump model of cell migration and proliferation

Cell invasion, characterised by moving fronts of cells, is an essential aspect of development, repair and disease. Typically, mathematical models of cell invasion are based on the Fisher-Kolmogorov equation. These traditional parabolic models cannot be used to represent experimental measurements of individual cell velocities within the invading population since they imply that information propagates with infinite speed. To overcome this limitation we study combined cell motility and proliferation based on a velocity-jump process where information propagates with finite speed. The model treats the total population of cells as two interacting subpopulations: a subpopulation of left-moving cells, L(x,t), and a subpopulation of right-moving cells, R(x,t). This leads to a system of hyperbolic partial differential equations that includes a turning rate, Λ⩾0, describing the rate at which individuals in the population change direction of movement. We present exact travelling wave solutions of the system of partial differential equations for the special case where Λ=0 and in the limit that Λ→∞. For intermediate turning rates, 0<Λ<∞, we analyse the travelling waves using the phase plane and we demonstrate a transition from smooth monotone travelling waves to smooth nonmonotone travelling waves as Λ decreases through a critical value Λcrit. We conclude by providing a qualitative comparison between the travelling wave solutions of our model and experimental observations of cell invasion. This comparison indicates that the small Λ limit produces results that are consistent with experimental observations.

Mathematical description of bacterial traveling pulses

The Keller-Segel system has been widely proposed as a model for bacterial waves driven by chemotactic processes. Current experiments on Escherichia coli have shown the precise structure of traveling pulses. We present here an alternative mathematical description of traveling pulses at the macroscopic scale. This modeling task is complemented with numerical simulations in accordance with the experimental observations. Our model is derived from an accurate kinetic description of the mesoscopic run-and-tumble process performed by bacteria. This can account for recent experimental observations with E. coli. Qualitative agreements include the asymmetry of the pulse and transition in the collective behaviour (clustered motion versus dispersion). In addition, we can capture quantitatively the traveling speed of the pulse as well as its characteristic length. This work opens several experimental and theoretical perspectives since coefficients at the macroscopic level are derived from considerations at the cellular scale. For instance, the particular response of a single cell to chemical cues turns out to have a strong effect on collective motion. Furthermore, the bottom-up scaling allows us to perform preliminary mathematical analysis and write efficient numerical schemes. This model is intended as a predictive tool for the investigation of bacterial collective motion.