Simulating human sleep spindle MEG and EEG from ion channel and circuit level dynamics

Realistic Subject-Specific Simulation of Resting State Scalp EEG Based on Physiological Model

Electroencephalography (EEG) recordings are widely used in neuroscience to identify healthy individual brain rhythms and to detect alterations associated with various brain diseases. However, understanding the cellular origins of scalp EEG signals and their spatiotemporal changes during the resting state (RS) in humans remains challenging, as cellular-level recordings are typically restricted to animal models. The objective of this study was to simulate individual-specific spatiotemporal features of RS EEG and measure the degree of similarity between real and simulated EEG. Using a physiologically grounded whole-brain computational model (based on known neuronal subtypes and their structural and functional connectivity) that simulates interregional cortical circuitry activity, realistic individual EEG recordings during RS of three healthy subjects were created. The model included interconnected neural mass modules simulating activities of different neuronal subtypes, including pyramidal cells and four types of GABAergic interneurons. High-definition EEG and source localization were used to delineate the cortical extent of alpha and beta-gamma rhythms. To evaluate the realism of the simulated EEG, we developed a similarity index based on cross-correlation analysis in the frequency domain across various bipolar channels respecting standard longitudinal montage. Alpha oscillations were produced by strengthening the somatostatin-pyramidal loop in posterior regions, while beta-gamma oscillations were generated by increasing the excitability of parvalbumin-interneurons on pyramidal neurons in anterior regions. The generation of realistic individual RS EEG rhythms represents a significant advance for research fields requiring data augmentation, including brain-computer interfaces and artificial intelligence training.

A topological method of generating action potentials and electroencephalography oscillations in a surface network

The brain is a source of continuous electrical activity, which includes one-dimensional voltage pulses (action potentials) that propagate along nerve fibres, transient localized oscillations and persistent surface oscillations in five distinct frequency bands. However, a unified theoretical framework for modelling these excitations is lacking. In this article, we provide such a framework by constructing a special surface network in which all observed brain-like signals, including surface oscillations, can be generated by topological means. Analytic expressions for all these excitations are found, and the values of the five frequency bands of surface oscillations are correctly predicted. It is shown how input signals of the system produce their own communication code to encode the information they carry and how the response output propagating signals produced carry this input information with them and can transfer it to the pathways they traverse as a non-transient topological memory structure of aligned spin-half protons. It is conjectured that the memory structure is located in the insulating sheaths of nerve fibres and is stable only if the pathways between the assembly of neurons, which represents a memory structure, include loops. The creation time and size of memory structures are estimated, and a memory-specific excitation frequency for a memory structure is identified and determined, which can be used to recall memories.

Sleep and circadian rhythms modeling: From hypothalamic regulatory networks to cortical dynamics and behavior

Sleep and circadian rhythms are regulated by dynamic physiologic processes that operate across multiple spatial and temporal scales. These include, but are not limited to, genetic oscillators, clearance of waste products from the brain, dynamic interplay among brain regions, and propagation of local dynamics across the cortex. The combination of these processes, modulated by environmental cues, such as light-dark cycles and work schedules, represents a complex multiscale system that regulates sleep-wake cycles and brain dynamics. Physiology-based mathematical models have successfully explained the mechanisms underpinning dynamics at specific scales and are a useful tool to investigate interactions across multiple scales. They can help answer questions such as how do electroencephalographic (EEG) features relate to subthalamic neuron activity? Or how are local cortical dynamics regulated by the homeostatic and circadian mechanisms? In this chapter, we review two types of models that are well-positioned to consider such interactions. Part I of the chapter focuses on the subthalamic sleep regulatory networks and a model of arousal dynamics capable of predicting sleep, circadian rhythms, and cognitive outputs. Part II presents a model of corticothalamic circuits, capable of predicting spatial and temporal EEG features. We then discuss existing approaches and unsolved challenges in developing unified multiscale models.

Python/NEURON code for simulating biophysically realistic thalamocortical dynamics during sleep

