A brief history of excitable map-based neurons and neural networks

Optimal input reverberation and homeostatic self-organization toward the edge of synchronization

Transient or partial synchronization can be used to do computations, although a fully synchronized network is sometimes related to the onset of epileptic seizures. Here, we propose a homeostatic mechanism that is capable of maintaining a neuronal network at the edge of a synchronization transition, thereby avoiding the harmful consequences of a fully synchronized network. We model neurons by maps since they are dynamically richer than integrate-and-fire models and more computationally efficient than conductance-based approaches. We first describe the synchronization phase transition of a dense network of neurons with different tonic spiking frequencies coupled by gap junctions. We show that at the transition critical point, inputs optimally reverberate through the network activity through transient synchronization. Then, we introduce a local homeostatic dynamic in the synaptic coupling and show that it produces a robust self-organization toward the edge of this phase transition. We discuss the potential biological consequences of this self-organization process, such as its relation to the Brain Criticality hypothesis, its input processing capacity, and how its malfunction could lead to pathological synchronization and the onset of seizure-like activity.

Self-tuned criticality: Controlling a neuron near its bifurcation point via temporal correlations

Previous work showed that the collective activity of large neuronal networks can be tamed to remain near its critical point by a feedback control that maximizes the temporal correlations of the mean-field fluctuations. Since such correlations behave similarly near instabilities across nonlinear dynamical systems, it is expected that the principle should control also low-dimensional dynamical systems exhibiting continuous or discontinuous bifurcations from fixed points to limit cycles. Here we present numerical evidence that the dynamics of a single neuron can be controlled in the vicinity of its bifurcation point. The approach is tested in two models: a two-dimensional generic excitable map and the paradigmatic FitzHugh-Nagumo neuron model. The results show that in both cases, the system can be self-tuned to its bifurcation point by modifying the control parameter according to the first coefficient of the autocorrelation function.

Variations of the spontaneous electrical activities of the neuronal networks imposed by the exposure of electromagnetic radiations using computational map-based modeling

The interaction between neurons in a neuronal network develops spontaneous electrical activities. But the effects of electromagnetic radiation on these activities have not yet been well explored. In this study, a ring of three coupled 1-dimensional Rulkov neurons and the generated electromagnetic field (EMF) are considered to investigate how the spontaneous activities might change regarding the EMF exposure. By employing the bifurcation analysis and time series, a comprehensive view of neuronal behavioral changes due to electromagnetic inductions is provided. The main findings of this study are as follows: 1) When a neuronal network is showing a spontaneous chaotic firing manner (without any external stimuli), a generated magnetic field inhibits this type of behavior. In fact, EMF completely eliminated the chaotic intrinsic behaviors of the neuronal loop. 2) When the network is exhibiting regular period-3 spiking patterns, the generated magnetic field changes its firing pattern to chaotic spiking, which is similar to epileptic seizures. 3) With weak synaptic connections, electromagnetic radiation inhibits and suppresses neuronal activities. 4) If the external magnetic flux has a high amplitude, it can change the shape of the induction current according to its shape 5) when there are weak synaptic connections in the network, a high-frequency external magnetic flux engenders high-frequency fluctuations in the membrane voltages. On the whole, electromagnetic radiation changes the pattern of the spontaneous activities of neuronal networks in the brain according to synaptic strengths and initial states of the neurons.

Spiral waves in a hybrid discrete excitable media with electromagnetic flux coupling

Though there are many neuron models based on differential equations, the complexity in realizing them into digital circuits is still a challenge. Hence, many new discrete neuron models have been recently proposed, which can be easily implemented in digital circuits. We consider the well-known FitzHugh-Nagumo model and derive the discrete version of the model considering the sigmoid type of recovery variable and electromagnetic flux coupling. We show the various time series plots confirming the existence of periodic and chaotic bursting as in differential equation type neuron models. Also, we have used the bifurcation plots, Lyapunov exponents, and frequency bifurcations to investigate the dynamics of the proposed discrete neuron model. Different topologies of networks like single, two, and three layers are considered to analyze the wave propagation phenomenon in the network. We introduce the concept of using energy levels of nodes to study the spiral wave existence and compare them with the spatiotemporal snapshots. Interestingly, the energy plots clearly show that when the energy level of nodes is different and distributed, the occurrence of the spiral waves is identified in the network.

