Modeling extracellular stimulation of retinal ganglion cells: theoretical and practical aspects

Avoidance of axonal stimulation with sinusoidal epiretinal stimulation

Objective.Neuromodulation, particularly electrical stimulation, necessitates high spatial resolution to achieve artificial vision with high acuity. In epiretinal implants, this is hindered by the undesired activation of distal axons. Here, we investigate focal and axonal activation of retinal ganglion cells (RGCs) in epiretinal configuration for different sinusoidal stimulation frequencies.Approach.RGC responses to epiretinal sinusoidal stimulation at frequencies between 40 and 100 Hz were tested inex-vivophotoreceptor degenerated (rd10) isolated retinae. Experiments were conducted using a high-density CMOS-based microelectrode array, which allows to localize RGC cell bodies and axons at high spatial resolution.Main results.We report current and charge density thresholds for focal and distal axon activation at stimulation frequencies of 40, 60, 80, and 100 Hz for an electrode size with an effective area of 0.01 mm2. Activation of distal axons is avoided up to a stimulation amplitude of 0.23µA (corresponding to 17.3µC cm-2) at 40 Hz and up to a stimulation amplitude of 0.28µA (14.8µC cm-2) at 60 Hz. The threshold ratio between focal and axonal activation increases from 1.1 for 100 Hz up to 1.6 for 60 Hz, while at 40 Hz stimulation frequency, almost no axonal responses were detected in the tested intensity range. With the use of synaptic blockers, we demonstrate the underlying direct activation mechanism of the ganglion cells. Finally, using high-resolution electrical imaging and label-free electrophysiological axon tracking, we demonstrate the extent of activation in axon bundles.Significance.Our results can be exploited to define a spatially selective stimulation strategy avoiding axonal activation in future retinal implants, thereby solving one of the major limitations of artificial vision. The results may be extended to other fields of neuroprosthetics to achieve selective focal electrical stimulation.

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.

Effects on Retinal Stimulation of the Geometry and the Insertion Location of Penetrating Electrodes

Retinal implants have been developed and implanted to restore vision from outer retinal degeneration, but their performance is still limited due to the poor spatial resolution. To improve the localization of stimulation, microelectrodes in various three-dimensional (3D) shapes have been investigated. In particular, computational simulation is crucial for optimizing the performance of a novel microelectrode design before actual fabrication. However, most previous studies have assumed a uniform conductivity for the entire retina without testing the effect of electrodes placement in different layers. In this study, we used the finite element method to simulate electric fields created by 3D microelectrodes of three different designs in a retina model with a stratified conductivity profile. The three electrode designs included two conventional shapes - a conical electrode (CE) and a pillar electrode (PE); we also proposed a novel structure of pillar electrode with an insulating wall (PEIW). A quantitative comparison of these designs shows the PEIW generates a stronger and more confined electric field with the same current injection, which is preferred for high-resolution retinal prostheses. Moreover, our results demonstrate both the magnitude and the shape of potential distribution generated by a penetrating electrode depend not only on the geometry, but also substantially on the insertion depth of the electrode. Although epiretinal insertions are mainly discussed, we also compared results for subretinal insertions. The results provide valuable insights for improving the spatial resolution of retinal implants using 3D penetrating microelectrodes and highlight the importance of considering the heterogeneity of conductivities in the retina.