A Continuum Neuronal Model for the Instigation and Propagation of Cortical Spreading Depression

Quantifying the relationship between spreading depolarization and perivascular cerebrospinal fluid flow

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.

Analysis of Network Models with Neuron-Astrocyte Interactions

Neural networks, composed of many neurons and governed by complex interactions between them, are a widely accepted formalism for modeling and exploring global dynamics and emergent properties in brain systems. In the past decades, experimental evidence of computationally relevant neuron-astrocyte interactions, as well as the astrocytic modulation of global neural dynamics, have accumulated. These findings motivated advances in computational glioscience and inspired several models integrating mechanisms of neuron-astrocyte interactions into the standard neural network formalism. These models were developed to study, for example, synchronization, information transfer, synaptic plasticity, and hyperexcitability, as well as classification tasks and hardware implementations. We here focus on network models of at least two neurons interacting bidirectionally with at least two astrocytes that include explicitly modeled astrocytic calcium dynamics. In this study, we analyze the evolution of these models and the biophysical, biochemical, cellular, and network mechanisms used to construct them. Based on our analysis, we propose how to systematically describe and categorize interaction schemes between cells in neuron-astrocyte networks. We additionally study the models in view of the existing experimental data and present future perspectives. Our analysis is an important first step towards understanding astrocytic contribution to brain functions. However, more advances are needed to collect comprehensive data about astrocyte morphology and physiology in vivo and to better integrate them in data-driven computational models. Broadening the discussion about theoretical approaches and expanding the computational tools is necessary to better understand astrocytes' roles in brain functions.

Accurate numerical simulation of electrodiffusion and water movement in brain tissue

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.

A tridomain model for potassium clearance in optic nerve of Necturus

Complex fluids flow in complex ways in complex structures. Transport of water and various organic and inorganic molecules in the central nervous system are important in a wide range of biological and medical processes. However, the exact driving mechanisms are often not known. In this work, we investigate flows induced by action potentials in an optic nerve as a prototype of the central nervous system. Different from traditional fluid dynamics problems, flows in biological tissues such as the central nervous system are coupled with ion transport. They are driven by osmosis created by concentration gradient of ionic solutions, which in turn influence the transport of ions. Our mathematical model is based on the known structural and biophysical properties of the experimental system used by the Harvard group Orkand et al. Asymptotic analysis and numerical computation show the significant role of water in convective ion transport. The full model (including water) and the electrodiffusion model (excluding water) are compared in detail to reveal an interesting interplay between water and ion transport. In the full model, convection due to water flow dominates inside the glial domain. This water flow in the glia contributes significantly to the spatial buffering of potassium in the extracellular space. Convection in the extracellular domain does not contribute significantly to spatial buffering. Electrodiffusion is the dominant mechanism for flows confined to the extracellular domain.

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

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

A computational study on the role of glutamate and NMDA receptors on cortical spreading depression using a multidomain electrodiffusion model

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.

An Algorithm for Automated, Noninvasive Detection of Cortical Spreading Depolarizations Based on EEG Simulations

OBJECTIVE

We present a novel signal processing algorithm for automated, noninvasive detection of cortical spreading depolarizations (CSDs) using electroencephalography (EEG) signals and validate the algorithm on simulated EEG signals. CSDs are waves of neurochemical changes that suppress the neuronal activity as they propagate across the brain's cortical surface. CSDs are believed to mediate secondary brain damage after brain trauma and cerebrovascular diseases like stroke. We address the following two key challenges in detecting CSDs from EEG signals: i) attenuation and loss of high spatial resolution information; and ii) cortical folds, which complicate tracking CSD waves.

METHODS

Our algorithm detects and tracks "wavefronts" of a CSD wave, and stitch together data across space and time to make a detection. To test our algorithm, we provide different models of CSD waves, including different widths of CSD suppressions and different patterns, and use them to simulate scalp EEG signals using head models of four subjects.

RESULTS AND CONCLUSION

Our results suggest that low-density EEG grids (40 electrodes) can detect CSD widths of 1.1 cm on average, while higher density EEG grids (340 electrodes) can detect CSD patterns as thin as 0.43 cm (less than minimum widths reported in prior works), among which single-gyrus CSDs are the hardest to detect because of their small suppression area.

SIGNIFICANCE