Altered communication dynamics reflect cognitive deficits in temporal lobe epilepsy

OBJECTIVE

Although temporal lobe epilepsy (TLE) is recognized as a system-level disorder, little work has investigated pathoconnectomics from a dynamic perspective. By leveraging computational simulations that quantify patterns of information flow across the connectome, we tested the hypothesis that network communication is abnormal in this condition, studied the interplay between hippocampal- and network-level disease effects, and assessed associations with cognition.

METHODS

We simulated signal spreading via a linear threshold model that temporally evolves on a structural graph derived from diffusion-weighted magnetic resonance imaging (MRI), comparing a homogeneous group of 31 patients with histologically proven hippocampal sclerosis to 31 age- and sex-matched healthy controls. We evaluated the modulatory effects of structural alterations of the neocortex and hippocampus on network dynamics. Furthermore, multivariate statistics addressed the relationship with cognitive parameters.

RESULTS

We observed a slowing of in- and out-spreading times across multiple areas bilaterally, indexing delayed information flow, with the strongest effects in ipsilateral frontotemporal regions, thalamus, and hippocampus. Effects were markedly reduced when controlling for hippocampal volume but not cortical thickness, underscoring the central role of the hippocampus in whole-brain disease expression. Multivariate analysis associated slower spreading time in frontoparietal, limbic, default mode, and subcortical networks with impairment across tasks tapping into sensorimotor, executive, memory, and verbal abilities.

SIGNIFICANCE

Moving beyond descriptions of static topology toward the formulation of brain dynamics, our work provides novel insight into structurally mediated network dysfunction and demonstrates that altered whole-brain communication dynamics contribute to common cognitive difficulties in TLE.

Spatially resolved dendritic integration: towards a functional classification of neurons

The vast tree-like dendritic structure of neurons allows them to receive and integrate input from many neurons. A wide variety of neuronal morphologies exist, however, their role in dendritic integration, and how it shapes the response of the neuron, is not yet fully understood. Here, we study the evolution and interactions of dendritic spikes in excitable neurons with complex real branch structures. We focus on dozens of digitally reconstructed illustrative neurons from the online repository NeuroMorpho.org, which contains over 130,000 neurons. Yet, our methods can be promptly extended to any other neuron. This approach allows us to estimate and map specific and heterogeneous patterns of activity observed across extensive dendritic trees with thousands of compartments. We propose a classification of neurons based on the location of the soma (centrality) and the number of branches connected to the soma. These are key topological factors in determining the neuron's energy consumption, firing rate, and the dynamic range, which quantifies the range in synaptic input rate that can be reliably encoded by the neuron's firing rate. Moreover, we find that bifurcations, the structural building blocks of complex dendrites, play a major role in increasing the dynamic range of neurons. Our results provide a better understanding of the effects of neuronal morphology in the diversity of neuronal dynamics and function.

Discrete Langevin machine: Bridging the gap between thermodynamic and neuromorphic systems

A formulation of Langevin dynamics for discrete systems is derived as a class of generic stochastic processes. The dynamics simplify for a two-state system and suggest a network architecture which is implemented by the Langevin machine. The Langevin machine represents a promising approach to compute successfully quantitative exact results of Boltzmann distributed systems by LIF neurons. Besides a detailed introduction of the dynamics, different simplified models of a neuromorphic hardware system are studied with respect to a control of emerging sources of errors.

Proper spatial heterogeneities expand the regime of scale-free behavior in a lattice of excitable elements

Signatures of criticality, such as power law scaling of observables, have been empirically found in a plethora of real-life settings, including biological systems. The presence of critical states is believed to have many functional advantages and is associated with optimal operational abilities. Typically, critical dynamics arises in the proximity of phase transition points between absorbing disordered states (subcriticality) and ordered active regimes (supercriticality) and requires a high degree of fine tuning to emerge, which is unlikely to occur in real biological systems. In the present study we propose a rather simple, and biologically relevant mechanism that profoundly expands the critical-like region. In particular, by means of numerical simulation we show that incorporating spatial heterogeneities into the square lattice of map-based excitable oscillators broadens the parameter space in which the distribution of excitation wave sizes follows closely a power law. Most importantly, this behavior is only observed if the spatial profile exhibits intermediate-sized patches with similar excitability levels, whereas for large and small spatial clusters only marginal widening of the critical state is detected. Furthermore, it turned out that the presence of spatial disorder in general amplifies the size of excitation waves, whereby the relatively highest contributions are observed in the proximity of the critical point. We argue that the reported mechanism is of particular importance for excitable systems with local interactions between individual elements.

