Turing instabilities and spatio-temporal chaos in ratio-dependent Holling-Tanner model

In this paper we consider a modified spatiotemporal ecological system originating from the temporal Holling-Tanner model, by incorporating diffusion terms. The original ODE system is studied for the stability of coexisting homogeneous steady-states. The modified PDE system is investigated in detail with both numerical and analytical approaches. Both the Turing and non-Turing patterns are examined for some fixed parametric values and some interesting results have been obtained for the prey and predator populations. Numerical simulation shows that either prey or predator population do not converge to any stationary state at any future time when parameter values are taken in the Turing-Hopf domain. Prey and predator populations exhibit spatiotemporal chaos resulting from temporal oscillation of both the population and spatial instability. With help of numerical simulations we have shown that Turing-Hopf bifurcation leads to onset of spatio-temporal chaos when predator's diffusivity is much higher compared to prey population. Our investigation reveals the fact that Hopf-bifurcation is essential for the onset of spatiotemporal chaos.     

Turing instabilities in general systems

We present necessary and sufficient conditions on the stability matrix of a general n(≥2)-dimensional reaction-diffusion system which guarantee that its uniform steady state can undergo a Turing bifurcation. The necessary (kinetic) condition, requiring that the system be composed of an unstable (or activator) and a stable (or inhibitor) subsystem, and the sufficient condition of sufficiently rapid inhibitor diffusion relative to the activator subsystem are established in three theorems which form the core of our results. Given the possibility that the unstable (activator) subsystem involves several species (dimensions), we present a classification of the analytically deduced Turing bifurcations into p (1 ≤p≤ (n− 1)) different classes. For n = 3 dimensions we illustrate numerically that two types of steady Turing pattern arise in one spatial dimension in a generic reaction-diffusion system. The results confirm the validity of an earlier conjecture [12] and they also characterise the class of so-called strongly stable matrices for which only necessary conditions have been known before [23, 24]. One of the main consequences of the present work is that biological morphogens, which have so far been expected to be single chemical species [1–9], may instead be composed of two or more interacting species forming an unstable subsystem. Received: 21 September 1999 / Revised version: 21 June 2000 / Published online: 24 November 2000     

A new necessary condition for Turing instabilities

Reactivity (a.k.a initial growth) is necessary for diffusion driven instability (Turing instability). Using a notion of common Lyapunov function we show that this necessary condition is a special case of a more powerful (i.e. tighter) necessary condition. Specifically, we show that if the linearised reaction matrix and the diffusion matrix share a common Lyapunov function, then Turing instability is not possible. The existence of common Lyapunov functions is readily checked using semi-definite programming. We apply this result to the Gierer-Meinhardt system modelling regenerative properties of Hydra, the Oregonator, to a host-parasite-hyperparasite system with diffusion and to a reaction-diffusion-chemotaxis model for a multi-species host-parasitoid community.     

Synthetic Turing protocells: vesicle self-reproduction through symmetry-breaking instabilities

The reproduction of a living cell requires a repeatable set of chemical events to be properly coordinated. Such events define a replication cycle, coupling the growth and shape change of the cell membrane with internal metabolic reactions. Although the logic of such process is determined by potentially simple physico-chemical laws, modelling of a full, self-maintained cell cycle is not trivial. Here we present a novel approach to the problem that makes use of so-called symmetry breaking instabilities as the engine of cell growth and division. It is shown that the process occurs as a consequence of the breaking of spatial symmetry and provides a reliable mechanism of vesicle growth and reproduction. Our model opens the possibility of a synthetic protocell lacking information but displaying self-reproduction under a very simple set of chemical reactions.     

Transient dynamics and pattern formation: reactivity is necessary for Turing instabilities.

The theory of spatial pattern formation via Turing bifurcations - wherein an equilibrium of a nonlinear system is asymptotically stable in the absence of dispersal but unstable in the presence of dispersal - plays an important role in biology, chemistry and physics. It is an asymptotic theory, concerned with the long-term behavior of perturbations. In contrast, the concept of reactivity describes the short-term transient behavior of perturbations to an asymptotically stable equilibrium. In this article we show that there is a connection between these two seemingly disparate concepts. In particular, we show that reactivity is necessary for Turing instability in multispecies systems of reaction-diffusion equations, integrodifference equations, coupled map lattices, and systems of ordinary differential equations.     

Propagation of Turing patterns in a plankton model

The paper is devoted to a reaction–diffusion system of equations describing phytoplankton and zooplankton distributions. Linear stability analysis of the model is carried out. Turing and Hopf stability boundaries are found. Emergence of two-dimensional spatial structures is illustrated by numerical simulations. Travelling waves between various stationary solutions are investigated. Transitions between homogeneous in space stationary solutions and Turing structures are studied.     

