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.

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.

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.

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.

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.

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.

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.

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.

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.

Effect of Stimulation Waveform on the Non-linear Entrainment of Cortical Alpha Oscillations.

In the past decade, there has been a surge of interest in using patterned brain stimulation to manipulate cortical oscillations, in both experimental and clinical settings. But the relationship between stimulation waveform and its impact on ongoing oscillations remains poorly understood and severely restrains the development of new paradigms. To address some aspects of this intricate problem, we combine computational and mathematical approaches, providing new insights into the influence of waveform of both low and high-frequency stimuli on synchronous neural activity. Using a cellular-based cortical microcircuit network model, we performed numerical simulations to test the influence of different waveforms on ongoing alpha oscillations, and derived a mean-field description of stimulation-driven dynamics to better understand the observed responses. Our analysis shows that high-frequency periodic stimulation translates into an effective transformation of the neurons' response function, leading to waveform-dependent changes in oscillatory dynamics and resting state activity. Moreover, we found that randomly fluctuating stimulation linearizes the neuron response function while constant input moves its activation threshold. Taken together, our findings establish a new theoretical framework in which stimulation waveforms impact neural systems at the population-scale through non-linear interactions.

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.

Optimal Model Parameter Estimation from EEG Power Spectrum Features Observed during General Anesthesia.

Mathematical modeling is a powerful tool that enables researchers to describe the experimentally observed dynamics of complex systems. Starting with a robust model including model parameters, it is necessary to choose an appropriate set of model parameters to reproduce experimental data. However, estimating an optimal solution of the inverse problem, i.e., finding a set of model parameters that yields the best possible fit to the experimental data, is a very challenging problem. In the present work, we use different optimization algorithms based on a frequentist approach, as well as Monte Carlo Markov Chain methods based on Bayesian inference techniques to solve the considered inverse problems. We first probe two case studies with synthetic data and study models described by a stochastic non-delayed linear second-order differential equation and a stochastic linear delay differential equation. In a third case study, a thalamo-cortical neural mass model is fitted to the EEG spectral power measured during general anesthesia induced by anesthetics propofol and desflurane. We show that the proposed neural mass model fits very well to the observed EEG power spectra, particularly to the power spectral peaks within δ - (0 - 4 Hz) and α - (8 - 13 Hz) frequency ranges. Furthermore, for each case study, we perform a practical identifiability analysis by estimating the confidence regions of the parameter estimates and interpret the corresponding correlation and sensitivity matrices. Our results indicate that estimating the model parameters from analytically computed spectral power, we are able to accurately estimate the unknown parameters while avoiding the computational costs due to numerical integration of the model equations.

Kernel Reconstruction for Delayed Neural Field Equations.

Understanding the neural field activity for realistic living systems is a challenging task in contemporary neuroscience. Neural fields have been studied and developed theoretically and numerically with considerable success over the past four decades. However, to make effective use of such models, we need to identify their constituents in practical systems. This includes the determination of model parameters and in particular the reconstruction of the underlying effective connectivity in biological tissues.In this work, we provide an integral equation approach to the reconstruction of the neural connectivity in the case where the neural activity is governed by a delay neural field equation. As preparation, we study the solution of the direct problem based on the Banach fixed-point theorem. Then we reformulate the inverse problem into a family of integral equations of the first kind. This equation will be vector valued when several neural activity trajectories are taken as input for the inverse problem. We employ spectral regularization techniques for its stable solution. A sensitivity analysis of the regularized kernel reconstruction with respect to the input signal u is carried out, investigating the Fréchet differentiability of the kernel with respect to the signal. Finally, we use numerical examples to show the feasibility of the approach for kernel reconstruction, including numerical sensitivity tests, which show that the integral equation approach is a very stable and promising approach for practical computational neuroscience.

Stochastic resonance mediates the state-dependent effect of periodic stimulation on cortical alpha oscillations.

