On the upper threshold phenomenon of extracellular neural stimulation

Neural activity of retinal ganglion cells under continuous, dynamically-modulated high frequency electrical stimulation

Objective.Current retinal prosthetics are limited in their ability to precisely control firing patterns of functionally distinct retinal ganglion cell (RGC) types. The aim of this study was to characterise RGC responses to continuous, kilohertz-frequency-varying stimulation to assess its utility in controlling RGC activity.Approach.We usedin vitropatch-clamp experiments to assess electrically-evoked ON and OFF RGC responses to frequency-varying pulse train sequences. In each sequence, the stimulation amplitude was kept constant while the stimulation frequency (0.5-10 kHz) was changed every 40 ms, in either a linearly increasing, linearly decreasing or randomised manner. The stimulation amplitude across sequences was increased from 10 to 300µA.Main results.We found that continuous stimulation without rest periods caused complex and irreproducible stimulus-response relationships, primarily due to strong stimulus-induced response adaptation and influence of the preceding stimulus frequency on the response to a subsequent stimulus. In addition, ON and OFF populations showed different sensitivities to continuous, frequency-varying pulse trains, with OFF cells generally exhibiting more dependency on frequency changes within a sequence. Finally, the ability to maintain spiking behaviour to continuous stimulation in RGCs significantly reduced over longer stimulation durations irrespective of the frequency order.Significance.This study represents an important step in advancing and understanding the utility of continuous frequency modulation in controlling functionally distinct RGCs. Our results indicate that continuous, kHz-frequency-varying stimulation sequences provide very limited control of RGC firing patterns due to inter-dependency between adjacent frequencies and generally, different RGC types do not display different frequency preferences under such stimulation conditions. For future stimulation strategies using kHz frequencies, careful consideration must be given to design appropriate pauses in stimulation, stimulation frequency order and the length of continuous stimulation duration.

Implications of Neural Plasticity in Retinal Prosthesis

Retinal degenerative diseases such as retinitis pigmentosa cause a progressive loss of photoreceptors that eventually prevents the affected person from perceiving visual sensations. The absence of a visual input produces a neural rewiring cascade that propagates along the visual system. This remodeling occurs first within the retina. Then, subsequent neuroplastic changes take place at higher visual centers in the brain, produced by either the abnormal neural encoding of the visual inputs delivered by the diseased retina or as the result of an adaptation to visual deprivation. While retinal implants can activate the surviving retinal neurons by delivering electric current, the unselective activation patterns of the different neural populations that exist in the retinal layers differ substantially from those in physiologic vision. Therefore, artificially induced neural patterns are being delivered to a brain that has already undergone important neural reconnections. Whether or not the modulation of this neural rewiring can improve the performance for retinal prostheses remains a critical question whose answer may be the enabler of improved functional artificial vision and more personalized neurorehabilitation strategies.

A finite element method framework to model extracellular neural stimulation

OBJECTIVE

Increasing complexity in extracellular stimulation experiments and neural implant design also requires realistic computer simulations capable of modeling the neural activity of nerve cells under the influence of an electrical stimulus. Classical model approaches are often based on simplifications, are not able to correctly calculate the electric field generated by complex electrode designs, and do not consider electrical effects of the cell on its surrounding. A more accurate approach is the finite element method (FEM), which provides necessary techniques to solve the Poisson equation for complex geometries under consideration of electrical tissue properties. Especially in situations where neurons experience large and non-symmetric extracellular potential gradients, a FEM solution that implements the cell membrane model can improve the computer simulation results. To investigate the response of neurons in an electric field generated by complex electrode designs, a FEM framework for extracellular stimulation was developed in COMSOL.

APPROACH

Methods to implement morphologically- and biophysically-detailed neurons including active Hodgkin-Huxley (HH) cell membrane dynamics as well as the stimulation setup are described in detail. Covered methods are (i) development of cell and electrode geometries including meshing strategies, (ii) assignment of physics for the conducting spaces and the realization of active electrodes, (iii) implementation of the HH model, and (iv) coupling of the physics to get a fully described model.

MAIN RESULTS

