Electrogenic properties of the Na⁺/K⁺ ATPase control transitions between normal and pathological brain states

Multistability of bursting rhythms in a half-center oscillator and the protective effects of synaptic inhibition

For animals to meet environmental challenges, the activity patterns of specialized oscillatory neural circuits, central pattern generators (CPGs), controlling rhythmic movements like breathing and locomotion, are adjusted by neuromodulation. As a representative example, the leech heartbeat is controlled by a CPG driven by two pairs of mutually inhibitory interneurons, heart interneuron (HN) half-center oscillators (HCO). Experiments and modeling indicate that neuromodulation of HCO navigates this CPG between dysfunctional regimes by employing a co-regulating inverted relation; reducing Na+/K+ pump current and increasing hyperpolarization-activated (h-) current. Simply reducing pump activity or increasing h-current leads to either seizure-like bursting or an asymmetric bursting dysfunctional regime, respectively. Here, we demonstrate through modeling that, alongside this coregulation path, a new bursting regime emerges. Both regimes fulfill the criteria for functional bursting activity. Although the cycle periods and burst durations of these patterns are roughly the same, the new one exhibits an intra-burst spike frequency that is twice as high as the other. This finding suggests that neuromodulation could introduce additional functional regimes with higher spike frequency, and thus more effective synaptic transmission to motor neurons. We found that this new regime co-exists with the original bursting. The HCO can be switched between them by a short pulse of excitatory or inhibitory conductance. In this domain of coexisting functional patterns, an isolated cell model exhibits only one regime, a severely dysfunctional plateau-containing, seizure-like activity. This aligns with widely reported notion that deficiency of inhibition can cause seizures and other dysfunctional neural activities. We show that along the coregulation path of neuromodulation, the high excitability of the single HNs induced by myomodulin is harnessed by mutually inhibitory synaptic interactions of the HCO into the functional bursting pattern.

Cholinergic modulation supports dynamic switching of resting state networks through selective DMN suppression

Brain activity during the resting state is widely used to examine brain organization, cognition and alterations in disease states. While it is known that neuromodulation and the state of alertness impact resting-state activity, neural mechanisms behind such modulation of resting-state activity are unknown. In this work, we used a computational model to demonstrate that change in excitability and recurrent connections, due to cholinergic modulation, impacts resting-state activity. The results of such modulation in the model match closely with experimental work on direct cholinergic modulation of Default Mode Network (DMN) in rodents. We further extended our study to the human connectome derived from diffusion-weighted MRI. In human resting-state simulations, an increase in cholinergic input resulted in a brain-wide reduction of functional connectivity. Furthermore, selective cholinergic modulation of DMN closely captured experimentally observed transitions between the baseline resting state and states with suppressed DMN fluctuations associated with attention to external tasks. Our study thus provides insight into potential neural mechanisms for the effects of cholinergic neuromodulation on resting-state activity and its dynamics.

Uncertainty quantification and sensitivity analysis of neuron models with ion concentration dynamics

This paper provides a comprehensive and computationally efficient case study for uncertainty quantification (UQ) and global sensitivity analysis (GSA) in a neuron model incorporating ion concentration dynamics. We address how challenges with UQ and GSA in this context can be approached and solved, including challenges related to computational cost, parameters affecting the system's resting state, and the presence of both fast and slow dynamics. Specifically, we analyze the electrodiffusive neuron-extracellular-glia (edNEG) model, which captures electrical potentials, ion concentrations (Na+, K+, Ca2+, and Cl-), and volume changes across six compartments. Our methodology includes a UQ procedure assessing the model's reliability and susceptibility to input uncertainty and a variance-based GSA identifying the most influential input parameters. To mitigate computational costs, we employ surrogate modeling techniques, optimized using efficient numerical integration methods. We propose a strategy for isolating parameters affecting the resting state and analyze the edNEG model dynamics under both physiological and pathological conditions. The influence of uncertain parameters on model outputs, particularly during spiking dynamics, is systematically explored. Rapid dynamics of membrane potentials necessitate a focus on informative spiking features, while slower variations in ion concentrations allow a meaningful study at each time point. Our study offers valuable guidelines for future UQ and GSA investigations on neuron models with ion concentration dynamics, contributing to the broader application of such models in computational neuroscience.

The Reciprocal Relationship between Sleep and Epilepsy

The relationship between sleep and epilepsy is bidirectional. Certain epilepsy syndromes predominantly or exclusively manifest during sleep, with seizures frequently originating from non-rapid eye movement (NREM) sleep. Interictal epileptiform discharges observed on electroencephalograms are most likely to be activated during the deep NREM sleep stage known as N3. Conversely, epileptiform discharges, anti-seizure medications (ASMs), as well as other anti-seizure therapies can exert detrimental effects on sleep architecture. Moreover, the co-occurrence of sleep disorders has the potential to exacerbate seizure control. Understating the relationship between sleep and epilepsy is crucial for healthcare providers. Addressing and managing sleep-related problems in individuals with epilepsy can potentially contribute to improved seizure control and overall well-being. At the same time, improving seizure control can improve sleep quality and quantity, thus further improving the health of individuals with epilepsy.

Bursting Dynamics Based on the Persistent Na+ and Na+/K+ Pump Currents: A Dynamic Clamp Approach

Life-supporting rhythmic motor functions like heart-beating in invertebrates and breathing in vertebrates require an indefatigable generation of a robust rhythm by specialized oscillatory circuits, central pattern generators (CPGs). These CPGs should be sufficiently flexible to adjust to environmental changes and behavioral goals. Continuous self-sustained operation of bursting neurons requires intracellular Na+ concentration to remain in a functional range and to have checks and balances of the Na+ fluxes met on a cycle-to-cycle basis during bursting. We hypothesize that at a high excitability state, the interaction of the Na+/K+ pump current, Ipump, and persistent Na+ current, INaP, produces a mechanism supporting functional bursting. INaP is a low voltage-activated inward current that initiates and supports the bursting phase. This current does not inactivate and is a significant source of Na+ influx. Ipump is an outward current activated by [Na+]i and is the major source of Na+ efflux. Both currents are active and counteract each other between and during bursts. We apply a combination of electrophysiology, computational modeling, and dynamic clamp to investigate the role of Ipump and INaP in the leech heartbeat CPG interneurons (HN neurons). Applying dynamic clamp to introduce additional Ipump and INaP into the dynamics of living synaptically isolated HN neurons in real time, we show that their joint increase produces transition into a new bursting regime characterized by higher spike frequency and larger amplitude of the membrane potential oscillations. Further increase of Ipump speeds up this rhythm by shortening burst duration (BD) and interburst interval (IBI).

From phenomenological to biophysical models of seizures

Epilepsy is a complex disease that requires various approaches for its study. In this short review, we discuss the contribution of theoretical and computational models. The review presents theoretical frameworks that underlie the understanding of certain seizure properties and their classification based on their dynamical properties at the onset and offset of seizures. Dynamical system tools are valuable resources in the study of seizures. By analyzing the complex, dynamic behavior of seizures, these tools can provide insights into seizure mechanisms and offer a framework for their classification. Additionally, computational models have high potential for clinical applications, as they can be used to develop more accurate diagnostic and personalized medicine tools. We discuss various modeling approaches that span different scales and levels, while also questioning the neurocentric view, and emphasize the importance of considering glial cells. Finally, we explore the epistemic value provided by this type of approach.

Non-trivial dynamics in a model of glial membrane voltage driven by open potassium pores

