Estimation of stimulus attenuation in cochlear implants

A Novel Algorithm to Derive Spread of Excitation Based on Deconvolution

OBJECTIVE

The width of the spread of excitation (SOE) curve has been widely thought to represent an estimate of SOE. Therefore, correlates between psychophysical parameters, such as pitch discrimination and speech perception, and the width of SOE curves, have long been investigated. However, to date, no relationships between these objective and subjective measurements have yet been determined. In a departure from the current thinking, the authors now propose that the SOE curve, recorded with forward masking, is the equivalent of a convolution operation. As such, deconvolution would be expected to retrieve the excitation areas attributable to either masker or probe, potentially more closely revealing the actual neural SOE. This study aimed to develop a new analytical tool with which to derive SOE using this principle.

DESIGN

Intraoperative SOE curve measurements of 16 subjects, implanted with an Advanced Bionics implant, were analyzed. Evoked compound action potential (ECAP)-based SOE curves were recorded on electrodes 3 to 16, using the forward masker paradigm, with variable masker. The measured SOE curves were then compared with predicted SOE curves, built by the convolution of basic excitation density profiles (EDPs). Predicted SOE curves were fitted to the measured SOEs by iterative adjustment of the EDPs for the masker and the probe.

RESULTS

It was possible to generate a good fit between the predicted and measured SOE curves, inclusive of their asymmetry. The rectangular EDP was of least value in terms of its ability to generate a good fit; smoother SOE curves were modeled using the exponential or Gaussian EDPs. In most subjects, the EDP width (i.e., the size of the excitation area) gradually changed from wide at the apex of the electrode array, to narrow at the base. A comparison of EDP widths to SOE curve widths, as calculated in the literature, revealed that the EDPs now provide a measure of the SOE that is qualitatively distinct from that provided using conventional methods.

CONCLUSIONS

This study shows that an eCAP-based SOE curve, measured with forward masking, can be treated as a convolution of EDPs for masker and probe. The poor fit achieved for the measured and modeled data using the rectangular EDP, emphasizes the requirement for a sloping excitation area to mimic actual SOE recordings. Our deconvolution method provides an explanation for the frequently observed asymmetry of SOE curves measured along the electrode array, as this is a consequence of a wider excitation area in the apical part of the cochlea, in the absence of any asymmetry in the actual EDP. In addition, broader apical EDPs underlie the higher eCAP amplitudes found for apical stimulation.

Analysis of a purely conductance-based stochastic nerve fibre model as applied to compound models of populations of human auditory nerve fibres used in cochlear implant simulations

The study presents the application of a purely conductance-based stochastic nerve fibre model to human auditory nerve fibres within finite element volume conduction models of a semi-generic head and user-specific cochleae. The stochastic, threshold and temporal characteristics of the human model are compared and successfully validated against physiological feline results with the application of a mono-polar, bi-phasic, cathodic first stimulus. Stochastic characteristics validated include: (i) the log(Relative Spread) versus log(fibre diameter) distribution for the discharge probability versus stimulus intensity plots and (ii) the required exponential membrane noise versus transmembrane voltage distribution. Intra-user, and to a lesser degree inter-user, comparisons are made with respect to threshold and dynamic range at short and long pulse widths for full versus degenerate single fibres as well as for populations of degenerate fibres of a single user having distributed and aligned somas with varying and equal diameters. Temporal characteristics validated through application of different stimulus pulse rates and different stimulus intensities include: (i) discharge rate, latency and latency standard deviation versus stimulus intensity, (ii) period histograms and (iii) interspike interval histograms. Although the stochastic population model does not reduce the modelled single deterministic fibre threshold, the simulated stochastic and temporal characteristics show that it could be used in future studies to model user-specific temporally encoded information, which influences the speech perception of CI users.

Development of a voltage-dependent current noise algorithm for conductance-based stochastic modelling of auditory nerve fibres

