Search

We present some existence and localization results for periodic solutions of impulsive first‐order coupled nonlinear systems of two equations, without requiring periodicity for the nonlinearities. The arguments are based on Schauder's fixed point theorem together with the upper and lower solution methods, where the upper and lower solutions are not necessarily well ordered. In addition, results on equi‐regulated functions are required for the impulsive analysis. An application to a Wilson‐Cowan system of two strongly coupled neurons illustrates one of the main results. [ABSTRACT FROM AUTHOR]

This work presents a fractional-order Wilson-Cowan model derivation under Caputo's formalism, considering an order of 0 < α ≤ 1 . To that end, we propose memory-dependent response functions and average neuronal excitation functions that permit us to naturally arrive at a fractional-order model that incorporates past dynamics into the description of synaptically coupled neuronal populations' activity. We then shift our focus on a particular example, aiming to analyze the fractional-order dynamics of the disinhibited cortex. This system mimics cortical activity observed during neurological disorders such as epileptic seizures, where an imbalance between excitation and inhibition is present, which allows brain dynamics to transition to a hyperexcited activity state. In the context of the first-order mathematical model, we recover traditional results showing a transition from a low-level activity state to a potentially pathological high-level activity state as an external factor modifies cortical inhibition. On the other hand, under the fractional-order formulation, we establish novel results showing that the system resists such transition as the order is decreased, permitting the possibility of staying in the low-activity state even with increased disinhibition. Furthermore, considering the memory index interpretation of the fractional-order model motivation here developed, our results establish that by increasing the memory index, the system becomes more resistant to transitioning towards the high-level activity state. That is, one possible effect of the memory index is to stabilize neuronal activity. Noticeably, this neuronal stabilizing effect is similar to homeostatic plasticity mechanisms. To summarize our results, we present a two-parameter structural portrait describing the system's dynamics dependent on a proposed disinhibition parameter and the order. We also explore numerical model simulations to validate our results. [ABSTRACT FROM AUTHOR]

Neural population modeling, including the role of neural attractors, is a promising tool for understanding many aspects of brain function. We propose a modeling framework to connect the abstract variables used in modeling to recent cellular-level estimates of the bioenergetic costs of different aspects of neural activity, measured in ATP consumed per second per neuron. Based on recent work, an empirical reference for brain ATP use for the awake resting brain was estimated as ~2 x 109 ATP/s-neuron across several mammalian species. The energetics framework was applied to the Wilson-Cowan (WC) model of two interacting populations of neurons, one excitatory (E) and one inhibitory (I). Attractors were considered to exhibit steady-state behavior and limit cycle behavior, both of which end when the excitatory stimulus ends, and sustained activity that persists after the stimulus ends. The energy cost of limit cycles, with oscillations much faster than the average neuronal firing rate of the population, is tracked more closely with the firing rate than the limit cycle frequency. Self-sustained firing driven by recurrent excitation, though, involves higher firing rates and a higher energy cost. As an example of a simple network in which each node is a WC model, a combination of three nodes can serve as a flexible circuit element that turns on with an oscillating output when input passes a threshold and then persists after the input ends (an "on-switch"), with moderate overall ATP use. The proposed framework can serve as a guide for anchoring neural population models to plausible bioenergetics requirements. [ABSTRACT FROM AUTHOR]

