Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Blog
    • Collections
    • Podcast
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • ABOUT
    • Overview
    • Editorial Board
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
  • SUBMIT

User menu

Search

  • Advanced search
eNeuro

eNeuro

Advanced Search

 

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Blog
    • Collections
    • Podcast
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • ABOUT
    • Overview
    • Editorial Board
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
  • SUBMIT
PreviousNext
Research ArticleResearch Article: New Research, Disorders of the Nervous System

Multiscale Computer Modeling of Spreading Depolarization in Brain Slices

Craig Kelley, Adam J. H. Newton, Sabina Hrabetova, Robert A. McDougal and William W Lytton
eNeuro 4 August 2022, 9 (4) ENEURO.0082-22.2022; DOI: https://doi.org/10.1523/ENEURO.0082-22.2022
Craig Kelley
1Program in Biomedical Engineering, SUNY Downstate Health Sciences University & NYU Tandon School of Engineering, Brooklyn, NY, 11203
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Adam J. H. Newton
2Department of Physiology and Pharmacology, SUNY Downstate Health Sciences University, Brooklyn, New York 11203
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sabina Hrabetova
3Department of Cell Biology, SUNY Downstate Health Sciences University, Brooklyn, New York 11203
4Robert F. Furchgott Center for Neural and Behavioral Science, SUNY Downstate Health Sciences University, Brooklyn, New York 11203
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Robert A. McDougal
5Department of Biostatistics, Yale University, New Haven, Connecticut 06513
6Yale Center for Medical Informatics, Yale University, New Haven, Connecticut 06513
7Program in Computational Biology and Bioinformatics, Yale University, New Haven, Connecticut 06513
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Robert A. McDougal
William W Lytton
2Department of Physiology and Pharmacology, SUNY Downstate Health Sciences University, Brooklyn, New York 11203
4Robert F. Furchgott Center for Neural and Behavioral Science, SUNY Downstate Health Sciences University, Brooklyn, New York 11203
8Department of Neurology, SUNY Downstate Health Sciences University, Brooklyn, New York 11203
9Department of Neurology, Kings County Hospital Center, Brooklyn, New York 11203
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Article Figures & Data

