Are buffers boring? Uniqueness and asymptotical stability of traveling wave fronts in the buffered bistable system

Traveling waves of calcium are widely observed under the condition that the free cytosolic calcium is buffered. Thus it is of physiological interest to determine how buffers affect the properties of calcium waves. Here we summarise and extend previous results on the existence, uniqueness and stability of traveling wave solutions of the buffered bistable equation, which is the simplest possible model of the upstroke of a calcium wave. Taken together, the results show that immobile buffers do not change the existence, uniqueness or stability of the traveling wave, while mobile buffers can eliminate a traveling wave. However, if a wave exists in the latter case, it remains unique and stable.      

Travelling waves of a delayed SIR epidemic model with nonlinear incidence rate and spatial diffusion

This paper is concerned with the existence of travelling waves to a SIR epidemic model with nonlinear incidence rate, spatial diffusion and time delay. By analyzing the corresponding characteristic equations, the local stability of a disease-free steady state and an endemic steady state to this system under homogeneous Neumann boundary conditions is discussed. By using the cross iteration method and the Schauder's fixed point theorem, we reduce the existence of travelling waves to the existence of a pair of upper-lower solutions. By constructing a pair of upper-lower solutions, we derive the existence of a travelling wave connecting the disease-free steady state and the endemic steady state. Numerical simulations are carried out to illustrate the main results.     

Aggregative effects for a reaction-advection equation

The time course of a single population subject to logistic growth and drift towards regions of increasing population density is modelled by a quasilinear differential equation of the first order. The stationary solutions and the travelling waves are investigated. The existence of travelling waves with compact support is used to prove, among other properties, that populations initially concentrated in a finite region preserve this property for all future times.Visiting A. von Humboldt Fellow at the Universität Tübingen, Lehrstuhl für Biomathematik, Auf der Morgenstelle 28, D-7406 Tübingen 1     

Homoclinic solutions in mechanical systems with small dissipation. Application to DNA dynamics

 This paper analyzes the problem of persistence of homoclinic solutions to perturbed systems of second order ODE's. These systems arise from PDE's, when considering solutions in the form of travelling waves. It is shown that homoclinic solutions persist in the presence of dissipation. Dissipation can be balanced by nonautonomous terms of compact support, which are controlled by a single parameter. This result is applied to prove the existence of torsional pulse-like travelling waves propagating along a nonelastic DNA molecule. In this case the energy is added to system by advancing the RNA polymerase. Received: 17 August 2000 / Revised version: 26 October 2001 / Published online: 14 March 2002     

On calcium-buffer dynamics within the excess buffer regime

In intracellular calcium signaling, calcium buffers has been recognized for their role in reshaping and localizing the calcium concentration profile in the vicinity of the channel, as well as reducing the effective diffusion of free calcium. In the presence of an excess of endogenous or exogenous buffers, linearization of the reaction-diffusion system describing the calcium-buffer dynamics has been instrumental in understanding the extent of the microdomain formation and in quantifying the apparent diffusion of the free calcium. In these linearized models, the conclusions are usually drawn from the steady-state solutions upon the opening of the channel. In this work, using the joint Laplace-Fourier method, we give an explicit integral transient solution, as well as, the long-time asymptotic behavior of the linearized calcium-buffer dynamics. The results confirm and emphasize the long stated intuitions on the diffusive character of the calcium-buffer dynamics. Numerical validations of our analytical results will be discussed.     

Travelling wave solutions of diffusive Lotka-Volterra equations

We establish the existence of travelling wave solutions for two reaction diffusion systems based on the Lotka-Volterra model for predator and prey interactions. For simplicity, we consider only 1 space dimension. The waves are of transition front type, analogous to the travelling wave solutions discussed by Fisher and Kolmogorov et al. for a scalar reaction diffusion equation. The waves discussed here are not necessarily monotone. For any speed c there is a travelling wave solution of transition front type. For one of the systems discussed here, there is a distinguished speed c^* dividing the waves into two types, waves of speed c < c^* being one type, waves of speed c ? c^* being of the other type. We present numerical evidence that for this system the wave of speed c^* is stable, and that c^* is an asymptotic speed of propagation in some sense. For the other system, waves of all speeds are in some sense stable. The proof of existence uses a shooting argument and a Lyapunov function. We also discuss some possible biological implications of the existence of these waves.     