Objectives: The exact mechanisms of deep brain stimulation (DBS) are still an active area of investigation, in spite of its clinical successes. This is due in part to the lack of understanding of the effects of stimulation on neuronal rhythms. Entrainment of brain oscillations has been hypothesised as a potential mechanism of neuromodulation. A better understanding of entrainment might further inform existing methods of continuous DBS, and help refine algorithms for adaptive methods. The purpose of this study is to develop and test a theoretical framework to predict entrainment of cortical rhythms to DBS across a wide range of stimulation parameters. Materials and Methods: We fit a model of interacting neural populations to selected features characterising PD patients' off-stimulation finely-tuned gamma rhythm recorded through electrocorticography. Using the fitted models, we predict basal ganglia DBS parameters that would result in 1:2 entrainment, a special case of sub-harmonic entrainment observed in patients and predicted by theory. Results: We show that the neural circuit models fitted to patient data exhibit 1:2 entrainment when stimulation is provided across a range of stimulation parameters. Furthermore, we verify key features of the region of 1:2 entrainment in the stimulation frequency/amplitude space with follow-up recordings from the same patients, such as the loss of 1:2 entrainment above certain stimulation amplitudes. Conclusion: Our results reveal that continuous, constant frequency DBS in patients may lead to nonlinear patterns of neuronal entrainment across stimulation parameters, and that these responses can be predicted by modelling. Should entrainment prove to be an important mechanism of therapeutic stimulation, our modelling framework may reduce the parameter space that clinicians must consider when programming devices for optimal benefit.

The exact mechanisms of deep brain stimulation (DBS) are still an active area of investigation, in spite of its clinical successes. This is due in part to the lack of understanding of the effects of stimulation on neuronal rhythms. Entrainment of brain oscillations has been hypothesised as a potential mechanism of neuromodulation. A better understanding of entrainment might further inform existing methods of continuous DBS, and help refine algorithms for adaptive methods. The purpose of this study is to develop and test a theoretical framework to predict entrainment of cortical rhythms to DBS across a wide range of stimulation parameters. We fit a model of interacting neural populations to selected features characterising PD patients' off-stimulation finely-tuned gamma rhythm recorded through electrocorticography. Using the fitted models, we predict basal ganglia DBS parameters that would result in 1:2 entrainment, a special case of sub-harmonic entrainment observed in patients and predicted by theory. We show that the neural circuit models fitted to patient data exhibit 1:2 entrainment when stimulation is provided across a range of stimulation parameters. Furthermore, we verify key features of the region of 1:2 entrainment in the stimulation frequency/amplitude space with follow-up recordings from the same patients, such as the loss of 1:2 entrainment above certain stimulation amplitudes. Our results reveal that continuous, constant frequency DBS in patients may lead to nonlinear patterns of neuronal entrainment across stimulation parameters, and that these responses can be predicted by modelling. Should entrainment prove to be an important mechanism of therapeutic stimulation, our modelling framework may reduce the parameter space that clinicians must consider when programming devices for optimal benefit. • Entrainment of brain rhythms at half stimulation frequency is predicted by theory. • Models fitted to off-stimulation patient data predict 1:2 synchronisation regions. • Observations of 1:2 entrainment in patients are in line with model predictions. • 1:2 entrainment is biased to low stimulation frequencies and lost at high amplitudes. [ABSTRACT FROM AUTHOR]

This work aims to study the interplay between the Wilson–Cowan model and connection matrices. These matrices describe cortical neural wiring, while Wilson–Cowan equations provide a dynamical description of neural interaction. We formulate Wilson–Cowan equations on locally compact Abelian groups. We show that the Cauchy problem is well posed. We then select a type of group that allows us to incorporate the experimental information provided by the connection matrices. We argue that the classical Wilson–Cowan model is incompatible with the small-world property. A necessary condition to have this property is that the Wilson–Cowan equations be formulated on a compact group. We propose a p-adic version of the Wilson–Cowan model, a hierarchical version in which the neurons are organized into an infinite rooted tree. We present several numerical simulations showing that the p-adic version matches the predictions of the classical version in relevant experiments. The p-adic version allows the incorporation of the connection matrices into the Wilson–Cowan model. We present several numerical simulations using a neural network model that incorporates a p-adic approximation of the connection matrix of the cat cortex. [ABSTRACT FROM AUTHOR]