Despite the molecular evidence that a nearly linear steady-state current-voltage relationship in mammalian astrocytes reflects a total current resulting from more than one differentially regulated K+ conductance, detailed ordinary differential equation (ODE) models of membrane voltage Vm are still lacking. Various experimental results reporting altered rectification of the major Kir currents in glia, dominated by Kir4.1, have motivated us to develop a detailed model of Vm dynamics incorporating the weaker potassium K2P-TREK1 current in addition to Kir4.1, and study the stability of the resting state Vr. The main question is whether, with the loss of monotonicity in glial I-V curve resulting from altered Kir rectification, the nominal resting state Vr remains stable, and the cell retains the trivial, potassium electrode behavior with Vm after EK. The minimal two-dimensional model of Vm near Vr showed that an N-shape deformed Kir I-V curve induces multistability of Vm in a model that incorporates K2P activation kinetics, and nonspecific K+ leak currents. More specifically, an asymmetrical, nonlinear decrease of outward Kir4.1 conductance, turning the channels into inward rectifiers, introduces instability of Vr. That happens through a robust bifurcation giving birth to a second, more depolarized stable resting state Vdr > -10 mV. Realistic recordings from electrographic seizures were used to perturb the model. Simulations of the model perturbed by constant current through gap junctions and seizure-like discharges as local field potentials led to depolarization and switching of Vm between the two stable states, in a downstate-upstate manner. In the event of prolonged depolarizations near Vdr, such catastrophic instability would affect all aspects of the glial function, from metabolic support to membrane transport, and practically all neuromodulatory roles assigned to glia.

Optogenetic Low-Frequency Stimulation of Principal Neurons, but Not Parvalbumin-Positive Interneurons, Prevents Generation of Ictal Discharges in Rodent Entorhinal Cortex in an In Vitro 4-Aminopyridine Model

Low-frequency electrical stimulation is used to treat some drug-resistant forms of epilepsy. Despite the effectiveness of the method in suppressing seizures, there is a considerable risk of side effects. An optogenetic approach allows the targeting of specific populations of neurons, which can increase the effectiveness and safety of low-frequency stimulation. In our study, we tested the efficacy of the suppression of ictal activity in entorhinal cortex slices in a 4-aminopyridine model with three variants of low-frequency light stimulation (LFLS): (1) activation of excitatory and inhibitory neurons (on Thy1-ChR2-YFP mice), (2) activation of inhibitory interneurons only (on PV-Cre mice after virus injection with channelrhodopsin2 gene), and (3) hyperpolarization of excitatory neurons (on Wistar rats after virus injection with archaerhodopsin gene). Only in the first variant did simultaneous LFLS of excitatory and inhibitory neurons replace ictal activity with interictal activity. We suggest that LFLS caused changes in the concentration gradients of K⁺ and Na⁺ cations across the neuron membrane, which activated Na-K pumping. According to the mathematical modeling, the increase in Na-K pump activity in neurons induced by LFLS led to an antiepileptic effect. Thus, a less specific and generalized optogenetic effect on entorhinal cortex neurons was more effective in suppressing ictal activity in the 4-aminopyridine model.

Relationship between chronic hypoxia and seizure susceptibility

Chronic hypobaric hypoxia in high‐altitude areas is closely related to the occurrence of many neurological diseases. Among these diseases, epilepsy is a common disease of the nervous system that is difficult to diagnose and treat, with a long treatment cycle. As of 2019, there were more than 70 million epilepsy patients worldwide, including 10 million in China. Studies have shown that chronic hypoxia promotes the occurrence and development of epilepsy, and elucidation of the relationship between chronic hypoxia and epilepsy is important for studying the pathogenesis of epilepsy and exploring the potential characteristics of epilepsy and new drug targets for epilepsy. In this article, we review the factors that may cause increased seizure susceptibility in chronic hypoxia and consider the potential relationship between chronic hypobaric hypoxia and seizure susceptibility in high‐altitude areas and prospects surrounding related research in the future.

Focal seizures are organized by feedback between neural activity and ion concentration changes

Human and animal EEG data demonstrate that focal seizures start with low-voltage fast activity, evolve into rhythmic burst discharges and are followed by a period of suppressed background activity. This suggests that processes with dynamics in the range of tens of seconds govern focal seizure evolution. We investigate the processes associated with seizure dynamics by complementing the Hodgkin-Huxley mathematical model with the physical laws that dictate ion movement and maintain ionic gradients. Our biophysically realistic computational model closely replicates the electrographic pattern of a typical human focal seizure characterized by low voltage fast activity onset, tonic phase, clonic phase and postictal suppression. Our study demonstrates, for the first time in silico, the potential mechanism of seizure initiation by inhibitory interneurons via the initial build-up of extracellular K+ due to intense interneuronal spiking. The model also identifies ionic mechanisms that may underlie a key feature in seizure dynamics, i.e., progressive slowing down of ictal discharges towards the end of seizure. Our model prediction of specific scaling of inter-burst intervals is confirmed by seizure data recorded in the whole guinea pig brain in vitro and in humans, suggesting that the observed termination pattern may hold across different species. Our results emphasize ionic dynamics as elementary processes behind seizure generation and indicate targets for new therapeutic strategies.

Start Me Up: How Can Surrounding Gangliosides Affect Sodium-Potassium ATPase Activity and Steer towards Pathological Ion Imbalance in Neurons?

Gangliosides, amphiphilic glycosphingolipids, tend to associate laterally with other membrane constituents and undergo extensive interactions with membrane proteins in cis or trans configurations. Studies of human diseases resulting from mutations in the ganglioside biosynthesis pathway and research on transgenic mice with the same mutations implicate gangliosides in the pathogenesis of epilepsy. Gangliosides are reported to affect the activity of the Na⁺/K⁺-ATPase, the ubiquitously expressed plasma membrane pump responsible for the stabilization of the resting membrane potential by hyperpolarization, firing up the action potential and ion homeostasis. Impaired Na⁺/K⁺-ATPase activity has also been hypothesized to cause seizures by several mechanisms. In this review we present different epileptic phenotypes that are caused by impaired activity of Na⁺/K⁺-ATPase or changed membrane ganglioside composition. We further discuss how gangliosides may influence Na⁺/K⁺-ATPase activity by acting as lipid sorting machinery providing the optimal stage for Na⁺/K⁺-ATPase function. By establishing a distinct lipid environment, together with other membrane lipids, gangliosides possibly modulate Na⁺/K⁺-ATPase activity and aid in “starting up” and “turning off” this vital pump. Therefore, structural changes of neuronal membranes caused by altered ganglioside composition can be a contributing factor leading to aberrant Na⁺/K⁺-ATPase activity and ion imbalance priming neurons for pathological firing.

Ictal wavefront propagation in slices and simulations with conductance-based refractory density model

The mechanisms determining ictal discharge (ID) propagation are still not clear. In the present study, we aimed to examine these mechanisms in animal and mathematical models of epileptiform activity. Using double-patch and extracellular potassium ion concentration recordings in rat hippocampal-cortical slices, we observed that IDs moved at a speed of about 1 mm/s or less. The mechanisms of such slow propagation have been studied with a mathematical, conductance-based refractory density (CBRD) model that describes the GABA- and glutamatergic neuronal populations’ interactions and ion dynamics in brain tissue. The modeling study reveals two main factors triggerring IDs: (i) increased interneuronal activity leading to chloride ion accumulation and a consequent depolarizing GABAergic effect and (ii) the elevation of extracellular potassium ion concentration. The local synaptic transmission followed by local potassium ion extrusion and GABA receptor-mediated chloride ion accumulation underlies the ID wavefront’s propagation. In contrast, potassium ion diffusion in the extracellular space is slower and does not affect ID’s speed. The short discharges, constituting the ID, propagate much faster than the ID front. The accumulation of sodium ions inside neurons due to their hyperactivity and glutamatergic currents boosts the Na⁺/K⁺ pump, which terminates the ID. Knowledge of the mechanism of ID generation and propagation contributes to the development of new treatments against epilepsy.

