A dual-process integrator-resonator model of the electrically stimulated human auditory nerve

What Is the Benefit of Ramped Pulse Shapes for Activating Auditory Cortex Neurons? An Electrophysiological Study in an Animal Model of Cochlear Implant

In all commercial cochlear implant (CI) devices, the activation of auditory nerve fibers is performed with rectangular pulses that have two phases of opposite polarity. Recently, several papers proposed that ramped pulse shapes could be an alternative shape for efficiently activating auditory nerve fibers. Here, we investigate whether ramped pulse shapes can activate auditory cortex (ACx) neurons in a more efficient way than the classical rectangular pulses. Guinea pigs were implanted with CI devices and responses of ACx neurons were tested with rectangular pulses and with four ramped pulse shapes, with a first-phase being either cathodic or anodic. The thresholds, i.e., the charge level necessary for obtaining significant cortical responses, were almost systematically lower with ramped pulses than with rectangular pulses. The maximal firing rate (FR) elicited by the ramped pulses was higher than with rectangular pulses. As the maximal FR occurred with lower charge levels, the dynamic range (between threshold and the maximal FR) was not modified. These effects were obtained with cathodic and anodic ramped pulses. By reducing the charge levels required to activate ACx neurons, the ramped pulse shapes should reduce charge consumption and should contribute to more battery-efficient CI devices in the future.

Temporal integration of short-duration pulse trains in cochlear implant listeners: Psychophysical and electrophysiological measurements

While electrically-evoked auditory brainstem response (eABR) thresholds for low-rate pulse trains correlate well with behavioral thresholds measured at the same rate, the correlation is much weaker with behavioral thresholds measured at high rates, such as used clinically. This implies that eABRs to low-rate stimuli cannot be reliably used for objective programming of threshold levels in cochlear implant (CI) users. Here, we investigate whether the use of bunched-up pulses (BUPS), consisting of groups of closely-spaced pulses may be used as an alternative stimulus. Experiment 1 measured psychophysical detection thresholds for several stimuli having a period of 32 ms in nine CI subjects implanted with a Med-EL device. The stimuli differed in the number of pulses present in each period (from 1 to 32), the pulse rate within period (1000 pps and as high as possible for BUPS) and the electrode location (apical or basal). The correlation between psychophysical thresholds obtained for a high-rate (1000 pps) clinical stimulus and for the BUPS stimuli increased as the number of pulses per period of BUPS increased from 1 to 32. This first psychophysical experiment suggests that the temporal processes affecting the threshold of clinical stimuli are also present for BUPS. Experiment 2 measured eABRs on the apical electrode of eight CI subjects for BUPS having 1, 2, 4, 8, 16 or 32 pulses per period. For most subjects, wave V was visible for BUPS having up to 16 pulses per period. The latency of wave V at threshold increased as a function of the number of pulses per period, suggesting that the eABR reflects the integration of multiple pulses at such low levels or that the neural response to each individual pulse increases along the sequence due to facilitation processes. There was also a strong within-subject correlation between electrophysiological and behavioral thresholds for the different BUPS stimuli. This demonstrates that the drop in behavioral threshold obtained when increasing the number of pulses per period of the BUPS can be measured electrophysiologically using eABRs. In contrast, the across-subject correlation between eABR thresholds for BUPS and clinical thresholds remained relatively weak and did not increase with the number of pulses per period. Implications of the use of BUPS for objective programming of CIs are discussed.

Effect of neural adaptation and degeneration on pulse-train ECAPs: A model study

Electrically evoked compound action potentials (eCAPs) are measurements of the auditory nerve's response to electrical stimulation. ECAP amplitudes during pulse trains can exhibit temporal alternations. The magnitude of this alternation tends to diminish over time during the stimulus. How this pattern relates to the temporal behavior of nerve fibers is not known. We hypothesized that the stochasticity, refractoriness, adaptation of the threshold and spike-times influence pulse-train eCAP responses. Thirty thousand auditory nerve fibers were modeled in a three-dimensional cochlear model incorporating pulse-shape effects, pulse-history effects, and stochasticity in the individual neural responses. ECAPs in response to pulse trains of different rates and amplitudes were modeled for fibers with different stochastic properties (by variation of the relative spread) and different temporal properties (by variation of the refractory periods, adaptation and latency). The model predicts alternation of peak amplitudes similar to available human data. In addition, the peak alternation was affected by changing the refractoriness, adaptation, and relative spread of auditory nerve fibers. As these parameters are related to factors such as the duration of deafness and neural survival, this study suggests that the eCAP pattern in response to pulse trains could be used to assess the underlying temporal and stochastic behavior of the auditory nerve. As these properties affect the nerve's response to pulse trains, they are of uttermost importance to sound perception with cochlear implants.

