Mathematical models for cell-matrix interactions during dermal wound healing

This paper contains a review of our recent work on the mathematical modeling of cell interaction with extracellular matrix components during the process of dermal wound healing. The models are of partial differential equation type and allow us to investigate in detail how various mechanochemical effects may be responsible for certain wound healing disorders such as fibrocontractive and fibroproliferative diseases. We also present a model for wound healing angiogenesis. The latter has several features in common with angiogenesis during cancer tumour growth and spread so a deeper understanding of the phenomenon in the context of wound healing may also help in the treatment of certain cancers

Macroscopic limits of individual-based models for motile cell populations with volume exclusion

Partial differential equation models are ubiquitous in studies of motile cell populations, giving a phenomenological description of events which can be analyzed and simulated using a wide range of existing tools. However, these models are seldom derived from individual cell behaviors and so it is difficult to accurately include biological hypotheses on this spatial scale. Moreover, studies which do attempt to link individual- and population-level behavior generally employ lattice-based frameworks in which the artifacts of lattice choice at the population level are unclear. In this work we derive limiting population-level descriptions of a motile cell population from an off-lattice, individual-based model (IBM) and investigate the effects of volume exclusion on the population-level dynamics. While motility with excluded volume in on-lattice IBMs can be accurately described by Fickian diffusion, we demonstrate that this is not the case off lattice. We show that the balance between two key parameters in the IBM (the distance moved in one step and the radius of an individual) determines whether volume exclusion results in enhanced or slowed diffusion. The magnitude of this effect is shown to increase with the number of cells and the rate of their movement. The method we describe is extendable to higher-dimensional and more complex systems and thereby provides a framework for deriving biologically realistic, continuum descriptions of motile populations

Diffusion-aggregation processes with mono-stable reaction terms

This paper analyses front propagation of the equation $u_\tau=[D(u)v_x]_x +f(v) \;\;\; \tau < 0, x \in \mathbb{R}$ where $f$ is a monostable (ie Fisher-type) nonlinear reaction term and $D(v)$ changes its sign once, from positive to negative values,in the interval $v \in[0,1]$ where the process is studied. This model equation accounts for simultaneous diffusive and aggregative behaviors of a population dynamic depending on the population density $v$ at time $\tau$ and position $x$. The existence of infinitely many travelling wave solutions is proven. These fronts are parametrized by their wave speed and monotonically connect the stationary states u = 0 and v = 1. In the degenerate case, i.e. when D(0) and/or D(1) = 0, sharp profiles appear, corresponding to the minimum wave speed. They also have new behaviors, in addition to those already observed in diffusive models, since they can be right compactly supported, left compactly supported, or both. The dynamics can exhibit, respectively, the phenomena of finite speed of propagation, finite speed of saturation, or both

Pattern formation in reaction diffusion models with spatially inhomogeneous diffusion coefficients

Reaction-diffusion models for biological pattern formation have been studied extensively in a variety of embryonic and ecological contexts. However, despite experimental evidence pointing to the existence of spatial inhomogeneities in various biological systems, most models have only been considered in a spatially homogeneous environment. The authors consider a two-chemical reaction-diffusion mechanism in one space dimension in which one of the diffusion coefficients depends explicitly on the spatial variable. The model is analysed in the case of a step function diffusion coefficient and the insight gained for this special case is used to discuss pattern generation for smoothly varying diffusion coefficients. The results show that spatial inhomogeneity may be an important biological pattern regulator, and possible applications of the model to chondrogenesis in the vertebrate limb are suggested

Unravelling the Turing bifurcation using spatially varying diffusion coefficients

The Turing bifurcation is the basic bifurcation generating spatial pattern, and lies at the heart of almost all mathematical models for patterning in biology and chemistry. In this paper the authors determine the structure of this bifurcation for two coupled reaction diffusion equations on a two-dimensional square spatial domain when the diffusion coefficients have a small explicit variation in space across the domain. In the case of homogeneous diffusivities, the Turing bifurcation is highly degenerate. Using a two variable perturbation method, the authors show that the small explicit spatial inhomogeneity splits the bifurcation into two separate primary and two separate secondary bifurcations, with all solution branches distinct. This splitting of the bifurcation is more effective than that given by making the domain slightly rectangular, and shows clearly the structure of the Turing bifurcation and the way in which the! var ious solution branches collapse together as the spatial variation is reduced. The authors determine the stability of the solution branches, which indicates that several new phenomena are introduced by the spatial variation, including stable subcritical striped patterns, and the possibility that stable stripes lose stability supercritically to give stable spotted patterns

Aggregative movement and front propagation for bi-stable population models

