Influence of Neuronal Morphology on the Shape of Extracellular Recordings With Microelectrode Arrays: A Finite Element Analysis

Finite-element modeling of neuromodulation via controlled delivery of potassium ions using conductive polymer-coated microelectrodes

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.

Multiformity of extracellular microelectrode recordings from Aδ neurons in the dorsal root ganglia: a computational modeling study

Microelectrodes serve as a fundamental tool in electrophysiology research throughout the nervous system, providing a means of exploring neural function with a high resolution of neural firing information. We constructed a hybrid computational model using the finite element method and multicompartment cable models to explore factors that contribute to extracellular voltage waveforms that are produced by sensory pseudounipolar neurons, specifically smaller A-type neurons, and that are recorded by microelectrodes in dorsal root ganglia. The finite element method model included a dorsal root ganglion, surrounding tissues, and a planar microelectrode array. We built a multicompartment neuron model with multiple trajectories of the glomerular initial segment found in many A-type sensory neurons. Our model replicated both the somatic intracellular voltage profile of Aδ low-threshold mechanoreceptor neurons and the unique extracellular voltage waveform shapes that are observed in experimental settings. Results from this model indicated that tortuous glomerular initial segment geometries can introduce distinct multiphasic properties into a neuron's recorded waveform. Our model also demonstrated how recording location relative to specific microanatomical components of these neurons, and recording distance from these components, can contribute to additional changes in the multiphasic characteristics and peak-to-peak voltage amplitude of the waveform. This knowledge may provide context for research employing microelectrode recordings of pseudounipolar neurons in sensory ganglia, including functional mapping and closed-loop neuromodulation. Furthermore, our simulations gave insight into the neurophysiology of pseudounipolar neurons by demonstrating how the glomerular initial segment aids in increasing the resistance of the stem axon and mitigating rebounding somatic action potentials.NEW & NOTEWORTHY We built a computational model of sensory neurons in the dorsal root ganglia to investigate factors that influence the extracellular waveforms recorded by microelectrodes. Our model demonstrates how the unique structure of these neurons can lead to diverse and often multiphasic waveform profiles depending on the location of the recording contact relative to microanatomical neural components. Our model also provides insight into the neurophysiological function of axon glomeruli that are often present in these neurons.

Positive and biphasic extracellular waveforms correspond to return currents and axonal spikes

Multiple biophysical mechanisms may generate non-negative extracellular waveforms during action potentials, but the origin and prevalence of positive spikes and biphasic spikes in the intact brain are unknown. Using extracellular recordings from densely-connected cortical networks in freely-moving mice, we find that a tenth of the waveforms are non-negative. Positive phases of non-negative spikes occur in synchrony or just before wider same-unit negative spikes. Narrow positive spikes occur in isolation in the white matter. Isolated biphasic spikes are narrower than negative spikes, occurring right after spikes of verified inhibitory units. In CA1, units with dominant non-negative spikes exhibit place fields, phase precession, and phase-locking to ripples. Thus, near-somatic narrow positive extracellular potentials correspond to return currents, and isolated non-negative spikes correspond to axonal potentials. Identifying non-negative extracellular waveforms that correspond to non-somatic compartments during spikes can enhance the understanding of physiological and pathological neural mechanisms in intact animals.