Intrinsic neural diversity quenches the dynamic volatility of neural networks.

Heterogeneity is the norm in biology. The brain is no different: Neuronal cell types are myriad, reflected through their cellular morphology, type, excitability, connectivity motifs, and ion channel distributions. While this biophysical diversity enriches neural systems' dynamical repertoire, it remains challenging to reconcile with the robustness and persistence of brain function over time (resilience). To better understand the relationship between excitability heterogeneity (variability in excitability within a population of neurons) and resilience, we analyzed both analytically and numerically a nonlinear sparse neural network with balanced excitatory and inhibitory connections evolving over long time scales. Homogeneous networks demonstrated increases in excitability, and strong firing rate correlations-signs of instability-in response to a slowly varying modulatory fluctuation. Excitability heterogeneity tuned network stability in a context-dependent way by restraining responses to modulatory challenges and limiting firing rate correlations, while enriching dynamics during states of low modulatory drive. Excitability heterogeneity was found to implement a homeostatic control mechanism enhancing network resilience to changes in population size, connection probability, strength and variability of synaptic weights, by quenching the volatility (i.e., its susceptibility to critical transitions) of its dynamics. Together, these results highlight the fundamental role played by cell-to-cell heterogeneity in the robustness of brain function in the face of change.

Selective control of synaptic plasticity in heterogeneous networks through transcranial alternating current stimulation (tACS).

Transcranial alternating current stimulation (tACS) represents a promising non-invasive treatment for an increasingly wide range of neurological and neuropsychiatric disorders. The ability to use periodically oscillating electric fields to non-invasively engage neural dynamics opens up the possibility of recruiting synaptic plasticity and to modulate brain function. However, despite consistent reports about tACS clinical effectiveness, strong state-dependence combined with the ubiquitous heterogeneity of cortical networks collectively results in high outcome variability. Introducing variations in intrinsic neuronal timescales, we explored how such heterogeneity influences stimulation-induced change in synaptic connectivity. We examined how spike timing dependent plasticity, at the level of cells, intra- and inter-laminar cortical networks, can be selectively and preferentially engaged by periodic stimulation. Using leaky integrate-and-fire neuron models, we analyzed cortical circuits comprised of multiple cell-types, alongside superficial multi-layered networks expressing distinct layer-specific timescales. Our results show that mismatch in neuronal timescales within and/or between cells-and the resulting variability in excitability, temporal integration properties and frequency tuning-enables selective and directional control on synaptic connectivity by tACS. Our work provides new vistas on how to recruit neural heterogeneity to guide brain plasticity using non-invasive stimulation paradigms.

Diversity-induced trivialization and resilience of neural dynamics.

Heterogeneity is omnipresent across all living systems. Diversity enriches the dynamical repertoire of these systems but remains challenging to reconcile with their manifest robustness and dynamical persistence over time, a fundamental feature called resilience. To better understand the mechanism underlying resilience in neural circuits, we considered a nonlinear network model, extracting the relationship between excitability heterogeneity and resilience. To measure resilience, we quantified the number of stationary states of this network, and how they are affected by various control parameters. We analyzed both analytically and numerically gradient and non-gradient systems modeled as non-linear sparse neural networks evolving over long time scales. Our analysis shows that neuronal heterogeneity quenches the number of stationary states while decreasing the susceptibility to bifurcations: a phenomenon known as trivialization. Heterogeneity was found to implement a homeostatic control mechanism enhancing network resilience to changes in network size and connection probability by quenching the system's dynamic volatility.

Arousal Fluctuations Govern Oscillatory Transitions Between Dominant [Formula: see text] and [Formula: see text] Occipital Activity During Eyes Open/Closed Conditions.

Arousal results in widespread activation of brain areas to increase their response in task and behavior relevant ways. Mediated by the Ascending Reticular Arousal System (ARAS), arousal-dependent inputs interact with neural circuitry to shape their dynamics. In the occipital cortex, such inputs may trigger shifts between dominant oscillations, where [Formula: see text] activity is replaced by [Formula: see text] activity, or vice versa. A salient example of this are spectral power alternations observed while eyes are opened and/or closed. These transitions closely follow fluctuations in arousal, suggesting a common origin. To better understand the mechanisms at play, we developed and analyzed a computational model composed of two modules: a thalamocortical feedback circuit coupled with a superficial cortical network. Upon activation by noise-like inputs originating from the ARAS, our model is able to demonstrate that noise-driven non-linear interactions mediate transitions in dominant peak frequency, resulting in the simultaneous suppression of [Formula: see text] limit cycle activity and the emergence of [Formula: see text] oscillations through coherence resonance. Reduction in input provoked the reverse effect - leading to anticorrelated transitions between [Formula: see text] and [Formula: see text] power. Taken together, these results shed a new light on how arousal shapes oscillatory brain activity.

