Pulsating fronts in periodically modulated neural field models

Does the brain behave like a (complex) network? I. Dynamics

Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network structure does not entail that the brain actually works as a network. Asking whether the brain behaves as a network means asking whether network properties count. From the viewpoint of neurophysiology and, possibly, of brain physics, the most substantial issues a network structure may be instrumental in addressing relate to the influence of network properties on brain dynamics and to whether these properties ultimately explain some aspects of brain function. Here, we address the dynamical implications of complex network, examining which aspects and scales of brain activity may be understood to genuinely behave as a network. To do so, we first define the meaning of networkness, and analyse some of its implications. We then examine ways in which brain anatomy and dynamics can be endowed with a network structure and discuss possible ways in which network structure may be shown to represent a genuine organisational principle of brain activity, rather than just a convenient description of its anatomy and dynamics.

MULTISCALE MOTION AND DEFORMATION OF BUMPS IN STOCHASTIC NEURAL FIELDS WITH DYNAMIC CONNECTIVITY

The distinct timescales of synaptic plasticity and neural activity dynamics play an important role in the brain's learning and memory systems. Activity-dependent plasticity reshapes neural circuit architecture, determining spontaneous and stimulus-encoding spatiotemporal patterns of neural activity. Neural activity bumps maintain short term memories of continuous parameter values, emerging in spatially organized models with short-range excitation and long-range inhibition. Previously, we demonstrated nonlinear Langevin equations derived using an interface method which accurately describe the dynamics of bumps in continuum neural fields with separate excitatory/inhibitory populations. Here we extend this analysis to incorporate effects of short term plasticity that dynamically modifies connectivity described by an integral kernel. Linear stability analysis adapted to these piecewise smooth models with Heaviside firing rates further indicates how plasticity shapes the bumps' local dynamics. Facilitation (depression), which strengthens (weakens) synaptic connectivity originating from active neurons, tends to increase (decrease) stability of bumps when acting on excitatory synapses. The relationship is inverted when plasticity acts on inhibitory synapses. Multiscale approximations of the stochastic dynamics of bumps perturbed by weak noise reveal that the plasticity variables evolve to slowly diffusing and blurred versions of their stationary profiles. Nonlinear Langevin equations associated with bump positions or interfaces coupled to slowly evolving projections of plasticity variables accurately describe how these smoothed synaptic efficacy profiles can tether or repel wandering bumps.

Distinct Excitatory and Inhibitory Bump Wandering in a Stochastic Neural Field

Localized persistent cortical neural activity is a validated neural substrate of parametric working memory. Such activity "bumps" represent the continuous location of a cue over several seconds. Pyramidal (excitatory (E )) and interneuronal (inhibitory (I )) subpopulations exhibit tuned bumps of activity, linking neural dynamics to behavioral inaccuracies observed in memory recall. However, many bump attractor models collapse these subpopulations into a single joint E /I (lateral inhibitory) population and do not consider the role of interpopulation neural architecture and noise correlations. Both factors have a high potential to impinge upon the stochastic dynamics of these bumps, ultimately shaping behavioral response variance. In our study, we consider a neural field model with separate E /I populations and leverage asymptotic analysis to derive a nonlinear Langevin system describing E /I bump interactions. While the E bump attracts the I bump, the I bump stabilizes but can also repel the E bump, which can result in prolonged relaxation dynamics when both bumps are perturbed. Furthermore, the structure of noise correlations within and between subpopulations strongly shapes the variance in bump position. Surprisingly, higher interpopulation correlations reduce variance.

Pattern formation in a 2-population homogenized neuronal network model

