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

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.

How does spreading depression spread? Physiology and modeling

Spreading depression (SD) is a wave phenomenon in gray matter tissue. Locally, it is characterized by massive redistribution of ions across cell membranes. As a consequence, there is sustained membrane depolarization and tissue polarization that depress any normal electrical activity. Despite these dramatic events, SD remains difficult to observe in humans noninvasively, which, for long, has slowed advances in this field. The growing appreciation of its clinical importance in migraine and stroke is therefore consistent with an increasing need for computational methods that tackle the complexity of the problem at multiple levels. In this review, we focus on mathematical tools to investigate the question of spread and its two complementary aspects: What are the physiological mechanisms and what is the spatial extent of SD in the cortex? This review discusses two types of models used to study these two questions, namely, Hodgkin-Huxley type and generic activator-inhibitor models, and the recent advances in techniques to link them.

Cortical hot spots and labyrinths: why cortical neuromodulation for episodic migraine with aura should be personalized

Stimulation protocols for medical devices should be rationally designed. For episodic migraine with aura we outline model-based design strategies toward preventive and acute therapies using stereotactic cortical neuromodulation. To this end, we regard a localized spreading depression (SD) wave segment as a central element in migraine pathophysiology. To describe nucleation and propagation features of the SD wave segment, we define the new concepts of cortical hot spots and labyrinths, respectively. In particular, we firstly focus exclusively on curvature-induced dynamical properties by studying a generic reaction-diffusion model of SD on the folded cortical surface. This surface is described with increasing level of details, including finally personalized simulations using patient's magnetic resonance imaging (MRI) scanner readings. At this stage, the only relevant factor that can modulate nucleation and propagation paths is the Gaussian curvature, which has the advantage of being rather readily accessible by MRI. We conclude with discussing further anatomical factors, such as areal, laminar, and cellular heterogeneity, that in addition to and in relation to Gaussian curvature determine the generalized concept of cortical hot spots and labyrinths as target structures for neuromodulation. Our numerical simulations suggest that these target structures are like fingerprints, they are individual features of each migraine sufferer. The goal in the future will be to provide individualized neural tissue simulations. These simulations should predict the clinical data and therefore can also serve as a test bed for exploring stereotactic cortical neuromodulation.

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.