The proposed algorithm is a first step toward noninvasive, automated detection of CSDs, which can help in reducing secondary brain damages.

A Bidomain Model for Lens Microcirculation

There exists a large body of research on the lens of the mammalian eye over the past several decades. The objective of this work is to provide a link between the most recent computational models and some of the pioneering work in the 1970s and 80s. We introduce a general nonelectroneutral model to study the microcirculation in the lens of the eye. It describes the steady-state relationships among ion fluxes, between water flow and electric field inside cells, and in the narrow extracellular spaces between cells in the lens. Using asymptotic analysis, we derive a simplified model based on physiological data and compare our results with those in the literature. We show that our simplified model can be reduced further to the first-generation models, whereas our full model is consistent with the most recent computational models. In addition, our simplified model captures in its equations the main features of the full computational models. Our results serve as a useful link intermediate between the computational models and the first-generation analytical models. Simplified models of this sort may be particularly helpful as the roles of similar osmotic pumps of microcirculation are examined in other tissues with narrow extracellular spaces, such as cardiac and skeletal muscle, liver, kidney, epithelia in general, and the narrow extracellular spaces of the central nervous system, the "brain." Simplified models may reveal the general functional plan of these systems before full computational models become feasible and specific.

Noise-sustained patterns in a model of volume-coupled neural tissue

Computational neuroscience operates on models based on several important paradigms. Among them is the assumption that coupling in neural ensembles is provided by chemical or electrical synapses. This assumption works well under normal conditions. However, there is a growing body of data that show the importance of other communication pathways caused by bi-directional transport of substances between the cells and the intercellular space. This type of interaction is called "volume transmission" and has not been rarely addressed in the model studies. The volume transmission pathway naturally appears in multidimensional quantitative models of cellular processes, but is not sufficiently represented at the level of lumped and computationally effective neural models. In this paper, we propose a simple model that allows one to study the features of volume transmission coupling at various spatial scales and taking into account various inhomogeneities. This model is obtained by the extension of the well-known FitzHugh-Nagumo system by the addition of the nonlinear terms and equations to describe, at a qualitative level, the release of potassium into the intercellular space, its diffusion, and the reverse effect on the neurons. The study of model dynamics in various spatial configurations has revealed a number of characteristic spatio-temporal types of behavior that include self-organizing bursting and phase-locked firing patterns, different scenarios of excitation spreading, noise-sustained target patterns, and long-living slow moving wave segments.

Turing-like structures in a functional model of cortical spreading depression

Cortical spreading depression (CSD) along with migraine waves and spreading depolarization events with stroke or injures are the front-line examples of extreme physiological behaviors of the brain cortex which manifest themselves via the onset and spreading of localized areas of neuronal hyperactivity followed by their depression. While much is known about the physiological pathways involved, the dynamical mechanisms of the formation and evolution of complex spatiotemporal patterns during CSD are still poorly understood, in spite of the number of modeling studies that have been already performed. Recently we have proposed a relatively simple mathematical model of cortical spreading depression which counts the effects of neurovascular coupling and cerebral blood flow redistribution during CSD. In the present study, we address the main dynamical consequences of newly included pathways, namely, the changes in the formation and propagation speed of the CSD front and the pattern formation features in two dimensions. Our most notable finding is that the combination of vascular-mediated spatial coupling with local regulatory mechanisms results in the formation of stationary Turing-like patterns during a CSD event.

A mathematical model of recurrent spreading depolarizations

A detailed biophysical model for a neuron/astrocyte network is developed in order to explore mechanisms responsible for the initiation and propagation of recurrent cortical spreading depolarizations. The model incorporates biophysical processes not considered in the earlier models. This includes a model for the Na⁺-glutamate transporter, which allows for a detailed description of reverse glutamate uptake. In particular, we consider the specific roles of elevated extracellular glutamate and K⁺ in the initiation, propagation and recurrence of spreading depolarizations.

Ionic mechanisms underlying history-dependence of conduction delay in an unmyelinated axon