Propagation of Turing patterns in a plankton model

The paper is devoted to a reaction-diffusion system of equations describing phytoplankton and zooplankton distributions. Linear stability analysis of the model is carried out. Turing and Hopf stability boundaries are found. Emergence of two-dimensional spatial structures is illustrated by numerical simulations. Travelling waves between various stationary solutions are investigated. Transitions between homogeneous in space stationary solutions and Turing structures are studied.     

A mathematical model of metabolic insulin signaling pathways

We develop a mathematical model that explicitly represents many of the known signaling components mediating translocation of the insulin-responsive glucose transporter GLUT4 to gain insight into the complexities of metabolic insulin signaling pathways. A novel mechanistic model of postreceptor events including phosphorylation of insulin receptor substrate-1, activation of phosphatidylinositol 3-kinase, and subsequent activation of downstream kinases Akt and protein kinase C-zeta is coupled with previously validated subsystem models of insulin receptor binding, receptor recycling, and GLUT4 translocation. A system of differential equations is defined by the structure of the model. Rate constants and model parameters are constrained by published experimental data. Model simulations of insulin dose-response experiments agree with published experimental data and also generate expected qualitative behaviors such as sequential signal amplification and increased sensitivity of downstream components. We examined the consequences of incorporating feedback pathways as well as representing pathological conditions, such as increased levels of protein tyrosine phosphatases, to illustrate the utility of our model for exploring molecular mechanisms. We conclude that mathematical modeling of signal transduction pathways is a useful approach for gaining insight into the complexities of metabolic insulin signaling.     

Calcium instabilities in mammalian cardiomyocyte networks

The degeneration of a regular heart rhythm into fibrillation (a chaotic or chaos-like sequence) can proceed via several classical routes described by nonlinear dynamics: period-doubling, quasiperiodicity, or intermittency. In this study, we experimentally examine one aspect of cardiac excitation dynamics, the long-term evolution of intracellular calcium signals in cultured cardiomyocyte networks subjected to increasingly faster pacing rates via field stimulation. In this spatially extended system, we observed alternans and higher-order periodicities, extra beats, and skipped beats or blocks. Calcium instabilities evolved nonmonotonically with the prevalence of phase-locking or Wenckebach rhythm, low-frequency magnitude modulations (signature of quasiperiodicity), and switches between patterns with occasional bursts (signature of intermittency), but period-doubling bifurcations were rare. Six ventricular-fibrillation-resembling episodes were pace-induced, for which significantly higher complexity was confirmed by approximate entropy calculations. The progressive destabilization of the heart rhythm by coexistent frequencies, seen in this study, can be related to theoretically predicted competition of control variables (voltage and calcium) at the single-cell level, or to competition of excitation and recovery at the cell network level. Optical maps of the response revealed multiple local spatiotemporal patterns, and the emergence of longer-period global rhythms as a result of wavebreak-induced reentries.     

Investigation of early events in Fc epsilon RI-mediated signaling using a detailed mathematical model

Aggregation of Fc epsilon RI on mast cells and basophils leads to autophosphorylation and increased activity of the cytosolic protein tyrosine kinase Syk. We investigated the roles of the Src kinase Lyn, the immunoreceptor tyrosine-based activation motifs (ITAMs) on the beta and gamma subunits of Fc epsilon RI, and Syk itself in the activation of Syk. Our approach was to build a detailed mathematical model of reactions involving Fc epsilon RI, Lyn, Syk, and a bivalent ligand that aggregates Fc(epsilon)RI. We applied the model to experiments in which covalently cross-linked IgE dimers stimulate rat basophilic leukemia cells. The model makes it possible to test the consistency of mechanistic assumptions with data that alone provide limited mechanistic insight. For example, the model helps sort out mechanisms that jointly control dephosphorylation of receptor subunits. In addition, interpreted in the context of the model, experimentally observed differences between the beta- and gamma-chains with respect to levels of phosphorylation and rates of dephosphorylation indicate that most cellular Syk, but only a small fraction of Lyn, is available to interact with receptors. We also show that although the beta ITAM acts to amplify signaling in experimental systems where its role has been investigated, there are conditions under which the beta ITAM will act as an inhibitor.     

A mathematical model of intercellular signaling during epithelial wound healing

Recent experiments monitoring the healing process of wounded epithelial monolayers have demonstrated the necessity of MAPK activation for coordinated cell movement after damage. This MAPK activity is characterized by two wave-like phenomena. One MAPK “wave” that originates immediately after injury, propagates deep into the cell sheet, away from the edge, and then rebounds back to the wound interface. After this initial MAPK activity has largely disappeared, a second MAPK front propagates slowly from the wound interface and also continues into the cell sheet, maintaining a sustained level of MAPK activity throughout the cell sheet. It has been suggested that the first wave is initiated by Reactive Oxygen Species (ROS) generated at the time of injury. In this work, we develop a minimal mathematical model that reproduces the observed behavior. The main ingredients of our model are a competition between ligand (e.g., Epithelial Growth Factor) and ROS for the activation of Epithelial Growth Factor Receptor, and a feedback loop between receptor occupancy and MAPK activation. We explore the mathematical properties of the model and look for traveling wave solutions consistent with the experimentally observed MAPK activity patterns.     