Effects of the relative timing of opposite-polarity pulses on loudness for cochlear implant listeners

The symmetric biphasic pulses used in contemporary cochlear implants (CIs) consist of both cathodic and anodic currents, which may stimulate different sites on spiral ganglion neurons and, potentially, interact with each other. The effect on the order of anodic and cathodic stimulation on loudness at short inter-pulse intervals (IPIs; 0-800 μs) is investigated. Pairs of opposite-polarity pseudomonophasic (PS) pulses were used and the amplitude of each pulse was manipulated independently. In experiment 1 the two PS pulses differed in their current level in order to elicit the same loudness when presented separately. Six users of the Advanced Bionics CI (Valencia, CA) loudness-ranked trains of the pulse pairs using a midpoint-comparison procedure. Stimuli with anodic-leading polarity were louder than those with cathodic-leading polarity for IPIs shorter than 400 μs. This effect was small-about 0.3 dB-but consistent across listeners. When the same procedure was repeated with both PS pulses having the same current level (experiment 2), anodic-leading stimuli were still louder than cathodic-leading stimuli at very short intervals. However, when using symmetric biphasic pulses (experiment 3) the effect disappeared at short intervals and reversed at long intervals. Possible peripheral sources of such polarity interactions are discussed.

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.

Modeled auditory nerve responses to amplitude modulated cochlear implant stimulation

Cochlear implants encode speech information by stimulating the auditory nerve with amplitude-modulated pulse trains. A computer model of the auditory nerve's response to electrical stimulation can be used to evaluate different approaches to improving CI patients' perception. In this paper a computationally efficient stochastic and adaptive auditory nerve model was used to investigate full nerve responses to amplitude-modulated electrical pulse trains. The model was validated for nerve responses to AM pulse trains via comparison with animal data. The influence of different parameters, such as adaptation and stochasticity, on long-term adaptation and modulation-following behavior was investigated. Responses to pulse trains with different pulse amplitudes, amplitude modulation frequencies, and modulation depths were modeled. Rate responses as well as period histograms, Vector Strength and the fundamental frequency were characterized in different time bins. The response alterations, including frequency following behavior, observed over the stimulus duration were similar to those seen in animal experiments. The tested model can be used to predict complete nerve responses to arbitrary input, and thus to different sound coding strategies.

Effect of Pulse Polarity on Thresholds and on Non-monotonic Loudness Growth in Cochlear Implant Users

Most cochlear implants (CIs) activate their electrodes non-simultaneously in order to eliminate electrical field interactions. However, the membrane of auditory nerve fibers needs time to return to its resting state, causing the probability of firing to a pulse to be affected by previous pulses. Here, we provide new evidence on the effect of pulse polarity and current level on these interactions. In experiment 1, detection thresholds and most comfortable levels (MCLs) were measured in CI users for 100-Hz pulse trains consisting of two consecutive biphasic pulses of the same or of opposite polarity. All combinations of polarities were studied: anodic-cathodic-anodic-cathodic (ACAC), CACA, ACCA, and CAAC. Thresholds were lower when the adjacent phases of the two pulses had the same polarity (ACCA and CAAC) than when they were different (ACAC and CACA). Some subjects showed a lower threshold for ACCA than for CAAC while others showed the opposite trend demonstrating that polarity sensitivity at threshold is genuine and subject- or electrode-dependent. In contrast, anodic (CAAC) pulses always showed a lower MCL than cathodic (ACCA) pulses, confirming previous reports. In experiments 2 and 3, the subjects compared the loudness of several pulse trains differing in current level separately for ACCA and CAAC. For 40 % of the electrodes tested, loudness grew non-monotonically as a function of current level for ACCA but never for CAAC. This finding may relate to a conduction block of the action potentials along the fibers induced by a strong hyperpolarization of their central processes. Further analysis showed that the electrodes showing a lower threshold for ACCA than for CAAC were more likely to yield a non-monotonic loudness growth. It is proposed that polarity sensitivity at threshold reflects the local neural health and that anodic asymmetric pulses should preferably be used to convey sound information while avoiding abnormal loudness percepts.

A Model of Electrically Stimulated Auditory Nerve Fiber Responses with Peripheral and Central Sites of Spike Generation

A computational model of cat auditory nerve fiber (ANF) responses to electrical stimulation is presented. The model assumes that (1) there exist at least two sites of spike generation along the ANF and (2) both an anodic (positive) and a cathodic (negative) charge in isolation can evoke a spike. A single ANF is modeled as a network of two exponential integrate-and-fire point-neuron models, referred to as peripheral and central axons of the ANF. The peripheral axon is excited by the cathodic charge, inhibited by the anodic charge, and exhibits longer spike latencies than the central axon; the central axon is excited by the anodic charge, inhibited by the cathodic charge, and exhibits shorter spike latencies than the peripheral axon. The model also includes subthreshold and suprathreshold adaptive feedback loops which continuously modify the membrane potential and can account for effects of facilitation, accommodation, refractoriness, and spike-rate adaptation in ANF. Although the model is parameterized using data for either single or paired pulse stimulation with monophasic rectangular pulses, it correctly predicts effects of various stimulus pulse shapes, stimulation pulse rates, and level on the neural response statistics. The model may serve as a framework to explore the effects of different stimulus parameters on psychophysical performance measured in cochlear implant listeners.

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.