During an epileptic seizure, neuronal excitation spreads across the brain tissue and is accompanied by significant changes in ionic concentrations. Ictal discharge front spreads at low speeds, less than 1 mm/s. Mechanisms underlying this phenomenon are not yet well understood. We study these mechanisms using electrophysiological recordings in brain slices and computer simulations. Our detailed biophysical model describing neuronal populations’ interaction, spatial propagation, and ionic dynamics reproduces the generation and propagation of spontaneously repeating ictal discharges. The simulations are consistent with our recordings of the electrical activity and the extracellular potassium ion concentration. We distinguished between the two alternative mechanisms of the ictal wavefront propagation: (i) the diffusion of potassium ions released from excited neurons, which depolarizes distant neurons and thus supports excitation, and (ii) the axonal spread of excitation followed by the local extracellular potassium ion accumulation that supports the excitation. Our simulations provide evidence in favor of the latter mechanism. Our experiment-based modeling contributes to a mathematical description of brain tissue functioning and potentially contributes to developing new treatments against epilepsy.

A unified physiological framework of transitions between seizures, sustained ictal activity and depolarization block at the single neuron level

The majority of seizures recorded in humans and experimental animal models can be described by a generic phenomenological mathematical model, the Epileptor. In this model, seizure-like events (SLEs) are driven by a slow variable and occur via saddle node (SN) and homoclinic bifurcations at seizure onset and offset, respectively. Here we investigated SLEs at the single cell level using a biophysically relevant neuron model including a slow/fast system of four equations. The two equations for the slow subsystem describe ion concentration variations and the two equations of the fast subsystem delineate the electrophysiological activities of the neuron. Using extracellular K⁺ as a slow variable, we report that SLEs with SN/homoclinic bifurcations can readily occur at the single cell level when extracellular K⁺ reaches a critical value. In patients and experimental models, seizures can also evolve into sustained ictal activity (SIA) and depolarization block (DB), activities which are also parts of the dynamic repertoire of the Epileptor. Increasing extracellular concentration of K⁺ in the model to values found during experimental status epilepticus and DB, we show that SIA and DB can also occur at the single cell level. Thus, seizures, SIA, and DB, which have been first identified as network events, can exist in a unified framework of a biophysical model at the single neuron level and exhibit similar dynamics as observed in the Epileptor.Author Summary: Epilepsy is a neurological disorder characterized by the occurrence of seizures. Seizures have been characterized in patients in experimental models at both macroscopic and microscopic scales using electrophysiological recordings. Experimental works allowed the establishment of a detailed taxonomy of seizures, which can be described by mathematical models. We can distinguish two main types of models. Phenomenological (generic) models have few parameters and variables and permit detailed dynamical studies often capturing a majority of activities observed in experimental conditions. But they also have abstract parameters, making biological interpretation difficult. Biophysical models, on the other hand, use a large number of variables and parameters due to the complexity of the biological systems they represent. Because of the multiplicity of solutions, it is difficult to extract general dynamical rules. In the present work, we integrate both approaches and reduce a detailed biophysical model to sufficiently low-dimensional equations, and thus maintaining the advantages of a generic model. We propose, at the single cell level, a unified framework of different pathological activities that are seizures, depolarization block, and sustained ictal activity.

A missense mutation converts the Na+,K+-ATPase into an ion channel and causes therapy-resistant epilepsy

The ion pump Na⁺,K⁺-ATPase is a critical determinant of neuronal excitability; however, its role in the etiology of diseases of the central nervous system (CNS) is largely unknown. We describe here the molecular phenotype of a Trp931Arg mutation of the Na⁺,K⁺-ATPase catalytic α1 subunit in an infant diagnosed with therapy-resistant lethal epilepsy. In addition to the pathological CNS phenotype, we also detected renal wasting of Mg²⁺. We found that membrane expression of the mutant α1 protein was low, and ion pumping activity was lost. Arginine insertion into membrane proteins can generate water-filled pores in the plasma membrane, and our molecular dynamic (MD) simulations of the principle states of Na⁺,K⁺-ATPase transport demonstrated massive water inflow into mutant α1 and destabilization of the ion-binding sites. MD simulations also indicated that a water pathway was created between the mutant arginine residue and the cytoplasm, and analysis of oocytes expressing mutant α1 detected a nonspecific cation current. Finally, neurons expressing mutant α1 were observed to be depolarized compared with neurons expressing wild-type protein, compatible with a lowered threshold for epileptic seizures. The results imply that Na⁺,K⁺-ATPase should be considered a neuronal locus minoris resistentia in diseases associated with epilepsy and with loss of plasma membrane integrity.

Interactions among diameter, myelination, and the Na/K pump affect axonal resilience to high-frequency spiking

The reliability of spike propagation in axons is determined by complex interactions among ionic currents, ion pumps, and morphological properties. We use compartment-based modeling to reveal that interactions of diameter, myelination, and the Na/K pump determine the reliability of high-frequency spike propagation. By acting as a “reservoir” of nodal Na⁺ influx, myelinated compartments efficiently increase propagation reliability. Although spike broadening was thought to oppose fast spiking, its effect on spike propagation is complicated, depending on the balance of Na⁺ channel inactivation gate recovery, Na⁺ influx, and axial charge. Our findings suggest that slow Na⁺ removal influences axonal resilience to high-frequency spike propagation and that different strategies may be required to overcome this constraint in different neurons.

Axons reliably conduct action potentials between neurons and/or other targets. Axons have widely variable diameters and can be myelinated or unmyelinated. Although the effect of these factors on propagation speed is well studied, how they constrain axonal resilience to high-frequency spiking is incompletely understood. Maximal firing frequencies range from ∼1 Hz to >300 Hz across neurons, but the process by which Na/K pumps counteract Na⁺ influx is slow, and the extent to which slow Na⁺ removal is compatible with high-frequency spiking is unclear. Modeling the process of Na⁺ removal shows that large-diameter axons are more resilient to high-frequency spikes than are small-diameter axons, because of their slow Na⁺ accumulation. In myelinated axons, the myelinated compartments between nodes of Ranvier act as a “reservoir” to slow Na⁺ accumulation and increase the reliability of axonal propagation. We now find that slowing the activation of K⁺ current can increase the Na⁺ influx rate, and the effect of minimizing the overlap between Na⁺ and K⁺ currents on spike propagation resilience depends on complex interactions among diameter, myelination, and the Na/K pump density. Our results suggest that, in neurons with different channel gating kinetic parameters, different strategies may be required to improve the reliability of axonal propagation.

Comodulation of h- and Na+/K+ Pump Currents Expands the Range of Functional Bursting in a Central Pattern Generator by Navigating between Dysfunctional Regimes