Front propagation for the aggregation-diffusion-reaction equation is investigated, where f is a bi-stable reaction-term and D(v) is a diffusion coefficient with changing sign, modeling aggregating-diffusing processes. We provide necessary and sufficient conditions for the existence of traveling wave solutions and classify them according to how or if they attain their equilibria at finite times. We also show that the dynamics can exhibit the phenomena of finite speed of propagation and/or finite speed of saturation

A mathematical model for the capillary endothelial cell-extracellular matrix interactions in wound-healing angiogenesis

Angiogenesis, the process by which new blood capillaries grow into a tissue from surrounding parent vessels, is a key event in dermal wound healing, malignant-tumour growth, and other pathologic conditions. In wound healing, new capillaries deliver vital metabolites such as amino acids and oxygen to the cells in the wound which are involved in a complex sequence of repair processes. The key cellular constituents of these new capillaries are endothelial cells: their interactions with soluble biochemical and insoluble extracellular matrix (ECM) proteins have been well documented recently, although the biological mechanisms underlying wound-healing angiogenesis are incompletely understood. Considerable recent research, including some continuum mathematical models, have focused on the interactions between endothelial cells and soluble regulators (such as growth factors). In this work, a similar modelling framework is used to investigate the roles of the insoluble ECM substrate, of which collagen is the predominant macromolecular protein. Our model consists of a partial differential equation for the endothelial-cell density (as a function of position and time) coupled to an ordinary differential equation for the ECM density. The ECM is assumed to regulate cell movement (both random and directed) and proliferation, whereas the cells synthesize and degrade the ECM. Analysis and numerical solutions of these equations highlights the roles of these processes in wound-healing angiogenesis. A nonstandard approximation analysis yields insight into the travel ling-wave structure of the system. The model is extended to two spatial dimensions (parallel and perpendicular to the plane of the skin), for which numerical simulations are presented. The model predicts that ECM-mediated random motility and cell proliferation are key processes which drive angiogenesis and that the details of the functional dependence of these processes on the ECM density, together with the rate of ECM remodelling, determine the qualitative nature of the angiogenic response. These predictions are experimentally testable, and they may lead towards a greater understanding of the biological mechanisms involved in wound-healing angiogenesis

Travelling waves in wound healing

We illustrate the role of travelling waves in wound healing by considering three different cases. Firstly, we review a model for surface wound healing in the cornea and focus on the speed of healing as a function of the application of growth factors. Secondly, we present a model for scar tissue formation in deep wounds and focus on the role of key chemicals in determining the quality of healing. Thirdly, we propose a model for excessive healing disorders and investigate how abnormal healing may be controlled

Rhythmic firing patterns in SCN: The role of circuit interactions

The suprachiasmatic nucleus (SCN) is believed to contain the main generator of circadian rhythmicity in mammals. In order to obtain further functional details of this, electrophysiological extracellular measurements in vitro were made. By means of an interspike interval distribution analysis, it is shown that there is a novel kind of neuronal firing pattern: the harmonic pattern. From these observations, we have developed a theoretical model based on possible filtering processes occurring during synaptic transmission. The model suffices to infer that regular ultradian oscillators could be an emergent property of circuit interactions of cells in the suprachiasmatic nucleus

Modeling chemotaxis reveals the role of reversed phosphotransfer and a bi-functional kinase-phosphatase

Understanding how multiple signals are integrated in living cells to produce a balanced response is a major challenge in biology. Two-component signal transduction pathways, such as bacterial chemotaxis, comprise histidine protein kinases (HPKs) and response regulators (RRs). These are used to sense and respond to changes in the environment. Rhodobacter sphaeroides has a complex chemosensory network with two signaling clusters, each containing a HPK, CheA. Here we demonstrate, using a mathematical model, how the outputs of the two signaling clusters may be integrated. We use our mathematical model supported by experimental data to predict that: (1) the main RR controlling flagellar rotation, CheY6, aided by its specific phosphatase, the bifunctional kinase CheA3, acts as a phosphate sink for the other RRs; and (2) a phosphorelay pathway involving CheB2 connects the cytoplasmic cluster kinase CheA3 with the polar localised kinase CheA2, and allows CheA3-P to phosphorylate non-cognate chemotaxis RRs. These two mechanisms enable the bifunctional kinase/phosphatase activity of CheA3 to integrate and tune the sensory output of each signaling cluster to produce a balanced response. The signal integration mechanisms identified here may be widely used by other bacteria, since like R. sphaeroides, over 50% of chemotactic bacteria have multiple cheA homologues and need to integrate signals from different sources