Self-sustained asynchronous irregular states and Up-Down states in thalamic, cortical and thalamocortical networks of nonlinear integrate-and-fire neurons

Emergent effects of synaptic connectivity on the dynamics of global and local slow waves in a large-scale thalamocortical network model of the human brain

Slow-wave sleep (SWS), characterized by slow oscillations (SO, <1Hz) of alternating active and silent states in the thalamocortical network, is a primary brain state during Non-Rapid Eye Movement (NREM) sleep. In the last two decades, the traditional view of SWS as a global and uniform whole-brain state has been challenged by a growing body of evidence indicating that SO can be local and can coexist with wake-like activity. However, the understanding of how global and local SO emerges from micro-scale neuron dynamics and network connectivity remains unclear. We developed a multi-scale, biophysically realistic human whole-brain thalamocortical network model capable of transitioning between the awake state and slow-wave sleep, and we investigated the role of connectivity in the spatio-temporal dynamics of sleep SO. We found that the overall strength and a relative balance between long and short-range synaptic connections determined the network state. Importantly, for a range of synaptic strengths, the model demonstrated complex mixed SO states, where periods of synchronized global slow-wave activity were intermittent with the periods of asynchronous local slow-waves. Increase of the overall synaptic strength led to synchronized global SO, while decrease of synaptic connectivity produced only local slow-waves that would not propagate beyond local area. These results were compared to human data to validate probable models of biophysically realistic SO. The model producing mixed states provided the best match to the spatial coherence profile and the functional connectivity estimated from human subjects. These findings shed light on how the spatio-temporal properties of SO emerge from local and global cortical connectivity and provide a framework for further exploring the mechanisms and functions of SWS in health and disease.

A Translaminar Spacetime Code Supports Touch-Evoked Traveling Waves

Linking sensory-evoked traveling waves to underlying circuit patterns is critical to understanding the neural basis of sensory perception. To form this link, we performed simultaneous electrophysiology and two-photon calcium imaging through transparent NeuroGrids and mapped touch-evoked cortical traveling waves and their underlying microcircuit dynamics. In awake mice, both passive and active whisker touch elicited traveling waves within and across barrels, with a fast early component followed by a variable late wave that lasted hundreds of milliseconds post-stimulus. Strikingly, late-wave dynamics were modulated by stimulus value and correlated with task performance. Mechanistically, the late wave component was i) modulated by motor feedback, ii) complemented by a sparse ensemble pattern across layer 2/3, which a balanced-state network model reconciled via inhibitory stabilization, and iii) aligned to regenerative Layer-5 apical dendritic Ca 2+ events. Our results reveal a translaminar spacetime pattern organized by cortical feedback in the sensory cortex that supports touch-evoked traveling waves.

GRAPHICAL ABSTRACT AND HIGHLIGHTS

Whisker touch evokes both early- and late-traveling waves in the barrel cortex over 100's of millisecondsReward reinforcement modulates wave dynamics Late wave emergence coincides with network sparsity in L23 and time-locked L5 dendritic Ca 2+ spikes Experimental and computational results link motor feedback to distinct translaminar spacetime patterns.

A mean-field model of gamma-frequency oscillations in networks of excitatory and inhibitory neurons

Gamma oscillations are widely seen in the cerebral cortex in different states of the wake-sleep cycle and are thought to play a role in sensory processing and cognition. Here, we study the emergence of gamma oscillations at two levels, in networks of spiking neurons, and a mean-field model. At the network level, we consider two different mechanisms to generate gamma oscillations and show that they are best seen if one takes into account the synaptic delay between neurons. At the mean-field level, we show that, by introducing delays, the mean-field can also produce gamma oscillations. The mean-field matches the mean activity of excitatory and inhibitory populations of the spiking network, as well as their oscillation frequencies, for both mechanisms. This mean-field model of gamma oscillations should be a useful tool to investigate large-scale interactions through gamma oscillations in the brain.

Sleep waves in a large-scale corticothalamic model constrained by activities intrinsic to neocortical networks and single thalamic neurons

AIM

Many biophysical and non-biophysical models have been able to reproduce the corticothalamic activities underlying different EEG sleep rhythms but none of them included the known ability of neocortical networks and single thalamic neurons to generate some of these waves intrinsically.

METHODS

We built a large-scale corticothalamic model with a high fidelity in anatomical connectivity consisting of a single cortical column and first- and higher-order thalamic nuclei. The model is constrained by different neocortical excitatory and inhibitory neuronal populations eliciting slow (<1 Hz) oscillations and by thalamic neurons generating sleep waves when isolated from the neocortex.

RESULTS

Our model faithfully reproduces all EEG sleep waves and the transition from a desynchronized EEG to spindles, slow (<1 Hz) oscillations, and delta waves by progressively increasing neuronal membrane hyperpolarization as it occurs in the intact brain. Moreover, our model shows that slow (<1 Hz) waves most often start in a small assembly of thalamocortical neurons though they can also originate in cortical layer 5. Moreover, the input of thalamocortical neurons increases the frequency of EEG slow (<1 Hz) waves compared to those generated by isolated cortical networks.

CONCLUSION

Our simulations challenge current mechanistic understanding of the temporal dynamics of sleep wave generation and suggest testable predictions.

Sleep and wake in a model of the thalamocortical system with Martinotti cells

The mechanisms leading to the alternation between active (UP) and silent (DOWN) states during sleep slow waves (SWs) remain poorly understood. Previous models have explained the transition to the DOWN state by a progressive failure of excitation because of the build-up of adaptation currents or synaptic depression. However, these models are at odds with recent studies suggesting a role for presynaptic inhibition by Martinotti cells (MaCs) in generating SWs. Here, we update a classical large-scale model of sleep SWs to include MaCs and propose a different mechanism for the generation of SWs. In the wake mode, the network exhibits irregular and selective activity with low firing rates (FRs). Following an increase in the strength of background inputs and a modulation of synaptic strength and potassium leak potential mimicking the reduced effect of acetylcholine during sleep, the network enters a sleep-like regime in which local increases of network activity trigger bursts of MaC activity, resulting in strong disfacilitation of the local network via presynaptic GABAB1a -type inhibition. This model replicates findings on slow wave activity (SWA) during sleep that challenge previous models, including low and skewed FRs that are comparable between the wake and sleep modes, higher synchrony of transitions to DOWN states than to UP states, the possibility of triggering SWs by optogenetic stimulation of MaCs, and the local dependence of SWA on synaptic strength. Overall, this work points to a role for presynaptic inhibition by MaCs in the generation of DOWN states during sleep.

High-Density Exploration of Activity States in a Multi-Area Brain Model

To simulate whole brain dynamics with only a few equations, biophysical, mesoscopic models of local neuron populations can be connected using empirical tractography data. The development of mesoscopic mean-field models of neural populations, in particular, the Adaptive Exponential (AdEx mean-field model), has successfully summarized neuron-scale phenomena leading to the emergence of global brain dynamics associated with conscious (asynchronous and rapid dynamics) and unconscious (synchronized slow-waves, with Up-and-Down state dynamics) brain states, based on biophysical mechanisms operating at cellular scales (e.g. neuromodulatory regulation of spike-frequency adaptation during sleep-wake cycles or anesthetics). Using the Virtual Brain (TVB) environment to connect mean-field AdEx models, we have previously simulated the general properties of brain states, playing on spike-frequency adaptation, but have not yet performed detailed analyses of other parameters possibly also regulating transitions in brain-scale dynamics between different brain states. We performed a dense grid parameter exploration of the TVB-AdEx model, making use of High Performance Computing. We report a remarkable robustness of the effect of adaptation to induce synchronized slow-wave activity. Moreover, the occurrence of slow waves is often paralleled with a closer relation between functional and structural connectivity. We find that hyperpolarization can also generate unconscious-like synchronized Up and Down states, which may be a mechanism underlying the action of anesthetics. We conclude that the TVB-AdEx model reveals large-scale properties identified experimentally in sleep and anesthesia.

