Synchronization and propagation of bursts in networks of coupled map neurons

Complex dynamics in a fractional order nephron pressure and flow regulation model

Cardiovascular diseases can be attributed to irregular blood pressure, which may be caused by malfunctioning kidneys that regulate blood pressure. Research has identified complex oscillations in the mechanisms used by the kidney to regulate blood pressure. This study uses established physiological knowledge and earlier autoregulation models to derive a fractional order nephron autoregulation model. The dynamical behaviour of the model is analyzed using bifurcation plots, revealing periodic oscillations, chaotic regions, and multistability. A lattice array of the model is used to study collective behaviour and demonstrates the presence of chimeras in the network. A ring network of the fractional order model is also considered, and a diffusion coupling strength is adopted. A basin of synchronization is derived, considering coupling strength, fractional order or number of neighbours as parameters, and measuring the strength of incoherence. Overall, the study provides valuable insights into the complex dynamics of the nephron autoregulation model and its potential implications for cardiovascular diseases.

A discrete Huber-Braun neuron model: from nodal properties to network performance

Many of the well-known neuron models are continuous time systems with complex mathematical definitions. Literatures have shown that a discrete mathematical model can effectively replicate the complete dynamical behaviour of a neuron with much reduced complexity. Hence, we propose a new discrete neuron model derived from the Huber-Braun neuron with two additional slow and subthreshold currents alongside the ion channel currents. We have also introduced temperature dependent ion channels to study its effects on the firing pattern of the neuron. With bifurcation and Lyapunov exponents we showed the chaotic and periodic regions of the discrete model. Further to study the complexity of the neuron model, we have used the sample entropy algorithm. Though the individual neuron analysis gives us an idea about the dynamical properties, it's the collective behaviour which decides the overall behavioural pattern of the neuron. Hence, we investigate the spatiotemporal behaviour of the discrete neuron model in single- and two-layer network. We have considered obstacle as an important factor which changes the excitability of the neurons in the network. Literatures have shown that spiral waves can play a positive role in breaking through quiescent areas of the brain as a pacemaker by creating a coherence resonance behaviour. Hence, we are interested in studying the induced spiral waves in the network. In this condition when an obstacle is introduced the wave propagation is disturbed and we could see multiple wave re-entry and spiral waves. In a two-layer network when the obstacle is considered only in one layer and stimulus applied to the layer having the obstacle, the wave re-entry is seen in both the layer though the other layer is not exposed to obstacle. But when both the layers are inserted with an obstacle and stimuli also applied to the layers, they behave like independent layers with no coupling effect. In a two-layer network, stimulus play an important role in spatiotemporal dynamics of the network.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11571-022-09806-1.

Intermittent evolution routes to the periodic or the chaotic orbits in Rulkov map

This paper concerns the intermittent evolution routes to the asymptotic regimes in the Rulkov map. That is, the windows with transient approximate periodic and transient chaotic behaviors occur alternatively before the system reaches the periodic or the chaotic orbits. Meanwhile, the evolution routes to chaotic orbits can be classified into different types according to the windows before reaching asymptotic chaotic states. In addition, the initial values can be regarded as a key factor affecting the asymptotic behaviors and the evolution routes. The effects of the initial values are given by parameter planes, bifurcation diagrams, and waveforms. In order to investigate whether the intermittent evolution routes can be learned by machine learning, some experiments are given to understanding the differences between the trajectories of the Rulkov map generated by the numerical simulations and predicted by the neural networks. These results show that there is about 60% accuracy rate of successfully predicting both the evolution routes and the asymptotic period-3 orbits using a three-layer feedforward neural network, while the bifurcation diagrams can be reconstructed using reservoir computing except a few parameter conditions.

Cluster burst synchronization in a scale-free network of inhibitory bursting neurons

