Restoring the encoding properties of a stochastic neuron model by an exogenous noise
- 1Department of Information Engineering, Electronics and Telecommunications, Sapienza University of Rome, Rome, Italy
- 2Italian Inter-University Center for the Study of Electromagnetic Fields and Biological Systems, Genova, Italy
Here we evaluate the possibility of improving the encoding properties of an impaired neuronal system by superimposing an exogenous noise to an external electric stimulation signal. The approach is based on the use of mathematical neuron models consisting of stochastic HH-like circuit, where the impairment of the endogenous presynaptic inputs is described as a subthreshold injected current and the exogenous stimulation signal is a sinusoidal voltage perturbation across the membrane. Our results indicate that a correlated Gaussian noise, added to the sinusoidal signal can significantly increase the encoding properties of the impaired system, through the Stochastic Resonance (SR) phenomenon. These results suggest that an exogenous noise, suitably tailored, could improve the efficacy of those stimulation techniques used in neuronal systems, where the presynaptic sensory neurons are impaired and have to be artificially bypassed.
Introduction
Different techniques for the stimulation of neuronal systems have been developed. Some are based on magnetic coupling, such as the Transcranial Magnetic Stimulation (TMS) (Corthout et al., 2001) and low-intensity magnetic stimulation (Di Lazzaro et al., 2013), while others use electric fields, such as Deep Brain Stimulation (DBS) (Okun et al., 2012; Paffi et al., 2013a), Functional Electric Stimulation (FES) of peripheral nerves (Peckham and Knutson, 2005), cochlear prostheses (Wilson et al., 1991; Clark, 2003), and Intracortical Microstimulation (ICMS) (Brock et al., 2013; Overstreet et al., 2013).
Despite great interest in such applications and the experimental activities to evaluate the effect of electromagnetic fields on single neurons and networks (Marchionni et al., 2006; Platano et al., 2007; Ahmed and Wieraszko, 2009; Moretti et al., 2013), the mechanisms of action are not clearly understood (Apollonio et al., 2013; Di Lazzaro et al., 2013) and the techniques are not yet optimized.
Theoretical studies to understand neuronal system functioning are based on biophysical models. At the single neuron level, a lot of work has been done using simple (Mainen et al., 1995; Rapp et al., 1996; Rinzel and Ermentrout, 1998) or augmented (Tateno et al., 1998; Pospischil et al., 2008) Hodgkin and Huxley (HH) descriptions (Hodgkin and Huxley, 1952), both under physiologic conditions and the action of exogenous stimulations (Mino et al., 2004; Giannì et al., 2005, 2006; Camera et al., 2012, 2013).
The HH model is a nonlinear active circuit, which behaves as an oscillator if the injected constant current (stimulation current) overcomes a threshold (Rinzel and Ermentrout, 1998). Such a current represents a physiologic stimulation from all the presynaptic neurons. Depending on this parameter, the HH model can display a stable resting state or/and a stable limit cycle (Hassard, 1978), corresponding to a periodic oscillation of the membrane voltage in the form of a spike train.
However, the HH model does not adequately take into account the stochastic behavior of neurons. Electrophysiology recordings have shown that actual neurons have an intrinsic stochastic behavior (Sigworth, 1980; Mainen and Sejnowski, 1995; Dorval and White, 2005), shown by the unreliable responses and non-deterministic current thresholds for firing. This is due to the noisy environment of the neuron, in particular to the intrinsic stochasticity of channel gating (channel noise) (White et al., 2000). The level of this noise decreases with the number of ionic channels and so depends on the channel density and the size of the neuron considered (Schneidman et al., 1998; White et al., 2000).
Several authors have provided stochastic neuron models (Schneidman et al., 1998; White et al., 2000) to study how channel noise may improve the encoding of a physiologic or artificial input stimulation (Schneidman et al., 1998; Jung and Shuai, 2001; Manwani et al., 2002; Mino et al., 2004; Giannì et al., 2005, 2006, 2007; Woo et al., 2010) with different frequency content (Bulsara et al., 1993; Longtin, 1993; Liu et al., 1999; Rudolph and Destexhe, 2001; Yu et al., 2001a; Giannì et al., 2005, 2006, 2007; Liberti et al., 2009a; Sengupta et al., 2013). This positive role of the endogenous noise in signal encoding and processing observed in neuronal models was often attributed to the Stochastic Resonance (SR) phenomenon (Gammaitoni, 1998; Moss et al., 2004; McDonnell and Abbott, 2009).
However the endogenous noise is essentially related to the type of neuron and is difficult to manipulate and control. From a biomedical perspective it is more interesting to use theoretical models to elucidate the role of an artificial noise externally applied.
In previous studies, we showed that a suitably tailored exogenous noise could increase firing activity and improve signal detection through the SR mechanism in compartmental models of neuronal systems with reduced levels of endogenous noise (Paffi et al., 2006, 2007, 2013b). Here we extend our idea to systems where the presynaptic stimulation was lowered due to impairment of sensory neurons that have to be artificially bypassed at different levels of the neuronal pathway toward the cortical region of sensorial processing. Examples are the cochlear prostheses, where an electrode inserted in the cochlea directly stimulates the fibers of the auditory nerve (Wilson et al., 1991; Clark, 2003), and ICMS to deliver sensory perceptions to the auditory or visual cortex (Brock et al., 2013; Overstreet et al., 2013).
To test our idea for the optimization of stimulation techniques, a simple and well-characterized HH neuron model is considered. The normal functioning is modeled with a suprathreshold input current, and pathologic conditions with a subthreshold presynaptic stimulation. Different kinds of neurons of the sensory pathway, characterized by different sizes and, hence, by different levels of endogenous noise, are accounted for by changing the number of ionic channels.
The first step is to demonstrate that the detectability of the exogenous stimulation signal is degraded in impaired sensory neurons as a function of the endogenous noise level. Then we want to show that a suitably tailored exogenous noise can partially restore the signal encoding in the impaired neurons, in agreement with the SR phenomenon.
The main aim of this paper is to show that the reduced encoding capability of pathologically understimulated neuronal systems can be improved using an exogenous noise, opening the way for prosthetic applicators delivering the exogenous stimulation signal and noise.
The paper is organized as follows. In Section Materials and Methods the stochastic neuron model is described (Section Stochastic Neuron Model), together with the methods for introducing the exogenous electric signal (Section Introduction of the Exogenous Signal) and noise (Section Introduction of the Exogenous Noise) and for evaluating neuron excitability and signal encoding (Section Observables). Results are presented and discussed in Section Results, without (Section Encoding Features of the Model) and with the exogenous noise (Section Role of the Exogenous Noise). Finally, in Section Discussion and Conclusions, results are discussed and conclusions are given.
Materials and Methods
Stochastic Neuron Model
We use a stochastic neuron model (Fitzhugh, 1965; Sigworth, 1980; Clay and DeFelice, 1983; Schneidman et al., 1998), derived from the HH model (Hodgkin and Huxley, 1952).
The equivalent scheme is shown by the parallel combination of five current branches, as in Figure 1.
Figure 1. Circuital representation of the neuron model. Cm is the specific membrane capacitance; gL is the leakage specific conductance; gNa and gk are Sodium and Potassium specific conductances, voltage (Vm) dependent; I0 is the specific stimulation current; VNa, VK, VL are the reversal potentials for Sodium, Potassium, and leakage currents, respectively; VES and Vnoise are the voltage perturbations, due to the exogenous electric signal and the exogenous Gaussian noise, superimposed to the physiological membrane voltage V; Vm is the total voltage between the intracellular and the extracellular space.
Sodium and Potassium currents were calculated using a channel-state-tracking algorithm (Mino et al., 2002) where the ionic channels are modeled as the combination of independent gating particles whose dynamics is well-described by Markov chains (Rubinstein, 1995; Mino et al., 2002).
The number of voltage-gated channels belonging to a given population determines the level of endogenous noise in the system (Schneidman et al., 1998; White et al., 2000). This decreases with the square root of the channel number and, for a given channel density, with the size of the neuron considered (Schneidman et al., 1998; White et al., 2000).
Therefore, the endogenous level of noise is typical of each kind of neuron and strongly varies with neuron size.
For example, the auditory fibers have a small diameter estimated between 1 and 9 μm in (Engstrom and Wersall, 1958), or between 1 and 5 μm in Gleich and Wilson (1993). The cell body is larger, with a cross-sectional area of 300 μm2 (Liberman and Oliver, 1984; Woo et al., 2010). At the brainstem in the cochlear nucleus, the bushy cells that receive inputs from the auditory fibers range between 300 and 1200 μm2 in area (Sento and Ryugo, 1989). Cortical neurons of Layer I have the soma size ranging from 7 to 17 μm in diameter (Hestrin and Armstrong, 1996), corresponding to an area from 154 to 900 μm2, if one assumes a spherical morphology.
In this work, membrane patches of 200, 300, and 600 μm2 have been considered, corresponding to NNa = 12,000, NNa = 18,000, and NNa = 36,000, respectively, for constant channel densities of 18 NK/μm2 and 60 NNa/μm2 (Hodgkin and Huxley, 1952).
A neuron of 200 μm2 could represent the largest axon in the auditory fiber (Engstrom and Wersall, 1958) or a small-sized neuron of Layer I in the cortex (Hestrin and Armstrong, 1996); one at 300 μm2 could be the smallest bushy cells in the cochlear nucleus (Sento and Ryugo, 1989) or a medium-sized cortical neuron of Layer I, and neurons of 600 μm2 represent a typical area of bushy cells and cortical neurons.
The current generator I0 in Figure 1 was set to different values: I0 = 2 μA/cm2, I0 = 4 μA/cm2, I0 = 7 μA/cm2, below or above the firing threshold (I0th = 6.3 μA/cm2) of the associated deterministic oscillator (Rinzel and Ermentrout, 1998).
Considering I0 as the total stimulation from the presynaptic neurons (presynaptic current), particularly from the sensory receptors, these conditions can be representative of a neuron normally stimulated by the sensory inputs (I0 = 7 μA/cm2) or a neuron where the receptors stimuli are slightly (I0 = 4 μA/cm2) or significantly (I0 = 2 μA/cm2) reduced. This is a typical impairment induced by aging, direct damage, or degenerative diseases of the sensory receptors such as the cells of the organ of Corti in the Cochlea (Ritter et al., 1981).
The neuron stochastic model used HH parameters (Hodgkin and Huxley, 1952) in the C++ environment using the Forward Euler integration method with time steps of 10 μs.
The output of the model is the time course of the voltage across the membrane V(t).
Introduction of the Exogenous Signal
There is consensus in literature that the effect of an exogenous magnetic or electric stimulation delivered by a coil or an implanted electrode is the creation of an electric field in the tissue that in turn induces a perturbation on the neuron membrane voltage (Foster and Schwan, 1986; Mino et al., 2004; Giannì et al., 2006; Merla et al., 2012). Therefore, unless a current is directly injected across the membrane, the interaction between the exogenous signal and the neuron membrane must be inserted as a voltage generator in series with the neuron circuital model (Figure 1), as already done in number of studies (Tsong and Astumian, 1987; Mino et al., 2004; Giannì et al., 2006; Woo et al., 2010; Paffi et al., 2013b). As previously discussed (Giannì et al., 2006), this additive voltage can describe the non-linear interaction between the exogenous signal and the neuron activity (Stodilka et al., 2011), since it induces a perturbation in the dynamics of the voltage-dependent Sodium and Potassium channels. The effects on the ionic currents are temporally integrated and reflected back on the membrane potential, showing a feedback interaction mechanism (Apollonio et al., 2000).
A weak deterministic sinusoidal signal, of amplitude VES = 500 μV was considered for frequencies between 10 and 500 Hz. The term “weak” means that it does not induce firing activity in the subthreshold neuron, provided no sources of stochasticity, endogenous (channel noise), or exogenous, are present.
Introduction of the Exogenous Noise
The endogenous noise, due to the neuron size, i.e., the number of ionic channels, is an intrinsic feature of each neuron type and cannot be artificially tuned or modulated. Conversely, a well-defined exogenous noise can be added to the sinusoidal signal and suitably tuned in terms of power and frequency content.
In the circuit scheme of Figure 1 the exogenous noise was modeled as a random voltage source whose level is given by, where ζ(t) is a Gaussian process with zero mean and unitary variance. Accordingly, the noise power D is measured in mV2 and was varied in the range [0.7-25] (mV2).
The well-known equation describing the current balance of the HH circuit in the presence of an exogenous signal and noise shown in Figure 1, becomes:
As shown by Paffi et al. (2013b) the voltage noise was filtered to obtain a Lorentzian Power Spectral Density (PSD):
The use of a Lorentzian behavior is a straightforward choice since, under the passive linear approximation (Steinmetz et al., 2000), the neuronal membrane behaves like a single pole filter with a time constant equal to the membrane capacitance divided by the total conductivity of the ionic channels in the patch (Rinzel and Ermentrout, 1998).
The cutoff angular frequency ωc = 2.5 × 103 rad/s of the Lorentzian filter was chosen on the basis of theoretical calculations (Paffi et al., 2013b) and numerical simulations (Paffi et al., 2007).
Figure 2 shows the features of the exogenous voltage noise in terms of the Gaussian distribution (Figure 2A) and correlation properties by the PSD (Figure 2B). The PSD was estimated as the Periodogram averaged over 10 runs of the Gaussian process 1 s long.
Figure 2. (A) Histogram showing the zero-mean Gaussian distribution of the exogenous voltage noise with power D = 2 mV2; the standard deviation (σ) of the distribution is equal to the square root of D. (B) Normalized Power Spectral Density of the exogenous voltage noise (blue line) compared with the theoretical Lorentzian spectrum with cutoff angular frequency ωc = 2.5 ×103 rad/s (red line).
Such noise can increase the firing activity of the subthreshold neurons as already demonstrated by Paffi et al. (2013b).
Observables
Due to the stochasticity of the model, the neuron properties, such as the number of spike per second or the frequency content of the spike sequence, can be calculated only as statistical values over a population of R runs of the model.
To quantify the excitability of the neuron, the average number of spikes per second was calculated over 100 V(t) traces 1 s long, together with the standard error.
To determine the “time encoding” of the neuron model, that is, the capability of encoding different input signals in the spike timing within the firing sequence, the Periodogram has been used as a spectral estimator. As a preliminary step to retain information only on the sequence of spikes (Giannì et al., 2005, 2006), disregarding their shape, the time course of the membrane voltage V(t) over 1 s, has been converted into a time series of Dirac pulses U(t), each corresponding to a spike, having height 100 mV as suggested in Gluckman et al. (1996); Levin and Miller (1996) and Yu et al. (2001a,b).
The PSD, averaged over R runs (R = 100 in the presence of the channel noise alone; R = 300 with the presence of channel and exogenous noise) of the signal U(t), has been calculated using the Fast Fourier Transform (FFT) algorithm. It is worth noting that in the presence of both noises the spectral estimator with R = 100 was not satisfactory due to high variance; so to calculate results of Section Role of the Exogenous Noise it was necessary to increase the number of runs from 100 to 300. For each frequency point, the standard error of the PSD was associated with the average value.
The signal to noise ratio (SNR), calculated as the ratio between the strength of the peak of the average spectrum at the forcing frequency (fs) and the background average spectrum around the same frequency (Figure 3) (Gluckman et al., 1996; Levin and Miller, 1996; Gammaitoni, 1998; Yu et al., 2001a,b; Giannì et al., 2006; Paffi et al., 2013b), was used to evaluate the signal detectability as a function of the signal frequency or the exogenous noise level. The background spectrum at the signal frequency was estimated as the average between the values assumed by the PSD 1 Hz before and 1 Hz after the signal frequency (Gammaitoni, 1998; Giannì et al., 2006; Paffi et al., 2013b). The standard error of the SNR was calculated using the standard errors propagation of the correlated variables PSD(fs), PSD(fs–1), PSD(fs + 1). The correlation coefficient of these adjacent samples was calculated to be around 0.9.
Figure 3. Average PSD with the associated standard error (R = 100) of the output sequence U(t) for the neuron model with patch area 300 μm2 and I0 = 4 μA/cm2 (A) or I0 = 7 μA/cm2 (B). The applied exogenous signal is a sinusoid at 150 Hz, 500 μV of amplitude.
Matlab functions have been used to extract the aforementioned observables from the model output.
Results
Encoding Features of the Model
It is known from literature that HH models exhibit a frequency sensitivity that depends on the model parameters, particularly on the constant input current I0 (Liu et al., 1999; Yu et al., 2001a; Giannì et al., 2006).
Here we examine and compare the encoding capability of the neuron model when I0 assumes values of 2, 4, and 7 μA/cm2 and the patch area, determining the endogenous noise, is equal to 200, 300, or 600 μm2.
As described in Section Stochastic Neuron Model, the exogenous signal was a deterministic sinusoid of amplitude 500 μV and frequency spanning from 10 to 500 Hz, applied as a voltage perturbation over the neuron membrane.
As already shown (Paffi et al., 2013b), the weak exogenous signal does not affect the neuron excitability, i.e., the number of spikes per second, independently of the frequency considered. However, this value significantly changes with the constant input current I0 and the endogenous noise level, as shown in Figure 4, where the number of spikes per second, averaged over the results obtained for all frequency values, are displayed together with the standard errors.
Figure 4. Mean number of spikes per second and standard error calculated over the frequency of the input sinusoidal signal vs. the input constant current I0 for patch areas: 200 (blue line), 300 (purple line), and 600 (orange line) μm2.
From Figure 4 it is evident that the standard errors on the number of spikes are very small with respect to variations due to I0 or to the patch areas, confirming that, when the applied signal is as low as 500 μV, the frequency encoding mechanism is not likely to occur.
As expected, the “impaired” neurons (I0 = 2 μA/cm2, I0 = 4 μA/cm2) are much less stimulated than the “healthy” one (I0 = 7 μA/cm2). This reduction in firing activity is likely to negatively affect the neuron encoding capability. In particular, the “severely impaired” neuron for the biggest patch area (600 μm2) is almost silent, losing any possibility of signal detection. For the “healthy” neuron the endogenous noise does not change the excitability that is completely determined by the input current; conversely, in subthreshold conditions, i.e., for the “impaired” neurons, the higher the endogenous noise level the higher the firing rate.
This shows that different neurons may be more or less sensitive to impairment depending on their channel noise. That noise has a beneficial effect on the firing activity of “impaired” neurons, suggesting a positive role of an exogenous noise.
If one considers the total power of the sequence U(t) 1 s long (Table 1), the same behavior is observed. The output power does not depend on the patch area for I0 = 7 μA/cm2 but increases up to 17 times for I0 = 2 μA/cm2 if the patch area decreases from 600 to 200 μm2.
Table 1. Total power (mV2) of the output sequence U(t) 1 s long calculated for the different values of patch areas and input currents.
Time encoding performances of the neuron model were measured using the SNR, as described in Section Observables. Indeed, if time encoding occurs, the PSD of U(t) will present a component at the signal frequency higher than the background level around the signal frequency, as in Figure 3 for a 150 Hz signal, a patch of 300 μm2, input currents I0 = 4 μA/cm2 (Figure 3A) and I0 = 7 μA/cm2 (Figure 3B).
Figure 3 shows that the exogenous signal can synchronize some spike events with its own frequency, leading to a frequency peak at 150 Hz emerging from the background PSD with a consequent SNR value greater than one.
The frequency sensitivity of the “healthy” neuron, from 10 to 500 Hz, is plotted in Figure 5 for the three patch areas. At first one can observe a bell-shaped behavior of the SNR, with a peak centered around 150 Hz. As the patch area decreases, i.e., the endogenous noise increases, the curves become smoothed and the maximum SNR decreases. This indicates that endogenous noise has a detrimental effect on the encoding capability in the healthy neuron.
Figure 5. Mean SNR and standard error (R = 100) vs. the signal frequency for I0 = 7 μA/cm2 and membrane patches of 200, 300, and 600 μm2; the signal is a sinusoid with amplitude VES = 500 μV and frequency ranging from 10 to 500 Hz.
Looking at the “impaired” neurons (Figure 6), the maximum of the curve shifts toward lower frequencies and the SNR at 150 Hz, where the “healthy” neuron exhibits the maximum sensitivity, significantly decreases (always below 3). For the “severely impaired” neuron and the largest size, the signal becomes undetectable (SNR almost equal to one with higher error bars), since the firing activity is almost completely suppressed. From a biophysical point of view, this can be interpreted as a severe worsening in neuron performances if the afferent inputs are lacking. Only the smallest neuron, for I0 = 4 μA/cm2 shows an encoding capability at 150 Hz similar to that of the “healthy” neuron with the same size, suggesting again a positive role of the noise in understimulated neurons.
Figure 6. Mean SNR and standard error (R = 100) vs. the signal frequency for I0 = 2 and 4 μA/cm2 and membrane patches of 200 μm2 (A), 300 μm2 (B), and 600 μm2 (C); the signal is a sinusoid with amplitude VES = 500 μV and frequency ranging from 10 to 500 Hz.
With an exogenous sinusoidal signal of 150 Hz and amplitude 500 μV, in Section Role of the Exogenous Noise we investigated the possibility of restoring degraded performances by adding a correlated noise to the system from the outside.
Role of the Exogenous Noise
As shown in Section Encoding Features of the Model, the effect of the impairment of the afferent stimulation is a drastic reduction of the neuron firing and signal encoding, especially for neurons with a lower endogenous noise, i.e., the larger ones.
From a biomedical perspective the question arises as to whether the exogenous noise can improve the firing activity and the encoding capability of the “impaired” neurons in terms of number of spikes per second and SNR. Therefore, we have added the correlated voltage noise described in Section Introduction of the Exogenous Noise to the neuron models in subthreshold conditions (I0 = 2 and 4 μA/cm2).
Among the six different conditions shown in Figure 6, we have not considered the patch of 600 μm2 with I0 = 2 μA/cm2, since the neuron has no residual activity, and the patch of 200 μm2 with I0 = 4 μA/cm2, since the encoding performances at 150 Hz are already comparable to those of “healthy” neurons with the same size.
The SNRs as a function of D are shown in Figure 7 for the “severely impaired” neuron (I0 = 2 μA/cm2) with area 200 (Figure 7A) and 300 μm2 (Figure 7B) and for the “impaired” neuron (I0 = 4 μA/cm2) with area 300 μm2 (Figure 7C) and 600 μm2 (Figure 7D).
Figure 7. Mean SNR and standard error (R = 300) as a function of the variance of the exogenous noise (D) for I0 = 2 μA/cm2 and membrane patches of 200 (A) and 300 μm2 (B), and for I0 = 4 μA/cm2 and membrane patches of 300 μm2 (C) and 600 μm2 (D); the signal is a sinusoid with amplitude VES = 500 μV and frequency f = 150 Hz; the exogenous noise is a zero mean Gaussian process with a Lorentzian spectrum (ωc = 2.5×103 rad/s).
Figure 7 shows that, for all the conditions considered, there exist noise levels that improve the SNR, showing the typical behavior of SR (Gammaitoni, 1998; Moss et al., 2004; McDonnell and Abbott, 2009). For each condition it is possible to identify an optimum noise level (Dopt), where the SNR is maximum, that depends on the membrane patch and the input current I0. In particular, for the same I0, Dopt is lower for the smaller patch area, showing that the higher the channel noise, the lower the exogenous noise to be supplied, in agreement with previous results (Schmid et al., 2001; Paffi et al., 2013b).
Table 2 summarizes, for each condition studied, the SNR values exhibited by the model at 150 Hz without the endogenous noise and with the optimum noise power (Dopt) reported in the last column. The table shows that the SNR increases in the presence of the optimum exogenous noise. This improvement is not significant for the neuron with I0 = 4 μA/cm2 and membrane patch 300 μm2, where the encoding performances without the exogenous noise (SNR = 2.51 ± 0.10) are still acceptable due to high levels of endogenous activity. In all the other cases, one may observe a significant increase in SNR by up to 57%.
Table 2. SNR values exhibited by the neuron model at 150 Hz without and with the exogenous noise at the optimum level Dopt.
Although the performances of the “impaired” neuron are not been completely restored, they are significantly improved, confirming the potentially beneficial effect of an exogenous noise according to the SR paradigm.
Another effect of the exogenous noise is a considerable increase of the number of spikes per second and, consequently, of the power associated with the spike sequence (U(t)). Unlike the SNR, these quantities exhibit a monotonic increase with D, approaching asymptotic values, as shown in Figure 8 for the number of spikes per second.
Figure 8. Mean number of spikes per second and standard error (R = 300) as a function of the variance of the exogenous noise (D) for I0 = 2 μA/cm2 and membrane patches of 200 μm2 (solid blue line) and 300 μm2 (solid purple line), and for I0 = 4 μA/cm2 and membrane patches of 300 μm2 (dashed purple line) and 600 μm2 (dashed orange line); the signal is a sinusoid with amplitude VES = 500 μV and frequency f = 150 Hz; the exogenous noise is a zero mean Gaussian process with a Lorentzian spectrum (ωc = 2.5×103 rad/s).
This means that the increase in firing activity does not necessary imply an improved encoding capability. Indeed, although a minimum number of spikes per second is required to efficiently encode a 150 Hz sinusoidal signal, if noise exceeds the optimum level, neuron activity is dominated by noise and less correlated with the signal, in agreement with the SR phenomenon (Gammaitoni, 1998; Moss et al., 2004; McDonnell and Abbott, 2009).
Table 3 shows the number of spikes per second and the power assumed by U(t) corresponding to Dopt, compared to their values in the absence of the exogenous noise. In all cases the number of spikes at the optimum noise level is between 28 and 38, with a power from 1.4 to 1.8 mV2, suggesting that these numbers of spikes per second can efficiently encode a sinusoid at 150 Hz in a subthreshold neuron. Thus, it is not surprising that the model with a 600 μm2 area and I0 = 4 μA/cm2 shows good encoding properties at 150 Hz (SNR = 3) (see Figure 6C) with 30 spikes per second (see Figure 4).
Table 3. Number of spikes per second and power of U(t) for the neuron model without and with the exogenous noise at the optimum level Dopt.
Discussion and Conclusions
In this work, moving from the original HH description, a stochastic neuron model has been developed and the presence of an exogenous signal and noise, representative of a possible electric or magnetic stimulation, has been added. The firing activity of the model has been studied for three different levels of channel noise, corresponding to different patch areas (200, 300, and 600 μm2) and for three input currents, representative of a “healthy” neuron (I0 = 7 μA/cm2), an “impaired” (I0 = 4 μA/cm2), neuron and a “severely impaired” neuron (I0 = 2 μA/cm2), where the afferent stimulation is reduced due to aging or degenerative diseases.
Results indicate that the “impaired” neurons are much less excited (less than 10 spikes/s for a patch of 600 μm2, I0 equal to 2 and 4 μA/cm2, and a patch of 300 μm2, I0 = 2 μA/cm2), suggesting reduced performances in signal encoding and processing.
The presence of the exogenous sinusoidal signal (VES = 500 μV; f = [10–500] Hz) does not significantly change the firing frequency (Figure 4), confirming that the neuron does not use the frequency encoding paradigm to sense such low-level alternate signals. On the contrary, the PSD of the spiking sequence U(t) reveals the presence of a peak corresponding to the signal frequency, suggesting a time encoding mechanism (Figure 3).
Results of the SNR, chosen as a measure of signal encoding, show a strong sensitivity to the signal frequency, as suggested by the bell-shaped curves of Figures 5, 6. For example, the “healthy” neuron does not sense signals at 10 or 500 Hz (SNR = 1), whereas, for a frequency of 150 Hz, the SNR can be as high as 6.8 (Figure 5). Such frequency sensitivity significantly depends on the stimulation current I0. This result confirms that the encoding capability of the neuron can be strongly altered by a decrease in the presynaptic input current.
The neuron size, and thus the endogenous noise, mainly affects the maximum SNR, indicating that the encoding capability is a function of the type of neuron (e.g., auditory fiber, bushy cell, cortical neuron). Interestingly, the endogenous noise reduces the encoding capability of the “healthy” neuron but facilitates the signal detection in “impaired neurons,” suggesting a similar behavior if the noise is delivered externally together with the stimulating sinusoidal signal.
An exogenous voltage noise, modeled as a zero-mean Gaussian process with a Lorentzian spectrum and variable power (D), has been added to the models of the “impaired” neurons. The SNR obtained as a function of the noise power exhibits a typical bell shaped behavior with a maximum value corresponding to a well-defined value of D (Dopt), which depends on the values considered for the I0 and patch area.
The exogenous noise at the optimum levels can significantly increase the SNR of the “impaired” neurons at 150 Hz (up to 58%, depending on the neuron size and the impairment level). Since the neuron model takes advantage of the noise to improve the detection of a weak sinusoidal input signal, the observed behavior can be attributed to the well-known SR phenomenon (Gammaitoni, 1998; Moss et al., 2004; McDonnell and Abbott, 2009) (Figure 7).
Beside the improvement in the SNR the exogenous noise induces an increased firing activity that for Dopt is characterized by 28–38 spikes/s.
These results are significant if considered as a proof of concept on how to use artificial exogenous noise to restore the functionalities of signal detection and processing in impaired neuronal systems.
This is a first step toward the optimization of specific biomedical applications such as cochlear prosthesis (Morse and Roper, 2000; Rattay, 2000; Stocks et al., 2002; Rubinstein and Hong, 2003) and ICMS (Overstreet et al., 2013).
Further developments of this work could be the optimization of the exogenous noise, in terms of the spectrum shape and/or the kind of stochastic process and a more accurate description of different neurons in terms of type and number of channels.
Finally, depending on the particular biomedical application, the typical waveform and amplitude of the exogenous voltage signal superimposed on the physiological transmembrane potential can be calculated using dosimetric (Liberti et al., 2007; Maggio et al., 2009, 2010; Paffi et al., 2013a,c, 2015) and microdosimetric techniques (Liberti et al., 2009b; Merla et al., 2010, 2011, 2012; Denzi et al., 2013, 2015), as described in the integrated methodology proposed by Apollonio et al. (2000, 2013).
Author Contributions
ML and FA developed the underlying concept of this study, with contributions from AP. AP and FC prepared computer-simulation codes and methods, and AP, FC, and ML carried out the analysis. All authors discussed results, interpreted data, and formulated findings. AP wrote the manuscript, with some contributions from FA and ML.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors want to thank Michael Repacholi for assistance in editing the final text.
References
Ahmed, Z., and Wieraszko, A. (2009). The influence of pulsed magnetic fields (PMFs) on nonsynaptic potentials recorded from the central and peripheral nervous systems in vitro. Bioelectromagnetics 30, 621–630. doi: 10.1002/bem.20516
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Apollonio, F., Liberti, M., d'Inzeo, G., and Tarricone, L. (2000). Integrated models for analysis of biological effects of EM fields used for mobile communications. IEEE Trans. Microw. Theory Tech. 48, 2082–2093. doi: 10.1109/22.884199
Apollonio, F., Liberti, M., Paffi, A., Merla, C., Marracino, P., Denzi, A., et al. (2013). Feasibility for microwaves energy to affect biological systems via nonthermal mechanisms: a systematic approach. Trans. Microw. Theory Tech. 61, 2031–2045. doi: 10.1109/TMTT.2013.2250298
Brock, A. A., Friedman, R. M., Fan, R. H., and Roe, A. W. (2013). Optical imaging of cortical networks via intracortical microstimulation. J. Neurophysiol. 110, 2670–2678. doi: 10.1152/jn.00879.2012
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Bulsara, A. R., Maren, A. J., and Schmera, G. (1993). Single effective neuron: Dendritic coupling effects and stochastic resonance. Biol. Cybern. 70, 145–156. doi: 10.1007/BF00200828
Camera, F., Paffi, A., Merla, C., Denzi, A., Apollonio, F., Marracino, P., et al. (2012). “Effects of nanosecond pulsed electric fields on the activity of a Hodgkin and Huxley neuron model,” in Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society (San Diego, CA: EMBS), 2567–2570.
Camera, F., Thomas, A. W., Paffi, A., d'Inzeo, G., Apollonio, F., Prato, F. S., et al. (2013). “Effects of pulsed magnetic field on neurons: Cnp signal silences a feed-forward network model,” in Proceedings of the 6th International IEEE EMBS Conference on Neural Engineering (San Diego, CA), 223–226.
Clark, G. M. (2003). Cochlear Implants: Fundamentals and Applications. New York, NY: Springer-Verlag.
Clay, J., and DeFelice, L. (1983). Relationship between membrane excitability and single channel open-close kinetics. Biophys. J. 42, 151–157. doi: 10.1016/S0006-3495(83)84381-1
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Corthout, E., Barker, A., and Cowey, A. (2001). Transcranial magnetic stimulation. Exp. Brain Res. 141, 128–132. doi: 10.1007/s002210100860
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Denzi, A., Merla, C., Camilleri, P., Paffi, A., d'Inzeo, G., Apollonio, F., et al. (2013). Microdosimetric study for nanosecond pulsed electric fields on a cell circuit model with nucleus. J. Membr. Biol. 246, 761–767. doi: 10.1007/s00232-013-9546-7
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Denzi, A., Merla, C., Palego, C., Paffi, A., Ning, Y., Multari, C., et al. (2015). Assessment of cytoplasm conductivity by nanosecond pulsed electric fields. IEEE Trans. Biomed. Eng. doi: 10.1109/TBME.2015.2399250. [Epub ahead of print].
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Di Lazzaro, V., Capone, F., Apollonio, F., Borea, P. A., Cadossi, R., Fassina, L., et al. (2013). A consensus panel review of central nervous system effects of the exposure to low-intensity extremely low-frequency magnetic fields. Brain Stimul. 6, 469–476. doi: 10.1016/j.brs.2013.01.004
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Dorval, A. D., and White, J. A. (2005). Channel noise is essential for perithreshold oscillations in entorhinal stellate neurons. J. Neurosci. 25, 10025–10028. doi: 10.1523/JNEUROSCI.3557-05.2005
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Engstrom, H., and Wersall, L. (1958). Structure and innervation of the inner ear sensory epithelia. Int. Rev. Cytol. 7, 535. doi: 10.1016/S0074-7696(08)62695-9
Fitzhugh, R. (1965). A kinetic model of the conductance changes in nerve membrane. J. Cell. Comp. Physiol. 66, 111–118. doi: 10.1002/jcp.1030660518
Foster, K. R., and Schwan, H. P. (1986). CRC Handbook of Biological Effects of Electromagnetic Fields. Boca Raton: CRC.
Gammaitoni, L. (1998). Stochastic resonance. Rev. Mod. Phys. 70, 223–287. doi: 10.1103/RevModPhys.70.223
Giannì, M., Maggio, F., Paffi, A., Liberti, M., Apollonio, F., and d'Inzeo, G. (2007). “Enhancement of EM signal detectability in a realistic model of feedforward neuronal network,” in 3rd International Conference on Neural Engineering, IEEE EMBS (Hawaii), 684–687.
Giannì, M., Paffi, A., Liberti, M., Apollonio, F., and d'Inzeo, G. (2005). Channel Noise may tune Electromagnetic fields detectability in neurons: stochastic resonance paradigm in a HH-like model. Wseas Trans. Commun. 4, 1406–1410.
Giannì, M., Paffi, A., Liberti, M., Apollonio, F., and d'Inzeo, G. (2006). Modeling electromagnetic fields detectability in a HH-like neuronal system: stochastic resonance and window behavior. Biol. Cybern. 94, 118–127. doi: 10.1007/s00422-005-0029-5
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Gleich, O., and Wilson, S. (1993). The diameters of guinea pig auditory nerve fibres: distribution and correlation with spontaneous rate. Hear. Res. 71, 69–79. doi: 10.1016/0378-5955(93)90022-S
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Gluckman, B. J., Netoff, T. I., Neel, E. J., Ditto, W. L., Spano, M. L., and Schiff, S. J. (1996). Stochastic resonance in a neuronal network from mammalian brain. Phys. Rev. Lett. 77, 4098–4101. doi: 10.1103/PhysRevLett.77.4098
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Hassard, B. D. (1978). Bifurcation of periodic solutions of the Hodgkin-Huxley model for the squid giant asson. J. Theor. Biol. 71, 401–420. doi: 10.1016/0022-5193(78)90168-6
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Hestrin, S., and Armstrong, W. E. (1996). Morphology and physiology of cortical neurons in layer I. J. Neurosci. 16, 5290–5300.
Hodgkin, A. L., and Huxley, A. F. (1952). A quantitative description of membrane current and its application to conduction and excitation in nerve. J. Physiol. 117, 500–544. doi: 10.1113/jphysiol.1952.sp004764
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Jung, P., and Shuai, J. W. (2001). Optimal sizes of ion channel clusters. Europhys. Lett. 56, 29–35. doi: 10.1209/epl/i2001-00483-y
Levin, J. E., and Miller, J. P. (1996). Broadband neural encoding in the cricket cercal sensory system enhanced by stochastic resonance. Nature 380, 165–168. doi: 10.1038/380165a0
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Liberman, M. C., and Oliver, M. E. (1984). Morphometry of intracellularly labeled neurons of the auditory nerve: correlations with functional properties. J. Comp. Neurol. 223, 163–176. doi: 10.1002/cne.902230203
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Liberti, M., Apollonio, F., Merla, C., and d'Inzeo, G. (2009b). Microdosimetry in the microwave range: a quantitative assessment at single cell level IEEE. Antennas Wirel. Propag. Lett. 8, 865–868. doi: 10.1109/LAWP.2009.2028045
Liberti, M., Apollonio, F., Paffi, A., Parazzini, M., Maggio, F., Novellino, T., et al. (2007). “Fundamental electrical quantities in deep brain stimulation: influence of domain dimensions and boundary conditions,” in Proceedings of the 29th IEEE EMBS Annual International Conference (Lyon), 6668–6671.
Liberti, M., Paffi, A., Maggio, F., De Angelis, A., Apollonio, F., and d'Inzeo, G. (2009a). “Channel noise enhances signal detectability in a model of acoustic neuron through the stochastic resonance paradigm,” in Proceedings of 31st IEEE EMBS Annual International Conference (Minneapolis, MN), 1525–1528.
Liu, F., Wang, J., and Wang, W. (1999). Frequency sensitivity in weak signal detection. Phys. Rev. E 59, 3453–3460. doi: 10.1103/PhysRevE.59.3453
Longtin, A. (1993). Stochastic resonance in neuron models. J. Stat. Phys. 70, 309–327. doi: 10.1007/BF01053970
Maggio, F., Liberti, M., Paffi, A., Apollonio, F., Parazzini, M., Ravazzani, P., et al. (2009). “A three dimensional electromagnetic model for the DBS application,” in Proceedings of the 4th International Conference. IEEE EMBS on Neural Engineering (Antalya), 22–25.
Maggio, F., Pasciuto, T., Paffi, A., Apollonio, F., Parazzini, M., Ravazzani, P., et al. (2010). “Micro vs macro electrode DBS stimulation: a dosimetric study,” in Engineering in Medicine and Biology Society (EMBC) Annual International Conference of the IEEE (Buenos Aires), 2057–2060.
Mainen, Z. F., Joerges, J., Huguenard, J., and Sejnowski, T. (1995). A model of spike initiation in neocortical pyramidal neurons. Neuron 15, 1427–1439. doi: 10.1016/0896-6273(95)90020-9
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Mainen, Z. F., and Sejnowski, T. (1995). Reliability of spike timing in neocortical neurons. Science 268, 1503–1508. doi: 10.1126/science.7770778
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Manwani, A., Steinmetz, P. N., and Koch, C. (2002). The impact of spike timing variability on the signal-encoding performance of neural spiking models. Neural Comput. 14, 347–367. doi: 10.1162/08997660252741158
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Marchionni, I., Paffi, A., Pellegrino, M., Liberti, M., Apollonio, F., Abeti, R., et al. (2006). Comparison between low-level 50 Hz and 900 MHz electromagnetic stimulation on single channel ionic currents and on firing frequency in dorsal root ganglion isolated neurons. Biochim. Biophys. Acta Biomembr. 1758, 597–605. doi: 10.1016/j.bbamem.2006.03.014
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
McDonnell, M. D., and Abbott, D. (2009). What is stochastic resonance? Definitions, misconceptions, debates, and its relevance to biology. PLoS Comput. Biol. 5:e1000348. doi: 10.1371/journal.pcbi.1000348
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Merla, C., Denzi, A., Paffi, A., Casciola, M., d'Inzeo, G., Apollonio, F., et al. (2012). Novel passive element circuits for microdosimetry of nanosecond pulsed electric fields. IEEE Trans. Biomed. Eng. 59, 2302–2311. doi: 10.1109/TBME.2012.2203133
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Merla, C., Liberti, M., Apollonio, F., Nervi, C., and d'Inzeo, G. (2010). A 3-D microdosimetric study on blood cells: a permittivity model of cell membrane and stochastic electromagnetic analysis. IEEE Trans. Microw. Theory Tech. 58, 691–698. doi: 10.1109/TMTT.2010.2040338
Merla, C., Paffi, A., Apollonio, F., Leveque, P., d'Inzeo, G., and Liberti, M. (2011). Microdosimetry for nanosecond pulsed electric field applications: a parametric study for a single cell. IEEE Trans. Biomed. Eng. 58, 1294–1302. doi: 10.1109/TBME.2010.2104150
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Mino, H., Rubinstein, J. T., Miller, C. A., and Abbas, P. J. (2004). Effects of electrode-to-fiber distance on temporal neural response with electrical stimulation. IEEE Trans. Biomed. Eng. 51, 13–20. doi: 10.1109/TBME.2003.820383
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Mino, H., Rubinstein, J. T., and White, J. A. (2002). Comparison of algorithms for the simulation of action potentials with stochastic sodium channels. Ann. Biom. Eng. 30, 578–583. doi: 10.1114/1.1475343
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Moretti, D., Garenne, A., Haro, E., Poulletier de Gannes, F., Lagroye, I., Veyret, B., et al. (2013). In-vitro exposure of neuronal networks to the GSM-1800 signal. Bioelectromagnetics 34, 571–578. doi: 10.1002/bem.21805
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Morse, R. P., and Roper, P. (2000). Enhanced coding in a cochlear-implant model using additive noise: aperiodic stochastic resonance with tuning. Phys. Rev. E 61, 5683–5692. doi: 10.1103/PhysRevE.61.5683
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Moss, F., Ward, L. M., and Sannita, W. G. (2004). Stochastic resonance and sensory information processing: a tutorial and review of application. Clin. Neurophysiol. 115, 267–281. doi: 10.1016/j.clinph.2003.09.014
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Okun, M. S., Gallo, B. V., Mandybur, G., Jagid, J., Foote, K. D., Revilla, F. J., et al. (2012). Subthalamic deep brain stimulation with a constant-current device in Parkinson's disease: an open-label randomised controlled trial. Lancet Neurol. 11, 140–149. doi: 10.1016/S1474-4422(11)70308-8
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Overstreet, C. K., Klein, J. D., and Helms Tillery, S. I. (2013). Computational modeling of direct neuronal recruitment during intracortical microstimulation in somatosensory cortex. J. Neural Eng. 10:066016. doi: 10.1088/1741-2560/10/6/066016
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Paffi, A., Apollonio, F., d'Inzeo, G., and Liberti, M. (2013b). Stochastic resonance induced by exogenous noise in a model of a neuronal network. Netw. Comput. Neural Syst. 24, 99–113. doi: 10.3109/0954898X.2013.793849
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Paffi, A., Apollonio, F., Puxeddu, M. G., Parazzini, M., d'Inzeo, G., Ravazzani, P., et al. (2013a). A numerical study to compare stimulations by intraoperative microelectrodes and chronic macroelectrodes in the DBS technique. Biomed Res. Int. 2013:262739. doi: 10.1155/2013/262739
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Paffi, A., Apollonio, F., Puxeddu, M. G., Parazzini, M., d'Inzeo, G., Ravazzani, P., et al. (2013c). “A dosimetric study comparing intra-operatory microelectrode and chronic macroelectrode in the DBS technique,” in Proceedings of the 6th International IEEE EMBS Conference on Neural Engineering (San Diego, CA), 1206–1209.
Paffi, A., Camera, F., Apollonio, F., d'Inzeo, G., and Liberti, M. (2015). Numerical characterization of intraoperative and chronic electrodes in deep brain stimulation. Front. Comput. Neurosci. 9:2. doi: 10.3389/fncom.2015.00002
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Paffi, A., Giannì, M., Maggio, F., Liberti, M., Apollonio, F., and d'Inzeo, G. (2007). “Effects of an exogenous noise on a realistic network model: encoding of an EM signal,” Proceedings of the 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS) (Lyon), 2404–2407.
Paffi, A., Liberti, M., Apollonio, F., Giann,ì, M., and d'Inzeo, G. (2006). “Effects of exogenous noise in a silent neuron model: firing induction and EM signal detection,” in Proceedings of the 28th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS) (New York, NY), 4183–4186.
Peckham, P. H., and Knutson, J. S. (2005). Functional electrical stimulation for neuromuscular applications. Annu. Rev. Biomed. Eng. 7, 327–360. doi: 10.1146/annurev.bioeng.6.040803.140103
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Platano, D., Mesirca, P., Paffi, A., Pellegrino, M., Liberti, M., Apollonio, F., et al. (2007). Acute exposure to 900 MHz CW and GSM-modulated radiofrequencies does not affect Ba2+ currents through voltage-gated calcium channels in rat cortical neurons. Bioelectromagnetics 28, 599–607. doi: 10.1002/bem.20345
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Pospischil, M., Toledo-Rodriguez, M., Monier, C., Piwkowska, Z., Bal, T., Frégnac, Y., et al. (2008). Minimal Hodgkin–Huxley type models for different classes of cortical and thalamic neurons. Biol. Cybern. 99, 427–441. doi: 10.1007/s00422-008-0263-8
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Rapp, M., Yarom, Y., and Segev, I. (1996). Modeling back propagating action potential in weakly excitable dendrites of neocortical pyramidal cells. Proc. Natl. Acad. Sci. U.S.A. 93, 11985–11990. doi: 10.1073/pnas.93.21.11985
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Rattay, F. (2000). Basics of hearing theory and noise in cochlear implants. Chaos Solitons Fractals 11, 1875–1884. doi: 10.1016/S0960-0779(99)00124-1
Rinzel, J., and Ermentrout, B. (1998). “Analysis of neural excitability and oscillations,” in Methods in Neuronal Modeling, eds C. Koch and I. Segev (Cambridge, MA: MIT Press), 135–169.
Ritter, J., Anniko, M., and Gerhardt, H. J. (1981). Some new aspects on damages in the organ of corti after pure tone exposure. Arch. Otorhinolaryngol. 232, 187–197. doi: 10.1007/BF00505037
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Rubinstein, J. T. (1995). Threshold Fluctuations in an N sodium channel model of the node of Ranvier. Biophys. J. 68, 779–785. doi: 10.1016/S0006-3495(95)80252-3
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Rubinstein, J. T., and Hong, R. (2003). Signal coding in cochlear implants: Exploiting stochastic effects of electrical stimulation. Ann. Otol. Rhinol. Laryngol. Suppl. 191, 14–19.
Rudolph, M., and Destexhe, A. (2001). Do neocortical pyramidal neurons display stochastic resonance? J. Comput. Neurosci. 11, 19–42. doi: 10.1023/A:1011200713411
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Schmid, G., Goychuk, I., and Hanggi, P. (2001). Stochastic resonance as a collective property of ion channel assemblies. Europhys. Lett. 56, 22–28. doi: 10.1209/epl/i2001-00482-6
Schneidman, E., Freedman, B., and Segev, I. (1998). Ion channel stochasticity may be critical in determining the reliability and precision of spike timing. Neural Comput. 10, 1679–1703. doi: 10.1162/089976698300017089
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Sengupta, B., Faisal, A. A., Laughlin, S. B., and Niven, J. E. (2013). The effect of cell size and channel density on neuronal information encoding and energy efficiency. J. Cereb. Blood Flow Metab. 33, 1465–1473. doi: 10.1038/jcbfm.2013.103
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Sento, S., and Ryugo, D. K. (1989). Endbulbs of held and spherical bushy cells in cats: morphological correlates with physiological properties. J. Comp. Neurol. 280, 553–562. doi: 10.1002/cne.902800406
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Sigworth, F. J. (1980). The variance of sodium current fluctuations at the node of Ranvier. J. Physiol. 307, 97–129. doi: 10.1113/jphysiol.1980.sp013426
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Steinmetz, P. N., Manwani, A., and Koch, C. (2000). Subthreshold voltage noise due to channel fluctuations in active neuronal membrane. J. Comput. Neurosci. 9, 133–148. doi: 10.1023/A:1008967807741
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Stocks, N. G., Allingham, D., and Morse, R. P. (2002). The application of suprathreshold stochastic resonance to cochlear implant coding. Fluct. Noise Lett. 3, L169–L181. doi: 10.1142/S0219477502000774
Stodilka, R. Z., Modolo, J., Prato, F. S., Robertson, J. A., Cook, C., Patrick, J., et al. (2011). Pulsed magnetic field exposure induces lasting changes in neural network dynamics. Neurocomputing 74, 2164–2175. doi: 10.1016/j.neucom.2011.01.025
Tateno, K., Hayashi, H., and Ishizuka, S. (1998). Complexity of spatiotemporal activity of a neural network model which depends on the degree of synchronization. Neural Netw. 11, 985–1003. doi: 10.1016/S0893-6080(98)00086-0
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Tsong, T. Y., and Astumian, R. D. (1987). Electroconformational coupling and membrane protein function. Prog. Biophys. Mol. Biol. 50, 1–45. doi: 10.1016/0079-6107(87)90002-2
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
White, J. A., Rubinstein, J. T., and Kay, A. R. (2000). Channel noise in neurons. Trends Neurosci. 23, 131–137. doi: 10.1016/S0166-2236(99)01521-0
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Wilson, B. S., Finley, C. C., Lawson, D. T., Wolford, R. D., Eddington, D. K., and Rabinowitz, W. M. (1991). Better speech recognition with cochlear implants. Nature 352, 236–238. doi: 10.1038/352236a0
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Woo, J., Miller, C. A., and Abbas, P. J. (2010). The dependence of auditory nerve rate adaptation on electric stimulus parameters, electrode position, and fiber diameter: a computer model study. J. Assoc. Res. Otolaryngol. 11, 283–296. doi: 10.1007/s10162-009-0199-2
PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar
Yu, Y., Liu, F., and Wang, W. (2001a). Frequency sensitivity in Hodgkin-Huxley systems. Biol. Cyber. 84, 227–235. doi: 10.1007/s004220000203
Keywords: single neuron, HH model, electric stimulation, exogenous noise, stochastic resonance, signal detection
Citation: Paffi A, Camera F, Apollonio F, d'Inzeo G and Liberti M (2015) Restoring the encoding properties of a stochastic neuron model by an exogenous noise. Front. Comput. Neurosci. 9:42. doi: 10.3389/fncom.2015.00042
Received: 02 August 2014; Accepted: 19 March 2015;
Published: 06 May 2015.
Edited by:
Misha Tsodyks, Weizmann Institute of Science, IsraelReviewed by:
Asaf Gal, Technion - Israel Institute of Technology, IsraelGianluigi Mongillo, Paris Descartes University, France
Lee DeVille, University of Illinois, USA
Copyright © 2015 Paffi, Camera, Apollonio, d'Inzeo and Liberti. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alessandra Paffi, Department of Information Engineering, Electronics and Telecommunications, Sapienza University of Rome, Via Eudossiana 18, 00184 Rome, Italy, paffi@die.uniroma1.it