Nav1.6 sodium channel model in globus pallidus neurons (Mercer et al. 2007)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:105385
Model files for the paper Mercer JN, Chan CS, Tkatch T, Held J, Surmeier DJ. Nav1.6 sodium channels are critical to pacemaking and fast spiking in globus pallidus neurons.,J Neurosci. 2007 Dec 5;27(49):13552-66.
Reference:
1 . Mercer JN, Chan CS, Tkatch T, Held J, Surmeier DJ (2007) Nav1.6 sodium channels are critical to pacemaking and fast spiking in globus pallidus neurons. J Neurosci 27:13552-66 [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): Globus pallidus neuron;
Channel(s): I Na,p; I Na,t; I K; I h; I K,Ca; I Sodium; I Calcium; I Potassium;
Gap Junctions:
Receptor(s):
Gene(s): Nav1.6 SCN8A;
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Action Potential Initiation; Action Potentials; Parkinson's;
Implementer(s): Held, Joshua [j-held at northwestern.edu];
Search NeuronDB for information about:  I Na,p; I Na,t; I K; I h; I K,Ca; I Sodium; I Calcium; I Potassium;
:KCNQ_GP.MOD
:
: 11/3/2003

NEURON {
        SUFFIX kcnq_gp
        USEION k READ ek WRITE ik
        RANGE ik, ek, g, gbar
        GLOBAL a0, b0, ah, bh, ac, bc, q10v
}

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

PARAMETER {
        gbar    = 1		(S/cm2)
        ek			(mV)
        a0      = .018		(/ms)	: parameters for alpha and beta
        b0      = .01		(/ms)
        ah      = -43.367	(mV)
        bh      = -43.367	(mV)
        ac      = 9.7054	(mV)
        bc      = -9.7054	(mV)
        q10v    = 3
        celsius			(degC)
}

ASSIGNED {
        v	(mV)
        g	(S/cm2)
        ik	(mA/cm2)
        alpha   (/ms)
        beta    (/ms)
}

STATE {
	c
	o
}

INITIAL {
    SOLVE kin STEADYSTATE sparse
}

BREAKPOINT {
        SOLVE kin METHOD sparse
        g = gbar*o
        ik = g*(v-ek) 
}

KINETIC kin {
        rates(v)
        ~ c <-> o       (alpha, beta)
        CONSERVE c + o = 1
}

PROCEDURE rates(v(mV)) {
        LOCAL qv
        qv = q10v^((celsius-22 (degC))/10 (degC))
        alpha = a0*qv / (1 + exp(-(v-ah)/ac))
        beta = b0*qv / (1 + exp(-(v-bh)/bc))
}

Loading data, please wait...