Figures

  • Tables
  • Movies
  • Extended Data
  • Figure 1.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 1.

    Multiscale model expanded. Tissue scale: a few of the 36 ⋅ 103 neurons (pink circles) embedded in the ECS of a brain slice submerged in a bath solution where ion and O2 concentrations were held constant. Glia are not explicitly modeled, but instead were represented as a field of sinks in every ECS voxel. Cell scale: each neuron had ion channels, 2 coexchangers; Na+/K+ pump (asterisk indicates ATP/O2 dependence). Ions were well mixed within each neuron (no intracellular diffusion). Protein scale: table (right) indicates species that control the activity of the intrinsic mechanisms in neurons and in glial field. Ion scale: ions diffused between ECS voxels by Fick’s law using diffusion coefficients in Table 1.

  • Figure 2.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 2.

    A small bolus of applied K+ initiates spreading K+ and depolarization waves in perfused slice. A, [K+]ECS averaged across slice depth (400 μm) at 4 time points during SD. B, Spike raster plot of 250 randomly selected neurons (of 36 ⋅103) during SD. Cells are ordered on y-axis by their radial distance from the center of the K+ bolus. Blank area under spikes represents region of spreading depression. Baseline values: [O2] = 0.1 mm; αECS = 0.2; λECS = 1.6; [Cl–]ECS = 130.0 mm; [Cl–]i = 6.0 mm.

  • Figure 3.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 3.

    Diameter (bottom, x-axis) and concentration (top, x-axis) of the K+ bolus had minor effects on K+ wave speed. Mean and SDs (n = 5) for K+ wave speed versus bolus diameter (solid line; bolus [K+]ECS = 70 mm; 5 random cell position initializations) and versus bolus concentration (dashed line; diameter, 200 μm).

  • Figure 4.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 4.

    Concentrations at 7 radial locations measured during SD in perfused 400 μm slice. A, Cell within the K+ bolus (37.5 μm) produced a single spike; cells farther out fired a burst. Cells remained in depolarization block for the remainder of the simulation (10 s). B–D, Intracellular ion concentrations. E–H, Extracellular O2 and ion concentrations in neighboring ECS voxels. Movie 1 shows extracellular ion and O2 concentrations across the slice, as well as neuronal spiking (white dots) from 250 neurons during the course of SD (most easily seen with slowed playback).

  • Figure 5.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 5.

    Hypoxia, propionate, and dynamic ECS increased SD speed principally through αECS reduction. A, Radial K+ wave position over time during SD in perfused (Movie 1), hypoxic, HSD (Movie 2), propionate conditions, and with dynamic changes in αECS . Hypoxia, propionate, and dynamic changes in αECS facilitated propagation. B, Radial position of SD wave represented by time to first spike in 126 selected cells at different distances from center. C, K+ wave speeds with individual parameter changes (Fig. 2). αECS had the greatest impact on SD speed over a physiologically plausible range (x-axis ranges: [O2] = 0.01–0.1 mm; αECS = 0.07–0.42; λECS = 1.4–2.0; [Cl–]ECS:[Cl–]I = 3.0:65.0–6.0:130.0 mm).

  • Figure 6.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 6.

    K+ wave propagation speed proportional to total neuronal surface area in slice. SD initiated in perfused and hypoxic slices (no reperfusion) by introducing a 100-μm-radius central 70 mm spherical K+ bolus. A–E, Effects of varying S/V of each cell while maintaining a constant βnrn (A); varying βnrn while keeping S/V constant, thus allowing Snrn to vary (B); varying cell density while keeping constant Snrn and volnrn, thus allowing βnrn to vary (C); and varying cell density while keeping βnrn and S/V constant, thus allowing Snrn and volnrn to vary (D). E, Pooled results: K+ wave speed increased linearly with total neuronal surface area in both perfused and hypoxic slices. Hypoxia increased wave speed across conditions (0 speed indicates no SD).

  • Figure 7.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 7.

    Slice thickness effects on SD propagation. K+ wave speed during SD in perfused and hypoxic slices of various thicknesses. SD could not be initiated in very thin (100 μm) slices when perfused but could in hypoxic slices. For perfused slices, K+ wave speed increased with slice thickness between 200 and 400 μm then saturated for slices of greater thickness. A similar pattern was observed in the hypoxic slices with consistently faster K+ wave speeds compared with in perfused slices.

  • Figure 8.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 8.

    Depth dependence of SD propagation in 400-μm-thick perfused slice. A, Spread of K+ wave through slice core (±50 μm from center). B, Wave of mildly elevated K+ reached the periphery (within 50 μm of bath) from the central bolus. C, Spread of SD through the slice core. Voltages in slice core (color map) with spike rasters for 120-cell subset overlaid in white. Neurons in the core showed typical SD voltage dynamics (bursting followed by depolarization block). D, Voltages and raster at periphery; cells show regular spiking patterns at 10–70 Hz (Vmemb averaged across all cells in 25 × 25 × 100 μm voxels at center or periphery).

Tables

  • Figures
  • Movies
  • Extended Data
    • View popup
    Table 1

    Diffusion coefficients and baseline concentrations for ions and O2 in perfused slice (Samson et al., 2003; Haynes, 2014; Wei et al., 2014)

    SpeciesD
    (⋅10–5 cm2/s)
    Intracellular
    concentration (mm)
    ECS
    concentration (mm)
    K+2.62140.03.5
    Na+ 1.7818.0144.0
    Cl–2.106.0130.0
    O23.300.10.1

