Tania Hanekom

Three-dimensional spiraling finite element model of the electrically stimulated cochlea.

Objective The objective of the article is to provide an accurate model of the human cochlea with which potential distributions and thus neural excitation patterns around cochlear implant electrodes can be determined. Improvements on previous models of the implanted cochlea are that this model 1) includes the spiral nature of the cochlea as well as many other anatomical details (and it is a model of the human cochlear rather than the guinea pig cochlea), and 2) facilitates modeling of different electrode geometries, array locations and electrode separations without changing the structure of the model. Design A three-dimensional spiraling finite element model of the human cochlea was created. The model incorporates the effect of neighboring canals and conduction along the fluid-filled canals of the cochlea. Potential distributions are used as inputs to a nerve fiber model to investigate auditory nerve excitation patterns around intracochlear electrode arrays. Results Potential distributions around intracochlear electrodes generated with the finite element model are presented. The effects of electrode separation, electrode geometry and array location on excitation threshold, excitation spread and ectopic excitation (i.e., excitation of nerve fibers at an undesirable location) are demonstrated. Conclusions The following conclusions should be considered preliminary, as their accuracy depends on the exactness of the underlying model. The spiraling geometry of the cochlea causes asymmetry in potential distributions. The location of electrodes along the length of the basilar membrane has a stronger influence on the site of excitation than the polarity of the leading phase of the stimulus. Array location is the primary parameter that controls excitation spread. Threshold currents and the effect of ongoing loss of peripheral dendrites on threshold currents can be limited by placing arrays close to the modiolus. Point electrode geometries are recommended above banded electrode geometries only when the array can be placed close to the modiolus. There is a tradeoff between array location and the degree of ectopic stimulation caused by a specific array location. Bimodal excitation patterns exist at comfortable stimulus intensities for longitudinal bipolar electrode configurations. It is shown that an electrode configuration with an electrode separation of approximately half that of the bipolar electrode separation of the Nucleus electrode can be used instead of radial and offset radial electrode configurations to create unimodal excitation patterns. The stimulation resolution of cochlear implant electrode arrays can potentially be improved by increasing the number of electrode contacts in an array.

Modelling encapsulation tissue around cochlear implant electrodes

The objective of the study was to explore the effect of electrode encapsulation by fibrous scar tissue on electrical potential distributions and auditory nerve fibre excitation patterns. A finite element model in combination with an auditory nerve fibre model was used to predict changes in threshold currents and auditory nerve fibre excitation patterns. The model showed that electrical potentials at the target nerve fibres and the electrode contacts changed in the presence of encapsulation tissue. This led to changes in threshold currents and spread of excitation. The effect of electrode encapsulation on threshold currents and spread of excitation depended on the thickness of the perilymph layer separating the fibrous tissue encapsulation and the electrode array, nerve fibre survival status, electrode geometry and configuration, and array location. Model results suggested that arrays located close to the modiolus were most sensitive to threshold changes caused by electrode encapsulation (changes were between −0.26 and 2.41 dB), whereas encapsulation of an electrode array had less effect on threshold currents when the array was located in a lateral position in the scala tympani (changes were between −0.64 and 1.5 dB). For medially located arrays, changes in the spread of excitation varied between an increase of 0.21 mm and a decrease of 0.33 mm along the length of the basilar membrane, and an increase of 0.18 mm and a decrease of 0.66 mm along the length of the basilar membrane were calculated for laterally located arrays.

A finite element model for describing the effect of muscle shortening on surface EMG

A finite-element model for the generation of single fiber action potentials in a muscle undergoing various degrees of fiber shortening is developed. The muscle is assumed fusiform with muscle fibers following a curvilinear path described by a Gaussian function. Different degrees of fiber shortening are simulated by changing the parameters of the fiber path and maintaining the volume of the muscle constant. The conductivity tensor is adapted to the muscle fiber orientation. In each point of the volume conductor, the conductivity of the muscle tissue in the direction of the fiber is larger than that in the transversal direction. Thus, the conductivity tensor changes point-by-point with fiber shortening, adapting to the fiber paths. An analytical derivation of the conductivity tensor is provided. The volume conductor is then studied with a finite-element approach using the analytically derived conductivity tensor. Representative simulations of single fiber action potentials with the muscle at different degrees of shortening are presented. It is shown that the geometrical changes in the muscle, which imply changes in the conductivity tensor, determine important variations in action potential shape, thus affecting its amplitude and frequency content. The model provides a new tool for interpreting surface EMG signal features with changes in muscle geometry, as it happens during dynamic contractions.

Effect of spatial filtering on crosstalk reduction in surface EMG recordings

