Activation and inhibition of retinal ganglion cells in response to epiretinal electrical stimulation: a computational modelling study

Patient-specific computational models of retinal prostheses

Retinal prostheses stimulate inner retinal neurons to create visual perception for blind patients. Implanted arrays have many small electrodes. Not all electrodes induce perception at the same stimulus amplitude, requiring clinicians to manually establish a visual perception threshold for each one. Phosphenes created by single-electrode stimuli can also vary in shape, size, and brightness. Computational models provide a tool to predict inter-electrode variability and automate device programming. In this study, we created statistical and patient-specific field-cable models to investigate inter-electrode variability across seven epiretinal prosthesis users. Our statistical analysis revealed that retinal thickness beneath the electrode correlated with perceptual threshold, with a significant fixed effect across participants. Electrode-retina distance and electrode impedance also correlated with perceptual threshold for some participants, but these effects varied by individual. We developed a novel method to construct patient-specific field-cable models from optical coherence tomography images. Predictions with these models significantly correlated with perceptual threshold for 80% of participants. Additionally, we demonstrated that patient-specific field-cable models could predict retinal activity and phosphene size. These computational models could be beneficial for determining optimal stimulation settings in silico, circumventing the trial-and-error testing of a large parameter space in clinic.

Patient-specific computational models of retinal prostheses

Retinal prostheses stimulate inner retinal neurons to create visual perception for blind patients. Implanted arrays have many small electrodes, which act as pixels. Not all electrodes induce perception at the same stimulus amplitude, requiring clinicians to manually establish a visual perception threshold for each one. Phosphenes created by single-electrode stimuli can also vary in shape, size, and brightness. Computational models provide a tool to predict inter-electrode variability and automate device programming. In this study, we created statistical and patient-specific field-cable models to investigate inter-electrode variability across seven epiretinal prosthesis users. Our statistical analysis revealed that retinal thickness beneath the electrode correlated with perceptual threshold, with a significant fixed effect across participants. Electrode-retina distance and electrode impedance also correlated with perceptual threshold for some participants, but these effects varied by individual. We developed a novel method to construct patient-specific field-cable models from optical coherence tomography images. Predictions with these models significantly correlated with perceptual threshold for 80% of participants. Additionally, we demonstrated that patient-specific field-cable models could predict retinal activity and phosphene size. These computational models could be beneficial for determining optimal stimulation settings in silico, circumventing the trial-and-error testing of a large parameter space in clinic.

A varying-radius cable equation for the modelling of impulse propagation in excitable fibres

In this study, we present a varying-radius cable equation for nerve fibres taking into account the varying diameter along the neuronal segments. Finite element neuronal models utilising the classical (fixed-radius) and varying-radius cable formulations were compared using simple and realistic morphologies under intra- and extracellular electrical stimulation protocols. We found that the use of the classical cable equation to model intracellular neural electrical stimulation exhibited an error of 17% in a passive resistive cable model with abrupt change in radius from 1 to 2 μm, when compared to the known analytical solution and varying-radius cable formulation. This error was observed to increase substantially using more realistic neuron morphologies and branching structures. In the case of extracellular stimulation however, the difference between the classical and varying-radius formulations was less pronounced, but we expect this difference will increase under more complex stimulation paradigms such as high-frequency stimulation. We conclude that for computational neuroscience applications, it is essential to use the varying-radius cable equation for accurate prediction of neuronal responses under electrical stimulation.

Somatic inhibition by microscopic magnetic stimulation

Electric currents can produce quick, reversible control of neural activity. Externally applied electric currents have been used in inhibiting certain ganglion cells in clinical practices. Via electromagnetic induction, a miniature-sized magnetic coil could provide focal stimulation to the ganglion neurons. Here we report that high-frequency stimulation with the miniature coil could reversibly block ganglion cell activity in marine mollusk Aplysia californica, regardless the firing frequency of the neurons, or concentration of potassium ions around the ganglion neurons. Presence of the ganglion sheath has minimal impact on the inhibitory effects of the coil. The inhibitory effect was local to the soma, and was sufficient in blocking the neuron's functional output. Biophysical modeling confirmed that the miniature coil induced a sufficient electric field in the vicinity of the targeted soma. Using a multi-compartment model of Aplysia ganglion neuron, we found that the high-frequency magnetic stimuli altered the ion channel dynamics that were essential for the sustained firing of action potentials in the soma. Results from this study produces several critical insights to further developing the miniature coil technology for neural control by targeting ganglion cells. The miniature coil provides an interesting neural modulation strategy in clinical applications and laboratory research.

Retinal Drug Delivery: Rethinking Outcomes for the Efficient Replication of Retinal Behavior

