Computer modelling of epilepsy

Small changes in synaptic gain lead to seizure-like activity in neuronal network at criticality

Epilepsy is a neurological disorder characterised by spontaneous recurrent seizures. The mechanisms by which multiple molecular and cellular changes lead to seizures is not well understood. Here, we study cortical seizure generation by simulating the activity of neuron groups in a network using the laminar cortex model. We identified a clear boundary between low-amplitude, asynchronous activity and high-amplitude, rhythmic activity, around which small changes in excitatory synaptic gain led to strong oscillatory activity. Neuron groups only responded significantly to stimulation around the boundary. The consequences of biophysical changes induced by epilepsy-related SCN1A mutations were also examined. Marked reduction in neuronal inhibition, as caused by mutations underlying Dravet syndrome, invariably led to strong neuronal firing, whereas small reductions in inhibition could cause significant changes when the network was poised close to the boundary. The study highlights the critical role of network dynamics in seizure genesis.

A Proposed Mechanism for Spontaneous Transitions between Interictal and Ictal Activity

Epileptic networks are characterized by two outputs: brief interictal spikes and rarer, more prolonged seizures. Although either output state is readily modeled in silico and induced experimentally, the transition mechanisms are unknown, in part because no models exhibit both output states spontaneously. In silico small-world neural networks were built using single-compartment neurons whose physiological parameters were derived from dual whole-cell recordings of pyramidal cells in organotypic hippocampal slice cultures that were generating spontaneous seizure-like activity. In silico, neurons were connected by abundant local synapses and rare long-distance synapses. Activity-dependent synaptic depression and gradual recovery delimited synchronous activity. Full synaptic recovery engendered interictal population spikes that spread via long-distance synapses. When synaptic recovery was incomplete, postsynaptic neurons required coincident activation of multiple presynaptic terminals to reach firing threshold. Only local connections were sufficiently dense to spread activity under these conditions. This coalesced network activity into traveling waves whose velocity varied with synaptic recovery. Seizures were comprised of sustained traveling waves that were similar to those recorded during experimental and human neocortical seizures. Sustained traveling waves occurred only when wave velocity, network dimensions, and the rate of synaptic recovery enabled wave reentry into previously depressed areas at precisely ictogenic levels of synaptic recovery. Wide-field, cellular-resolution GCamP7b calcium imaging demonstrated similar initial patterns of activation in the hippocampus, although the anatomical distribution of traveling waves of synaptic activation was altered by the pattern of synaptic connectivity in the organotypic hippocampal cultures.SIGNIFICANCE STATEMENT When computerized distributed neural network models are required to generate both features of epileptic networks (i.e., spontaneous interictal population spikes and seizures), the network structure is substantially constrained. These constraints provide important new hypotheses regarding the nature of epileptic networks and mechanisms of seizure onset.

Computer modelling of connectivity change suggests epileptogenesis mechanisms in idiopathic generalised epilepsy

Patients with idiopathic generalised epilepsy (IGE) typically have normal conventional magnetic resonance imaging (MRI), hence diagnosis based on MRI is challenging. Anatomical abnormalities underlying brain dysfunctions in IGE are unclear and their relation to the pathomechanisms of epileptogenesis is poorly understood.

In this study, we applied connectometry, an advanced quantitative neuroimaging technique for investigating localised changes in white-matter tissues in vivo. Analysing white matter structures of 32 subjects we incorporated our in vivo findings in a computational model of seizure dynamics to suggest a plausible mechanism of epileptogenesis.

Patients with IGE have significant bilateral alterations in major white-matter fascicles. In the cingulum, fornix, and superior longitudinal fasciculus, tract integrity is compromised, whereas in specific parts of tracts between thalamus and the precentral gyrus, tract integrity is enhanced in patients. Combining these alterations in a logistic regression model, we computed the decision boundary that discriminated patients and controls. The computational model, informed with the findings on the tract abnormalities, specifically highlighted the importance of enhanced cortico-reticular connections along with impaired cortico-cortical connections in inducing pathological seizure-like dynamics.

We emphasise taking directionality of brain connectivity into consideration towards understanding the pathological mechanisms; this is possible by combining neuroimaging and computational modelling. Our imaging evidence of structural alterations suggest the loss of cortico-cortical and enhancement of cortico-thalamic fibre integrity in IGE. We further suggest that impaired connectivity from cortical regions to the thalamic reticular nucleus offers a therapeutic target for selectively modifying the brain circuit for reversing the mechanisms leading to epileptogenesis.

•

Significant focal alterations along major white-matter fascicles in IGE patients are characterised.

•

Increased white matter integrity found in thalamo-cortical connections.

•

Decreased white matter integrity found in cortico-cortical connections.

•

Disease mechanism is investigated by combining the neuroimaging findings with a dynamical model of seizure activity.

•

Model implicates cortical projections to the thalamic reticular nucleus in IGE.

Trkb-IP3 Pathway Mediating Neuroprotection in Rat Hippocampal Neuronal Cell Culture Following Induction of Kainic Acid

BACKGROUND

Following brain injury, development of hippocampal sclerosis often led to the temporal lobe epilepsy which is sometimes resistant to common anti-epileptic drugs. Cellular and molecular changes underlying epileptogenesis in animal models were studied, however, the underlying mechanisms of kainic acid (KA) mediated neuronal damage in rat hippocampal neuron cell culture alone has not been elucidated yet.

METHODS

Embryonic day 18 (E-18) rat hippocampus neurons were cultured with poly-L-lysine coated glass coverslips. Following optimisation, KA (0.5 μM), a chemoconvulsant agent, was administered at three different time-points (30, 60 and 90 min) to induce seizure in rat hippocampal neuronal cell culture. We examined cell viability, neurite outgrowth density and immunoreactivity of the hippocampus neuron culture by measuring brain derived neurotrophic factor (BDNF), γ-amino butyric acid A (GABAA) subunit α-1 (GABRA1), tyrosine receptor kinase B (TrkB), and inositol trisphosphate receptor (IP3R/IP3) levels.

RESULTS

The results revealed significantly decreased and increased immunoreactivity changes in TrkB (a BDNF receptor) and IP3R, respectively, at 60 min time point.

CONCLUSION

The current findings suggest that TrkB and IP3 could have a neuroprotective role which could be a potential pharmacological target for anti-epilepsy drugs.

NMDA-receptor antibodies alter cortical microcircuit dynamics

