Comparing Epileptiform Behavior of Mesoscale Detailed Models and Population Models of Neocortex

Multi-scale modelling of the epileptic brain: advantages of computational therapy exploration

Objective: Epilepsy is a complex disease spanning across multiple scales, from ion channels in neurons to neuronal circuits across the entire brain. Over the past decades, computational models have been used to describe the pathophysiological activity of the epileptic brain from different aspects. Traditionally, each computational model can aid in optimizing therapeutic interventions, therefore, providing a particular view to design strategies for treating epilepsy. As a result, most studies are concerned with generating specific models of the epileptic brain that can help us understand the certain machinery of the pathological state. Those specific models vary in complexity and biological accuracy, with system-level models often lacking biological details.Approach: Here, we review various types of computational model of epilepsy and discuss their potential for different therapeutic approaches and scenarios, including drug discovery, surgical strategies, brain stimulation, and seizure prediction. We propose that we need to consider an integrated approach with a unified modelling framework across multiple scales to understand the epileptic brain. Our proposal is based on the recent increase in computational power, which has opened up the possibility of unifying those specific epileptic models into simulations with an unprecedented level of detail.Main results: A multi-scale epilepsy model can bridge the gap between biologically detailed models, used to address molecular and cellular questions, and brain-wide models based on abstract models which can account for complex neurological and behavioural observations.Significance: With these efforts, we move toward the next generation of epileptic brain models capable of connecting cellular features, such as ion channel properties, with standard clinical measures such as seizure severity.

Cross-scale effects of neural interactions during human neocortical seizure activity

Small-scale neuronal networks may impose widespread effects on large network dynamics. To unravel this relationship, we analyzed eight multiscale recordings of spontaneous seizures from four patients with epilepsy. During seizures, multiunit spike activity organizes into a submillimeter-sized wavefront, and this activity correlates significantly with low-frequency rhythms from electrocorticographic recordings across a 10-cm-sized neocortical network. Notably, this correlation effect is specific to the ictal wavefront and is absent interictally or from action potential activity outside the wavefront territory. To examine the multiscale interactions, we created a model using a multiscale, nonlinear system and found evidence for a dual role for feedforward inhibition in seizures: while inhibition at the wavefront fails, allowing seizure propagation, feedforward inhibition of the surrounding centimeter-scale networks is activated via long-range excitatory connections. Bifurcation analysis revealed that distinct dynamical pathways for seizure termination depend on the surrounding inhibition strength. Using our model, we found that the mesoscopic, local wavefront acts as the forcing term of the ictal process, while the macroscopic, centimeter-sized network modulates the oscillatory seizure activity.

DYNAMICS OF COLLECTIVE MULTI-STABILITY IN MODELS OF MULTI-UNIT NEURONAL SYSTEMS

In this study, we investigate the correspondence between dynamic patterns of behavior in two types of computational models of neuronal activity. The first model type is the realistic neuronal model; the second model type is the phenomenological or analytical model. In the simplest model set-up of two interconnected units, we define a parameter space for both types of systems where their behavior is similar. Next we expand the analytical model to two sets of 90 fully interconnected units with some overlap, which can display multi-stable behavior. This system can be in three classes of states: (i) a class consisting of a single resting state, where all units of a set are in steady state, (ii) a class consisting of multiple preserving states, where subsets of the units of a set participate in limit cycle, and (iii) a class consisting of a single saturated state, where all units of a set are recruited in a global limit cycle. In the third and final part of the work, we demonstrate that phase synchronization of units can be detected by a single output unit.

Large-scale modeling of epileptic seizures: scaling properties of two parallel neuronal network simulation algorithms

Our limited understanding of the relationship between the behavior of individual neurons and large neuronal networks is an important limitation in current epilepsy research and may be one of the main causes of our inadequate ability to treat it. Addressing this problem directly via experiments is impossibly complex; thus, we have been developing and studying medium-large-scale simulations of detailed neuronal networks to guide us. Flexibility in the connection schemas and a complete description of the cortical tissue seem necessary for this purpose. In this paper we examine some of the basic issues encountered in these multiscale simulations. We have determined the detailed behavior of two such simulators on parallel computer systems. The observed memory and computation-time scaling behavior for a distributed memory implementation were very good over the range studied, both in terms of network sizes (2,000 to 400,000 neurons) and processor pool sizes (1 to 256 processors). Our simulations required between a few megabytes and about 150 gigabytes of RAM and lasted between a few minutes and about a week, well within the capability of most multinode clusters. Therefore, simulations of epileptic seizures on networks with millions of cells should be feasible on current supercomputers.

Modeling Neural Activity

This paper provides an overview of different types of models for studying activity of nerve cells and their networks with a special emphasis on neural oscillations. One part describes the neuronal models based on the Hodgkin and Huxley formalism first described in the 1950s. It is discussed how further simplifications of this formalism enable mathematical analysis of the process of neural excitability. The focus of the paper’s second component is on network activity. Understanding network function is one of the important frontiers remaining in neuroscience. At present, experimental techniques can only provide global recordings or samples of the activity of the huge networks that form the nervous system. Models in neuroscience can therefore play a critical role by providing a framework for integration of necessarily incomplete datasets, thereby providing insight into the mechanisms of neural function. Network models can either explicitly contain individual network nodes that model the neurons, or they can be based on representations of compound population activity. The latter approach was pioneered by Wilson and Cowan in the 1970s. Finally I provide an overview and discuss how network models are employed in the study of neuronal network pathology such as epilepsy.

Analysis of stability and bifurcations of fixed points and periodic solutions of a lumped model of neocortex with two delays

A lumped model of neural activity in neocortex is studied to identify regions of multi-stability of both steady states and periodic solutions. Presence of both steady states and periodic solutions is considered to correspond with epileptogenesis. The model, which consists of two delay differential equations with two fixed time lags is mainly studied for its dependency on varying connection strength between populations. Equilibria are identified, and using linear stability analysis, all transitions are determined under which both trivial and non-trivial fixed points lose stability. Periodic solutions arising at some of these bifurcations are numerically studied with a two-parameter bifurcation analysis.

Vulnerability to paroxysmal oscillations in delayed neural networks: a basis for nocturnal frontal lobe epilepsy?

Resonance can occur in bistable dynamical systems due to the interplay between noise and delay (τ) in the absence of a periodic input. We investigate resonance in a two-neuron model with mutual time-delayed inhibitory feedback. For appropriate choices of the parameters and inputs three fixed-point attractors co-exist: two are stable and one is unstable. In the absence of noise, delay-induced transient oscillations (referred to herein as DITOs) arise whenever the initial function is tuned sufficiently close to the unstable fixed-point. In the presence of noisy perturbations, DITOs arise spontaneously. Since the correlation time for the stationary dynamics is ∼τ, we approximated a higher order Markov process by a three-state Markov chain model by rescaling time as t → 2sτ, identifying the states based on whether the sub-intervals were completely confined to one basin of attraction (the two stable attractors) or straddled the separatrix, and then determining the transition probability matrix empirically. The resultant Markov chain model captured the switching behaviors including the statistical properties of the DITOs. Our observations indicate that time-delayed and noisy bistable dynamical systems are prone to generate DITOs as switches between the two attractors occur. Bistable systems arise transiently in situations when one attractor is gradually replaced by another. This may explain, for example, why seizures in certain epileptic syndromes tend to occur as sleep stages change.