We consider a scale-free network of inhibitory Hindmarsh-Rose (HR) bursting neurons, and make a computational study on coupling-induced cluster burst synchronization by varying the average coupling strength J 0 . For sufficiently small J 0 , non-cluster desynchronized states exist. However, when passing a critical point J c ∗ ( ≃ 0.16 ) , the whole population is segregated into 3 clusters via a constructive role of synaptic inhibition to stimulate dynamical clustering between individual burstings, and thus 3-cluster desynchronized states appear. As J 0 is further increased and passes a lower threshold J l ∗ ( ≃ 0.78 ) , a transition to 3-cluster burst synchronization occurs due to another constructive role of synaptic inhibition to favor population synchronization. In this case, HR neurons in each cluster make burstings every 3rd cycle of the instantaneous burst rate R w ( t ) of the whole population, and exhibit burst synchronization. However, as J 0 passes an intermediate threshold J m ∗ ( ≃ 5.2 ) , HR neurons fire burstings intermittently at a 4th cycle of R w ( t ) via burst skipping rather than at its 3rd cycle, and hence they begin to make intermittent hoppings between the 3 clusters. Due to such intermittent intercluster hoppings via burst skippings, the 3 clusters become broken up (i.e., the 3 clusters are integrated into a single one). However, in spite of such break-up (i.e., disappearance) of the 3-cluster states, (non-cluster) burst synchronization persists in the whole population, which is well visualized in the raster plot of burst onset times where bursting stripes (composed of burst onset times and indicating burst synchronization) appear successively. With further increase in J 0 , intercluster hoppings are intensified, and bursting stripes also become dispersed more and more due to a destructive role of synaptic inhibition to spoil the burst synchronization. Eventually, when passing a higher threshold J h ∗ ( ≃ 17.8 ) a transition to desynchronization occurs via complete overlap between the bursting stripes. Finally, we also investigate the effects of stochastic noise on both 3-cluster burst synchronization and intercluster hoppings.

Synchronization of Rulkov neuron networks coupled by excitatory and inhibitory chemical synapses

This paper takes into account a neuron network model in which the excitatory and the inhibitory Rulkov neurons interact each other through excitatory and inhibitory chemical coupling, respectively. Firstly, for two or more identical or non-identical Rulkov neurons, the existence conditions of the synchronization manifold of the fixed points are investigated, which have received less attention over the past decades. Secondly, the master stability equation of the arbitrarily connected neuron network under the existence conditions of the synchronization manifold is discussed. Thirdly, taking three identical Rulkov neurons as an example, some new results are presented: (1) topological structures that can make the synchronization manifold exist are given, (2) the stability of synchronization when different parameters change is discussed, and (3) the roles of the control parameters, the ratio, as well as the size of the coupling strength and sigmoid function are analyzed. Finally, for the chemical coupling between two non-identical neurons, the transversal system is given and the effect of two coupling strengths on synchronization is analyzed.

Burst synchronization in a scale-free neuronal network with inhibitory spike-timing-dependent plasticity

We are concerned about burst synchronization (BS), related to neural information processes in health and disease, in the Barabási-Albert scale-free network (SFN) composed of inhibitory bursting Hindmarsh-Rose neurons. This inhibitory neuronal population has adaptive dynamic synaptic strengths governed by the inhibitory spike-timing-dependent plasticity (iSTDP). In previous works without considering iSTDP, BS was found to appear in a range of noise intensities for fixed synaptic inhibition strengths. In contrast, in our present work, we take into consideration iSTDP and investigate its effect on BS by varying the noise intensity. Our new main result is to find occurrence of a Matthew effect in inhibitory synaptic plasticity: good BS gets better via LTD, while bad BS get worse via LTP. This kind of Matthew effect in inhibitory synaptic plasticity is in contrast to that in excitatory synaptic plasticity where good (bad) synchronization gets better (worse) via LTP (LTD). We note that, due to inhibition, the roles of LTD and LTP in inhibitory synaptic plasticity are reversed in comparison with those in excitatory synaptic plasticity. Moreover, emergences of LTD and LTP of synaptic inhibition strengths are intensively investigated via a microscopic method based on the distributions of time delays between the pre- and the post-synaptic burst onset times. Finally, in the presence of iSTDP we investigate the effects of network architecture on BS by varying the symmetric attachment degree l ∗ and the asymmetry parameter Δ l in the SFN.

