Tonic firing in substantia gelatinosa neurons (Melnick et al 2004)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:62285
Ionic conductances underlying excitability in tonically firing neurons (TFNs) from substantia gelatinosa (SG) were studied by the patch-clamp method in rat spinal cord slices. ... Suppression of Ca2+ and KCA currents ... did not abolish the basic pattern of tonic firing, indicating that it was generated by voltage-gated Na+ and K+ currents. ... on the basis of present data, we created a model of TFN and showed that Na+ and KDR currents are sufficient to generate a basic pattern of tonic firing. It is concluded that the balanced contribution of all ionic conductances described here is important for generation and modulation of tonic firing in SG neurons. See paper for more and details.
Reference:
1 . Melnick IV, Santos SF, Szokol K, Szûcs P, Safronov BV (2004) Ionic basis of tonic firing in spinal substantia gelatinosa neurons of rat. J Neurophysiol 91:646-55 [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):
Channel(s): I Na,t; I A; I K;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Activity Patterns;
Implementer(s): Safronov, Boris [safronov at ibmc.up.pt];
Search NeuronDB for information about:  I Na,t; I A; I K;
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

Loading data, please wait...