Application of a neuroelectric model to electrocutaneous sensory sensitivity: parameter variation study

Avoiding nerve stimulation in irreversible electroporation: a numerical modeling study

Electroporation based treatments consist in applying one or multiple high voltage pulses to the tissues to be treated. As an undesired side effect, these pulses cause electrical stimulation of excitable tissues such as nerves and muscles. This increases the complexity of the treatments and may pose a risk to the patient. To minimize electrical stimulation during electroporation based treatments, it has been proposed to replace the commonly used monopolar pulses by bursts of short bipolar pulses. In the present study, we have numerically analyzed the rationale for such approach. We have compared different pulsing protocols in terms of their electroporation efficacy and their capability of triggering action potentials in nerves. For that, we have developed a modeling framework that combines numerical models of nerve fibers and experimental data on irreversible electroporation. Our results indicate that, by replacing the conventional relatively long monopolar pulses by bursts of short bipolar pulses, it is possible to ablate a large tissue region without triggering action potentials in a nearby nerve. Our models indicate that this is possible because, as the pulse length of these bipolar pulses is reduced, the stimulation thresholds raise faster than the irreversible electroporation thresholds. We propose that this different dependence on the pulse length is due to the fact that transmembrane charging for nerve fibers is much slower than that of cells treated by electroporation because of their geometrical differences.

A μm-Scale Computational Model of Magnetic Neural Stimulation in Multifascicular Peripheral Nerves

There has been recurring interest in using magnetic neural stimulation for implantable localized stimulation. However, the large stimulation voltages and energies necessary to evoke neuronal activity have tempered this interest. To investigate the potential of magnetic stimulation as a viable methodology and to provide the ability to investigate novel coil designs that can result in lower stimulation threshold voltages and energies, there is a need for a model that accurately predicts the magnetic field-tissue interaction that results in neuronal stimulation. In this study, we provide a computational framework to accurately estimate the stimulation threshold and have validated the model with in vivo magnetic stimulation experiments. To make such predictions, we developed a micrometer-resolution anatomically driven computational model of rat sciatic nerve and quantified the effect of tissue heterogeneity (i.e., fascicular organization, axon distribution, and density) and axonal membrane capacitance on the resulting threshold. Using the multiresolution impedance method, we computed the spatial-temporal distribution of the induced electric field in the nerve and applied this field to a Frankenhaeuser-Huxley axon model in NEURON to simulate the nonlinear mechanisms of the membrane channels. The computational model developed predicts the stimulation thresholds for four magnetic coil designs with different geometrical parameters within the 95% confidence interval (experiments count = 4) of measured in vivo stimulation thresholds for the rat sciatic nerve.

A model for transcutaneous current stimulation: simulations and experiments

Complex nerve models have been developed for describing the generation of action potentials in humans. Such nerve models have primarily been used to model implantable electrical stimulation systems, where the stimulation electrodes are close to the nerve (near-field). To address if these nerve models can also be used to model transcutaneous electrical stimulation (TES) (far-field), we have developed a TES model that comprises a volume conductor and different previously published non-linear nerve models. The volume conductor models the resistive and capacitive properties of electrodes, electrode-skin interface, skin, fat, muscle, and bone. The non-linear nerve models were used to conclude from the potential field within the volume conductor on nerve activation. A comparison of simulated and experimentally measured chronaxie values (a measure for the excitability of nerves) and muscle twitch forces on human volunteers allowed us to conclude that some of the published nerve models can be used in TES models. The presented TES model provides a first step to more extensive model implementations for TES in which e.g., multi-array electrode configurations can be tested.

Spatial relationships in electrostimulation: application to electromagnetic field standards

Spatial relationships in electrostimulation are examined with a myelinated nerve model. Excitation thresholds are determined for: 1) terminated axon within a locally constant electric field; 2) bent axon within a locally constant field; and 3) a field of finite extent over the affected axon. For application to electromagnetic field standards, the minimum excitation threshold of the in situ electric field applies to a straight, terminated axon; a field measurement averaging distance of 5 mm is recommended.

The effect of electrode configuration and duration of deafness on threshold and selectivity of responses to intracochlear electrical stimulation

This report examines the effects of intracochlear electrode configuration and mode of stimulation (bipolar or monopolar) on neural threshold and spatial selectivity in the inferior colliculus (IC) of the cat. Single and multiunit IC recordings were made in three groups of animals; acutely deafened adults (controls), neonatally deafened animals studied at 6 to 18 months of age and neonatally deafened cats studied at 2.5 to 6.5 years. Response thresholds were plotted versus IC depth to measure the spatial distribution of responses. The response selectivity for each stimulating configuration was defined as the width of the resulting spatial tuning curve (STC) measured at 6 dB above threshold. Spiral ganglion cell (SG) survival was examined histologically in all neonatally deafened animals and correlated with physiological results. Animals studied at less than 1.5 years had SG densities of 23.5%-64.4% of normal (mean=42.7%) while animals studied at greater than 2.5 years had densities of 5.1%-18.3% of normal (mean=9.9%). Electrophysiological results include the following. (1) Monopolar thresholds were 7-8 dB lower than bipolar thresholds in the same animals. (2) Varying the configuration of bipolar contacts (measured as radial, offset radial and longitudinal pairs) did not systematically affect IC threshold in either controls or short-term neonatally deafened animals. In contrast, the long-term neonatally deafened animals showed a difference in threshold with each configuration. (3) The spatial distributions (Q(6 dB)) of responses to bipolar stimulation were approximately 40% more restricted than those for monopolar stimulation. (4) The spatial selectivity of neonatally deafened animals studied at ages up to 1.5 years was equal to that of control animals with normal auditory experience. However, selectivity was degraded in the older animals. (5) Selectivity was decreased in some animals with the longitudinal bipolar configuration and multiple response peaks were seen in several cases using this stimulus configuration.

