Dynamics of excitable elements with time-delayed coupling

Calcium waves in a grid of clustered channels with synchronous IP3 binding and unbinding

Calcium signals in cells occur at multiple spatial scales and variable temporal duration. However, a physical explanation for transitions between long-lasting global oscillations and localized short-term elevations (puffs) of cytoplasmic Ca²⁺ is still lacking. Here we introduce a phenomenological, coarse-grained model for the calcium variable, which is represented by ordinary differential equations. Due to its small number of parameters, and its simplicity, this model allows us to numerically study the interplay of multi-scale calcium concentrations with stochastic ion channel gating dynamics even in larger systems. We apply this model to a single cluster of inositol trisphosphate (IP ₃) receptor channels and find further evidence for the results presented in earlier work: a single cluster may be capable of producing different calcium release types, where long-lasting events are accompanied by unbinding of IP ₃ from the receptor (Rückl et al., PLoS Comput. Biol. 11, e1003965 (2015)). Finally, we show the practicability of the model in a grid of 64 clusters which is computationally intractable with previous high-resolution models. Here long-lasting events can lead to synchronized oscillations and waves, while short events stay localized. The frequency of calcium releases as well as their coherence can thereby be regulated by the amplitude of IP ₃ stimulation. Finally the model allows for a new explanation of oscillating [IP ₃], which is not based on metabolic production and degradation of IP ₃.

Analytical and experimental study of two delay-coupled excitable units

We investigate the onset of time-periodic oscillations for a system of two identical delay-coupled excitable (nonoscillatory) units. We first analyze these solutions by using asymptotic methods. The oscillations are described as relaxation oscillations exhibiting successive slow and fast changes. The analysis highlights the determinant role of the delay during the fast transition layers. We then study experimentally a system of two coupled electronic circuits that is modeled mathematically by the same delay differential equations. We obtain quantitative agreements between analytical and experimental bifurcation diagrams.

Control of bursting behavior in neurons by autaptic modulation

Firing properties of biological neurons have long been recognized to be determined by extrinsic synaptic afferents that neurons receive and intrinsic ionic mechanisms that neurons possess, however, previous researches have also demonstrated that firing behavior of single neurons can be modulated by the neurons themselves, realized by the autapses. Thus in this investigation, we argued that how autaptic modulations shape the bursting behavior of biological neurons. We considered the issue from the following two aspects: autaptic-excitation and -inhibition. Our results suggested that for autaptic-excitation, under the condition of relatively weak stimulus, regular bursting was more incline to occur when the autaptic strength was weak, while regular spiking was more likely to appear when the autaptic strength was strong. However, larger stimulus would diminish the portion of bursting, but increase the portion of spiking. For autaptic-inhibition, under relatively weak stimulus, a wide range of regular bursting emerges when the autaptic strength was small, but when stronger stimulus were applied, the range of regular bursting shrinked into a small region. Meanwhile, we observed that synaptic delays have no obvious effects in the case of autaptic-excitation, while a subtle effect of synaptic delays was observed in the case of autaptic-inhibition. These results showed that bursting behavior of neurons could be controlled and modulated by the autaptic mechanisms that biological neurons intrinsically possess, and the final results may further promote the understanding in the generation of various neuronal firing patterns.

Traveling echo waves in an array of excitable elements with time-delayed coupling

We study time-delayed coupling in one-dimensional arrays of discrete excitable FitzHugh-Nagumo elements. For small time delays in an excitatory coupling, noise-triggered spikes undergo a slight prolongation. For larger delay times, the coupling to adjacent elements leads to a splitting of spikes and an "echo" effect that produces dampened series of spikes. For very large delay times, a transition to a stationary regime appears, as the spikes are sustained by mutual excitations, resulting in antiphase synchronization of all elements. We characterize several features exhibited by the system, such as speed and duration of pulses as well as bifurcations leading to the stationary regime.

Spontaneous spiking in an autaptic Hodgkin-Huxley setup

The effect of intrinsic channel noise is investigated for the dynamic response of a neuronal cell with a delayed feedback loop. The loop is based on the so-called autapse phenomenon in which dendrites establish connections not only to neighboring cells but also to its own axon. The biophysical modeling is achieved in terms of a stochastic Hodgkin-Huxley model containing such a built in delayed feedback. The fluctuations stem from intrinsic channel noise, being caused by the stochastic nature of the gating dynamics of ion channels. The influence of the delayed stimulus is systematically analyzed with respect to the coupling parameter and the delay time in terms of the interspike interval histograms and the average interspike interval. The delayed feedback manifests itself in the occurrence of bursting and a rich multimodal interspike interval distribution, exhibiting a delay-induced reduction in the spontaneous spiking activity at characteristic frequencies. Moreover, a specific frequency-locking mechanism is detected for the mean interspike interval.

Calcium domains around single and clustered IP3 receptors and their modulation by buffers

We study Ca(2+) release through single and clustered IP(3) receptor channels on the ER membrane under presence of buffer proteins. Our computational scheme couples reaction-diffusion equations and a Markovian channel model and allows our investigating the effects of buffer proteins on local calcium concentrations and channel gating. We find transient and stationary elevations of calcium concentrations around active channels and show how they determine release amplitude. Transient calcium domains occur after closing of isolated channels and constitute an important part of the channel's feedback. They cause repeated openings (bursts) and mediate increased release due to Ca(2+) buffering by immobile proteins. Stationary domains occur during prolonged activity of clustered channels, where the spatial proximity of IP(3)Rs produces a distinct [Ca(2+)] scale (0.5-10 microM), which is smaller than channel pore concentrations (>100 microM) but larger than transient levels. While immobile buffer affects transient levels only, mobile buffers in general reduce both transient and stationary domains, giving rise to Ca(2+) evacuation and biphasic modulation of release amplitude. Our findings explain recent experiments in oocytes and provide a general framework for the understanding of calcium signals.

Puff-wave transition in an inhomogeneous model for calcium signals

In many cell types, calcium ion channels on the endoplasmic reticulum membrane occur in a clustered distribution. The channels generate either localized puffs, each comprising channels of only one cluster, or global calcium waves. In this work we model the calcium system as a two-dimensional lattice of active elements distributed regularly in an otherwise passive space. We address an important feature of the puff-wave transition, which is the difference in lifetime of puffs at a few hundred milliseconds and long-lived global waves with periods of several seconds. We show that such a lifetime difference between puffs and waves can be understood with strongly reduced ordinary differential equations modified by a time-scale factor that takes into account the coupling strength of active and passive regions determined by the Ca2+ diffusion coefficient. Furthermore, we show that the point model can also describe very well the dependence of Ca2+ oscillation characteristics on the cluster-cluster distance in the case of large diffusivity.

Interplay of time-delayed feedback control and temporally correlated noise in excitable systems

The interplay of time-delayed feedback and temporally correlated coloured noise in a single and two coupled excitable systems is studied in the framework of the FitzHugh-Nagumo (FHN) model. By using coloured noise instead of white noise, the noise correlation time is introduced as an additional time scale. We show that in a single FHN system the major time scale of oscillations is strongly influenced by the noise correlation time, which in turn affects the maxima of coherence with respect to the delay time. In two coupled FHN systems, coloured noise input to one subsystem influences coherence resonance and stochastic synchronization of both subsystems. Application of delayed feedback control to the coloured noise-driven subsystem is shown to change coherence and time scales of noise-induced oscillations in both systems, and to enhance or suppress stochastic synchronization under certain conditions.