Mechanisms underlying pathological cortical bursts during metabolic depletion

Cortical activity depends upon a continuous supply of oxygen and other metabolic resources. Perinatal disruption of oxygen availability is a common clinical scenario in neonatal intensive care units, and a leading cause of lifelong disability. Pathological patterns of brain activity including burst suppression and seizures are a hallmark of the recovery period, yet the mechanisms by which these patterns arise remain poorly understood. Here, we use computational modeling of coupled metabolic-neuronal activity to explore the mechanisms by which oxygen depletion generates pathological brain activity. We find that restricting oxygen supply drives transitions from normal activity to several pathological activity patterns (isoelectric, burst suppression, and seizures), depending on the potassium supply. Trajectories through parameter space track key features of clinical electrophysiology recordings and reveal how infants with good recovery outcomes track toward normal parameter values, whereas the parameter values for infants with poor outcomes dwell around the pathological values. These findings open avenues for studying and monitoring the metabolically challenged infant brain, and deepen our understanding of the link between neuronal and metabolic activity.

M-current modulation of cortical slow oscillations: Network dynamics and computational modeling

The slow oscillation is a synchronized network activity expressed by the cortical network in slow wave sleep and under anesthesia. Waking up requires a transition from this synchronized brain state to a desynchronized one. Cholinergic innervation is critical for the transition from slow-wave-sleep to wakefulness, and muscarinic action is largely exerted through the muscarinic-sensitive potassium current (M-current) block. We investigated the dynamical impact of blocking the M-current on slow oscillations, both in cortical slices and in a cortical network computational model. Blocking M-current resulted in an elongation of Up states (by four times) and in a significant firing rate increase, reflecting an increased network excitability, albeit no epileptiform discharges occurred. These effects were replicated in a biophysical cortical model, where a parametric reduction of the M-current resulted in a progressive elongation of Up states and firing rate. All neurons, and not only those modeled with M-current, increased their firing rates due to network recurrency. Further increases in excitability induced even longer Up states, approaching the microarousals described in the transition towards wakefulness. Our results bridge an ionic current with network modulation, providing a mechanistic insight into network dynamics of awakening.

Assessing brain state and anesthesia level with two-photon calcium signals

Brain states, such as wake, sleep, or different depths of anesthesia are usually assessed using electrophysiological techniques, such as the local field potential (LFP) or the electroencephalogram (EEG), which are ideal signals for detecting activity patterns such as asynchronous or oscillatory activities. However, it is technically challenging to have these types of measures during calcium imaging recordings such as two-photon or wide-field techniques. Here, using simultaneous two-photon and LFP measurements, we demonstrate that despite the slower dynamics of the calcium signal, there is a high correlation between the LFP and two-photon signals taken from the neuropil outside neuronal somata. Moreover, we find the calcium signal to be systematically delayed from the LFP signal, and we use a model to show that the delay between the two signals is due to the physical distance between the recording sites. These results suggest that calcium signals alone can be used to detect activity patterns such as slow oscillations and ultimately assess the brain state and level of anesthesia.

A comprehensive neural simulation of slow-wave sleep and highly responsive wakefulness dynamics

Hallmarks of neural dynamics during healthy human brain states span spatial scales from neuromodulators acting on microscopic ion channels to macroscopic changes in communication between brain regions. Developing a scale-integrated understanding of neural dynamics has therefore remained challenging. Here, we perform the integration across scales using mean-field modeling of Adaptive Exponential (AdEx) neurons, explicitly incorporating intrinsic properties of excitatory and inhibitory neurons. The model was run using The Virtual Brain (TVB) simulator, and is open-access in EBRAINS. We report that when AdEx mean-field neural populations are connected via structural tracts defined by the human connectome, macroscopic dynamics resembling human brain activity emerge. Importantly, the model can qualitatively and quantitatively account for properties of empirically observed spontaneous and stimulus-evoked dynamics in space, time, phase, and frequency domains. Large-scale properties of cortical dynamics are shown to emerge from both microscopic-scale adaptation that control transitions between wake-like to sleep-like activity, and the organization of the human structural connectome; together, they shape the spatial extent of synchrony and phase coherence across brain regions consistent with the propagation of sleep-like spontaneous traveling waves at intermediate scales. Remarkably, the model also reproduces brain-wide, enhanced responsiveness and capacity to encode information particularly during wake-like states, as quantified using the perturbational complexity index. The model was run using The Virtual Brain (TVB) simulator, and is open-access in EBRAINS. This approach not only provides a scale-integrated understanding of brain states and their underlying mechanisms, but also open access tools to investigate brain responsiveness, toward producing a more unified, formal understanding of experimental data from conscious and unconscious states, as well as their associated pathologies.

Creation of Neuronal Ensembles and Cell-Specific Homeostatic Plasticity through Chronic Sparse Optogenetic Stimulation

Cortical computations emerge from the dynamics of neurons embedded in complex cortical circuits. Within these circuits, neuronal ensembles, which represent subnetworks with shared functional connectivity, emerge in an experience-dependent manner. Here we induced ensembles in ex vivo cortical circuits from mice of either sex by differentially activating subpopulations through chronic optogenetic stimulation. We observed a decrease in voltage correlation, and importantly a synaptic decoupling between the stimulated and nonstimulated populations. We also observed a decrease in firing rate during Up-states in the stimulated population. These ensemble-specific changes were accompanied by decreases in intrinsic excitability in the stimulated population, and a decrease in connectivity between stimulated and nonstimulated pyramidal neurons. By incorporating the empirically observed changes in intrinsic excitability and connectivity into a spiking neural network model, we were able to demonstrate that changes in both intrinsic excitability and connectivity accounted for the decreased firing rate, but only changes in connectivity accounted for the observed decorrelation. Our findings help ascertain the mechanisms underlying the ability of chronic patterned stimulation to create ensembles within cortical circuits and, importantly, show that while Up-states are a global network-wide phenomenon, functionally distinct ensembles can preserve their identity during Up-states through differential firing rates and correlations.SIGNIFICANCE STATEMENT The connectivity and activity patterns of local cortical circuits are shaped by experience. This experience-dependent reorganization of cortical circuits is driven by complex interactions between different local learning rules, external input, and reciprocal feedback between many distinct brain areas. Here we used an ex vivo approach to demonstrate how simple forms of chronic external stimulation can shape local cortical circuits in terms of their correlated activity and functional connectivity. The absence of feedback between different brain areas and full control of external input allowed for a tractable system to study the underlying mechanisms and development of a computational model. Results show that differential stimulation of subpopulations of neurons significantly reshapes cortical circuits and forms subnetworks referred to as neuronal ensembles.

Computational models of Idling brain activity for memory processing

Studying the underlying neural mechanisms of cognitive functions of the brain is one of the central questions in modern biology. Moreover, it has significantly impacted the development of novel technologies in artificial intelligence. Spontaneous activity is a unique feature of the brain and is currently lacking in many artificially constructed intelligent machines. Spontaneous activity may represent the brain's idling states, which are internally driven by neuronal networks and possibly participate in offline processing during awake, sleep, and resting states. Evidence is accumulating that the brain's spontaneous activity is not mere noise but part of the mechanisms to process information about previous experiences. A bunch of literature has shown how previous sensory and behavioral experiences influence the subsequent patterns of brain activity with various methods in various animals. It seems, however, that the patterns of neural activity and their computational roles differ significantly from area to area and from function to function. In this article, I review the various forms of the brain's spontaneous activity, especially those observed during memory processing, and some attempts to model the generation mechanisms and computational roles of such activities.