This study presents the development of an alternative noise current term and novel voltage-dependent current noise algorithm for conductance-based stochastic auditory nerve fibre (ANF) models. ANFs are known to have significant variance in threshold stimulus which affects temporal characteristics such as latency. This variance is primarily caused by the stochastic behaviour or microscopic fluctuations of the node of Ranvier's voltage-dependent sodium channels of which the intensity is a function of membrane voltage. Though easy to implement and low in computational cost, existing current noise models have two deficiencies: it is independent of membrane voltage, and it is unable to inherently determine the noise intensity required to produce in vivo measured discharge probability functions. The proposed algorithm overcomes these deficiencies while maintaining its low computational cost and ease of implementation compared to other conductance and Markovian-based stochastic models. The algorithm is applied to a Hodgkin-Huxley-based compartmental cat ANF model and validated via comparison of the threshold probability and latency distributions to measured cat ANF data. Simulation results show the algorithm's adherence to in vivo stochastic fibre characteristics such as an exponential relationship between the membrane noise and transmembrane voltage, a negative linear relationship between the log of the relative spread of the discharge probability and the log of the fibre diameter and a decrease in latency with an increase in stimulus intensity.

Stimulation strategies and electrode design in computational models of the electrically stimulated cochlea: An overview of existing literature

Since the 1970s, computational modeling has been used to investigate the fundamental mechanisms of cochlear implant stimulation. Lumped parameter models and analytical models have been used to simulate cochlear potentials, as well as three-dimensional volume conduction models based on the Finite Difference, Finite Element, and Boundary Element methods. Additionally, in order to simulate neural responses, several of these cochlear models have been combined with nerve models, which were either simple activation functions or active nerve fiber models of the cochlear auditory neurons. This review paper will present an overview of the ways in which these computational models have been employed to study different stimulation strategies and electrode designs. Research into stimulation strategies has concentrated mainly on multipolar stimulation as a means of achieving current focussing and current steering, while modeling work on electrode design has been chiefly concerned with finding the optimal position and insertion depth of the electrode array. Finally, the present and future of computational modeling of the electrically stimulated cochlea is discussed.

Electrical stimulation characteristics of denervated orbicularis oculi muscle

This research is to study the electrical stimulation characteristics of orbicularis oculi muscle and the characteristics of the mechanical contraction. We observed the stimulus current diffusion regularity and its relationship with mechanical contraction in the orbicularis oculi muscle using an electrode gathering line. Under different stimulus intensities of 2 or 4 mA, the closer the recording electrodes were to the stimulating electrode, the larger was the amplitude. When the recording electrode and stimulating electrode distance increased, the amplitude declined linearly with decreasing function. In addition, current conduction across the muscle fiber was studied. Under different stimulus intensities of 2 or 4 mA, it was found that the closer the recording electrodes were to the stimulating electrode, the larger was the amplitude. When the recording electrode and stimulating electrode distance increased, the amplitude declined linearly with decreasing function. The transverse current reached a maximum 4 mA range, and increasing the current intensity did not increase the propagation range. Under different stimulation intensities, the larger the stimulus intensity, the greater is the potential change and the faster is the attenuation. Longitudinal current, even in the range of 6 mm, can still record electrical activity. While a transverse current diffuser has a maximum range of 4 mm, increasing the current intensity does not increase the propagation range.

Effects of electrode position on spatiotemporal auditory nerve fiber responses: a 3D computational model study

A cochlear implant (CI) is an auditory prosthesis that enables hearing by providing electrical stimuli through an electrode array. It has been previously established that the electrode position can influence CI performance. Thus, electrode position should be considered in order to achieve better CI results. This paper describes how the electrode position influences the auditory nerve fiber (ANF) response to either a single pulse or low- (250 pulses/s) and high-rate (5,000 pulses/s) pulse-trains using a computational model. The field potential in the cochlea was calculated using a three-dimensional finite-element model, and the ANF response was simulated using a biophysical ANF model. The effects were evaluated in terms of the dynamic range, stochasticity, and spike excitation pattern. The relative spread, threshold, jitter, and initiated node were analyzed for single-pulse response; and the dynamic range, threshold, initiated node, and interspike interval were analyzed for pulse-train stimuli responses. Electrode position was found to significantly affect the spatiotemporal pattern of the ANF response, and this effect was significantly dependent on the stimulus rate. We believe that these modeling results can provide guidance regarding perimodiolar and lateral insertion of CIs in clinical settings and help understand CI performance.

