Spike frequency adaptation in spinal sensory neurones (Melnick et al 2004)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:62284
Using tight-seal recordings from rat spinal cord slices, intracellular labelling and computer simulation, we analysed the mechanisms of spike frequency adaptation in substantia gelatinosa (SG) neurones. Adapting-firing neurones (AFNs) generated short bursts of spikes during sustained depolarization and were mostly found in lateral SG. ... Ca2 + -dependent conductances do not contribute to adapting firing. Transient KA current was small and completely inactivated at resting potential suggesting that adapting firing was mainly generated by voltage-gated Na+ and delayed-rectifier K+ (KDR ) currents. ... Computer simulation has further revealed that down-regulation of Na+ conductance represents an effective mechanism for the induction of firing adaptation. It is suggested that the cell-specific regulation of Na+ channel expression can be an important factor underlying the diversity of firing patterns in SG neurones. See paper for more and details.
Reference:
1 . Melnick IV, Santos SF, Safronov BV (2004) Mechanism of spike frequency adaptation in substantia gelatinosa neurones of rat. J Physiol 559:383-95 [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;
Brain Region(s)/Organism:
Cell Type(s):
Channel(s):
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Short-term Synaptic Plasticity; Spike Frequency Adaptation;
Implementer(s): Safronov, Boris [safronov at ibmc.up.pt];
TITLE HH sodium channel
: Hodgkin - Huxley squid sodium channel

: The model used in Safronov et al. 2000
:
: 5/17/2017  Revised by N.T. Carnevale for the sake of conceptual clarity
: and to facilitate attributed reuse.
: In this version, the reference temperature is 23 deg C
: and the value assigned to celsius is the actual operating temperature
: in degrees celsius.

NEURON {
	SUFFIX SS
	USEION na READ ena WRITE ina
	RANGE gnabar, ina
	GLOBAL inf
}

UNITS {
	(mA) = (milliamp)
	(mV) = (millivolt)
}

INDEPENDENT {t FROM 0 TO 1 WITH 1 (ms)}
PARAMETER {
	v (mV)
:	celsius = 6.3	(degC)
	celsius = 23 (degC) : actual operating temperature
	dt (ms)
	gnabar=0 (mho/cm2) <0,1e9>
	ena = 53 (mV)
}
STATE {
	m h
}
ASSIGNED {
	ina (mA/cm2)
	inf[2]
}
LOCAL	fac[2]

INITIAL {
	rate(v*1(/mV))
	m = inf[0]
	h = inf[1]
}

BREAKPOINT {
	SOLVE states
	ina = gnabar*m*m*m*(v - ena)
}

PROCEDURE states() {	: exact when v held constant
	rate(v*1(/mV))
	m = m + fac[0]*(inf[0] - m)
	h = h + fac[1]*(inf[1] - h)
	VERBATIM
	return 0;
	ENDVERBATIM
}

UNITSOFF
FUNCTION alp(v(mV),i) { LOCAL a,b,c,q10 :rest = -70  order m,h
	v = v :convert to hh convention
:	q10 = 3^((celsius - 6.3)/10)
	q10 = 3^((celsius - 23)/10) : actual reference temperature
	if (i==0) {
		alp = q10* 1 *.182*expM1(-v - 45, 9)
	}else if (i==1){
		alp = q10*.024*expM1(-v - 50, 5)


	}
}

FUNCTION bet(v,i) { LOCAL a,b,c,q10 :rest = -70  order m,h
	v = v 
:	q10 = 3^((celsius - 6.3)/10)
	q10 = 3^((celsius - 23)/10) : actual reference temperature
	if (i==0) {
		bet = q10* 1 *.124*expM1(v + 45, 9)
	}else if (i==1){
		bet = q10*.0091*expM1(v + 75, 5)
	}
}

FUNCTION expM1(x,y) {
	if (fabs(x/y) < 1e-6) {
		expM1 = y*(1 - x/y/2)
	}else{
		expM1 = x/(exp(x/y) - 1)
	}
}

PROCEDURE rate(v) {LOCAL a, b, tau :rest = -70
	TABLE inf, fac DEPEND dt, celsius FROM -150 TO 100 WITH 200
	FROM i=0 TO 1 {
		a = alp(v,i)  b=bet(v,i)
		tau = 1/(a + b)
		if (i==0) {		
		inf[i] = a/(a+b)
	}else if (i==1) {
		inf[i] = 1/(1+exp((v+75)/9))
	} 
		fac[i] = (1 - exp(-dt/tau))
	}
}
UNITSON