Simplest relationship between local field potential and intracellular signals in layered neural tissue

Seeing beyond the spikes: reconstructing the complete spatiotemporal membrane potential distribution from paired intra- and extracellular recordings

Although electrophysiologists have been recording intracellular neural activity routinely ever since the ground-breaking work of Hodgkin and Huxley, and extracellular multichannel electrodes have also been used frequently and extensively, a practical experimental method to track changes in membrane potential along a complete single neuron is still lacking. Instead of obtaining multiple intracellular measurements on the same neuron, we propose an alternative method by combining single-channel somatic patch-clamp and multichannel extracellular potential recordings. In this work, we show that it is possible to reconstruct the complete spatiotemporal distribution of the membrane potential of a single neuron with the spatial resolution of an extracellular probe during action potential generation. Moreover, the reconstruction of the membrane potential allows us to distinguish between the two major but previously hidden components of the current source density (CSD) distribution: the resistive and the capacitive currents. This distinction provides a clue to the clear interpretation of the CSD analysis, because the resistive component corresponds to transmembrane ionic currents (all the synaptic, voltage-sensitive and passive currents), whereas capacitive currents are considered to be the main contributors of counter-currents. We validate our model-based reconstruction approach on simulations and demonstrate its application to experimental data obtained in vitro via paired extracellular and intracellular recordings from a single pyramidal cell of the rat hippocampus. In perspective, the estimation of the spatial distribution of resistive membrane currents makes it possible to distiguish between active and passive sinks and sources of the CSD map and the localization of the synaptic input currents, which make the neuron fire. KEY POINTS: A new computational method is introduced to calculate the unbiased current source density distribution on a single neuron with known morphology. The relationship between extracellular and intracellular electric potential is determined via mathematical formalism, and a new reconstruction method is applied to reveal the full spatiotemporal distribution of the membrane potential and the resistive and capacitive current components. The new reconstruction method was validated on simulations. Simultaneous and colocalized whole-cell patch-clamp and multichannel silicon probe recordings were performed from the same pyramidal neuron in the rat hippocampal CA1 region, in vitro. The method was applied in experimental measurements and returned precise and distinctive characteristics of various intracellular phenomena, such as action potential generation, signal back-propagation and the initial dendritic depolarization preceding the somatic action potential.

HH model based smart deep brain stimulator to detect, predict and control epilepsy using machine learning algorithm

BACKGROUND

Epilepsy is the most common neurological disorder in the world. To control epilepsy, deep brain stimulation is one of the widely accepted treatment techniques. However, conventional deep brain stimulation technique provides continuous stimulation without optimizing the stimulation parameters, resulting in adverse side effects and unexpected death. Hence, understanding the dynamic behavior of brain neural networks at a cellular level is required for patient-specific epilepsy treatment. Considering the underlying mechanism of a single neuronal shift in the brain neural network, computational model-based techniques have a new face for healthcare, which aims to develop effective medical devices for preclinical investigations.

NEW METHOD

This paper discusses the design of a Smart Deep Brain Stimulator (SDBS) using the Hodgkin-Huxley (HH) conductance-based cellular model of brain neurons to automatically detect, predict and regulate epilepsy against patient-specific conditions. Epileptic activity is simulated as a spike train of action potential due to sodium and potassium channel conductance variations in the single-neuron HH model. The proposed SDBS consists of three components:- i) seizure detection using bagging and boosting-based ensemble machine learning classifiers, ii) channel conductance prediction using Long Short Term Memory-Recurrent Neural Network (LSTM-RNN) based Deep Neural Network (DNN) for updating model parameters of brain neuron, and iii) model-based intelligent control of epileptic seizure with Nonlinear Autoregressive Moving Average-L2 (NARMA-L2) Controller and Nonlinear Model Predictive Controller (NMPC).

RESULTS

For effective treatment, improving the overall accuracy and efficiency of SDBS is essential. For epilepsy detection, the ensemble bagging machine learning algorithm provides better accuracy of 92.7% compared to the ensemble boosting algorithm. LSTM-RNN deep neural network model with four layers predicts the variations in channel conductance with Root Mean Square Error (RMSE) of 0.00568 and 0.009081 for sodium and potassium channel conductance, respectively. From the closed-loop performances of SDBS with an intelligent control scheme, it is observed that SDBS with NMPC provides efficient and accurate stimulation with minimum energy consumption. From a stability point of view, SDBS with NMPC provides better stability than SDBS with NARMA-L2 Controller.