Two-dimensional determination of the cellular Ca2+ binding in bovine chromaffin cells.

The spatiotemporal profile of intracellular calcium signals is determined by the flux of calcium ions across different biological membranes as well as by the diffusional mobility of calcium and different calcium buffers in the cell. To arrive at a quantitative understanding of the determinants of these signals, one needs to dissociate the flux contribution from the redistribution and buffering of calcium. Since the cytosol can be heterogeneous with respect to its calcium buffering property, it is essential to assess this property in a spatially resolved manner. In this paper we report on two different methods to estimate the cellular calcium binding of bovine adrenal chromaffin cells. In the first method, we use voltage-dependent calcium channels as a source to generate calcium gradients in the cytosol. Using imaging techniques, we monitor the dissipation of these gradients to estimate local apparent calcium diffusion coefficients and, from these, local calcium binding ratios. This approach requires a very high signal-to-noise ratio of the calcium measurement and can be used when well-defined calcium gradients can be generated throughout the cell. In the second method, we overcome these problems by using calcium-loaded DM-nitrophen as a light-dependent calcium source to homogeneously and quantitatively release calcium in the cytosol. By measuring [Ca2+] directly before and after the photorelease process and knowing the total amount of calcium being released photolytically, we get an estimate of the fraction of calcium ions which does not appear as free calcium and hence must be bound to either the indicator dye or the endogenous calcium buffer. This finally results in a two-dimensional map of the distribution of the immobile endogenous calcium buffer. We did not observe significant variations of the cellular calcium binding at a spatial resolution of approximately 2 micron. Furthermore, the calcium binding is not reduced by increasing the resting [Ca2+] to levels as high as 1.1 microM. This is indicative of a low calcium affinity of the corresponding buffers and is in agreement with a recent report on the affinity of these buffers (Xu, T., M. Naraghi, H. Kang, and E. Neher. 1997. Biophys. J. 73:532-545). In contrast to the homogeneous distribution of the calcium buffers, the apparant calcium diffusion coefficient did show inhomogeneities, which can be attributed to restricted diffusion at the nuclear envelope and to rim effects at the cell membrane.     

A Simple Mathematical Model Inspired by the Purkinje Cells: From Delayed Travelling Waves to Fractional Diffusion

Recently, several experiments have demonstrated the existence of fractional diffusion in the neuronal transmission occurring in the Purkinje cells, whose malfunctioning is known to be related to the lack of voluntary coordination and the appearance of tremors. Also, a classical mathematical feature is that (fractional) parabolic equations possess smoothing effects, in contrast with the case of hyperbolic equations, which typically exhibit shocks and discontinuities. In this paper, we show how a simple toy-model of a highly ramified structure, somehow inspired by that of the Purkinje cells, may produce a fractional diffusion via the superposition of travelling waves that solve a hyperbolic equation. This could suggest that the high ramification of the Purkinje cells might have provided an evolutionary advantage of “smoothing” the transmission of signals and avoiding shock propagations (at the price of slowing a bit such transmission). Although an experimental confirmation of the possibility of such evolutionary advantage goes well beyond the goals of this paper, we think that it is intriguing, as a mathematical counterpart, to consider the time fractional diffusion as arising from the superposition of delayed travelling waves in highly ramified transmission media. The case of a travelling concave parabola with sufficiently small curvature is explicitly computed. The new link that we propose between time fractional diffusion and hyperbolic equation also provides a novelty with respect to the usual paradigm relating time fractional diffusion with parabolic equations in the limit. This paper is written in such a way as to be of interest to both biologists and mathematician alike. In order to accomplish this aim, both complete explanations of the objects considered and detailed lists of references are provided.     