Recently, autoantibodies against NMDA receptors (NMDARs) were identified as a major cause of autoimmune encephalitis. They cause abnormalities in brain function often associated with significant changes in patients’ brain dynamics. Here we use computational modeling to identify how NMDAR dysfunction causes abnormalities in brain dynamics using patient EEGs and local field potential recordings in a mouse model of NMDAR-Ab encephalitis. NMDAR autoantibodies cause a specific shift in excitatory coupling within cortical circuits that places the circuits closer to pathological transitions between dynamic brain states. Because of the proximity to these phase transitions, otherwise benign fluctuations in neuronal coupling cause abnormal EEG responses in the presence of the antibodies. Our modeling results thus explain fluctuating abnormalities in brain dynamics observed in patients.

NMDA-receptor antibodies (NMDAR-Abs) cause an autoimmune encephalitis with a diverse range of EEG abnormalities. NMDAR-Abs are believed to disrupt receptor function, but how blocking this excitatory synaptic receptor can lead to paroxysmal EEG abnormalities—or even seizures—is poorly understood. Here we show that NMDAR-Abs change intrinsic cortical connections and neuronal population dynamics to alter the spectral composition of spontaneous EEG activity and predispose brain dynamics to paroxysmal abnormalities. Based on local field potential recordings in a mouse model, we first validate a dynamic causal model of NMDAR-Ab effects on cortical microcircuitry. Using this model, we then identify the key synaptic parameters that best explain EEG paroxysms in pediatric patients with NMDAR-Ab encephalitis. Finally, we use the mouse model to show that NMDAR-Ab–related changes render microcircuitry critically susceptible to overt EEG paroxysms when these key parameters are changed, even though the same parameter fluctuations are tolerated in the in silico model of the control condition. These findings offer mechanistic insights into circuit-level dysfunction induced by NMDAR-Ab.

Seizure evolution can be characterized as path through synaptic gain space of a neural mass model

Physiologically based models could facilitate better understanding of mechanisms underlying epileptic seizures. In this paper, we attempt to reveal the dynamic evolution of intracranial EEG activity during epileptic seizures based on synaptic gain identification procedure of a neural mass model. The distribution of average excitatory, slow and fast inhibitory synaptic gain in the parameter space and their temporal evolution, i.e., the path through the model parameter space, were analyzed in thirty seizures from ten temporal lobe epileptic patients. Results showed that the synaptic gain values located roughly on a plane before seizure onset, dispersed during seizure and returned to the plane when seizure terminated. Cluster analysis was performed on seizure paths and demonstrated consistency in synaptic gain evolution across different seizures from the individual patient. Furthermore, two patient groups were identified, each one corresponding to a specific synaptic gain evolution in the parameter space during a seizure. Results were validated by a bootstrapping approach based on comparison with random paths. The differences in the path revealed variations in EEG dynamics for patients despite showing identical seizure onset pattern. Our approach may have the potential to classify the epileptic patients into subgroups based on different mechanisms revealed by subtle changes in synaptic gains and further enable more robust decisions regarding treatment strategy.

Adaptation and Inhibition Control Pathological Synchronization in a Model of Focal Epileptic Seizure

Pharmacoresistant epilepsy is a common neurological disorder in which increased neuronal intrinsic excitability and synaptic excitation lead to pathologically synchronous behavior in the brain. In the majority of experimental and theoretical epilepsy models, epilepsy is associated with reduced inhibition in the pathological neural circuits, yet effects of intrinsic excitability are usually not explicitly analyzed. Here we present a novel neural mass model that includes intrinsic excitability in the form of spike-frequency adaptation in the excitatory population. We validated our model using local field potential (LFP) data recorded from human hippocampal/subicular slices. We found that synaptic conductances and slow adaptation in the excitatory population both play essential roles for generating seizures and pre-ictal oscillations. Using bifurcation analysis, we found that transitions towards seizure and back to the resting state take place via Andronov-Hopf bifurcations. These simulations therefore suggest that single neuron adaptation as well as synaptic inhibition are responsible for orchestrating seizure dynamics and transition towards the epileptic state.

Calcium imaging and dynamic causal modelling reveal brain-wide changes in effective connectivity and synaptic dynamics during epileptic seizures

Pathophysiological explanations of epilepsy typically focus on either the micro/mesoscale (e.g. excitation-inhibition imbalance), or on the macroscale (e.g. network architecture). Linking abnormalities across spatial scales remains difficult, partly because of technical limitations in measuring neuronal signatures concurrently at the scales involved. Here we use light sheet imaging of the larval zebrafish brain during acute epileptic seizure induced with pentylenetetrazole. Spectral changes of spontaneous neuronal activity during the seizure are then modelled using neural mass models, allowing Bayesian inference on changes in effective network connectivity and their underlying synaptic dynamics. This dynamic causal modelling of seizures in the zebrafish brain reveals concurrent changes in synaptic coupling at macro- and mesoscale. Fluctuations of both synaptic connection strength and their temporal dynamics are required to explain observed seizure patterns. These findings highlight distinct changes in local (intrinsic) and long-range (extrinsic) synaptic transmission dynamics as a possible seizure pathomechanism and illustrate how our Bayesian model inversion approach can be used to link existing neural mass models of seizure activity and novel experimental methods.

Dependence of absence seizure dynamics on physiological parameter evolution

A neural field model of the corticothalamic system is applied to investigate the temporal and spectral characteristics of absence seizures in the presence of a temporally varying connection strength between the cerebral cortex and thalamus. Increasing connection strength drives the system into an absence seizure-like state once a threshold is passed and a supercritical Hopf bifurcation occurs. The dynamics and spectral characteristics of the resulting model seizures are explored as functions of maximum connection strength, time above threshold, and the rate at which the connection strength increases (ramp rate). Our results enable spectral and temporal characteristics of seizures to be related to changes in the underlying physiological evolution of connections via nonlinear dynamics and neural field theory. Spectral analysis reveals that the power of the harmonics and the duration of the oscillations increase as the maximum connection strength and the time above threshold increase. It is also found that the time to reach the stable limit-cycle seizure oscillation from the instability threshold decreases with the square root of the ramp rate.

Activity-induced spontaneous spikes in GABAergic neurons suppress seizure discharges: an implication of computational modeling

Background

Epilepsy, a prevalent neurological disorder, appears self-termination. The endogenous mechanism for seizure self-termination remains to be addressed in order to develop new strategies for epilepsy treatment. We aim to examine the role of activity-induced spontaneous spikes at GABAergic neurons as an endogenous mechanism in the seizure self-termination.

Methods and Results

Neuronal spikes were induced by depolarization pulses at cortical GABAergic neurons from temporal lobe epilepsy patients and mice, in which some of these neurons fired activity-induced spontaneous spikes. Neural networks including excitatory and inhibitory neurons were computationally constructed, and their functional properties were based on our studies from whole-cell recordings. With the changes in the portion and excitability of inhibitory neurons that generated activity-induced spontaneous spike, the efficacies to suppress synchronous seizure activity were analyzed, such as its onset time, decay slope and spike frequency. The increases in the proportion and excitability of inhibitory neurons that generated activity-induced spontaneous spikes effectively suppressed seizure activity in neural networks. These factors synergistically strengthened the efficacy of seizure activity suppression.