Central pattern generators (CPGs), specialized oscillatory neuronal networks controlling rhythmic motor behaviors such as breathing and locomotion, must adjust their patterns of activity to a variable environment and changing behavioral goals. Neuromodulation adjusts these patterns by orchestrating changes in multiple ionic currents. In the medicinal leech, the endogenous neuromodulator myomodulin speeds up the heartbeat CPG by reducing the electrogenic Na+/K+ pump current and increasing h-current in pairs of mutually inhibitory leech heart interneurons (HNs), which form half-center oscillators (HN HCOs). Here we investigate whether the comodulation of two currents could have advantages over a single current in the control of functional bursting patterns of a CPG. We use a conductance-based biophysical model of an HN HCO to explain the experimental effects of myomodulin. We demonstrate that, in the model, comodulation of the Na+/K+ pump current and h-current expands the range of functional bursting activity by avoiding transitions into nonfunctional regimes, such as asymmetric bursting and plateau-containing seizure-like activity. We validate the model by finding parameters that reproduce temporal bursting characteristics matching experimental recordings from HN HCOs under control, three different myomodulin concentrations, and Cs+ treated conditions. The matching cases are located along the border of an asymmetric regime away from the border with more dangerous seizure-like activity. We found a simple comodulation mechanism with an inverse relation between the pump and h-currents makes a good fit of the matching cases and comprises a general mechanism for the robust and flexible control of oscillatory neuronal networks.SIGNIFICANCE STATEMENT Rhythm-generating neuronal circuits adjust their oscillatory patterns to accommodate a changing environment through neuromodulation. In different species, chemical messengers participating in such processes may target two or more membrane currents. In medicinal leeches, the neuromodulator myomodulin speeds up the heartbeat central pattern generator by reducing Na+/K+ pump current and increasing h-current. In a computational model, we show that this comodulation expands the range of central pattern generator's functional activity by navigating the circuit between dysfunctional regimes resulting in a much wider range of cycle period. This control would not be attainable by modulating only one current, emphasizing the synergy of combined effects. Given the prevalence of h-current and Na+/K+ pump current in neurons, similar comodulation mechanisms may exist across species.

An electrodiffusive neuron-extracellular-glia model for exploring the genesis of slow potentials in the brain

Within the computational neuroscience community, there has been a focus on simulating the electrical activity of neurons, while other components of brain tissue, such as glia cells and the extracellular space, are often neglected. Standard models of extracellular potentials are based on a combination of multicompartmental models describing neural electrodynamics and volume conductor theory. Such models cannot be used to simulate the slow components of extracellular potentials, which depend on ion concentration dynamics, and the effect that this has on extracellular diffusion potentials and glial buffering currents. We here present the electrodiffusive neuron-extracellular-glia (edNEG) model, which we believe is the first model to combine compartmental neuron modeling with an electrodiffusive framework for intra- and extracellular ion concentration dynamics in a local piece of neuro-glial brain tissue. The edNEG model (i) keeps track of all intraneuronal, intraglial, and extracellular ion concentrations and electrical potentials, (ii) accounts for action potentials and dendritic calcium spikes in neurons, (iii) contains a neuronal and glial homeostatic machinery that gives physiologically realistic ion concentration dynamics, (iv) accounts for electrodiffusive transmembrane, intracellular, and extracellular ionic movements, and (v) accounts for glial and neuronal swelling caused by osmotic transmembrane pressure gradients. The edNEG model accounts for the concentration-dependent effects on ECS potentials that the standard models neglect. Using the edNEG model, we analyze these effects by splitting the extracellular potential into three components: one due to neural sink/source configurations, one due to glial sink/source configurations, and one due to extracellular diffusive currents. Through a series of simulations, we analyze the roles played by the various components and how they interact in generating the total slow potential. We conclude that the three components are of comparable magnitude and that the stimulus conditions determine which of the components that dominate.

Activity-mediated accumulation of potassium induces a switch in firing pattern and neuronal excitability type

During normal neuronal activity, ionic concentration gradients across a neuron’s membrane are often assumed to be stable. Prolonged spiking activity, however, can reduce transmembrane gradients and affect voltage dynamics. Based on mathematical modeling, we investigated the impact of neuronal activity on ionic concentrations and, consequently, the dynamics of action potential generation. We find that intense spiking activity on the order of a second suffices to induce changes in ionic reversal potentials and to consistently induce a switch from a regular to an intermittent firing mode. This transition is caused by a qualitative alteration in the system’s voltage dynamics, mathematically corresponding to a co-dimension-two bifurcation from a saddle-node on invariant cycle (SNIC) to a homoclinic orbit bifurcation (HOM). Our electrophysiological recordings in mouse cortical pyramidal neurons confirm the changes in action potential dynamics predicted by the models: (i) activity-dependent increases in intracellular sodium concentration directly reduce action potential amplitudes, an effect typically attributed solely to sodium channel inactivation; (ii) extracellular potassium accumulation switches action potential generation from tonic firing to intermittently interrupted output. Thus, individual neurons may respond very differently to the same input stimuli, depending on their recent patterns of activity and/or the current brain-state.

Ionic concentrations in the brain are not constant. We show that during intense neuronal activity, they can change on the order of seconds and even switch neuronal spiking patterns under identical stimulation from a regular firing mode to an intermittently interrupted one. Triggered by an accumulation of extracellular potassium, such a transition is caused by a specific, qualitative change in of the neuronal voltage dynamics—a so-called bifurcation—which affects crucial features of action-potential generation and bears consequences for how information is encoded and how neurons behave together in the network. Also, changes in intracellular sodium can induce measurable effects, like a reduction of spike amplitude that occurs independently of the fast amplitude effects attributed to sodium channel inactivation. Taken together, our results demonstrate that a neuron can respond very differently to the same stimulus, depending on its previous activity or the current brain state. This finding may be particularly relevant when other regulatory mechanisms of ionic homeostasis are challenged, for example, during pathological states of glial impairment or oxygen deprivation. Finally, categorization of cortical neurons as intrinsically bursting or regular spiking may be biased by the ionic concentrations at the time of the observation, highlighting the non-static nature of neuronal dynamics.

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

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

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

Glial Modulation of Electrical Rhythms in a Neuroglial Network Model of Epilepsy

OBJECTIVE

An important EEG-based biomarker for epilepsy is the phase-amplitude cross-frequency coupling (PAC) of electrical rhythms; however, the underlying pathways of these pathologic markers are not always clear. Since glial cells have been shown to play an active role in neuroglial networks, it is likely that some of these PAC markers are modulated via glial effects.

METHODS

We developed a 4-unit hybrid model of a neuroglial network, consisting of 16 sub-units, that combines a mechanistic representation of neurons with an oscillator-based Cognitive Rhythm Generator (CRG) representation of glial cells-astrocytes and microglia. The model output was compared with recorded generalized tonic-clonic patient data, both in terms of PAC features, and state classification using an unsupervised hidden Markov model (HMM).

RESULTS

The neuroglial model output showed PAC features similar to those observed in epileptic seizures. These generated PAC features were able to accurately identify spontaneous epileptiform discharges (SEDs) as seizure-like states, as well as a postictal-like state following the long-duration SED, when applied to the HMM machine learning algorithm trained on patient data. The evolution profile of the maximal PAC during the SED compared well with patient data, showing similar association with the duration of the postictal state.

CONCLUSION

The hybrid neuroglial network model was able to generate PAC features similar to those observed in ictal and postictal epileptic states, which has been used for state classification and postictal state duration prediction.

SIGNIFICANCE

Since PAC biomarkers are important for epilepsy research and postictal state duration has been linked with risk of sudden unexplained death in epilepsy, this model suggests glial synaptic effects as potential targets for further analysis and treatment.

Dynamics of a neuron-glia system: the occurrence of seizures and the influence of electroconvulsive stimuli : A mathematical and numerical study

