Phase-space structure of a thermoreceptor

Computational and experimental modulation of a noisy chaotic neuronal system

In this work, we study the interplay between chaos and noise in neuronal state transitions involving period doubling cascades. Our approach involves the implementation of a neuronal mathematical model under the action of neuromodulatory input, with and without noise, as well as equivalent experimental work on a biological neuron in the stomatogastric ganglion of the crab Cancer borealis. Our simulations show typical transitions between tonic and bursting regimes that are mediated by chaos and period doubling cascades. While this transition is less evident when intrinsic noise is present in the model, the noisy computational output displays features akin to our experimental results. The differences and similarities observed in the computational and experimental approaches are discussed.

Obstacle induced spiral waves in a multilayered Huber-Braun (HB) neuron model

Various dynamical properties of four-dimensional mammalian cold receptor model have been discussed widely in the literature considering noise and temperature as important parameters of discussion. Though various spiking and bursting behaviors of the neuron under various noise and temperature conditions studied for a single neuron, no much discussions have been done on the collective behavior. We investigate the collective behavior of these temperature dependent stochastic neurons and unlike the neuron models when forced by periodic external force there is no wave reentry or spiral waves in the network. Hence, we introduce obstacle in the network and depending on the orientation and size of the introduced obstacle, we could show their effects on the wave reentry in the network. Various significant discussions are produced in this paper to confirm that obstacles placed parallel to the wave entry affects the excitability of the tissues significantly compared to those obstacles place perpendicular. We could also show that those obstacles which are lesser in dimensions doesn't affect the excitabilities and hence doesn't contribute for wave reentry. We introduce a new technique to identify wave reentry and spiral waves using the period of individual nodes is proposed. This technique could help us identify even the lowest of excitability change which cannot be seen when using spatiotemporal snapshots.

Front-propagation in bacterial inter-colony communication

Many bacterial species exchange signaling molecules to coordinate population-wide responses. For this process, known as quorum sensing, the concentration of the respective molecules is crucial. Here, we consider the interaction between spatially distributed bacterial colonies so that the spreading of the signaling molecules in space becomes important. The exponential growth of the signal-producing populations and the corresponding increase in signaling molecule production result in an exponential concentration profile that spreads with uniform speed. The theoretical predictions are supported by experiments with different strains of the soil bacterium Sinorhizobium meliloti that display fluorescence when either producing or responding to the signaling molecules.

Modeling effects of neural fluctuations and inter-scale interactions

One of the greatest challenges to science, in particular, to neuroscience, is to understand how processes at different levels of organization are related to each other. In connection with this problem is the question of the functional significance of fluctuations, noise, and chaos. This paper deals with three related issues: (1) how processes at different organizational levels of neural systems might be related, (2) the functional significance of non-linear neurodynamics, including oscillations, chaos, and noise, and (3) how computational models can serve as useful tools in elucidating these types of issues. In order to capture and describe phenomena at different micro (molecular), meso (cellular), and macro (network) scales, the computational models need to be of appropriate complexity making use of available experimental data. I exemplify by two major types of computational models, those of Hans Braun and colleagues and those of my own group, which both aim at bridging gaps between different levels of neural systems. In particular, the constructive role of noise and chaos in such systems is modelled and related to functions, such as sensation, perception, learning/memory, decision making, and transitions between different (un-)conscious states. While there is, in general, a focus on upward causation, I will also discuss downward causation, where higher level activity may affect the activity at lower levels, which should be a condition for any functional role of consciousness and free will, often considered to be problematic to science.

Activity patterns with silent states in a heterogeneous network of gap-junction coupled Huber-Braun model neurons

