On stochastic diff. eq. models for ion channel noise in Hodgkin-Huxley neurons (Goldwyn et al. 2010)

 Download zip file 
Help downloading and running models
Accession:128502
" ... We analyze three SDE models that have been proposed as approximations to the Markov chain model: one that describes the states of the ion channels and two that describe the states of the ion channel subunits. We show that the former channel-based approach can capture the distribution of channel noise and its effect on spiking in a Hodgkin-Huxley neuron model to a degree not previously demonstrated, but the latter two subunit-based approaches cannot. ..."
Reference:
1 . Goldwyn JH, Imennov NS, Famulare M, Shea-Brown E (2011) Stochastic differential equation models for ion channel noise in Hodgkin-Huxley neurons. Phys Rev E Stat Nonlin Soft Matter Phys 83:041908 [PubMed]
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:
Cell Type(s): Squid axon;
Channel(s): I Sodium; I Potassium;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: FORTRAN;
Model Concept(s): Ion Channel Kinetics; Action Potentials; Methods; Noise Sensitivity;
Implementer(s): Goldwyn, Joshua [jhgoldwyn at gmail.com];
Search NeuronDB for information about:  I Sodium; I Potassium;
!***********************************************************************************
MODULE ParameterModule

IMPLICIT NONE

SAVE

! HH Parameters (Hodgkin and Huxley, 1952)
DOUBLE PRECISION, PARAMETER :: C = 1 ! muF /cm^2
DOUBLE PRECISION, PARAMETER :: gNa = 120 ! mS/cm^2
DOUBLE PRECISION, PARAMETER :: ENa = 115 !  mV  ! changed from 120 to 115 because of apparent typo in Izhikevich book
DOUBLE PRECISION, PARAMETER :: DNa = 60. ! Channels/mu m^2 Channel Density
DOUBLE PRECISION, PARAMETER :: gK = 36 ! mS/cm^2
DOUBLE PRECISION, PARAMETER :: EK = -12 !  mV
DOUBLE PRECISION, PARAMETER :: DK = 18. ! Channels/mu m^2 Channel Density
DOUBLE PRECISION, PARAMETER :: gL = 0.3 ! mS / cm^2
DOUBLE PRECISION, PARAMETER :: EL = 10.6 ! mV
DOUBLE PRECISION, PARAMETER :: Erest = 0.
DOUBLE PRECISION, PARAMETER :: oned = 1.0
DOUBLE PRECISION, PARAMETER :: spikethreshold = 60. ! Threshold for spike count (mV)
DOUBLE PRECISION, PARAMETER :: minisi = 2.0 ! Minimum ISI for spike count (ms)
DOUBLE PRECISION, PARAMETER :: maxisi = 2000. ! Maximum ISI to kill program (ms)
DOUBLE PRECISION, PARAMETER :: PI=3.141592653589793238462643	! PI (Find a better way to compute this?)      

END MODULE ParameterModule
!***********************************************************************************

Loading data, please wait...