Showing 1–19 of 19 results for author: Painter, K J

Phenotypic switching mechanisms determine the structure of cell migration into extracellular matrix under the `go-or-grow' hypothesis

Authors: Rebecca M. Crossley, Kevin J. Painter, Tommaso Lorenzi, Philip K. Maini, Ruth E. Baker

Abstract: A fundamental feature of collective cell migration is phenotypic heterogeneity which, for example, influences tumour progression and relapse. While current mathematical models often consider discrete phenotypic structuring of the cell population, in-line with the `go-or-grow' hypothesis \cite{hatzikirou2012go, stepien2018traveling}, they regularly overlook the role that the environment may play in… ▽ More A fundamental feature of collective cell migration is phenotypic heterogeneity which, for example, influences tumour progression and relapse. While current mathematical models often consider discrete phenotypic structuring of the cell population, in-line with the `go-or-grow' hypothesis \cite{hatzikirou2012go, stepien2018traveling}, they regularly overlook the role that the environment may play in determining the cells' phenotype during migration. Comparing a previously studied volume-filling model for a homogeneous population of generalist cells that can proliferate, move and degrade extracellular matrix (ECM) \cite{crossley2023travelling} to a novel model for a heterogeneous population comprising two distinct sub-populations of specialist cells that can either move and degrade ECM or proliferate, this study explores how different hypothetical phenotypic switching mechanisms affect the speed and structure of the invading cell populations. Through a continuum model derived from its individual-based counterpart, insights into the influence of the ECM and the impact of phenotypic switching on migrating cell populations emerge. Notably, specialist cell populations that cannot switch phenotype show reduced invasiveness compared to generalist cell populations, while implementing different forms of switching significantly alters the structure of migrating cell fronts. This key result suggests that the structure of an invading cell population could be used to infer the underlying mechanisms governing phenotypic switching. △ Less

Submitted 10 June, 2024; v1 submitted 14 January, 2024; originally announced January 2024.

Comments: 35 pages, 12 figures

Derivation and travelling wave analysis of phenotype-structured haptotaxis models of cancer invasion

Authors: Tommaso Lorenzi, Fiona R. Macfarlane, Kevin J. Painter