We study pattern formation in a 2-population homogenized neural field model of the Hopfield type in one spatial dimension with periodic microstructure. The connectivity functions are periodically modulated in both the synaptic footprint and in the spatial scale. It is shown that the nonlocal synaptic interactions promote a finite band width instability. The stability method relies on a sequence of wave-number dependent invariants of [Formula: see text]-stability matrices representing the sequence of Fourier-transformed linearized evolution equations for the perturbation imposed on the homogeneous background. The generic picture of the instability structure consists of a finite set of well-separated gain bands. In the shallow firing rate regime the nonlinear development of the instability is determined by means of the translational invariant model with connectivity kernels replaced with the corresponding period averaged connectivity functions. In the steep firing rate regime the pattern formation process depends sensitively on the spatial localization of the connectivity kernels: For strongly localized kernels this process is determined by the translational invariant model with period averaged connectivity kernels, whereas in the complementary regime of weak and moderate localization requires the homogenized model as a starting point for the analysis. We follow the development of the instability numerically into the nonlinear regime for both steep and shallow firing rate functions when the connectivity kernels are modeled by means of an exponentially decaying function. We also study the pattern forming process numerically as a function of the heterogeneity parameters in four different regimes ranging from the weakly modulated case to the strongly heterogeneous case. For the weakly modulated regime, we observe that stable spatial oscillations are formed in the steep firing rate regime, whereas we get spatiotemporal oscillations in the shallow regime of the firing rate functions.

Bumps and oscillons in networks of spiking neurons

We study localized patterns in an exact mean-field description of a spatially extended network of quadratic integrate-and-fire neurons. We investigate conditions for the existence and stability of localized solutions, so-called bumps, and give an analytic estimate for the parameter range, where these solutions exist in parameter space, when one or more microscopic network parameters are varied. We develop Galerkin methods for the model equations, which enable numerical bifurcation analysis of stationary and time-periodic spatially extended solutions. We study the emergence of patterns composed of multiple bumps, which are arranged in a snake-and-ladder bifurcation structure if a homogeneous or heterogeneous synaptic kernel is suitably chosen. Furthermore, we examine time-periodic, spatially localized solutions (oscillons) in the presence of external forcing, and in autonomous, recurrently coupled excitatory and inhibitory networks. In both cases, we observe period-doubling cascades leading to chaotic oscillations.

Alternating chimeras in networks of ephaptically coupled bursting neurons

The distinctive phenomenon of the chimera state has been explored in neuronal systems under a variety of different network topologies during the last decade. Nevertheless, in all the works, the neurons are presumed to interact with each other directly with the help of synapses only. But, the influence of ephaptic coupling, particularly magnetic flux across the membrane, is mostly unexplored and should essentially be dealt with during the emergence of collective electrical activities and propagation of signals among the neurons in a network. Through this article, we report the development of an emerging dynamical state, namely, the alternating chimera, in a network of identical neuronal systems induced by an external electromagnetic field. Owing to this interaction scenario, the nonlinear neuronal oscillators are coupled indirectly via electromagnetic induction with magnetic flux, through which neurons communicate in spite of the absence of physical connections among them. The evolution of each neuron, here, is described by the three-dimensional Hindmarsh-Rose dynamics. We demonstrate that the presence of such non-locally and globally interacting external environments induces a stationary alternating chimera pattern in the ensemble of neurons, whereas in the local coupling limit, the network exhibits a transient chimera state whenever the local dynamics of the neurons is of the chaotic square-wave bursting type. For periodic square-wave bursting of the neurons, a similar qualitative phenomenon has been witnessed with the exception of the disappearance of cluster states for non-local and global interactions. Besides these observations, we advance our work while providing confirmation of the findings for neuronal ensembles exhibiting plateau bursting dynamics and also put forward the fact that the plateau pattern actually favors the alternating chimera more than others. These results may deliver better interpretations for different aspects of synchronization appearing in a network of neurons through field coupling that also relaxes the prerequisite of synaptic connectivity for realizing the chimera state in neuronal networks.

Synaptic efficacy shapes resource limitations in working memory

Working memory (WM) is limited in its temporal length and capacity. Classic conceptions of WM capacity assume the system possesses a finite number of slots, but recent evidence suggests WM may be a continuous resource. Resource models typically assume there is no hard upper bound on the number of items that can be stored, but WM fidelity decreases with the number of items. We analyze a neural field model of multi-item WM that associates each item with the location of a bump in a finite spatial domain, considering items that span a one-dimensional continuous feature space. Our analysis relates the neural architecture of the network to accumulated errors and capacity limitations arising during the delay period of a multi-item WM task. Networks with stronger synapses support wider bumps that interact more, whereas networks with weaker synapses support narrower bumps that are more susceptible to noise perturbations. There is an optimal synaptic strength that both limits bump interaction events and the effects of noise perturbations. This optimum shifts to weaker synapses as the number of items stored in the network is increased. Our model not only provides a circuit-based explanation for WM capacity, but also speaks to how capacity relates to the arrangement of stored items in a feature space.