Several implementation examples are briefly presented: (i) a full FEM implementation of a HH model cell stimulated with a honeycomb electrode, (ii) the electric field of a cochlear electrode placed inside the cochlea, and (iii) a proof of concept implementation of a detailed double-cable cell membrane model for myelinated nerve fibers.

SIGNIFICANCE

The presented concepts and methods provide basic and advanced techniques to realize a full FEM framework for innovative studies of neural excitation in response to extracellular stimulation.

Block Phenomena During Electric Micro-Stimulation of Pyramidal Cells and Retinal Ganglion Cells

Electric micro-stimulation of the nervous system is a means to restore various body functions. The stimulus amplitude necessary to generate action potentials, the lower threshold (LT), is well characterized for many neuronal populations. However, electric overstimulation above an upper threshold (UT) prevents action potential generation and therefore hinders optimal neuro-rehabilitation. Previous studies demonstrated the impact of the UT in micro-stimulation of retinal ganglion cells (RGCs). The observed phenomenon is mostly explained by (i) reversed sodium ion flow in the soma membrane, and (ii) anodal surround block that hinders spike conduction in strongly hyperpolarized regions of the axon at high stimulus intensities. However, up to now, no detailed study of the nature of these phenomena has been presented, particularly for different cell types. Here, we present computational analyses of LT and UT for layer 5 pyramidal cells (PCs) as well as alpha RGCs. Model neurons were stimulated in close vicinity to the cell body and LTs and UTs as well as the ratio UT/LT were compared. Aside from a simple point source electrode and monophasic stimuli also realistic electrode and pulse configurations were examined. The analysis showed: (i) in RGCs, the soma contributed to action potential initiation and block for small electrode distances, whereas in PCs the soma played no role in LTs or UTs. (ii) In both cell types, action potential always initiated within the axon initial segment at LT. (iii) In contrast to a complete block of spike conductance at UT that occurred in RGCs, an incomplete block of spiking appeared in PC axon collaterals. (iv) PC axon collateral arrangement influenced UTs but had small impact on LTs. (v) Population responses of RGCs change from circular regions of activation to ring-shaped patterns for increasing stimulus amplitude. A better understanding of the stimulation window that can reliably activate target neurons will benefit the future development of neuroprostheses.

Comparison of responses of visual cortical neurons in the mouse to intraocular and extraocular electric stimulation of the retina

Retinal implants offered the promise of restoring functional vision to the blind via the delivery of electrical stimulation to the retina. To enhance the efficacy of these devices, stimulation should elicit neural responses that are similar to the responses that occur naturally in the retina as these have the best chance of carrying a robust signal to visual cortex. A corollary of this is that the responses that arise in visual cortical neurons can be used to compare the effectiveness of different stimulation strategies in the retina. Here, we studied how visual cortical neurons in the mouse respond to monophasic cathodal and anodal electric stimulation delivered via a wire electrode positioned on the outer surface of the eye (extraocular) or within the vitreous cavity of the eye (intraocular). Responses of visual cortical neurons were recorded from primary visual cortex on the side contralateral to the stimulatated eye. For both stimulation modalities, response patterns consisted of a brief burst of spikes followed by a 400-500 ms period of inhibition. Both modalities also elicited stronger responses to cathodal stimuli (vs. anodal). The preferential sensitivity to cathodal stimuli is similar to that of epiretinal stimulation (anodal stimuli are more effective with sub-retinal stimulation) suggest the extraocular approach mirrors epiretinal stimulation. Extraocular stimulation also showed some response characteristics that were different from those observed in the retina, e.g., at very strong amplitudes, cathodal and anodal stimulation produced similar responses.

Mediating Retinal Ganglion Cell Spike Rates Using High-Frequency Electrical Stimulation