Frequency-dependent plasticity refers to changes in synaptic strength in response to different stimulation frequencies. Resonance is a factor known to be of importance in such frequency dependence, however, the role of neural noise in the process remains elusive. Considering the brain is an inherently noisy system, understanding its effects may prove beneficial in shaping therapeutic interventions based on non-invasive brain stimulation protocols. The Wilson-Cowan (WC) model is a well-established model to describe the average dynamics of neural populations and has been shown to exhibit bistability in the presence of noise. However, the important question of how the different stable regimes in the WC model can affect synaptic plasticity when cortical populations interact has not yet been addressed. Therefore, we investigated plasticity dynamics in a WC-based model of interacting neural populations coupled with activity-dependent synapses in which a periodic stimulation was applied in the presence of noise of controlled intensity. The results indicate that for a narrow range of the noise variance, synaptic strength can be optimized. In particular, there is a regime of noise intensity for which synaptic strength presents a triple-stable state. Regulating noise intensity affects the probability that the system chooses one of the stable states, thereby controlling plasticity. These results suggest that noise is a highly influential factor in determining the outcome of plasticity induced by stimulation.

Fifty years ago, Wilson and Cowan developed a mathematical model to describe the activity of neural populations. In this seminal work, they divided the cells in three groups: active, sensitive and refractory, and obtained a dynamical system to describe the evolution of the average firing rates of the populations. In the present work, we investigate the impact of the often neglected refractory state and show that taking it into account can introduce new dynamics. Starting from a continuous-time Markov chain, we perform a rigorous derivation of a mean-field model that includes the refractory fractions of populations as dynamical variables. Then, we perform bifurcation analysis to explain the occurrence of periodic solutions in cases where the classical Wilson–Cowan does not predict oscillations. We also show that our mean-field model is able to predict chaotic behavior in the dynamics of networks with as little as two populations. [ABSTRACT FROM AUTHOR]

This paper attempts to explore and compare the regulatory mechanisms of optogenetic stimulation (OS), deep brain stimulation (DBS) and electromagnetic induction on epilepsy. Based on the Wilson–Cowan model, we first demonstrate that the external input received by excitatory and inhibitory neural populations can induce rich dynamic bifurcation behaviors such as Hopf bifurcation, and make the system exhibit epileptic and normal states. Then, both OS and DBS are shown to be effective in controlling the epileptic state to a normal low-level state, and the stimulus parameters have a broad effective range. However, electromagnetic induction cannot directly control epilepsy to this desired state, even if it can significantly reduce the oscillation frequency of neural populations. One main difference worth noting is that the high spatiotemporal specificity of OS allows it to target inhibitory neuronal populations, whereas DBS and electromagnetic induction can only stimulate excitatory as well as inhibitory neuronal populations together. Next, the propagation behavior of epilepsy is explored under a typical three-node feedback loop structure. An increase in coupling strength accelerates and exacerbates epileptic activity in other brain regions. Finally, OS and DBS applied to the epileptic focus play similar positive roles in controlling the behavior of the area of seizure propagation, while electromagnetic induction still only achieves unsatisfactory effects. It is hoped that these dynamical results can provide insights into the treatment of epilepsy as well as other neurological disorders. [ABSTRACT FROM AUTHOR]

This work aims to study the interplay between the Wilson–Cowan model and connection matrices. These matrices describe cortical neural wiring, while Wilson–Cowan equations provide a dynamical description of neural interaction. We formulate Wilson–Cowan equations on locally compact Abelian groups. We show that the Cauchy problem is well posed. We then select a type of group that allows us to incorporate the experimental information provided by the connection matrices. We argue that the classical Wilson–Cowan model is incompatible with the small-world property. A necessary condition to have this property is that the Wilson–Cowan equations be formulated on a compact group. We propose a p-adic version of the Wilson–Cowan model, a hierarchical version in which the neurons are organized into an infinite rooted tree. We present several numerical simulations showing that the p-adic version matches the predictions of the classical version in relevant experiments. The p-adic version allows the incorporation of the connection matrices into the Wilson–Cowan model. We present several numerical simulations using a neural network model that incorporates a p-adic approximation of the connection matrix of the cat cortex.