In this paper, we investigate the dynamics of a neuron-glia cell system and the underlying mechanism for the occurrence of seizures. For our mathematical and numerical investigation of the cell model we will use bifurcation analysis and some computational methods. It turns out that an increase of the potassium concentration in the reservoir is one trigger for seizures and is related to a torus bifurcation. In addition, we will study potassium dynamics of the model by considering a reduced version and we will show how both mechanisms are linked to each other. Moreover, the reduction of the potassium leak current will also induce seizures. Our study will show that an enhancement of the extracellular potassium concentration, which influences the Nernst potential of the potassium current, may lead to seizures. Furthermore, we will show that an external forcing term (e.g. electroshocks as unidirectional rectangular pulses also known as electroconvulsive therapy) will establish seizures similar to the unforced system with the increased extracellular potassium concentration. To this end, we describe the unidirectional rectangular pulses as an autonomous system of ordinary differential equations. These approaches will explain the appearance of seizures in the cellular model. Moreover, seizures, as they are measured by electroencephalography (EEG), spread on the macro-scale (cm). Therefore, we extend the cell model with a suitable homogenised monodomain model, propose a set of (numerical) experiment to complement the bifurcation analysis performed on the single-cell model. Based on these experiments, we introduce a bidomain model for a more realistic modelling of white and grey matter of the brain. Performing similar (numerical) experiment as for the monodomain model leads to a suitable comparison of both models. The individual cell model, with its seizures explained in terms of a torus bifurcation, extends directly to corresponding results in both the monodomain and bidomain models where the neural firing spreads almost synchronous through the domain as fast traveling waves, for physiologically relevant paramenters.

An electrodiffusive, ion conserving Pinsky-Rinzel model with homeostatic mechanisms

In most neuronal models, ion concentrations are assumed to be constant, and effects of concentration variations on ionic reversal potentials, or of ionic diffusion on electrical potentials are not accounted for. Here, we present the electrodiffusive Pinsky-Rinzel (edPR) model, which we believe is the first multicompartmental neuron model that accounts for electrodiffusive ion concentration dynamics in a way that ensures a biophysically consistent relationship between ion concentrations, electrical charge, and electrical potentials in both the intra- and extracellular space. The edPR model is an expanded version of the two-compartment Pinsky-Rinzel (PR) model of a hippocampal CA3 neuron. Unlike the PR model, the edPR model includes homeostatic mechanisms and ion-specific leakage currents, and keeps track of all ion concentrations (Na⁺, K⁺, Ca²⁺, and Cl⁻), electrical potentials, and electrical conductivities in the intra- and extracellular space. The edPR model reproduces the membrane potential dynamics of the PR model for moderate firing activity. For higher activity levels, or when homeostatic mechanisms are impaired, the homeostatic mechanisms fail in maintaining ion concentrations close to baseline, and the edPR model diverges from the PR model as it accounts for effects of concentration changes on neuronal firing. We envision that the edPR model will be useful for the field in three main ways. Firstly, as it relaxes commonly made modeling assumptions, the edPR model can be used to test the validity of these assumptions under various firing conditions, as we show here for a few selected cases. Secondly, the edPR model should supplement the PR model when simulating scenarios where ion concentrations are expected to vary over time. Thirdly, being applicable to conditions with failed homeostasis, the edPR model opens up for simulating a range of pathological conditions, such as spreading depression or epilepsy.

Neurons generate their electrical signals by letting ions pass through their membranes. Despite this fact, most models of neurons apply the simplifying assumption that ion concentrations remain effectively constant during neural activity. This assumption is often quite good, as neurons contain a set of homeostatic mechanisms that make sure that ion concentrations vary quite little under normal circumstances. However, under some conditions, these mechanisms can fail, and ion concentrations can vary quite dramatically. Standard models are thus not able to simulate such conditions. Here, we present what to our knowledge is the first multicompartmental neuron model that accounts for ion concentration variations in a way that ensures complete and consistent ion concentration and charge conservation. In this work, we use the model to explore under which activity conditions the ion concentration variations become important for predicting the neurodynamics. We expect the model to be of great value for the field of neuroscience, as it can be used to simulate a range of pathological conditions, such as spreading depression or epilepsy, which are associated with large changes in extracellular ion concentrations.

Ionic and synaptic mechanisms of seizure generation and epileptogenesis

The biophysical mechanisms underlying epileptogenesis and the generation of seizures remain to be better understood. Among many factors triggering epileptogenesis are traumatic brain injury breaking normal synaptic homeostasis and genetic mutations disrupting ionic concentration homeostasis. Impairments in these mechanisms, as seen in various brain diseases, may push the brain network to a pathological state characterized by increased susceptibility to unprovoked seizures. Here, we review recent computational studies exploring the roles of ionic concentration dynamics in the generation, maintenance, and termination of seizures. We further discuss how ionic and synaptic homeostatic mechanisms may give rise to conditions which prime brain networks to exhibit recurrent spontaneous seizures and epilepsy.

Regulation of Neuronal Na+/K+-ATPase by Specific Protein Kinases and Protein Phosphatases

The Na⁺/K⁺-ATPase (NKA) is a ubiquitous membrane-bound enzyme responsible for generating and maintaining the Na⁺ and K⁺ electrochemical gradients across the plasmalemma of living cells. Numerous studies in non-neuronal tissues have shown that this transport mechanism is reversibly regulated by phosphorylation/dephosphorylation of the catalytic α subunit and/or associated proteins. In neurons, Na⁺/K⁺ transport by NKA is essential for almost all neuronal operations, consuming up to two-thirds of the neuron's energy expenditure. However, little is known about its cellular regulatory mechanisms. Here we have used an electrophysiological approach to monitor NKA transport activity in male rat hippocampal neurons in situ We report that this activity is regulated by a balance between serine/threonine phosphorylation and dephosphorylation. Phosphorylation by the protein kinases PKG and PKC inhibits NKA activity, whereas dephosphorylation by the protein phosphatases PP-1 and PP-2B (calcineurin) reverses this effect. Given that these kinases and phosphatases serve as downstream effectors in key neuronal signaling pathways, they may mediate the coupling of primary messengers, such as neurotransmitters, hormones, and growth factors, to the NKAs, through which multiple brain functions can be regulated or dysregulated.SIGNIFICANCE STATEMENT The Na⁺/K⁺-ATPase (NKA), known as the "Na⁺ pump," is a ubiquitous membrane-bound enzyme responsible for generating and maintaining the Na⁺ and K⁺ electrochemical gradients across the plasma membrane of living cells. In neurons, as in most types of cells, the NKA generates the negative resting membrane potential, which is the basis for almost all aspects of cellular function. Here we used an electrophysiological approach to monitor physiological NKA transport activity in single hippocampal pyramidal cells in situ We have found that neuronal NKA activity is oppositely regulated by phosphorylation and dephosphorylation, and we have identified the main protein kinases and phosphatases mediating this regulation. This fundamental form of NKA regulation likely plays a role in multiple brain functions.

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.

Germline De Novo Mutations in ATP1A1 Cause Renal Hypomagnesemia, Refractory Seizures, and Intellectual Disability

Over the last decades, a growing spectrum of monogenic disorders of human magnesium homeostasis has been clinically characterized, and genetic studies in affected individuals have identified important molecular components of cellular and epithelial magnesium transport. Here, we describe three infants who are from non-consanguineous families and who presented with a disease phenotype consisting of generalized seizures in infancy, severe hypomagnesemia, and renal magnesium wasting. Seizures persisted despite magnesium supplementation and were associated with significant intellectual disability. Whole-exome sequencing and conventional Sanger sequencing identified heterozygous de novo mutations in the catalytic Na⁺, K⁺-ATPase α1 subunit (ATP1A1). Functional characterization of mutant Na⁺, K⁺-ATPase α1 subunits in heterologous expression systems revealed not only a loss of Na⁺, K⁺-ATPase function but also abnormal cation permeabilities, which led to membrane depolarization and possibly aggravated the effect of the loss of physiological pump activity. These findings underline the indispensable role of the α1 isoform of the Na⁺, K⁺-ATPase for renal-tubular magnesium handling and cellular ion homeostasis, as well as maintenance of physiologic neuronal activity.