Conclusion

Our study supports a notion that activity-induced spontaneous spikes in GABAergic neurons may be an endogenous mechanism for seizure self-termination. A potential therapeutic strategy for epilepsy is to upregulate the cortical inhibitory neurons that generate activity-induced spontaneous spikes.

Credibility, Replicability, and Reproducibility in Simulation for Biomedicine and Clinical Applications in Neuroscience

Modeling and simulation in computational neuroscience is currently a research enterprise to better understand neural systems. It is not yet directly applicable to the problems of patients with brain disease. To be used for clinical applications, there must not only be considerable progress in the field but also a concerted effort to use best practices in order to demonstrate model credibility to regulatory bodies, to clinics and hospitals, to doctors, and to patients. In doing this for neuroscience, we can learn lessons from long-standing practices in other areas of simulation (aircraft, computer chips), from software engineering, and from other biomedical disciplines. In this manuscript, we introduce some basic concepts that will be important in the development of credible clinical neuroscience models: reproducibility and replicability; verification and validation; model configuration; and procedures and processes for credible mechanistic multiscale modeling. We also discuss how garnering strong community involvement can promote model credibility. Finally, in addition to direct usage with patients, we note the potential for simulation usage in the area of Simulation-Based Medical Education, an area which to date has been primarily reliant on physical models (mannequins) and scenario-based simulations rather than on numerical simulations.

Analytically determining frequency and amplitude of spontaneous alpha oscillation in Jansen's neural mass model using the describing function method

Spontaneous alpha oscillations are a ubiquitous phenomenon in the brain and play a key role in neural information processing and various cognitive functions. Jansen's neural mass model (NMM) was initially proposed to study the origin of alpha oscillations. Most of previous studies of the spontaneous alpha oscillations in the NMM were conducted using numerical methods. In this study, we aim to propose an analytical approach using the describing function method to elucidate the spontaneous alpha oscillation mechanism in the NMM. First, the sigmoid nonlinear function in the NMM is approximated by its describing function, allowing us to reformulate the NMM and derive its standard form composed of one nonlinear part and one linear part. Second, by conducting a theoretical analysis, we can assess whether or not the spontaneous alpha oscillation would occur in the NMM and, furthermore, accurately determine its amplitude and frequency. The results reveal analytically that the interaction between linearity and nonlinearity of the NMM plays a key role in generating the spontaneous alpha oscillations. Furthermore, strong nonlinearity and large linear strength are required to generate the spontaneous alpha oscillations.

Modeling pathogenesis and treatment response in childhood absence epilepsy

OBJECTIVE

Childhood absence epilepsy (CAE) is a genetic generalized epilepsy syndrome with polygenic inheritance, with genes for γ-aminobutyric acid (GABA) receptors and T-type calcium channels implicated in the disorder. Previous studies of T-type calcium channel electrophysiology have shown genetic changes and medications have multiple effects. The aim of this study was to use an established thalamocortical computer model to determine how T-type calcium channels work in concert with cortical excitability to contribute to pathogenesis and treatment response in CAE.

METHODS

The model is comprised of cortical pyramidal, cortical inhibitory, thalamocortical relay, and thalamic reticular single-compartment neurons, implemented with Hodgkin-Huxley model ion channels and connected by AMPA, GABA_A , and GABA_B synapses. Network behavior was simulated for different combinations of T-type calcium channel conductance, inactivation time, steady state activation/inactivation shift, and cortical GABA_A conductance.

RESULTS

Decreasing cortical GABA_A conductance and increasing T-type calcium channel conductance converted spindle to spike and wave oscillations; smaller changes were required if both were changed in concert. In contrast, left shift of steady state voltage activation/inactivation did not lead to spike and wave oscillations, whereas right shift reduced network propensity for oscillations of any type.

SIGNIFICANCE

These results provide a window into mechanisms underlying polygenic inheritance in CAE, as well as a mechanism for treatment effects and failures mediated by these channels. Although the model is a simplification of the human thalamocortical network, it serves as a useful starting point for predicting the implications of ion channel electrophysiology in polygenic epilepsy such as CAE.

Multiscale modeling in the clinic: diseases of the brain and nervous system

Computational neuroscience is a field that traces its origins to the efforts of Hodgkin and Huxley, who pioneered quantitative analysis of electrical activity in the nervous system. While also continuing as an independent field, computational neuroscience has combined with computational systems biology, and neural multiscale modeling arose as one offshoot. This consolidation has added electrical, graphical, dynamical system, learning theory, artificial intelligence and neural network viewpoints with the microscale of cellular biology (neuronal and glial), mesoscales of vascular, immunological and neuronal networks, on up to macroscales of cognition and behavior. The complexity of linkages that produces pathophysiology in neurological, neurosurgical and psychiatric disease will require multiscale modeling to provide understanding that exceeds what is possible with statistical analysis or highly simplified models: how to bring together pharmacotherapeutics with neurostimulation, how to personalize therapies, how to combine novel therapies with neurorehabilitation, how to interlace periodic diagnostic updates with frequent reevaluation of therapy, how to understand a physical disease that manifests as a disease of the mind. Multiscale modeling will also help to extend the usefulness of animal models of human diseases in neuroscience, where the disconnects between clinical and animal phenomenology are particularly pronounced. Here we cover areas of particular interest for clinical application of these new modeling neurotechnologies, including epilepsy, traumatic brain injury, ischemic disease, neurorehabilitation, drug addiction, schizophrenia and neurostimulation.

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.

Sample-based approach can outperform the classical dynamical analysis - experimental confirmation of the basin stability method

In this paper we show the first broad experimental confirmation of the basin stability approach. The basin stability is one of the sample-based approach methods for analysis of the complex, multidimensional dynamical systems. We show that investigated method is a reliable tool for the analysis of dynamical systems and we prove that it has a significant advantages which make it appropriate for many applications in which classical analysis methods are difficult to apply. We study theoretically and experimentally the dynamics of a forced double pendulum. We examine the ranges of stability for nine different solutions of the system in a two parameter space, namely the amplitude and the frequency of excitation. We apply the path-following and the extended basin stability methods (Brzeski et al., Meccanica 51(11), 2016) and we verify obtained theoretical results in experimental investigations. Comparison of the presented results show that the sample-based approach offers comparable precision to the classical method of analysis. However, it is much simpler to apply and can be used despite the type of dynamical system and its dimensions. Moreover, the sample-based approach has some unique advantages and can be applied without the precise knowledge of parameter values.