Movies

  • Figures
  • Tables
  • Extended Data
  • Movie 1.

    Extracellular ion and O2 concentrations across the slice averaged over depth, as well as neuronal spiking (white dots) from 250 neurons during the course of SD in a perfused slice. The spread of spiking and the K+ wave can be seen in real time. We recommend downloading the file and using slower playback to visualize the spread of hypoxia.

  • Movie 2.

    Extracellular ion and O2 concentrations across the slice averaged over depth, as well as neuronal spiking (white dots) from 250 neurons during the course of hypoxic SD-like depolarization. The spread of spiking and the K+ wave can be seen in real time. We recommend downloading the file and using slower playback to visualize the delay between initiating the simulation and the spread of the K+ wave.

Extended Data

  • Figures
  • Tables
  • Movies
  • Extended Data 1

    Simulation code: our tissue-scale model of SD in brain slices is available on GitHub. We used the NEURON simulator reaction-diffusion framework to implement embed thousands of neurons (based on the model from Wei et al., 2014) in the extracellular space of a brain slice, which is itself embedded in an bath solution. We initiated SD in the slice by elevating extracellular K+ in a spherical region at the center of the slice. The effects of hypoxia and propionate on the slice were modeled by appropriate changes to the volume fraction and tortuosity of the extracellular space and oxygen/chloride concentrations. Users need to install NEURON, and we recommend using MPI to parallelize simulations. Download Extended Data 1, ZIP file.

Back to top

In this issue

eneuro: 9 (4)
eNeuro
Vol. 9, Issue 4
July/August 2022
  • Table of Contents
  • Index by author
  • Ed Board (PDF)
Email

Thank you for sharing this eNeuro article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Multiscale Computer Modeling of Spreading Depolarization in Brain Slices
(Your Name) has forwarded a page to you from eNeuro
(Your Name) thought you would be interested in this article in eNeuro.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Multiscale Computer Modeling of Spreading Depolarization in Brain Slices
Craig Kelley, Adam J. H. Newton, Sabina Hrabetova, Robert A. McDougal, William W Lytton
eNeuro 4 August 2022, 9 (4) ENEURO.0082-22.2022; DOI: 10.1523/ENEURO.0082-22.2022

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Share
Multiscale Computer Modeling of Spreading Depolarization in Brain Slices
Craig Kelley, Adam J. H. Newton, Sabina Hrabetova, Robert A. McDougal, William W Lytton
eNeuro 4 August 2022, 9 (4) ENEURO.0082-22.2022; DOI: 10.1523/ENEURO.0082-22.2022
del.icio.us logo Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Significance Statement
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Acknowledgments
    • Footnotes
    • References
    • Synthesis
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • extracellular space
  • hypoxia
  • NEURON
  • reaction-diffusion
  • spreading depolarization
  • spreading depression

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Article: New Research

  • Opponent Learning with Different Representations in the Cortico-Basal Ganglia Circuits
  • Nonspiking Interneurons in the Drosophila Antennal Lobe Exhibit Spatially Restricted Activity
  • Pattern of Driver-Like Input onto Neurons of the Mouse Ventral Lateral Geniculate Nucleus
Show more Research Article: New Research

Disorders of the Nervous System

  • Brain FNDC5/irisin expression in patients and mouse models of major depression
  • Increased physiological GDNF levels have no effect on dopamine neuron protection and restoration in a proteasome inhibition mouse model of Parkinson's disease
  • Microglial Expression of the Wnt Signaling Modulator DKK2 Differs between Human Alzheimer’s Disease Brains and Mouse Neurodegeneration Models
Show more Disorders of the Nervous System

Subjects

  • Disorders of the Nervous System

  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Latest Articles
  • Issue Archive
  • Blog
  • Browse by Topic

Information

  • For Authors
  • For the Media

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
  • Feedback
(eNeuro logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
eNeuro eISSN: 2373-2822

The ideas and opinions expressed in eNeuro do not necessarily reflect those of SfN or the eNeuro Editorial Board. Publication of an advertisement or other product mention in eNeuro should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in eNeuro.