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

 Download zip file 
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: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 GABA 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 GABA cell; I Na,p; I p,q; I A; I K; I M; I K,Ca; I Sodium; I Calcium; I Potassium;
TITLE BK calcium-activated potassium current
: Calcium activated K channel.

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 {
	(molar) = (1/liter)
}

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


INDEPENDENT {t FROM 0 TO 1 WITH 1 (ms)}

NEURON {
	SUFFIX KC3cvode
	USEION ca READ cai
	USEION k WRITE ik
	RANGE gkbar,gk,ik
}


PARAMETER {
	celsius=37	(degC)
	v		(mV)
	gkbar=.08	(mho/cm2)	: Maximum Permeability
	cai = .04e-3	(mM)
	ekcvode  = -85	(mV)
      mon = 1
	zon = 1	
}


ASSIGNED {
	ik		(mA/cm2)
	minf
	zinf
	gk
	tau 
	q10 
	alpha 
	beta 
	sum 
}

STATE {	m z }		: fraction of open channels

BREAKPOINT {
	SOLVE state METHOD cnexp
:	gk = gkbar*m*z*z
	ik = gkbar*m*z*z*(v - ekcvode)
}

:UNITSOFF

INITIAL {
	
	m = minf
	z = zinf
}


DERIVATIVE state { :callable from hoc
	
	alpha = 400/(cai*1000)
	beta = 0.11/exp((v-35)/14.9)
	
	zinf = 1/(1+alpha)
      tau= 10
	z' = zon * (zinf-z)/tau      
	
	minf = 7.5/(7.5+beta)
      tau= 1/(7.5 + beta)
	m' = mon * (minf-m)/tau      

}

:UNITSON

Loading data, please wait...