Do calcium buffers always slow down the propagation of calcium waves?

Calcium buffers are large proteins that act as binding sites for free cytosolic calcium. Since a large fraction of cytosolic calcium is bound to calcium buffers, calcium waves are widely observed under the condition that free cytosolic calcium is heavily buffered. In addition, all physiological buffered excitable systems contain multiple buffers with different affinities. It is thus important to understand the properties of waves in excitable systems with the inclusion of buffers. There is an ongoing controversy about whether or not the addition of calcium buffers into the system always slows down the propagation of calcium waves. To solve this controversy, we incorporate the buffering effect into the generic excitable system, the FitzHugh–Nagumo model, to get the buffered FitzHugh–Nagumo model, and then to study the effect of the added buffer with large diffusivity on traveling waves of such a model in one spatial dimension. We can find a critical dissociation constant ( $K=K(a)$ ) characterized by system excitability parameter $a$ such that calcium buffers can be classified into two types: weak buffers ( $K\in (K(a),\infty )$ ) and strong buffers ( $K\in (0,K(a))$ ). We analytically show that the addition of weak buffers or strong buffers but with its total concentration $b_0^1$ below some critical total concentration $b_{0,c}^1$ into the system can generate a traveling wave of the resulting system which propagates faster than that of the origin system, provided that the diffusivity $D_1$ of the added buffers is sufficiently large. Further, the magnitude of the wave speed of traveling waves of the resulting system is proportional to $\sqrt{D_1}$ as $D_1\rightarrow \infty $ . In contrast, the addition of strong buffers with the total concentration $b_0^1>b_{0,c}^1$ into the system may not be able to support the formation of a biologically acceptable wave provided that the diffusivity $D_1$ of the added buffers is sufficiently large.     

A continuum model for surface contraction waves

The periodic travelling waves which appear on some animal eggs after fertilization are considered here. These are thought to be caused by a calcium initiated calcium release on the surface, causing calcium waves. A continuum model is developed where the cell is treated as a small viscous droplet with a surface contamination. When a periodic source of surfactant acts at one pole and propagates down the cell surface to the opposite pole, the drop responds by forming constriction rings which move from pole to pole.     

一类具有分布时滞的扩散种群模型行波解的存在性

研究了一类具有分布时滞的扩散种群模型行波解的存在性,证明了当平均时滞充分小时,方程具有连接两个平衡点的单调行波解.     

Travelling Waves in Hyperbolic Chemotaxis Equations

Mathematical models of bacterial populations are often written as systems of partial differential equations for the densities of bacteria and concentrations of extracellular (signal) chemicals. This approach has been employed since the seminal work of Keller and Segel in the 1970s (Keller and Segel, J. Theor. Biol. 30:235–248, 1971). The system has been shown to permit travelling wave solutions which correspond to travelling band formation in bacterial colonies, yet only under specific criteria, such as a singularity in the chemotactic sensitivity function as the signal approaches zero. Such a singularity generates infinite macroscopic velocities which are biologically unrealistic. In this paper, we formulate a model that takes into consideration relevant details of the intracellular processes while avoiding the singularity in the chemotactic sensitivity. We prove the global existence of solutions and then show the existence of travelling wave solutions both numerically and analytically.     

Existence and stability of periodic travelling wave solutions to Nagumo's nerve equation

Summary Nagumo's nerve conduction equation has a one-parameter family of spatially periodic travelling wave solutions. First, we prove the existence of these solutions by using a topological method. (There are some exceptional cases in which this method cannot be applied in showing the existence.) A periodic travelling wave solution corresponds to a closed orbit of a third-order dynamical system. The Poincaré index of the closed orbit is determined as a direct consequence of the proof of the existence. Second, we prove that the periodic travelling wave solution is unstable if the Poincaré index of the corresponding closed orbit is + 1. By using this result, together with the result of the author's previous paper, it is concluded that the slow periodic travelling wave solutions are always unstable. Third, we consider the stability of the fast periodic travelling wave solutions. We denote by L(c) the spatial period of the travelling wave solution with the propagation speed c. It is shown that the fast solution is unstable if its period is close to L_min, the minimum of L(c).     