Recent retinal studies have directed more attention to sophisticated stimulation strategies based on high-frequency (>1.0 kHz) electrical stimulation (HFS). In these studies, each retinal ganglion cell (RGC) type demonstrated a characteristic stimulus-strength-dependent response to HFS, offering the intriguing possibility of focally targeting retinal neurons to provide useful visual information by retinal prosthetics. Ionic mechanisms are known to affect the responses of electrogenic cells during electrical stimulation. However, how these mechanisms affect RGC responses is not well understood at present, particularly when applying HFS. Here, we investigate this issue via an in silico model of the RGC. We calibrate and validate the model using an in vitro retinal preparation. An RGC model based on accurate biophysics and realistic representation of cell morphology, was used to investigate how RGCs respond to HFS. The model was able to closely replicate the stimulus-strength-dependent suppression of RGC action potentials observed experimentally. Our results suggest that spike inhibition during HFS is due to local membrane hyperpolarization caused by outward membrane currents near the stimulus electrode. In addition, the extent of HFS-induced inhibition can be largely altered by the intrinsic properties of the inward sodium current. Finally, stimulus-strength-dependent suppression can be modulated by a wide range of stimulation frequencies, under generalized electrode placement conditions. In vitro experiments verified the computational modeling data. This modeling and experimental approach can be extended to further our understanding on the effects of novel stimulus strategies by simulating RGC stimulus-response profiles over a wider range of stimulation frequencies and electrode locations than have previously been explored.

Analysis of upper threshold mechanisms of spherical neurons during extracellular stimulation

Exceeding a certain stimulation strength can prevent the generation of somatic action potentials, as has been demonstrated in vitro with extracellularly stimulated dorsal root ganglion cells as well as retinal ganglion cells. This phenomenon, termed upper threshold, is currently thought to be a consequence of sodium current reversal in strongly depolarized regions. Here we analyze the contribution of membrane kinetics, using spherical model neurons that are stimulated externally with a microelectrode, in more detail. During extracellular pulse application, the electric field depolarizes one part and hyperpolarizes the other part of the cell. Strong transmembrane currents are generated only in the active depolarized region, changing the overall polarization level. The asymmetric membrane voltage distribution caused by the stimulus strongly influences the cell's behavior during and even after the stimulus. Effects on membrane voltage and transmembrane currents during and after the stimulus are shown and discussed in detail. Aside from the sodium current reversal, two more key mechanisms were identified in causing the upper threshold: strong potassium currents and inactivation of sodium channels. The contributions of the mechanisms involved strongly depend on cell properties, stimulus parameters, and other factors such as temperature. The conclusions presented here are based on several retinal ganglion cell models of the Fohlmeister group, a model with original Hodgkin-Huxley membrane, and a pyramidal cell model. NEW & NOTEWORTHY The upper threshold phenomenon in extracellular stimulation is analyzed in detail for spherical cells. Three main mechanisms were identified that prevent the generation of action potentials at high stimulation strengths: 1) strong potassium currents, 2) inactivating sodium ion channels, and 3) sodium current reversal. Ion channel kinetics in retinal ganglion cells, pyramidal cells, and the original Hodgkin-Huxley model were investigated under the influence of an extracellular stimulus.

Upper stimulation threshold for retinal ganglion cell activation

OBJECTIVE

The existence of an upper threshold in electrically stimulated retinal ganglion cells (RGCs) is of interest because of its relevance to the development of visual prosthetic devices, which are designed to restore partial sight to blind patients. The upper threshold is defined as the stimulation level above which no action potentials (direct spikes) can be elicited in electrically stimulated retina.

APPROACH

We collected and analyzed in vitro recordings from rat RGCs in response to extracellular biphasic (anodic-cathodic) pulse stimulation of varying amplitudes and pulse durations. Such responses were also simulated using a multicompartment model.

MAIN RESULTS

We identified the individual cell variability in response to stimulation and the phenomenon known as upper threshold in all but one of the recorded cells (n  =  20/21). We found that the latencies of spike responses relative to stimulus amplitude had a characteristic U-shape. In silico, we showed that the upper threshold phenomenon was observed only in the soma. For all tested biphasic pulse durations, electrode positions, and pulse amplitudes above lower threshold, a propagating action potential was observed in the distal axon. For amplitudes above the somatic upper threshold, the axonal action potential back-propagated in the direction of the soma, but the soma's low level of hyperpolarization prevented action potential generation in the soma itself.