In weakly coupled neural oscillator networks describing brain dynamics, the coupling delay is often distributed. We present a theoretical framework to calculate the phase response curve of distributed-delay induced limit cycles with infinite-dimensional phase space. Extending previous works, in which non-delayed or discrete-delay systems were investigated, we develop analytical results for phase response curves of oscillatory systems with distributed delay using Gaussian and log-normal delay distributions. We determine the scalar product and normalization condition for the linearized adjoint of the system required for the calculation of the phase response curve. As a paradigmatic example, we apply our technique to the Wilson–Cowan oscillator model of excitatory and inhibitory neuronal populations under the two delay distributions. We calculate and compare the phase response curves for the Gaussian and log-normal delay distributions. The phase response curves obtained from our adjoint calculations match those compiled by the direct perturbation method, thereby proving that the theory of weakly coupled oscillators can be applied successfully for distributed-delay-induced limit cycles. We further use the obtained phase response curves to derive phase interaction functions and determine the possible phase locked states of multiple inter-coupled populations to illuminate different synchronization scenarios. In numerical simulations, we show that the coupling delay distribution can impact the stability of the synchronization between inter-coupled gamma-oscillatory networks. [ABSTRACT FROM AUTHOR]

The Wilson-Cowan model can emulate gamma oscillations, and thus is extensively used to research the generation of gamma oscillations closely related to cognitive functions. Previous studies have revealed that excitatory and inhibitory inputs to the model can modulate its gamma oscillations. Inhibitory and excitatory self-feedback loops are important structural features of the model, however, its functional role in the regulation of gamma oscillations in the model is still unclear. In the present study, bifurcation analysis and spectrum analysis are employed to elucidate the regulating mechanism of gamma oscillations underlined by the inhibitory and excitatory self-feedback loops, especially how the two self-feedback loops cooperate to generate the gamma oscillations and regulate the oscillation frequency. The present results reveal that, on one hand, the inhibitory self-feedback loop is not conducive to the generation of gamma oscillations, and increased inhibitory self-feedback strength facilitates the enhancement of the oscillation frequency. On the other hand, the excitatory self-feedback loop promotes the generation of gamma oscillations, and increased excitatory self-feedback strength leads to the decrease of oscillation frequency. Finally, theoretical analysis is conducted to provide explain on how the two self-feedback loops play a crucial role in the generation and regulation of neural oscillations in the model. To sum up, Inhibitory and excitatory self-feedback loops play a complementary role in generating and regulating the gamma oscillation in Wilson-Cowan model, and cooperate to bidirectionally regulate the gamma-oscillation frequency in a more flexible manner. These results might provide testable hypotheses for future experimental research. [ABSTRACT FROM AUTHOR]

The Wilson-Cowan model can emulate gamma oscillations, and thus is extensively used to research the generation of gamma oscillations closely related to cognitive functions. Previous studies have revealed that excitatory and inhibitory inputs to the model can modulate its gamma oscillations. Inhibitory and excitatory self-feedback loops are important structural features of the model, however, its functional role in the regulation of gamma oscillations in the model is still unclear. In the present study, bifurcation analysis and spectrum analysis are employed to elucidate the regulating mechanism of gamma oscillations underlined by the inhibitory and excitatory self-feedback loops, especially how the two self-feedback loops cooperate to generate the gamma oscillations and regulate the oscillation frequency. The present results reveal that, on one hand, the inhibitory self-feedback loop is not conducive to the generation of gamma oscillations, and increased inhibitory self-feedback strength facilitates the enhancement of the oscillation frequency. On the other hand, the excitatory self-feedback loop promotes the generation of gamma oscillations, and increased excitatory self-feedback strength leads to the decrease of oscillation frequency. Finally, theoretical analysis is conducted to provide explain on how the two self-feedback loops play a crucial role in the generation and regulation of neural oscillations in the model. To sum up, Inhibitory and excitatory self-feedback loops play a complementary role in generating and regulating the gamma oscillation in Wilson-Cowan model, and cooperate to bidirectionally regulate the gamma-oscillation frequency in a more flexible manner. These results might provide testable hypotheses for future experimental research.

