A simplified cerebellar Purkinje neuron (the PPR model) (Brown et al. 2011)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:126637
These models were implemented in NEURON by Sherry-Ann Brown in the laboratory of Leslie M. Loew. The files reproduce Figures 2c-f from Brown et al, 2011 "Virtual NEURON: a Strategy For Merged Biochemical and Electrophysiological Modeling".
Reference:
1 . Brown SA, Moraru II, Schaff JC, Loew LM (2011) Virtual NEURON: a strategy for merged biochemical and electrophysiological modeling. J Comput Neurosci 31(2):385-400 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell; Dendrite;
Brain Region(s)/Organism: Cerebellum;
Cell Type(s): Cerebellum Purkinje cell;
Channel(s): I Na,p; I p,q; I A; I K; I M; I K,Ca; I Sodium; I Calcium; I Potassium;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Dendritic Action Potentials; Simplified Models; Active Dendrites; Influence of Dendritic Geometry; Detailed Neuronal Models; Intrinsic plasticity; Methods; Synaptic Integration;
Implementer(s): Brown, Sherry-Ann [sabrown at student.uchc.edu];
Search NeuronDB for information about:  Cerebellum Purkinje cell; I Na,p; I p,q; I A; I K; I M; I K,Ca; I Sodium; I Calcium; I Potassium;
TITLE slowly inactivating K current
 
COMMENT
This is adapted for CVODE (by Brown et al 2010) from the original file published in "An Active Membrane Model of the Cerebellar Purkinje Cell 1. Simulation of Current Clamp in Slice".
ENDCOMMENT
 
UNITS {
        (mA) = (milliamp)
        (mV) = (millivolt)
}
 
NEURON {
        SUFFIX KDcvode
        USEION k  WRITE ik
        RANGE  gkbar, ik, gk
} 
 
INDEPENDENT {t FROM 0 TO 1 WITH 1 (ms)}
 
PARAMETER {
        v (mV)
        celsius = 37 (degC)
        mon = 1
	  hon = 1
        gkbar = .0045 (mho/cm2)
        ekcvode = -85 (mV)
}
 
STATE {
        m h
}
 
ASSIGNED {
        ik (mA/cm2)
        gk minf hinf tau q10 alpha beta sum 
}
 
BREAKPOINT {
        SOLVE state METHOD cnexp
        gk = gkbar * m*h
	  ik = gk* (v-ekcvode)
}
 
UNITSOFF
 
INITIAL {
	
	m = minf
	h = hinf
}


DERIVATIVE state {  :Computes rate and other constants at current v.
                      :Call once from HOC to initialize inf at resting v.
        
        q10 = 3^((celsius - 37)/10)
        
                :"m" potassium activation system

        alpha = 8.5/(1+exp((v+17)/(-12.5)))
        beta =  35/(1+exp((v+99)/14.5))
        sum = alpha + beta
        minf = alpha/sum
	  tau= 10/(q10 *sum)
        m' = mon * (minf-m)/tau      

                :"h" potassium inactivation system

        alpha = 0.0015/(1+exp((v+89)/8))
        beta = 0.0055/(1+exp((v+83)/(-8)))
        sum = alpha + beta
        hinf = alpha/sum
	  tau = 1/(q10 * sum *1.6)
     	  h' = hon * (hinf-h)/tau
}

 
UNITSON


Loading data, please wait...