The retina is a highly organized structure that is considered to be "an approachable part of the brain." It is attracting the interest of development scientists, as it provides a model neurovascular system. Over the last few years, we have been witnessing significant development in the knowledge of the mechanisms that induce the shape of the retinal vascular system, as well as knowledge of disease processes that lead to retina degeneration. Knowledge and understanding of how our vision works are crucial to creating a hardware-adaptive computational model that can replicate retinal behavior. The neuronal system is nonlinear and very intricate. It is thus instrumental to have a clear view of the neurophysiological and neuroanatomic processes and to take into account the underlying principles that govern the process of hardware transformation to produce an appropriate model that can be mapped to a physical device. The mechanistic and integrated computational models have enormous potential toward helping to understand disease mechanisms and to explain the associations identified in large model-free data sets. The approach used is modulated and based on different models of drug administration, including the geometry of the eye. This work aimed to review the recently used mathematical models to map a directed retinal network.

A Three-Dimensional Microelectrode Array to Generate Virtual Electrodes for Epiretinal Prosthesis Based on a Modeling Study

Despite many advances in the development of retinal prostheses, clinical reports show that current retinal prosthesis subjects can only perceive prosthetic vision with poor visual acuity. A possible approach for improving visual acuity is to produce virtual electrodes (VEs) through electric field modulation. Generating controllable and localized VEs is a crucial factor in effectively improving the perceptive resolution of the retinal prostheses. In this paper, we aimed to design a microelectrode array (MEA) that can produce converged and controllable VEs by current steering stimulation strategies. Through computational modeling, we designed a three-dimensional concentric ring–disc MEA and evaluated its performance with different stimulation strategies. Our simulation results showed that electrode–retina distance (ERD) and inter-electrode distance (IED) can dramatically affect the distribution of electric field. Also the converged VEs could be produced when the parameters of the three-dimensional MEA were appropriately set. VE sites can be controlled by manipulating the proportion of current on each adjacent electrode in a current steering group (CSG). In addition, spatial localization of electrical stimulation can be greatly improved under quasi-monopolar (QMP) stimulation. This study may provide support for future application of VEs in epiretinal prosthesis for potentially increasing the visual acuity of prosthetic vision.

A Patient-Specific Computational Framework for the Argus II Implant

Goal

Retinal prosthesis performance is limited by the variability of elicited phosphenes. The stimulating electrode's position with respect to retinal ganglion cells (RGCs) affects both perceptual threshold and phosphene shape. We created a modeling framework incorporating patient-specific anatomy and electrode location to investigate RGC activation and predict inter-electrode differences for one Argus II user.

Methods

We used ocular imaging to build a three-dimensional finite element model characterizing retinal morphology and implant placement. To predict the neural response to stimulation, we coupled electric fields with multi-compartment cable models of RGCs. We evaluated our model predictions by comparing them to patient-reported perceptual threshold measurements.

Results

Our model was validated by the ability to replicate clinical impedance and threshold values, along with known neurophysiological trends. Inter-electrode threshold differences in silico correlated with in vivo results.

Conclusions

We developed a patient-specific retinal stimulation framework to quantitatively predict RGC activation and better explain phosphene variations.

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.

A continuum model of electrical stimulation of multi-compartmental retinal ganglion cells

A continuum multi-domain model of electrical stimulation of the retina is presented. Each point in the retinal ganglion cell layer could be thought of as representing a single cell, whose biophysics is described using a four-compartment formulation incorporating varying ion channel expressions in the soma, axon initial segment, dendrites and axon. Our continuum model was validated against a discrete morphologically-realistic OFF RGC model, using intra- and extra-cellular electrical stimulation scenarios. Simulations from the continuum model reproduced the same results as that of the discrete model. Our continuum model is the first multi-domain model to represent all main RGC compartments, not just the soma. Moreover, we demonstrated that this model allows the investigation of axonal activation which has been observed to influence the perception of phosphenes.

Progress in artificial vision through suprachoroidal retinal implants

Retinal implants have proven their ability to restore visual sensation to people with degenerative retinopathy, characterized by photoreceptor cell death and the retina's inability to sense light. Retinal bionics operate by electrically stimulating the surviving neurons in the retina, thus triggering the transfer of visual sensory information to the brain. Suprachoroidal implants were first investigated in Australia in the 1950s. In this approach, the neuromodulation hardware is positioned between the sclera and the choroid, thus providing significant surgical and safety benefits for patients, with the potential to maintain residual vision combined with the artificial input from the device. Here we review the latest advances and state of the art devices for suprachoroidal prostheses, highlight future technologies and discuss challenges and perspectives towards improved rehabilitation of vision.

Increasing Electrical Stimulation Efficacy in Degenerated Retina: Stimulus Waveform Design in a Multiscale Computational Model

A computational model of electrical stimulation of the retina is proposed for investigating current waveforms used in prosthetic devices for restoring partial vision lost to retinal degenerative diseases. The model framework combines a connectome-based neural network model characterized by accurate morphological and synaptic properties with an admittance method model of bulk tissue and prosthetic electronics. In this model, the retina was computationally "degenerated," considering cellular death and anatomical changes that occur early in disease, as well as altered neural behavior that develops throughout the neurodegeneration and is likely interfering with current attempts at restoring vision. A resulting analysis of stimulation range and threshold of ON ganglion cells within the retina that are either healthy or in beginning stages of degeneration is presented for currently used stimulation waveforms, and an asymmetric biphasic current stimulation for subduing spontaneous firing to allow increased control over ganglion cell firing patterns in degenerated retina is proposed. Results show that stimulation thresholds of retinal ganglion cells do not notably vary after beginning stages of retina degeneration. In addition, simulation of proposed asymmetric waveforms showed the ability to enhance the control of ganglion cell firing via electrical stimulation.