Axonal conduction velocity can change substantially during ongoing activity, thus modifying spike interval structures and, potentially, temporal coding. We used a biophysical model to unmask mechanisms underlying the history-dependence of conduction. The model replicates activity in the unmyelinated axon of the crustacean stomatogastric pyloric dilator neuron. At the timescale of a single burst, conduction delay has a non-monotonic relationship with instantaneous frequency, which depends on the gating rates of the fast voltage-gated Na⁺ current. At the slower timescale of minutes, the mean value and variability of conduction delay increase. These effects are because of hyperpolarization of the baseline membrane potential by the Na⁺/K⁺ pump, balanced by an h-current, both of which affect the gating of the Na⁺ current. We explore the mechanisms of history-dependence of conduction delay in axons and develop an empirical equation that accurately predicts this history-dependence, both in the model and in experimental measurements.

Neuroprotective Role of Gap Junctions in a Neuron Astrocyte Network Model

A detailed biophysical model for a neuron/astrocyte network is developed to explore mechanisms responsible for the initiation and propagation of cortical spreading depolarizations and the role of astrocytes in maintaining ion homeostasis, thereby preventing these pathological waves. Simulations of the model illustrate how properties of spreading depolarizations, such as wave speed and duration of depolarization, depend on several factors, including the neuron and astrocyte Na(+)-K(+) ATPase pump strengths. In particular, we consider the neuroprotective role of astrocyte gap junction coupling. The model demonstrates that a syncytium of electrically coupled astrocytes can maintain a physiological membrane potential in the presence of an elevated extracellular K(+) concentration and efficiently distribute the excess K(+) across the syncytium. This provides an effective neuroprotective mechanism for delaying or preventing the initiation of spreading depolarizations.

Effects of Glia in a Triphasic Continuum Model of Cortical Spreading Depression

Cortical spreading depression (SD) is a spreading disruption in brain ionic homeostasis during which neurons experience complete and prolonged depolarizations. SD is generally believed to be the physiological substrate of migraine aura and is associated with many other brain pathologies. Here, we perform simulations with a model of SD treating brain tissue as a triphasic continuum of neurons, glia and the extracellular space. A thermodynamically consistent incorporation of the major biophysical effects, including ionic electrodiffusion and osmotic water flow, allows for the computation of important physiological variables including the extracellular voltage (DC) shift. A systematic parameter study reveals that glia can act as both a disperser and buffer of potassium in SD propagation. Furthermore, we show that the timing of the DC shift with respect to extracellular [Formula: see text] rise is highly dependent on glial parameters, a result with implications for the identification of the propagating mechanism of SD.

Large extracellular space leads to neuronal susceptibility to ischemic injury in a Na+/K+ pumps-dependent manner

The extent of anoxic depolarization (AD), the initial electrophysiological event during ischemia, determines the degree of brain region-specific neuronal damage. Neurons in higher brain regions exhibiting nonreversible, strong AD are more susceptible to ischemic injury as compared to cells in lower brain regions that exhibit reversible, weak AD. While the contrasting ADs in different brain regions in response to oxygen-glucose deprivation (OGD) is well established, the mechanism leading to such differences is not clear. Here we use computational modeling to elucidate the mechanism behind the brain region-specific recovery from AD. Our extended Hodgkin-Huxley (HH) framework consisting of neural spiking dynamics, processes of ion accumulation, and ion homeostatic mechanisms unveils that glial-vascular K(+) clearance and Na(+)/K(+)-exchange pumps are key to the cell's recovery from AD. Our phase space analysis reveals that the large extracellular space in the upper brain regions leads to impaired Na(+)/K(+)-exchange pumps so that they function at lower than normal capacity and are unable to bring the cell out of AD after oxygen and glucose is restored.

Computational modeling of neurostimulation in brain diseases

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

Dynamics from seconds to hours in Hodgkin-Huxley model with time-dependent ion concentrations and buffer reservoirs

The classical Hodgkin-Huxley (HH) model neglects the time-dependence of ion concentrations in spiking dynamics. The dynamics is therefore limited to a time scale of milliseconds, which is determined by the membrane capacitance multiplied by the resistance of the ion channels, and by the gating time constants. We study slow dynamics in an extended HH framework that includes time-dependent ion concentrations, pumps, and buffers. Fluxes across the neuronal membrane change intra- and extracellular ion concentrations, whereby the latter can also change through contact to reservoirs in the surroundings. Ion gain and loss of the system is identified as a bifurcation parameter whose essential importance was not realized in earlier studies. Our systematic study of the bifurcation structure and thus the phase space structure helps to understand activation and inhibition of a new excitability in ion homeostasis which emerges in such extended models. Also modulatory mechanisms that regulate the spiking rate can be explained by bifurcations. The dynamics on three distinct slow times scales is determined by the cell volume-to-surface-area ratio and the membrane permeability (seconds), the buffer time constants (tens of seconds), and the slower backward buffering (minutes to hours). The modulatory dynamics and the newly emerging excitable dynamics corresponds to pathological conditions observed in epileptiform burst activity, and spreading depression in migraine aura and stroke, respectively.