Increasing the selectivity of the detection system in surface electromyography (EMG) is beneficial in the collection of information of a specific portion of the investigated muscle and to reduce the contribution of undesired components, such as non-propagating components (due to generation or end-of-fibre effects) or crosstalk from nearby muscles. A comparison of the ability of different spatial filters to reduce the amount of crosstalk in surface EMG measurements was conducted in this paper using simulated signals. It focused on the influence of different properties of the muscle anatomy (changing subcutaneous layer thickness, skin conductivity, fibre length) and detection system (single, double and normal double differential, with two inter-electrode distances - IED) on the amount of crosstalk present in the measurements. A cylindrical multilayer (skin, subcutaneous tissue, muscle, bone) analytical model was used to simulate single fibre action potentials (SFAPs). Fibres were grouped together in motor units (MUs) and motor unit action potentials (MUAPs) were obtained by adding the SFAPs of the corresponding fibres. Interference surface EMG signals were obtained, modelling the recruitment of MUs and rate coding. The average rectified value (ARV) and mean frequency (MNF) content of the EMG signals were studied and used as a basis for determining the selectivity of each spatial filter. From these results it was found that the selectivity of each spatial filter varies depending on the transversal location of the measurement electrodes and on the anatomy. An increase in skin conductivity favourably affects the selectivity of normal double differential filters as does an increase in subcutaneous layer thickness. An increase in IED decreases the selectivity of all the analysed filters.

Can subject-specific single-fibre electrically evoked auditory brainstem response data be predicted from a model?

This article investigates whether prediction of subject-specific physiological data is viable through an individualised computational model of a cochlear implant. Subject-specific predictions could be particularly useful to assess and quantify the peripheral factors that cause inter-subject variations in perception. The results of such model predictions could potentially be translated to clinical application through optimisation of mapping parameters for individual users, since parameters that affect perception would be reflected in the model structure and parameters. A method to create a subject-specific computational model of a guinea pig with a cochlear implant is presented. The objectives of the study are to develop a method to construct subject-specific models considering translation of the method to in vivo human models and to assess the effectiveness of subject-specific models to predict peripheral neural excitation on subject level. Neural excitation patterns predicted by the model are compared with single-fibre electrically evoked auditory brainstem responses obtained from the inferior colliculus in the same animal. Results indicate that the model can predict threshold frequency location, spatial spread of bipolar and tripolar stimulation and electrode thresholds relative to one another where electrodes are located in different cochlear structures. Absolute thresholds and spatial spread using monopolar stimulation are not predicted accurately. Improvements to the model should address this.

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.

Estimation of stimulus attenuation in cochlear implants

Neural excitation profile widths at the neural level, for monopolar stimulation with Nucleus straight and contour arrays respectively, were simulated using a combined volume-conduction-neural model. The electrically evoked compound action potential profile widths at the electrode array level were calculated with a simple approximation method employing stimulus attenuation inside the cochlear duct, and the results compared to profile width data from literature. The objective of the article is to develop a simple method to estimate stimulus attenuation values by calculating the values that best fit the modelled excitation profile widths to the measured evoked compound action potential profile widths. Results indicate that the modelled excitation profile widths decrease with increasing stimulus attenuation. However, fitting of modelled excitation profile widths to measured evoked compound action potential profile widths show that different stimulus attenuation values are needed for different stimulation levels. It is suggested that the proposed simple model can provide an estimate of stimulus attenuation by calculating the value of the parameter that produces the best fit to experimental data in specific human subjects.

Constructing a three-dimensional electrical model of a living cochlear implant user's cochlea.

BACKGROUND Hearing performance varies greatly among users of cochlear implants. Current three-dimensional cochlear models that predict the electrical fields inside a stimulated cochlea and their effect on neural excitation are generally based on a generic human or guinea pig cochlear shape that does not take inter-user morphological variations into account. This precludes prediction of user-specific performance. AIMS The aim of this study is to develop a model of the implanted cochlea of a specific living human individual and to assess if the inclusion of morphological variations in cochlear models affects predicted outcomes significantly. METHODS Five three-dimensional electric volume conduction models of the implanted cochleae of individual living users were constructed from standard CT scan data. These models were embedded in head models that include monopolar return electrodes in accurate anatomic positions. Potential distributions and neural excitation patterns were predicted for each of the models. RESULTS Modeled potential distributions and neural excitation profiles (threshold amplitudes, center frequencies, and bandwidths) are affected by user-specific cochlear morphology and electrode placement within the cochlea. CONCLUSIONS This work suggests that the use of user-specific models is indicated when more detailed analysis is required than what is available from generic models. Copyright

Three-dimensional models of cochlear implants: A review of their development and how they could support management and maintenance of cochlear implant performance.

ABSTRACT Three-dimensional (3D) computational modeling of the auditory periphery forms an integral part of modern-day research in cochlear implants (CIs). These models consist of a volume conduction description of implanted stimulation electrodes and the current distribution around these, coupled with auditory nerve fiber models. Cochlear neural activation patterns can then be predicted for a given input stimulus. The objective of this article is to present the context of 3D modeling within the field of CIs, the different models, and approaches to models that have been developed over the years, as well as the applications and potential applications of these models. The process of development of 3D models is discussed, and the article places specific emphasis on the complementary roles of generic models and user-specific models, as the latter is important for translation of these models into clinical application.

Initial results from a model of ephaptic excitation in the electrically excited peripheral auditory nervous system.

This article reports on a study performed to investigate the occurrence and effect of ephaptic excitation in electrical stimulation of the auditory system. The objective of the study was to quantify the influence of ephaptic excitation on nerve stimulation and determine whether it is a necessary factor in neuromodeling. It was shown with a simple model that ephaptic excitation could be important at stimulus intensities close to threshold. The results show that the contribution of ephaptic excitation is significant up to at least 6-7dB above threshold. Cochlear implant patients normally have a small dynamic range (average of 7dB), indicating that the ephaptic effect might be important in models of the implanted cochlea.