A mathematical model of a network of nearest neighbor gap-junction coupled neurons has been used to examine the impact of neuronal heterogeneity on the networks' activity during increasing coupling strength. Heterogeneity has been introduced by Huber-Braun model neurons with randomization of the temperature as a scaling factor. This leads to neurons of an enormous diversity of impulse pattern, including burst discharges, chaotic activity, and two different types of tonic firing-all of them experimentally observed in the peripheral as well as central nervous system. When the network is composed of all these types of neurons, randomly selected, a particular phenomenon can be observed. At a certain coupling strength, the network goes into a completely silent state. Examination of voltage traces and inter-spike intervals of individual neurons suggests that all neurons, irrespective of their original pattern, go through a well-known bifurcation scenario, resembling those of single neurons especially on external current injection. All the originally spontaneously firing neurons can achieve constant membrane potentials at which all intrinsic and gap-junction currents are balanced. With limited diversity, i.e., taking out neurons of specific patterns from the lower and upper temperature range, spontaneous firing can be reinstalled with further increasing coupling strength, especially when the tonic firing regimes are missing. Reinstalled firing develops from slowly increasing subthreshold oscillations leading to tonic firing activity with already fairly well synchronized action potentials, while the subthreshold potentials can still be significantly different. Full in phase synchronization is not achieved. Additional studies are needed elucidating the underlying mechanisms and the conditions under which such particular transitions can appear.

Detection of nonstationary transition to synchronized states of a neural network using recurrence analyses

We study the stability of asymptotic states displayed by a complex neural network. We focus on the loss of stability of a stationary state of networks using recurrence quantifiers as tools to diagnose local and global stabilities as well as the multistability of a coupled neural network. Numerical simulations of a neural network composed of 1024 neurons in a small-world connection scheme are performed using the model of Braun et al. [Int. J. Bifurcation Chaos 08, 881 (1998)IJBEE40218-127410.1142/S0218127498000681], which is a modified model from the Hodgkin-Huxley model [J. Phys. 117, 500 (1952)]. To validate the analyses, the results are compared with those produced by Kuramoto's order parameter [Chemical Oscillations, Waves, and Turbulence (Springer-Verlag, Berlin Heidelberg, 1984)]. We show that recurrence tools making use of just integrated signals provided by the networks, such as local field potential (LFP) (LFP signals) or mean field values bring new results on the understanding of neural behavior occurring before the synchronization states. In particular we show the occurrence of different stationary and nonstationarity asymptotic states.

Chimera states in a Hodgkin-Huxley model of thermally sensitive neurons

Chimera states occur when identically coupled groups of nonlinear oscillators exhibit radically different dynamics, with one group exhibiting synchronized oscillations and the other desynchronized behavior. This dynamical phenomenon has recently been studied in computational models and demonstrated experimentally in mechanical, optical, and chemical systems. The theoretical basis of these states is currently under active investigation. Chimera behavior is of particular relevance in the context of neural synchronization, given the phenomenon of unihemispheric sleep and the recent observation of asymmetric sleep in human patients with sleep apnea. The similarity of neural chimera states to neural "bump" states, which have been suggested as a model for working memory and visual orientation tuning in the cortex, adds to their interest as objects of study. Chimera states have been demonstrated in the FitzHugh-Nagumo model of excitable cells and in the Hindmarsh-Rose neural model. Here, we demonstrate chimera states and chimera-like behaviors in a Hodgkin-Huxley-type model of thermally sensitive neurons both in a system with Abrams-Strogatz (mean field) coupling and in a system with Kuramoto (distance-dependent) coupling. We map the regions of parameter space for which chimera behavior occurs in each of the two coupling schemes.

Phase synchronization of coupled bursting neurons and the generalized Kuramoto model

Bursting neurons fire rapid sequences of action potential spikes followed by a quiescent period. The basic dynamical mechanism of bursting is the slow currents that modulate a fast spiking activity caused by rapid ionic currents. Minimal models of bursting neurons must include both effects. We considered one of these models and its relation with a generalized Kuramoto model, thanks to the definition of a geometrical phase for bursting and a corresponding frequency. We considered neuronal networks with different connection topologies and investigated the transition from a non-synchronized to a partially phase-synchronized state as the coupling strength is varied. The numerically determined critical coupling strength value for this transition to occur is compared with theoretical results valid for the generalized Kuramoto model.