COMPARISON WITH EXISTING METHOD

The proposed SDBS is designed to generate accurate stimulation pulses for epilepsy patients with specific conditions depending on the neuronal activity of a single neuron. Moreover, it will also adapt to the dynamic condition of epilepsy patients. The existing deep brain stimulator continuously provides stimulation pulses without adapting to the patient's conditions.

CONCLUSION

The proposed SDBS could provide patient-specific treatment based on sodium/potassium channel conductance variations of brain neurons. It will help increase the use of deep brain stimulation techniques and reduce sudden death. Furthermore, the proposed technique will be extended to neural network models with larger neuronal populations to improve the practical feasibility.

A strategy for mapping biophysical to abstract neuronal network models applied to primary visual cortex

A fundamental challenge for the theoretical study of neuronal networks is to make the link between complex biophysical models based directly on experimental data, to progressively simpler mathematical models that allow the derivation of general operating principles. We present a strategy that successively maps a relatively detailed biophysical population model, comprising conductance-based Hodgkin-Huxley type neuron models with connectivity rules derived from anatomical data, to various representations with fewer parameters, finishing with a firing rate network model that permits analysis. We apply this methodology to primary visual cortex of higher mammals, focusing on the functional property of stimulus orientation selectivity of receptive fields of individual neurons. The mapping produces compact expressions for the parameters of the abstract model that clearly identify the impact of specific electrophysiological and anatomical parameters on the analytical results, in particular as manifested by specific functional signatures of visual cortex, including input-output sharpening, conductance invariance, virtual rotation and the tilt after effect. Importantly, qualitative differences between model behaviours point out consequences of various simplifications. The strategy may be applied to other neuronal systems with appropriate modifications.

A hierarchy of theoretical approaches to study a neuronal network depends on a tradeoff between biological fidelity and mathematical tractibility. Biophysically-detailed models consider cellular mechanisms and anatomically defined synaptic circuits, but are often too complex to reveal insights into fundamental principles. In contrast, increasingly abstract reduced models facilitate analytical insights. To better ground the latter to the underlying biology, we describe a systematic procedure to move across the model hierarchy that allows understanding how changes in biological parameters—physiological, pathophysiological, or because of new data—impact the behaviour of the network. We apply this approach to mammalian primary visual cortex, and examine how the different models in the hierarchy reproduce functional signatures of this area, in particular the tuning of neurons to the orientation of a visual stimulus. Our work provides a navigation of the complex parameter space of neural network models faithful to biology, as well as highlighting how simplifications made for mathematical convenience can fundamentally change their behaviour.

Refractory density model of cortical direction selectivity: Lagged-nonlagged, transient-sustained, and On-Off thalamic neuron-based mechanisms and intracortical amplification

A biophysically detailed description of the mechanisms of the primary vision is still being developed. We have incorporated a simplified, filter-based description of retino-thalamic visual signal processing into the detailed, conductance-based refractory density description of the neuronal population activity of the primary visual cortex. We compared four mechanisms of the direction selectivity (DS), three of them being based on asymmetrical projections of different types of thalamic neurons to the cortex, distinguishing between (i) lagged and nonlagged, (ii) transient and sustained, and (iii) On and Off neurons. The fourth mechanism implies a lack of subcortical bias and is an epiphenomenon of intracortical interactions between orientation columns. The simulations of the cortical response to moving gratings have verified that first three mechanisms provide DS to an extent compared with experimental data and that the biophysical model realistically reproduces characteristics of the visual cortex activity, such as membrane potential, firing rate, and synaptic conductances. The proposed model reveals the difference between the mechanisms of both the intact and the silenced cortex, favoring the second mechanism. In the fourth case, DS is weaker but significant; it completely vanishes in the silenced cortex.DS in the On-Off mechanism derives from the nonlinear interactions within the orientation map. Results of simulations can help to identify a prevailing mechanism of DS in V1. This is a step towards a comprehensive biophysical modeling of the primary visual system in the frameworks of the population rate coding concept.