Connectivity concepts in neuronal network modeling

Sustainable research on computational models of neuronal networks requires published models to be understandable, reproducible, and extendable. Missing details or ambiguities about mathematical concepts and assumptions, algorithmic implementations, or parameterizations hinder progress. Such flaws are unfortunately frequent and one reason is a lack of readily applicable standards and tools for model description. Our work aims to advance complete and concise descriptions of network connectivity but also to guide the implementation of connection routines in simulation software and neuromorphic hardware systems. We first review models made available by the computational neuroscience community in the repositories ModelDB and Open Source Brain, and investigate the corresponding connectivity structures and their descriptions in both manuscript and code. The review comprises the connectivity of networks with diverse levels of neuroanatomical detail and exposes how connectivity is abstracted in existing description languages and simulator interfaces. We find that a substantial proportion of the published descriptions of connectivity is ambiguous. Based on this review, we derive a set of connectivity concepts for deterministically and probabilistically connected networks and also address networks embedded in metric space. Beside these mathematical and textual guidelines, we propose a unified graphical notation for network diagrams to facilitate an intuitive understanding of network properties. Examples of representative network models demonstrate the practical use of the ideas. We hope that the proposed standardizations will contribute to unambiguous descriptions and reproducible implementations of neuronal network connectivity in computational neuroscience.

Not so spontaneous: Multi-dimensional representations of behaviors and context in sensory areas

Sensory areas are spontaneously active in the absence of sensory stimuli. This spontaneous activity has long been studied; however, its functional role remains largely unknown. Recent advances in technology, allowing large-scale neural recordings in the awake and behaving animal, have transformed our understanding of spontaneous activity. Studies using these recordings have discovered high-dimensional spontaneous activity patterns, correlation between spontaneous activity and behavior, and dissimilarity between spontaneous and sensory-driven activity patterns. These findings are supported by evidence from developing animals, where a transition toward these characteristics is observed as the circuit matures, as well as by evidence from mature animals across species. These newly revealed characteristics call for the formulation of a new role for spontaneous activity in neural sensory computation.

State-dependent activity dynamics of hypothalamic stress effector neurons

The stress response necessitates an immediate boost in vital physiological functions from their homeostatic operation to an elevated emergency response. However, the neural mechanisms underlying this state-dependent change remain largely unknown. Using a combination of in vivo and ex vivo electrophysiology with computational modeling, we report that corticotropin releasing hormone (CRH) neurons in the paraventricular nucleus of the hypothalamus (PVN), the effector neurons of hormonal stress response, rapidly transition between distinct activity states through recurrent inhibition. Specifically, in vivo optrode recording shows that under non-stress conditions, CRH_PVN neurons often fire with rhythmic brief bursts (RB), which, somewhat counterintuitively, constrains firing rate due to long (~2 s) interburst intervals. Stressful stimuli rapidly switch RB to continuous single spiking (SS), permitting a large increase in firing rate. A spiking network model shows that recurrent inhibition can control this activity-state switch, and more broadly the gain of spiking responses to excitatory inputs. In biological CRH_PVN neurons ex vivo, the injection of whole-cell currents derived from our computational model recreates the in vivo-like switch between RB and SS, providing direct evidence that physiologically relevant network inputs enable state-dependent computation in single neurons. Together, we present a novel mechanism for state-dependent activity dynamics in CRH_PVN neurons.

Waveform detection by deep learning reveals multi-area spindles that are selectively modulated by memory load

Sleep is generally considered to be a state of large-scale synchrony across thalamus and neocortex; however, recent work has challenged this idea by reporting isolated sleep rhythms such as slow oscillations and spindles. What is the spatial scale of sleep rhythms? To answer this question, we adapted deep learning algorithms initially developed for detecting earthquakes and gravitational waves in high-noise settings for analysis of neural recordings in sleep. We then studied sleep spindles in non-human primate electrocorticography (ECoG), human electroencephalogram (EEG), and clinical intracranial electroencephalogram (iEEG) recordings in the human. Within each recording type, we find widespread spindles occur much more frequently than previously reported. We then analyzed the spatiotemporal patterns of these large-scale, multi-area spindles and, in the EEG recordings, how spindle patterns change following a visual memory task. Our results reveal a potential role for widespread, multi-area spindles in consolidation of memories in networks widely distributed across primate cortex.

The brain processes memories as we sleep, generating rhythms of electrical activity called ‘sleep spindles’. Sleep spindles were long thought to be a state where the entire brain was fully synchronized by this rhythm. This was based on EEG recordings, short for electroencephalogram, a technique that uses electrodes on the scalp to measure electrical activity in the outermost layer of the brain, the cortex. But more recent intracranial recordings of people undergoing brain surgery have challenged this idea and suggested that sleep spindles may not be a state of global brain synchronization, but rather localised to specific areas.

Mofrad et al. sought to clarify the extent to which spindles co-occur at multiple sites in the brain, which could shed light on how networks of neurons coordinate memory storage during sleep. To analyse highly variable brain wave recordings, Mofrad et al. adapted deep learning algorithms initially developed for detecting earthquakes and gravitational waves. The resulting algorithm, designed to more sensitively detect spindles amongst other brain activity, was then applied to a range of sleep recordings from humans and macaque monkeys.

The analyses revealed that widespread and complex patterns of spindle rhythms, spanning multiple areas in the cortex of the brain, actually appear much more frequently than previously thought. This finding was consistent across all the recordings analysed, even recordings under the skull, which provide the clearest window into brain circuits. Further analyses found that these multi-area spindles occurred more often in sleep after people had completed tasks that required holding many visual scenes in memory, as opposed to control conditions with fewer visual scenes.

In summary, Mofrad et al. show that neuroscientists had previously not appreciated the complex and dynamic patterns in this sleep rhythm. These patterns in sleep spindles may be able to adapt based on the demands needed for memory storage, and this will be the subject of future work. Moreover, the findings support the idea that sleep spindles help coordinate the consolidation of memories in brain circuits that stretch across the cortex. Understanding this mechanism may provide insights into how memory falters in aging and sleep-related diseases, such as Alzheimer’s disease. Lastly, the algorithm developed by Mofrad et al. stands to be a useful tool for analysing other rhythmic waveforms in noisy recordings.

Information transmission in recurrent networks: Consequences of network noise for synchronous and asynchronous signal encoding

Information about natural time-dependent stimuli encoded by the sensory periphery or communication between cortical networks may span a large frequency range or be localized to a smaller frequency band. Biological systems have been shown to multiplex such disparate broadband and narrow-band signals and then discriminate them in later populations by employing either an integration (low-pass) or coincidence detection (bandpass) encoding strategy. Analytical expressions have been developed for both encoding methods in feedforward populations of uncoupled neurons and confirm that the integration of a population's output low-pass filters the information, whereas synchronous output encodes less information overall and retains signal information in a selected frequency band. The present study extends the theory to recurrent networks and shows that recurrence may sharpen the synchronous bandpass filter. The frequency of the pass band is significantly influenced by the synaptic strengths, especially for inhibition-dominated networks. Synchronous information transfer is also increased when network models take into account heterogeneity that arises from the stochastic distribution of the synaptic weights.

Cortical oscillations support sampling-based computations in spiking neural networks

Being permanently confronted with an uncertain world, brains have faced evolutionary pressure to represent this uncertainty in order to respond appropriately. Often, this requires visiting multiple interpretations of the available information or multiple solutions to an encountered problem. This gives rise to the so-called mixing problem: since all of these "valid" states represent powerful attractors, but between themselves can be very dissimilar, switching between such states can be difficult. We propose that cortical oscillations can be effectively used to overcome this challenge. By acting as an effective temperature, background spiking activity modulates exploration. Rhythmic changes induced by cortical oscillations can then be interpreted as a form of simulated tempering. We provide a rigorous mathematical discussion of this link and study some of its phenomenological implications in computer simulations. This identifies a new computational role of cortical oscillations and connects them to various phenomena in the brain, such as sampling-based probabilistic inference, memory replay, multisensory cue combination, and place cell flickering.