Synchronization of bursting Hodgkin-Huxley-type neurons in clustered networks

We considered a clustered network of bursting neurons described by the Huber-Braun model. In the upper level of the network we used the connectivity matrix of the cat cerebral cortex network, and in the lower level each cortex area (or cluster) is modelled as a small-world network. There are two different coupling strengths related to inter- and intracluster dynamics. Each bursting cycle is composed of a quiescent period followed by a rapid chaotic sequence of spikes, and we defined a geometric phase which enables us to investigate the onset of synchronized bursting, as the state in which the neuron start bursting at the same time, whereas their spikes may remain uncorrelated. The bursting synchronization of a clustered network has been investigated using an order parameter and the average field of the network in order to identify regimes in which each cluster may display synchronized behavior, whereas the overall network does not. We introduce quantifiers to evaluate the relative contribution of each cluster in the partial synchronized behavior of the whole network. Our main finding is that we typically observe in the clustered network not a complete phase synchronized regime but instead a complex pattern of partial phase synchronization in which different cortical areas may be internally synchronized at distinct phase values, hence they are not externally synchronized, unless the coupling strengths are too large.

Temperature-dependent bursting pattern analysis by modified Plant model

Many electrophysiological properties of neuron including firing rates and rhythmical oscillation change in response to a temperature variation, but the mechanism underlying these correlations remains unverified. In this study, we analyzed various action potential (AP) parameters of bursting pacemaker neurons in the abdominal ganglion of Aplysia juliana to examine whether or not bursting patterns are altered in response to temperature change. Here we found that the inter-burst interval, burst duration, and number of spike during burst decreased as temperature increased. On the other hand, the numbers of bursts per minute and numbers of spikes per minute increased and then decreased, but interspike interval during burst firstly decreased and then increased. We also tested the reproducibility of temperature-dependent changes in bursting patterns and AP parameters. Finally we performed computational simulations of these phenomena by using a modified Plant model composed of equations with temperature-dependent scaling factors to mathematically clarify the temperature-dependent changes of bursting patterns in burst-firing neurons. Taken together, we found that the modified Plant model could trace the ionic mechanism underlying the temperature-dependent change in bursting pattern from experiments with bursting pacemaker neurons in the abdominal ganglia of Aplysia juliana.

Control of bursting synchronization in networks of Hodgkin-Huxley-type neurons with chemical synapses

Thermally sensitive neurons present bursting activity for certain temperature ranges, characterized by fast repetitive spiking of action potential followed by a short quiescent period. Synchronization of bursting activity is possible in networks of coupled neurons, and it is sometimes an undesirable feature. Control procedures can suppress totally or partially this collective behavior, with potential applications in deep-brain stimulation techniques. We investigate the control of bursting synchronization in small-world networks of Hodgkin-Huxley-type thermally sensitive neurons with chemical synapses through two different strategies. One is the application of an external time-periodic electrical signal and another consists of a time-delayed feedback signal. We consider the effectiveness of both strategies in terms of protocols of applications suitable to be applied by pacemakers.

Temperature-dependent stochastic dynamics of the Huber-Braun neuron model

The response of a four-dimensional mammalian cold receptor model to different implementations of noise is studied across a wide temperature range. It is observed that for noisy activation kinetics, the parameter range decomposes into two regions in which the system reacts qualitatively completely different to small perturbations through noise, and these regions are separated by a homoclinic bifurcation. Noise implemented as an additional current yields a substantially different system response at low temperature values, while the response at high temperatures is comparable to activation-kinetic noise. We elucidate how this phenomenon can be understood in terms of state space dynamics and gives quantitative results on the statistics of interspike interval distributions across the relevant parameter range.

Noise-induced precursors of tonic-to-bursting transitions in hypothalamic neurons and in a conductance-based model