Bottom-up approach to torus bifurcation in neuron models

We study the quasi-periodicity phenomena occurring at the transition between tonic spiking and bursting activities in exemplary biologically plausible Hodgkin-Huxley type models of individual cells and reduced phenomenological models with slow and fast dynamics. Using the geometric slow-fast dissection and the parameter continuation approach, we show that the transition is due to either the torus bifurcation or the period-doubling bifurcation of a stable periodic orbit on the 2D slow-motion manifold near a characteristic fold. Various torus bifurcations including stable and saddle torus-canards, resonant tori, the co-existence of nested tori, and the torus breakdown leading to the onset of complex and bistable dynamics in such systems are examined too.

Origin of slow spontaneous resting-state neuronal fluctuations in brain networks

Resting- or baseline-state low-frequency (0.01-0.2 Hz) brain activity is observed in fMRI, EEG, and local field potential recordings. These fluctuations were found to be correlated across brain regions and are thought to reflect neuronal activity fluctuations between functionally connected areas of the brain. However, the origin of these infra-slow resting-state fluctuations remains unknown. Here, using a detailed computational model of the brain network, we show that spontaneous infra-slow (<0.05 Hz) activity could originate due to the ion concentration dynamics. The computational model implemented dynamics for intra- and extracellular K⁺ and Na⁺ and intracellular Cl^- ions, Na⁺/K⁺ exchange pump, and KCC2 cotransporter. In the network model simulating resting awake-like brain state, we observed infra-slow fluctuations in the extracellular K⁺ concentration, Na⁺/K⁺ pump activation, firing rate of neurons, and local field potentials. Holding K⁺ concentration constant prevented generation of the infra-slow fluctuations. The amplitude and peak frequency of this activity were modulated by the Na⁺/K⁺ pump, AMPA/GABA synaptic currents, and glial properties. Further, in a large-scale network with long-range connections based on CoCoMac connectivity data, the infra-slow fluctuations became synchronized among remote clusters similar to the resting-state activity observed in vivo. Overall, our study proposes that ion concentration dynamics mediated by neuronal and glial activity may contribute to the generation of very slow spontaneous fluctuations of brain activity that are reported as the resting-state fluctuations in fMRI and EEG recordings.

Minimal model of interictal and ictal discharges "Epileptor-2"

Seizures occur in a recurrent manner with intermittent states of interictal and ictal discharges (IIDs and IDs). The transitions to and from IDs are determined by a set of processes, including synaptic interaction and ionic dynamics. Although mathematical models of separate types of epileptic discharges have been developed, modeling the transitions between states remains a challenge. A simple generic mathematical model of seizure dynamics (Epileptor) has recently been proposed by Jirsa et al. (2014); however, it is formulated in terms of abstract variables. In this paper, a minimal population-type model of IIDs and IDs is proposed that is as simple to use as the Epileptor, but the suggested model attributes physical meaning to the variables. The model is expressed in ordinary differential equations for extracellular potassium and intracellular sodium concentrations, membrane potential, and short-term synaptic depression variables. A quadratic integrate-and-fire model driven by the population input current is used to reproduce spike trains in a representative neuron. In simulations, potassium accumulation governs the transition from the silent state to the state of an ID. Each ID is composed of clustered IID-like events. The sodium accumulates during discharge and activates the sodium-potassium pump, which terminates the ID by restoring the potassium gradient and thus polarizing the neuronal membranes. The whole-cell and cell-attached recordings of a 4-AP-based in vitro model of epilepsy confirmed the primary model assumptions and predictions. The mathematical analysis revealed that the IID-like events are large-amplitude stochastic oscillations, which in the case of ID generation are controlled by slow oscillations of ionic concentrations. The IDs originate in the conditions of elevated potassium concentrations in a bath solution via a saddle-node-on-invariant-circle-like bifurcation for a non-smooth dynamical system. By providing a minimal biophysical description of ionic dynamics and network interactions, the model may serve as a hierarchical base from a simple to more complex modeling of seizures.

In pathological conditions of epilepsy, the functioning of the neural network crucially depends on the ionic concentrations inside and outside neurons. A number of factors that affect neuronal activity is large. That is why the development of a minimal model that reproduces typical seizures could structure further experimental and analytical studies of the pathological mechanisms. Here, on a base of known biophysical models, we present a simple population-type model that includes only four principal variables, the extracellular potassium concentration, the intracellular sodium concentration, the membrane potential and the synaptic resource diminishing due to short-term synaptic depression. A simple modeled neuron is used as an observer of the population activity. We validate the model assumptions with in vitro experiments. Our model reproduces ictal and interictal events, where the latter result in bursts of spikes in single neurons, and the former represent the cluster of spike bursts. Mathematical analysis reveals that the bursts are spontaneous large-amplitude oscillations, which may cluster after a saddle-node on invariant circle bifurcation in the pro-epileptic conditions. Our consideration has significant bearing in understanding pathological neuronal network dynamics.

Computational model of brain-stem circuit for state-dependent control of hypoglossal motoneurons

In patients with obstructive sleep apnea (OSA), the pharyngeal muscles become relaxed during sleep, which leads to a partial or complete closure of upper airway. Experimental studies suggest that withdrawal of noradrenergic and serotonergic drives importantly contributes to depression of hypoglossal motoneurons and, therefore, may contribute to OSA pathophysiology; however, specific cellular and synaptic mechanisms remain unknown. In this new study, we developed a biophysical network model to test the hypothesis that, to explain experimental observations, the neuronal network for monoaminergic control of excitability of hypoglossal motoneurons needs to include excitatory and inhibitory perihypoglossal interneurons that mediate noradrenergic and serotonergic drives to hypoglossal motoneurons. In the model, the state-dependent activation of the hypoglossal motoneurons was in qualitative agreement with in vivo data during simulated rapid eye movement (REM) and non-REM sleep. The model was applied to test the mechanisms of action of noradrenergic and serotonergic drugs during REM sleep as observed in vivo. We conclude that the proposed minimal neuronal circuit is sufficient to explain in vivo data and supports the hypothesis that perihypoglossal interneurons may mediate state-dependent monoaminergic drive to hypoglossal motoneurons. The population of the hypothesized perihypoglossal interneurons may serve as novel targets for pharmacological treatment of OSA. NEW & NOTEWORTHY In vivo studies suggest that during rapid eye movement sleep, withdrawal of noradrenergic and serotonergic drives critically contributes to depression of hypoglossal motoneurons (HMs), which innervate the tongue muscles. By means of a biophysical model, which is consistent with a broad range of empirical data, we demonstrate that the neuronal network controlling the excitability of HMs needs to include excitatory and inhibitory interneurons that mediate noradrenergic and serotonergic drives to HMs.

A new water-soluble polythiophene derivative as a probe for real-time monitoring adenosine 5'-triphosphatase activity in lysosome of living cells

Detection of the adenosine 5'-triphosphatase (ATPase) activity in lysosome of living cells is of great importance for clinical diagnosis of many related diseases, including cancer. In this work, a new water-soluble polythiophene derivative named ZnPT bearing both quaternary ammonium salt groups and dipicolylamine-Zn²⁺ (DPA-Zn²⁺) complexes in its side chain, was designed and synthesized for this propose. The probe mainly localized to lysosome with good biocompatibility and membrane penetration. The real-time, continuous, direct, and label-free assays were achieved through a fluorescence "turn-on" mode by taking advantages of the reaction specificity of ATPase with ATP and the high binding selectivity of ZnPT toward ATP substrate over its hydrolysis product (ADP). This well designed strategy should provide a facile and effective way for investigating ATPase-relevant biological processes.

