Simulating ion channel noise in an auditory brainstem neuron model (Schmerl & McDonnell 2013)

 Download zip file 
Help downloading and running models
" ... Here we demonstrate that biophysical models of channel noise can give rise to two kinds of recently discovered stochastic facilitation effects in a Hodgkin-Huxley-like model of auditory brainstem neurons. The first, known as slope-based stochastic resonance (SBSR), enables phasic neurons to emit action potentials that can encode the slope of inputs that vary slowly relative to key time constants in the model. The second, known as inverse stochastic resonance (ISR), occurs in tonically firing neurons when small levels of noise inhibit tonic firing and replace it with burstlike dynamics. ..." Preprint available at
1 . Schmerl BA, McDonnell MD (2013) Channel noise induced stochastic facilitation in an auditory brainstem neuron model Physical Review E 88:052722 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell; Channel/Receptor;
Brain Region(s)/Organism: Auditory brainstem;
Cell Type(s): Cochlear nucleus bushy GLU cell; CN stellate cell; Ventral cochlear nucleus T stellate (chopper) neuron; Hodgkin-Huxley neuron;
Channel(s): I h; I Sodium; I Potassium;
Gap Junctions:
Gene(s): Kv1.1 KCNA1; Kv3.1 KCNC1;
Simulation Environment: MATLAB;
Model Concept(s): Bursting; Ion Channel Kinetics; Action Potentials; Methods; Noise Sensitivity; Bifurcation; Audition;
Search NeuronDB for information about:  Cochlear nucleus bushy GLU cell; I h; I Sodium; I Potassium;
%Rothman and Manis Type II model
ICs = [-63.6269,0.5122,0.6618,0.0077,0.0011,0.0251,0.4430,0.1458];

C = 12;

gP = 45;
gR = 40; 
glk = 4; 

EP = -70;
ER = -43; 
Elk = -65;

gs = [gNa,gLT,gN,gP,gR,glk];
Es = [ENa,ELT,EN,EP,ER,Elk];

NumChannelTypes = 5;
ActivationVarsPerChannel = [2,2,1,1,1];
NumActivationVars = sum(ActivationVarsPerChannel);
NumGatesPerActivationVariable{1} = [3,1];
NumGatesPerActivationVariable{2} = [4,1];
NumGatesPerActivationVariable{3} = [2];
NumGatesPerActivationVariable{4} = [1];
NumGatesPerActivationVariable{5} = [1];

SpikeThreshold = -35;