The dynamics of neurons is characterized by a variety of different spiking patterns in response to external stimuli. One of the most important transitions in neuronal response patterns is the transition from tonic firing to burst discharges, i.e., when the neuronal activity changes from single spikes to the grouping of spikes. An increased number of interspike-interval sequences of specific temporal correlations was detected in anticipation of temperature induced tonic-to-bursting transitions in both, experimental impulse recordings from hypothalamic brain slices and numerical simulations of a stochastic model. Analysis of the modelling data elucidates that the appearance of such patterns can be related to particular system dynamics in the vicinity of the period-doubling bifurcation. It leads to a nonlinear response on de- and hyperpolarizing perturbations introduced by noise. This explains why such particular patterns can be found as reliable precursors of the neurons' transition to burst discharges.

On the role of subthreshold currents in the Huber-Braun cold receptor model

We study the role of the strength of subthreshold currents in a four-dimensional Hodgkin-Huxley-type model of mammalian cold receptors. Since a total diminution of subthreshold activity corresponds to a decomposition of the model into a slow, subthreshold, and a fast, spiking subsystem, we first elucidate their respective dynamics separately and draw conclusions about their role for the generation of different spiking patterns. These results motivate a numerical bifurcation analysis of the effect of varying the strength of subthreshold currents, which is done by varying a suitable control parameter. We work out the key mechanisms which can be attributed to subthreshold activity and furthermore elucidate the dynamical backbone of different activity patterns generated by this model.

Dynamics of globally delay-coupled neurons displaying subthreshold oscillations

We study an ensemble of neurons that are coupled through their time-delayed collective mean field. The individual neuron is modelled using a Hodgkin-Huxley-type conductance model with parameters chosen such that the uncoupled neuron displays autonomous subthreshold oscillations of the membrane potential. We find that the ensemble generates a rich variety of oscillatory activities that are mainly controlled by two time scales: the natural period of oscillation at the single neuron level and the delay time of the global coupling. When the neuronal oscillations are synchronized, they can be either in-phase or out-of-phase. The phase-shifted activity is interpreted as the result of a phase-flip bifurcation, also occurring in a set of globally delay-coupled limit cycle oscillators. At the bifurcation point, there is a transition from in-phase to out-of-phase (or vice versa) synchronized oscillations, which is accompanied by an abrupt change in the common oscillation frequency. This phase-flip bifurcation was recently investigated in two mutually delay-coupled oscillators and can play a role in the mechanisms by which the neurons switch among different firing patterns.

Interplay of subthreshold activity, time-delayed feedback, and noise on neuronal firing patterns

Feedback connections and noise are ubiquitous features of neuronal networks and affect in a determinant way the patterns of neural activity. Here we study how the subthreshold dynamics of a neuron interacts with time-delayed feedback and noise. We use a Hodgkin-Huxley-type model of a thermoreceptor neuron and assume the feedback to be linear, corresponding effectively to a recurrent electrical connection via gap junctions. This type of feedback can model electrical autapses, which connect the terminal fibers of a neuron's axon with dendrites from the same neuron. Thus the delay in the feedback loop is due basically to the axonal propagation time. We chose model parameters for which the neuron displays, in the absence of feedback and noise, only subthreshold oscillations. These oscillations, however, take the neuron close to the firing threshold, such that small perturbations can drive it above the level for generation of action potentials. The resulting interplay between weak delayed feedback, noise, and the subthreshold intrinsic activity is nontrivial. For negative feedback, depending on the delay, the firing rate can be lower than in the noise-free situation. This is due to the fact that noise inhibits feedback-induced spikes by driving the neuronal oscillations away from the firing threshold. For positive feedback, there are regions of delay values where the noise-induced spikes are inhibited by the feedback; in this case, it is the feedback that drives the neuronal oscillations away from the threshold. Our study contributes to a better understanding of the role of electrical self-connections in the presence of noise and subthreshold activity.

Propagation effects of current and conductance noise in a model neuron with subthreshold oscillations