Slow waves form expanding, memory-rich mesostates steered by local excitability in fading anesthesia

In the arousal process, the brain restores its integrative activity from the synchronized state of slow wave activity (SWA). The mechanisms underpinning this state transition remain, however, to be elucidated. Here we simultaneously probed neuronal assemblies throughout the whole cortex with micro-electrocorticographic recordings in mice. We investigated the progressive shaping of propagating SWA at different levels of isoflurane. We found a form of memory of the wavefront shapes at deep anesthesia, tightly alternating posterior-anterior-posterior patterns. At low isoflurane, metastable patterns propagated in more directions, reflecting an increased complexity. The wandering across these mesostates progressively increased its randomness, as predicted by simulations of a network of spiking neurons, and confirmed in our experimental data. The complexity increase is explained by the elevated excitability of local assemblies with no modifications of the network connectivity. These results shed new light on the functional reorganization of the cortical network as anesthesia fades out.

•

Complexity of isoflurane-induced slow waves reliably determines anesthesia level

•

In deep anesthesia, the propagation strictly alternates between front-back-front patterns

•

In light anesthesia, there is a continuum of directions and faster propagation

•

Local excitability underpins the cortical reorganization in fading anesthesia

Critical current for giant Fano factor in neural models with bistable firing dynamics and implications for signal transmission

Bistability in the firing rate is a prominent feature in different types of neurons as well as in neural networks. We show that for a constant input below a critical value, such bistability can lead to a giant spike-count diffusion. We study the transmission of a periodic signal and demonstrate that close to the critical bias current, the signal-to-noise ratio suffers a sharp increase, an effect that can be traced back to the giant diffusion and large Fano factor.

Cholinergic modulation of hierarchical inhibitory control over cortical resting state dynamics: Local circuit modeling of schizophrenia-related hypofrontality

Nicotinic acetylcholine receptors (nAChRs) modulate the cholinergic drive to a hierarchy of inhibitory neurons in the superficial layers of the PFC, critical to cognitive processes. It has been shown that genetic deletions of the various types of nAChRs impact the properties of ultra-slow transitions between high and low PFC activity states in mice during quiet wakefulness. The impact characteristics depend on specific interneuron populations expressing the manipulated receptor subtype. In addition, recent data indicate that a genetic mutation of the α5 nAChR subunit, located on vasoactive intestinal polypeptide (VIP) inhibitory neurons, the rs16969968 single nucleotide polymorphism (α5 SNP), plays a key role in the hypofrontality observed in schizophrenia patients carrying the SNP. Data also indicate that chronic nicotine application to α5 SNP mice relieves the hypofrontality. We developed a computational model to show that the activity patterns recorded in the genetically modified mice can be explained by changes in the dynamics of the local PFC circuit. Notably, our model shows that these altered PFC circuit dynamics are due to changes in the stability structure of the activity states. We identify how this stability structure is differentially modulated by cholinergic inputs to the parvalbumin (PV), somatostatin (SOM) or the VIP inhibitory populations. Our model uncovers that a change in amplitude, but not duration of the high activity states can account for the lowered pyramidal (PYR) population firing rates recorded in α5 SNP mice. We demonstrate how nicotine-induced desensitization and upregulation of the β2 nAChRs located on SOM interneurons, as opposed to the activation of α5 nAChRs located on VIP interneurons, is sufficient to explain the nicotine-induced activity normalization in α5 SNP mice. The model further implies that subsequent nicotine withdrawal may exacerbate the hypofrontality over and beyond one caused by the SNP mutation.

•

Prefrontal cortex shows ultra-slow alterations between low and high activity states at rest. This activity is characteristically decreased in schizophrenia patients.

•

Model identifies local circuit origin of hypofrontality associated with schizophrenia and a5 nicotinic receptor malfunction.

•

Decrease in PFC VIP-interneuron excitability drives decrease in high-activity-state stability and overall hypofrontality.

•

Model shows desensitization/upregulation of SOM-expressed β2-NAChRs drive nicotine-induced renormalization of PFC activity.

Spatial synaptic modulation through IP3 diffusion triggered by ECB: a computational study with an astrocyte-neurons model

Recently, functional interactions between neurons and astrocytes have been steadily clarified. In particular, the effects of presynaptic depolarization-induced suppression of excitation (DSE) through endocannabinoid (ECB) and endocannabinoids-mediated synaptic potentiation (eSP) by an astrocyte have been used as an evidence of global heterosynaptic modulation. However, the mechanism of occurrence of spatial modulation in a neural network remains unknown. Although the Ca²⁺ density in astrocytes is strongly related to eSP through ECB, the mechanism of the rise in the ECB receptor in Ca²⁺ remains unclear. Since Ca²⁺ is closely related to inositol-1,4,5-trisphosphate (IP₃), it is believed that the released IP₃ affects Ca²⁺ in astrocytes that receive ECB. Therefore, this study approximately showed the spatial distribution of DSE or eSP with astrocyte-neuron computational models, assuming that the IP₃ caused by ECB is transmitted in an astrocyte. The results showed doughnut-shaped DSE, eSP, and DSE regions from the central ECB released points to the surroundings. They suggested that IP₃ diffusion plays an important role in mediating this synaptic modulation.

Spontaneous traveling waves naturally emerge from horizontal fiber time delays and travel through locally asynchronous-irregular states

Studies of sensory-evoked neuronal responses often focus on mean spike rates, with fluctuations treated as internally-generated noise. However, fluctuations of spontaneous activity, often organized as traveling waves, shape stimulus-evoked responses and perceptual sensitivity. The mechanisms underlying these waves are unknown. Further, it is unclear whether waves are consistent with the low rate and weakly correlated “asynchronous-irregular” dynamics observed in cortical recordings. Here, we describe a large-scale computational model with topographically-organized connectivity and conduction delays relevant to biological scales. We find that spontaneous traveling waves are a general property of these networks. The traveling waves that occur in the model are sparse, with only a small fraction of neurons participating in any individual wave. Consequently, they do not induce measurable spike correlations and remain consistent with locally asynchronous irregular states. Further, by modulating local network state, they can shape responses to incoming inputs as observed in vivo.

Spontaneous traveling cortical waves shape neural responses. Using a large-scale computational model, the authors show that transmission delays shape locally asynchronous spiking dynamics into traveling waves without inducing correlations and boost responses to external input, as observed in vivo.

Nonlinear collision between propagating waves in mouse somatosensory cortex

How does cellular organization shape the spatio-temporal patterns of activity in the cortex while processing sensory information? After measuring the propagation of activity in the mouse primary somatosensory cortex (S1) in response to single whisker deflections with Voltage Sensitive Dye (VSD) imaging, we developed a two dimensional model of S1. We designed an inference method to reconstruct model parameters from VSD data, revealing that a spatially heterogeneous organization of synaptic strengths between pyramidal neurons in S1 is likely to be responsible for the heterogeneous spatio-temporal patterns of activity measured experimentally. The model shows that, for strong enough excitatory cortical interactions, whisker deflections generate a propagating wave in S1. Finally, we report that two consecutive stimuli activating different spatial locations in S1 generate two waves which collide sub-linearly, giving rise to a suppressive wave. In the inferred model, the suppressive wave is explained by a lower sensitivity to external perturbations of neural networks during activated states.

Frequency Selectivity of Neural Circuits With Heterogeneous Discrete Transmission Delays