Coupling conditions for globally stable and robust synchrony of chaotic systems

We propose a set of general coupling conditions to select a coupling profile (a set of coupling matrices) from the linear flow matrix of dynamical systems for realizing global stability of complete synchronization (CS) in identical systems and robustness to parameter perturbation. The coupling matrices define the coupling links between any two oscillators in a network that consists of a conventional diffusive coupling link (self-coupling link) as well as a cross-coupling link. The addition of a selective cross-coupling link in particular plays constructive roles that ensure the global stability of synchrony and furthermore enables robustness of synchrony against small to nonsmall parameter perturbation. We elaborate the general conditions for the selection of coupling profiles for two coupled systems, three- and four-node network motifs analytically as well as numerically using benchmark models, the Lorenz system, the Hindmarsh-Rose neuron model, the Shimizu-Morioka laser model, the Rössler system, and a Sprott system. The role of the cross-coupling link is, particularly, exemplified with an example of a larger network, where it saves the network from a breakdown of synchrony against large parameter perturbation in any node. The perturbed node in the network transits from CS to generalized synchronization (GS) when all the other nodes remain in CS. The GS is manifested by an amplified response of the perturbed node in a coherent state.

Computer modeling of epilepsy: opportunities for drug discovery

Analysis of the brain as a dynamical system can assist drug development for dynamical diseases such as epilepsy. The pathological trajectories that make up a seizure differ significantly from the physiological trajectories of normal brain function. These trajectories depend on parameters - conductances and time constants of ion channels and synapses - that can be modified by drugs. Drug development will benefit by taking account of the way in which multiple parameters - multiple drug targets - produce trajectory alterations. This may lead us to reconsider potential benefits of multi-target polypharmacy, of drug cocktails, and of so-called "dirty drugs" (drugs with activity at multiple locations).

Predicting epileptic seizures from scalp EEG based on attractor state analysis

BACKGROUND AND OBJECTIVE

Epilepsy is the second most common disease of the brain. Epilepsy makes it difficult for patients to live a normal life because it is difficult to predict when seizures will occur. In this regard, if seizures could be predicted a reasonable period of time before their occurrence, epilepsy patients could take precautions against them and improve their safety and quality of life. In this paper, we investigate a novel seizure precursor based on attractor state analysis for seizure prediction.

METHODS

We analyze the transition process from normal to seizure attractor state and investigate a precursor phenomenon seen before reaching the seizure attractor state. From the result of an analysis, we define a quantified spectral measure in scalp EEG for seizure prediction. From scalp EEG recordings, the Fourier coefficients of six EEG frequency bands are extracted, and the defined spectral measure is computed based on the coefficients for each half-overlapped 20-second-long window. The computed spectral measure is applied to seizure prediction using a low-complexity methodology.

RESULTS

Within scalp EEG, we identified an early-warning indicator before an epileptic seizure occurs. Getting closer to the bifurcation point that triggers the transition from normal to seizure state, the power spectral density of low frequency bands of the perturbation of an attractor in the EEG, showed a relative increase. A low-complexity seizure prediction algorithm using this feature was evaluated, using ∼583h of scalp EEG in which 143 seizures in 16 patients were recorded. With the test dataset, the proposed method showed high sensitivity (86.67%) with a false prediction rate of 0.367h^-1 and average prediction time of 45.3min.

CONCLUSIONS

A novel seizure prediction method using scalp EEG, based on attractor state analysis, shows potential for application with real epilepsy patients. This is the first study in which the seizure-precursor phenomenon of an epileptic seizure is investigated based on attractor-based analysis of the macroscopic dynamics of the brain. With the scalp EEG, we first propose use of a spectral feature identified for seizure prediction, in which the dynamics of an attractor are excluded, and only the perturbation dynamics from the attractor are considered.

Mechanisms underlying different onset patterns of focal seizures

Focal seizures are episodes of pathological brain activity that appear to arise from a localised area of the brain. The onset patterns of focal seizure activity have been studied intensively, and they have largely been distinguished into two types-low amplitude fast oscillations (LAF), or high amplitude spikes (HAS). Here we explore whether these two patterns arise from fundamentally different mechanisms. Here, we use a previously established computational model of neocortical tissue, and validate it as an adequate model using clinical recordings of focal seizures. We then reproduce the two onset patterns in their most defining properties and investigate the possible mechanisms underlying the different focal seizure onset patterns in the model. We show that the two patterns are associated with different mechanisms at the spatial scale of a single ECoG electrode. The LAF onset is initiated by independent patches of localised activity, which slowly invade the surrounding tissue and coalesce over time. In contrast, the HAS onset is a global, systemic transition to a coexisting seizure state triggered by a local event. We find that such a global transition is enabled by an increase in the excitability of the "healthy" surrounding tissue, which by itself does not generate seizures, but can support seizure activity when incited. In our simulations, the difference in surrounding tissue excitability also offers a simple explanation of the clinically reported difference in surgical outcomes. Finally, we demonstrate in the model how changes in tissue excitability could be elucidated, in principle, using active stimulation. Taken together, our modelling results suggest that the excitability of the tissue surrounding the seizure core may play a determining role in the seizure onset pattern, as well as in the surgical outcome.

Control of Absence Seizures by the Thalamic Feed-Forward Inhibition

As a subtype of idiopathic generalized epilepsies, absence epilepsy is believed to be caused by pathological interactions within the corticothalamic (CT) system. Using a biophysical mean-field model of the CT system, we demonstrate here that the feed-forward inhibition (FFI) in thalamus, i.e., the pathway from the cerebral cortex (Ctx) to the thalamic reticular nucleus (TRN) and then to the specific relay nuclei (SRN) of thalamus that are also directly driven by the Ctx, may participate in controlling absence seizures. In particular, we show that increasing the excitatory Ctx-TRN coupling strength can significantly suppress typical electrical activities during absence seizures. Further, investigation demonstrates that the GABA_A- and GABA_B-mediated inhibitions in the TRN-SRN pathway perform combination roles in the regulation of absence seizures. Overall, these results may provide an insightful mechanistic understanding of how the thalamic FFI serves as an intrinsic regulator contributing to the control of absence seizures.

Understanding Epileptiform After-Discharges as Rhythmic Oscillatory Transients

