Regulation of motoneuron excitability by KCNQ/Kv7 modulators (Lombardo & Harrington 2016)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:217882
" ... Computer simulations confirmed that pharmacological enhancement of KCNQ/Kv7 channel (M current) activity decreases excitability and also suggested that the effects of inhibition of KCNQ/Kv7 channels on the excitability of spinal MNs do not depend on a direct effect in these neurons but likely on spinal cord synaptic partners. These results indicate that KCNQ/Kv7 channels have a fundamental role in the modulation of the excitability of spinal MNs acting both in these neurons and in their local presynaptic partners. ..."
Reference:
1 . Lombardo J, Harrington MA (2016) Nonreciprocal mechanisms in up- and downregulation of spinal motoneuron excitability by modulators of KCNQ/Kv7 channels. J Neurophysiol 116:2114-2124 [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: Spinal motoneuron;
Cell Type(s): Spinal cord lumbar motor neuron alpha ACh cell;
Channel(s): I Potassium; I K; I Na,t; I M;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Axonal Action Potentials;
Implementer(s): Lombardo, Joseph [josslomb at gmail.com];
Search NeuronDB for information about:  Spinal cord lumbar motor neuron alpha ACh cell; I Na,t; I K; I M; I Potassium;
/
LombardoHarrington2016
readme.html
Gfluctdv.mod *
ghchan.mod *
kca2.mod *
kdrRL.mod *
Km.mod *
kv1_gp.mod
L_Ca.mod *
mAHP.mod *
na3n.mod
na3rp.mod *
naps.mod *
napsi.mod *
buttons.png
ctrl.png
FR3cablepas.hoc
FRMot3dendNaHH.hoc
gKm0.png
GraphicsKmModulators.hoc
ModifiedFRMotoneuron.hoc
mosinit.hoc
retigabine.png
standard_0.hoc
Tools.ses
XE991.png
                            
:   KV1_GP.MOD
:
:   Kv1.2 channel model using HH-type activation/inactivation
:
:   3/2003
:   Josh Held

NEURON {
	SUFFIX kv1_gp
	USEION k READ ek WRITE ik
	RANGE th, tm, ik, hinf, minf, g, gbar
	GLOBAL p
	GLOBAL vhm, vcm
	GLOBAL vhh, vch
	GLOBAL Cth, vhth, ath, bth, th0
	GLOBAL tm0, Ctm, vhtm, vctm
	GLOBAL th90
	GLOBAL Cq10
}

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

PARAMETER {
	gbar = 1	(S/cm2)
	ek		(mV)
	vhm = -27	(mV)
	vhh = -33.477	(mV)
	vcm = -16	(mV)
	vch = 21.5	(mV)
	Cth = 548.67	(ms)
	vhth = -0.956	(mV)
	ath = 29.013	(mV)
	bth= 100	(mV)
	th0 = 779	(ms)
	tm0 = 3.4	(ms)
	Ctm = 89.2	(ms)
	vhtm = -34.3	(mV)
	vctm = 30.1	(mV)
	p = 0.004
	celsius		(degC)
	Cq10 = 3
}

ASSIGNED {
	v	(mV)
	minf
	hinf
	tm	(ms)
	th	(ms)
	th90	(ms)
	ik	(mA/cm2)
	g	(S/cm2)
}

STATE {
	m
	h
}

BREAKPOINT {
	SOLVE states METHOD cnexp
	g = gbar * (m^2) * h 
	ik = g * (v - ek) 
}

DERIVATIVE states{
	values()
	m' = (minf - m)/tm
	h' = (hinf - h)/th
}

INITIAL {
	values()
	m = minf
	h = hinf
}

PROCEDURE values() {LOCAL q10
	q10 = Cq10^((celsius-22 (degC))/10 (degC))
	minf = 1/(1 + exp((v - vhm)/vcm))
	tm = (1/q10)*(tm0 + Ctm*exp(-((v-vhtm)/vctm)^2))

	hinf = (1-p)/(1 + exp((v - vhh)/vch)) + p
	th = (1/q10)*((Cth/(exp((v-vhth)/ath) + exp(-(v-vhth)/bth))) + th0)
	th90 = (1/q10)*((Cth/(exp((-90-vhth)/ath) + exp(-(-90-vhth)/bth))) + th0)
}






Loading data, please wait...