On a mathematical model for body tissue inflammation

In the present paper I will try to prove the mathematical validity of a model on the localized bacterial infection for tissue inflammation. This model was proposed by Lauffenburger and Kennedy [3], and it describes the inflammatory response to bacterial invasion of body tissue. I prove the mathematical validity of the model by means of a positivity theorem, an existence theorem and a uniqueness theorem. In spite of the apparent simplicity of the problem, the solution requires a delicate set of techniques. It seems very difficult to extend these techniques to a model in more than one dimension.     

A discrete mathematical model applied to genetic regulation and metabolic networks

This paper describes the use of a discrete mathematical model to represent the basic mechanisms of regulation of the bacteria E. coli in batch fermentation. The specific phenomena studied were the changes in metabolism and genetic regulation when the bacteria use three different carbon substrates (glucose, glycerol, and acetate). The model correctly predicts the behavior of E. coli vis-à-vis substrate mixtures. In a mixture of glucose, glycerol, and acetate, it prefers glucose, then glycerol, and finally acetate. The model included 67 nodes; 28 were genes, 20 enzymes, and 19 regulators/biochemical compounds. The model represents both the genetic regulation and metabolic networks in an inrtegrated form, which is how they function biologically. This is one of the first attempts to include both of these networks in one model. Previously, discrete mathematical models were used only to describe genetic regulation networks. The study of the network dynamics generated 8 (2(3)) fixed points, one for each nutrient configuration (substrate mixture) in the medium. The fixed points of the discrete model reflect the phenotypes described. Gene expression and the patterns of the metabolic fluxes generated are described accurately. The activation of the gene regulation network depends basically on the presence of glucose and glycerol. The model predicts the behavior when mixed carbon sources are utilized as well as when there is no carbon source present. Fictitious jokers (Joker1, Joker2, and Repressor SdhC) had to be created to control 12 genes whose regulation mechanism is unknown, since glycerol and glucose do not act directly on the genes. The approach presented in this paper is particularly useful to investigate potential unknown gene regulation mechanisms; such a novel approach can also be used to describe other gene regulation situations such as the comparison between non-recombinant and recombinant yeast strain, producing recombinant proteins, presently under investigation in our group.     

Scrapie transmission in Britain: a recipe for a mathematical model

Responses to an anonymous postal survey concerning scrapie are analysed. Risk factors associated with farms that have had scrapie are identified as size, geographical region, lambing practices and holding of certain breeds. Further analysis of farms that have scrapie only in bought-in animals reveals that such farms tend to breed a smaller proportion of their replacement animals than farms without scrapie. Farms that have had scrapie in home-bred animals have attributes associated with breeding many animals: large numbers of rams bought, few ewes bought, and many animals that are home-bred. The demography of British sheep farms as described by size, breeds, purchasing behaviour, age structure and proportion of animals that are home-bred is summarized. British farms with scrapie reveal certain special features: they have more sheep that are found dead, more elderly ewes and more cases of scab.     

Analysis of unstable behavior in a mathematical model for erythropoiesis

We consider an age-structured model that describes the regulation of erythropoiesis through the negative feedback loop between erythropoietin and hemoglobin. This model is reduced to a system of two ordinary differential equations with two constant delays for which we show existence of a unique steady state. We determine all instances at which this steady state loses stability via a Hopf bifurcation through a theoretical bifurcation analysis establishing analytical expressions for the scenarios in which they arise. We show examples of supercritical Hopf bifurcations for parameter values estimated according to physiological values for humans found in the literature and present numerical simulations in agreement with the theoretical analysis. We provide a strategy for parameter estimation to match empirical measurements and predict dynamics in experimental settings, and compare existing data on hemoglobin oscillation in rabbits with predictions of our model.     

Analysis of a mathematical model for transport of I-albumin

The hypotheses and results given are motivated by the study of the distribution of albumin in man which represents a class of delay-differential systems. The approach used is to study the behavior of the solutions of nonlinear delay-differential systems with variable coefficients under the assumptions of continuity and boundedness of coefficients. The criterions are conditions on the roots of a certain “quasi-polynomial”, i.e., a polynomial in a variable and exponential of that variable. These criterions bear a resemblance to the ones in the case of constant coefficients and retardations and are applicable to this case also. The method is based on Lyapunov type functional with appropriate properties.