Electro-cortical activity in patients with epilepsy may show abnormal rhythmic transients in response to stimulation. Even when using the same stimulation parameters in the same patient, wide variability in the duration of transient response has been reported. These transients have long been considered important for the mapping of the excitability levels in the epileptic brain but their dynamic mechanism is still not well understood. To investigate the occurrence of abnormal transients dynamically, we use a thalamo-cortical neural population model of epileptic spike-wave activity and study the interaction between slow and fast subsystems. In a reduced version of the thalamo-cortical model, slow wave oscillations arise from a fold of cycles (FoC) bifurcation. This marks the onset of a region of bistability between a high amplitude oscillatory rhythm and the background state. In vicinity of the bistability in parameter space, the model has excitable dynamics, showing prolonged rhythmic transients in response to suprathreshold pulse stimulation. We analyse the state space geometry of the bistable and excitable states, and find that the rhythmic transient arises when the impending FoC bifurcation deforms the state space and creates an area of locally reduced attraction to the fixed point. This area essentially allows trajectories to dwell there before escaping to the stable steady state, thus creating rhythmic transients. In the full thalamo-cortical model, we find a similar FoC bifurcation structure. Based on the analysis, we propose an explanation of why stimulation induced epileptiform activity may vary between trials, and predict how the variability could be related to ongoing oscillatory background activity. We compare our dynamic mechanism with other mechanisms (such as a slow parameter change) to generate excitable transients, and we discuss the proposed excitability mechanism in the context of stimulation responses in the epileptic cortex.

Basin stability measure of different steady states in coupled oscillators

In this report, we investigate the stabilization of saddle fixed points in coupled oscillators where individual oscillators exhibit the saddle fixed points. The coupled oscillators may have two structurally different types of suppressed states, namely amplitude death and oscillation death. The stabilization of saddle equilibrium point refers to the amplitude death state where oscillations are ceased and all the oscillators converge to the single stable steady state via inverse pitchfork bifurcation. Due to multistability features of oscillation death states, linear stability theory fails to analyze the stability of such states analytically, so we quantify all the states by basin stability measurement which is an universal nonlocal nonlinear concept and it interplays with the volume of basins of attractions. We also observe multi-clustered oscillation death states in a random network and measure them using basin stability framework. To explore such phenomena we choose a network of coupled Duffing-Holmes and Lorenz oscillators which are interacting through mean-field coupling. We investigate how basin stability for different steady states depends on mean-field density and coupling strength. We also analytically derive stability conditions for different steady states and confirm by rigorous bifurcation analysis.

Computer modeling for pharmacological treatments for dystonia

Dystonia is a movement disorder that produces involuntary muscle contractions. Current pharmacological treatments are of limited efficacy. Dystonia, like epilepsy is a disorder involving excessive activty of motor areas including motor cortex and several causal gene mutations have been identified. In order to evaluate potential novel agents for multitarget therapy for dystonia, we have developed a computer model of cortex that includes some of the complex array of molecular interactions that, along with membrane ion channels, control cell excitability.

Synchronization dynamics on the picosecond time scale in coupled Josephson junction neurons

Conventional digital computation is rapidly approaching physical limits for speed and energy dissipation. Here we fabricate and test a simple neuromorphic circuit that models neuronal somas, axons, and synapses with superconducting Josephson junctions. The circuit models two mutually coupled excitatory neurons. In some regions of parameter space the neurons are desynchronized. In others, the Josephson neurons synchronize in one of two states, in-phase or antiphase. An experimental alteration of the delay and strength of the connecting synapses can toggle the system back and forth in a phase-flip bifurcation. Firing synchronization states are calculated >70 000 times faster than conventional digital approaches. With their speed and low energy dissipation (10^{-17}J/spike), this set of proof-of-concept experiments establishes Josephson junction neurons as a viable approach for improvements in neuronal computation as well as applications in neuromorphic computing.

KCC2 downregulation facilitates epileptic seizures

GABAA receptor-mediated inhibition depends on the maintenance of low level intracellular [Cl-] concentration, which in adult depends on neuron specific K+-Cl- cotransporter-2 (KCC2). Previous studies have shown that KCC2 was downregulated in both epileptic patients and various epileptic animal models. However, the temporal relationship between KCC2 downregulation and seizure induction is unclear yet. In this study, we explored the temporal relationship and the influence of KCC2 downregulation on seizure induction. Significant downregulation of plasma membrane KCC2 was directly associated with severe (Racine Score III and above) behavioral seizures in vivo, and occurred before epileptiform bursting activities in vitro induced by convulsant. Overexpression of KCC2 using KCC2 plasmid effectively enhanced resistance to convulsant-induced epileptiform bursting activities in vitro. Furthermore, suppression of membrane KCC2 expression, using shRNAKCC2 plasmid in vitro and shRNAKCC2 containing lentivirus in vivo, induced spontaneous epileptiform bursting activities in vitro and Racine III seizure behaviors accompanied by epileptic EEG in vivo. Our findings novelly demonstrated that altered expression of KCC2 is not the consequence of seizure occurrence but likely is the contributing factor.

Analysis and application of neuronal network controllability and observability

Controllability and observability analyses are important prerequisite for designing suitable neural control strategy, which can help lower the efforts required to control and observe the system dynamics. First, 3-neuron motifs including the excitatory motif, the inhibitory motif, and the mixed motif are constructed to investigate the effects of single neuron and synaptic dynamics on network controllability (observability). Simulation results demonstrate that for networks with the same topological structure, the controllability (observability) of the node always changes if the properties of neurons and synaptic coupling strengths vary. Besides, the inhibitory networks are more controllable (observable) than the excitatory networks when the coupling strengths are the same. Then, the numerically determined controllability results of 3-neuron excitatory motifs are generalized to the desynchronization control of the modular motif network. The control energy and neuronal synchrony measure indexes are used to quantify the controllability of each node in the modular network. The best driver node obtained in this way is the same as the deduced one from motif analysis.

Temporally correlated fluctuations drive epileptiform dynamics

Macroscopic models of brain networks typically incorporate assumptions regarding the characteristics of afferent noise, which is used to represent input from distal brain regions or ongoing fluctuations in non-modelled parts of the brain. Such inputs are often modelled by Gaussian white noise which has a flat power spectrum. In contrast, macroscopic fluctuations in the brain typically follow a 1/fb spectrum. It is therefore important to understand the effect on brain dynamics of deviations from the assumption of white noise. In particular, we wish to understand the role that noise might play in eliciting aberrant rhythms in the epileptic brain. To address this question we study the response of a neural mass model to driving by stochastic, temporally correlated input. We characterise the model in terms of whether it generates "healthy" or "epileptiform" dynamics and observe which of these dynamics predominate under different choices of temporal correlation and amplitude of an Ornstein-Uhlenbeck process. We find that certain temporal correlations are prone to eliciting epileptiform dynamics, and that these correlations produce noise with maximal power in the δ and θ bands. Crucially, these are rhythms that are found to be enhanced prior to seizures in humans and animal models of epilepsy. In order to understand why these rhythms can generate epileptiform dynamics, we analyse the response of the model to sinusoidal driving and explain how the bifurcation structure of the model gives rise to these findings. Our results provide insight into how ongoing fluctuations in brain dynamics can facilitate the onset and propagation of epileptiform rhythms in brain networks. Furthermore, we highlight the need to combine large-scale models with noise of a variety of different types in order to understand brain (dys-)function.