Neurons are connected to other neurons by axons and dendrites that conduct signals with finite velocities, resulting in delays between the firing of a neuron and the arrival of the resultant impulse at other neurons. Since delays greatly complicate the analytical treatment and interpretation of models, they are usually neglected or taken to be uniform, leading to a lack in the comprehension of the effects of delays in neural systems. This letter shows that heterogeneous transmission delays make small groups of neurons respond selectively to inputs with differing frequency spectra. By studying a single integrate-and-fire neuron receiving correlated time-shifted inputs, it is shown how the frequency response is linked to both the strengths and delay times of the afferent connections. The results show that incorporating delays alters the functioning of neural networks, and changes the effect that neural connections and synaptic strengths have.

Mutual Interactions between Brain States and Alzheimer's Disease Pathology: A Focus on Gamma and Slow Oscillations

Brain state varies from moment to moment. While brain state can be defined by ongoing neuronal population activity, such as neuronal oscillations, this is tightly coupled with certain behavioural or vigilant states. In recent decades, abnormalities in brain state have been recognised as biomarkers of various brain diseases and disorders. Intriguingly, accumulating evidence also demonstrates mutual interactions between brain states and disease pathologies: while abnormalities in brain state arise during disease progression, manipulations of brain state can modify disease pathology, suggesting a therapeutic potential. In this review, by focusing on Alzheimer's disease (AD), the most common form of dementia, we provide an overview of how brain states change in AD patients and mouse models, and how controlling brain states can modify AD pathology. Specifically, we summarise the relationship between AD and changes in gamma and slow oscillations. As pathological changes in these oscillations correlate with AD pathology, manipulations of either gamma or slow oscillations can modify AD pathology in mouse models. We argue that neuromodulation approaches to target brain states are a promising non-pharmacological intervention for neurodegenerative diseases.

Chaos in memory function of sleep: A nonlinear dynamical analysis in thalamocortical study

Studying the dynamical behaviors of neuronal models may help in better understanding of real nervous system. In addition, it can help researchers to understand some specific phenomena in neuronal system. The thalamocortical network is made of neurons in the thalamus and cortex. In it, the memory function is consolidated in sleep by creating up and down state oscillations (1 Hz) and fast (13-17 Hz) - slow (8-12 Hz) spindles. Recently, a nonlinear biological model for up-down oscillations and fast-slow spindles of the thalamocortical network has been proposed. In this research, the power spectral for the fast-slow spindle of the model is extracted. Dynamical properties of the model, such as the bifurcation diagrams, and attractors are investigated. The results show that the variation of the synaptic power between the excitatory neurons of the cortex and the reticular neurons in the thalamus changes the spindles' activity. According to previous experimental findings, it is an essential rule for consolidating the memory function during sleep. It is also pointed out that when the fast-slow spindles of the brain increase, the dynamics of the thalamocortical system tend to chaos.

A novel density-based neural mass model for simulating neuronal network dynamics with conductance-based synapses and membrane current adaptation

Despite its success in understanding brain rhythms, the neural mass model, as a low-dimensional mean-field network model, is phenomenological in nature, so that it cannot replicate some of rich repertoire of responses seen in real neuronal tissues. Here, using a colored-synapse population density method, we derived a novel neural mass model, termed density-based neural mass model (dNMM), as the mean-field description of network dynamics of adaptive exponential integrate-and-fire (aEIF) neurons, in which two critical neuronal features, i.e., voltage-dependent conductance-based synaptic interactions and adaptation of firing rate responses, were included. Our results showed that the dNMM was capable of correctly estimating firing rate responses of a neuronal population of aEIF neurons receiving stationary or time-varying excitatory and inhibitory inputs. Finally, it was also able to quantitatively describe the effect of spike-frequency adaptation in the generation of asynchronous irregular activity of excitatory-inhibitory cortical networks. We conclude that in terms of its biological reality and calculation efficiency, the dNMM is a suitable candidate to build significantly large-scale network models involving multiple brain areas, where the neuronal population is the smallest dynamic unit.

Thalamocortical Spectral Transmission Relies on Balanced Input Strengths

The thalamus is a key element of sensory transmission in the brain, as it gates and selects sensory streams through a modulation of its internal activity. A preponderant role in these functions is played by its internal activity in the alpha range ([8–14] Hz), but the mechanism underlying this process is not completely understood. In particular, how do thalamocortical connections convey stimulus driven information selectively over the back-ground of thalamic internally generated activity? Here we investigate this issue with a spiking network model of feedforward connectivity between thalamus and primary sensory cortex reproducing the local field potential of both areas. We found that in a feedforward network, thalamic oscillations in the alpha range do not entrain cortical activity for two reasons: (i) alpha range oscillations are weaker in neurons projecting to the cortex, (ii) the gamma resonance dynamics of cortical networks hampers oscillations over the 10–20 Hz range thus weakening alpha range oscillations. This latter mechanism depends on the balance of the strength of thalamocortical connections toward excitatory and inhibitory neurons in the cortex. Our results highlight the relevance of corticothalamic feedback to sustain alpha range oscillations and pave the way toward an integrated understanding of the sensory streams traveling between the periphery and the cortex.

Attractor competition enriches cortical dynamics during awakening from anesthesia

Slow oscillations (≲ 1 Hz), a hallmark of slow-wave sleep and deep anesthesia across species, arise from spatiotemporal patterns of activity whose complexity increases as wakefulness is approached and cognitive functions emerge. The arousal process constitutes an open window to the unknown mechanisms underlying the emergence of such dynamical richness in awake cortical networks. Here, we investigate the changes in network dynamics as anesthesia fades out in the rat visual cortex. Starting from deep anesthesia, slow oscillations gradually increase their frequency, eventually expressing maximum regularity. This stage is followed by the abrupt onset of an infra-slow (~0.2 Hz) alternation between sleep-like oscillations and activated states. A population rate model reproduces this transition driven by an increased excitability that brings it to periodically cross a critical point. Based on our model, dynamical richness emerges as a competition between two metastable attractor states, a conclusion strongly supported by the data.

Growth rules for the repair of Asynchronous Irregular neuronal networks after peripheral lesions

Several homeostatic mechanisms enable the brain to maintain desired levels of neuronal activity. One of these, homeostatic structural plasticity, has been reported to restore activity in networks disrupted by peripheral lesions by altering their neuronal connectivity. While multiple lesion experiments have studied the changes in neurite morphology that underlie modifications of synapses in these networks, the underlying mechanisms that drive these changes are yet to be explained. Evidence suggests that neuronal activity modulates neurite morphology and may stimulate neurites to selective sprout or retract to restore network activity levels. We developed a new spiking network model of peripheral lesioning and accurately reproduced the characteristics of network repair after deafferentation that are reported in experiments to study the activity dependent growth regimes of neurites. To ensure that our simulations closely resemble the behaviour of networks in the brain, we model deafferentation in a biologically realistic balanced network model that exhibits low frequency Asynchronous Irregular (AI) activity as observed in cerebral cortex. Our simulation results indicate that the re-establishment of activity in neurons both within and outside the deprived region, the Lesion Projection Zone (LPZ), requires opposite activity dependent growth rules for excitatory and inhibitory post-synaptic elements. Analysis of these growth regimes indicates that they also contribute to the maintenance of activity levels in individual neurons. Furthermore, in our model, the directional formation of synapses that is observed in experiments requires that pre-synaptic excitatory and inhibitory elements also follow opposite growth rules. Lastly, we observe that our proposed structural plasticity growth rules and the inhibitory synaptic plasticity mechanism that also balances our AI network both contribute to the restoration of the network to pre-deafferentation stable activity levels.

