Spontaneous firing caused by stochastic channel gating (Chow, White 1996)

 Download zip file 
Help downloading and running models
NEURON implementation of model of stochastic channel gating, resulting in spontaneous firing. Qualitatively reproduces the phenomena described in the reference.
1 . Chow CC, White JA (1996) Spontaneous action potentials due to channel fluctuations. Biophys J 71:3013-21 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Channel/Receptor;
Brain Region(s)/Organism:
Cell Type(s):
Channel(s): I K; I Sodium;
Gap Junctions:
Simulation Environment: NEURON;
Model Concept(s): Action Potential Initiation; Activity Patterns; Ion Channel Kinetics; Temporal Pattern Generation; Axonal Action Potentials; Action Potentials;
Implementer(s): Carnevale, Ted [Ted.Carnevale at Yale.edu]; Donnelly, David [david.donnelly at yale.edu];
Search NeuronDB for information about:  I K; I Sodium;
This NEURON implementation qualitatively reproduces the phenomena described in 
Chow, C.C. and White, J.A.
Spontaneous action potentials due to channel fluctuations.
Biophysical Journal 71:3013-3021, 1996

An exact reproduction of their figures is not possible for the following reasons:
--differences in the pseudorandom number generators
--differences in numerical integration methods
--differences in the algorithms used to govern channel state 

When first loaded, this code automatically emulates part of Figure 2.
With minor modifications, it can be used to generate data for Figures 1, 3, and 4; 
however, note that this will require prolonged, repeated simulations.

On a 2 GHz PC, ~1600 seconds are needed to simulate 30 seconds of activity 
in a 100 um2 patch of membrane.


////////// Specification of Biological Properties //////////

load_file("cell100.hoc")	// basic passive model
load_file("channels.ses")	// voltage-gated stochastic channels 
                                // implemented with the Channel Builder

Ionic        Conductance  channel  Single channel
conductance  density      density  conductance
gna          0.12 S/cm2   60/um2   2e-5 uS
gk           0.036 S/cm2  18/um2   2e-5 uS

Area = 100 // um2

// argument should be desired area in um2
proc setarea() {
	L = Area/diam/PI
	DDINa[0].Nsingle = 60 * Area
	ddhhkp[0].Nsingle = 18 * Area


xpanel("Set Area", 0)
xvalue("Area (um2)","Area", 1,"setarea()", 0, 1 )

////////// Simulation Control and Graphical Displays //////////

load_file("rig.ses")	// run control and v vs. t

// Emulate Fig. 2


print "Emulating Fig. 2"
print "This takes about 100 seconds on a 2 GHz PC."


Loading data, please wait...