Effect of spike-timing-dependent plasticity on stochastic burst synchronization in a scale-free neuronal network

We consider an excitatory population of subthreshold Izhikevich neurons which cannot fire spontaneously without noise. As the coupling strength passes a threshold, individual neurons exhibit noise-induced burstings. This neuronal population has adaptive dynamic synaptic strengths governed by the spike-timing-dependent plasticity (STDP). However, STDP was not considered in previous works on stochastic burst synchronization (SBS) between noise-induced burstings of sub-threshold neurons. Here, we study the effect of additive STDP on SBS by varying the noise intensity D in the Barabási-Albert scale-free network (SFN). One of our main findings is a Matthew effect in synaptic plasticity which occurs due to a positive feedback process. Good burst synchronization (with higher bursting measure) gets better via long-term potentiation (LTP) of synaptic strengths, while bad burst synchronization (with lower bursting measure) gets worse via long-term depression (LTD). Consequently, a step-like rapid transition to SBS occurs by changing D, in contrast to a relatively smooth transition in the absence of STDP. We also investigate the effects of network architecture on SBS by varying the symmetric attachment degree [Formula: see text] and the asymmetry parameter [Formula: see text] in the SFN, and Matthew effects are also found to occur by varying [Formula: see text] and [Formula: see text]. Furthermore, emergences of LTP and LTD of synaptic strengths are investigated in details via our own microscopic methods based on both the distributions of time delays between the burst onset times of the pre- and the post-synaptic neurons and the pair-correlations between the pre- and the post-synaptic instantaneous individual burst rates (IIBRs). Finally, a multiplicative STDP case (depending on states) with soft bounds is also investigated in comparison with the additive STDP case (independent of states) with hard bounds. Due to the soft bounds, a Matthew effect with some quantitative differences is also found to occur for the case of multiplicative STDP.

Frequency-domain order parameters for the burst and spike synchronization transitions of bursting neurons

We are interested in characterization of synchronization transitions of bursting neurons in the frequency domain. Instantaneous population firing rate (IPFR) [Formula: see text], which is directly obtained from the raster plot of neural spikes, is often used as a realistic collective quantity describing population activities in both the computational and the experimental neuroscience. For the case of spiking neurons, a realistic time-domain order parameter, based on [Formula: see text], was introduced in our recent work to characterize the spike synchronization transition. Unlike the case of spiking neurons, the IPFR [Formula: see text] of bursting neurons exhibits population behaviors with both the slow bursting and the fast spiking timescales. For our aim, we decompose the IPFR [Formula: see text] into the instantaneous population bursting rate [Formula: see text] (describing the bursting behavior) and the instantaneous population spike rate [Formula: see text] (describing the spiking behavior) via frequency filtering, and extend the realistic order parameter to the case of bursting neurons. Thus, we develop the frequency-domain bursting and spiking order parameters which are just the bursting and spiking "coherence factors" [Formula: see text] and [Formula: see text] of the bursting and spiking peaks in the power spectral densities of [Formula: see text] and [Formula: see text] (i.e., "signal to noise" ratio of the spectral peak height and its relative width). Through calculation of [Formula: see text] and [Formula: see text], we obtain the bursting and spiking thresholds beyond which the burst and spike synchronizations break up, respectively. Consequently, it is shown in explicit examples that the frequency-domain bursting and spiking order parameters may be usefully used for characterization of the bursting and the spiking transitions, respectively.

Local and global synchronization transitions induced by time delays in small-world neuronal networks with chemical synapses

