Recurrent discharge in a reduced model of cat spinal motoneuron (Balbi et al, 2013)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:151443
Following a distal stimulation of a motor fibre, only a fraction of spinal motoneurons are able to produce a re-excitation of the initial segment leading to an orthodromically conducted action potential, known as recurrent discharge. In order to show the reciprocal interplay of the axonal initial segment and the soma leading to recurrent discharge in detail, a reduced model of a cat spinal motoneuron was developed.
Reference:
1 . Balbi P, Martinoia S, Colombo R, Massobrio P (2014) Modelling recurrent discharge in the spinal a-motoneuron: reappraisal of the F wave. Clin Neurophysiol 125:427-9 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell;
Brain Region(s)/Organism:
Cell Type(s): Spinal cord lumbar motor neuron alpha ACh cell;
Channel(s): I Na,p; I Na,t; I K; I K,Ca;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Recurrent Discharge;
Implementer(s): Balbi, Pietro [piero.balbi at fsm.it];
Search NeuronDB for information about:  Spinal cord lumbar motor neuron alpha ACh cell; I Na,p; I Na,t; I K; I K,Ca;
TITLE Kdr.mod   delayed rectifier potassium channel
 
COMMENT
This is the original Hodgkin-Huxley treatment for the set of potassium channel found
in the squid giant axon membrane.
Some parameters have been changed to correspond to McIntyre and Grill (2002) "Extracellular
stimulation of central neurons"
Author: Balbi
ENDCOMMENT
 
UNITS {
        (mA) = (milliamp)
        (mV) = (millivolt)
		(S) = (siemens)
}
 
NEURON {
        SUFFIX Kdr
        USEION k READ ek WRITE ik
        RANGE gkmax, gk
        RANGE ninf, ntau
}
 
PARAMETER { : value for soma
        gkmax = .5 (S/cm2)	  <0,1e9>
}
 
STATE {
        n
}
 
ASSIGNED {
        v (mV)
        ek (mV)

		gk (S/cm2)
        ik (mA/cm2)
        ninf
		ntau (ms)
}
 
BREAKPOINT {
        SOLVE states METHOD cnexp
        gk = gkmax*n*n*n*n
		ik = gk*(v - ek)
} 
 
INITIAL {
	rates(v)
	n = ninf
}

DERIVATIVE states {
        rates(v)
        n' = (ninf-n)/ntau
}
 
PROCEDURE rates(v(mV)) {  
		:Call once from HOC to initialize inf at resting v.
        LOCAL  alpha, beta, sum

UNITSOFF
                :"n" potassium activation system
        ntau = 5/(exp((50+v)/40)+exp(-(50+v)/50))
        ninf = 1/(1+exp((38+v)/-15))
}
 
FUNCTION vtrap(x,y) {  :Traps for 0 in denominator of rate eqns.
        if (fabs(x/y) < 1e-6) {
                vtrap = y*(1 - x/y/2)
        }else{
                vtrap = x/(exp(x/y) - 1)
        }
}
 
UNITSON

Loading data, please wait...