Abstract: We formulate haptotaxis models of cancer invasion wherein the infiltrating cancer cells can occupy a spectrum of states in phenotype space, ranging from `fully mesenchymal' to `fully epithelial'. The more mesenchymal cells are those that display stronger haptotaxis responses and have greater capacity to modify the extracellular matrix (ECM) through enhanced secretion of matrix-degrading enzymes (M… ▽ More We formulate haptotaxis models of cancer invasion wherein the infiltrating cancer cells can occupy a spectrum of states in phenotype space, ranging from `fully mesenchymal' to `fully epithelial'. The more mesenchymal cells are those that display stronger haptotaxis responses and have greater capacity to modify the extracellular matrix (ECM) through enhanced secretion of matrix-degrading enzymes (MDEs). However, as a trade-off, they have lower proliferative capacity than the more epithelial cells. The framework is multiscale in that we start with an individual-based model that tracks the dynamics of single cells, which is based on a branching random walk over a lattice representing both physical and phenotype space. We formally derive the corresponding continuum model, which takes the form of a coupled system comprising a partial integro-differential equation for the local cell population density function, a partial differential equation for the MDE concentration and an infinite-dimensional ordinary differential equation for the ECM density. Despite the intricacy of the model, we show, through formal asymptotic techniques, that for certain parameter regimes it is possible to carry out a detailed travelling wave analysis and obtain invading fronts with spatial structuring of phenotypes. Precisely, the most mesenchymal cells dominate the leading edge of the invasion wave and the most epithelial (and most proliferative) dominate the rear, representing a bulk tumour population. As such, the model recapitulates similar observations into a front to back structuring of invasion waves into leader-type and follower-type cells, witnessed in an increasing number of experimental studies over recent years. △ Less

Submitted 29 January, 2024; v1 submitted 30 October, 2023; originally announced October 2023.

A biased random walk approach for modeling the collective chemotaxis of neural crest cells

Authors: Viktoria Freingruber, Kevin J. Painter, Mariya Ptashnyk, Linus Schumacher

Abstract: Collective cell migration is a multicellular phenomenon that arises in various biological contexts, including cancer and embryo development. "Collectiveness" can be promoted by cell-cell interactions such as co-attraction and contact inhibition of locomotion. These mechanisms act on cell polarity, pivotal for directed cell motility, through influencing the intracellular dynamics of small GTPases s… ▽ More Collective cell migration is a multicellular phenomenon that arises in various biological contexts, including cancer and embryo development. "Collectiveness" can be promoted by cell-cell interactions such as co-attraction and contact inhibition of locomotion. These mechanisms act on cell polarity, pivotal for directed cell motility, through influencing the intracellular dynamics of small GTPases such as Rac1. To model these dynamics we introduce a biased random walk model, where the bias depends on the internal state of Rac1, and the Rac1 state is influenced by cell-cell interactions and chemoattractive cues. In an extensive simulation study we demonstrate and explain the scope and applicability of the introduced model in various scenarios. The use of a biased random walk model allows for the derivation of a corresponding partial differential equation for the cell density while still maintaining a certain level of intracellular detail from the individual based setting. △ Less

Submitted 2 October, 2023; originally announced October 2023.

Distinguishing between long-transient and asymptotic states in a biological aggregation model

Authors: Jonathan R. Potts, Kevin J. Painter

Abstract: Aggregations are emergent features common to many biological systems. Mathematical models to understand their emergence are consequently widespread, with the aggregation-diffusion equation being a prime example. Here we study the aggregation-diffusion equation with linear diffusion. This equation is known to support solutions that involve both single and multiple aggregations. However, numerical e… ▽ More Aggregations are emergent features common to many biological systems. Mathematical models to understand their emergence are consequently widespread, with the aggregation-diffusion equation being a prime example. Here we study the aggregation-diffusion equation with linear diffusion. This equation is known to support solutions that involve both single and multiple aggregations. However, numerical evidence suggests that the latter, which we term `multi-peaked solutions' may often be long-transient solutions rather than asymptotic steady states. We develop a novel technique for distinguishing between long transients and asymptotic steady states via an energy minimisation approach. The technique involves first approximating our study equation using a limiting process and a moment closure procedure. We then analyse local minimum energy states of this approximate system, hypothesising that these will correspond to asymptotic patterns in the aggregation-diffusion equation. Finally, we verify our hypotheses through numerical investigation, showing that our approximate analytic technique gives good predictions as to whether a state is asymptotic or a long transient. Overall, we find that almost all twin-peaked, and by extension multi-peaked, solutions are transient, except for some very special cases. We demonstrate numerically that these transients can be arbitrarily long-lived, depending on the parameters of the system. △ Less

Submitted 27 September, 2023; originally announced September 2023.

Biological Modelling with Nonlocal Advection Diffusion Equations

Authors: Kevin J Painter, Thomas Hillen, Jonathan R Potts

Abstract: The employment of nonlocal PDE models to describe biological aggregation and other phenomena has gained considerable traction in recent years. For cell populations, these methods grant a means of accommodating essential elements such as cell adhesion, critical to the development and structure of tissues. For animals, they can be used to describe how the nearby presence of conspecifics and/or heter… ▽ More The employment of nonlocal PDE models to describe biological aggregation and other phenomena has gained considerable traction in recent years. For cell populations, these methods grant a means of accommodating essential elements such as cell adhesion, critical to the development and structure of tissues. For animals, they can be used to describe how the nearby presence of conspecifics and/or heterospecifics influence movement behaviour. In this review, we will focus on classes of biological movement models in which the advective (or directed) component to motion is governed by an integral term that accounts for how the surrounding distribution(s) of the population(s) impact on a member's movement. We recount the fundamental motivation for these models: the intrinsic capacity of cell populations to self-organise and spatially sort within tissues; the wide-ranging tendency of animals towards spatial structuring, from the formations of herds and swarms to territorial segregation. We examine the derivation of these models from an individual level, illustrating in the process methods that allow models to be connected to data. We explore a growing analytical literature, including methods of stability and bifurcation analysis, and existence results. We conclude with a short section that lays out some future challenges and connections to the modelling of sociological phenomena including opinion dynamics. △ Less