Suppression of underlying neuronal fluctuations mediates EEG slowing during general anaesthesia.

The physiological mechanisms by which anaesthetic drugs modulate oscillatory brain activity remain poorly understood. Combining human data, mathematical and computational analysis of both spiking and mean-field models, we investigated the spectral dynamics of encephalographic (EEG) beta-alpha oscillations, observed in human patients undergoing general anaesthesia. The effect of anaesthetics can be modelled as a reduction of neural fluctuation intensity, and/or an increase in inhibitory synaptic gain in the thalamo-cortical circuit. Unlike previous work, which suggested the primary importance of gamma-amino-butryic-acid (GABA) augmentation in causing a shift to low EEG frequencies, our analysis demonstrates that a non-linear transition, triggered by a simple decrease in neural fluctuation intensity, is sufficient to explain the clinically-observed appearance - and subsequent slowing - of the beta-alpha narrowband EEG peak. In our model, increased synaptic inhibition alone, did not correlate with the clinically-observed encephalographic spectral changes, but did cause the anaesthetic-induced decrease in neuronal firing rate. Taken together, our results show that such a non-linear transition results in functional fragmentation of cortical and thalamic populations; highly correlated intra-population dynamics triggered by anaesthesia decouple and isolate neural populations. Our results are able to parsimoniously unify and replicate the observed anaesthetic effects on both the EEG spectra and inter-regional connectivity, and further highlight the importance of neural activity fluctuations in the genesis of altered brain states.

Breakdown of local information processing may underlie isoflurane anesthesia effects.

The disruption of coupling between brain areas has been suggested as the mechanism underlying loss of consciousness in anesthesia. This hypothesis has been tested previously by measuring the information transfer between brain areas, and by taking reduced information transfer as a proxy for decoupling. Yet, information transfer is a function of the amount of information available in the information source-such that transfer decreases even for unchanged coupling when less source information is available. Therefore, we reconsidered past interpretations of reduced information transfer as a sign of decoupling, and asked whether impaired local information processing leads to a loss of information transfer. An important prediction of this alternative hypothesis is that changes in locally available information (signal entropy) should be at least as pronounced as changes in information transfer. We tested this prediction by recording local field potentials in two ferrets after administration of isoflurane in concentrations of 0.0%, 0.5%, and 1.0%. We found strong decreases in the source entropy under isoflurane in area V1 and the prefrontal cortex (PFC)-as predicted by our alternative hypothesis. The decrease in source entropy was stronger in PFC compared to V1. Information transfer between V1 and PFC was reduced bidirectionally, but with a stronger decrease from PFC to V1. This links the stronger decrease in information transfer to the stronger decrease in source entropy-suggesting reduced source entropy reduces information transfer. This conclusion fits the observation that the synaptic targets of isoflurane are located in local cortical circuits rather than on the synapses formed by interareal axonal projections. Thus, changes in information transfer under isoflurane seem to be a consequence of changes in local processing more than of decoupling between brain areas. We suggest that source entropy changes must be considered whenever interpreting changes in information transfer as decoupling.

Shaping Intrinsic Neural Oscillations with Periodic Stimulation.