Neural field model of memory-guided search

Many organisms can remember locations they have previously visited during a search. Visual search experiments have shown exploration is guided away from these locations, reducing redundancies in the search path before finding a hidden target. We develop and analyze a two-layer neural field model that encodes positional information during a search task. A position-encoding layer sustains a bump attractor corresponding to the searching agent's current location, and search is modeled by velocity input that propagates the bump. A memory layer sustains persistent activity bounded by a wave front, whose edges expand in response to excitatory input from the position layer. Search can then be biased in response to remembered locations, influencing velocity inputs to the position layer. Asymptotic techniques are used to reduce the dynamics of our model to a low-dimensional system of equations that track the bump position and front boundary. Performance is compared for different target-finding tasks.

Spatiotemporal canards in neural field equations

Canards are special solutions to ordinary differential equations that follow invariant repelling slow manifolds for long time intervals. In realistic biophysical single-cell models, canards are responsible for several complex neural rhythms observed experimentally, but their existence and role in spatially extended systems is largely unexplored. We identify and describe a type of coherent structure in which a spatial pattern displays temporal canard behavior. Using interfacial dynamics and geometric singular perturbation theory, we classify spatiotemporal canards and give conditions for the existence of folded-saddle and folded-node canards. We find that spatiotemporal canards are robust to changes in the synaptic connectivity and firing rate. The theory correctly predicts the existence of spatiotemporal canards with octahedral symmetry in a neural field model posed on the unit sphere.

Neural Field Models with Threshold Noise

The original neural field model of Wilson and Cowan is often interpreted as the averaged behaviour of a network of switch like neural elements with a distribution of switch thresholds, giving rise to the classic sigmoidal population firing-rate function so prevalent in large scale neuronal modelling. In this paper we explore the effects of such threshold noise without recourse to averaging and show that spatial correlations can have a strong effect on the behaviour of waves and patterns in continuum models. Moreover, for a prescribed spatial covariance function we explore the differences in behaviour that can emerge when the underlying stationary distribution is changed from Gaussian to non-Gaussian. For travelling front solutions, in a system with exponentially decaying spatial interactions, we make use of an interface approach to calculate the instantaneous wave speed analytically as a series expansion in the noise strength. From this we find that, for weak noise, the spatially averaged speed depends only on the choice of covariance function and not on the shape of the stationary distribution. For a system with a Mexican-hat spatial connectivity we further find that noise can induce localised bump solutions, and using an interface stability argument show that there can be multiple stable solution branches.

Sensory feedback in a bump attractor model of path integration

Mammalian spatial navigation systems utilize several different sensory information channels. This information is converted into a neural code that represents the animal's current position in space by engaging place cell, grid cell, and head direction cell networks. In particular, sensory landmark (allothetic) cues can be utilized in concert with an animal's knowledge of its own velocity (idiothetic) cues to generate a more accurate representation of position than path integration provides on its own (Battaglia et al. The Journal of Neuroscience 24(19):4541-4550 (2004)). We develop a computational model that merges path integration with feedback from external sensory cues that provide a reliable representation of spatial position along an annular track. Starting with a continuous bump attractor model, we explore the impact of synaptic spatial asymmetry and heterogeneity, which disrupt the position code of the path integration process. We use asymptotic analysis to reduce the bump attractor model to a single scalar equation whose potential represents the impact of asymmetry and heterogeneity. Such imperfections cause errors to build up when the network performs path integration, but these errors can be corrected by an external control signal representing the effects of sensory cues. We demonstrate that there is an optimal strength and decay rate of the control signal when cues appear either periodically or randomly. A similar analysis is performed when errors in path integration arise from dynamic noise fluctuations. Again, there is an optimal strength and decay of discrete control that minimizes the path integration error.

Asymmetry in neural fields: a spatiotemporal encoding mechanism