Unidirectional Synchronization of Hodgkin-Huxley Neurons With Prescribed Performance Under Transcranial Magneto-Acoustical Simulation

This paper exploits the unidirectional synchronization dynamics of two Hodgkin-Huxley (HH) neurons under transcranial magneto-acoustical stimulation (TMAS). The major purpose is to explore a control scheme to make the spiking modes of the neural potentials stimulated by TMAS achieve synchronization states under the feedback input. For this purpose, an adaptive neural controller, which makes the neurons satisfy the prescribed master-slaver synchronization performance, is designed by introducing a tracking error into Lyapunov analysis. Under the proposed control scheme, the slaver neuron can not only overcome the model uncertainties and the difficulties brought by prescribed performance, but also track the spiking patterns of the master neuron. Finally, the simulations are implemented to demonstrate the effectiveness of the proposed controller, that is, the TMAS induced synchronization states of the HH neuron system can achieve the prescribed performance under the proposed controller.

Observability and synchronization of neuron models

Observability is the property that enables recovering the state of a dynamical system from a reduced number of measured variables. In high-dimensional systems, it is therefore important to make sure that the variable recorded to perform the analysis conveys good observability of the system dynamics. The observability of a network of neuron models depends nontrivially on the observability of the node dynamics and on the topology of the network. The aim of this paper is twofold. First, to perform a study of observability using four well-known neuron models by computing three different observability coefficients. This not only clarifies observability properties of the models but also shows the limitations of applicability of each type of coefficients in the context of such models. Second, to study the emergence of phase synchronization in networks composed of neuron models. This is done performing multivariate singular spectrum analysis which, to the best of the authors' knowledge, has not been used in the context of networks of neuron models. It is shown that it is possible to detect phase synchronization: (i) without having to measure all the state variables, but only one (that provides greatest observability) from each node and (ii) without having to estimate the phase.

A Mathematical Model for Storage and Recall of Images using Targeted Synchronization of Coupled Maps

We propose a mathematical model for storage and recall of images using coupled maps. We start by theoretically investigating targeted synchronization in coupled map systems wherein only a desired (partial) subset of the maps is made to synchronize. A simple method is introduced to specify coupling coefficients such that targeted synchronization is ensured. The principle of this method is extended to storage/recall of images using coupled Rulkov maps. The process of adjusting coupling coefficients between Rulkov maps (often used to model neurons) for the purpose of storing a desired image mimics the process of adjusting synaptic strengths between neurons to store memories. Our method uses both synchronisation and synaptic weight modification, as the human brain is thought to do. The stored image can be recalled by providing an initial random pattern to the dynamical system. The storage and recall of the standard image of Lena is explicitly demonstrated.

Transient sequences in a hypernetwork generated by an adaptive network of spiking neurons

We propose a model of an adaptive network of spiking neurons that gives rise to a hypernetwork of its dynamic states at the upper level of description. Left to itself, the network exhibits a sequence of transient clustering which relates to a traffic in the hypernetwork in the form of a random walk. Receiving inputs the system is able to generate reproducible sequences corresponding to stimulus-specific paths in the hypernetwork. We illustrate these basic notions by a simple network of discrete-time spiking neurons together with its FPGA realization and analyse their properties.This article is part of the themed issue 'Mathematical methods in medicine: neuroscience, cardiology and pathology'.

Phase diagrams and dynamics of a computationally efficient map-based neuron model

We introduce a new map-based neuron model derived from the dynamical perceptron family that has the best compromise between computational efficiency, analytical tractability, reduced parameter space and many dynamical behaviors. We calculate bifurcation and phase diagrams analytically and computationally that underpins a rich repertoire of autonomous and excitable dynamical behaviors. We report the existence of a new regime of cardiac spikes corresponding to nonchaotic aperiodic behavior. We compare the features of our model to standard neuron models currently available in the literature.

A modular architecture for transparent computation in recurrent neural networks