UNLABELLED: Rhythmic brain activity plays an important role in neural processing and behavior. Features of these oscillations, including amplitude, phase, and spectrum, can be influenced by internal states (e.g., shifts in arousal, attention or cognitive ability) or external stimulation. Electromagnetic stimulation techniques such as transcranial magnetic stimulation, transcranial direct current stimulation, and transcranial alternating current stimulation are used increasingly in both research and clinical settings. Currently, the mechanisms whereby time-dependent external stimuli influence population-scale oscillations remain poorly understood. Here, we provide computational insights regarding the mapping between periodic pulsatile stimulation parameters such as amplitude and frequency and the response dynamics of recurrent, nonlinear spiking neural networks. Using a cortical model built of excitatory and inhibitory neurons, we explored a wide range of stimulation intensities and frequencies systematically. Our results suggest that rhythmic stimulation can form the basis of a control paradigm in which one can manipulate the intrinsic oscillatory properties of driven networks via a plurality of input-driven mechanisms. Our results show that, in addition to resonance and entrainment, nonlinear acceleration is involved in shaping the rhythmic response of our modeled network. Such nonlinear acceleration of spontaneous and synchronous oscillatory activity in a neural network occurs in regimes of intense, high-frequency rhythmic stimulation. These results open new perspectives on the manipulation of synchronous neural activity for basic and clinical research. SIGNIFICANCE STATEMENT: Oscillatory activity is widely recognized as a core mechanism for information transmission within and between brain circuits. Noninvasive stimulation methods can shape this activity, something that is increasingly capitalized upon in basic research and clinical practice. Here, we provide computational insights on the mechanistic bases for such effects. Our results show that rhythmic stimulation forms the basis of a control paradigm in which one can manipulate the intrinsic oscillatory properties of driven networks via a plurality of input-driven mechanisms. In addition to resonance and entrainment, nonlinear acceleration is involved in shaping the rhythmic response of our modeled network, particularly in regimes of high-frequency rhythmic stimulation. These results open new perspectives on the manipulation of synchronous neural activity for basic and clinical research.

Stimulus statistics shape oscillations in nonlinear recurrent neural networks.

Rhythmic activity plays a central role in neural computations and brain functions ranging from homeostasis to attention, as well as in neurological and neuropsychiatric disorders. Despite this pervasiveness, little is known about the mechanisms whereby the frequency and power of oscillatory activity are modulated, and how they reflect the inputs received by neurons. Numerous studies have reported input-dependent fluctuations in peak frequency and power (as well as couplings across these features). However, it remains unresolved what mediates these spectral shifts among neural populations. Extending previous findings regarding stochastic nonlinear systems and experimental observations, we provide analytical insights regarding oscillatory responses of neural populations to stimulation from either endogenous or exogenous origins. Using a deceptively simple yet sparse and randomly connected network of neurons, we show how spiking inputs can reliably modulate the peak frequency and power expressed by synchronous neural populations without any changes in circuitry. Our results reveal that a generic, non-nonlinear and input-induced mechanism can robustly mediate these spectral fluctuations, and thus provide a framework in which inputs to the neurons bidirectionally regulate both the frequency and power expressed by synchronous populations. Theoretical and computational analysis of the ensuing spectral fluctuations was found to reflect the underlying dynamics of the input stimuli driving the neurons. Our results provide insights regarding a generic mechanism supporting spectral transitions observed across cortical networks and spanning multiple frequency bands.

Awake vs. anesthetized: layer-specific sensory processing in visual cortex and functional connectivity between cortical areas.

During general anesthesia, global brain activity and behavioral state are profoundly altered. Yet it remains mostly unknown how anesthetics alter sensory processing across cortical layers and modulate functional cortico-cortical connectivity. To address this gap in knowledge of the micro- and mesoscale effects of anesthetics on sensory processing in the cortical microcircuit, we recorded multiunit activity and local field potential in awake and anesthetized ferrets (Mustela putoris furo) during sensory stimulation. To understand how anesthetics alter sensory processing in a primary sensory area and the representation of sensory input in higher-order association areas, we studied the local sensory responses and long-range functional connectivity of primary visual cortex (V1) and prefrontal cortex (PFC). Isoflurane combined with xylazine provided general anesthesia for all anesthetized recordings. We found that anesthetics altered the duration of sensory-evoked responses, disrupted the response dynamics across cortical layers, suppressed both multimodal interactions in V1 and sensory responses in PFC, and reduced functional cortico-cortical connectivity between V1 and PFC. Together, the present findings demonstrate altered sensory responses and impaired functional network connectivity during anesthesia at the level of multiunit activity and local field potential across cortical layers.

How the cortico-thalamic feedback affects the EEG power spectrum over frontal and occipital regions during propofol-induced sedation.

Increasing concentrations of the anaesthetic agent propofol initially induces sedation before achieving full general anaesthesia. During this state of anaesthesia, the observed specific changes in electroencephalographic (EEG) rhythms comprise increased activity in the δ- (0.5-4 Hz) and α- (8-13 Hz) frequency bands over the frontal region, but increased δ- and decreased α-activity over the occipital region. It is known that the cortex, the thalamus, and the thalamo-cortical feedback loop contribute to some degree to the propofol-induced changes in the EEG power spectrum. However the precise role of each structure to the dynamics of the EEG is unknown. In this paper we apply a thalamo-cortical neuronal population model to reproduce the power spectrum changes in EEG during propofol-induced anaesthesia sedation. The model reproduces the power spectrum features observed experimentally both in frontal and occipital electrodes. Moreover, a detailed analysis of the model indicates the importance of multiple resting states in brain activity. The work suggests that the α-activity originates from the cortico-thalamic relay interaction, whereas the emergence of δ-activity results from the full cortico-reticular-relay-cortical feedback loop with a prominent enforced thalamic reticular-relay interaction. This model suggests an important role for synaptic GABAergic receptors at relay neurons and, more generally, for the thalamus in the generation of both the δ- and the α- EEG patterns that are seen during propofol anaesthesia sedation.