Neural field models have been successfully applied to model diverse brain mechanisms like visual attention, motor control, and memory. Most theoretical and modeling works have focused on the study of the dynamics of such systems under variations in neural connectivity, mainly symmetric connectivity among neurons. However, less attention has been given to the emerging properties of neuron populations when neural connectivity is asymmetric, although asymmetric activity propagation has been observed in cortical tissue. Here we explore the dynamics of neural fields with asymmetric connectivity and show, in the case of front propagation, that it can bias the population to follow a certain trajectory with higher activation. We find that asymmetry relates linearly to the input speed when the input is spatially localized, and this relation holds for different kernels and input shapes. To illustrate the behavior of asymmetric connectivity, we present an application: standard video sequences of human motion were encoded using the asymmetric neural field and compared to computer vision techniques. Overall, our results indicate that asymmetric neural fields are a competitive approach for spatiotemporal encoding with two main advantages: online classification and distributed operation.

Propagation of CaMKII translocation waves in heterogeneous spiny dendrites

CaMKII (Ca²⁺-calmodulin-dependent protein kinase II) is a key regulator of glutamatergic synapses and plays an essential role in many forms of synaptic plasticity. It has recently been observed experimentally that stimulating a local region of dendrite not only induces the local translocation of CaMKII from the dendritic shaft to synaptic targets within spines, but also initiates a wave of CaMKII translocation that spreads distally through the dendrite with an average speed of order 1 μm/s. We have previously developed a simple reaction-diffusion model of CaMKII translocation waves that can account for the observed wavespeed and predicts wave propagation failure if the density of spines is too high. A major simplification of our previous model was to treat the distribution of spines as spatially uniform. However, there are at least two sources of heterogeneity in the spine distribution that occur on two different spatial scales. First, spines are discrete entities that are joined to a dendritic branch via a thin spine neck of submicron radius, resulting in spatial variations in spine density at the micron level. The second source of heterogeneity occurs on a much longer length scale and reflects the experimental observation that there is a slow proximal to distal variation in the density of spines. In this paper, we analyze how both sources of heterogeneity modulate the speed of CaMKII translocation waves along a spiny dendrite. We adapt methods from the study of the spread of biological invasions in heterogeneous environments, including homogenization theory of pulsating fronts and Hamilton-Jacobi dynamics of sharp interfaces.

Response of traveling waves to transient inputs in neural fields

We analyze the effects of transient stimulation on traveling waves in neural field equations. Neural fields are modeled as integro-differential equations whose convolution term represents the synaptic connections of a spatially extended neuronal network. The adjoint of the linearized wave equation can be used to identify how a particular input will shift the location of a traveling wave. This wave response function is analogous to the phase response curve of limit cycle oscillators. For traveling fronts in an excitatory network, the sign of the shift depends solely on the sign of the transient input. A complementary estimate of the effective shift is derived using an equation for the time-dependent speed of the perturbed front. Traveling pulses are analyzed in an asymmetric lateral inhibitory network and they can be advanced or delayed, depending on the position of spatially localized transient inputs. We also develop bounds on the amplitude of transient input necessary to terminate traveling pulses, based on the global bifurcation structure of the neural field.

The mechanisms for compression and reflection of cortical waves

Waves are common in cortical networks and may be important for carrying information about a stimulus from one local circuit to another. In a recent study of visually evoked waves in rat cortex, compression and reflection of waves are observed as the activation passes from visual areas V1 to V2. The authors of this study apply bicuculline (BMI) and demonstrate that the reflection disappears. They conclude that inhibition plays a major role in compression and reflection. We present several models for propagating waves in heterogeneous media and show that the velocity and thus compression depends weakly on inhibition. We propose that the main site of action of BMI with respect to wave propagation is on the threshold for firing which we suggest is related to action on potassium channels. We combine numerical and analytic methods to explore both compression and reflection in an excitable system with synaptic coupling.

Bifurcations of lurching waves in a thalamic neuronal network

We consider a two-layer, one-dimensional lattice of neurons; one layer consists of excitatory thalamocortical neurons, while the other is comprised of inhibitory reticular thalamic neurons. Such networks are known to support "lurching" waves, for which propagation does not appear smooth, but rather progresses in a saltatory fashion; these waves can be characterized by different spatial widths (different numbers of neurons active at the same time). We show that these lurching waves are fixed points of appropriately defined Poincaré maps, and follow these fixed points as parameters are varied. In this way, we are able to explain observed transitions in behavior, and, in particular, to show how branches with different spatial widths are linked with each other. Our computer-assisted analysis is quite general and could be applied to other spatially extended systems which exhibit this non-trivial form of wave propagation.