An accumulating body of evidence suggests that our brain can compensate for peripheral lesions by adaptive rewiring of its neuronal circuitry. The underlying process, structural plasticity, can modify the connectivity of neuronal networks in the brain, thus affecting their function. To better understand the mechanisms of structural plasticity in the brain, we have developed a novel model of peripheral lesions and the resulting activity-dependent rewiring in a simplified balanced cortical network model that exhibits biologically realistic Asynchronous Irregular (AI) activity. In order to accurately reproduce the directionality and course of network rewiring after injury that is observed in peripheral lesion experiments, we derive activity dependent growth rules for different synaptic elements: dendritic and axonal contacts. Our simulation results suggest that excitatory and inhibitory synaptic elements have to react to changes in neuronal activity in opposite ways. We show that these rules result in a homeostatic stabilisation of activity in individual neurons. In our simulations, both synaptic and structural plasticity mechanisms contribute to network repair. Furthermore, our simulations indicate that while activity is restored in neurons deprived by the peripheral lesion, the temporal firing characteristics of the network may not be retained by the rewiring process.

Polynomial, piecewise-Linear, Step (PLS): A Simple, Scalable, and Efficient Framework for Modeling Neurons

Biological neurons can be modeled with different levels of biophysical/biochemical details. The accuracy with which a model reflects the actual physiological processes and ultimately the information function of a neuron, can range from very detailed to a schematic phenomenological representation. This range exists due to the common problem: one needs to find an optimal trade-off between the level of details needed to capture the necessary information processing in a neuron and the computational load needed to compute 1 s of model time. An increase in modeled network size or model-time, for which the solution should be obtained, makes this trade-off pivotal in model development. Numerical simulations become incredibly challenging when an extensive network with a detailed representation of each neuron needs to be modeled over a long time interval to study slow evolving processes, e.g., development of the thalamocortical circuits. Here we suggest a simple, powerful and flexible approach in which we approximate the right-hand sides of differential equations by combinations of functions from three families: Polynomial, piecewise-Linear, Step (PLS). To obtain a single coherent framework, we provide four core principles in which PLS functions should be combined. We show the rationale behind each of the core principles. Two examples illustrate how to build a conductance-based or phenomenological model using the PLS-framework. We use the first example as a benchmark on three different computational platforms: CPU, GPU, and mobile system-on-chip devices. We show that the PLS-framework speeds up computations without increasing the memory footprint and maintains high model fidelity comparable to the fully-computed model or with lookup-table approximation. We are convinced that the full range of neuron models: from biophysical to phenomenological and even to abstract models, may benefit from using the PLS-framework.

Emergence of Neuronal Synchronisation in Coupled Areas

One of the most fundamental questions in the field of neuroscience is the emergence of synchronous behaviour in the brain, such as phase, anti-phase, and shift-phase synchronisation. In this work, we investigate how the connectivity between brain areas can influence the phase angle and the neuronal synchronisation. To do this, we consider brain areas connected by means of excitatory and inhibitory synapses, in which the neuron dynamics is given by the adaptive exponential integrate-and-fire model. Our simulations suggest that excitatory and inhibitory connections from one area to another play a crucial role in the emergence of these types of synchronisation. Thus, in the case of unidirectional interaction, we observe that the phase angles of the neurons in the receiver area depend on the excitatory and inhibitory synapses which arrive from the sender area. Moreover, when the neurons in the sender area are synchronised, the phase angle variability of the receiver area can be reduced for some conductance values between the areas. For bidirectional interactions, we find that phase and anti-phase synchronisation can emerge due to excitatory and inhibitory connections. We also verify, for a strong inhibitory-to-excitatory interaction, the existence of silent neuronal activities, namely a large number of excitatory neurons that remain in silence for a long time.

A Potential Mechanism of Sodium Channel Mediating the General Anesthesia Induced by Propofol

General anesthesia has revolutionized healthcare over the past 200 years and continues to show advancements. However, many phenomena induced by general anesthetics including paradoxical excitation are still poorly understood. Voltage-gated sodium channels (Na_V) were believed to be one of the proteins targeted during general anesthesia. Based on electrophysiological measurements before and after propofol treatments of different concentrations, we mathematically modified the Hodgkin–Huxley sodium channel formulations and constructed a thalamocortical model to investigate the potential roles of Na_V. The ion channels of individual neurons were modeled using the Hodgkin–Huxley type equations. The enhancement of propofol-induced GABAa current was simulated by increasing the maximal conductance and the time-constant of decay. Electroencephalogram (EEG) was evaluated as the post-synaptic potential from pyramidal (PY) cells. We found that a left shift in activation of Na_V was induced primarily by a low concentration of propofol (0.3–10 μM), while a left shift in inactivation of Na_V was induced by an increasing concentration (0.3–30 μM). Mathematical simulation indicated that a left shift of Na_V activation produced a Hopf bifurcation, leading to cell oscillations. Left shift of Na_V activation around a value of 5.5 mV in the thalamocortical models suppressed normal bursting of thalamocortical (TC) cells by triggering its chaotic oscillations. This led to irregular spiking of PY cells and an increased frequency in EEG readings. This observation suggests a mechanism leading to paradoxical excitation during general anesthesia. While a left shift in inactivation led to light hyperpolarization in individual cells, it inhibited the activity of the thalamocortical model after a certain depth of anesthesia. This finding implies that high doses of propofol inhibit the network partly by accelerating Na_V toward inactivation. Additionally, this result explains why the application of sodium channel blockers decreases the requirement for general anesthetics. Our study provides an insight into the roles that Na_V plays in the mechanism of general anesthesia. Since the activation and inactivation of Na_V are structurally independent, it should be possible to avoid side effects by state-dependent binding to the Na_V to achieve precision medicine in the future.

Simple models including energy and spike constraints reproduce complex activity patterns and metabolic disruptions

In this work, we introduce new phenomenological neuronal models (eLIF and mAdExp) that account for energy supply and demand in the cell as well as the inactivation of spike generation how these interact with subthreshold and spiking dynamics. Including these constraints, the new models reproduce a broad range of biologically-relevant behaviors that are identified to be crucial in many neurological disorders, but were not captured by commonly used phenomenological models. Because of their low dimensionality eLIF and mAdExp open the possibility of future large-scale simulations for more realistic studies of brain circuits involved in neuronal disorders. The new models enable both more accurate modeling and the possibility to study energy-associated disorders over the whole time-course of disease progression instead of only comparing the initially healthy status with the final diseased state. These models, therefore, provide new theoretical and computational methods to assess the opportunities of early diagnostics and the potential of energy-centered approaches to improve therapies.

Neurons, even “at rest”, require a constant supply of energy to function. They cannot sustain high-frequency activity over long periods because of regulatory mechanisms, such as adaptation or sodium channels inactivation, and metabolic limitations. These limitations are especially severe in many neuronal disorders, where energy can become insufficient and make the neuronal response change drastically, leading to increased burstiness, network oscillations, or seizures. Capturing such behaviors and impact of energy constraints on them is an essential prerequisite to study disorders such as Parkinson’s disease and epilepsy. However, energy and spiking constraints are not present in any of the standard neuronal models used in computational neuroscience. Here we introduce models that provide a simple and scalable way to account for these features, enabling large-scale theoretical and computational studies of neurological disorders and activity patterns that could not be captured by previously used models. These models provide a way to study energy-associated disorders over the whole time-course of disease progression, and they enable a better assessment of energy-centered approaches to improve therapies.

Influence of Autapses on Synchronization in Neural Networks With Chemical Synapses

A great deal of research has been devoted on the investigation of neural dynamics in various network topologies. However, only a few studies have focused on the influence of autapses, synapses from a neuron onto itself via closed loops, on neural synchronization. Here, we build a random network with adaptive exponential integrate-and-fire neurons coupled with chemical synapses, equipped with autapses, to study the effect of the latter on synchronous behavior. We consider time delay in the conductance of the pre-synaptic neuron for excitatory and inhibitory connections. Interestingly, in neural networks consisting of both excitatory and inhibitory neurons, we uncover that synchronous behavior depends on their synapse type. Our results provide evidence on the synchronous and desynchronous activities that emerge in random neural networks with chemical, inhibitory and excitatory synapses where neurons are equipped with autapses.