Reduced Efficacy of the KCC2 Cotransporter Promotes Epileptic Oscillations in a Subiculum Network Model

Pharmacoresistant epilepsy is a chronic neurological condition in which a basal brain hyperexcitability results in paroxysmal hypersynchronous neuronal discharges. Human temporal lobe epilepsy has been associated with dysfunction or loss of the potassium-chloride cotransporter KCC2 in a subset of pyramidal cells in the subiculum, a key structure generating epileptic activities. KCC2 regulates intraneuronal chloride and extracellular potassium levels by extruding both ions. Absence of effective KCC2 may alter the dynamics of chloride and potassium levels during repeated activation of GABAergic synapses due to interneuron activity. In turn, such GABAergic stress may itself affect Cl- regulation. Such changes in ionic homeostasis may switch GABAergic signaling from inhibitory to excitatory in affected pyramidal cells and also increase neuronal excitability. Possibly these changes contribute to periodic bursting in pyramidal cells, an essential component in the onset of ictal epileptic events. We tested this hypothesis with a computational model of a subicular network with realistic connectivity. The pyramidal cell model explicitly incorporated the cotransporter KCC2 and its effects on the internal/external chloride and potassium levels. Our network model suggested the loss of KCC2 in a critical number of pyramidal cells increased external potassium and intracellular chloride concentrations leading to seizure-like field potential oscillations. These oscillations included transient discharges leading to ictal-like field events with frequency spectra as in vitro Restoration of KCC2 function suppressed seizure activity and thus may present a useful therapeutic option. These simulations therefore suggest that reduced KCC2 cotransporter activity alone may underlie the generation of ictal discharges.

SIGNIFICANCE STATEMENT

Ion regulation in the brain is a major determinant of neural excitability. Intracellular chloride in neurons, a partial determinant of the resting potential and the inhibitory reversal potentials, is regulated together with extracellular potassium via kation chloride cotransporters. During temporal lobe epilepsy, the homeostatic regulation of intracellular chloride is impaired in pyramidal cells, yet how this dysregulation may lead to seizures has not been explored. Using a realistic neural network model describing ion mechanisms, we show that chloride homeostasis pathology provokes seizure activity analogous to recordings from epileptogenic brain tissue. We show that there is a critical percentage of pathological cells required for seizure initiation. Our model predicts that restoration of the chloride homeostasis in pyramidal cells could be a viable antiepileptic strategy.

Particle swarm optimization algorithm based parameters estimation and control of epileptiform spikes in a neural mass model

This paper proposes an epilepsy detection and closed-loop control strategy based on Particle Swarm Optimization (PSO) algorithm. The proposed strategy can effectively suppress the epileptic spikes in neural mass models, where the epileptiform spikes are recognized as the biomarkers of transitions from the normal (interictal) activity to the seizure (ictal) activity. In addition, the PSO algorithm shows capabilities of accurate estimation for the time evolution of key model parameters and practical detection for all the epileptic spikes. The estimation effects of unmeasurable parameters are improved significantly compared with unscented Kalman filter. When the estimated excitatory-inhibitory ratio exceeds a threshold value, the epileptiform spikes can be inhibited immediately by adopting the proportion-integration controller. Besides, numerical simulations are carried out to illustrate the effectiveness of the proposed method as well as the potential value for the model-based early seizure detection and closed-loop control treatment design.

Suppressing epileptic activity in a neural mass model using a closed-loop proportional-integral controller

Closed-loop control is a promising deep brain stimulation (DBS) strategy that could be used to suppress high-amplitude epileptic activity. However, there are currently no analytical approaches to determine the stimulation parameters for effective and safe treatment protocols. Proportional-integral (PI) control is the most extensively used closed-loop control scheme in the field of control engineering because of its simple implementation and perfect performance. In this study, we took Jansen's neural mass model (NMM) as a test bed to develop a PI-type closed-loop controller for suppressing epileptic activity. A graphical stability analysis method was employed to determine the stabilizing region of the PI controller in the control parameter space, which provided a theoretical guideline for the choice of the PI control parameters. Furthermore, we established the relationship between the parameters of the PI controller and the parameters of the NMM in the form of a stabilizing region, which provided insights into the mechanisms that may suppress epileptic activity in the NMM. The simulation results demonstrated the validity and effectiveness of the proposed closed-loop PI control scheme.

Role of Sodium Channels in Epilepsy

Voltage-gated sodium channels (VGSCs) are fundamentally important for the generation and coordinated transmission of action potentials throughout the nervous system. It is, therefore, unsurprising that they have been shown to play a central role in the genesis and alleviation of epilepsy. Genetic studies on patients with epilepsy have identified more than 700 mutations among the genes that encode for VGSCs attesting to their role in pathogenesis. Further, many common antiepileptic drugs act on VGSCs to suppress seizure activity. Here, we present an account of the role of VGSCs in epilepsy, both through their pathogenic dysfunction and as targets for pharmacotherapy.

Astrocytes as new targets to improve cognitive functions

Astrocytes are now viewed as key elements of brain wiring as well as neuronal communication. Indeed, they not only bridge the gap between metabolic supplies by blood vessels and neurons, but also allow fine control of neurotransmission by providing appropriate signaling molecules and insulation through a tight enwrapping of synapses. Recognition that astroglia is essential to neuronal communication is nevertheless fairly recent and the large body of evidence dissecting such role has focused on the synaptic level by identifying neuro- and gliotransmitters uptaken and released at synaptic or extrasynaptic sites. Yet, more integrated research deciphering the impact of astroglial functions on neuronal network activity have led to the reasonable assumption that the role of astrocytes in supervising synaptic activity translates in influencing neuronal processing and cognitive functions. Several investigations using recent genetic tools now support this notion by showing that inactivating or boosting astroglial function directly affects cognitive abilities. Accordingly, brain diseases resulting in impaired cognitive functions have seen their physiopathological mechanisms revisited in light of this primary protagonist of brain processing. We here provide a review of the current knowledge on the role of astrocytes in cognition and in several brain diseases including neurodegenerative disorders, psychiatric illnesses, as well as other conditions such as epilepsy. Potential astroglial therapeutic targets are also discussed.