Detecting event-related recurrences by symbolic analysis: applications to human language processing.

Quasi-stationarity is ubiquitous in complex dynamical systems. In brain dynamics, there is ample evidence that event-related potentials (ERPs) reflect such quasi-stationary states. In order to detect them from time series, several segmentation techniques have been proposed. In this study, we elaborate a recent approach for detecting quasi-stationary states as recurrence domains by means of recurrence analysis and subsequent symbolization methods. We address two pertinent problems of contemporary recurrence analysis: optimizing the size of recurrence neighbourhoods and identifying symbols from different realizations for sequence alignment. As possible solutions for these problems, we suggest a maximum entropy criterion and a Hausdorff clustering algorithm. The resulting recurrence domains for single-subject ERPs are obtained as partition cells reflecting quasi-stationary brain states.

Study of GABAergic extra-synaptic tonic inhibition in single neurons and neural populations by traversing neural scales: application to propofol-induced anaesthesia.

Anaesthetic agents are known to affect extra-synaptic GABAergic receptors, which induce tonic inhibitory currents. Since these receptors are very sensitive to small concentrations of agents, they are supposed to play an important role in the underlying neural mechanism of general anaesthesia. Moreover anaesthetic agents modulate the encephalographic activity (EEG) of subjects and hence show an effect on neural populations. To understand better the tonic inhibition effect in single neurons on neural populations and hence how it affects the EEG, the work considers single neurons and neural populations in a steady-state and studies numerically and analytically the modulation of their firing rate and nonlinear gain with respect to different levels of tonic inhibition. We consider populations of both type-I (Leaky Integrate-and-Fire model) and type-II (Morris-Lecar model) neurons. To bridge the single neuron description to the population description analytically, a recently proposed statistical approach is employed which allows to derive new analytical expressions for the population firing rate for type-I neurons. In addition, the work shows the derivation of a novel transfer function for type-I neurons as considered in neural mass models and studies briefly the interaction of synaptic and extra-synaptic inhibition. We reveal a strong subtractive and divisive effect of tonic inhibition in type-I neurons, i.e. a shift of the firing rate to higher excitation levels accompanied by a change of the nonlinear gain. Tonic inhibition shortens the excitation window of type-II neurons and their populations while maintaining the nonlinear gain. The gained results are interpreted in the context of recent experimental findings under propofol-induced anaesthesia.

Anesthesia differentially modulates spontaneous network dynamics by cortical area and layer.

Anesthesia is widely used in medicine and research to achieve altered states of consciousness and cognition. Whereas changes to macroscopic cortical activity patterns by anesthesia measured at the spatial resolution of electroencephalography have been widely studied, modulation of mesoscopic and microscopic network dynamics by anesthesia remain poorly understood. To address this gap in knowledge, we recorded spontaneous mesoscopic (local field potential) and microscopic (multiunit activity) network dynamics in primary visual cortex (V1) and prefrontal cortex (PFC) of awake and isoflurane anesthetized ferrets (Mustela putoris furo). This approach allowed for examination of activity as a function of cortical area, cortical layer, and anesthetic depth with much higher spatial and temporal resolution than in previous studies. We hypothesized that a primary sensory area and an association cortical area would exhibit different patterns of network modulation by anesthesia due to their different functional roles. Indeed, we found effects specific to cortical area and cortical layer. V1 exhibited minimal changes in rhythmic structure with anesthesia but differential modulation of input layer IV. In contrast, anesthesia profoundly altered spectral power in PFC, with more uniform modulation across cortical layers. Our results demonstrate that anesthesia modulates spontaneous cortical activity in an area- and layer-specific manner. These finding provide the basis for 1) refining anesthesia monitoring algorithms, 2) reevaluating the large number of systems neuroscience studies performed in anesthetized animals, and 3) increasing our understanding of differential dynamics across cortical layers and areas.

Detecting recurrence domains of dynamical systems by symbolic dynamics.