Minimizing activation of overlying axons with epiretinal stimulation: The role of fiber orientation and electrode configuration

Currently, a challenge in electrical stimulation of the retina with a visual prosthesis (bionic eye) is to excite only the cells lying directly under the electrode in the ganglion cell layer, while avoiding excitation of axon bundles that pass over the surface of the retina in the nerve fiber layer. Stimulation of overlying axons results in irregular visual percepts, limiting perceptual efficacy. This research explores how differences in fiber orientation between the nerve fiber layer and ganglion cell layer leads to differences in the electrical activation of the axon initial segment and axons of passage. Approach. Axons of passage of retinal ganglion cells in the nerve fiber layer are characterized by a narrow distribution of fiber orientations, causing highly anisotropic spread of applied current. In contrast, proximal axons in the ganglion cell layer have a wider distribution of orientations. A four-layer computational model of epiretinal extracellular stimulation that captures the effect of neurite orientation in anisotropic tissue has been developed using a volume conductor model known as the cellular composite model. Simulations are conducted to investigate the interaction of neural tissue orientation, stimulating electrode configuration, and stimulation pulse duration and amplitude. Main results. Our model shows that simultaneous stimulation with multiple electrodes aligned with the nerve fiber layer can be used to achieve selective activation of axon initial segments rather than passing fibers. This result can be achieved while reducing required stimulus charge density and with only modest increases in the spread of activation in the ganglion cell layer, and is shown to extend to the general case of arbitrary electrode array positioning and arbitrary target volume. Significance. These results elucidate a strategy for more targeted stimulation of retinal ganglion cells with experimentally-relevant multi-electrode geometries and achievable stimulation requirements.

Activation of ganglion cells and axon bundles using epiretinal electrical stimulation

Epiretinal prostheses for treating blindness activate axon bundles, causing large, arc-shaped visual percepts that limit the quality of artificial vision. Improving the function of epiretinal prostheses therefore requires understanding and avoiding axon bundle activation. This study introduces a method to detect axon bundle activation on the basis of its electrical signature and uses the method to test whether epiretinal stimulation can directly elicit spikes in individual retinal ganglion cells without activating nearby axon bundles. Combined electrical stimulation and recording from isolated primate retina were performed using a custom multielectrode system (512 electrodes, 10-μm diameter, 60-μm pitch). Axon bundle signals were identified by their bidirectional propagation, speed, and increasing amplitude as a function of stimulation current. The threshold for bundle activation varied across electrodes and retinas, and was in the same range as the threshold for activating retinal ganglion cells near their somas. In the peripheral retina, 45% of electrodes that activated individual ganglion cells (17% of all electrodes) did so without activating bundles. This permitted selective activation of 21% of recorded ganglion cells (7% of expected ganglion cells) over the array. In one recording in the central retina, 75% of electrodes that activated individual ganglion cells (16% of all electrodes) did so without activating bundles. The ability to selectively activate a subset of retinal ganglion cells without axon bundles suggests a possible novel architecture for future epiretinal prostheses.NEW & NOTEWORTHY Large-scale multielectrode recording and stimulation were used to test how selectively retinal ganglion cells can be electrically activated without activating axon bundles. A novel method was developed to identify axon activation on the basis of its unique electrical signature and was used to find that a subset of ganglion cells can be activated at single-cell, single-spike resolution without producing bundle activity in peripheral and central retina. These findings have implications for the development of advanced retinal prostheses.

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.

On the computation of a retina resistivity profile for applications in multi-scale modeling of electrical stimulation and absorption

This study proposes a methodology for computationally estimating resistive properties of tissue in multi-scale computational models, used for studying the interaction of electromagnetic fields with neural tissue, with applications to both dosimetry and neuroprosthetics. Traditionally, models at bulk tissue- and cellular-level scales are solved independently, linking resulting voltage from existing resistive tissue-scale models as extracellular sources to cellular models. This allows for solving the effects that external electric fields have on cellular activity. There are two major limitations to this approach: first, the resistive properties of the tissue need to be chosen, of which there are contradicting measurements in literature; second, the measurements of resistivity themselves may be inaccurate, leading to the mentioned contradicting results found across different studies. Our proposed methodology allows for constructing computed resistivity profiles using knowledge of only the neural morphology within the multi-scale model, resulting in a practical implementation of the effective medium theory; this bypasses concerns regarding the choice of resistive properties and accuracy of measurement setups. A multi-scale model of retina is constructed with an external electrode to serve as a test bench for analyzing existing and resulting resistivity profiles, and validation is presented through the reconstruction of a published resistivity profile of retina tissue. Results include a computed resistivity profile of retina tissue for use with a retina multi-scale model used to analyze effects of external electric fields on neural activity.