We have examined the effects of current and conductance noise in a single-neuron model which can generate a variety of physiologically important impulse patterns. Current noise enters the membrane equation directly while conductance noise is propagated through the activation variables. Additive Gaussian white noise which is implemented as conductance noise appears in the voltage equations as an additive and a multiplicative term. Moreover, the originally white noise is turned into colored noise. The noise correlation time is a function of the system's control parameters which may explain the different effects of current and conductance noise in different dynamic states. We have found the most significant, qualitative differences between different noise implementations in a pacemaker-like, tonic firing regime at the transition to chaotic burst discharges. This reflects a dynamic state of high physiological relevance.

Frustration, drift, and antiphase coupling in a neural array

Synchronization among neurons is critical for many processes in the nervous system, ranging from the processing of sensory information to the onset of pathological conditions such as epilepsy. Here, we study synchronization in an array of neurons, each modeled by a set of nonlinear ordinary differential equations. We find that an array of 20x20 coupled neurons undergoes a series of alternating low and high synchronization states, as measured by phase-locking and frequency entrainment, as the coupling constant is tuned. The role of long-range connections in inducing "small-world networks" has recently been of great interest in many physical and biological problems. Since long-range connections do exist in the brain, we investigated the role of such connections in our neural array. Introducing a biologically realistic percentage of long-range connections has no significant effect on synchronization. We find that it is rather the type of coupling and the total number of connections that determine the synchronization state of the array. We also show that some coupling conditions can lead to frustration in the system, resulting from an inability to simultaneously satisfy conflicting phase requirements. This frustration leads to a drift in the overall behavior of the network, which may offer an explanation for transitions between different types of neural oscillations observed experimentally.

Impulse pattern in bi-directionally coupled model neurons of different dynamics

The effects of bi-directional gap junction coupling of two model neurons with subthreshold oscillations have been examined when the individual neurons are operating at different dynamical states either in the tonic or bursting firing mode. Our simulations indicate that intermediate coupling strengths mostly lead to highly variable, often chaotic impulse patterns whereas transition to completely synchronized activity at high coupling strengths is generally going along with transitions to regular limit cycle activity. The synchronized activity pattern, however, can be completely different from the original pattern of the uncoupled neurons.

Zero-lag long-range synchronization via dynamical relaying

We show that isochronous synchronization between two delay-coupled oscillators can be achieved by relaying the dynamics via a third mediating element, which surprisingly lags behind the synchronized outer elements. The zero-lag synchronization thus obtained is robust over a considerable parameter range. We substantiate our claims with experimental and numerical evidence of such synchronization solutions in a chain of three coupled semiconductor lasers with long interelement coupling delays. The generality of the mechanism is validated in a neuronal model with the same coupling architecture. Thus, our results show that zero-lag synchronized chaotic dynamical states can occur over long distances through relaying, without restriction by the amount of delay.

Ordering spatiotemporal chaos in complex thermosensitive neuron networks

We have studied the effect of random long-range connections in chaotic thermosensitive neuron networks with each neuron being capable of exhibiting diverse bursting behaviors, and found stochastic synchronization and optimal spatiotemporal patterns. For a given coupling strength, the chaotic burst-firings of the neurons become more and more synchronized as the number of random connections (or randomness) is increased and, rather, the most pronounced spatiotemporal pattern appears for an optimal randomness. As the coupling strength is increased, the optimal randomness shifts towards a smaller strength. This result shows that random long-range connections can tame the chaos in the neural networks and make the neurons more effectively reach synchronization. Since the model studied can be used to account for hypothalamic neurons of dogfish, catfish, etc., this result may reflect the significant role of random connections in transferring biological information.

Influence of time-delayed feedback in the firing pattern of thermally sensitive neurons