Submitted 26 July, 2023; originally announced July 2023.

MSC Class: 92B05; 92C17; 92D40; 92C15

On the study of slow-fast dynamics, when the fast process has multiple invariant measures

Authors: B. D. Goddard, M. Ottobre, K. J. Painter, I. Souttar

Abstract: Motivated by applications to mathematical biology, we study the averaging problem for slow-fast systems, {\em in the case in which the fast dynamics is a stochastic process with multiple invariant measures}. We consider both the case in which the fast process is decoupled from the slow process and the case in which the two components are fully coupled. We work in the setting in which the slow proc… ▽ More Motivated by applications to mathematical biology, we study the averaging problem for slow-fast systems, {\em in the case in which the fast dynamics is a stochastic process with multiple invariant measures}. We consider both the case in which the fast process is decoupled from the slow process and the case in which the two components are fully coupled. We work in the setting in which the slow process evolves according to an Ordinary Differential Equation (ODE) and the fast process is a continuous time Markov Process with finite state space and show that, in this setting, the limiting (averaged) dynamics can be described as a random ODE (that is, an ODE with random coefficients.) Keywords. Multiscale methods, Processes with multiple equilibria, Averaging, Collective Navigation, Interacting Piecewise Deterministic Markov Processes. △ Less

Submitted 15 August, 2023; v1 submitted 8 May, 2023; originally announced May 2023.

Comments: 24 pages

MSC Class: 34C29; 34D05; 34E10; 34F05; 37H10; 60K35

Phenotype switching in chemotaxis aggregation models controls the spontaneous emergence of large densities

Authors: Kevin J Painter, Michael Winkler

Abstract: We consider a phenotype-switching chemotaxis model for aggregation, in which a chemotactic population is capable of switching back and forth between a chemotaxing state (performing chemotactic movement) and a secreting state (producing the attractant). We show that the switching rate provides a powerful mechanism for controlling the densities of spontaneously emerging aggregates. Specifically, in… ▽ More We consider a phenotype-switching chemotaxis model for aggregation, in which a chemotactic population is capable of switching back and forth between a chemotaxing state (performing chemotactic movement) and a secreting state (producing the attractant). We show that the switching rate provides a powerful mechanism for controlling the densities of spontaneously emerging aggregates. Specifically, in two- and three-dimensional settings it is shown that when both switching rates coincide and are suitably large, then the densities of both the chemotaxing and the secreting population will exceed any prescribed level at some points in the considered domain. This is complemented by two results asserting the absence of such aggregation phenomena in corresponding scenarios in which one of the switching rates remains within some bounded interval. △ Less

Submitted 10 October, 2022; originally announced October 2022.

MSC Class: 35B44 (primary); 35B36; 35B45; 35K55; 92C17 (secondary)

The impact of phenotypic heterogeneity on chemotactic self-organisation

Authors: Fiona R Macfarlane, Tommaso Lorenzi, Kevin J Painter

Abstract: The capacity to aggregate through chemosensitive movement forms a paradigm of self-organisation, with examples spanning cellular and animal systems. A basic mechanism assumes a phenotypically homogeneous population that secretes its own attractant, with the well known system introduced more than five decades ago by Keller and Segel proving resolutely popular in modelling studies. The typical assum… ▽ More The capacity to aggregate through chemosensitive movement forms a paradigm of self-organisation, with examples spanning cellular and animal systems. A basic mechanism assumes a phenotypically homogeneous population that secretes its own attractant, with the well known system introduced more than five decades ago by Keller and Segel proving resolutely popular in modelling studies. The typical assumption of population phenotypic homogeneity, however, often lies at odds with the heterogeneity of natural systems, where populations may comprise distinct phenotypes that vary according to their chemotactic ability, attractant secretion, {\it etc}. To initiate an understanding into how this diversity can impact on autoaggregation, we propose a simple extension to the classical Keller and Segel model, in which the population is divided into two distinct phenotypes: those performing chemotaxis and those producing attractant. Using a combination of linear stability analysis and numerical simulations, we demonstrate that switching between these phenotypic states alters the capacity of a population to self-aggregate. Further, we show that switching based on the local environment (population density or chemoattractant level) leads to diverse patterning and provides a route through which a population can effectively curb the size and density of an aggregate. We discuss the results in the context of real world examples of chemotactic aggregation, as well as theoretical aspects of the model such as global existence and blow-up of solutions. △ Less