Phenomenological modelling of electrically stimulated auditory nerve fibers: A review

Auditory nerve fibers (ANFs) play a crucial role in hearing by encoding and transporting the synaptic input from inner hair cells into afferent spiking information for higher stages of the auditory system. If the inner hair cells are degenerated, cochlear implants may restore hearing by directly stimulating the ANFs. The response of an ANF is affected by several characteristics of the electrical stimulus and of the ANF, and neurophysiological measurements are needed to know how the ANF responds to a particular stimulus. However, recording from individual nerve fibers in humans is not feasible and obtaining compound neural or psychophysical responses is often time-consuming. This motivates the design and use of models to estimate the ANF response to the electrical stimulation. Phenomenological models reproduce the ANF response based on a simplified description of ANF functionality and on a limited parameter space by not directly describing detailed biophysical mechanisms. Here, we give an overview of phenomenological models published to date and demonstrate how different modeling approaches can account for the diverse phenomena affecting the ANF response. To highlight the success achieved in designing such models, we also describe a number of applications of phenomenological models to predict percepts of cochlear implant listeners.

A Phenomenological Model of the Electrically Stimulated Auditory Nerve Fiber: Temporal and Biphasic Response Properties

We present a phenomenological model of electrically stimulated auditory nerve fibers (ANFs). The model reproduces the probabilistic and temporal properties of the ANF response to both monophasic and biphasic stimuli, in isolation. The main contribution of the model lies in its ability to reproduce statistics of the ANF response (mean latency, jitter, and firing probability) under both monophasic and cathodic-anodic biphasic stimulation, without changing the model's parameters. The response statistics of the model depend on stimulus level and duration of the stimulating pulse, reproducing trends observed in the ANF. In the case of biphasic stimulation, the model reproduces the effects of pseudomonophasic pulse shapes and also the dependence on the interphase gap (IPG) of the stimulus pulse, an effect that is quantitatively reproduced. The model is fitted to ANF data using a procedure that uniquely determines each model parameter. It is thus possible to rapidly parameterize a large population of neurons to reproduce a given set of response statistic distributions. Our work extends the stochastic leaky integrate and fire (SLIF) neuron, a well-studied phenomenological model of the electrically stimulated neuron. We extend the SLIF neuron so as to produce a realistic latency distribution by delaying the moment of spiking. During this delay, spiking may be abolished by anodic current. By this means, the probability of the model neuron responding to a stimulus is reduced when a trailing phase of opposite polarity is introduced. By introducing a minimum wait period that must elapse before a spike may be emitted, the model is able to reproduce the differences in the threshold level observed in the ANF for monophasic and biphasic stimuli. Thus, the ANF response to a large variety of pulse shapes are reproduced correctly by this model.

A point process framework for modeling electrical stimulation of the auditory nerve

Model-based studies of responses of auditory nerve fibers to electrical stimulation can provide insight into the functioning of cochlear implants. Ideally, these studies can identify limitations in sound processing strategies and lead to improved methods for providing sound information to cochlear implant users. To accomplish this, models must accurately describe spiking activity while avoiding excessive complexity that would preclude large-scale simulations of populations of auditory nerve fibers and obscure insight into the mechanisms that influence neural encoding of sound information. In this spirit, we develop a point process model of individual auditory nerve fibers that provides a compact and accurate description of neural responses to electric stimulation. Inspired by the framework of generalized linear models, the proposed model consists of a cascade of linear and nonlinear stages. We show how each of these stages can be associated with biophysical mechanisms and related to models of neuronal dynamics. Moreover, we derive a semianalytical procedure that uniquely determines each parameter in the model on the basis of fundamental statistics from recordings of single fiber responses to electric stimulation, including threshold, relative spread, jitter, and chronaxie. The model also accounts for refractory and summation effects that influence the responses of auditory nerve fibers to high pulse rate stimulation. Throughout, we compare model predictions to published physiological data of response to high and low pulse rate stimulation. We find that the model, although constructed to fit data from single and paired pulse experiments, can accurately predict responses to unmodulated and modulated pulse train stimuli. We close by performing an ideal observer analysis of simulated spike trains in response to sinusoidally amplitude modulated stimuli and find that carrier pulse rate does not affect modulation detection thresholds.