We explore the dynamics of a Hodgkin-Huxley-type model for thermally sensitive neurons that exhibit intrinsic oscillatory activity. The model is modified to include a feedback loop that is represented by two parameters: the synaptic strength and the transmission delay time. We analyze the dynamics of the neuron depending on the temperature, the synaptic strength, and the delay time. We find parameter regions where the effect of the recurrent connexion is excitatory, inducing spikes or trains of spikes, and regions where it is inhibitory, reducing or eliminating completely the spiking behavior. We characterize the complex interplay of the intrinsic dynamics of the neuron with the recurrent feedback input and a noisy input.

Bifurcating neuron: computation and learning

The ability of bifurcating processing units and their networks to rapidly switch between different dynamic modes has been used in recent research efforts to model new computational properties of neural systems. In this spirit, we devise a bifurcating neuron based on control of chaos collapsing to a period-3 orbit in the dynamics of a quadratic logistic map (QLM). Proposed QLM3 neuron is constructed with the third iterate of QLM and uses an external input, which governs its dynamics. The input shifts the neuron's dynamics from chaos to one of the stable fixed points. This way the inputs from certain ranges (clusters) are mapped to stable fixed points, while the rest of the inputs is mapped to chaotic or periodic output dynamics. It has been shown that QLM3 neuron is able to learn a specific mapping by adaptively adjusting its bifurcation parameter, the idea of which is based on the principles of parametric control of logistic maps [Proceedings of the International Symposium on Nonlinear Theory and its Applications (NOLTA'97), Honolulu, HI, 1997; Proceedings of SPIE, 2000]. Learning algorithm for the bifurcation parameter is proposed, which employs the error gradient descent method.

Oscillations, resonances and noise: basis of flexible neuronal pattern generation

Modulation of neuronal impulse pattern is examined by means of a simplified Hodgkin-Huxley type computer model which refers to experimental recordings of cold receptor discharges. This model essentially consists of two potentially oscillating subsystems: a spike generator and a subthreshold oscillator. With addition of noise the model successfully mimics the major types of experimentally recorded impulse patterns and thereby elucidate different resonance behaviors. (1) There is a range of rhythmic spiking or bursting where the spike generator is strongly coupled to the subthreshold oscillator. (2) There is a pacemaker activity of more complex interactions where the spike generator has overtaken part of the control. (3) There is a situation where the two subsystems are decoupled and only resonate with the help of noise.

Anticipating the response of excitable systems driven by random forcing

We study the regime of anticipated synchronization in unidirectionally coupled model neurons subject to a common external aperiodic forcing that makes their behavior unpredictable. We show numerically and by analog hardware electronic circuits that, under appropriate coupling conditions, the pulses fired by the slave neuron anticipate (i.e., predict) the pulses fired by the master neuron. This anticipated synchronization occurs even when the common external forcing is white noise.

Noise-induced impulse pattern modifications at different dynamical period-one situations in a computer model of temperature encoding

We used a minimal Hodgkin-Huxley type model of cold receptor discharges to examine how noise interferes with the non-linear dynamics of the ionic mechanisms of neuronal stimulus encoding. The model is based on the assumption that spike-generation depends on subthreshold oscillations. With physiologically plausible temperature scaling, it passes through different impulse patterns which, with addition of noise, are in excellent agreement with real experimental data. The interval distributions of purely deterministic simulations, however, exhibit considerable differences compared to the noisy simulations especially at the bifurcations of deterministically period-one discharges. We, therefore, analyzed the effects of noise in different situations of deterministically regular period-one discharges: (1) at high-temperatures near the transition to subthreshold oscillations and to burst discharges, and (2) at low-temperatures close to and more far away from the bifurcations to chaotic dynamics. The data suggest that addition of noise can considerably extend the dynamical behavior of the system with coexistence of different dynamical situations at deterministically fixed parameter constellations. Apart from well-described coexistence of spike-generating and subthreshold oscillations also mixtures of tonic and bursting patterns can be seen and even transitions to unstable period-one orbits seem to appear. The data indicate that cooperative effects between low- and high-dimensional dynamics have to be considered as qualitatively important factors in neuronal encoding.