We propose an algorithm for the detection of recurrence domains of complex dynamical systems from time series. Our approach exploits the characteristic checkerboard texture of recurrence domains exhibited in recurrence plots. In phase space, recurrence plots yield intersecting balls around sampling points that could be merged into cells of a phase space partition. We construct this partition by a rewriting grammar applied to the symbolic dynamics of time indices. A maximum entropy principle defines the optimal size of intersecting balls. The final application to high-dimensional brain signals yields an optimal symbolic recurrence plot revealing functional components of the signal.

Arousal system stimulation and anesthetic state alter visuoparietal connectivity.

Cortical information processing is under the precise control of the ascending arousal system (AAS). Anesthesia suppresses cortical arousal that can be mitigated by exogenous stimulation of the AAS. The question remains to what extent cortical information processing is regained by AAS stimulation. We investigate the effect of electrical stimulation of the nucleus Pontis Oralis (PnO), a distinct source of ascending AAS projections, on cortical functional connectivity (FC) and information storage at mild, moderate, and deep anesthesia. Local field potentials (LFPs) recorded previously in the secondary visual cortex (V2) and the adjacent parietal association cortex (PtA) in chronically instrumented unrestrained rats. We hypothesized that PnO stimulation would induce electrocortical arousal accompanied by enhanced FC and active information storage (AIS) implying improved information processing. In fact, stimulation reduced FC in slow oscillations (0.3-2.5 Hz) at low anesthetic level and increased FC at high anesthetic level. These effects were augmented following stimulation suggesting stimulus-induced plasticity. The observed opposite stimulation-anesthetic impact was less clear in the γ-band activity (30-70 Hz). In addition, FC in slow oscillations was more sensitive to stimulation and anesthetic level than FC in γ-band activity which exhibited a rather constant spatial FC structure that was symmetric between specific, topographically related sites in V2 and PtA. Invariant networks were defined as a set of strongly connected electrode channels, which were invariant to experimental conditions. In invariant networks, stimulation decreased AIS and increasing anesthetic level increased AIS. Conversely, in non-invariant (complement) networks, stimulation did not affect AIS at low anesthetic level but increased it at high anesthetic level. The results suggest that arousal stimulation alters cortical FC and information storage as a function of anesthetic level with a prolonged effect beyond the duration of stimulation. The findings help better understand how the arousal system may influence information processing in cortical networks at different levels of anesthesia.

Delayed closed-loop neurostimulation for the treatment of pathological brain rhythms in mental disorders: a computational study.

Mental disorders are among the top most demanding challenges in world-wide health. A large number of mental disorders exhibit pathological rhythms, which serve as the disorders characteristic biomarkers. These rhythms are the targets for neurostimulation techniques. Open-loop neurostimulation employs stimulation protocols, which are rather independent of the patients health and brain state in the moment of treatment. Most alternative closed-loop stimulation protocols consider real-time brain activity observations but appear as adaptive open-loop protocols, where e.g., pre-defined stimulation sets in if observations fulfil pre-defined criteria. The present theoretical work proposes a fully-adaptive closed-loop neurostimulation setup, that tunes the brain activities power spectral density (PSD) according to a user-defined PSD. The utilized brain model is non-parametric and estimated from the observations via magnitude fitting in a pre-stimulus setup phase. Moreover, the algorithm takes into account possible conduction delays in the feedback connection between observation and stimulation electrode. All involved features are illustrated on pathological α- and γ-rhythms known from psychosis. To this end, we simulate numerically a linear neural population brain model and a non-linear cortico-thalamic feedback loop model recently derived to explain brain activity in psychosis.

Reduced dynamics for delayed systems with harmonic or stochastic forcing.

The analysis of nonlinear delay-differential equations (DDEs) subjected to external forcing is difficult due to the infinite dimensionality of the space in which they evolve. To simplify the analysis of such systems, the present work develops a non-homogeneous center manifold (CM) reduction scheme, which allows the derivation of a time-dependent order parameter equation in finite dimension. This differential equation captures the major dynamical features of the delayed system. The forcing is assumed to be small compared to the amplitude of the autonomous system, in order to cause only small variations of the fixed points and of the autonomous CM. The time-dependent CM is shown to satisfy a non-homogeneous partial differential equation. We first briefly review CM theory for DDEs. Then we show, for the general scalar case, how an ansatz that separates the CM into one for the autonomous problem plus an additional time-dependent order-two correction leads to satisfying results. The paper then details the application to a transcritical bifurcation subjected to single or multiple periodic forcings. The validity limits of the reduction scheme are also highlighted. Finally, we characterize the specific case of additive stochastic driving of the transcritical bifurcation, where additive white noise shifts the mode of the probability density function of the state variable to larger amplitudes.