This paper attempts to explore and compare the regulatory mechanisms of optogenetic stimulation (OS), deep brain stimulation (DBS) and electromagnetic induction on epilepsy. Based on the Wilson–Cowan model, we first demonstrate that the external input received by excitatory and inhibitory neural populations can induce rich dynamic bifurcation behaviors such as Hopf bifurcation, and make the system exhibit epileptic and normal states. Then, both OS and DBS are shown to be effective in controlling the epileptic state to a normal low-level state, and the stimulus parameters have a broad effective range. However, electromagnetic induction cannot directly control epilepsy to this desired state, even if it can significantly reduce the oscillation frequency of neural populations. One main difference worth noting is that the high spatiotemporal specificity of OS allows it to target inhibitory neuronal populations, whereas DBS and electromagnetic induction can only stimulate excitatory as well as inhibitory neuronal populations together. Next, the propagation behavior of epilepsy is explored under a typical three-node feedback loop structure. An increase in coupling strength accelerates and exacerbates epileptic activity in other brain regions. Finally, OS and DBS applied to the epileptic focus play similar positive roles in controlling the behavior of the area of seizure propagation, while electromagnetic induction still only achieves unsatisfactory effects. It is hoped that these dynamical results can provide insights into the treatment of epilepsy as well as other neurological disorders.

• We characterise effects of parameter heterogeneity in a neuromorphic directed chain. • Long-range links trigger the emergence of synchronisation and phase coherence. • Local topological condition leads to strong amplification even with backward links. Neuromorphic networks can be described in terms of coarse-grained variables, where emergent sustained behaviours spontaneously arise if stochasticity is properly taken into account. For example it has been recently found that a directed linear chain of connected patch of neurons amplifies an input signal, also tuning its characteristic frequency. Here we study a generalization of such a simple model, introducing heterogeneity and variability in the parameter space and long-range interactions, breaking, in turn, the preferential direction of information transmission of a directed chain. On one hand, enlarging the region of parameters leads to a more complex state space that we analytically characterise; moreover, we explicitly link the strength distribution of the non-local interactions with the frequency distribution of the network oscillations. On the other hand, we found that adding long-range interactions can cause the onset of novel phenomena, as coherent and synchronous oscillations among all the interacting units, which can also coexist with the amplification of the signal. [ABSTRACT FROM AUTHOR]

Neocortical circuits exhibit a rich dynamic repertoire, and their ability to achieve entrainment (adjustment of their frequency to match the input frequency) is thought to support many cognitive functions and indicate functional flexibility. Although previous studies have explored the influence of various circuit properties on this phenomenon, the role of divisive gain modulation (or divisive inhibition) is unknown. This gain control mechanism is thought to be delivered mainly by the soma-targeting interneurons in neocortical microcircuits. In this study, we use a neural mass model of the neocortical microcircuit (extended Wilson–Cowan model) featuring both soma-targeting and dendritetargeting interneuronal subpopulations to investigate the role of divisive gain modulation in entrainment. Our results demonstrate that the presence of divisive inhibition in the microcircuit, as delivered by the soma-targeting interneurons, enables its entrainment to a wider range of input frequencies. Divisive inhibition also promotes a faster entrainment, with the microcircuit needing less time to converge to the fully entrained state. We suggest that divisive inhibition, working alongside subtractive inhibition, allows for more adaptive oscillatory responses in neocortical circuits and, thus, supports healthy brain functioning. NEW & NOTEWORTHY We introduce a computational neocortical microcircuit model that features two inhibitory neural populations, with one providing subtractive and the other divisive inhibition to the excitatory population. We demonstrate that divisive inhibition widens the range of input frequencies to which the microcircuit can become entrained and diminishes the time needed to reach full entrainment. We suggest that divisive inhibition enables more adaptive oscillatory activity, with important implications for both normal and pathological brain function. [ABSTRACT FROM AUTHOR]