Understanding the function of sleep and its associated neural rhythms is an important goal in neuroscience. While many theoretical models of neural dynamics during sleep exist, few include the effects of neuromodulators on sleep oscillations and describe transitions between sleep and wake states or different sleep stages. Here, we started with a C++-based thalamocortical network model that describes characteristic thalamic and cortical oscillations specific to sleep. This model, which includes a biophysically realistic description of intrinsic and synaptic channels, allows for testing the effects of different neuromodulators, intrinsic cell properties, and synaptic connectivity on neural dynamics during sleep. We present a complete reimplementation of this previously-published sleep model in the standardized NEURON/Python framework, making it more accessible to the wider scientific community.

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 (SOs, <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 mechanisms by which global and local SOs arise from micro-scale neuronal dynamics and network connectivity remain poorly understood. We developed a multi-scale, biophysically realistic human whole-brain thalamocortical network model capable of transitioning between the awake state and SWS, 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. An increase in the overall synaptic strength led to synchronized global SO, while a decrease in synaptic connectivity produced only local slow-waves that would not propagate beyond local areas. 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.

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.

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.

Multielectrode Cortical Stimulation Selectively Induces Unidirectional Wave Propagation of Excitatory Neuronal Activity in Biophysical Neural Model

Cortical stimulation is emerging as an experimental tool in basic research and a promising therapy for a range of neuropsychiatric conditions. As multielectrode arrays enter clinical practice, the possibility of using spatiotemporal patterns of electrical stimulation to induce desired physiological patterns has become theoretically possible, but in practice can only be implemented by trial-and-error because of a lack of predictive models. Experimental evidence increasingly establishes traveling waves as fundamental to cortical information-processing, but we lack an understanding of how to control wave properties despite rapidly improving technologies. This study uses a hybrid biophysical-anatomical and neural-computational model to predict and understand how a simple pattern of cortical surface stimulation could induce directional traveling waves via asymmetric activation of inhibitory interneurons. We found that pyramidal cells and basket cells are highly activated by the anodal electrode and minimally activated by the cathodal electrodes, while Martinotti cells are moderately activated by both electrodes but exhibit a slight preference for cathodal stimulation. Network model simulations found that this asymmetrical activation results in a traveling wave in superficial excitatory cells that propagates unidirectionally away from the electrode array. Our study reveals how asymmetric electrical stimulation can easily facilitate traveling waves by relying on two distinct types of inhibitory interneuron activity to shape and sustain the spatiotemporal dynamics of endogenous local circuit mechanisms.SIGNIFICANCE STATEMENT Electrical brain stimulation is becoming increasingly useful to probe the workings of brain and to treat a variety of neuropsychiatric disorders. However, stimulation is currently performed in a trial-and-error fashion as there are no methods to predict how different electrode arrangements and stimulation paradigms will affect brain functioning. In this study, we demonstrate a hybrid modeling approach, which makes experimentally testable predictions that bridge the gap between the microscale effects of multielectrode stimulation and the resultant circuit dynamics at the mesoscale. Our results show how custom stimulation paradigms can induce predictable, persistent changes in brain activity, which has the potential to restore normal brain function and become a powerful therapy for neurological and psychiatric conditions.

Neuronal Circuit-Based Computer Modeling as a Phenotypic Strategy for CNS R&D

With the success rate of drugs for CNS indications at an all-time low, new approaches are needed to turn the tide of failed clinical trials. This paper reviews the history of CNS drug Discovery over the last 60 years and proposes a new paradigm based on the lessons learned. The initial wave of successful therapeutics discovered using careful clinical observations was followed by an emphasis on a phenotypic target-agnostic approach, often leading to successful drugs with a rich pharmacology. The subsequent introduction of molecular biology and the focus on a target-driven strategy has largely dominated drug discovery efforts over the last 30 years, but has not increased the probability of success, because these highly selective molecules are unlikely to address the complex pathological phenotypes of most CNS disorders. In many cases, reliance on preclinical animal models has lacked robust translational power. We argue that Quantitative Systems Pharmacology (QSP), a mechanism-based computer model of biological processes informed by preclinical knowledge and enhanced by neuroimaging and clinical data could be a new powerful knowledge generator engine and paradigm for rational polypharmacy. Progress in the academic discipline of computational neurosciences, allows one to model the effect of pathology and therapeutic interventions on neuronal circuit firing activity that can relate to clinical phenotypes, driven by complex properties of specific brain region activation states. The model is validated by optimizing the correlation between relevant emergent properties of these neuronal circuits and historical clinical and imaging datasets. A rationally designed polypharmacy target profile will be discovered using reverse engineering and sensitivity analysis. Small molecules will be identified using a combination of Artificial Intelligence methods and computational modeling, tested subsequently in heterologous cellular systems with human targets. Animal models will be used to establish target engagement and for ADME-Tox, with the QSP approach complemented by in vivo preclinical models that can be further refined to increase predictive validity. The QSP platform can also mitigate the variability in clinical trials with the concept of virtual patients. Because the QSP platform integrates knowledge from a wide variety of sources in an actionable simulation, it offers the possibility of substantially improving the success rate of CNS R&D programs while, at the same time, reducing both cost and the number of animals.