Mathematical Model Insights into EEG Origin under Transcranial Direct Current Stimulation (tDCS) in the Context of Psychosis.

Schizophrenia is a psychotic disease that develops progressively over years with a transition from prodromal to psychotic state associated with a disruption in brain activity. Transcranial Direct Current Stimulation (tDCS), known to alleviate pharmaco-resistant symptoms in patients suffering from schizophrenia, promises to prevent such a psychotic transition. To understand better how tDCS affects brain activity, we propose a neural cortico-thalamo-cortical (CTC) circuit model involving the Ascending Reticular Arousal System (ARAS) that permits to describe major impact features of tDCS, such as excitability for short-duration stimulation and electroencephalography (EEG) power modulation for long-duration stimulation. To this end, the mathematical model relates stimulus duration and Long-Term Plasticity (LTP) effect, in addition to describing the temporal LTP decay after stimulus offset. This new relation promises to optimize future stimulation protocols. Moreover, we reproduce successfully EEG-power modulation under tDCS in a ketamine-induced psychosis model and confirm the N-methyl-d-aspartate (NMDA) receptor hypofunction hypothesis in the etiopathophysiology of schizophrenia. The model description points to an important role of the ARAS and the δ-rhythm synchronicity in CTC circuit in early-stage psychosis.

Digital Addiction and Sleep.

In 2020, the World Health Organization formally recognized addiction to digital technology (connected devices) as a worldwide problem, where excessive online activity and internet use lead to inability to manage time, energy, and attention during daytime and produce disturbed sleep patterns or insomnia during nighttime. Recent studies have shown that the problem has increased in magnitude worldwide during the COVID-19 pandemic. The extent to which dysfunctional sleep is a consequence of altered motivation, memory function, mood, diet, and other lifestyle variables or results from excess of blue-light exposure when looking at digital device screens for long hours at day and night is one of many still unresolved questions. This article offers a narrative overview of some of the most recent literature on this topic. The analysis provided offers a conceptual basis for understanding digital addiction as one of the major reasons why people, and adolescents in particular, sleep less and less well in the digital age. It discusses definitions as well as mechanistic model accounts in context. Digital addiction is identified as functionally equivalent to all addictions, characterized by the compulsive, habitual, and uncontrolled use of digital devices and an excessively repeated engagement in a particular online behavior. Once the urge to be online has become uncontrollable, it is always accompanied by severe sleep loss, emotional distress, depression, and memory dysfunction. In extreme cases, it may lead to suicide. The syndrome has been linked to the known chronic effects of all drugs, producing disturbances in cellular and molecular mechanisms of the GABAergic and glutamatergic neurotransmitter systems. Dopamine and serotonin synaptic plasticity, essential for impulse control, memory, and sleep function, are measurably altered. The full spectrum of behavioral symptoms in digital addicts include eating disorders and withdrawal from outdoor and social life. Evidence pointing towards dysfunctional melatonin and vitamin D metabolism in digital addicts should be taken into account for carving out perspectives for treatment. The conclusions offer a holistic account for digital addiction, where sleep deficit is one of the key factors.

Neurostimulation stabilizes spiking neural networks by disrupting seizure-like oscillatory transitions.

An improved understanding of the mechanisms underlying neuromodulatory approaches to mitigate seizure onset is needed to identify clinical targets for the treatment of epilepsy. Using a Wilson-Cowan-motivated network of inhibitory and excitatory populations, we examined the role played by intrinsic and extrinsic stimuli on the network's predisposition to sudden transitions into oscillatory dynamics, similar to the transition to the seizure state. Our joint computational and mathematical analyses revealed that such stimuli, be they noisy or periodic in nature, exert a stabilizing influence on network responses, disrupting the development of such oscillations. Based on a combination of numerical simulations and mean-field analyses, our results suggest that high variance and/or high frequency stimulation waveforms can prevent multi-stability, a mathematical harbinger of sudden changes in network dynamics. By tuning the neurons' responses to input, stimuli stabilize network dynamics away from these transitions. Furthermore, our research shows that such stabilization of neural activity occurs through a selective recruitment of inhibitory cells, providing a theoretical undergird for the known key role these cells play in both the healthy and diseased brain. Taken together, these findings provide new vistas on neuromodulatory approaches to stabilize neural microcircuit activity.