• Feasibility of conducting cognitive experiments during brain gamma oscillations. • Optogenetic stimulation enhances low gamma power in the visual cortex. • Visual stimulation has differing effects on gamma power in different layers of the visual cortex. Studying brain functions and activity during gamma oscillations can be a challenge because it requires careful planning to create the necessary conditions for a controlled experiment. Such an experiment consists of placing the brain into a gamma state and investigating cognitive processing with a careful design. Cortical oscillations in the gamma frequency range (30–80 Hz) play an essential role in a variety of cognitive processes, including visual processing and cognition. The present study aims to investigate the effects of a visual stimulus on the primary visual cortex under gamma oscillations. Specifically, we sought to explore the behavior of gamma oscillations triggered by optogenetic stimulation in the II and IV layers of the visual cortex, both with and without concurrent visual stimulation. Our results show that optogenetic stimulation increases the power of gamma oscillation in both layers of the visual cortex. However, the combined stimuli resulted in a reduction of gamma power in layer II and an increase and reinforcement in gamma power in layer IV. Modelling the results with the Wilson-Cowan model suggests changes in the input of the excitatory population due to the combined stimuli. In addition, our analysis of the data using the Lempel-Ziv complexity method supports our interpretations from the modeling. Thus, our results suggest that optogenetic stimulation enhances low gamma power in both layers of the visual cortex, while simultaneous visual stimulation has differing effects on the two layers, reducing gamma power in layer II and increasing it in layer IV. [ABSTRACT FROM AUTHOR]

Gamma oscillations are readily observed in a variety of brain regions during both waking and sleeping states. Computational models of gamma oscillations typically involve simulations of large networks of synaptically coupled spiking units. These networks can exhibit strongly synchronized gamma behavior, whereby neurons fire in near synchrony on every cycle, or weakly modulated gamma behavior, corresponding to stochastic, sparse firing of the individual units on each cycle of the population gamma rhythm. These spiking models offer valuable biophysical descriptions of gamma oscillations; however, because they involve many individual neuronal units they are limited in their ability to communicate general network-level dynamics. Here we demonstrate that few-variable firing rate models with established synaptic timescales can account for both strongly synchronized and weakly modulated gamma oscillations. These models go beyond the classical formulations of rate models by including at least two dynamic variables per population: firing rate and synaptic activation. The models’ flexibility to capture the broad range of gamma behavior depends directly on the timescales that represent recruitment of the excitatory and inhibitory firing rates. In particular, we find that weakly modulated gamma oscillations occur robustly when the recruitment timescale of inhibition is faster than that of excitation. We present our findings by using an extended Wilson-Cowan model and a rate model derived from a network of quadratic integrate-and-fire neurons. These biophysical rate models capture the range of weakly modulated and coherent gamma oscillations observed in spiking network models, while additionally allowing for greater tractability and systems analysis. [ABSTRACT FROM AUTHOR]

In this paper we prove the propagation of chaos property for an ensemble of interacting neurons subject to independent Brownian noise. The propagation of chaos property means that in the large network size limit, the neurons behave as if they are probabilistically independent. The model for the internal dynamics of the neurons is taken to be that of Wilson and Cowan, and we consider there to be multiple different populations. The synaptic connections are modeled with a nonlinear "electrical" model. The nonlinearity of the synaptic connections means that our model lies outside the scope of classical propagation of chaos results. We obtain the propagation of chaos result by taking advantage of the fact that the mean-field equations are Gaussian, which allows us to use Borell's Inequality to prove that its tails decay exponentially. [ABSTRACT FROM AUTHOR]