Deterministic characteristics of spontaneous activity detected by multi-fractal analysis in a spiking neural network with long-tailed distributions of synaptic weights

Cortical neural networks maintain autonomous electrical activity called spontaneous activity that represents the brain's dynamic internal state even in the absence of sensory stimuli. The spatio-temporal complexity of spontaneous activity is strongly related to perceptual, learning, and cognitive brain functions; multi-fractal analysis can be utilized to evaluate the complexity of spontaneous activity. Recent studies have shown that the deterministic dynamic behavior of spontaneous activity especially reflects the topological neural network characteristics and changes of neural network structures. However, it remains unclear whether multi-fractal analysis, recently widely utilized for neural activity, is effective for detecting the complexity of the deterministic dynamic process. To verify this point, we focused on the log-normal distribution of excitatory postsynaptic potentials (EPSPs) to evaluate the multi-fractality of spontaneous activity in a spiking neural network with a log-normal distribution of EPSPs. We found that the spiking activities exhibited multi-fractal characteristics. Moreover, to investigate the presence of a deterministic process in the spiking activity, we conducted a surrogate data analysis against the time-series of spiking activity. The results showed that the spontaneous spiking activity included the deterministic dynamic behavior. Overall, the combination of multi-fractal analysis and surrogate data analysis can detect deterministic complex neural activity. The multi-fractal analysis of neural activity used in this study could be widely utilized for brain modeling and evaluation methods for signals obtained by neuroimaging modalities.

Neuronal spike-rate adaptation supports working memory in language processing

Language processing involves the ability to store and integrate pieces of information in working memory over short periods of time. According to the dominant view, information is maintained through sustained, elevated neural activity. Other work has argued that short-term synaptic facilitation can serve as a substrate of memory. Here we propose an account where memory is supported by intrinsic plasticity that downregulates neuronal firing rates. Single neuron responses are dependent on experience, and we show through simulations that these adaptive changes in excitability provide memory on timescales ranging from milliseconds to seconds. On this account, spiking activity writes information into coupled dynamic variables that control adaptation and move at slower timescales than the membrane potential. From these variables, information is continuously read back into the active membrane state for processing. This neuronal memory mechanism does not rely on persistent activity, excitatory feedback, or synaptic plasticity for storage. Instead, information is maintained in adaptive conductances that reduce firing rates and can be accessed directly without cued retrieval. Memory span is systematically related to both the time constant of adaptation and baseline levels of neuronal excitability. Interference effects within memory arise when adaptation is long lasting. We demonstrate that this mechanism is sensitive to context and serial order which makes it suitable for temporal integration in sequence processing within the language domain. We also show that it enables the binding of linguistic features over time within dynamic memory registers. This work provides a step toward a computational neurobiology of language.

A kernel-based method to calculate local field potentials from networks of spiking neurons

BACKGROUND

The local field potential (LFP) is usually calculated from current sources arising from transmembrane currents, in particular in asymmetric cellular morphologies such as pyramidal neurons.

NEW METHOD

Here, we adopt a different point of view and relate the spiking of neurons to the LFP through efferent synaptic connections and provide a method to calculate LFPs.

RESULTS

We show that the so-called unitary LFPs (uLFP) provide the key to such a calculation. We show experimental measurements and simulations of uLFPs in neocortex and hippocampus, for both excitatory and inhibitory neurons. We fit a "kernel" function to measurements of uLFPs, and we estimate its spatial and temporal spread by using simulations of morphologically detailed reconstructions of hippocampal pyramidal neurons. Assuming that LFPs are the sum of uLFPs generated by every neuron in the network, the LFP generated by excitatory and inhibitory neurons can be calculated by convolving the trains of action potentials with the kernels estimated from uLFPs. This provides a method to calculate the LFP from networks of spiking neurons, even for point neurons for which the LFP is not easily defined. We show examples of LFPs calculated from networks of point neurons and compare to the LFP calculated from synaptic currents.

CONCLUSIONS

The kernel-based method provides a practical way to calculate LFPs from networks of point neurons.

Experimental and Computational Study on Motor Control and Recovery After Stroke: Toward a Constructive Loop Between Experimental and Virtual Embodied Neuroscience

Being able to replicate real experiments with computational simulations is a unique opportunity to refine and validate models with experimental data and redesign the experiments based on simulations. However, since it is technically demanding to model all components of an experiment, traditional approaches to modeling reduce the experimental setups as much as possible. In this study, our goal is to replicate all the relevant features of an experiment on motor control and motor rehabilitation after stroke. To this aim, we propose an approach that allows continuous integration of new experimental data into a computational modeling framework. First, results show that we could reproduce experimental object displacement with high accuracy via the simulated embodiment in the virtual world by feeding a spinal cord model with experimental registration of the cortical activity. Second, by using computational models of multiple granularities, our preliminary results show the possibility of simulating several features of the brain after stroke, from the local alteration in neuronal activity to long-range connectivity remodeling. Finally, strategies are proposed to merge the two pipelines. We further suggest that additional models could be integrated into the framework thanks to the versatility of the proposed approach, thus allowing many researchers to achieve continuously improved experimental design.

Deep brain stimulation-guided optogenetic rescue of parkinsonian symptoms

Deep brain stimulation (DBS) of the subthalamic nucleus is a symptomatic treatment of Parkinson’s disease but benefits only to a minority of patients due to stringent eligibility criteria. To investigate new targets for less invasive therapies, we aimed at elucidating key mechanisms supporting deep brain stimulation efficiency. Here, using in vivo electrophysiology, optogenetics, behavioral tasks and mathematical modeling, we found that subthalamic stimulation normalizes pathological hyperactivity of motor cortex pyramidal cells, while concurrently activating somatostatin and inhibiting parvalbumin interneurons. In vivo opto-activation of cortical somatostatin interneurons alleviates motor symptoms in a parkinsonian mouse model. A computational model highlights that a decrease in pyramidal neuron activity induced by DBS or by a stimulation of cortical somatostatin interneurons can restore information processing capabilities. Overall, these results demonstrate that activation of cortical somatostatin interneurons may constitute a less invasive alternative than subthalamic stimulation.

Deep brain stimulation (DBS) is a symptomatic treatment of Parkinson’s disease (PD) that benefits only a minority of patients. Here, the authors show that activation of cortical somatostatin interneurons alleviates motor symptoms in a mouse model of PD and may constitute a less invasive alternative than DBS.

Biologically Relevant Dynamical Behaviors Realized in an Ultra-Compact Neuron Model

We demonstrate a variety of biologically relevant dynamical behaviors building on a recently introduced ultra-compact neuron (UCN) model. We provide the detailed circuits which all share a common basic block that realizes the leaky-integrate-and-fire (LIF) spiking behavior. All circuits have a small number of active components and the basic block has only three, two transistors and a silicon controlled rectifier (SCR). We also demonstrate that numerical simulations can faithfully represent the variety of spiking behavior and can be used for further exploration of dynamical behaviors. Taking Izhikevich’s set of biologically relevant behaviors as a reference, our work demonstrates that a circuit of a LIF neuron model can be used as a basis to implement a large variety of relevant spiking patterns. These behaviors may be useful to construct neural networks that can capture complex brain dynamics or may also be useful for artificial intelligence applications. Our UCN model can therefore be considered the electronic circuit counterpart of Izhikevich’s (2003) mathematical neuron model, sharing its two seemingly contradicting features, extreme simplicity and rich dynamical behavior.

Biophysically grounded mean-field models of neural populations under electrical stimulation