A Phenomenological Synapse Model for Asynchronous Neurotransmitter Release

Neurons communicate with each other via synapses. Action potentials cause release of neurotransmitters at the axon terminal. Typically, this neurotransmitter release is tightly time-locked to the arrival of an action potential and is thus called synchronous release. However, neurotransmitter release is stochastic and the rate of release of small quanta of neurotransmitters can be considerably elevated even long after the ceasing of spiking activity, leading to asynchronous release of neurotransmitters. Such asynchronous release varies for tissue and neuron types and has been shown recently to be pronounced in fast-spiking neurons. Notably, it was found that asynchronous release is enhanced in human epileptic tissue implicating a possibly important role in generating abnormal neural activity. Current neural network models for simulating and studying neural activity virtually only consider synchronous release and ignore asynchronous transmitter release. Here, we develop a phenomenological model for asynchronous neurotransmitter release, which, on one hand, captures the fundamental features of the asynchronous release process, and, on the other hand, is simple enough to be incorporated in large-size network simulations. Our proposed model is based on the well-known equations for short-term dynamical synaptic interactions and includes an additional stochastic term for modeling asynchronous release. We use experimental data obtained from inhibitory fast-spiking synapses of human epileptic tissue to fit the model parameters, and demonstrate that our model reproduces the characteristics of realistic asynchronous transmitter release.

Critical Roles of the Direct GABAergic Pallido-cortical Pathway in Controlling Absence Seizures

The basal ganglia (BG), serving as an intermediate bridge between the cerebral cortex and thalamus, are believed to play crucial roles in controlling absence seizure activities generated by the pathological corticothalamic system. Inspired by recent experiments, here we systematically investigate the contribution of a novel identified GABAergic pallido-cortical pathway, projecting from the globus pallidus externa (GPe) in the BG to the cerebral cortex, to the control of absence seizures. By computational modelling, we find that both increasing the activation of GPe neurons and enhancing the coupling strength of the inhibitory pallido-cortical pathway can suppress the bilaterally synchronous 2-4 Hz spike and wave discharges (SWDs) during absence seizures. Appropriate tuning of several GPe-related pathways may also trigger the SWD suppression, through modulating the activation level of GPe neurons. Furthermore, we show that the previously discovered bidirectional control of absence seizures due to the competition between other two BG output pathways also exists in our established model. Importantly, such bidirectional control is shaped by the coupling strength of this direct GABAergic pallido-cortical pathway. Our work suggests that the novel identified pallido-cortical pathway has a functional role in controlling absence seizures and the presented results might provide testable hypotheses for future experimental studies.

Computational modeling of neurostimulation in brain diseases

Neurostimulation as a therapeutic tool has been developed and used for a range of different diseases such as Parkinson's disease, epilepsy, and migraine. However, it is not known why the efficacy of the stimulation varies dramatically across patients or why some patients suffer from severe side effects. This is largely due to the lack of mechanistic understanding of neurostimulation. Hence, theoretical computational approaches to address this issue are in demand. This chapter provides a review of mechanistic computational modeling of brain stimulation. In particular, we will focus on brain diseases, where mechanistic models (e.g., neural population models or detailed neuronal models) have been used to bridge the gap between cellular-level processes of affected neural circuits and the symptomatic expression of disease dynamics. We show how such models have been, and can be, used to investigate the effects of neurostimulation in the diseased brain. We argue that these models are crucial for the mechanistic understanding of the effect of stimulation, allowing for a rational design of stimulation protocols. Based on mechanistic models, we argue that the development of closed-loop stimulation is essential in order to avoid inference with healthy ongoing brain activity. Furthermore, patient-specific data, such as neuroanatomic information and connectivity profiles obtainable from neuroimaging, can be readily incorporated to address the clinical issue of variability in efficacy between subjects. We conclude that mechanistic computational models can and should play a key role in the rational design of effective, fully integrated, patient-specific therapeutic brain stimulation.

Anatomy and physiology of the thick-tufted layer 5 pyramidal neuron

The thick-tufted layer 5 (TTL5) pyramidal neuron is one of the most extensively studied neuron types in the mammalian neocortex and has become a benchmark for understanding information processing in excitatory neurons. By virtue of having the widest local axonal and dendritic arborization, the TTL5 neuron encompasses various local neocortical neurons and thereby defines the dimensions of neocortical microcircuitry. The TTL5 neuron integrates input across all neocortical layers and is the principal output pathway funneling information flow to subcortical structures. Several studies over the past decades have investigated the anatomy, physiology, synaptology, and pathophysiology of the TTL5 neuron. This review summarizes key discoveries and identifies potential avenues of research to facilitate an integrated and unifying understanding on the role of a central neuron in the neocortex.

Impaired dendritic inhibition leads to epileptic activity in a computer model of CA3

Temporal lobe epilepsy (TLE) is a common type of epilepsy with hippocampus as the usual site of origin. The CA3 subfield of hippocampus is reported to have a low epileptic threshold and hence initiates the disorder in patients with TLE. This study computationally investigates how impaired dendritic inhibition of pyramidal cells in the vulnerable CA3 subfield leads to generation of epileptic activity. A model of CA3 subfield consisting of 800 pyramidal cells, 200 basket cells (BC) and 200 Oriens-Lacunosum Moleculare (O-LM) interneurons was used. The dendritic inhibition provided by O-LM interneurons is reported to be selectively impaired in some TLEs. A step-wise approach is taken to investigate how alterations in network connectivity lead to generation of epileptic patterns. Initially, dendritic inhibition alone was reduced, followed by an increase in the external inputs received at the distal dendrites of pyramidal cells, and finally additional changes were made at the synapses between all neurons in the network. In the first case, when the dendritic inhibition of pyramidal cells alone was reduced, the local field potential activity changed from a theta-modulated gamma pattern to a prominently gamma frequency pattern. In the second case, in addition to this reduction of dendritic inhibition, with a simultaneous large increase in the external excitatory inputs received by pyramidal cells, the basket cells entered a state of depolarization block, causing the network to generate a typical ictal activity pattern. In the third case, when the dendritic inhibition onto the pyramidal cells was reduced and changes were simultaneously made in synaptic connectivity between all neurons in the network, the basket cells were again observed to enter depolarization block. In the third case, impairment of dendritic inhibition required to generate an ictal activity pattern was lesser than the two previous cases. Moreover, the ictal like activity began earlier in the third case. Hence, our study suggests that greater synaptic plasticity occurring in the whole network due to increase in reception of external excitatory inputs (due to impaired dendritic inhibition) makes the network more susceptible to generation of epileptic activity.

Computational study of synchrony in fields and microclusters of ephaptically coupled neurons

