Travelling wave analysis of cellular invasion into surrounding tissues

Abstract

Single-species reaction-diffusion equations, such as the Fisher-KPP and Porous-Fisher equations, support travelling wave solutions that are often interpreted as simple mathematical models of biological invasion. Such travelling wave solutions are thought to play a role in various applications including development, wound healing and malignant invasion. One criticism of these single-species equations is that they do not explicitly describe interactions between the invading population and the surrounding environment. In this work we study a reaction-diffusion equation that describes malignant invasion which has been used to interpret experimental measurements describing the invasion of malignant melanoma cells into surrounding human skin tissues [1]. This model explicitly describes how the population of cancer cells degrade the surrounding tissues, thereby creating free space into which the cancer cells migrate and proliferate to form an invasion wave of malignant tissue that is coupled to a retreating wave of skin tissue. We analyse travelling wave solutions of this model using a combination of numerical simulation, phase plane analysis and perturbation techniques. Our analysis shows that the travelling wave solutions involve a range of very interesting properties that resemble certain well-established features of both the Fisher-KPP and Porous-Fisher equations, as well as a range of novel properties that can be thought of as extensions of these well-studied single-species equations. Of particular interest is that travelling wave solutions of the invasion model are very well approximated by trajectories in the Fisher-KPP phase plane that are normally disregarded. This observation establishes a previously unnoticed link between coupled multi-species reaction diffusion models of invasion and a different class of models of invasion that involve moving boundary problems. Matlab software to implement all calculations is available at GitHub.

Travelling-wave analysis of a model of tumour invasion with degenerate, cross-dependent diffusion

“…In particular, the existence of TWS for the Gatenby-Gawlinski model has been largely supported by a combination of numerical and analytical results [4,9,13,21,22,31]. This also holds for a simplified model of invasion by Browning et al [1,8]. However, key results were recently proved by Gallay and Mascia [11] for a reduced version of the Gatenby-Gawlinski model: they showed the existence of a form of weak TWS for any positive wave speed, c > 0.…”