Electrical stimulation of neural systems is a key tool for understanding neural dynamics and ultimately for developing clinical treatments. Many applications of electrical stimulation affect large populations of neurons. However, computational models of large networks of spiking neurons are inherently hard to simulate and analyze. We evaluate a reduced mean-field model of excitatory and inhibitory adaptive exponential integrate-and-fire (AdEx) neurons which can be used to efficiently study the effects of electrical stimulation on large neural populations. The rich dynamical properties of this basic cortical model are described in detail and validated using large network simulations. Bifurcation diagrams reflecting the network's state reveal asynchronous up- and down-states, bistable regimes, and oscillatory regions corresponding to fast excitation-inhibition and slow excitation-adaptation feedback loops. The biophysical parameters of the AdEx neuron can be coupled to an electric field with realistic field strengths which then can be propagated up to the population description. We show how on the edge of bifurcation, direct electrical inputs cause network state transitions, such as turning on and off oscillations of the population rate. Oscillatory input can frequency-entrain and phase-lock endogenous oscillations. Relatively weak electric field strengths on the order of 1 V/m are able to produce these effects, indicating that field effects are strongly amplified in the network. The effects of time-varying external stimulation are well-predicted by the mean-field model, further underpinning the utility of low-dimensional neural mass models.

Cortical Responses to Input From Distant Areas are Modulated by Local Spontaneous Alpha/Beta Oscillations

Any given area in human cortex may receive input from multiple, functionally heterogeneous areas, potentially representing different processing threads. Alpha (8-13 Hz) and beta oscillations (13-20 Hz) have been hypothesized by other investigators to gate local cortical processing, but their influence on cortical responses to input from other cortical areas is unknown. To study this, we measured the effect of local oscillatory power and phase on cortical responses elicited by single-pulse electrical stimulation (SPES) at distant cortical sites, in awake human subjects implanted with intracranial electrodes for epilepsy surgery. In 4 out of 5 subjects, the amplitudes of corticocortical evoked potentials (CCEPs) elicited by distant SPES were reproducibly modulated by the power, but not the phase, of local oscillations in alpha and beta frequencies. Specifically, CCEP amplitudes were higher when average oscillatory power just before distant SPES (-110 to -10 ms) was high. This effect was observed in only a subset (0-33%) of sites with CCEPs and, like the CCEPs themselves, varied with stimulation at different distant sites. Our results suggest that although alpha and beta oscillations may gate local processing, they may also enhance the responsiveness of cortex to input from distant cortical sites.

Cortical state transitions and stimulus response evolve along stiff and sloppy parameter dimensions, respectively

Previous research showed that spontaneous neuronal activity presents sloppiness: the collective behavior is strongly determined by a small number of parameter combinations, defined as ‘stiff’ dimensions, while it is insensitive to many others (‘sloppy’ dimensions). Here, we analyzed neural population activity from the auditory cortex of anesthetized rats while the brain spontaneously transited through different synchronized and desynchronized states and intermittently received sensory inputs. We showed that cortical state transitions were determined by changes in stiff parameters associated with the activity of a core of neurons with low responses to stimuli and high centrality within the observed network. In contrast, stimulus-evoked responses evolved along sloppy dimensions associated with the activity of neurons with low centrality and displaying large ongoing and stimulus-evoked fluctuations without affecting the integrity of the network. Our results shed light on the interplay among stability, flexibility, and responsiveness of neuronal collective dynamics during intrinsic and induced activity.

Cholinergic Switch between Two Types of Slow Waves in Cerebral Cortex

Sleep slow waves are known to participate in memory consolidation, yet slow waves occurring under anesthesia present no positive effects on memory. Here, we shed light onto this paradox, based on a combination of extracellular recordings in vivo, in vitro, and computational models. We find two types of slow waves, based on analyzing the temporal patterns of successive slow-wave events. The first type is consistently observed in natural slow-wave sleep, while the second is shown to be ubiquitous under anesthesia. Network models of spiking neurons predict that the two slow wave types emerge due to a different gain on inhibitory versus excitatory cells and that different levels of spike-frequency adaptation in excitatory cells can account for dynamical distinctions between the two types. This prediction was tested in vitro by varying adaptation strength using an agonist of acetylcholine receptors, which demonstrated a neuromodulatory switch between the two types of slow waves. Finally, we show that the first type of slow-wave dynamics is more sensitive to external stimuli, which can explain how slow waves in sleep and anesthesia differentially affect memory consolidation, as well as provide a link between slow-wave dynamics and memory diseases.

Bridging Single Neuron Dynamics to Global Brain States

Biological neural networks produce information backgrounds of multi-scale spontaneous activity that become more complex in brain states displaying higher capacities for cognition, for instance, attentive awake versus asleep or anesthetized states. Here, we review brain state-dependent mechanisms spanning ion channel currents (microscale) to the dynamics of brain-wide, distributed, transient functional assemblies (macroscale). Not unlike how microscopic interactions between molecules underlie structures formed in macroscopic states of matter, using statistical physics, the dynamics of microscopic neural phenomena can be linked to macroscopic brain dynamics through mesoscopic scales. Beyond spontaneous dynamics, it is observed that stimuli evoke collapses of complexity, most remarkable over high dimensional, asynchronous, irregular background dynamics during consciousness. In contrast, complexity may not be further collapsed beyond synchrony and regularity characteristic of unconscious spontaneous activity. We propose that increased dimensionality of spontaneous dynamics during conscious states supports responsiveness, enhancing neural networks' emergent capacity to robustly encode information over multiple scales.

Emergence of spatially periodic diffusive waves in small-world neuronal networks

It has been observed in experiment that the anatomical structure of neuronal networks in the brain possesses the feature of small-world networks. Yet how the small-world structure affects network dynamics remains to be fully clarified. Here we study the dynamics of a class of small-world networks consisting of pulse-coupled integrate-and-fire (I&F) neurons. Under stochastic Poisson drive, we find that the activity of the entire network resembles diffusive waves. To understand its underlying mechanism, we analyze the simplified regular-lattice network consisting of firing-rate-based neurons as an approximation to the original I&F small-world network. We demonstrate both analytically and numerically that, with strongly coupled connections, in the absence of noise, the activity of the firing-rate-based regular-lattice network spatially forms a static grating pattern that corresponds to the spatial distribution of the firing rate observed in the I&F small-world neuronal network. We further show that the spatial grating pattern with different phases comprise the continuous attractor of both the I&F small-world and firing-rate-based regular-lattice network dynamics. In the presence of input noise, the activity of both networks is perturbed along the continuous attractor, which gives rise to the diffusive waves. Our numerical simulations and theoretical analysis may potentially provide insights into the understanding of the generation of wave patterns observed in cortical networks.

Jensen's force and the statistical mechanics of cortical asynchronous states

Cortical networks are shaped by the combined action of excitatory and inhibitory interactions. Among other important functions, inhibition solves the problem of the all-or-none type of response that comes about in purely excitatory networks, allowing the network to operate in regimes of moderate or low activity, between quiescent and saturated regimes. Here, we elucidate a noise-induced effect that we call “Jensen’s force” –stemming from the combined effect of excitation/inhibition balance and network sparsity– which is responsible for generating a phase of self-sustained low activity in excitation-inhibition networks. The uncovered phase reproduces the main empirically-observed features of cortical networks in the so-called asynchronous state, characterized by low, un-correlated and highly-irregular activity. The parsimonious model analyzed here allows us to resolve a number of long-standing issues, such as proving that activity can be self-sustained even in the complete absence of external stimuli or driving. The simplicity of our approach allows for a deep understanding of asynchronous states and of the phase transitions to other standard phases it exhibits, opening the door to reconcile, asynchronous-state and critical-state hypotheses, putting them within a unified framework. We argue that Jensen’s forces are measurable experimentally and might be relevant in contexts beyond neuroscience.