Local positive feedback by calcium in the propagation of intracellular calcium waves.

In many types of eukaryotic cells, the activation of surface receptors leads to the production of inositol 1,4,5-trisphosphate and calcium release from intracellular stores. Calcium release can occur in complex spatial patterns, including waves of release that traverse the cytoplasm. Fluorescence video microscopy was used to view calcium waves in single mouse neuroblastoma cells. The propagation of calcium waves was slowed by buffers that bind calcium quickly, such as BAPTA, but not by a buffer with slower on-rate, EGTA. This shows that a key feedback event in wave propagation is rapid diffusion of calcium occurring locally on a scale of < 1 micron. The length-speed product of wavefronts was used to determine that calcium acting in feedback diffuses at nearly the rate expected for free diffusion in aqueous solution. In cytoplasm, which contains immobile Ca2+ buffers, this rate of diffusion occurs only in the first 0.2 ms after release, within 0.4 micron of a Ca2+ release channel mouth. Calcium diffusion from an open channel to neighboring release sites is, therefore, a rate-determining regenerative step in calcium wave propagation. The theoretical limitations of the wave front analysis are discussed.     

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.     

From periodic travelling waves to travelling fronts in the spike-diffuse-spike model of dendritic waves

In the vertebrate brain excitatory synaptic contacts typically occur on the tiny evaginations of neuron dendritic surface known as dendritic spines. There is clear evidence that spine heads are endowed with voltage-dependent excitable channels and that action potentials invade spines. Computational models are being increasingly used to gain insight into the functional significance of a spine with an excitable membrane. The spike-diffuse-spike (SDS) model is one such model that admits to a relatively straightforward mathematical analysis. In this paper we demonstrate that not only can the SDS model support solitary travelling pulses, already observed numerically in more detailed biophysical models, but that it has periodic travelling wave solutions. The exact mathematical treatment of periodic travelling waves in the SDS model is used, within a kinematic framework, to predict the existence of connections between two periodic spike trains of different interspike interval. The associated wave front in the sequence of interspike intervals travels with a constant velocity without degradation of shape, and might therefore be used for the robust encoding of information.     

Effects of rapid buffers on Ca2+ diffusion and Ca2+ oscillations.

Based on realistic mechanisms of Ca2+ buffering that include both stationary and mobile buffers, we derive and investigate models of Ca2+ diffusion in the presence of rapid buffers. We obtain a single transport equation for Ca2+ that contains the effects caused by both stationary and mobile buffers. For stationary buffers alone, we obtain an expression for the effective diffusion constant of Ca2+ that depends on local Ca2+ concentrations. Mobile buffers, such as fura-2, BAPTA, or small endogenous proteins, give rise to a transport equation that is no longer strictly diffusive. Calculations are presented to show that these effects can modify greatly the manner and rate at which Ca2+ diffuses in cells, and we compare these results with recent measurements by Allbritton et al. (1992). As a prelude to work on Ca2+ waves, we use a simplified version of our model of the activation and inhibition of the IP3 receptor Ca2+ channel in the ER membrane to illustrate the way in which Ca2+ buffering can affect both the amplitude and existence of Ca2+ oscillations.     

Travelling waves in a tissue interaction model for skin pattern formation

Tissue interaction plays a major role in many morphogenetic processes, particularly those associated with skin organ primordia. We examine travelling wave solutions in a tissue interaction model for skin pattern formation which is firmly based on the known biology. From a phase space analysis we conjecture the existence of travelling waves with specific wave speeds. Subsequently, analytical approximations to the wave profiles are derived using perturbation methods. We then show numerically that such travelling wave solutions do exist and that they are in good agreement with our analytical results. Finally, the biological implications of our analysis are discussed.