Propensity for Bistability of Bursting and Silence in the Leech Heart Interneuron

The coexistence of neuronal activity regimes has been reported under normal and pathological conditions. Such multistability could enhance the flexibility of the nervous system and has many implications for motor control, memory, and decision making. Multistability is commonly promoted by neuromodulation targeting specific membrane ionic currents. Here, we investigated how modulation of different ionic currents could affect the neuronal propensity for bistability. We considered a leech heart interneuron model. It exhibits bistability of bursting and silence in a narrow range of the leak current parameters, conductance (g_leak) and reversal potential (E_leak). We assessed the propensity for bistability of the model by using bifurcation diagrams. On the diagram (g_leak, E_leak), we mapped bursting and silent regimes. For the canonical value of E_leak we determined the range of g_leak which supported the bistability. We use this range as an index of propensity for bistability. We investigated how this index was affected by alterations of ionic currents. We systematically changed their conductances, one at a time, and built corresponding bifurcation diagrams in parameter planes of the maximal conductance of a given current and the leak conductance. We found that conductance of only one current substantially affected the index of propensity; the increase of the maximal conductance of the hyperpolarization-activated cationic current increased the propensity index. The second conductance with the strongest effect was the conductance of the low-threshold fast Ca²⁺ current; its reduction increased the propensity index although the effect was about two times smaller in magnitude. Analyzing the model with both changes applied simultaneously, we found that the diagram (g_leak, E_leak) showed a progressively expanded area of bistability of bursting and silence.

Role of KCC2-dependent potassium efflux in 4-Aminopyridine-induced Epileptiform synchronization

A balance between excitation and inhibition is necessary to maintain stable brain network dynamics. Traditionally, seizure activity is believed to arise from the breakdown of this delicate balance in favor of excitation with loss of inhibition. Surprisingly, recent experimental evidence suggests that this conventional view may be limited, and that inhibition plays a prominent role in the development of epileptiform synchronization. Here, we explored the role of the KCC2 co-transporter in the onset of inhibitory network-induced seizures. Our experiments in acute mouse brain slices, of either sex, revealed that optogenetic stimulation of either parvalbumin- or somatostatin-expressing interneurons induced ictal discharges in rodent entorhinal cortex during 4-aminopyridine application. These data point to a proconvulsive role of GABAA receptor signaling that is independent of the inhibitory input location (i.e., dendritic vs. somatic). We developed a biophysically realistic network model implementing dynamics of ion concentrations to explore the mechanisms leading to inhibitory network-induced seizures. In agreement with experimental results, we found that stimulation of the inhibitory interneurons induced seizure-like activity in a network with reduced potassium A-current. Our model predicts that interneuron stimulation triggered an increase of interneuron firing, which was accompanied by an increase in the intracellular chloride concentration and a subsequent KCC2-dependent gradual accumulation of the extracellular potassium promoting epileptiform ictal activity. When the KCC2 activity was reduced, stimulation of the interneurons was no longer able to induce ictal events. Overall, our study provides evidence for a proconvulsive role of GABAA receptor signaling that depends on the involvement of the KCC2 co-transporter.

Chaotic Dynamics Mediate Brain State Transitions, Driven by Changes in Extracellular Ion Concentrations

Previous studies have suggested that changes in extracellular ion concentrations initiate the transition from an activity state that characterizes sleep in cortical neurons to states that characterize wakefulness. However, because neuronal activity and extracellular ion concentrations are interdependent, isolating their unique roles during sleep-wake transitions is not possible in vivo. Here, we extend the Averaged-Neuron model and demonstrate that, although changes in extracellular ion concentrations occur concurrently, decreasing the conductance of calcium-dependent potassium channels initiates the transition from sleep to wakefulness. We find that sleep is governed by stable, self-sustained oscillations in neuronal firing patterns, whereas the quiet awake state and active awake state are both governed by irregular oscillations and chaotic dynamics; transitions between these separable awake states are prompted by ionic changes. Although waking is indicative of a shift from stable to chaotic neuronal firing patterns, we illustrate that the properties of chaotic dynamics ensure that the transition between states is smooth and robust to noise.

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.

Low-to-High Cross-Frequency Coupling in the Electrical Rhythms as Biomarker for Hyperexcitable Neuroglial Networks of the Brain

OBJECTIVE

One of the features used in the study of hyperexcitablility is high-frequency oscillations (HFOs, >80 Hz). HFOs have been reported in the electrical rhythms of the brain's neuroglial networks under physiological and pathological conditions. Cross-frequency coupling (CFC) of HFOs with low-frequency rhythms was used to identify pathologic HFOs in the epileptogenic zones of epileptic patients and as a biomarker for the severity of seizure-like events in genetically modified rodent models. We describe a model to replicate reported CFC features extracted from recorded local field potentials (LFPs) representing network properties.

METHODS

This study deals with a four-unit neuroglial cellular network model where each unit incorporates pyramidal cells, interneurons, and astrocytes. Three different pathways of hyperexcitability generation-Na - ATPase pump, glial potassium clearance, and potassium afterhyperpolarization channel-were used to generate LFPs. Changes in excitability, average spontaneous electrical discharge (SED) duration, and CFC were then measured and analyzed.

RESULTS

Each parameter caused an increase in network excitability and the consequent lengthening of the SED duration. Short SEDs showed CFC between HFOs and theta oscillations (4-8 Hz), but in longer SEDs the low frequency changed to the delta range (1-4 Hz).

CONCLUSION

Longer duration SEDs exhibit CFC features similar to those reported by our team.

SIGNIFICANCE

First, Identifying the exponential relationship between network excitability and SED durations; second, highlighting the importance of glia in hyperexcitability (as they relate to extracellular potassium); and third, elucidation of the biophysical basis for CFC coupling features.

Sodium Pumps Mediate Activity-Dependent Changes in Mammalian Motor Networks

Ubiquitously expressed sodium pumps are best known for maintaining the ionic gradients and resting membrane potential required for generating action potentials. However, activity- and state-dependent changes in pump activity can also influence neuronal firing and regulate rhythmic network output. Here we demonstrate that changes in sodium pump activity regulate locomotor networks in the spinal cord of neonatal mice. The sodium pump inhibitor, ouabain, increased the frequency and decreased the amplitude of drug-induced locomotor bursting, effects that were dependent on the presence of the neuromodulator dopamine. Conversely, activating the pump with the sodium ionophore monensin decreased burst frequency. When more "natural" locomotor output was evoked using dorsal-root stimulation, ouabain increased burst frequency and extended locomotor episode duration, whereas monensin slowed and shortened episodes. Decreasing the time between dorsal-root stimulation, and therefore interepisode interval, also shortened and slowed activity, suggesting that pump activity encodes information about past network output and contributes to feedforward control of subsequent locomotor bouts. Using whole-cell patch-clamp recordings from spinal motoneurons and interneurons, we describe a long-duration (∼60 s), activity-dependent, TTX- and ouabain-sensitive, hyperpolarization (∼5 mV), which is mediated by spike-dependent increases in pump activity. The duration of this dynamic pump potential is enhanced by dopamine. Our results therefore reveal sodium pumps as dynamic regulators of mammalian spinal motor networks that can also be affected by neuromodulatory systems. Given the involvement of sodium pumps in movement disorders, such as amyotrophic lateral sclerosis and rapid-onset dystonia parkinsonism, knowledge of their contribution to motor network regulation also has considerable clinical importance.

SIGNIFICANCE STATEMENT