Effects of time delay on the local and global synchronization in small-world neuronal networks with chemical synapses are investigated in this paper. Numerical results show that, for both excitatory and inhibitory coupling types, the information transmission delay can always induce synchronization transitions of spiking neurons in small-world networks. In particular, regions of in-phase and out-of-phase synchronization of connected neurons emerge intermittently as the synaptic delay increases. For excitatory coupling, all transitions to spiking synchronization occur approximately at integer multiples of the firing period of individual neurons; while for inhibitory coupling, these transitions appear at the odd multiples of the half of the firing period of neurons. More importantly, the local synchronization transition is more profound than the global synchronization transition, depending on the type of coupling synapse. For excitatory synapses, the local in-phase synchronization observed for some values of the delay also occur at a global scale; while for inhibitory ones, this synchronization, observed at the local scale, disappears at a global scale. Furthermore, the small-world structure can also affect the phase synchronization of neuronal networks. It is demonstrated that increasing the rewiring probability can always improve the global synchronization of neuronal activity, but has little effect on the local synchronization of neighboring neurons.

Noise-induced burst and spike synchronizations in an inhibitory small-world network of subthreshold bursting neurons

We are interested in noise-induced firings of subthreshold neurons which may be used for encoding environmental stimuli. Noise-induced population synchronization was previously studied only for the case of global coupling, unlike the case of subthreshold spiking neurons. Hence, we investigate the effect of complex network architecture on noise-induced synchronization in an inhibitory population of subthreshold bursting Hindmarsh-Rose neurons. For modeling complex synaptic connectivity, we consider the Watts-Strogatz small-world network which interpolates between regular lattice and random network via rewiring, and investigate the effect of small-world connectivity on emergence of noise-induced population synchronization. Thus, noise-induced burst synchronization (synchrony on the slow bursting time scale) and spike synchronization (synchrony on the fast spike time scale) are found to appear in a synchronized region of the [Formula: see text]-[Formula: see text] plane ([Formula: see text]: synaptic inhibition strength and [Formula: see text]: noise intensity). As the rewiring probability [Formula: see text] is decreased from 1 (random network) to 0 (regular lattice), the region of spike synchronization shrinks rapidly in the [Formula: see text]-[Formula: see text] plane, while the region of the burst synchronization decreases slowly. We separate the slow bursting and the fast spiking time scales via frequency filtering, and characterize the noise-induced burst and spike synchronizations by employing realistic order parameters and statistical-mechanical measures introduced in our recent work. Thus, the bursting and spiking thresholds for the burst and spike synchronization transitions are determined in terms of the bursting and spiking order parameters, respectively. Furthermore, we also measure the degrees of burst and spike synchronizations in terms of the statistical-mechanical bursting and spiking measures, respectively.

Universal dynamical properties preclude standard clustering in a large class of biochemical data

MOTIVATION

Clustering of chemical and biochemical data based on observed features is a central cognitive step in the analysis of chemical substances, in particular in combinatorial chemistry, or of complex biochemical reaction networks. Often, for reasons unknown to the researcher, this step produces disappointing results. Once the sources of the problem are known, improved clustering methods might revitalize the statistical approach of compound and reaction search and analysis. Here, we present a generic mechanism that may be at the origin of many clustering difficulties.

RESULTS

The variety of dynamical behaviors that can be exhibited by complex biochemical reactions on variation of the system parameters are fundamental system fingerprints. In parameter space, shrimp-like or swallow-tail structures separate parameter sets that lead to stable periodic dynamical behavior from those leading to irregular behavior. We work out the genericity of this phenomenon and demonstrate novel examples for their occurrence in realistic models of biophysics. Although we elucidate the phenomenon by considering the emergence of periodicity in dependence on system parameters in a low-dimensional parameter space, the conclusions from our simple setting are shown to continue to be valid for features in a higher-dimensional feature space, as long as the feature-generating mechanism is not too extreme and the dimension of this space is not too high compared with the amount of available data.

AVAILABILITY AND IMPLEMENTATION