A major mechanism that underlies tuning of cortical neurons to the direction of a moving stimulus is still debated. Considering the visual cortex structured with orientation-selective columns, we have realized and compared in our biophysically detailed mathematical model four hypothetical mechanisms of the direction selectivity (DS) known from experiments. The present model accomplishes our previous model that was tuned to experimental data on excitability in slices and reproduces orientation tuning effects in vivo. In simulations, we have found that the convergence of inputs from so-called transient and sustained (or lagged and nonlagged) thalamic neurons in the cortex provides an initial bias for DS, whereas cortical interactions amplify the tuning. In the absence of any bias, DS emerges as an epiphenomenon of the orientation map. In the case of a biased convergence of On- and Off- thalamic inputs, DS emerges with the help of the intracortical interactions on the orientation map, also. Thus, we have proposed a comprehensive description of the primary vision and revealed characteristic features of different mechanisms of DS in the visual cortex with columnar structure.

Generation of Sharp Wave-Ripple Events by Disinhibition

Sharp wave-ripple complexes (SWRs) are hippocampal network phenomena involved in memory consolidation. To date, the mechanisms underlying their occurrence remain obscure. Here, we show how the interactions between pyramidal cells, parvalbumin-positive (PV⁺) basket cells, and an unidentified class of anti-SWR interneurons can contribute to the initiation and termination of SWRs. Using a biophysically constrained model of a network of spiking neurons and a rate-model approximation, we demonstrate that SWRs emerge as a result of the competition between two interneuron populations and the resulting disinhibition of pyramidal cells. Our models explain how the activation of pyramidal cells or PV⁺ cells can trigger SWRs, as shown in vitro, and suggests that PV⁺ cell-mediated short-term synaptic depression influences the experimentally reported dynamics of SWR events. Furthermore, we predict that the silencing of anti-SWR interneurons can trigger SWRs. These results broaden our understanding of the microcircuits supporting the generation of memory-related network dynamics.

SIGNIFICANCE STATEMENT The hippocampus is a part of the mammalian brain that is crucial for episodic memories. During periods of sleep and inactive waking, the extracellular activity of the hippocampus is dominated by sharp wave-ripple events (SWRs), which have been shown to be important for memory consolidation. The mechanisms regulating the emergence of these events are still unclear. We developed a computational model to study the emergence of SWRs and to explain the roles of different cell types in regulating them. The model accounts for several previously unexplained features of SWRs and thus advances the understanding of memory-related dynamics.

KCC2 antagonism and gabaergic synchronization in the entorhinal cortex in the absence of ionotropic glutamatergic receptor signalling

γ-Aminobutyric acid (GABA), which is released by interneurons, plays an active role in generating interictal epileptiform spikes during blockade of ionotropic glutamatergic signalling, but it remains unclear whether and how the K+-Cl- cotransporter 2 (KCC2) influences these paroxysmal events. Therefore, we employed tetrode recordings in the in vitro rat entorhinal cortex (EC) to analyze the effects of the KCC2 antagonist VU0463271 on 4-aminopyridine (4AP)-induced interictal spikes that were pharmacologically isolated by applying ionotropic glutamatergic receptor antagonists. After the addition of VU0463271, these interictal spikes continued to occur at similar rates as in control (i.e., during application of 4AP with ionotropic glutamatergic receptor antagonists) but were smaller and shorter. Despite the absence of ionotropic glutamatergic receptor signalling, both interneurons and principal cells increased their firing during interictal spikes. Moreover, we found that KCC2 antagonism increased interneuron firing but decreased principal cell firing during the interictal spike rising phase; in contrast, during the falling phase, interneuron firing decreased in the presence of VU0463271 while no change was observed in principal cell firing. Overall, our results show that KCC2 antagonism enhances interneuron excitability at the onset of interictal spikes generated by the EC neuronal networks during blockade of ionotropic glutamatergic transmission but disrupts later neuronal recruitment.

Mind the last spike - firing rate models for mesoscopic populations of spiking neurons