The sodium pump is ubiquitously expressed and responsible for at least half of total brain energy consumption. The pumps maintain ionic gradients and the resting membrane potential of neurons, but increasing evidence suggests that activity- and state-dependent changes in pump activity also influence neuronal firing. Here we demonstrate that changes in sodium pump activity regulate locomotor output in the spinal cord of neonatal mice. We describe a sodium pump-mediated afterhyperpolarization in spinal neurons, mediated by spike-dependent increases in pump activity, which is affected by dopamine. Understanding how sodium pumps contribute to network regulation and are targeted by neuromodulators, including dopamine, has clinical relevance due to the role of the sodium pump in diseases, including amyotrophic lateral sclerosis, parkinsonism, epilepsy, and hemiplegic migraine.

Brain sodium MRI in human epilepsy: Disturbances of ionic homeostasis reflect the organization of pathological regions

In light of technical advancements supporting exploration of MR signals other than ¹H, sodium (²³Na) has received attention as a marker of ionic homeostasis and cell viability. Here, we evaluate for the first time the possibility that ²³Na-MRI is sensitive to pathological processes occurring in human epilepsy. A normative sample of 27 controls was used to normalize regions of interest (ROIs) from 1424 unique brain locales on quantitative ²³Na-MRI and high-resolution ¹H-MPRAGE images. ROIs were based on intracerebral electrodes in ten patients undergoing epileptic network mapping. The stereo-EEG gold standard was used to define regions as belonging to primarily epileptogenic, secondarily irritative and to non-involved regions. Estimates of total sodium concentration (TSC) on ²³Na-MRI and cerebrospinal fluid (CSF) on ¹H imaging were extracted for each patient ROI, and normalized against the same region in controls. ROIs with disproportionate CSF contributions (Z_CSF≥1.96) were excluded. TSC levels were found to be elevated in patients relative to controls except in one patient, who suffered non-convulsive seizures during the scan, in whom we found reduced TSC levels. In the remaining patients, an ANOVA (F₁₁₀₀= 12.37, p<0.0001) revealed a highly significant effect of clinically-defined zones (F₁₁₀₀= 11.13, p<0.0001), with higher normalized TSC in the epileptogenic zone relative to both secondarily irritative (F₁₁₀₀= 11, p=0.0009) and non-involved regions (F₁₁₀₀= 17.8, p<0.0001). We provide the first non-invasive, in vivo evidence of a chronic TSC elevation alongside Z_CSF levels within the normative range, associated with the epileptogenic region even during the interictal period in human epilepsy, and the possibility of reduced TSC levels due to seizure. In line with modified homeostatic mechanisms in epilepsy - including altered mechanisms underlying ionic gating, clearance and exchange - we provide the first indication of ²³Na-MRI as an assay of altered sodium concentrations occurring in epilepsy associated with the organization of clinically relevant divisions of pathological cortex.

Sodium pump regulation of locomotor control circuits

Sodium pumps are ubiquitously expressed membrane proteins that extrude three Na⁺ ions in exchange for two K⁺ ions, using ATP as an energy source. Recent studies have illuminated additional, dynamic roles for sodium pumps in regulating the excitability of neuronal networks in an activity-dependent fashion. We review their role in a novel form of short-term memory within rhythmic locomotor networks. The data we review derives mainly from recent studies on Xenopus tadpoles and neonatal mice. The role and underlying mechanisms of pump action broadly match previously published data from an invertebrate, the Drosophila larva. We therefore propose a highly conserved mechanism by which sodium pump activity increases following a bout of locomotion. This results in an ultraslow afterhyperpolarization (usAHP) of the membrane potential that lasts around 1 min, but which only occurs in around half the network neurons. This usAHP in turn alters network excitability so that network output is reduced in a locomotor interval-dependent manner. The pumps therefore confer on spinal locomotor networks a temporary memory trace of recent network performance.

Na(+)/K(+) pump interacts with the h-current to control bursting activity in central pattern generator neurons of leeches

The dynamics of different ionic currents shape the bursting activity of neurons and networks that control motor output. Despite being ubiquitous in all animal cells, the contribution of the Na(+)/K(+) pump current to such bursting activity has not been well studied. We used monensin, a Na(+)/H(+) antiporter, to examine the role of the pump on the bursting activity of oscillator heart interneurons in leeches. When we stimulated the pump with monensin, the period of these neurons decreased significantly, an effect that was prevented or reversed when the h-current was blocked by Cs(+). The decreased period could also occur if the pump was inhibited with strophanthidin or K(+)-free saline. Our monensin results were reproduced in model, which explains the pump's contributions to bursting activity based on Na(+) dynamics. Our results indicate that a dynamically oscillating pump current that interacts with the h-current can regulate the bursting activity of neurons and networks.

Evaluation of cortical local field potential diffusion in stereotactic electro-encephalography recordings: A glimpse on white matter signal

While there is a strong interest in meso-scale field potential recording using intracranial electroencephalography with penetrating depth electrodes (i.e. stereotactic EEG or S-EEG) in humans, the signal recorded in the white matter remains ignored. White matter is generally considered electrically neutral and often included in the reference montage. Moreover, re-referencing electrophysiological data is a critical preprocessing choice that could drastically impact signal content and consequently the results of any given analysis. In the present stereotactic electroencephalography study, we first illustrate empirically the consequences of commonly used references (subdermal, white matter, global average, local montage) on inter-electrode signal correlation. Since most of these reference montages incorporate white matter signal, we next consider the difference between signals recorded in cortical gray matter and white matter. Our results reveal that electrode contacts located in the white matter record a mixture of activity, with part arising from the volume conduction (zero time delay) of activity from nearby gray matter. Furthermore, our analysis shows that white matter signal may be correlated with distant gray matter signal. While residual passive electrical spread from nearby matter may account for this relationship, our results suggest the possibility that this long distance correlation arises from the white matter fiber tracts themselves (i.e. activity from distant gray matter traveling along axonal fibers with time lag larger than zero); yet definitive conclusions about the origin of the white matter signal would require further experimental substantiation. By characterizing the properties of signals recorded in white matter and in gray matter, this study illustrates the importance of including anatomical prior knowledge when analyzing S-EEG data.

Ion dynamics during seizures

Changes in membrane voltage brought about by ion fluxes through voltage and transmitter-gated channels represent the basis of neural activity. As such, electrochemical gradients across the membrane determine the direction and driving force for the flow of ions and are therefore crucial in setting the properties of synaptic transmission and signal propagation. Ion concentration gradients are established by a variety of mechanisms, including specialized transporter proteins. However, transmembrane gradients can be affected by ionic fluxes through channels during periods of elevated neural activity, which in turn are predicted to influence the properties of on-going synaptic transmission. Such activity-induced changes to ion concentration gradients are a feature of both physiological and pathological neural processes. An epileptic seizure is an example of severely perturbed neural activity, which is accompanied by pronounced changes in intracellular and extracellular ion concentrations. Appreciating the factors that contribute to these ion dynamics is critical if we are to understand how a seizure event evolves and is sustained and terminated by neural tissue. Indeed, this issue is of significant clinical importance as status epilepticus—a type of seizure that does not stop of its own accord—is a life-threatening medical emergency. In this review we explore how the transmembrane concentration gradient of the six major ions (K⁺, Na⁺, Cl⁻, Ca²⁺, H⁺and HCO3−) is altered during an epileptic seizure. We will first examine each ion individually, before describing how multiple interacting mechanisms between ions might contribute to concentration changes and whether these act to prolong or terminate epileptic activity. In doing so, we will consider how the availability of experimental techniques has both advanced and restricted our ability to study these phenomena.