For online versions of super-paramagnetic clustering see http://stoop.ini.uzh.ch/research/clustering.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Exact synchronization bound for coupled time-delay systems

We obtain an exact bound for synchronization in coupled time-delay systems using the generalized Halanay inequality for the general case of time-dependent delay, coupling, and coefficients. Furthermore, we show that the same analysis is applicable to both uni- and bidirectionally coupled time-delay systems with an appropriate evolution equation for their synchronization manifold, which can also be defined for different types of synchronization. The exact synchronization bound assures an exponential stabilization of the synchronization manifold which is crucial for applications. The analytical synchronization bound is independent of the nature of the modulation and can be applied to any time-delay system satisfying a Lipschitz condition. The analytical results are corroborated numerically using the Ikeda system.

Chaotic phase synchronization in a modular neuronal network of small-world subnetworks

We investigate the onset of chaotic phase synchronization of bursting oscillators in a modular neuronal network of small-world subnetworks. A transition to mutual phase synchronization takes place on the bursting time scale of coupled oscillators, while on the spiking time scale, they behave asynchronously. It is shown that this bursting synchronization transition can be induced not only by the variations of inter- and intra-coupling strengths but also by changing the probability of random links between different subnetworks. We also analyze the effect of external chaotic phase synchronization of bursting behavior in this clustered network by an external time-periodic signal applied to a single neuron. Simulation results demonstrate a frequency locking tongue in the driving parameter plane, where bursting synchronization is maintained, even with the external driving. The width of this synchronization region increases with the signal amplitude and the number of driven neurons but decreases rapidly with the network size. Considering that the synchronization of bursting neurons is thought to play a key role in some pathological conditions, the presented results could have important implications for the role of externally applied driving signal in controlling bursting activity in neuronal ensembles.

Chaotic phase synchronization in small-world networks of bursting neurons

We investigate the chaotic phase synchronization in a system of coupled bursting neurons in small-world networks. A transition to mutual phase synchronization takes place on the bursting time scale of coupled oscillators, while on the spiking time scale, they behave asynchronously. It is shown that phase synchronization is largely facilitated by a large fraction of shortcuts, but saturates when it exceeds a critical value. We also study the external chaotic phase synchronization of bursting oscillators in the small-world network by a periodic driving signal applied to a single neuron. It is demonstrated that there exists an optimal small-world topology, resulting in the largest peak value of frequency locking interval in the parameter plane, where bursting synchronization is maintained, even with the external driving. The width of this interval increases with the driving amplitude, but decrease rapidly with the network size. We infer that the externally applied driving parameters outside the frequency locking region can effectively suppress pathologically synchronized rhythms of bursting neurons in the brain.

Delay-induced intermittent transition of synchronization in neuronal networks with hybrid synapses

We study the dependence of synchronization transitions in scale-free networks of bursting neurons with hybrid synapses on the information transmission delay and the probability of inhibitory synapses. It is shown that, irrespective of the probability of inhibitory synapses, the delay always plays a subtle role during synchronization transition of the scale-free neuronal networks. In particular, regions of irregular and regular propagating excitatory fronts appear intermittently as the delay increases. These delay-induced synchronization transitions are manifested as well-expressed minima in the measure for spatiotemporal synchrony. In addition, it is found that, for smaller and larger probability of inhibitory synapses, intermittent synchronization transition is relatively profound, while for the moderate probability of inhibitory synapses, synchronization transition seems less profound. More interestingly, it is found that as the probability of inhibitory synapses is large, regions of synchronization are upscattering.

Hybrid discrete-time neural networks