The dominant modeling framework for understanding cortical computations are heuristic firing rate models. Despite their success, these models fall short to capture spike synchronization effects, to link to biophysical parameters and to describe finite-size fluctuations. In this opinion article, we propose that the refractory density method (RDM), also known as age-structured population dynamics or quasi-renewal theory, yields a powerful theoretical framework to build rate-based models for mesoscopic neural populations from realistic neuron dynamics at the microscopic level. We review recent advances achieved by the RDM to obtain efficient population density equations for networks of generalized integrate-and-fire (GIF) neurons - a class of neuron models that has been successfully fitted to various cell types. The theory not only predicts the nonstationary dynamics of large populations of neurons but also permits an extension to finite-size populations and a systematic reduction to low-dimensional rate dynamics. The new types of rate models will allow a re-examination of models of cortical computations under biological constraints.

Mathematical model of Na-K-Cl homeostasis in ictal and interictal discharges

Despite big experimental data on the phenomena and mechanisms of the generation of ictal and interictal discharges (IDs and IIDs), mathematical models that can describe the synaptic interactions of neurons and the ionic dynamics in biophysical detail are not well-established. Based on experimental recordings of combined hippocampal-entorhinal cortex slices from rats in a high-potassium and a low-magnesium solution containing 4-aminopyridine as well as previous observations of similar experimental models, this type of mathematical model has been developed. The model describes neuronal excitation through the application of the conductance-based refractory density approach for three neuronal populations: two populations of glutamatergic neurons with hyperpolarizing and depolarizing GABAergic synapses and one GABAergic population. The ionic dynamics account for the contributions of voltage-gated and synaptic channels, active and passive transporters, and diffusion. The relatively slow dynamics of potassium, chloride, and sodium ion concentrations determine the transitions from pure GABAergic IIDs to IDs and GABA-glutamatergic IIDs. The model reproduces different types of IIDs, including those initiated by interneurons; repetitive IDs; tonic and bursting modes of an ID composed of clustered IID-like events. The simulations revealed contributions from different ionic channels to the ion concentration dynamics before and during ID generation. The proposed model is a step forward to an optimal mathematical description of the mechanisms of epileptic discharges.

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.

Computational model of interictal discharges triggered by interneurons

Interictal discharges (IIDs) are abnormal waveforms registered in the periods before or between seizures. IIDs that are initiated by GABAergic interneurons have not been mathematically modeled yet. In the present study, a mathematical model that describes the mechanisms of these discharges is proposed. The model is based on the experimental recordings of IIDs in pyramidal neurons of the rat entorhinal cortex and estimations of synaptic conductances during IIDs. IIDs were induced in cortico-hippocampal slices by applying an extracellular solution with 4-aminopyridine, high potassium, and low magnesium concentrations. Two different types of IIDs initiated by interneurons were observed. The first type of IID (IID1) was pure GABAergic. The second type of IID (IID2) was induced by GABAergic excitation and maintained by recurrent interactions of both GABA- and glutamatergic neuronal populations. The model employed the conductance-based refractory density (CBRD) approach, which accurately approximates the firing rate of a population of similar Hodgkin-Huxley-like neurons. The model of coupled excitatory and inhibitory populations includes AMPA, NMDA, and GABA-receptor-mediated synapses and gap junctions. These neurons receive both arbitrary deterministic input and individual colored Gaussian noise. Both types of IIDs were successfully reproduced in the model by setting two different depolarized levels for GABA-mediated current reversal potential. It was revealed that short-term synaptic depression is a crucial factor in ceasing each of the discharges, and it also determines their durations and frequencies.

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.

Detection of transient synchrony across oscillating receptors by the central electrosensory system of mormyrid fish

Recently, we reported evidence for a novel mechanism of peripheral sensory coding based on oscillatory synchrony. Spontaneously oscillating electroreceptors in weakly electric fish (Mormyridae) respond to electrosensory stimuli with a phase reset that results in transient synchrony across the receptor population (Baker et al., 2015). Here, we asked whether the central electrosensory system actually detects the occurrence of synchronous oscillations among receptors. We found that electrosensory stimulation elicited evoked potentials in the midbrain exterolateral nucleus at a short latency following receptor synchronization. Frequency tuning in the midbrain resembled peripheral frequency tuning, which matches the intrinsic oscillation frequencies of the receptors. These frequencies are lower than those in individual conspecific signals, and instead match those found in collective signals produced by groups of conspecifics. Our results provide further support for a novel mechanism for sensory coding based on the detection of oscillatory synchrony among peripheral receptors.

DOI:http://dx.doi.org/10.7554/eLife.16851.001