Computation is classically studied in terms of automata, formal languages and algorithms; yet, the relation between neural dynamics and symbolic representations and operations is still unclear in traditional eliminative connectionism. Therefore, we suggest a unique perspective on this central issue, to which we would like to refer as transparent connectionism, by proposing accounts of how symbolic computation can be implemented in neural substrates. In this study we first introduce a new model of dynamics on a symbolic space, the versatile shift, showing that it supports the real-time simulation of a range of automata. We then show that the Gödelization of versatile shifts defines nonlinear dynamical automata, dynamical systems evolving on a vectorial space. Finally, we present a mapping between nonlinear dynamical automata and recurrent artificial neural networks. The mapping defines an architecture characterized by its granular modularity, where data, symbolic operations and their control are not only distinguishable in activation space, but also spatially localizable in the network itself, while maintaining a distributed encoding of symbolic representations. The resulting networks simulate automata in real-time and are programmed directly, in the absence of network training. To discuss the unique characteristics of the architecture and their consequences, we present two examples: (i) the design of a Central Pattern Generator from a finite-state locomotive controller, and (ii) the creation of a network simulating a system of interactive automata that supports the parsing of garden-path sentences as investigated in psycholinguistics experiments.

Inverse period-doubling bifurcations determine complex structure of bursting in a one-dimensional non-autonomous map

We propose a simple one-dimensional non-autonomous map, in which some novel bursting patterns (e.g., "fold/double inverse flip" bursting, "fold/multiple inverse flip" bursting, and "fold/a cascade of inverse flip" bursting) can be observed. Typically, these bursting patterns exhibit complex structures containing a chain of inverse period-doubling bifurcations. The active states related to these bursting can be period-2(n) (n = 1, 2, 3,…) attractors or chaotic attractors, which may evolve to quiescence by a chain of inverse period-doubling bifurcations when the slow excitation decreases through period-doubling bifurcation points of the map. This accounts for the complex inverse period-doubling bifurcation structures observed in bursting patterns. Our findings enrich the possible routes to bursting as well as the underlying mechanisms of bursting.

The transfer function of neuron spike

The mathematical modeling of neuronal signals is a relevant problem in neuroscience. The complexity of the neuron behavior, however, makes this problem a particularly difficult task. Here, we propose a discrete-time linear time-invariant (LTI) model with a rational function in order to represent the neuronal spike detected by an electrode located in the surroundings of the nerve cell. The model is presented as a cascade association of two subsystems: one that generates an action potential from an input stimulus, and one that represents the medium between the cell and the electrode. The suggested approach employs system identification and signal processing concepts, and is dissociated from any considerations about the biophysical processes of the neuronal cell, providing a low-complexity alternative to model the neuronal spike. The model is validated by using in vivo experimental readings of intracellular and extracellular signals. A computational simulation of the model is presented in order to assess its proximity to the neuronal signal and to observe the variability of the estimated parameters. The implications of the results are discussed in the context of spike sorting.

Towards a complex system understanding of bipolar disorder: A map based model of a complex winnerless competition

Bipolar disorder is characterized by repeated erratic episodes of mania and depression, which can be understood as pathological complex system behavior involving cognitive, affective and psychomotor disturbance. In order to illuminate dynamical aspects of the longitudinal course of the illness, we propose here a novel complex model based on the notion of competition between recurrent maps, which mathematically represent the dynamics of activation in excitatory (Glutamatergic) and inhibitory (GABAergic) pathways. We assume that manic and depressive states can be considered stable sub attractors of a dynamical system through which the mood trajectory moves. The model provides a theoretical framework which can account for a number of complex phenomena of bipolar disorder, including intermittent transition between the two poles of the disorder, rapid and ultra-rapid cycling of episodes and manicogenic effects of antidepressants.

Critical avalanches and subsampling in map-based neural networks coupled with noisy synapses

Many different kinds of noise are experimentally observed in the brain. Among them, we study a model of noisy chemical synapse and obtain critical avalanches for the spatiotemporal activity of the neural network. Neurons and synapses are modeled by dynamical maps. We discuss the relevant neuronal and synaptic properties to achieve the critical state. We verify that networks of functionally excitable neurons with fast synapses present power-law avalanches, due to rebound spiking dynamics. We also discuss the measuring of neuronal avalanches by subsampling our data, shedding light on the experimental search for self-organized criticality in neural networks.