Gamma rhythm refers to oscillatory neural activity between 30 and 80 Hz, induced in visual cortex by stimuli such as iso-luminant hues or gratings. The power and peak frequency of gamma depend on the properties of the stimulus such as size and contrast. Gamma waveform is typically arch-shaped, with narrow troughs and broad peaks, and can be replicated in a self-oscillating Wilson-Cowan (WC) model operating in an appropriate regime. However, oscillations in this model are infinitely long, unlike physiological gamma that occurs in short bursts. Further, unlike the model, gamma is faster after stimulus onset and slows down over time. Here, we first characterized gamma burst duration in local field potential data recorded from two monkeys as they viewed full screen iso-luminant hues. We then added different types of noise in the inputs to the WC model and tested how that affected duration and temporal dynamics of gamma. While the model failed with the often-used Poisson noise, Ornstein-Uhlenbeck noise applied to both the excitatory and the inhibitory populations replicated the duration and slowing of gamma and replicated the shape and stimulus dependencies. Thus, the temporal dynamics of gamma oscillations put constraints on the type and properties of underlying neural noise., (© The Author(s) 2023. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Abstract Modeling and interpreting (partial) synchronous neural activity can be a challenge. We illustrate this by deriving the phase dynamics of two seminal neural mass models: the Wilson-Cowan firing rate model and the voltage-based Freeman model. We established that the phase dynamics of these models differed qualitatively due to an attractive coupling in the first and a repulsive coupling in the latter. Using empirical structural connectivity matrices, we determined that the two dynamics cover the functional connectivity observed in resting state activity. We further searched for two pivotal dynamical features that have been reported in many experimental studies: (1) a partial phase synchrony with a possibility of a transition towards either a desynchronized or a (fully) synchronized state; (2) long-term autocorrelations indicative of a scale-free temporal dynamics of phase synchronization. Only the Freeman phase model exhibited scale-free behavior. Its repulsive coupling, however, let the individual phases disperse and did not allow for a transition into a synchronized state. The Wilson-Cowan phase model, by contrast, could switch into a (partially) synchronized state, but it did not generate long-term correlations although being located close to the onset of synchronization, i.e. in its critical regime. That is, the phase-reduced models can display one of the two dynamical features, but not both. Highlights • Networks of neural masses are reduced to phase oscillator models. • Freeman- and Wilson-Cowan-based phase models can resemble empirically observed functional connectivity. • Neither a network of Freeman nor of Wilson-Cowan phase models can. • Capture both partial synchronization and scale-free temporal correlations. [ABSTRACT FROM AUTHOR]

We present a detailed computational model of interacting neuronal populations that mimic the hatchling Xenopus tadpole nervous system. The model includes four sensory pathways, integrators of sensory information, and a central pattern generator (CPG) network. Sensory pathways of different modalities receive inputs from an “environment”; these inputs are then processed and integrated to select the most appropriate locomotor action. The CPG populations execute the selected action, generating output in motor neuron populations. Thus, the model describes a detailed and biologically plausible chain of information processing from external signals to sensors, sensory pathways, integration and decision-making, action selection and execution and finally, generation of appropriate motor activity and behaviour. We show how the model produces appropriate behaviours in response to a selected scenario, which consists of a sequence of “environmental” signals. These behaviours might be relatively complex due to noisy sensory pathways and the possibility of spontaneous actions. [ABSTRACT FROM AUTHOR]

The inhibitory restraint necessary to suppress aberrant activity can fail when inhibitory neurons cease to generate action potentials as they enter depolarization block. We investigate possible bifurcation structures that arise at the onset of seizure-like activity resulting from depolarization block in inhibitory neurons. Networks of conductance-based excitatory and inhibitory neurons are simulated to characterize different types of transitions to the seizure state, and a mean field model is developed to verify the generality of the observed phenomena of excitatory-inhibitory dynamics. Specifically, the inhibitory population's activation function in the Wilson-Cowan model is modified to be non-monotonic to reflect that inhibitory neurons enter depolarization block given strong input. We find that a physiological state and a seizure state can coexist, where the seizure state is characterized by high excitatory and low inhibitory firing rate. Bifurcation analysis of the mean field model reveals that a transition to the seizure state may occur via a saddle-node bifurcation or a homoclinic bifurcation. We explain the hysteresis observed in network simulations using these two bifurcation types. We also demonstrate that extracellular potassium concentration affects the depolarization block threshold; the consequent changes in bifurcation structure enable the network to produce the tonic to clonic phase transition observed in biological epileptic networks. [ABSTRACT FROM AUTHOR]

Much progress has been made in uncovering the computational capabilities of spiking neural networks. However, spiking neurons will always be more expensive to simulate compared to rate neurons because of the inherent disparity in time scales - the spike duration time is much shorter than the inter-spike time, which is much shorter than any learning time scale. In numerical analysis, this is a classic stiff problem. Spiking neurons are also much more difficult to study analytically. One possible approach to making spiking networks more tractable is to augment mean field activity models with some information about spiking correlations. For example, such a generalized activity model could carry information about spiking rates and correlations between spikes self-consistently. Here, we will show how this can be accomplished by constructing a complete formal probabilistic description of the network and then expanding around a small parameter such as the inverse of the number of neurons in the network. The mean field theory of the system gives a rate-like description. The first order terms in the perturbation expansion keep track of covariances.

The Wilson-Cowan neural field equations describe the dynamical behavior of a 1-D continuum of excitatory and inhibitory cortical neural aggregates, using a pair of coupled integro-differential equations. Here we use bifurcation theory and small-noise linear stochastics to study the range of a phase transitions--sudden qualitative changes in the state of a dynamical system emerging from a bifurcation-- accessible to the Wilson-Cowan network. Specifically, we examine saddle-node, Hopf, Turing, and Turing-Hopf instabilities. We introduce stochasticity by adding small-amplitude spatio-temporal white noise, and analyze the resulting subthreshold fluctuations using an Ornstein-Uhlenbeck linearization. This analysis predicts divergent changes in correlation and spectral characteristics of neural activity during close approach to bifurcation from below. We validate these theoretical predictions using numerical simulations. The results demonstrate the role of noise in the emergence of critically slowed precursors in both space and time, and suggest that these earlywarning signals are a universal feature of a neural system close to bifurcation. In particular, these precursor signals are likely to have neurobiological significance as early warnings of impending state change in the cortex. We support this claim with an analysis of the in vitro local field potentials recorded from slices of mouse-brain tissue. We show that in the period leading up to emergence of spontaneous seizurelike events, the mouse field potentials show a characteristic spectral focusing toward lower frequencies concomitant with a growth in fluctuation variance, consistent with critical slowing near a bifurcation point. This observation of biological criticality has clear implications regarding the feasibility of seizure prediction. [ABSTRACT FROM AUTHOR]

Although it is known that two coupled Wilson-Cowan models with reciprocal connections induce aperiodic oscillations, little attention has been paid to the dynamical mechanism for such oscillations so far. In this study, we aim to elucidate the fundamental mechanism to induce the aperiodic oscillations in the coupled model. First, aperiodic oscillations observed are investigated for the case when the connections are unidirectional and when the input signal is a periodic oscillation. By the phase portrait analysis, we determine that the aperiodic oscillations are caused by periodically forced state transitions between a stable equilibrium and a stable limit cycle attractors around the saddle-node and saddle separatrix loop bifurcation points. It is revealed that the dynamical mechanism where the state crosses over the saddle-node and saddle separatrix loop bifurcations significantly contributes to the occurrence of chaotic oscillations forced by a periodic input. In addition, this mechanism can also give rise to chaotic oscillations in reciprocally connected Wilson-Cowan models. These results suggest that the dynamic attractor transition underlies chaotic behaviors in two coupled Wilson-Cowan oscillators. [ABSTRACT FROM AUTHOR]

Much progress has been made in uncovering the computational capabilities of spiking neural networks. However, spiking neurons will always be more expensive to simulate compared to rate neurons because of the inherent disparity in time scales-the spike duration time is much shorter than the inter-spike time, which is much shorter than any learning time scale. In numerical analysis, this is a classic stiff problem. Spiking neurons are also much more difficult to study analytically. One possible approach to making spiking networks more tractable is to augment mean field activity models with some information about spiking correlations. For example, such a generalized activity model could carry information about spiking rates and correlations between spikes self-consistently. Here, we will show how this can be accomplished by constructing a complete formal probabilistic description of the network and then expanding around a small parameter such as the inverse of the number of neurons in the network. The mean field theory of the system gives a rate-like description. The first order terms in the perturbation expansion keep track of covariances [ABSTRACT FROM AUTHOR]