Linking a genetic defect in migraine to spreading depression in a computational model

Familial hemiplegic migraine (FHM) is a rare subtype of migraine with aura. A mutation causing FHM type 3 (FHM3) has been identified in SCN1A encoding the Nav1.1 Na⁺ channel. This genetic defect affects the inactivation gate. While the Na⁺ tail currents following voltage steps are consistent with both hyperexcitability and hypoexcitability, in this computational study, we investigate functional consequences beyond these isolated events. Our extended Hodgkin–Huxley framework establishes a connection between genotype and cellular phenotype, i.e., the pathophysiological dynamics that spans over multiple time scales and is relevant to migraine with aura. In particular, we investigate the dynamical repertoire from normal spiking (milliseconds) to spreading depression and anoxic depolarization (tens of seconds) and show that FHM3 mutations render gray matter tissue more vulnerable to spreading depression despite opposing effects associated with action potential generation. We conclude that the classification in terms of hypoexcitability vs. hyperexcitability is too simple a scheme. Our mathematical analysis provides further basic insight into also previously discussed criticisms against this scheme based on psychophysical and clinical data.

Bistable dynamics underlying excitability of ion homeostasis in neuron models

When neurons fire action potentials, dissipation of free energy is usually not directly considered, because the change in free energy is often negligible compared to the immense reservoir stored in neural transmembrane ion gradients and the long-term energy requirements are met through chemical energy, i.e., metabolism. However, these gradients can temporarily nearly vanish in neurological diseases, such as migraine and stroke, and in traumatic brain injury from concussions to severe injuries. We study biophysical neuron models based on the Hodgkin-Huxley (HH) formalism extended to include time-dependent ion concentrations inside and outside the cell and metabolic energy-driven pumps. We reveal the basic mechanism of a state of free energy-starvation (FES) with bifurcation analyses showing that ion dynamics is for a large range of pump rates bistable without contact to an ion bath. This is interpreted as a threshold reduction of a new fundamental mechanism of ionic excitability that causes a long-lasting but transient FES as observed in pathological states. We can in particular conclude that a coupling of extracellular ion concentrations to a large glial-vascular bath can take a role as an inhibitory mechanism crucial in ion homeostasis, while the Na⁺/K⁺ pumps alone are insufficient to recover from FES. Our results provide the missing link between the HH formalism and activator-inhibitor models that have been successfully used for modeling migraine phenotypes, and therefore will allow us to validate the hypothesis that migraine symptoms are explained by disturbed function in ion channel subunits, Na⁺/K⁺ pumps, and other proteins that regulate ion homeostasis.

Mathematical approaches to modeling of cortical spreading depression

Migraine with aura (MwA) is a debilitating disease that afflicts about 25%-30% of migraine sufferers. During MwA, a visual illusion propagates in the visual field, then disappears, and is followed by a sustained headache. MwA was conjectured by Lashley to be related to some neurological phenomenon. A few years later, Leão observed electrophysiological waves in the brain that are now known as cortical spreading depression (CSD). CSD waves were soon conjectured to be the neurological phenomenon underlying MwA that had been suggested by Lashley. However, the confirmation of the link between MwA and CSD was not made until 2001 by Hadjikhani et al. [Proc. Natl. Acad. Sci. U.S.A. 98, 4687-4692 (2001)] using functional MRI techniques. Despite the fact that CSD has been studied continuously since its discovery in 1944, our detailed understandings of the interactions between the mechanisms underlying CSD waves have remained elusive. The connection between MwA and CSD makes the understanding of CSD even more compelling and urgent. In addition to all of the information gleaned from the many experimental studies on CSD since its discovery, mathematical modeling studies provide a general and in some sense more precise alternative method for exploring a variety of mechanisms, which may be important to develop a comprehensive picture of the diverse mechanisms leading to CSD wave instigation and propagation. Some of the mechanisms that are believed to be important include ion diffusion, membrane ionic currents, osmotic effects, spatial buffering, neurotransmitter substances, gap junctions, metabolic pumps, and synaptic connections. Discrete and continuum models of CSD consist of coupled nonlinear differential equations for the ion concentrations. In this review of the current quantitative understanding of CSD, we focus on these modeling paradigms and various mechanisms that are felt to be important for CSD.