A model of myelinated nerve fibres for electrical prosthesis design

Starting with the spatially extended non-linear node model (Reilly et al., 1985), which incorporates Frankenhaeuser-Huxley non-linearities at each of several nodes in a row, a model is developed to describe many aspects of the behaviour of mammalian nerve fibres in a quantitative way. By taking into account the effects of temperature and by introducing a realistic nerve morphology, a good fit is obtained between the shape, duration and conduction velocity of simulated and in vivo action potentials in mammals. The resulting model correctly predicts the influence of physiological variations of body temperature on various aspects of nerve behaviour. It is shown that the absolute refractory period predicted by the model is within physiological ranges. Both in vivo and in the model, the spike amplitude and the spike conduction velocity are reduced in the relative refractory period. It is concluded that single-node models (although widely used) cannot replace this multiple nonlinear node model, as the stimulus repetition rates that can be followed by the simulated nerve fibre are limited by impulse conduction properties, rather than by the frequency following behaviour of a single node.

A quantitative approach to modeling mammalian myelinated nerve fibers for electrical prosthesis design

This paper presents an upgraded cable model of mammalian myelinated nerve fibers in an extracellularly applied field. The kinetics of the nodes is based upon voltage clamp data in rat motor fibers at 37 degrees C, while the resting membrane potential is computed with the Goldman equation. The resulting spike shape, conduction velocity, strength/duration behavior, and absolute and relative refractory period are in good quantitative agreement with published experimental data in mammals at normal body temperature and at 20 degrees C. Results at intermediate temperatures however, suggest that the widely used concept of a constant Q10 for the rate constants is invalid. In addition, the model generates realistic abortive spikes towards the end of the absolute refractory period and it can describe the consequences of repetitive firing. The results stress the advantages of a multiple nonlinear node model even if only time aspects of nerve behavior are under study. It turned out, that the model presented here describes in vivo neural properties relevant for electrical prosthesis design better than previous models in literature.

Axon termination conditions for electrical stimulation

The cable model for electrical stimulation near the terminal of a passive fiber is derived for excitation by an arbitrary, time-varying, applied extracellular field. Unless the termination impedance is comparable to that of mammalian node of Ranvier, the end-conditions require the longitudinal intracellular current at the fiber terminal to be negligibly small. This requirement substantially alters the membrane potential profile from that obtained with a fiber of infinite length. Stimulation near the end of a fiber may result in lower thresholds and may reverse the anodal/cathodal threshold ratio obtained with stimulation in the mid-portion of the fiber. Chronaxie for stimulation near the terminal may be much smaller than at a distance from the terminal and the strength-duration curve may be nonmonotonic. These differences may have significant implications for any application of electrical stimulation where fiber terminations may play a role in the excitatory process.

Magnetic field excitation of peripheral nerves and the heart: a comparison of thresholds

Time-varying magnetic fields can theoretically excite the heart or peripheral nerves. Relative excitation thresholds of nerve and heart are compared using ellipsoidal representations of the human torso. Relative magnetic thresholds depend on the excitability of nerve and cardiac tissue, the geometric positions of anatomical features, stimulus waveform features and the direction and spatial distribution of the incident magnetic field. Minimum electric field thresholds for excitation of nerve and heart do not differ greatly. Nevertheless, nerve and heart magnetic thresholds may be disparate because of factors related to body geometry and the dependence of excitability on the stimulus waveform.

Electrotactile and vibrotactile displays for sensory substitution systems

Sensory substitution systems provide their users with environmental information through a human sensory channel (eye, ear, or skin) different from that normally used, or with the information processed in some useful way. We review the methods used to present visual, auditory, and modified tactile information to the skin. First, we discuss present and potential future applications of sensory substitution, including tactile vision substitution (TVS), tactile auditory substitution, and remote tactile sensing or feedback (teletouch). Next, we review the relevant sensory physiology of the skin, including both the mechanisms of normal touch and the mechanisms and sensations associated with electrical stimulation of the skin using surface electrodes (electrotactile (also called electrocutaneous) stimulation). We briefly summarize the information-processing ability of the tactile sense and its relevance to sensory substitution. Finally, we discuss the limitations of current tactile display technologies and suggest areas requiring further research for sensory substitution systems to become more practical.