Brain stimulation can be used to engage and modulate rhythmic activity in brain networks. However, the outcomes of brain stimulation are shaped by behavioral states and endogenous fluctuations in brain activity. To better understand how this intrinsic oscillatory activity controls the susceptibility of the brain to stimulation, we analyzed a computational model of the thalamo-cortical system in two distinct states (rest and task-engaged) to identify the mechanisms by which endogenous alpha oscillations (8Hz-12Hz) are modulated by periodic stimulation. Our analysis shows that the different responses to stimulation observed experimentally in these brain states can be explained by a passage through a bifurcation combined with stochastic resonance - a mechanism by which irregular fluctuations amplify the response of a nonlinear system to weak periodic signals. Indeed, our findings suggest that modulation of brain oscillations is best achieved in states of low endogenous rhythmic activity, and that irregular state-dependent fluctuations in thalamic inputs shape the susceptibility of cortical population to periodic stimulation.

Anesthetic action on the transmission delay between cortex and thalamus explains the beta-buzz observed under propofol anesthesia.

In recent years, more and more surgeries under general anesthesia have been performed with the assistance of electroencephalogram (EEG) monitors. An increase in anesthetic concentration leads to characteristic changes in the power spectra of the EEG. Although tracking the anesthetic-induced changes in EEG rhythms can be employed to estimate the depth of anesthesia, their precise underlying mechanisms are still unknown. A prominent feature in the EEG of some patients is the emergence of a strong power peak in the ß-frequency band, which moves to the α-frequency band while increasing the anesthetic concentration. This feature is called the beta-buzz. In the present study, we use a thalamo-cortical neural population feedback model to reproduce observed characteristic features in frontal EEG power obtained experimentally during propofol general anesthesia, such as this beta-buzz. First, we find that the spectral power peak in the α- and δ-frequency ranges depend on the decay rate constant of excitatory and inhibitory synapses, but the anesthetic action on synapses does not explain the beta-buzz. Moreover, considering the action of propofol on the transmission delay between cortex and thalamus, the model reveals that the beta-buzz may result from a prolongation of the transmission delay by increasing propofol concentration. A corresponding relationship between transmission delay and anesthetic blood concentration is derived. Finally, an analytical stability study demonstrates that increasing propofol concentration moves the systems resting state towards its stability threshold.

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.

Dynamic Control of Synchronous Activity in Networks of Spiking Neurons.

Oscillatory brain activity is believed to play a central role in neural coding. Accumulating evidence shows that features of these oscillations are highly dynamic: power, frequency and phase fluctuate alongside changes in behavior and task demands. The role and mechanism supporting this variability is however poorly understood. We here analyze a network of recurrently connected spiking neurons with time delay displaying stable synchronous dynamics. Using mean-field and stability analyses, we investigate the influence of dynamic inputs on the frequency of firing rate oscillations. We show that afferent noise, mimicking inputs to the neurons, causes smoothing of the system's response function, displacing equilibria and altering the stability of oscillatory states. Our analysis further shows that these noise-induced changes cause a shift of the peak frequency of synchronous oscillations that scales with input intensity, leading the network towards critical states. We lastly discuss the extension of these principles to periodic stimulation, in which externally applied driving signals can trigger analogous phenomena. Our results reveal one possible mechanism involved in shaping oscillatory activity in the brain and associated control principles.

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.

Neural field simulator: two-dimensional spatio-temporal dynamics involving finite transmission speed.

Neural Field models (NFM) play an important role in the understanding of neural population dynamics on a mesoscopic spatial and temporal scale. Their numerical simulation is an essential element in the analysis of their spatio-temporal dynamics. The simulation tool described in this work considers scalar spatially homogeneous neural fields taking into account a finite axonal transmission speed and synaptic temporal derivatives of first and second order. A text-based interface offers complete control of field parameters and several approaches are used to accelerate simulations. A graphical output utilizes video hardware acceleration to display running output with reduced computational hindrance compared to simulators that are exclusively software-based. Diverse applications of the tool demonstrate breather oscillations, static and dynamic Turing patterns and activity spreading with finite propagation speed. The simulator is open source to allow tailoring of code and this is presented with an extension use case.