Submitted 23 September, 2022; v1 submitted 29 June, 2022; originally announced June 2022.

Non-mean-field Vicsek-type models for collective behaviour

Authors: P. Buttà, B. Goddard, T. M. Hodgson, M. Ottobre, K. J. Painter

Abstract: We consider interacting particle dynamics with Vicsek type interactions, and their macroscopic PDE limit, in the non-mean-field regime; that is, we consider the case in which each particle/agent in the system interacts only with a prescribed subset of the particles in the system (for example, those within a certain distance). In this non-mean-field regime the influence between agents (i.e. the int… ▽ More We consider interacting particle dynamics with Vicsek type interactions, and their macroscopic PDE limit, in the non-mean-field regime; that is, we consider the case in which each particle/agent in the system interacts only with a prescribed subset of the particles in the system (for example, those within a certain distance). In this non-mean-field regime the influence between agents (i.e. the interaction term) can be normalised either by the total number of agents in the system (\textit{global scaling}) or by the number of agents with which the particle is effectively interacting (\textit{local scaling}). We compare the behaviour of the globally scaled and the locally scaled systems in many respects, considering for each scaling both the PDE and the corresponding particle model. In particular we observe that both the locally and globally scaled particle system exhibit pattern formation (i.e. formation of travelling-wave-like solutions) within certain parameter regimes, and generally display similar dynamics. The same is not true of the corresponding PDE models. Indeed, while both PDE models have multiple stationary states, for the globally scaled PDE such (space-homogeneous) equilibria are unstable for certain parameter regimes, with the instability leading to travelling wave solutions, while they are always stable for the locally scaled one, which never produces travelling waves. This observation is based on a careful numerical study of the model, supported by further analysis. △ Less

Submitted 17 May, 2022; originally announced June 2022.

Comments: 44 pages

Trade-offs between chemotaxis and proliferation shape the phenotypic structuring of invading waves

Authors: Tommaso Lorenzi, Kevin J Painter

Abstract: Chemotaxis-driven invasions have been proposed across a broad spectrum of biological processes, from cancer to ecology. The influential system of equations introduced by Keller and Segel has proven a popular choice in the modelling of such phenomena, but in its original form restricts to a homogeneous population. To account for the possibility of phenotypic heterogeneity, we extend to the case of… ▽ More Chemotaxis-driven invasions have been proposed across a broad spectrum of biological processes, from cancer to ecology. The influential system of equations introduced by Keller and Segel has proven a popular choice in the modelling of such phenomena, but in its original form restricts to a homogeneous population. To account for the possibility of phenotypic heterogeneity, we extend to the case of a population continuously structured across space, time and phenotype, where the latter determines variation in chemotactic responsiveness, proliferation rate, and the level of chemical environment modulation. The extended model considered here comprises a non-local partial differential equation for the local phenotype distribution of cells which is coupled, through an integral term, with a differential equation for the concentration of an attractant, which is sensed and degraded by the cells. In the framework of this model, we concentrate on a chemotaxis/proliferation trade-off scenario, where the cell phenotypes span a spectrum of states from highly-chemotactic but minimally-proliferative to minimally-chemotactic but highly-proliferative. Using a combination of numerical simulation and formal asymptotic analysis, we explore the properties of travelling-wave solutions. The results of our study demonstrate how incorporating phenotypic heterogeneity may lead to a highly-structured wave profile, where cells in different phenotypic states dominate different spatial positions across the invading wave, and clarify how the phenotypic structuring of the wave can be shaped by trade-offs between chemotaxis and proliferation. △ Less

Submitted 25 November, 2021; v1 submitted 8 September, 2021; originally announced September 2021.

Comments: 35

MSC Class: 92C17; 35Q92; 35R09; 35B40; 35C07

Leadership through influence: what mechanisms allow leaders to steer a swarm?

Authors: Sara Bernardi, Raluca Eftimie, Kevin J. Painter

Abstract: Collective migration of cells and animals often relies on a specialised set of "leaders", whose role is to steer a population of naive followers towards some target. We formulate a continuous model to understand the dynamics and structure of such groups, splitting a population into separate follower and leader types with distinct orientation responses. We incorporate "leader influence" via three p… ▽ More Collective migration of cells and animals often relies on a specialised set of "leaders", whose role is to steer a population of naive followers towards some target. We formulate a continuous model to understand the dynamics and structure of such groups, splitting a population into separate follower and leader types with distinct orientation responses. We incorporate "leader influence" via three principal mechanisms: a bias in the orientation of leaders according to the destination, distinct speeds of movement and distinct levels of conspicuousness. Using a combination of analysis and numerical computation on a sequence of models of increasing complexity, we assess the extent to which leaders successfully shepherd the swarm. While all three mechanisms can lead to a successfully steered swarm, parameter regime is crucial with non successful choices generating a variety of unsuccessful attempts, including movement away from the target, swarm splitting or swarm dispersal. △ Less

Submitted 9 February, 2021; originally announced February 2021.

Direction-Dependent Turning Leads to Anisotropic Diffusion and Persistence

Authors: Nadia Loy, Thomas Hillen, Kevin John Painter

Abstract: Cells and organisms follow aligned structures in their environment, a process that can generate persistent migration paths. Kinetic transport equations are a popular modelling tool for describing biological movements at the mesoscopic level, yet their formulations usually assume a constant turning rate. Here we relax this simplification, extending to include a turning rate that varies according to… ▽ More Cells and organisms follow aligned structures in their environment, a process that can generate persistent migration paths. Kinetic transport equations are a popular modelling tool for describing biological movements at the mesoscopic level, yet their formulations usually assume a constant turning rate. Here we relax this simplification, extending to include a turning rate that varies according to the anisotropy of a heterogeneous environment. We extend known methods of parabolic and hyperbolic scaling and apply the results to cell movement on micro-patterned domains. We show that inclusion of orientation dependence in the turning rate can lead to persistence of motion in an otherwise fully symmetric environment, and generate enhanced diffusion in structured domains. △ Less

Submitted 12 January, 2021; originally announced January 2021.

Macroscopic descriptions of follower-leader systems

Authors: Sara Bernardi, Gissell Estrada-Rodriguez, Heiko Gimperlein, Kevin J. Painter

Abstract: The fundamental derivation of macroscopic model equations to describe swarms based on microscopic movement laws and mathematical analyses into their self-organisation capabilities remains a challenge from the perspective of both modelling and analysis. In this paper we clarify relevant continuous macroscopic model equations that describe follower-leader interactions for a swarm where these two pop… ▽ More The fundamental derivation of macroscopic model equations to describe swarms based on microscopic movement laws and mathematical analyses into their self-organisation capabilities remains a challenge from the perspective of both modelling and analysis. In this paper we clarify relevant continuous macroscopic model equations that describe follower-leader interactions for a swarm where these two populations are fixed. We study the behaviour of the swarm over long and short time scales to shed light on the number of leaders needed to initiate swarm movement, according to the homogeneous or inhomogeneous nature of the interaction (alignment) kernel. The results indicate the crucial role played by the interaction kernel to model transient behaviour. △ Less

Submitted 18 October, 2021; v1 submitted 1 August, 2019; originally announced August 2019.

Comments: 22 pages, accepted for publication in Kinetic and Related Models

Journal ref: Kinetic and Related Models 14 (2021), 981-1002

Mathematical models for chemotaxis and their applications in self-organisation phenomena

Authors: K. J. Painter

Abstract: Chemotaxis is a fundamental guidance mechanism of cells and organisms, responsible for attracting microbes to food, embryonic cells into developing tissues, immune cells to infection sites, animals towards potential mates, and mathematicians into biology. The Patlak-Keller-Segel (PKS) system forms part of the bedrock of mathematical biology, a go-to-choice for modellers and analysts alike. For the… ▽ More Chemotaxis is a fundamental guidance mechanism of cells and organisms, responsible for attracting microbes to food, embryonic cells into developing tissues, immune cells to infection sites, animals towards potential mates, and mathematicians into biology. The Patlak-Keller-Segel (PKS) system forms part of the bedrock of mathematical biology, a go-to-choice for modellers and analysts alike. For the former it is simple yet recapitulates numerous phenomena; the latter are attracted to these rich dynamics. Here I review the adoption of PKS systems when explaining self-organisation processes. I consider their foundation, returning to the initial efforts of Patlak and Keller and Segel, and briefly describe their patterning properties. Applications of PKS systems are considered in their diverse areas, including microbiology, development, immunology, cancer, ecology and crime. In each case a historical perspective is provided on the evidence for chemotactic behaviour, followed by a review of modelling efforts; a compendium of the models is included as an Appendix. Finally, a half-serious/half-tongue-in-cheek model is developed to explain how cliques form in academia. Assumptions in which scholars alter their research line according to available problems leads to clustering of academics and the formation of "hot" research topics. △ Less

Submitted 25 June, 2018; v1 submitted 22 June, 2018; originally announced June 2018.

Comments: 35 pages, 8 figures, Submitted to Journal of Theoretical Biology

MSC Class: 92C17

Space-time fractional diffusion in cell movement models with delay

Authors: Gissell Estrada-Rodriguez, Heiko Gimperlein, Kevin J. Painter, Jakub Stocek

Abstract: The movement of organisms and cells can be governed by occasional long distance runs, according to an approximate Lévy walk. For T cells migrating through chronically-infected brain tissue, runs are further interrupted by long pauses, and the aim here is to clarify the form of continuous model equations which describe such movements. Starting from a microscopic velocity-jump model based on experim… ▽ More The movement of organisms and cells can be governed by occasional long distance runs, according to an approximate Lévy walk. For T cells migrating through chronically-infected brain tissue, runs are further interrupted by long pauses, and the aim here is to clarify the form of continuous model equations which describe such movements. Starting from a microscopic velocity-jump model based on experimental observations, we include power-law distributions of run and waiting times and investigate the relevant parabolic limit from a kinetic equation for resting and moving individuals. In biologically relevant regimes we derive nonlocal diffusion equations, including fractional Laplacians in space and fractional time derivatives. Its analysis and numerical experiments shed light on how the searching strategy, and the impact from chemokinesis responses to chemokines, shorten the average time taken to find rare targets in the absence of direct guidance information such as chemotaxis. △ Less

Submitted 29 October, 2018; v1 submitted 23 February, 2018; originally announced February 2018.

Comments: 25 pages, 8 figures, Mathematical Models and Methods in Applied Sciences (2019)

Journal ref: Mathematical Models and Methods in Applied Sciences 29 (2019), 65-88

Fractional Patlak-Keller-Segel equations for chemotactic superdiffusion

Authors: Gissell Estrada-Rodriguez, Heiko Gimperlein, Kevin J. Painter

Abstract: The long range movement of certain organisms in the presence of a chemoattractant can be governed by long distance runs, according to an approximate Levy distribution. This article clarifies the form of biologically relevant model equations: We derive Patlak-Keller-Segel-like equations involving nonlocal, fractional Laplacians from a microscopic model for cell movement. Starting from a power-law d… ▽ More The long range movement of certain organisms in the presence of a chemoattractant can be governed by long distance runs, according to an approximate Levy distribution. This article clarifies the form of biologically relevant model equations: We derive Patlak-Keller-Segel-like equations involving nonlocal, fractional Laplacians from a microscopic model for cell movement. Starting from a power-law distribution of run times, we derive a kinetic equation in which the collision term takes into account the long range behaviour of the individuals. A fractional chemotactic equation is obtained in a biologically relevant regime. Apart from chemotaxis, our work has implications for biological diffusion in numerous processes. △ Less

Submitted 24 January, 2018; v1 submitted 9 August, 2017; originally announced August 2017.

Comments: 20 pages, 4 figures, to appear in SIAM Journal on Applied Mathematics

Journal ref: SIAM Journal on Applied Mathematics 78 (2018), 1155-1173

Spatio-temporal Models of Lymphangiogenesis in Wound Healing

Authors: Arianna Bianchi, Kevin J. Painter, Jonathan A. Sherratt

Abstract: Several studies suggest that one possible cause of impaired wound healing is failed or insufficient lymphangiogenesis, that is the formation of new lymphatic capillaries. Although many mathematical models have been developed to describe the formation of blood capillaries (angiogenesis), very few have been proposed for the regeneration of the lymphatic network. Lymphangiogenesis is a markedly diffe… ▽ More Several studies suggest that one possible cause of impaired wound healing is failed or insufficient lymphangiogenesis, that is the formation of new lymphatic capillaries. Although many mathematical models have been developed to describe the formation of blood capillaries (angiogenesis), very few have been proposed for the regeneration of the lymphatic network. Lymphangiogenesis is a markedly different process from angiogenesis, occurring at different times and in response to different chemical stimuli. Two main hypotheses have been proposed: 1) lymphatic capillaries sprout from existing interrupted ones at the edge of the wound in analogy to the blood angiogenesis case; 2) lymphatic endothelial cells first pool in the wound region following the lymph flow and then, once sufficiently populated, start to form a network. Here we present two PDE models describing lymphangiogenesis according to these two different hypotheses. Further, we include the effect of advection due to interstitial flow and lymph flow coming from open capillaries. The variables represent different cell densities and growth factor concentrations, and where possible the parameters are estimated from biological data. The models are then solved numerically and the results are compared with the available biological literature. △ Less