Bistability and resonance in the periodically stimulated Hodgkin-Huxley model with noise

We describe general characteristics of the Hodgkin-Huxley neuron's response to a periodic train of short current pulses with Gaussian noise. The deterministic neuron is bistable for antiresonant frequencies. When the stimuli arrive at the resonant frequency the firing rate is a continuous function of the current amplitude I(0) and scales as (I(0)-I(th))(1/2), characteristic of a saddle-node bifurcation at the threshold I(th). Intervals of continuous irregular response alternate with integer mode-locked regions with bistable excitation edge. There is an even-all multimodal transition between the 2 : 1 and 3 : 1 states in the vicinity of the main resonance, which is analogous to the odd-all transition discovered earlier in the high-frequency regime. For I(0)<I(th) and small noise the firing rate has a maximum at the resonant frequency. For larger noise and subthreshold stimulation the maximum firing rate initially shifts toward lower frequencies, then returns to higher frequencies in the limit of large noise. The stochastic coherence antiresonance, defined as a simultaneous occurrence of (i) the maximum of the coefficient of variation and (ii) the minimum of the firing rate vs the noise intensity, occurs over a wide range of parameter values, including monostable regions. Results of this work can be verified experimentally.

A multi-timescale adaptive threshold model for the SAI tactile afferent to predict response to mechanical vibration

The Leaky Integrate and Fire (LIF) model of a neuron is one of the best known models for a spiking neuron. A current limitation of the LIF model is that it may not accurately reproduce the dynamics of an action potential. There have recently been some studies suggesting that a LIF coupled with a multi-timescale adaptive threshold (MAT) may increase LIF's accuracy in predicting spikes in cortical neurons. We propose a mechanotransduction process coupled with a LIF model with multi-timescale adaptive threshold to model slowly adapting type I (SAI) mechanoreceptor in monkey's glabrous skin. In order to test the performance of the model, the spike timings predicted by this MAT model are compared with neural data. We also test a fixed threshold variant of the model by comparing its outcome with the neural data. Initial results indicate that the MAT model predicts spike timings better than a fixed threshold LIF model only.

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.

Modelled temperature-dependent excitability behaviour of a generalised human peripheral sensory nerve fibre

The objective of this study was to determine if a recently developed human Ranvier node model, which is based on a modified version of the Hodgkin-Huxley model, could predict the excitability behaviour in human peripheral sensory nerve fibres with diameters ranging from 5.0 to 15.0 microm. The Ranvier node model was extended to include a persistent sodium current and was incorporated into a generalised single cable nerve fibre model. Parameter temperature dependence was included. All calculations were performed in Matlab. Sensory nerve fibre excitability behaviour characteristics predicted by the new nerve fibre model at different temperatures and fibre diameters compared well with measured data. Absolute refractory periods deviated from measured data, while relative refractory periods were similar to measured data. Conduction velocities showed both fibre diameter and temperature dependence and were underestimated in fibres thinner than 12.5 microm. Calculated strength-duration time constants ranged from 128.5 to 183.0 micros at 37 degrees C over the studied nerve fibre diameter range, with chronaxie times about 30% shorter than strength-duration time constants. Chronaxie times exhibited temperature dependence, with values overestimated by a factor 5 at temperatures lower than body temperature. Possible explanations include the deviated absolute refractory period trend and inclusion of a nodal strangulation relationship.

Alternative pulse shapes in electrical hearing

Cochlear implants (CIs) stimulate the auditory nerve with trains of symmetric biphasic (BI) pulses. We review studies showing that more efficient stimulation can be achieved by modifying these pulses by (1) increasing the inter-phase gap (IPG) between the two phases of each pulse, thereby delaying the recovery of charge, (2) increasing the duration and decreasing the amplitude of one phase - so-called "pseudomonophasic (PS)" waveforms, and (3) combining the pseudomonophasic stimulus with an IPG in a "delayed pseudomonophasic" waveform (PS_IPG). These efficiency gains, measured using changes in threshold and loudness, occur at a wide range of pulse rates, including those commonly used in current CI systems. In monopolar mode, dynamic ranges are larger for PS and for long-IPG pulse shapes than for BI pulses, but this increase in DR is not accompanied by a higher number of discriminable loudness steps, and hence, in a better coding of loudness. Moreover, waveforms with relatively short and long interphase gaps do not yield different patterns of excitation despite the relatively large differences in threshold. Two important findings are that, contrary to data obtained in animal experiments, anodic currents are more effective than cathodic stimulation for human CI patients and that the thresholds decrease with increases in IPG over a much longer time course (more than 3 ms) than for animals. In this review it is discussed how these alternative pulse shapes may be beneficial in terms of reducing power consumption and channel interactions, which issues remain to be addressed, and how models contribute to guiding our research.