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

Model-based studies of auditory nerve responses 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 auditory nerve spiking 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 the auditory nerve 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 semi-analytical 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 and explain differences in auditory nerve responses to high and low pulse rate stimulation. 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.Comment: 1 title page, 27 manuscript pages, 14 figures, 1 table, 1 appendi

Evaluation of motor neuron excitability by CMAP scanning with modulated current

It is important to have better evaluation and understanding of the motor neuron physiology, with the goal to early and objectively diagnose and treat patients with neurodegenerative pathologies. The Compound Muscle Action Potential (CMAP) scan is a non-invasive diagnosis technique for neurodegenerative pathologies, such as ALS, and enables a quick analysis of the muscle action potentials in response to motor nerve stimulation. This work aims to study the influence of different pulse modulated waveforms in peripheral nerve excitability by CMAP scan technique on healthy subjects. A total of 13 healthy subjects were submitted to the same test. The stimuli were applied in the medium nerve on the right wrist and electromyography signal collected on the Abductor Pollicis Brevis (APB) muscle surface on the right thumb. Stimulation was performed with an increasing intensities range from 4 to 30 mA, with varying steps, 3 stimuli per step. The procedure was repeated 4 times per subject, each repetition using a different single pulse stimulation waveform: monophasic square, monophasic triangular, monophasic quadratic and biphasic square. Results were retrieved from the averaging of the stimuli on each current intensity step. The square pulse needs less current intensity to generate the same response amplitude regarding the other waves and presents a more steep curve slope and this effect is gradually decreasing for the triangular and quadratic pulse,respectively, being the difference even more evident regarding the biphasic pulse. The control of the waveform stimulation pulse allows varying the stimulusresponse curve slope

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

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

Paired associative stimulation targeting the tibialis anterior muscle using either mono or biphasic transcranial magnetic stimulation

Paired associative stimulation (PAS) protocols induce plastic changes within the motor cortex. The objectives of this study were to investigate PAS effects targeting the tibialis anterior (TA) muscle using a biphasic transcranial magnetic stimulation (TMS) pulse form and, to determine whether a reduced intensity of this pulse would lead to significant changes as has been reported for hand muscles using a monophasic TMS pulse. Three interventions were investigated: (1) suprathreshold PAbi-PAS (n = 11); (2) suprathreshold PAmono-PAS (n = 11) where PAS was applied using a biphasic or monophasic pulse form at 120% resting motor threshold (RMT); (3) subthreshold PAbi-PAS (n = 10) where PAS was applied as for (1) at 95% active motor threshold (AMT). The peak-to-peak motor evoked potentials (MEPs) were quantified prior to, immediately following, and 30 min after the cessation of the intervention. TA MEP size increased significantly for all interventions immediately post (61% for suprathreshold PAbi-PAS, 83% for suprathreshold PAmono-PAS, 55% for subthreshold PAbi-PAS) and 30 min after the cessation of the intervention (123% for suprathreshold PAbi-PAS, 105% for suprathreshold PAmono-PAS, 80% for subthreshold PAbi-PAS. PAS using a biphasic pulse form at subthreshold intensities induces similar effects to conventional PAS

Comparison of Multi-Compartment Cable Models of Human Auditory Nerve Fibers

Background: Multi-compartment cable models of auditory nerve fibers have been developed to assist in the improvement of cochlear implants. With the advancement of computational technology and the results obtained from in vivo and in vitro experiments, these models have evolved to incorporate a considerable degree of morphological and physiological details. They have also been combined with three-dimensional volume conduction models of the cochlea to simulate neural responses to electrical stimulation. However, no specific rules have been provided on choosing the appropriate cable model, and most models adopted in recent studies were chosen without a specific reason or by inheritance. Methods: Three of the most cited biophysical multi-compartment cable models of the human auditory nerve, i.e., Rattay et al. (2001b), Briaire and Frijns (2005), and Smit et al. (2010), were implemented in this study. Several properties of single fibers were compared among the three models, including threshold, conduction velocity, action potential shape, latency, refractory properties, as well as stochastic and temporal behaviors. Experimental results regarding these properties were also included as a reference for comparison. Results: For monophasic single-pulse stimulation, the ratio of anodic vs. cathodic thresholds in all models was within the experimental range despite a much larger ratio in the model by Briaire and Frijns. For biphasic pulse-train stimulation, thresholds as a function of both pulse rate and pulse duration differed between the models, but none matched the experimental observations even coarsely. Similarly, for all other properties including the conduction velocity, action potential shape, and latency, the models presented different outcomes and not all of them fell within the range observed in experiments. Conclusions: While all three models presented similar values in certain single fiber properties to those obtained in experiments, none matched all experimental observations satisfactorily. In particular, the adaptation and temporal integration behaviors were completely missing in all models. Further extensions and analyses are required to explain and simulate realistic auditory nerve fiber responses to electrical stimulation

Sensitivity study of neuronal excitation and cathodal blocking thresholds of myelinated axons for percutaneous auricular vagus nerve stimulation

Objective: Excitation of myelinated nerve fibers is investigated by means of numerical simulations, for the application of percutaneous auricular vagus nerve stimulation (pVNS). High sensitivity to axon diameter is of interest regarding the goal of targeting thicker fibers. Methods: Excitation and blocking thresholds for different pulse types, phase durations, axon depths, axon-electrode distances, temperatures and axon diameters are investigated. The used model consists of a 50 mm long axon and a centrally located needle electrode in a layered medium representing the auricle. Neuronal excitation is simulated using the Frankenhaeuser-Huxley equations for all combinations of parameter values. Results and conclusion: Multiple modes and locations of excitation along the axon were observed, depending on the pulse type and amplitude. When increasing the axon-electrode distance from 1 mm to 2 mm, sensitivity of thresholds to axon depth decreased with ca. 50%, while sensitivity to axon-electrode distance, axon diameter and phase duration each increased with ca. 15% to 20%, except from monophasic anodal pulses, showing a 45% decrease for axon-electrode distance. These trends for axon diameter and axon-electrode distance allow for more selective stimulation of thicker target fibers using monophasic anodal pulses at higher axon-electrode distances. Cathodal monophasic pulses did not perform well due to blocking of the thicker fibers, which was only rarely seen for other pulse types. Significance: Sensitivities of stimulation thresholds to these parameters by numerical simulation reveal how the stimulation parameters can be changed in order to increase therapeutic effect and comfort during pVNS by enabling more selective stimulation

Sensitivity analysis of a numerical model for percutaneous auricular vagus nerve stimulation

Background: Less-invasive percutaneous stimulation of the auricular branch of the vagus nerve (pVNS) gained importance as a possible nonpharmacological treatment for various diseases. The objective is to perform a sensitivity analysis of a realistic numerical model of pVNS and to investigate the effects of the model parameters on the excitation threshold for single and bundled axons. Methods: Sim4Life electrostatic solver and neural tissue models were combined for electromagnetic and neural simulation. The numerical model consisted of a high-resolution model of a human ear, blood vessels, nerves, and three needle electrodes. Investigated parameters include the axon diameter and number, model temperature, ear conductivity, and electrodes' penetration depth and position. Results: The electric field distribution was evaluated. Model temperature and ear conductivity are the non-influential parameters. Axons fiber diameter and the electrodes' penetration depth are the most influential parameters with a maximum threshold voltage sensitivity of 32 mV for each 1 mu m change in the axon diameter and 38 mV for each 0.1 mm change in the electrodes' penetration depth. Conclusions: The established sensitivity analysis allows the identification of the influential and the non-influential parameters with a sensitivity quantification. Results suggest that the electrodes' penetration depth is the most influential parameter