Submitted 15 April, 2016; originally announced April 2016.

Comments: 29 pages, 9 Figures, 6 Tables (39 figure files in total)

A Mathematical Model for Lymphangiogenesis in Normal and Diabetic Wounds

Authors: Arianna Bianchi, Kevin J. Painter, Jonathan A. Sherratt

Abstract: Several studies suggest that one possible cause of impaired wound healing is failed or insufficient lymphangiogenesis, that is the formation of new lymphatic capillaries. Although many mathematical models have been developed to describe the formation of blood capillaries (angiogenesis) very few have been proposed for the regeneration of the lymphatic network. Moreover, lymphangiogenesis is markedl… ▽ More Several studies suggest that one possible cause of impaired wound healing is failed or insufficient lymphangiogenesis, that is the formation of new lymphatic capillaries. Although many mathematical models have been developed to describe the formation of blood capillaries (angiogenesis) very few have been proposed for the regeneration of the lymphatic network. Moreover, lymphangiogenesis is markedly distinct from angiogenesis, occurring at different times and in a different manner. Here a model of five ordinary differential equations is presented to describe the formation of lymphatic capillaries following a skin wound. The variables represent different cell densities and growth factor concentrations, and where possible the parameters are estimated from experimental and clinical data. The system is then solved numerically and the results are compared with the available biological literature. Finally, a parameter sensitivity analysis of the model is taken as a starting point for suggesting new therapeutic approaches targeting the enhancement of lymphangiogenesis in diabetic wounds. The work provides a deeper understanding of the phenomenon in question, clarifying the main factors involved. In particular, the balance between TGF-$βべーた$ and VEGF levels, rather than their absolute values, is identified as crucial to effective lymphangiogenesis. In addition, the results indicate lowering the macrophage-mediated activation of TGF-$βべーた$ and increasing the basal lymphatic endothelial cell growth rate, \emph{inter alia}, as potential treatments. It is hoped the findings of this paper may be considered in the development of future experiments investigating novel lymphangiogenic therapies. △ Less

Submitted 2 January, 2015; originally announced January 2015.

MSC Class: 92B05

Travelling waves in hybrid chemotaxis models

Authors: Benjamin Franz, Chuan Xue, Kevin J. Painter, Radek Erban

Abstract: 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 b… ▽ More 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. △ Less

Submitted 11 February, 2013; originally announced February 2013.

Comments: 24 pages, submitted to Bulletin of Mathematical Biology