SIGNIFICANCE

An upper threshold observed in the soma does not prevent spike conductance in the axon.

Closed-Loop Efficient Searching of Optimal Electrical Stimulation Parameters for Preferential Excitation of Retinal Ganglion Cells

The ability for visual prostheses to preferentially activate functionally-distinct retinal ganglion cells (RGCs) is important for improving visual perception. This study investigates the use of high frequency stimulation (HFS) to elicit RGC activation, using a closed-loop algorithm to search for optimal stimulation parameters for preferential ON and OFF RGC activation, resembling natural physiological neural encoding in response to visual stimuli. We evaluated the performance of a wide range of electrical stimulation amplitudes and frequencies on RGC responses in vitro using murine retinal preparations. It was possible to preferentially excite either ON or OFF RGCs by adjusting amplitudes and frequencies in HFS. ON RGCs can be preferentially activated at relatively higher stimulation amplitudes (>150 μA) and frequencies (2-6.25 kHz) while OFF RGCs are activated by lower stimulation amplitudes (40-90 μA) across all tested frequencies (1-6.25 kHz). These stimuli also showed great promise in eliciting RGC responses that parallel natural RGC encoding: ON RGCs exhibited an increase in spiking activity during electrical stimulation while OFF RGCs exhibited decreased spiking activity, given the same stimulation amplitude. In conjunction with the in vitro studies, in silico simulations indicated that optimal HFS parameters could be rapidly identified in practice, whilst sampling spiking activity of relevant neuronal subtypes. This closed-loop approach represents a step forward in modulating stimulation parameters to achieve appropriate neural encoding in retinal prostheses, advancing control over RGC subtypes activated by electrical stimulation.

Impact of Electrode Position on the Elicitation of Sodium Spikes in Retinal Bipolar Cells

Bipolar cells of the magnocellular pathway in the primate retina can generate action potentials because they have an axonal segment with high sodium channel density, comparable to the sodium channel band in retinal ganglion cells or pyramidal cells. The similarity between the non-human primate and the human retina is of interest for the research on retinal implants for the blind, and especially, the conditions to elicit sodium spikes in bipolar cells using extracellular stimulation. A comparison of excitation characteristics of three model neurons, a bipolar cell, a retinal ganglion cell, and a cortical pyramidal cell, demonstrates the similarities and differences regarding stimulation with microelectrodes. Moving a microelectrode parallel to the axon of a neuron commonly allows to generate spikes for every position - and this rule holds both for cathodic and anodic pulses. However, for the simulated bipolar cell anodic pulses cannot generate sodium spikes directly. Further, there is only a small region for electrode placing where extracellular cathodic stimulation causes direct spike initiation in the sodium channel band. For all other positions, a sodium spike can only be generated by antidromic current flow originating from strongly depolarized terminals.

High-amplitude electrical stimulation can reduce elicited neuronal activity in visual prosthesis

Retinal electrostimulation is promising a successful therapy to restore functional vision. However, a narrow stimulating current range exists between retinal neuron excitation and inhibition which may lead to misperformance of visual prostheses. As the conveyance of representation of complex visual scenes may require neighbouring electrodes to be activated simultaneously, electric field summation may contribute to reach this inhibitory threshold. This study used three approaches to assess the implications of relatively high stimulating conditions in visual prostheses: (1) in vivo, using a suprachoroidal prosthesis implanted in a feline model, (2) in vitro through electrostimulation of murine retinal preparations, and (3) in silico by computing the response of a population of retinal ganglion cells. Inhibitory stimulating conditions led to diminished cortical activity in the cat. Stimulus-response relationships showed non-monotonic profiles to increasing stimulating current. This was observed in vitro and in silico as the combined response of groups of neurons (close to the stimulating electrode) being inhibited at certain stimulating amplitudes, whilst other groups (far from the stimulating electrode) being recruited. These findings may explain the halo-like phosphene shapes reported in clinical trials and suggest that simultaneous stimulation in retinal prostheses is limited by the inhibitory threshold of the retinal ganglion cells.