1
|
Li Y, Gu H, Qi C. Uncommon and common roles of inhibitory interneuron and autapse and their cooperations to induce or eliminate epileptiform firing of pyramidal neuron. Cogn Neurodyn 2025; 19:59. [PMID: 40206239 PMCID: PMC11977076 DOI: 10.1007/s11571-025-10243-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2025] [Revised: 02/27/2025] [Accepted: 03/13/2025] [Indexed: 04/11/2025] Open
Abstract
Different from the common role of inhibitory modulations to suppress firing activities, uncommon roles of inhibitory modulations are observed in recent experiments. For instance, inhibitory autapse can enhance spiking frequency of interneuron, and inhibitory interneuron can enhance spiking of pyramidal neuron to epileptiform firing with high membrane potential and extracellular potassium concentration, presenting possible novel etiology of brain diseases and challenge to excitation-inhibition balance. In the present paper, the uncommon roles, the common roles, and their cooperations are studied in a computation model. Firstly, the inhibitory interneuron with fast instead of slow decay synaptic current plays an uncommon role, and the complex process for the uncommon role is obtained. Compared with slow decay, the fast decay inhibitory synaptic current is strong enough to induce silence with low membrane potential, resulting in long silence and high level of extracellular potassium concentration when firing recovers, initiating positive feedback between firing and potassium concentration to induce the epileptiform firing. Secondly, inhibitory autaptic current with fast rather than slow decay plays an uncommon role to enhance spiking frequency of interneuron. Autaptic current with slow decay causes weak potassium current during downstroke of action potential to induce spike advanced. Finally, different cooperations between the common and uncommon roles of interneuron and autapse are obtained. Especially, fast autapse with great uncommon role can reverse the common role of interneuron, which can induce spiking to the epileptiform firing, and slow autapse with great common role can reverse the uncommon role of interneuron, which can change the epileptiform firing to spiking for the normal state. These findings present explanations to the uncommon roles of inhibitory modulations and multiple feasible measures to modulate the epileptiform firing and brain diseases.
Collapse
Affiliation(s)
- Yuye Li
- College of Mathematics and Computer Science, Chifeng University, Chifeng, 024000 China
| | - Huaguang Gu
- School of Aerospace Engineering and Applied Mechanics, Tongji University, Shanghai, 200092 China
| | - Changsheng Qi
- College of Chemistry and Life Science, Chifeng University, Chifeng, 024000 China
| |
Collapse
|
2
|
Norby JH, Hummel D, Ricks N, Rolston J, Rahimpour S, Cowan R, Voipio J, Trevelyan AJ, Smith EH, Parrish RR. Rodent and human seizures demonstrate a dynamic interplay with spreading depolarizations. Neurobiol Dis 2025; 211:106937. [PMID: 40348201 PMCID: PMC12086139 DOI: 10.1016/j.nbd.2025.106937] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2025] [Revised: 04/17/2025] [Accepted: 04/29/2025] [Indexed: 05/14/2025] Open
Abstract
Seizure termination has been linked to spreading depolarizations (SDs) in experimental models of epilepsy, and SDs have been suggested to protect against seizures. However, the precise mechanisms of seizure-associated SDs remain unclear. Additionally, the co-occurrence of SDs with human seizures remains debated. In this study, we found that SDs are a prominent feature following ictal events in both human clinical recordings and in a rodent model of ictogenesis. Approximately one-third of rodent seizure-like events (SLEs) associated with SDs, while all human seizures analyzed associated with propagating infraslow shifts, indicative of SDs. In rodents, SDs clustered towards the end of ictal events, resulting in significantly shorter SLEs and delayed onset of subsequent SLEs. Interestingly, SLEs with SDs displayed significantly more low gamma activity during ictal events than SLEs that did not end in SDs. Furthermore, we found no significant correlation between [K+]o levels and the likelihood of SLEs ending in SDs, questioning the role of [K+]o in SD induction during seizures. Interestingly, the human data demonstrate clear SD propagation during seizures and show that SDs appear and propagate in multiple brain regions simultaneously with ictal events. Collectively, these results indicate that SDs are a hallmark of ictal activity associated with seizure termination. Furthermore, these findings provide unique insight into the neuronal dynamics that promote SD induction by showing that increased low gamma activity during SLEs is more predictive of SD induction than [K+]o levels. Taken together, these findings provide rationale for further exploration of SDs to prematurely terminate life-threatening seizures.
Collapse
Affiliation(s)
- Jacob H Norby
- Department of Cell Biology and Physiology, Brigham Young University, Provo, UT, USA; Neuroscience Center, Brigham Young University, Provo, UT, USA
| | - Daniel Hummel
- Department of Cell Biology and Physiology, Brigham Young University, Provo, UT, USA; Neuroscience Center, Brigham Young University, Provo, UT, USA
| | - Noah Ricks
- Department of Cell Biology and Physiology, Brigham Young University, Provo, UT, USA
| | | | - Shervin Rahimpour
- Department of Neurosurgery, University of Utah, Salt Lake City 84112, UT, USA
| | - Rhiannon Cowan
- Department of Neurosurgery, University of Utah, Salt Lake City 84112, UT, USA
| | - Juha Voipio
- Faculty of Biological and Environmental Sciences, Molecular and Integrative Biosciences, University of Helsinki, Helsinki 00014, Finland
| | - Andrew J Trevelyan
- Newcastle University Biosciences Institute, Medical School, Framlington Place, Newcastle upon Tyne NE2 4HH, UK
| | - Elliot H Smith
- Department of Neurosurgery, University of Utah, Salt Lake City 84112, UT, USA.
| | - R Ryley Parrish
- Department of Cell Biology and Physiology, Brigham Young University, Provo, UT, USA; Neuroscience Center, Brigham Young University, Provo, UT, USA.
| |
Collapse
|
3
|
Ho 何鎮宇 ECY, Newton AJH, Urdapilleta E, Dura-Bernal S, Truccolo W. Downmodulation of Potassium Conductances Induces Epileptic Seizures in Cortical Network Models Via Multiple Synergistic Factors. J Neurosci 2025; 45:e1909232025. [PMID: 39880680 PMCID: PMC11949479 DOI: 10.1523/jneurosci.1909-23.2025] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2023] [Revised: 01/17/2025] [Accepted: 01/22/2025] [Indexed: 01/31/2025] Open
Abstract
Voltage-gated potassium conductances g K play a critical role not only in normal neural function, but also in many neurological disorders and related therapeutic interventions. In particular, in an important animal model of epileptic seizures, 4-aminopyridine (4-AP) administration is thought to induce seizures by reducing g K in cortex and other brain areas. Interestingly, 4-AP has also been useful in the treatment of neurological disorders such as multiple sclerosis and spinal cord injury, where it is thought to improve action potential propagation in axonal fibers. Here, we examined g K downmodulation in biophysical models of cortical networks that included different neuron types organized in layers, potassium diffusion in interstitial and larger extracellular spaces, and glial buffering. Our findings are fourfold. First, g K downmodulation in pyramidal and fast-spiking inhibitory interneurons led to differential effects, making the latter much more likely to enter depolarization block. Second, both neuron types showed an increase in the duration and amplitude of action potentials, with more pronounced effects in pyramidal neurons. Third, a sufficiently strong g K reduction dramatically increased network synchrony, resulting in seizure-like dynamics. Fourth, we hypothesized that broader action potentials were likely to not only improve their propagation, as in 4-AP therapeutic uses, but also to increase synaptic coupling. Notably, graded-synapses incorporating this effect further amplified network synchronization and seizure-like dynamics. Overall, our findings elucidate different effects that g K downmodulation may have in cortical networks, explaining its potential role in both pathological neural dynamics and therapeutic applications.
Collapse
Affiliation(s)
- Ernest C Y Ho 何鎮宇
- Department of Neuroscience, Brown University, Providence, Rhode Island 02912
| | - Adam J H Newton
- Department of Physiology and Pharmacology, State University of New York (SUNY), Downstate Health Sciences University, Brooklyn, New York 11203
| | - Eugenio Urdapilleta
- Centro Atómico Bariloche and Instituto Balseiro, Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Universidad Nacional de Cuyo, R8402AGP Bariloche, Río Negro, Argentina
| | - Salvador Dura-Bernal
- Department of Physiology and Pharmacology, State University of New York (SUNY), Downstate Health Sciences University, Brooklyn, New York 11203
- Center for Biomedical Imaging and Neuromodulation, Nathan S. Kline Institute for Psychiatric Research, Orangeburg, New York 10962
| | - Wilson Truccolo
- Department of Neuroscience, Brown University, Providence, Rhode Island 02912
- Carney Institute for Brain Science, Brown University, Providence, Rhode Island 02912
| |
Collapse
|
4
|
Sato TR, Ode KL, Kinoshita FL, Ueda HR. A design principle for neuronal firing with up-down oscillation through Na + dynamics. iScience 2025; 28:111904. [PMID: 40028276 PMCID: PMC11869597 DOI: 10.1016/j.isci.2025.111904] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2024] [Revised: 12/10/2024] [Accepted: 01/23/2025] [Indexed: 03/05/2025] Open
Abstract
Nonrapid eye movement sleep is characterized by high-amplitude and low-frequency electroencephalography signals. These signals are thought to be produced by the synchronized activity of cortical neurons, demonstrating the alternating bursting (up) and resting (down) states. Here, such an activity is referred to as up-down oscillation (UDO). Previously, we discussed the importance of the Ca2+-dependent hyperpolarization pathway in the generation of UDO by simulating neuronal activity based on the Hodgkin-Huxley-type model. We herein focus on intracellular Na+ dynamics. The Na+-centered model indicates that the activation of voltage-gated Na+ channels leads to intracellular Na+ accumulation, which in turn activates Na+-dependent K+ (KNa) channels or Na+/K+ ATPases, resulting in the down state. Activation kinetics of voltage-gated Na+ channels are important in shaping the UDO firing. Therefore, our model demonstrates that voltage-gated Na+ and KNa channels or Na+/K+ ATPases are candidate pathways for UDO induction.
Collapse
Affiliation(s)
- Tomohide R. Sato
- Department of Systems Pharmacology, Graduate School of Medicine, The University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo 113-0033, Japan
| | - Koji L. Ode
- Department of Systems Pharmacology, Graduate School of Medicine, The University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo 113-0033, Japan
| | - Fukuaki L. Kinoshita
- Department of Neurology, Graduate School of Medicine, Osaka University, 2-2 Yamadaoka, Suita, Osaka 565-0871, Japan
- Laboratory for Synthetic Biology, RIKEN Center for Biosystems Dynamics Research, 1-3 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Hiroki R. Ueda
- Department of Systems Pharmacology, Graduate School of Medicine, The University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo 113-0033, Japan
- Laboratory for Synthetic Biology, RIKEN Center for Biosystems Dynamics Research, 1-3 Yamadaoka, Suita, Osaka 565-0871, Japan
- Department of Systems Biology, Institute of Life Science, Kurume University, 67 Asahimachi, Kurume, Fukuoka 830-0011, Japan
| |
Collapse
|
5
|
Pietrobon D, Brennan KC. Mechanisms underlying CSD initiation implicated by genetic mouse models of migraine. J Headache Pain 2025; 26:17. [PMID: 39871148 PMCID: PMC11773941 DOI: 10.1186/s10194-025-01948-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2024] [Accepted: 01/06/2025] [Indexed: 01/29/2025] Open
Abstract
A key unanswered question in migraine neurobiology concerns the mechanisms that make the brain of migraineurs susceptible to cortical spreading depression (CSD, a spreading depolarization that underlies migraine aura and may trigger the migraine pain mechanisms). Important insights into this question can be obtained by studying the mechanisms of facilitation of CSD initiation in genetic mouse models of the disease. These models, all generated from families with hereditary migraine, allow the investigation of the functional consequences of disease-causing mutations at the molecular, cellular, synaptic and neural circuit levels. In this review, after describing the available genetic mouse models of migraine, which all share increased susceptibility to experimentally induced CSD, we will discuss the functional alterations in their cerebral cortex and the mechanisms underlying the facilitation of CSD initiation in their cortex, as well as the insights that these mechanisms may give into the mechanisms of initiation of spontaneous CSDs in migraine.
Collapse
Affiliation(s)
- Daniela Pietrobon
- Department of Biomedical Sciences and Padova Neuroscience Center, University of Padova, Via Ugo Bassi 58, 35131, Padua, Italy.
| | - K C Brennan
- Department of Neurology, University of Utah, 383 Colorow Drive, Salt Lake City, UT, 84108, USA.
| |
Collapse
|
6
|
Liu Z, De Schutter E, Li Y. GABA-Induced Seizure-Like Events Caused by Multi-ionic Interactive Dynamics. eNeuro 2024; 11:ENEURO.0308-24.2024. [PMID: 39443111 PMCID: PMC11524612 DOI: 10.1523/eneuro.0308-24.2024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2024] [Accepted: 09/17/2024] [Indexed: 10/25/2024] Open
Abstract
Experimental evidence showed that an increase in intracellular chloride concentration [Formula: see text] caused by gamma-aminobutyric acid (GABA) input can promote epileptic firing activity, but the actual mechanisms remain elusive. Here in this theoretical work, we show that influx of chloride and concomitant bicarbonate ion [Formula: see text] efflux upon GABA receptor activation can induce epileptic firing activity by transition of GABA from inhibition to excitation. We analyzed the intrinsic property of neuron firing states as a function of [Formula: see text] We found that as [Formula: see text] increases, the system exhibits a saddle-node bifurcation, above which the neuron exhibits a spectrum of intensive firing, periodic bursting interrupted by depolarization block (DB) state, and eventually a stable DB through a Hopf bifurcation. We demonstrate that only GABA stimuli together with [Formula: see text] efflux can switch GABA's effect to excitation which leads to a series of seizure-like events (SLEs). Exposure to a low [Formula: see text] can drive neurons with high concentrations of [Formula: see text] downward to lower levels of [Formula: see text], during which it could also trigger SLEs depending on the exchange rate with the bath. Our analysis and simulation results show how the competition between GABA stimuli-induced accumulation of [Formula: see text] and [Formula: see text] application-induced decrease of [Formula: see text] regulates the neuron firing activity, which helps to understand the fundamental ionic dynamics of SLE.
Collapse
Affiliation(s)
- Zichao Liu
- School of Systems Science, Beijing Normal University, Beijing 100875, China
| | - Erik De Schutter
- Computational Neuroscience Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| | - Yinyun Li
- School of Systems Science, Beijing Normal University, Beijing 100875, China
- Computational Neuroscience Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| |
Collapse
|
7
|
Behbood M, Lemaire L, Schleimer JH, Schreiber S. The Na+/K+-ATPase generically enables deterministic bursting in class I neurons by shearing the spike-onset bifurcation structure. PLoS Comput Biol 2024; 20:e1011751. [PMID: 39133755 PMCID: PMC11383233 DOI: 10.1371/journal.pcbi.1011751] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2023] [Revised: 09/09/2024] [Accepted: 08/02/2024] [Indexed: 09/10/2024] Open
Abstract
Slow brain rhythms, for example during slow-wave sleep or pathological conditions like seizures and spreading depolarization, can be accompanied by oscillations in extracellular potassium concentration. Such slow brain rhythms typically have a lower frequency than tonic action-potential firing. They are assumed to arise from network-level mechanisms, involving synaptic interactions and delays, or from intrinsically bursting neurons. Neuronal burst generation is commonly attributed to ion channels with slow kinetics. Here, we explore an alternative mechanism generically available to all neurons with class I excitability. It is based on the interplay of fast-spiking voltage dynamics with a one-dimensional slow dynamics of the extracellular potassium concentration, mediated by the activity of the Na+/K+-ATPase. We use bifurcation analysis of the complete system as well as the slow-fast method to reveal that this coupling suffices to generate a hysteresis loop organized around a bistable region that emerges from a saddle-node loop bifurcation-a common feature of class I excitable neurons. Depending on the strength of the Na+/K+-ATPase, bursts are generated from pump-induced shearing the bifurcation structure, spiking is tonic, or cells are silenced via depolarization block. We suggest that transitions between these dynamics can result from disturbances in extracellular potassium regulation, such as glial malfunction or hypoxia affecting the Na+/K+-ATPase activity. The identified minimal mechanistic model outlining the sodium-potassium pump's generic contribution to burst dynamics can, therefore, contribute to a better mechanistic understanding of pathologies such as epilepsy syndromes and, potentially, inform therapeutic strategies.
Collapse
Affiliation(s)
- Mahraz Behbood
- Institute for Theoretical Biology, Department of Biology, Humboldt-Universität zu Berlin, Berlin, Germany
- Bernstein Center for Computational Neuroscience, Berlin, Germany
| | - Louisiane Lemaire
- Institute for Theoretical Biology, Department of Biology, Humboldt-Universität zu Berlin, Berlin, Germany
- Bernstein Center for Computational Neuroscience, Berlin, Germany
| | - Jan-Hendrik Schleimer
- Institute for Theoretical Biology, Department of Biology, Humboldt-Universität zu Berlin, Berlin, Germany
- Bernstein Center for Computational Neuroscience, Berlin, Germany
| | - Susanne Schreiber
- Institute for Theoretical Biology, Department of Biology, Humboldt-Universität zu Berlin, Berlin, Germany
- Bernstein Center for Computational Neuroscience, Berlin, Germany
| |
Collapse
|
8
|
Schoknecht K, Maechler M, Wallach I, Dreier JP, Liotta A, Berndt N. Isoflurane lowers the cerebral metabolic rate of oxygen and prevents hypoxia during cortical spreading depolarization in vitro: An integrative experimental and modeling study. J Cereb Blood Flow Metab 2024; 44:1000-1012. [PMID: 38140913 PMCID: PMC11318408 DOI: 10.1177/0271678x231222306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/05/2023] [Revised: 10/23/2023] [Accepted: 11/29/2023] [Indexed: 12/24/2023]
Abstract
Cortical spreading depolarization (SD) imposes a massive increase in energy demand and therefore evolves as a target for treatment following acute brain injuries. Anesthetics are empirically used to reduce energy metabolism in critical brain conditions, yet their effect on metabolism during SD remains largely unknown. We investigated oxidative metabolism during SD in brain slices from Wistar rats. Extracellular potassium ([K+]o), local field potential and partial tissue oxygen pressure (ptiO2) were measured simultaneously. The cerebral metabolic rate of oxygen (CMRO2) was calculated using a reaction-diffusion model. By that, we tested the effect of clinically relevant concentrations of isoflurane on CMRO2 during SD and modeled tissue oxygenation for different capillary pO2 values. During SD, CMRO2 increased 2.7-fold, resulting in transient hypoxia in the slice core. Isoflurane decreased CMRO2, reduced peak [K+]o, and prolonged [K+]o clearance, which indicates reduced synaptic transmission and sodium-potassium ATPase inhibition. Modeling tissue oxygenation during SD illustrates the need for increased capillary pO2 levels to prevent hypoxia. In the absence thereof, isoflurane could improve tissue oxygenation by lowering CMRO2. Therefore, isoflurane is a promising candidate for pre-clinical studies on neuronal survival in conditions involving SD.
Collapse
Affiliation(s)
- Karl Schoknecht
- Carl-Ludwig-Institute of Physiology, Medical Faculty, Leipzig University, Leipzig, Germany
| | - Mathilde Maechler
- Department of Anesthesiology and Intensive Care, Charité – Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany
- Institute of Neurophysiology, Charité – Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany
- Institute of Computer-Assisted Cardiovascular Medicine, Deutsches Herzzentrum der Charité (DHZC), Berlin, Germany
| | - Iwona Wallach
- Institute of Computer-Assisted Cardiovascular Medicine, Deutsches Herzzentrum der Charité (DHZC), Berlin, Germany
- Charité – Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany
| | - Jens P Dreier
- Centre for Stroke Research Berlin, Charité – Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany
- Department of Experimental Neurology, Charité – Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany
- Department of Neurology, Charité – Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany
- Bernstein Centre for Computational Neuroscience Berlin, Berlin, Germany
- Einstein Centre for Neurosciences Berlin, Berlin, Germany
| | - Agustin Liotta
- Department of Anesthesiology and Intensive Care, Charité – Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany
- Institute of Neurophysiology, Charité – Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany
- Institute of Computer-Assisted Cardiovascular Medicine, Deutsches Herzzentrum der Charité (DHZC), Berlin, Germany
- Department of Experimental Neurology, Charité – Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany
- Institute of Health at Charité – Universitätsmedizin Berlin, Berlin
- Neuroscience Research Center, Charité – Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany
| | - Nikolaus Berndt
- Institute of Computer-Assisted Cardiovascular Medicine, Deutsches Herzzentrum der Charité (DHZC), Berlin, Germany
- Charité – Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany
- German Institute of Human Nutrition Potsdam-Rehbruecke (DIfE), Department of Molecular Toxicology, Nuthetal, Germany
| |
Collapse
|
9
|
Verardo C, Mele LJ, Selmi L, Palestri P. Finite-element modeling of neuromodulation via controlled delivery of potassium ions using conductive polymer-coated microelectrodes. J Neural Eng 2024; 21:026002. [PMID: 38306702 DOI: 10.1088/1741-2552/ad2581] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2023] [Accepted: 02/02/2024] [Indexed: 02/04/2024]
Abstract
Objective. The controlled delivery of potassium is an interesting neuromodulation modality, being potassium ions involved in shaping neuron excitability, synaptic transmission, network synchronization, and playing a key role in pathological conditions like epilepsy and spreading depression. Despite many successful examples of pre-clinical devices able to influence the extracellular potassium concentration, computational frameworks capturing the corresponding impact on neuronal activity are still missing.Approach. We present a finite-element model describing a PEDOT:PSS-coated microelectrode (herein, simplyionic actuator) able to release potassium and thus modulate the activity of a cortical neuron in anin-vitro-like setting. The dynamics of ions in the ionic actuator, the neural membrane, and the cellular fluids are solved self-consistently.Main results. We showcase the capability of the model to describe on a physical basis the modulation of the intrinsic excitability of the cell and of the synaptic transmission following the electro-ionic stimulation produced by the actuator. We consider three case studies for the ionic actuator with different levels of selectivity to potassium: ideal selectivity, no selectivity, and selectivity achieved by embedding ionophores in the polymer.Significance. This work is the first step toward a comprehensive computational framework aimed to investigate novel neuromodulation devices targeting specific ionic species, as well as to optimize their design and performance, in terms of the induced modulation of neural activity.
Collapse
Affiliation(s)
- Claudio Verardo
- Polytechnic Department of Engineering and Architecture, Università degli Studi di Udine, Udine, Italy
- BioRobotics Institute and Department of Excellence in Robotics and AI, Scuola Superiore Sant'Anna, Pisa, Italy
| | - Leandro Julian Mele
- Polytechnic Department of Engineering and Architecture, Università degli Studi di Udine, Udine, Italy
- Department of Materials Science and Engineering, Stanford University, Stanford, CA, United States of America
| | - Luca Selmi
- Department of Engineering "Enzo Ferrari", Università degli Studi di Modena e Reggio Emilia, Modena, Italy
| | - Pierpaolo Palestri
- Polytechnic Department of Engineering and Architecture, Università degli Studi di Udine, Udine, Italy
- Department of Engineering "Enzo Ferrari", Università degli Studi di Modena e Reggio Emilia, Modena, Italy
| |
Collapse
|
10
|
Mukherjee S, Mirzaee M, Tithof J. Quantifying the relationship between spreading depolarization and perivascular cerebrospinal fluid flow. Sci Rep 2023; 13:12405. [PMID: 37524734 PMCID: PMC10390554 DOI: 10.1038/s41598-023-38938-5] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2023] [Accepted: 07/17/2023] [Indexed: 08/02/2023] Open
Abstract
Recent studies have linked spreading depolarization (SD, an electro-chemical wave in the brain following stroke, migraine, traumatic brain injury, and more) with increase in cerebrospinal fluid (CSF) flow through the perivascular spaces (PVSs, annular channels lining the brain vasculature). We develop a novel computational model that couples SD and CSF flow. We first use high order numerical simulations to solve a system of physiologically realistic reaction-diffusion equations which govern the spatiotemporal dynamics of ions in the extracellular and intracellular spaces of the brain cortex during SD. We then couple the SD wave with a 1D CSF flow model that captures the change in cross-sectional area, pressure, and volume flow rate through the PVSs. The coupling is modelled using an empirical relationship between the excess potassium ion concentration in the extracellular space following SD and the vessel radius. We find that the CSF volumetric flow rate depends intricately on the length and width of the PVS, as well as the vessel radius and the angle of incidence of the SD wave. We derive analytical expressions for pressure and volumetric flow rates of CSF through the PVS for a given SD wave and quantify CSF flow variations when two SD waves collide. Our numerical approach is very general and could be extended in the future to obtain novel, quantitative insights into how CSF flow in the brain couples with slow waves, functional hyperemia, seizures, or externally applied neural stimulations.
Collapse
Affiliation(s)
- Saikat Mukherjee
- Department of Mechanical Engineering, University of Minnesota, Minneapolis, MN, 55455, USA.
- Department of Mechanical Engineering, Iowa State University, Ames, IA, 50011, USA.
| | - Mahsa Mirzaee
- Department of Mechanical Engineering, University of Minnesota, Minneapolis, MN, 55455, USA
| | - Jeffrey Tithof
- Department of Mechanical Engineering, University of Minnesota, Minneapolis, MN, 55455, USA
| |
Collapse
|
11
|
Øyehaug L. Slow ion concentration oscillations and multiple states in neuron-glia interaction-insights gained from reduced mathematical models. FRONTIERS IN NETWORK PHYSIOLOGY 2023; 3:1189118. [PMID: 37284003 PMCID: PMC10241345 DOI: 10.3389/fnetp.2023.1189118] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/18/2023] [Accepted: 04/28/2023] [Indexed: 06/08/2023]
Abstract
When potassium in the extracellular space separating neurons and glia reaches sufficient levels, neurons may fire spontaneous action potentials or even become inactivated due to membrane depolarisation, which, in turn, may lead to increased extracellular potassium levels. Under certain circumstances, this chain of events may trigger periodic bursts of neuronal activity. In the present study, reduced neuron-glia models are applied to explore the relationship between bursting behaviour and ion concentration dynamics. These reduced models are built based on a previously developed neuron-glia model, in which channel-mediated neuronal sodium and potassium currents are replaced by a function of neuronal sodium and extracellular potassium concentrations. Simulated dynamics of the resulting two reduced models display features that are qualitatively similar to those of the existing neuron-glia model. Bifurcation analyses of the reduced models show rich and interesting dynamics that include the existence of Hopf bifurcations between which the models exhibit slow ion concentration oscillations for a wide range of parameter values. The study demonstrates that even very simple models can provide insights of possible relevance to complex phenomena.
Collapse
|
12
|
Hanssen KØ, Grødem S, Fyhn M, Hafting T, Einevoll GT, Ness TV, Halnes G. Responses in fast-spiking interneuron firing rates to parameter variations associated with degradation of perineuronal nets. J Comput Neurosci 2023; 51:283-298. [PMID: 37058180 PMCID: PMC10182141 DOI: 10.1007/s10827-023-00849-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2022] [Revised: 03/01/2023] [Accepted: 03/10/2023] [Indexed: 04/15/2023]
Abstract
The perineuronal nets (PNNs) are sugar coated protein structures that encapsulate certain neurons in the brain, such as parvalbumin positive (PV) inhibitory neurons. As PNNs are theorized to act as a barrier to ion transport, they may effectively increase the membrane charge-separation distance, thereby affecting the membrane capacitance. Tewari et al. (2018) found that degradation of PNNs induced a 25%-50% increase in membrane capacitance [Formula: see text] and a reduction in the firing rates of PV-cells. In the current work, we explore how changes in [Formula: see text] affects the firing rate in a selection of computational neuron models, ranging in complexity from a single compartment Hodgkin-Huxley model to morphologically detailed PV-neuron models. In all models, an increased [Formula: see text] lead to reduced firing, but the experimentally reported increase in [Formula: see text] was not alone sufficient to explain the experimentally reported reduction in firing rate. We therefore hypothesized that PNN degradation in the experiments affected not only [Formula: see text], but also ionic reversal potentials and ion channel conductances. In simulations, we explored how various model parameters affected the firing rate of the model neurons, and identified which parameter variations in addition to [Formula: see text] that are most likely candidates for explaining the experimentally reported reduction in firing rate.
Collapse
Affiliation(s)
- Kine Ødegård Hanssen
- Department of Physics, University of Oslo, Oslo, Norway.
- Centre for Integrative Neuroplasticity, University of Oslo, Oslo, Norway.
| | - Sverre Grødem
- Centre for Integrative Neuroplasticity, University of Oslo, Oslo, Norway
- Department of Biosciences, University of Oslo, Oslo, Norway
| | - Marianne Fyhn
- Centre for Integrative Neuroplasticity, University of Oslo, Oslo, Norway
- Department of Biosciences, University of Oslo, Oslo, Norway
| | - Torkel Hafting
- Centre for Integrative Neuroplasticity, University of Oslo, Oslo, Norway
- Institute of Basic Medical Sciences, University of Oslo, Oslo, Norway
| | - Gaute T Einevoll
- Department of Physics, University of Oslo, Oslo, Norway
- Centre for Integrative Neuroplasticity, University of Oslo, Oslo, Norway
- Department of Physics, Norwegian University of Life Sciences, Ås, Norway
| | - Torbjørn Vefferstad Ness
- Centre for Integrative Neuroplasticity, University of Oslo, Oslo, Norway
- Department of Physics, Norwegian University of Life Sciences, Ås, Norway
| | - Geir Halnes
- Centre for Integrative Neuroplasticity, University of Oslo, Oslo, Norway
- Department of Physics, Norwegian University of Life Sciences, Ås, Norway
| |
Collapse
|
13
|
Tripathi R, Gluckman BJ. Development of Mechanistic Neural Mass (mNM) Models that Link Physiology to Mean-Field Dynamics. FRONTIERS IN NETWORK PHYSIOLOGY 2022; 2:911090. [PMID: 36876035 PMCID: PMC9980379 DOI: 10.3389/fnetp.2022.911090] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
Abstract
Brain rhythms emerge from the mean-field activity of networks of neurons. There have been many efforts to build mathematical and computational embodiments in the form of discrete cell-group activities-termed neural masses-to understand in particular the origins of evoked potentials, intrinsic patterns of activities such as theta, regulation of sleep, Parkinson's disease related dynamics, and mimic seizure dynamics. As originally utilized, standard neural masses convert input through a sigmoidal function to a firing rate, and firing rate through a synaptic alpha function to other masses. Here we define a process to build mechanistic neural masses (mNMs) as mean-field models of microscopic membrane-type (Hodgkin Huxley type) models of different neuron types that duplicate the stability, firing rate, and associated bifurcations as function of relevant slow variables - such as extracellular potassium - and synaptic current; and whose output is both firing rate and impact on the slow variables - such as transmembrane potassium flux. Small networks composed of just excitatory and inhibitory mNMs demonstrate expected dynamical states including firing, runaway excitation and depolarization block, and these transitions change in biologically observed ways with changes in extracellular potassium and excitatory-inhibitory balance.
Collapse
Affiliation(s)
- Richa Tripathi
- Center for Neural Engineering, The Pennsylvania State University, University Park, PA, United States.,Indian Institute of Technology Gandhinagar, Gandhinagar, India.,Center for Advanced Systems Understanding (CASUS), HZDR, Görlitz, Germany
| | - Bruce J Gluckman
- Center for Neural Engineering, The Pennsylvania State University, University Park, PA, United States.,Departments of Engineering Science and Mechanics, Biomedical Engineering, The Pennsylvania State University, University Park, PA, United States.,Department of Neurosurgery, College of Medicine, The Pennsylvania State University, Hershey, PA, United States
| |
Collapse
|
14
|
Gentiletti D, de Curtis M, Gnatkovsky V, Suffczynski P. Focal seizures are organized by feedback between neural activity and ion concentration changes. eLife 2022; 11:68541. [PMID: 35916367 PMCID: PMC9377802 DOI: 10.7554/elife.68541] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Accepted: 07/12/2022] [Indexed: 11/13/2022] Open
Abstract
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.
Collapse
|
15
|
Neuroelectric Mechanisms of Delayed Cerebral Ischemia after Aneurysmal Subarachnoid Hemorrhage. Int J Mol Sci 2022; 23:ijms23063102. [PMID: 35328523 PMCID: PMC8951073 DOI: 10.3390/ijms23063102] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2022] [Revised: 03/11/2022] [Accepted: 03/11/2022] [Indexed: 12/16/2022] Open
Abstract
Delayed cerebral ischemia (DCI) remains a challenging but very important condition, because DCI is preventable and treatable for improving functional outcomes after aneurysmal subarachnoid hemorrhage (SAH). The pathologies underlying DCI are multifactorial. Classical approaches to DCI focus exclusively on preventing and treating the reduction of blood flow supply. However, recently, glutamate-mediated neuroelectric disruptions, such as excitotoxicity, cortical spreading depolarization and seizures, and epileptiform discharges, have been reported to occur in high frequencies in association with DCI development after SAH. Each of the neuroelectric disruptions can trigger the other, which augments metabolic demand. If increased metabolic demand exceeds the impaired blood supply, the mismatch leads to relative ischemia, resulting in DCI. The neuroelectric disruption also induces inverted vasoconstrictive neurovascular coupling in compromised brain tissues after SAH, causing DCI. Although glutamates and the receptors may play central roles in the development of excitotoxicity, cortical spreading ischemia and epileptic activity-related events, more studies are needed to clarify the pathophysiology and to develop novel therapeutic strategies for preventing or treating neuroelectric disruption-related DCI after SAH. This article reviews the recent advancement in research on neuroelectric disruption after SAH.
Collapse
|
16
|
Ellingsrud AJ, Dukefoss DB, Enger R, Halnes G, Pettersen K, Rognes ME. Validating a Computational Framework for Ionic Electrodiffusion with Cortical Spreading Depression as a Case Study. eNeuro 2022; 9:ENEURO.0408-21.2022. [PMID: 35365505 PMCID: PMC9045477 DOI: 10.1523/eneuro.0408-21.2022] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2021] [Revised: 02/21/2022] [Accepted: 03/12/2022] [Indexed: 11/21/2022] Open
Abstract
Cortical spreading depression (CSD) is a wave of pronounced depolarization of brain tissue accompanied by substantial shifts in ionic concentrations and cellular swelling. Here, we validate a computational framework for modeling electrical potentials, ionic movement, and cellular swelling in brain tissue during CSD. We consider different model variations representing wild-type (WT) or knock-out/knock-down mice and systematically compare the numerical results with reports from a selection of experimental studies. We find that the data for several CSD hallmarks obtained computationally, including wave propagation speed, direct current shift duration, peak in extracellular K+ concentration as well as a pronounced shrinkage of extracellular space (ECS) are well in line with what has previously been observed experimentally. Further, we assess how key model parameters including cellular diffusivity, structural ratios, membrane water and/or K+ permeabilities affect the set of CSD characteristics.
Collapse
Affiliation(s)
- Ada J Ellingsrud
- Department for Numerical Analysis and Scientific Computing, Simula Research Laboratory, Oslo 0164, Norway
| | - Didrik B Dukefoss
- Letten Centre, Division of Anatomy, Department of Molecular Medicine, Institute of Basic Medical Sciences, University of Oslo, Oslo 0317, Norway
| | - Rune Enger
- Letten Centre, Division of Anatomy, Department of Molecular Medicine, Institute of Basic Medical Sciences, University of Oslo, Oslo 0317, Norway
| | - Geir Halnes
- CINPLA, Faculty of Mathematics and Natural Sciences, University of Oslo, Oslo 0316, Norway
- Institute of Physics, Faculty of Science and Technology, Norwegian University of Life Sciences, Ås 1432, Norway
| | - Klas Pettersen
- NORA, Faculty of Mathematics and Natural Sciences, University of Oslo, Oslo 0316, Norway
| | - Marie E Rognes
- Department for Numerical Analysis and Scientific Computing, Simula Research Laboratory, Oslo 0164, Norway
- Department of Mathematics, Faculty of Mathematics and Natural Sciences, University of Bergen, Bergen 5020, Norway
| |
Collapse
|
17
|
Membrane electrical properties of mouse hippocampal CA1 pyramidal neurons during strong inputs. Biophys J 2022; 121:644-657. [PMID: 34999132 PMCID: PMC8873947 DOI: 10.1016/j.bpj.2022.01.002] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2021] [Revised: 09/21/2021] [Accepted: 01/05/2022] [Indexed: 11/24/2022] Open
Abstract
In this work, we highlight an electrophysiological feature often observed in recordings from mouse CA1 pyramidal cells that has so far been ignored by experimentalists and modelers. It consists of a large and dynamic increase in the depolarization baseline (i.e., the minimum value of the membrane potential between successive action potentials during a sustained input) in response to strong somatic current injections. Such an increase can directly affect neurotransmitter release properties and, more generally, the efficacy of synaptic transmission. However, it cannot be explained by any currently available conductance-based computational model. Here we present a model addressing this issue, demonstrating that experimental recordings can be reproduced by assuming that an input current modifies, in a time-dependent manner, the electrical and permeability properties of the neuron membrane by shifting the ionic reversal potentials and channel kinetics. For this reason, we propose that any detailed model of ion channel kinetics for neurons exhibiting this characteristic should be adapted to correctly represent the response and the synaptic integration process during strong and sustained inputs.
Collapse
|
18
|
A unified physiological framework of transitions between seizures, sustained ictal activity and depolarization block at the single neuron level. J Comput Neurosci 2022; 50:33-49. [PMID: 35031915 PMCID: PMC8818009 DOI: 10.1007/s10827-022-00811-1] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2021] [Revised: 11/10/2021] [Accepted: 01/03/2022] [Indexed: 10/29/2022]
Abstract
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.
Collapse
|
19
|
Morris CE, Wheeler JJ, Joos B. The Donnan-dominated resting state of skeletal muscle fibers contributes to resilience and longevity in dystrophic fibers. J Gen Physiol 2022; 154:212743. [PMID: 34731883 PMCID: PMC8570295 DOI: 10.1085/jgp.202112914] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2021] [Accepted: 09/30/2021] [Indexed: 11/28/2022] Open
Abstract
Duchenne muscular dystrophy (DMD) is an X-linked dystrophin-minus muscle-wasting disease. Ion homeostasis in skeletal muscle fibers underperforms as DMD progresses. But though DMD renders these excitable cells intolerant of exertion, sodium overloaded, depolarized, and spontaneously contractile, they can survive for several decades. We show computationally that underpinning this longevity is a strikingly frugal, robust Pump-Leak/Donnan (P-L/D) ion homeostatic process. Unlike neurons, which operate with a costly “Pump-Leak–dominated” ion homeostatic steady state, skeletal muscle fibers operate with a low-cost “Donnan-dominated” ion homeostatic steady state that combines a large chloride permeability with an exceptionally small sodium permeability. Simultaneously, this combination keeps fiber excitability low and minimizes pump expenditures. As mechanically active, long-lived multinucleate cells, skeletal muscle fibers have evolved to handle overexertion, sarcolemmal tears, ischemic bouts, etc.; the frugality of their Donnan dominated steady state lets them maintain the outsized pump reserves that make them resilient during these inevitable transient emergencies. Here, P-L/D model variants challenged with DMD-type insult/injury (low pump-strength, overstimulation, leaky Nav and cation channels) show how chronic “nonosmotic” sodium overload (observed in DMD patients) develops. Profoundly severe DMD ion homeostatic insult/injury causes spontaneous firing (and, consequently, unwanted excitation–contraction coupling) that elicits cytotoxic swelling. Therefore, boosting operational pump-strength and/or diminishing sodium and cation channel leaks should help extend DMD fiber longevity.
Collapse
Affiliation(s)
- Catherine E Morris
- Neuroscience, Ottawa Hospital Research Institute, Ottawa, Canada.,Center for Neural Dynamics, University of Ottawa, Ottawa, Canada
| | | | - Béla Joos
- Center for Neural Dynamics, University of Ottawa, Ottawa, Canada.,Department of Physics, University of Ottawa, Ottawa, Canada
| |
Collapse
|
20
|
Jafarian A, Wykes RC. Impact of DC-Coupled Electrophysiological Recordings for Translational Neuroscience: Case Study of Tracking Neural Dynamics in Rodent Models of Seizures. Front Comput Neurosci 2022; 16:900063. [PMID: 35936824 PMCID: PMC9351053 DOI: 10.3389/fncom.2022.900063] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2022] [Accepted: 06/15/2022] [Indexed: 11/29/2022] Open
Abstract
We propose that to fully understand biological mechanisms underlying pathological brain activity with transitions (e.g., into and out of seizures), wide-bandwidth electrophysiological recordings are important. We demonstrate the importance of ultraslow potential shifts and infraslow oscillations for reliable tracking of synaptic physiology, within a neural mass model, from brain recordings that undergo pathological phase transitions. We use wide-bandwidth data (direct current (DC) to high-frequency activity), recorded using epidural and penetrating graphene micro-transistor arrays in a rodent model of acute seizures. Using this technological approach, we capture the dynamics of infraslow changes that contribute to seizure initiation (active pre-seizure DC shifts) and progression (passive DC shifts). By employing a continuous-discrete unscented Kalman filter, we track biological mechanisms from full-bandwidth data with and without active pre-seizure DC shifts during paroxysmal transitions. We then apply the same methodological approach for tracking the same parameters after application of high-pass-filtering >0.3Hz to both data sets. This approach reveals that ultraslow potential shifts play a fundamental role in the transition to seizure, and the use of high-pass-filtered data results in the loss of key information in regard to seizure onset and termination dynamics.
Collapse
Affiliation(s)
- Amirhossein Jafarian
- Department of Clinical Neurosciences and Cambridge University Hospitals NHS Trust, University of Cambridge, Cambridge, United Kingdom
| | - Rob C Wykes
- Department of Clinical and Experimental Epilepsy, UCL Queen Square Institute of Neurology, London, United Kingdom.,Nanomedicine Lab, University of Manchester, Manchester, United Kingdom
| |
Collapse
|
21
|
Ellingsrud AJ, Boullé N, Farrell PE, Rognes ME. Accurate numerical simulation of electrodiffusion and water movement in brain tissue. MATHEMATICAL MEDICINE AND BIOLOGY-A JOURNAL OF THE IMA 2021; 38:516-551. [PMID: 34791309 DOI: 10.1093/imammb/dqab016] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/12/2021] [Revised: 08/09/2021] [Accepted: 10/19/2021] [Indexed: 12/19/2022]
Abstract
Mathematical modelling of ionic electrodiffusion and water movement is emerging as a powerful avenue of investigation to provide a new physiological insight into brain homeostasis. However, in order to provide solid answers and resolve controversies, the accuracy of the predictions is essential. Ionic electrodiffusion models typically comprise non-trivial systems of non-linear and highly coupled partial and ordinary differential equations that govern phenomena on disparate time scales. Here, we study numerical challenges related to approximating these systems. We consider a homogenized model for electrodiffusion and osmosis in brain tissue and present and evaluate different associated finite element-based splitting schemes in terms of their numerical properties, including accuracy, convergence and computational efficiency for both idealized scenarios and for the physiologically relevant setting of cortical spreading depression (CSD). We find that the schemes display optimal convergence rates in space for problems with smooth manufactured solutions. However, the physiological CSD setting is challenging: we find that the accurate computation of CSD wave characteristics (wave speed and wave width) requires a very fine spatial and fine temporal resolution.
Collapse
Affiliation(s)
| | - Nicolas Boullé
- Mathematical Institute, University of Oxford, Oxford, UK
| | | | | |
Collapse
|
22
|
Amakhin DV, Soboleva EB, Chizhov AV, Zaitsev AV. Insertion of Calcium-Permeable AMPA Receptors during Epileptiform Activity In Vitro Modulates Excitability of Principal Neurons in the Rat Entorhinal Cortex. Int J Mol Sci 2021; 22:12174. [PMID: 34830051 PMCID: PMC8621524 DOI: 10.3390/ijms222212174] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2021] [Revised: 11/05/2021] [Accepted: 11/07/2021] [Indexed: 12/19/2022] Open
Abstract
Epileptic activity leads to rapid insertion of calcium-permeable α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (CP-AMPARs) into the synapses of cortical and hippocampal glutamatergic neurons, which generally do not express them. The physiological significance of this process is not yet fully understood; however, it is usually assumed to be a pathological process that augments epileptic activity. Using whole-cell patch-clamp recordings in rat entorhinal cortex slices, we demonstrate that the timing of epileptiform discharges, induced by 4-aminopyridine and gabazine, is determined by the shunting effect of Ca2+-dependent slow conductance, mediated predominantly by K+-channels. The blockade of CP-AMPARs by IEM-1460 eliminates this extra conductance and consequently increases the rate of discharge generation. The blockade of NMDARs reduced the additional conductance to a lesser extent than the blockade of CP-AMPARs, indicating that CP-AMPARs are a more significant source of intracellular Ca2+. The study's main findings were implemented in a mathematical model, which reproduces the shunting effect of activity-dependent conductance on the generation of discharges. The obtained results suggest that the expression of CP-AMPARs in principal neurons reduces the discharge generation rate and may be considered as a protective mechanism.
Collapse
Affiliation(s)
- Dmitry V. Amakhin
- Sechenov Institute of Evolutionary Physiology and Biochemistry, Toreza Prospekt 44, 194223 Saint Petersburg, Russia; (D.V.A.); (E.B.S.); (A.V.C.)
| | - Elena B. Soboleva
- Sechenov Institute of Evolutionary Physiology and Biochemistry, Toreza Prospekt 44, 194223 Saint Petersburg, Russia; (D.V.A.); (E.B.S.); (A.V.C.)
| | - Anton V. Chizhov
- Sechenov Institute of Evolutionary Physiology and Biochemistry, Toreza Prospekt 44, 194223 Saint Petersburg, Russia; (D.V.A.); (E.B.S.); (A.V.C.)
- Ioffe Institute, Russian Academy of Sciences, Polytekhnicheskaya 26, 194021 Saint Petersburg, Russia
| | - Aleksey V. Zaitsev
- Sechenov Institute of Evolutionary Physiology and Biochemistry, Toreza Prospekt 44, 194223 Saint Petersburg, Russia; (D.V.A.); (E.B.S.); (A.V.C.)
| |
Collapse
|
23
|
Poelzing S, Weinberg SH, Keener JP. Initiation and entrainment of multicellular automaticity via diffusion limited extracellular domains. Biophys J 2021; 120:5279-5294. [PMID: 34757078 DOI: 10.1016/j.bpj.2021.10.034] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2021] [Revised: 09/12/2021] [Accepted: 10/26/2021] [Indexed: 01/07/2023] Open
Abstract
Electrically excitable cells often spontaneously and synchronously depolarize in vitro and in vivo preparations. It remains unclear how cells entrain and autorhythmically activate above the intrinsic mean activation frequency of isolated cells with or without pacemaking mechanisms. Recent studies suggest that cyclic ion accumulation and depletion in diffusion-limited extracellular volumes modulate electrophysiology by ephaptic mechanisms (nongap junction or synaptic coupling). This report explores how potassium accumulation and depletion in a restricted extracellular domain induces spontaneous action potentials in two different computational models of excitable cells without gap junctional coupling: Hodgkin-Huxley and Luo-Rudy. Importantly, neither model will spontaneously activate on its own without external stimuli. Simulations demonstrate that cells sharing a diffusion-limited extracellular compartment can become autorhythmic and entrained despite intercellular electrical heterogeneity. Autorhythmic frequency is modulated by the cleft volume and potassium fluxes through the cleft. Additionally, inexcitable cells can suppress or induce autorhythmic activity in an excitable cell via a shared cleft. Diffusion-limited shared clefts can also entrain repolarization. Critically, this model predicts a mechanism by which diffusion-limited shared clefts can initiate, entrain, and modulate multicellular automaticity in the absence of gap junctions.
Collapse
Affiliation(s)
- Steven Poelzing
- Fralin Biomedical Research Institute at Virginia Tech Carilion, Center for Heart and Reparative Medicine, and the Department of Biomedical Engineering and Mechanics, Virginia Polytechnic Institute and State University, Roanoke, Virginia.
| | - Seth H Weinberg
- Department of Biomedical Engineering, Davis Heart and Lung Research Institute, and the Wexner Medical Center, The Ohio State University, Columbus, Ohio
| | - James P Keener
- Department of Mathematics, University of Utah, Salt Lake City, Utah
| |
Collapse
|
24
|
Lemaire L, Desroches M, Krupa M, Pizzamiglio L, Scalmani P, Mantegazza M. Modeling NaV1.1/SCN1A sodium channel mutations in a microcircuit with realistic ion concentration dynamics suggests differential GABAergic mechanisms leading to hyperexcitability in epilepsy and hemiplegic migraine. PLoS Comput Biol 2021; 17:e1009239. [PMID: 34314446 PMCID: PMC8345895 DOI: 10.1371/journal.pcbi.1009239] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2021] [Revised: 08/06/2021] [Accepted: 07/02/2021] [Indexed: 11/19/2022] Open
Abstract
Loss of function mutations of SCN1A, the gene coding for the voltage-gated sodium channel NaV1.1, cause different types of epilepsy, whereas gain of function mutations cause sporadic and familial hemiplegic migraine type 3 (FHM-3). However, it is not clear yet how these opposite effects can induce paroxysmal pathological activities involving neuronal networks’ hyperexcitability that are specific of epilepsy (seizures) or migraine (cortical spreading depolarization, CSD). To better understand differential mechanisms leading to the initiation of these pathological activities, we used a two-neuron conductance-based model of interconnected GABAergic and pyramidal glutamatergic neurons, in which we incorporated ionic concentration dynamics in both neurons. We modeled FHM-3 mutations by increasing the persistent sodium current in the interneuron and epileptogenic mutations by decreasing the sodium conductance in the interneuron. Therefore, we studied both FHM-3 and epileptogenic mutations within the same framework, modifying only two parameters. In our model, the key effect of gain of function FHM-3 mutations is ion fluxes modification at each action potential (in particular the larger activation of voltage-gated potassium channels induced by the NaV1.1 gain of function), and the resulting CSD-triggering extracellular potassium accumulation, which is not caused only by modifications of firing frequency. Loss of function epileptogenic mutations, on the other hand, increase GABAergic neurons’ susceptibility to depolarization block, without major modifications of firing frequency before it. Our modeling results connect qualitatively to experimental data: potassium accumulation in the case of FHM-3 mutations and facilitated depolarization block of the GABAergic neuron in the case of epileptogenic mutations. Both these effects can lead to pyramidal neuron hyperexcitability, inducing in the migraine condition depolarization block of both the GABAergic and the pyramidal neuron. Overall, our findings suggest different mechanisms of network hyperexcitability for migraine and epileptogenic NaV1.1 mutations, implying that the modifications of firing frequency may not be the only relevant pathological mechanism. The voltage-gated sodium channel NaV1.1 is a major target of human mutations implicated in different pathologies. In particular, mutations identified in certain types of epilepsy cause loss of function of the channel, whereas mutations identified in certain types of migraine (in which spreading depolarizations of the cortical circuits of the brain are involved) cause instead gain of function. Here, we study dysfunctions induced by these differential effects in a two-neuron (GABAergic and pyramidal) conductance-based model with dynamic ion concentrations. We obtain results that can be related to experimental findings in both situations. Namely, extracellular potassium accumulation induced by the activity of the GABAergic neuron in the case of CSD, and higher propensity of the GABAergic neuron to depolarization block in the epileptogenic scenario, without significant modifications of its firing frequency prior to it. Both scenarios can induce hyperexcitability of the pyramidal neuron, leading in the migraine condition to depolarization block of both the GABAergic and the pyramidal neuron. Our results are successfully confronted to experimental data and suggest that modification of firing frequency is not the only key mechanism in these pathologies of neuronal excitability.
Collapse
Affiliation(s)
- Louisiane Lemaire
- Inria Sophia Antipolis Méditerranée Research Centre, MathNeuro Team, Valbonne-Sophia Antipolis, France
- Université Côte d’Azur, Nice, France
- * E-mail: (LL); (MM)
| | - Mathieu Desroches
- Inria Sophia Antipolis Méditerranée Research Centre, MathNeuro Team, Valbonne-Sophia Antipolis, France
- Université Côte d’Azur, Nice, France
| | - Martin Krupa
- Inria Sophia Antipolis Méditerranée Research Centre, MathNeuro Team, Valbonne-Sophia Antipolis, France
- Université Côte d’Azur, Laboratoire Jean-Alexandre Dieudonné, Nice, France
| | - Lara Pizzamiglio
- Université Côte d’Azur, Institute of Molecular and Cellular Pharmacology (IPMC), Valbonne-Sophia Antipolis, France
- CNRS UMR7275, Institute of Molecular and Cellular Pharmacology (IPMC), Valbonne-Sophia Antipolis, France
| | - Paolo Scalmani
- U.O. VII Clinical and Experimental Epileptology, Foundation IRCCS Neurological Institute Carlo Besta, Milan, Italy
| | - Massimo Mantegazza
- Université Côte d’Azur, Institute of Molecular and Cellular Pharmacology (IPMC), Valbonne-Sophia Antipolis, France
- CNRS UMR7275, Institute of Molecular and Cellular Pharmacology (IPMC), Valbonne-Sophia Antipolis, France
- Inserm, Valbonne-Sophia Antipolis, France
- * E-mail: (LL); (MM)
| |
Collapse
|
25
|
Sætra MJ, Einevoll GT, Halnes G. An electrodiffusive neuron-extracellular-glia model for exploring the genesis of slow potentials in the brain. PLoS Comput Biol 2021; 17:e1008143. [PMID: 34270543 PMCID: PMC8318289 DOI: 10.1371/journal.pcbi.1008143] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2020] [Revised: 07/28/2021] [Accepted: 06/28/2021] [Indexed: 11/29/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Marte J. Sætra
- Department of Numerical Analysis and Scientific Computing, Simula Research Laboratory, Oslo, Norway
| | - Gaute T. Einevoll
- Centre for Integrative Neuroplasticity, University of Oslo, Oslo, Norway
- Department of Physics, University of Oslo, Oslo, Norway
- Department of Physics, Norwegian University of Life Sciences, Ås, Norway
| | - Geir Halnes
- Centre for Integrative Neuroplasticity, University of Oslo, Oslo, Norway
- Department of Physics, Norwegian University of Life Sciences, Ås, Norway
| |
Collapse
|
26
|
Kalia M, Meijer HGE, van Gils SA, van Putten MJAM, Rose CR. Ion dynamics at the energy-deprived tripartite synapse. PLoS Comput Biol 2021; 17:e1009019. [PMID: 34143772 PMCID: PMC8244923 DOI: 10.1371/journal.pcbi.1009019] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2020] [Revised: 06/30/2021] [Accepted: 04/28/2021] [Indexed: 01/09/2023] Open
Abstract
The anatomical and functional organization of neurons and astrocytes at 'tripartite synapses' is essential for reliable neurotransmission, which critically depends on ATP. In low energy conditions, synaptic transmission fails, accompanied by a breakdown of ion gradients, changes in membrane potentials and cell swelling. The resulting cellular damage and cell death are causal to the often devastating consequences of an ischemic stroke. The severity of ischemic damage depends on the age and the brain region in which a stroke occurs, but the reasons for this differential vulnerability are far from understood. In the present study, we address this question by developing a comprehensive biophysical model of a glutamatergic synapse to identify key determinants of synaptic failure during energy deprivation. Our model is based on fundamental biophysical principles, includes dynamics of the most relevant ions, i.e., Na+, K+, Ca2+, Cl- and glutamate, and is calibrated with experimental data. It confirms the critical role of the Na+/K+-ATPase in maintaining ion gradients, membrane potentials and cell volumes. Our simulations demonstrate that the system exhibits two stable states, one physiological and one pathological. During energy deprivation, the physiological state may disappear, forcing a transit to the pathological state, which can be reverted when blocking voltage-gated Na+ and K+ channels. Our model predicts that the transition to the pathological state is favoured if the extracellular space fraction is small. A reduction in the extracellular space volume fraction, as, e.g. observed with ageing, will thus promote the brain's susceptibility to ischemic damage. Our work provides new insights into the brain's ability to recover from energy deprivation, with translational relevance for diagnosis and treatment of ischemic strokes.
Collapse
Affiliation(s)
- Manu Kalia
- Applied Analysis, Department of Applied Mathematics, University of Twente, Enschede, The Netherlands
- * E-mail:
| | - Hil G. E. Meijer
- Applied Analysis, Department of Applied Mathematics, University of Twente, Enschede, The Netherlands
| | - Stephan A. van Gils
- Applied Analysis, Department of Applied Mathematics, University of Twente, Enschede, The Netherlands
| | | | - Christine R. Rose
- Institute of Neurobiology, Heinrich Heine University Düsseldorf, Düsseldorf, Germany
| |
Collapse
|
27
|
Contreras SA, Schleimer JH, Gulledge AT, Schreiber S. Activity-mediated accumulation of potassium induces a switch in firing pattern and neuronal excitability type. PLoS Comput Biol 2021; 17:e1008510. [PMID: 34043638 PMCID: PMC8205125 DOI: 10.1371/journal.pcbi.1008510] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2020] [Revised: 06/15/2021] [Accepted: 04/16/2021] [Indexed: 01/30/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Susana Andrea Contreras
- Institute for Theoretical Biology, Humboldt-University of Berlin, Berlin, Germany
- Bernstein Center for Computational Neuroscience Berlin, Berlin, Germany
| | - Jan-Hendrik Schleimer
- Institute for Theoretical Biology, Humboldt-University of Berlin, Berlin, Germany
- Bernstein Center for Computational Neuroscience Berlin, Berlin, Germany
| | - Allan T. Gulledge
- Molecular and Systems Biology, Geisel School of Medicine at Dartmouth College, Hanover, New Hampshire, United States of America
| | - Susanne Schreiber
- Institute for Theoretical Biology, Humboldt-University of Berlin, Berlin, Germany
- Bernstein Center for Computational Neuroscience Berlin, Berlin, Germany
- * E-mail:
| |
Collapse
|
28
|
Liu X, Yang X, Ma Y, Liu J, Shi D, Schipper D. Construction of a nano-rectangular Zn-Nd complex with near-infrared luminescent response towards metal ions. CHINESE CHEM LETT 2021. [DOI: 10.1016/j.cclet.2020.03.016] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
29
|
Sambath K, Liu X, Wan Z, Hutnik L, Belfield KD, Zhang Y. Potassium Ion Fluorescence Probes: Structures, Properties and Bioimaging. CHEMPHOTOCHEM 2020. [DOI: 10.1002/cptc.202000236] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Karthik Sambath
- Department of Chemistry and Environmental Science College of Science and Liberal Arts New Jersey Institute of Technology 323 Martin Luther King Jr. Blvd. Newark NJ 07102 USA
| | - Xiangshan Liu
- Department of Chemistry and Environmental Science College of Science and Liberal Arts New Jersey Institute of Technology 323 Martin Luther King Jr. Blvd. Newark NJ 07102 USA
| | - Zhaoxiong Wan
- Department of Chemistry and Environmental Science College of Science and Liberal Arts New Jersey Institute of Technology 323 Martin Luther King Jr. Blvd. Newark NJ 07102 USA
| | - Lauren Hutnik
- Department of Chemistry and Environmental Science College of Science and Liberal Arts New Jersey Institute of Technology 323 Martin Luther King Jr. Blvd. Newark NJ 07102 USA
| | - Kevin D. Belfield
- Department of Chemistry and Environmental Science College of Science and Liberal Arts New Jersey Institute of Technology 323 Martin Luther King Jr. Blvd. Newark NJ 07102 USA
| | - Yuanwei Zhang
- Department of Chemistry and Environmental Science College of Science and Liberal Arts New Jersey Institute of Technology 323 Martin Luther King Jr. Blvd. Newark NJ 07102 USA
| |
Collapse
|
30
|
An electrodiffusive, ion conserving Pinsky-Rinzel model with homeostatic mechanisms. PLoS Comput Biol 2020; 16:e1007661. [PMID: 32348299 PMCID: PMC7213750 DOI: 10.1371/journal.pcbi.1007661] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2020] [Revised: 05/11/2020] [Accepted: 04/07/2020] [Indexed: 02/05/2023] Open
Abstract
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+, Ca2+, 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.
Collapse
|
31
|
Chizhov AV, Sanin AE. A simple model of epileptic seizure propagation: Potassium diffusion versus axo-dendritic spread. PLoS One 2020; 15:e0230787. [PMID: 32275724 PMCID: PMC7147746 DOI: 10.1371/journal.pone.0230787] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2019] [Accepted: 03/08/2020] [Indexed: 12/29/2022] Open
Abstract
The mechanisms of epileptic discharge generation and spread are not yet fully known. A recently proposed simple biophysical model of interictal and ictal discharges, Epileptor-2, reproduces well the main features of neuronal excitation and ionic dynamics during discharge generation. In order to distinguish between two hypothesized mechanisms of discharge propagation, we extend the model to the case of two-dimensional propagation along the cortical neural tissue. The first mechanism is based on extracellular potassium diffusion, and the second is the propagation of spikes and postsynaptic signals along axons and dendrites. Our simulations show that potassium diffusion is too slow to reproduce an experimentally observed speed of ictal wavefront propagation (tenths of mm/s). By contrast, the synaptic mechanism predicts well the speed and synchronization of the pre-ictal bursts before the ictal front and the afterdischarges in the ictal core. Though this fact diminishes the role of diffusion and electrodiffusion, the model nevertheless highlights the role of potassium extrusion during neuronal excitation, which provides a positive feedback that changes at the ictal wavefront the balance of excitation versus inhibition in favor of excitation. This finding may help to find a target for a treatment to prevent seizure propagation.
Collapse
Affiliation(s)
- Anton V. Chizhov
- Computational Physics Laboratory, Ioffe Institute, Saint Petersburg, Russia
- Laboratory of Molecular Mechanisms of Neural Interactions, Sechenov Institute of Evolutionary Physiology and Biochemistry of the Russian Academy of Sciences, Saint Petersburg, Russia
| | - Aleksei E. Sanin
- Laboratory of Molecular Mechanisms of Neural Interactions, Sechenov Institute of Evolutionary Physiology and Biochemistry of the Russian Academy of Sciences, Saint Petersburg, Russia
- Blue Brain Project, École polytechnique fédérale de Lausanne (EPFL), Campus Biotech, Geneva, Switzerland
| |
Collapse
|
32
|
Ellingsrud AJ, Solbrå A, Einevoll GT, Halnes G, Rognes ME. Finite Element Simulation of Ionic Electrodiffusion in Cellular Geometries. Front Neuroinform 2020; 14:11. [PMID: 32269519 PMCID: PMC7109287 DOI: 10.3389/fninf.2020.00011] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2019] [Accepted: 03/05/2020] [Indexed: 12/31/2022] Open
Abstract
Mathematical models for excitable cells are commonly based on cable theory, which considers a homogenized domain and spatially constant ionic concentrations. Although such models provide valuable insight, the effect of altered ion concentrations or detailed cell morphology on the electrical potentials cannot be captured. In this paper, we discuss an alternative approach to detailed modeling of electrodiffusion in neural tissue. The mathematical model describes the distribution and evolution of ion concentrations in a geometrically-explicit representation of the intra- and extracellular domains. As a combination of the electroneutral Kirchhoff-Nernst-Planck (KNP) model and the Extracellular-Membrane-Intracellular (EMI) framework, we refer to this model as the KNP-EMI model. Here, we introduce and numerically evaluate a new, finite element-based numerical scheme for the KNP-EMI model, capable of efficiently and flexibly handling geometries of arbitrary dimension and arbitrary polynomial degree. Moreover, we compare the electrical potentials predicted by the KNP-EMI and EMI models. Finally, we study ephaptic coupling induced in an unmyelinated axon bundle and demonstrate how the KNP-EMI framework can give new insights in this setting.
Collapse
Affiliation(s)
- Ada J. Ellingsrud
- Department for Scientific Computing and Numerical Analysis, Simula Research Laboratory, Oslo, Norway
| | - Andreas Solbrå
- Centre for Integrative Neuroplasticity, University of Oslo, Oslo, Norway
- Department of Physics, University of Oslo, Oslo, Norway
| | - Gaute T. Einevoll
- Centre for Integrative Neuroplasticity, University of Oslo, Oslo, Norway
- Department of Physics, University of Oslo, Oslo, Norway
- Faculty of Science and Technology, Norwegian University of Life Sciences, Ås, Norway
| | - Geir Halnes
- Centre for Integrative Neuroplasticity, University of Oslo, Oslo, Norway
- Faculty of Science and Technology, Norwegian University of Life Sciences, Ås, Norway
| | - Marie E. Rognes
- Department for Scientific Computing and Numerical Analysis, Simula Research Laboratory, Oslo, Norway
| |
Collapse
|
33
|
El Houssaini K, Bernard C, Jirsa VK. The Epileptor Model: A Systematic Mathematical Analysis Linked to the Dynamics of Seizures, Refractory Status Epilepticus, and Depolarization Block. eNeuro 2020; 7:ENEURO.0485-18.2019. [PMID: 32066612 PMCID: PMC7096539 DOI: 10.1523/eneuro.0485-18.2019] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2018] [Revised: 10/10/2019] [Accepted: 10/28/2019] [Indexed: 11/21/2022] Open
Abstract
One characteristic of epilepsy is the variety of mechanisms leading to the epileptic state, which are still largely unknown. Refractory status epilepticus (RSE) and depolarization block (DB) are other pathological brain activities linked to epilepsy, whose patterns are different and whose mechanisms remain poorly understood. In epileptogenic network modeling, the Epileptor is a generic phenomenological model that has been recently developed to describe the dynamics of seizures. Here, we performed a detailed qualitative analysis of the Epileptor model based on dynamical systems theory and bifurcation analysis, and investigate the dynamic evolution of "normal" activity toward seizures and to the pathological RSE and DB states. The mechanisms of the transition between states are called bifurcations. Our detailed analysis demonstrates that the generic model undergoes different bifurcation types at seizure offset, when varying some selected parameters. We show that the pathological and normal activities can coexist within the same model under some conditions, and demonstrate that there are many pathways leading to and away from these activities. We here archive systematically all behaviors and dynamic regimes of the Epileptor model to serve as a resource in the development of patient-specific brain network models, and more generally in epilepsy research.
Collapse
Affiliation(s)
- Kenza El Houssaini
- Aix Marseille University, INSERM, INS, Institut de Neurosciences des Systèmes, 13005 Marseille, France
| | - Christophe Bernard
- Aix Marseille University, INSERM, INS, Institut de Neurosciences des Systèmes, 13005 Marseille, France
| | - Viktor K Jirsa
- Aix Marseille University, INSERM, INS, Institut de Neurosciences des Systèmes, 13005 Marseille, France
| |
Collapse
|
34
|
Cui MR, Chen LX, Li XL, Xu JJ, Chen HY. NIR Remote-Controlled "Lock-Unlock" Nanosystem for Imaging Potassium Ions in Living Cells. Anal Chem 2020; 92:4558-4565. [PMID: 32066238 DOI: 10.1021/acs.analchem.9b05820] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Despite great achievements in sensitive and selective detection of important biomolecules in living cells, it is still challenging to develop smart and controllable sensing nanodevices for cellular studies that can be activated at desired time in target sites. To address this issue, we have constructed a remote-controlled "lock-unlock" nanosystem for visual analysis of endogenous potassium ions (K+), which employed a dual-stranded aptamer precursor (DSAP) as recognition molecules, SiO2 based gold nanoshells (AuNS) as nanocarriers, and near-infrared ray (NIR) as the remotely applied stimulus. With the well-designed and activatable DSAP-AuNS, the deficiencies of traditional aptamer-based sensors have been successfully overcome, and the undesired response during transport has been avoided, especially in complex physiological microenvironments. While triggered by NIR, the increased local temperature of AuNS induced the dehybridiztion of DSAP, realized the "lock-unlock" switch of the DSAP-AuNS nanosystem, activated the binding capability of aptamer, and then monitored intracellular K+ via the change of fluorescence signal. This DSAP-AuNS nanosystem not only allows us to visualize endogenous ions in living cells at a desired time but also paves the way for fabricating temporal controllable nanodevices for cellular studies.
Collapse
Affiliation(s)
- Mei-Rong Cui
- State Key Laboratory of Analytical Chemistry for Life Science, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P.R. China
| | - Li-Xian Chen
- State Key Laboratory of Analytical Chemistry for Life Science, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P.R. China
| | - Xiang-Ling Li
- State Key Laboratory of Analytical Chemistry for Life Science, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P.R. China.,College of Life Science and Pharmaceutical Engineering, Nanjing Tech University, Nanjing 211816, P.R. China
| | - Jing-Juan Xu
- State Key Laboratory of Analytical Chemistry for Life Science, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P.R. China
| | - Hong-Yuan Chen
- State Key Laboratory of Analytical Chemistry for Life Science, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P.R. China
| |
Collapse
|
35
|
Capo-Rangel G, Gerardo-Giorda L, Somersalo E, Calvetti D. Metabolism plays a central role in the cortical spreading depression: Evidence from a mathematical model. J Theor Biol 2020; 486:110093. [PMID: 31778711 DOI: 10.1016/j.jtbi.2019.110093] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2019] [Revised: 10/28/2019] [Accepted: 11/23/2019] [Indexed: 11/24/2022]
Abstract
The slow propagating waves of strong depolarization of neural cells characterizing cortical spreading depression, or depolarization, (SD) are known to break cerebral homeostasis and induce significant hemodynamic and electro-metabolic alterations. Mathematical models of cortical spreading depression found in the literature tend to focus on the changes occurring at the electrophysiological level rather than on the ensuing metabolic changes. In this paper, we propose a novel mathematical model which is able to simulate the coupled electrophysiology and metabolism dynamics of SD events, including the swelling of neurons and astrocytes and the concomitant shrinkage of extracellular space. The simulations show that the metabolic coupling leads to spontaneous repetitions of the SD events, which the electrophysiological model alone is not capable to produce. The model predictions, which corroborate experimental findings from the literature, show a strong disruption in metabolism accompanying each wave of spreading depression in the form of a sharp decrease of glucose and oxygen concentrations, with a simultaneous increase in lactate concentration which, in turn, delays the clearing of excess potassium in extracellular space. Our model suggests that the depletion of glucose and oxygen concentration is more pronounced in astrocyte than neuron, in line with the partitioning of the energetic cost of potassium clearing. The model suggests that the repeated SD events are electro-metabolic oscillations that cannot be explained by the electrophysiology alone. The model highlights the crucial role of astrocytes in cleaning the excess potassium flooding extracellular space during a spreading depression event: further, if the ratio of glial/neuron density increases, the frequency of cortical SD events decreases, and the peak potassium concentration in extracellular space is lower than with equal volume fractions.
Collapse
Affiliation(s)
| | | | - E Somersalo
- Basque Center for Applied Mathematics, Spain
| | - D Calvetti
- Department of Mathematics, Applied Mathematics and Statistics, Case Western Reserve University, Ohio.
| |
Collapse
|
36
|
Tuttle A, Riera Diaz J, Mori Y. A computational study on the role of glutamate and NMDA receptors on cortical spreading depression using a multidomain electrodiffusion model. PLoS Comput Biol 2019; 15:e1007455. [PMID: 31790388 PMCID: PMC6907880 DOI: 10.1371/journal.pcbi.1007455] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2019] [Revised: 12/12/2019] [Accepted: 10/02/2019] [Indexed: 11/25/2022] Open
Abstract
Cortical spreading depression (SD) is a spreading disruption of ionic homeostasis in the brain during which neurons experience complete and prolonged depolarizations. SD is the basis of migraine aura and is increasingly associated with many other brain pathologies. Here, we study the role of glutamate and NMDA receptor dynamics in the context of an ionic electrodiffusion model. We perform simulations in one (1D) and two (2D) spatial dimension. Our 1D simulations reproduce the "inverted saddle" shape of the extracellular voltage signal for the first time. Our simulations suggest that SD propagation depends on two overlapping mechanisms; one dependent on extracellular glutamate diffusion and NMDA receptors and the other dependent on extracellular potassium diffusion and persistent sodium channel conductance. In 2D simulations, we study the dynamics of spiral waves. We study the properties of the spiral waves in relation to the planar 1D wave, and also compute the energy expenditure associated with the recurrent SD spirals.
Collapse
Affiliation(s)
- Austin Tuttle
- School of Mathematics, University of Minnesota, Minneapolis, Minnesota, United States of America
| | - Jorge Riera Diaz
- Department of Biomedical Engineering, Florida International University, Miami, Florida, United States of America
| | - Yoichiro Mori
- School of Mathematics, University of Minnesota, Minneapolis, Minnesota, United States of America
- Department of Mathematics, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
- Department of Biology, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
| |
Collapse
|
37
|
Abstract
Leukodystrophies are genetically determined disorders affecting the white matter of the central nervous system. The combination of MRI pattern recognition and next-generation sequencing for the definition of novel disease entities has recently demonstrated that many leukodystrophies are due to the primary involvement and/or mutations in genes selectively expressed by cell types other than the oligodendrocytes, the myelin-forming cells in the brain. This has led to a new definition of leukodystrophies as genetic white matter disorders resulting from the involvement of any white matter structural component. As a result, the research has shifted its main focus from oligodendrocytes to other types of neuroglia. Astrocytes are the housekeeping cells of the nervous system, responsible for maintaining homeostasis and normal brain physiology and to orchestrate repair upon injury. Several lines of evidence show that astrocytic interactions with the other white matter cellular constituents play a primary pathophysiologic role in many leukodystrophies. These are thus now classified as astrocytopathies. This chapter addresses how the crosstalk between astrocytes, other glial cells, axons and non-neural cells are essential for the integrity and maintenance of the white matter in health. It also addresses the current knowledge of the cellular pathomechanisms of astrocytic leukodystrophies, and specifically Alexander disease, vanishing white matter, megalencephalic leukoencephalopathy with subcortical cysts and Aicardi-Goutière Syndrome.
Collapse
Affiliation(s)
- M S Jorge
- Department of Pathology, Free University Medical Centre, Amsterdam, The Netherlands
| | - Marianna Bugiani
- Department of Pathology, Free University Medical Centre, Amsterdam, The Netherlands.
| |
Collapse
|
38
|
Toglia P, Ullah G. Mitochondrial dysfunction and role in spreading depolarization and seizure. J Comput Neurosci 2019; 47:91-108. [PMID: 31506806 DOI: 10.1007/s10827-019-00724-6] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2018] [Revised: 03/12/2019] [Accepted: 07/26/2019] [Indexed: 11/24/2022]
Abstract
The effect of pathological phenomena such as epileptic seizures and spreading depolarization (SD) on mitochondria and the potential feedback of mitochondrial dysfunction into the dynamics of those phenomena are complex and difficult to study experimentally due to the simultaneous changes in many variables governing neuronal behavior. By combining a model that accounts for a wide range of neuronal behaviors including seizures, normoxic SD, and hypoxic SD (HSD), together with a detailed model of mitochondrial function and intracellular Ca2+ dynamics, we investigate mitochondrial dysfunction and its potential role in recovery of the neuron from seizures, HSD, and SD. Our results demonstrate that HSD leads to the collapse of mitochondrial membrane potential and cellular ATP levels that recover only when normal oxygen supply is restored. Mitochondrial organic phosphate and pH gradients determine the strength of the depolarization block during HSD and SD, how quickly the cell enters the depolarization block when the oxygen supply is disrupted or potassium in the bath solution is raised beyond the physiological value, and how fast the cell recovers from SD and HSD when normal potassium concentration and oxygen supply are restored. Although not as dramatic as phosphate and pH gradients, mitochondrial Ca2+ uptake has a similar effect on neuronal behavior during these conditions.
Collapse
Affiliation(s)
- Patrick Toglia
- Department of Physics, University of South Florida, 4202 E. Fowler Ave., Tampa, FL, 33620, USA
| | - Ghanim Ullah
- Department of Physics, University of South Florida, 4202 E. Fowler Ave., Tampa, FL, 33620, USA.
| |
Collapse
|
39
|
Kenny A, Plank MJ, David T. The effects of cerebral curvature on cortical spreading depression. J Theor Biol 2019; 472:11-26. [PMID: 30978351 DOI: 10.1016/j.jtbi.2019.04.006] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2019] [Revised: 04/05/2019] [Accepted: 04/08/2019] [Indexed: 10/27/2022]
Abstract
Neuronal activity evokes a localised increase in cerebral blood flow through neurovascular coupling (NVC), a communication system between a group of cells known as a neurovascular unit (NVU). Dysfunctional NVC can lead to pathologies such as cortical spreading depression (CSD), characterised by a slowly propagating wave of neuronal depolarisation and high extracellular potassium (K+) levels. CSD is associated with several neurological disorders such as migraine, stroke, and traumatic brain injury. Insight into the spatial dynamics of CSD in humans is mainly deduced from animal experiments on the smooth lissencephalic brain (in particular murine experiments), however the human cortex is gyrencephalic (highly folded) and is considered likely to exhibit different and more complex patterns of CSD. In this study a large scale numerical NVC model of multiple NVUs is coupled to a vascular tree simulating a two-dimensional cerebral tissue slice. This model is extended with a spatial Gaussian curvature mapping that can simulate the highly folded nature of the human cortex. For a flat surface comparable to a lissencephalic cortex the model can simulate propagating waves of high extracellular K+ travelling radially outwards from a stimulated area at approximately 6.7 mm/min, corresponding well with multiple experimental results. The high K+ concentration induces a corresponding wave of vasoconstriction (with decreased blood flow) then slight vasodilation, achieved through cellular communication within the NVU. The BOLD response decreases below baseline by approximately 10% followed by an increase of 1%. For a surface with spatially varied curvature comparable to a section of gyrencephalic cortex, areas of positive Gaussian curvature inhibit wave propagation due to decreased extracellular diffusion rate. Whereas areas of negative curvature promote propagation. Consequently extracellular K+ is observed travelling as wave segments (as opposed to radial waves) through flat or negatively curved "valleys" corresponding to folds (sulci) in the cortex. If the wave size (defined as the activated area of high K+ concentration) is too small or diffusion rate too low then wave segments can cease propagation. If the diffusion rate is high enough the wave segments can grow from open ends forming loose spiral waves. These results may provide some insight into the differences seen between human and animal experiments.
Collapse
Affiliation(s)
- Allanah Kenny
- Department of Mechanical Engineering, University of Canterbury, New Zealand.
| | - Michael J Plank
- School of Mathematics and Statistics and Te Punaha Matatini, University of Canterbury, New Zealand
| | - Tim David
- Department of Mechanical Engineering, University of Canterbury, New Zealand
| |
Collapse
|
40
|
González OC, Krishnan GP, Timofeev I, Bazhenov M. Ionic and synaptic mechanisms of seizure generation and epileptogenesis. Neurobiol Dis 2019; 130:104485. [PMID: 31150792 DOI: 10.1016/j.nbd.2019.104485] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2019] [Revised: 05/23/2019] [Accepted: 05/27/2019] [Indexed: 01/09/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Oscar C González
- Neurosciences Graduate Program, University of California, San Diego, CA 92093, United States of America; Department of Medicine, University of California, San Diego, CA 92093, United States of America
| | - Giri P Krishnan
- Department of Medicine, University of California, San Diego, CA 92093, United States of America
| | - Igor Timofeev
- Centre de recherche de l'Institut universitaire en santé mentale de Québec (CRIUSMQ), 2601 de la Canardière, Québec, QC, Canada; Department of Psychiatry and Neuroscience, Université Laval, Québec, QC, Canada
| | - Maxim Bazhenov
- Neurosciences Graduate Program, University of California, San Diego, CA 92093, United States of America; Department of Medicine, University of California, San Diego, CA 92093, United States of America.
| |
Collapse
|
41
|
Chizhov AV, Amakhin DV, Zaitsev AV. Mathematical model of Na-K-Cl homeostasis in ictal and interictal discharges. PLoS One 2019; 14:e0213904. [PMID: 30875397 PMCID: PMC6420042 DOI: 10.1371/journal.pone.0213904] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2018] [Accepted: 03/04/2019] [Indexed: 12/22/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Anton V. Chizhov
- Computational Physics Laboratory, Ioffe Institute, Saint Petersburg, Russia
- Laboratory of Molecular Mechanisms of Neural Interactions, Sechenov Institute of Evolutionary Physiology and Biochemistry of the Russian Academy of Sciences, Saint Petersburg, Russia
- * E-mail:
| | - Dmitry V. Amakhin
- Laboratory of Molecular Mechanisms of Neural Interactions, Sechenov Institute of Evolutionary Physiology and Biochemistry of the Russian Academy of Sciences, Saint Petersburg, Russia
| | - Aleksey V. Zaitsev
- Laboratory of Molecular Mechanisms of Neural Interactions, Sechenov Institute of Evolutionary Physiology and Biochemistry of the Russian Academy of Sciences, Saint Petersburg, Russia
| |
Collapse
|
42
|
Hart JL, Gremaud PA, David T. Global Sensitivity Analysis of High-Dimensional Neuroscience Models: An Example of Neurovascular Coupling. Bull Math Biol 2019; 81:1805-1828. [PMID: 30820832 DOI: 10.1007/s11538-019-00578-0] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2018] [Accepted: 02/05/2019] [Indexed: 11/29/2022]
Abstract
The complexity and size of state-of-the-art cell models have significantly increased in part due to the requirement that these models possess complex cellular functions which are thought-but not necessarily proven-to be important. Modern cell models often involve hundreds of parameters; the values of these parameters come, more often than not, from animal experiments whose relationship to the human physiology is weak with very little information on the errors in these measurements. The concomitant uncertainties in parameter values result in uncertainties in the model outputs or quantities of interest (QoIs). Global sensitivity analysis (GSA) aims at apportioning to individual parameters (or sets of parameters) their relative contribution to output uncertainty thereby introducing a measure of influence or importance of said parameters. New GSA approaches are required to deal with increased model size and complexity; a three-stage methodology consisting of screening (dimension reduction), surrogate modeling, and computing Sobol' indices, is presented. The methodology is used to analyze a physiologically validated numerical model of neurovascular coupling which possess 160 uncertain parameters. The sensitivity analysis investigates three quantities of interest, the average value of [Formula: see text] in the extracellular space, the average volumetric flow rate through the perfusing vessel, and the minimum value of the actin/myosin complex in the smooth muscle cell. GSA provides a measure of the influence of each parameter, for each of the three QoIs, giving insight into areas of possible physiological dysfunction and areas of further investigation.
Collapse
Affiliation(s)
- J L Hart
- Department of Mathematics, North Carolina State University, Raleigh, NC, USA
| | - P A Gremaud
- Department of Mathematics, North Carolina State University, Raleigh, NC, USA
| | - T David
- Department of Mechanical Engineering, University of Canterbury, Christchurch, New Zealand.
| |
Collapse
|
43
|
Liu J, Yang X, Zhu T, Yuan B, Chen H, Shi D, Schipper D, Jones RA. Construction of a crystalline 14-metal Zn–Nd rectangular nanocluster with a dual-emissive response towards metal ions. RSC Adv 2019; 9:40017-40022. [PMID: 35541392 PMCID: PMC9076200 DOI: 10.1039/c9ra09409a] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2019] [Accepted: 11/18/2019] [Indexed: 12/01/2022] Open
Abstract
A crystalline 14-metal Zn–Nd cluster [Zn6Nd8L2(OAc)20(O)2(NO3)4(OC2H5)4] (1) was constructed using a flexible Schiff base ligand with long-chain (CH2)2O(CH2)2O(CH2)2 backbone. 1 exhibits a nanoscale rectangular structure (10 × 14 × 18 Å). Besides the visible ligand-centered emission, 1 displays typical NIR luminescence of Nd3+. Interestingly, 1 shows both a ligand-centered and lanthanide fluorescence response towards metal ions, especially to K+ and Co2+ at the ppm level. A crystalline 14-metal Zn–Nd nanocluster was constructed, and it shows ligand-centered and lanthanide fluorescence response towards metal ions. ![]()
Collapse
Affiliation(s)
- Jieni Liu
- Zhejiang Key Laboratory of Carbon Materials
- College of Chemistry and Materials Engineering
- Wenzhou University
- Wenzhou 325035
- China
| | - Xiaoping Yang
- Zhejiang Key Laboratory of Carbon Materials
- College of Chemistry and Materials Engineering
- Wenzhou University
- Wenzhou 325035
- China
| | - Ting Zhu
- Zhejiang Key Laboratory of Carbon Materials
- College of Chemistry and Materials Engineering
- Wenzhou University
- Wenzhou 325035
- China
| | - Bichen Yuan
- Zhejiang Key Laboratory of Carbon Materials
- College of Chemistry and Materials Engineering
- Wenzhou University
- Wenzhou 325035
- China
| | - Hongfen Chen
- Zhejiang Key Laboratory of Carbon Materials
- College of Chemistry and Materials Engineering
- Wenzhou University
- Wenzhou 325035
- China
| | - Dongliang Shi
- Zhejiang Key Laboratory of Carbon Materials
- College of Chemistry and Materials Engineering
- Wenzhou University
- Wenzhou 325035
- China
| | - Desmond Schipper
- Department of Chemistry and Biochemistry
- The University of Texas at Austin
- Austin
- USA
| | - Richard A. Jones
- Department of Chemistry and Biochemistry
- The University of Texas at Austin
- Austin
- USA
| |
Collapse
|
44
|
Macro scale modelling of cortical spreading depression and the role of astrocytic gap junctions. J Theor Biol 2018; 458:78-91. [DOI: 10.1016/j.jtbi.2018.09.006] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2018] [Revised: 08/08/2018] [Accepted: 09/07/2018] [Indexed: 12/22/2022]
|
45
|
Düsterwald KM, Currin CB, Burman RJ, Akerman CJ, Kay AR, Raimondo JV. Biophysical models reveal the relative importance of transporter proteins and impermeant anions in chloride homeostasis. eLife 2018; 7:39575. [PMID: 30260315 PMCID: PMC6200395 DOI: 10.7554/elife.39575] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2018] [Accepted: 09/24/2018] [Indexed: 11/17/2022] Open
Abstract
Fast synaptic inhibition in the nervous system depends on the transmembrane flux of Cl- ions based on the neuronal Cl- driving force. Established theories regarding the determinants of Cl- driving force have recently been questioned. Here, we present biophysical models of Cl- homeostasis using the pump-leak model. Using numerical and novel analytic solutions, we demonstrate that the Na+/K+-ATPase, ion conductances, impermeant anions, electrodiffusion, water fluxes and cation-chloride cotransporters (CCCs) play roles in setting the Cl- driving force. Our models, together with experimental validation, show that while impermeant anions can contribute to setting [Cl-]i in neurons, they have a negligible effect on the driving force for Cl- locally and cell-wide. In contrast, we demonstrate that CCCs are well-suited for modulating Cl- driving force and hence inhibitory signaling in neurons. Our findings reconcile recent experimental findings and provide a framework for understanding the interplay of different chloride regulatory processes in neurons. Cells called neurons in the brain communicate by triggering or inhibiting electrical activity in other neurons. To inhibit electrical activity, a signal from one neuron usually triggers specific receptors on the second neuron to open, which allows particles called chloride ions to flow into or out of the neuron. The force that moves chloride ions (the so-called ‘chloride driving force’) depends on two main factors. Firstly, chloride ions, like other particles, tend to move from an area where they are plentiful to areas where they are less abundant. Secondly, chloride ions are negatively charged and are therefore attracted to areas where the net charge (determined by the mix of positively and negatively charged particles) is more positive than their current position. It was previously believed that a group of proteins known as CCCs, which transport chloride ions and positive ions together across the membranes surrounding cells, sets the chloride driving force. However, it has recently been suggested that negatively charged ions that are unable to cross the membrane (or ‘impermeant anions’ for short) may set the driving force instead by contributing to the net charge across the membrane. Düsterwald et al. used a computational model of the neuron to explore these two possibilities. In the simulations, altering the activity of the CCCs led to big changes in the chloride driving force. Changing the levels of impermeant anions altered the volume of cells, but did not drive changes in the chloride driving force. This was because the flow of chloride ions across the membrane led to a compensatory change in the net charge across the membrane. Düsterwald et al. then used an experimental technique called patch-clamping in mice and rats to confirm the model’s predictions. Defects in controlling the chloride driving force in brain cells have been linked with epilepsy, stroke and other neurological diseases. Therefore, a better knowledge of these mechanisms may in future help to identify the best targets for drugs to treat such conditions.
Collapse
Affiliation(s)
- Kira M Düsterwald
- Division of Cell Biology, Department of Human Biology, Faculty of Health Sciences, University of Cape Town, Cape Town, South Africa.,Neuroscience Institute, University of Cape Town, Cape Town, South Africa.,Institute of Infectious Disease and Molecular Medicine, University of Cape Town, Cape Town, South Africa
| | - Christopher B Currin
- Division of Cell Biology, Department of Human Biology, Faculty of Health Sciences, University of Cape Town, Cape Town, South Africa.,Neuroscience Institute, University of Cape Town, Cape Town, South Africa.,Institute of Infectious Disease and Molecular Medicine, University of Cape Town, Cape Town, South Africa
| | - Richard J Burman
- Division of Cell Biology, Department of Human Biology, Faculty of Health Sciences, University of Cape Town, Cape Town, South Africa.,Neuroscience Institute, University of Cape Town, Cape Town, South Africa.,Institute of Infectious Disease and Molecular Medicine, University of Cape Town, Cape Town, South Africa
| | - Colin J Akerman
- Department of Pharmacology, University of Oxford, Oxford, United Kingdom
| | - Alan R Kay
- Department of Biology, University of Iowa, Iowa City Iowa, United States
| | - Joseph V Raimondo
- Division of Cell Biology, Department of Human Biology, Faculty of Health Sciences, University of Cape Town, Cape Town, South Africa.,Neuroscience Institute, University of Cape Town, Cape Town, South Africa.,Institute of Infectious Disease and Molecular Medicine, University of Cape Town, Cape Town, South Africa
| |
Collapse
|
46
|
Buchin A, Kerr CC, Huberfeld G, Miles R, Gutkin B. Adaptation and Inhibition Control Pathological Synchronization in a Model of Focal Epileptic Seizure. eNeuro 2018; 5:ENEURO.0019-18.2018. [PMID: 30302390 PMCID: PMC6173584 DOI: 10.1523/eneuro.0019-18.2018] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2018] [Revised: 06/07/2018] [Accepted: 06/07/2018] [Indexed: 01/12/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Anatoly Buchin
- University of Washington, Department of Physiology and Biophysics (United States, Seattle), 1959 NE Pacific St, 98195
| | - Cliff C. Kerr
- University of Sydney, School of Physics (Australia, Sydney), Physics Rd, NSW 2006
| | - Gilles Huberfeld
- Sorbonne Université-UPMC, Pitié-Salpêtrière Hô, Neurophysiology Department (France, Paris), 47-83 Boulevard de l’Hôpital, 75013
- Institut national de la santé et de la recherche médicale Unit 1129 “Infantile Epilepsies and Brain Plasticity”, Paris Descartes University, Sorbonne Paris Cité University group, (France, Paris), 149 rue de Sévres 75015
| | - Richard Miles
- Brain and Spine Institute, Cortex and Epilepsie Group (France, Paris), 47 Boulevard Hôpital, 75013
| | - Boris Gutkin
- Paris Sciences & Lettres Research University, Laboratoire des Neurosciences Cognitives, Group for Neural Theory (France, Paris), 29, rue d'Ulm, 75005 France
- National Research University Higher School of Economics, Center for Cognition and Decision Making (Russia, Moscow), 20 Myasnitskaya, 109316
| |
Collapse
|
47
|
Yao C, He Z, Nakano T, Shuai J. Spiking patterns of a neuron model to stimulus: Rich dynamics and oxygen's role. CHAOS (WOODBURY, N.Y.) 2018; 28:083112. [PMID: 30180647 DOI: 10.1063/1.5018707] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/08/2017] [Accepted: 07/23/2018] [Indexed: 06/08/2023]
Abstract
Neuronal spiking patterns, which are of fundamental importance for the understanding of information processing in neural systems, can be generated in response to different stimuli. We here investigate in detail the stimulus-induced spiking patterns in a biologically plausible neuron model in which the oxygen concentration and the dynamical concentrations of potassium, sodium, and chloride are considered. Various types of spiking patterns can be induced by the different external potassium accumulations in response to the stimulus, including two different types of epileptic seizure (SZ) and spreading depression (SD) states, two different mixed states of SD and SZ, SZ state with multi-burst, and tonic firing behaviors. Interestingly, we show that these rich spiking patterns can also be induced by the current stimulus with a low oxygen concentration. Furthermore, we reveal that the stimulus can induce two different phase transitions from the SD state to the SZ state according to the phase transition theory, which results in the different electrical activities. All these findings may provide insight into information processing in neural systems.
Collapse
Affiliation(s)
- Chenggui Yao
- Department of Mathematics, Shaoxing University, Shaoxing 312000, China
| | - Zhiwei He
- Department of Mathematics, Shaoxing University, Shaoxing 312000, China
| | - Tadashi Nakano
- Graduate School of Frontier Biosciences, Osaka University, 1-1 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Jianwei Shuai
- Department of Physics, Xiamen University, Xiamen 361005, People's Republic of China
| |
Collapse
|
48
|
Newton AJH, McDougal RA, Hines ML, Lytton WW. Using NEURON for Reaction-Diffusion Modeling of Extracellular Dynamics. Front Neuroinform 2018; 12:41. [PMID: 30042670 PMCID: PMC6049079 DOI: 10.3389/fninf.2018.00041] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2018] [Accepted: 06/12/2018] [Indexed: 11/13/2022] Open
Abstract
Development of credible clinically-relevant brain simulations has been slowed due to a focus on electrophysiology in computational neuroscience, neglecting the multiscale whole-tissue modeling approach used for simulation in most other organ systems. We have now begun to extend the NEURON simulation platform in this direction by adding extracellular modeling. The extracellular medium of neural tissue is an active medium of neuromodulators, ions, inflammatory cells, oxygen, NO and other gases, with additional physiological, pharmacological and pathological agents. These extracellular agents influence, and are influenced by, cellular electrophysiology, and cellular chemophysiology-the complex internal cellular milieu of second-messenger signaling and cascades. NEURON's extracellular reaction-diffusion is supported by an intuitive Python-based where/who/what command sequence, derived from that used for intracellular reaction diffusion, to support coarse-grained macroscopic extracellular models. This simulation specification separates the expression of the conceptual model and parameters from the underlying numerical methods. In the volume-averaging approach used, the macroscopic model of tissue is characterized by free volume fraction-the proportion of space in which species are able to diffuse, and tortuosity-the average increase in path length due to obstacles. These tissue characteristics can be defined within particular spatial regions, enabling the modeler to account for regional differences, due either to intrinsic organization, particularly gray vs. white matter, or to pathology such as edema. We illustrate simulation development using spreading depression, a pathological phenomenon thought to play roles in migraine, epilepsy and stroke. Simulation results were verified against analytic results and against the extracellular portion of the simulation run under FiPy. The creation of this NEURON interface provides a pathway for interoperability that can be used to automatically export this class of models into complex intracellular/extracellular simulations and future cross-simulator standardization.
Collapse
Affiliation(s)
- Adam J. H. Newton
- Department of Neuroscience, Yale University, New Haven, CT, United States
- SUNY Downstate Medical Center, The State University of New York, New York, NY, United States
| | - Robert A. McDougal
- Department of Neuroscience, Yale University, New Haven, CT, United States
- Center for Medical Informatics, Yale University, New Haven, CT, United States
| | - Michael L. Hines
- Department of Neuroscience, Yale University, New Haven, CT, United States
| | - William W. Lytton
- SUNY Downstate Medical Center, The State University of New York, New York, NY, United States
- Neurology, Kings County Hospital Center, Brooklyn, NY, United States
| |
Collapse
|
49
|
Integrated models of neurovascular coupling and BOLD signals: Responses for varying neural activations. Neuroimage 2018. [DOI: 10.1016/j.neuroimage.2018.03.010] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
|
50
|
Chizhov AV, Zefirov AV, Amakhin DV, Smirnova EY, Zaitsev AV. Minimal model of interictal and ictal discharges "Epileptor-2". PLoS Comput Biol 2018; 14:e1006186. [PMID: 29851959 PMCID: PMC6005638 DOI: 10.1371/journal.pcbi.1006186] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2017] [Revised: 06/18/2018] [Accepted: 05/09/2018] [Indexed: 12/01/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Anton V. Chizhov
- Laboratory of Molecular Mechanisms of Neural Interactions, Sechenov Institute of Evolutionary Physiology and Biochemistry of the Russian Academy of Sciences, Saint Petersburg, Russia
- Computational Physics Laboratory, Ioffe Institute, Saint Petersburg, Russia
- * E-mail:
| | - Artyom V. Zefirov
- Laboratory of Molecular Mechanisms of Neural Interactions, Sechenov Institute of Evolutionary Physiology and Biochemistry of the Russian Academy of Sciences, Saint Petersburg, Russia
- Computational Physics Laboratory, Ioffe Institute, Saint Petersburg, Russia
| | - Dmitry V. Amakhin
- Laboratory of Molecular Mechanisms of Neural Interactions, Sechenov Institute of Evolutionary Physiology and Biochemistry of the Russian Academy of Sciences, Saint Petersburg, Russia
| | - Elena Yu. Smirnova
- Laboratory of Molecular Mechanisms of Neural Interactions, Sechenov Institute of Evolutionary Physiology and Biochemistry of the Russian Academy of Sciences, Saint Petersburg, Russia
- Computational Physics Laboratory, Ioffe Institute, Saint Petersburg, Russia
| | - Aleksey V. Zaitsev
- Laboratory of Molecular Mechanisms of Neural Interactions, Sechenov Institute of Evolutionary Physiology and Biochemistry of the Russian Academy of Sciences, Saint Petersburg, Russia
- Institute of Experimental Medicine, Almazov National Medical Research Centre, Saint Petersburg, Russia
| |
Collapse
|