Neuronal hypersynchrony is implicated in epilepsy and other diseases. The low-frequency, spatially averaged electric fields from many thousands of neurons have been shown to promote synchrony. It remains unclear whether highly transient, spatially localized electric fields from single action potentials (ephaptic coupling) significantly affect spike timing of neighboring cells and in consequence, population synchrony. In this study, we simulated the extracellular potentials and the resulting coupling between neurons in the NEURON environment and generalized their connection rules to create an oscillator network model of a sheet of ephaptically coupled neurons. With the use of both models, we explained several aspects of epileptiform behavior not previously modeled by synaptically coupled networks. Importantly, reduction of neuron spacing induced synchronization via single-spike ephaptic coupling, agreeing with seizure suppression seen clinically and in vitro via extracellular volume adjustment. Further reduction of neuron spacing yielded locally synchronized clusters, providing a mechanism for recent in vitro observations of localized neuronal synchrony in the absence of synaptic and gap-junction coupling.

A biologically constrained, mathematical model of cortical wave propagation preceding seizure termination

Epilepsy--the condition of recurrent, unprovoked seizures--manifests in brain voltage activity with characteristic spatiotemporal patterns. These patterns include stereotyped semi-rhythmic activity produced by aggregate neuronal populations, and organized spatiotemporal phenomena, including waves. To assess these spatiotemporal patterns, we develop a mathematical model consistent with the observed neuronal population activity and determine analytically the parameter configurations that support traveling wave solutions. We then utilize high-density local field potential data recorded in vivo from human cortex preceding seizure termination from three patients to constrain the model parameters, and propose basic mechanisms that contribute to the observed traveling waves. We conclude that a relatively simple and abstract mathematical model consisting of localized interactions between excitatory cells with slow adaptation captures the quantitative features of wave propagation observed in the human local field potential preceding seizure termination.

A portion of inhibitory neurons in human temporal lobe epilepsy are functionally upregulated: an endogenous mechanism for seizure termination

MAIN PROBLEM

Epilepsy is one of the more common neurological disorders. The medication is often ineffective to the patients suffering from intractable temporal lobe epilepsy (TLE). As their seizures are usually self-terminated, the elucidation of the mechanism underlying endogenous seizure termination will help to find a new strategy for epilepsy treatment. We aim to examine the role of inhibitory interneurons in endogenous seizure termination in TLE patients.

METHODS

Whole-cell recordings were conducted on inhibitory interneurons in seizure-onset cortices of intractable TLE patients and the temporal lobe cortices of nonseizure individuals. The intrinsic property of the inhibitory interneurons and the strength of their GABAergic synaptic outputs were measured. The quantitative data were introduced into the computer-simulated neuronal networks to figure out a role of these inhibitory units in the seizure termination.

RESULTS

In addition to functional downregulation, a portion of inhibitory interneurons in seizure-onset cortices were upregulated in encoding the spikes and controlling their postsynaptic neurons. A patch-like upregulation of inhibitory neurons in the local network facilitated seizure termination. The upregulations of both inhibitory neurons and their output synapses synergistically shortened seizure duration, attenuated seizure strength, and terminated seizure propagation.

CONCLUSION

Automatic seizure termination is likely due to the fact that a portion of the inhibitory neurons and synapses are upregulated in the seizure-onset cortices. This mechanism may create novel therapeutic strategies to treat intractable epilepsy, such as the simultaneous upregulation of cortical inhibitory neurons and their output synapses.

Activation of extrasynaptic GABA(A) receptors inhibits cyclothiazide-induced epileptiform activity in hippocampal CA1 neurons

Extrasynaptic GABA(A) receptors (GABA(A)Rs)-mediated tonic inhibition is reported to involve in the pathogenesis of epilepsy. In this study, we used cyclothiazide (CTZ)-induced in vitro brain slice seizure model to explore the effect of selective activation of extrasynaptic GABA(A)Rs by 4,5,6,7-tetrahydroisoxazolo[5,4-c] pyridine-3-ol (THIP) on the CTZ-induced epileptiform activity in hippocampal neurons. Perfusion with CTZ dose-dependently induced multiple epileptiform peaks of evoked population spikes (PSs) in CA1 pyramidal neurons, and treatment with THIP (5 μmol/L) significantly reduced the multiple PS peaks induced by CTZ stimulation. Western blot showed that the δ-subunit of the GABA(A)R, an extrasynaptic specific GABA(A)R subunit, was also significantly down-regulated in the cell membrane 2 h after CTZ treatment. Our results suggest that the CTZ-induced epileptiform activity in hippocampal CA1 neurons is suppressed by the activation of extrasynaptic GABA(A)Rs, and further support the hypothesis that tonic inhibition mediated by extrasynaptic GABA(A)Rs plays a prominent role in seizure generation.

Systems biology, complexity, and the impact on antiepileptic drug discovery

The number of available anticonvulsant drugs increased in the period spanning over more than a century, amounting to the current panoply of nearly two dozen so-called antiepileptic drugs (AEDs). However, none of them actually prevents/reduces the post-brain insult development of epilepsy in man, and in no less than a third of patients with epilepsy, the seizures are not drug-controlled. Plausibly, the enduring limitation of AEDs' efficacy derives from the insufficient understanding of epileptic pathology. This review pinpoints the unbalanced reductionism of the analytic approaches that overlook the intrinsic complexity of epilepsy and of the drug resistance in epilepsy as the core conceptual flaw hampering the discovery of truly antiepileptogenic drugs. A rising awareness of the complexity of epileptic pathology is, however, brought about by the emergence of nonreductionist systems biology (SB) that considers the networks of interactions underlying the normal organismic functions and of SB-based systems (network) pharmacology that aims to restore pathological networks. By now, the systems pharmacology approaches of AED discovery are fairly meager, but their forthcoming development is both a necessity and a realistic prospect, explored in this review.

Structural connectivity based whole brain modelling in epilepsy

Epilepsy is a neurological condition characterised by the recurrence of seizures. During seizures multiple brain areas can behave abnormally. Rather than considering each abnormal area in isolation, one can consider them as an interconnected functional 'network'. Recently, there has been a shift in emphasis to consider epilepsy as a disorder involving more widespread functional brain networks than perhaps was previously thought. The basis for these functional networks is proposed to be the static structural brain network established through the connectivity of the white matter. Additionally, it has also been argued that time varying aspects of epilepsy are of crucial importance and as such computational models of these dynamical properties have recently advanced. We describe how dynamic computer models can be combined with static human in vivo connectivity obtained through diffusion weighted magnetic resonance imaging. We predict that in future the use of these two methods in concert will lead to predictions for optimal surgery and brain stimulation sites for epilepsy and other neurological disorders.