Hybrid dynamical systems combine evolution equations with state transitions. When the evolution equations are discrete-time (also called map-based), the result is a hybrid discrete-time system. A class of biological neural network models that has recently received some attention falls within this category: map-based neuron models connected by means of fast threshold modulation (FTM). FTM is a connection scheme that aims to mimic the switching dynamics of a neuron subject to synaptic inputs. The dynamic equations of the neuron adopt different forms according to the state (either firing or not firing) and type (excitatory or inhibitory) of their presynaptic neighbours. Therefore, the mathematical model of one such network is a combination of discrete-time evolution equations with transitions between states, constituting a hybrid discrete-time (map-based) neural network. In this paper, we review previous work within the context of these models, exemplifying useful techniques to analyse them. Typical map-based neuron models are low-dimensional and amenable to phase-plane analysis. In bursting models, fast-slow decomposition can be used to reduce dimensionality further, so that the dynamics of a pair of connected neurons can be easily understood. We also discuss a model that includes electrical synapses in addition to chemical synapses with FTM. Furthermore, we describe how master stability functions can predict the stability of synchronized states in these networks. The main results are extended to larger map-based neural networks.

Spatial-temporal dynamics of chaotic behavior in cultured hippocampal networks

Using multiple nonlinear techniques, we revealed the existence of chaos in the spontaneous activity of neuronal networks in vitro. The spatial-temporal dynamics of these networks indicated that emergent transition between chaotic behavior and superburst occurred periodically in low-frequency oscillations. An analysis of network-wide activity indicated that chaos was synchronized among different sites. Moreover, we found that the degree of chaos increased as the number of active sites in the network increased during long-term development (over three months in vitro). The chaotic behavior of the dissociated networks had similar spatial-temporal characteristics (rapid transition, periodicity, and synchronization) as the intact brain; however, the degree of chaos depended on the number of active sites at the mesoscopic level. This work could provide insight into neural coding and neurocybernetics.

Bursting regimes in map-based neuron models coupled through fast threshold modulation

A system consisting of two map-based neurons coupled through reciprocal excitatory or inhibitory chemical synapses is discussed. After a brief explanation of the basic mechanism behind generation and synchronization of bursts, parameter space is explored to determine less obvious but biologically meaningful regimes and effects. Among them, we show how excitatory synapses without any delays may induce antiphase synchronization; that a synapse may change its character from excitatory to inhibitory and vice versa by changing its conductance, without any change in reversal potential; and that small variations in the synaptic threshold may result in drastic changes in the synchronization of spikes within bursts. Finally we show how the synchronization effects found in the two-neuron system carry over to larger networks.

Sensitivity versus resonance in two-dimensional spiking-bursting neuron models

Through phase plane analysis of a class of two-dimensional spiking and bursting neuron models, covering some of the most popular map-based neuron models, we show that there exists a trade-off between the sensitivity of the neuron to steady external stimulation and its resonance properties, and how this trade-off may be tuned by the neutral or asymptotic character of the slow variable. Implications of the results for the suprathreshold behavior of the neurons, both by themselves and as part of networks, are presented in different regimes of interest, such as the excitable, regular spiking, and bursting regimes. These results establish a consistent link between single-neuron parameters and resulting network dynamics, and will hopefully be useful as a guide for modeling.

Patterns in inhibitory networks of simple map neurons

We study the dynamics of networks of inhibitory map-based bursting neurons. Linear analysis allows us to understand how the patterns of bursting are determined by network topology and how they depend on the strength of synaptic connections, when inhibition is balanced. Two kinds of patterns are found depending on the symmetry of the network: slow cyclic patterns riding on subthreshold oscillations where almost all neurons contribute bursts in a sparse manner and fast patterns of bursts in which only one of two mutually exclusive groups of neurons take part. We also discuss the properties of the neuron model that underlie the described phenomena, comment on the limitations of the technique of analysis, and point to some possible ways to overcome them.

Sparse repulsive coupling enhances synchronization in complex networks

Through the last years, different strategies to enhance synchronization in complex networks have been proposed. In this work, we show that synchronization of nonidentical dynamical units that are attractively coupled in a small-world network is strongly improved by just making phase-repulsive a tiny fraction of the couplings. By a purely topological analysis that does not depend on the dynamical model, we link the emerging dynamical behavior with the structural properties of the sparsely coupled repulsive network.