Enabling Electric Field Model of Microscopically Realistic Brain

Abstract

Modeling brain stimulation at the microscopic scale may reveal new paradigms for various stimulation modalities. We present the largest map to date of extracellular electric field distributions within a layer L2/L3 mouse primary visual cortex brain sample. This was enabled by the automated analysis of serial section electron microscopy images with improved handling of image defects, covering a volume of 250 x 140 x 90 m3. The map was obtained by applying a uniform brain stimulation electric field at three different polarizations and accurately computing microscopic field perturbations using the boundary element fast multipole method. We used the map to identify the effect of microscopic field perturbations on the activation thresholds of individual neurons. Previous relevant studies modeled a macroscopically homogeneous cortical volume. Our result shows that the microscopic field perturbations - an electric field spatial noise with a mean value of zero - only modestly influence the macroscopically predicted stimulation field strengths necessary for neuronal activation. The thresholds do not change by more than 10% on average. Under the stated limitations and assumptions of our method, this result justifies the conventional theory of "invisible" neurons embedded in a macroscopic brain model for transcranial magnetic and transcranial electrical stimulation. However, our result is solely sample-specific and largely neglects the effect of the microcapillary network. Furthermore, we only considered the uniform impressed field and a single- pulse stimulation time course. Significance statementThis study is arguably the first attempt to model brain stimulation at the microscopic scale, enabled by automated analysis of modern scanning electron microscopy images of the brain. It concentrates on modeling microscopic perturbations of the extracellular electric field caused by the physical cell structure and is applicable to any type of brain stimulation. Data availability statementPost-processed cell CAD models (383, stl format), microcapillary CAD models (34, stl format), post-processed neuron morphologies (267, swc format), extracellular electric field and potential distributions at different polarizations (267x3, MATLAB format), *.ses projects files for biophysical modeling with Neuron software (267x2, Neuron format), and computed neuron activating thresholds at different conditions (267x8, Excel tables, without the sample polarization correction from Section 2.8) are made available online through BossDB, a volumetric open-source database for 3D and 4D neuroscience data.