Towards dynamical network biomarkers in neuromodulation of episodic migraine

Computational methods have complemented experimental and clinical neurosciences and led to improvements in our understanding of the nervous systems in health and disease. In parallel, neuromodulation in form of electric and magnetic stimulation is gaining increasing acceptance in chronic and intractable diseases. In this paper, we firstly explore the relevant state of the art in fusion of both developments towards translational computational neuroscience. Then, we propose a strategy to employ the new theoretical concept of dynamical network biomarkers (DNB) in episodic manifestations of chronic disorders. In particular, as a first example, we introduce the use of computational models in migraine and illustrate on the basis of this example the potential of DNB as early-warning signals for neuromodulation in episodic migraine.

A mathematical model of the metabolic and perfusion effects on cortical spreading depression

Cortical spreading depression (CSD) is a slow-moving ionic and metabolic disturbance that propagates in cortical brain tissue. In addition to massive cellular depolarizations, CSD also involves significant changes in perfusion and metabolism-aspects of CSD that had not been modeled and are important to traumatic brain injury, subarachnoid hemorrhage, stroke, and migraine. In this study, we develop a mathematical model for CSD where we focus on modeling the features essential to understanding the implications of neurovascular coupling during CSD. In our model, the sodium-potassium-ATPase, mainly responsible for ionic homeostasis and active during CSD, operates at a rate that is dependent on the supply of oxygen. The supply of oxygen is determined by modeling blood flow through a lumped vascular tree with an effective local vessel radius that is controlled by the extracellular potassium concentration. We show that during CSD, the metabolic demands of the cortex exceed the physiological limits placed on oxygen delivery, regardless of vascular constriction or dilation. However, vasoconstriction and vasodilation play important roles in the propagation of CSD and its recovery. Our model replicates the qualitative and quantitative behavior of CSD--vasoconstriction, oxygen depletion, extracellular potassium elevation, prolonged depolarization--found in experimental studies. We predict faster, longer duration CSD in vivo than in vitro due to the contribution of the vasculature. Our results also help explain some of the variability of CSD between species and even within the same animal. These results have clinical and translational implications, as they allow for more precise in vitro, in vivo, and in silico exploration of a phenomenon broadly relevant to neurological disease.

Diffusing substances during spreading depolarization: analytical expressions for propagation speed, triggering, and concentration time courses

Spreading depolarization (SD) is an important phenomenon in stroke and migraine. However, the processes underlying the propagation of SD are still poorly understood, and an elementary model that is both physiological and quantitative is lacking. We show that, during the onset and propagation of SD, the concentration time courses of excitatory substances such as potassium and glutamate can be described with a reaction-diffusion equation. This equation contains four physiological parameters: (1) a concentration threshold for excitation; (2) a release rate; (3) a removal rate; and (4) an effective diffusion constant. Solving this equation yields expressions for the propagation velocity, concentration time courses, and the minimum stimulus that can trigger SD. This framework allows for analyzing experimental results in terms of these four parameters. The derived time courses are validated with measurements of potassium in rat brain tissue.

Self-terminating wave patterns and self-organized pacemakers in a phenomenological model of spreading depression

A simple reaction-diffusion model of spreading depression (SD) is presented. Its local dynamics are governed by two activator and two inhibitor variables that provide an extremely simplified description of the mutual interaction between the neurons and extracellular space. This interaction is realized by the substances in the extracellular space that are increasing excitability of the neurons that have released them and are diffusing to the neighboring neurons, thereby spreading this excitation. Typical dynamic patterns of simulated activity are presented. The focus is laid on the case where response of the extracellular medium is relatively fast, and retracting waves, spiral-shaped waves, and autonomous pacemakers are observed, which is in good agreement with experimental observations. The underlying mechanisms are found to be related to switching between the local bi-stable, excitable, and self-sustained dynamics in the simulated medium. This article is part of a Special Issue entitled: Neural Coding.