Comparative studies on the implantation of nano platinum black and pure platinum electrodes in the rabbit orbicularis oculi muscle

OBJECTIVE

To study the interactive influence of implanted nano platinum black electrodes (as compared with pure platinum electrodes) on rabbit orbicularis oculi muscle morphology and function.

METHODS

The influence of the two types of electrode on the orbicularis oculi muscle was monitored in a rabbit model of facial paralysis. Biological electric current and exciting current were administered to biological tissue, and morphological and functional changes were identified. Changes in orbicularis oculi muscle contraction, electrode configuration and performance associated with long-term electrical stimulation were observed over 28 days of implantation.

RESULTS

The nano platinum black electrode was superior to the pure platinum electrode in the following aspects: morphology and functionality, electrical excitation function of the orbicularis oculi muscle (as assessed by electromyography), muscle contraction function and biological tissue changes. Furthermore, the nano platinum black electrode features had good stability.

CONCLUSION

Microelectrode surface modification with nano platinum black can effectively increase the microelectrode surface area and improve electrode performance, and is associated with good tissue compatibility.

Morphometric classification and spatial organization of spiral ganglion neurons in the human cochlea: consequences for single fiber response to electrical stimulation

The unique, unmyelinated perikarya of spiral ganglion cells (SGCs) in the human cochlea are often arranged in functional units covered by common satellite glial cells. This micro anatomical peculiarity presents a crucial barrier for an action potential (AP) travelling from the sensory receptors to the brain. Confocal microscopy was used to acquire systematically volumetric data on perikarya and corresponding nuclei in their full dimension along the cochlea of two individuals.

Four populations of SGCs within the human inner ear of two different specimens were identified using agglomerative hierarchical clustering, contrary to the present distinction of two groups of SGCs. Furthermore, we found evidence of a spatial arrangement of perikarya and their accordant nuclei along the cochlea spiral. In this arrangement, the most uniform sizes of cell bodies are located in the middle turn, which represents the majority of phonational frequencies.

Since single-cell recordings from other mammalians may not be representative to humans and human SGCs are not accessible for physiological measurements, computer simulation has been used to quantify the effect of varying soma size on single neuron response to electrical micro stimulation. Results show that temporal parameters of the spiking pattern are affected by the size of the cell body. Cathodic stimulation was found to induce stronger variations of spikes while also leading to the lowest thresholds and longest latencies. Therefore, anodic stimulation leads to a more uniform excitation profile among SGCs with different cell body size.

Threshold predictions of different pulse shapes using a human auditory nerve fibre model containing persistent sodium and slow potassium currents

The ability of a human auditory nerve fibre computational model to predict threshold differences for biphasic, pseudomonophasic and alternating monophasic waveforms was investigated. The effect of increasing the interphase gap, interpulse interval and pulse rate on thresholds was also simulated. Simulations were performed for both anodic-first and cathodic-first stimuli. Results indicated that the model correctly predicted threshold reductions for pseudomonophasic compared to biphasic waveforms, although reduction for alternating monophasic waveforms was underestimated. Threshold reductions were more pronounced for cathodic-first stimuli compared to anodic-first stimuli. Reversal of the phases in pseudomonophasic stimuli suggested a threshold reduction for anodic-first stimuli, but a threshold increase in cathodic-first stimuli. Inclusion of the persistent sodium and slow potassium currents in the model resulted in a reasonably accurate prediction of the non-monotonic threshold behaviour for pulse rates higher than 1000 pps. However, the model did not correctly predict the threshold changes observed for low pulse rate biphasic and alternating monophasic waveforms. It was suggested that these results could in part be explained by the difference in the refractory periods between real and simulated auditory nerve fibres, but also by the lack of representation of stochasticity observed in real auditory nerve fibres in our auditory nerve model.