A set of reduced models of layer 5 pyramidal neurons (Bahl et al. 2012)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:146026
These are the NEURON files for 10 different models of a reduced L5 pyramidal neuron. The parameters were obtained by automatically fitting the models to experimental data using a multi objective evolutionary search strategy. Details on the algorithm can be found at http://www.g-node.org/emoo and in Bahl et al. (2012).
Reference:
1 . Bahl A, Stemmler MB, Herz AV, Roth A (2012) Automated optimization of a reduced layer 5 pyramidal cell model based on experimental data. J Neurosci Methods 210:22-34 [PubMed]
Citations  Citation Browser
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:
Cell Type(s): Neocortex U1 L5B pyramidal pyramidal tract GLU cell;
Channel(s): I Na,p; I Na,t; I K; I M; I h; I K,Ca; I Calcium; I A, slow;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Action Potential Initiation; Parameter Fitting; Simplified Models; Active Dendrites; Detailed Neuronal Models; Action Potentials; Methods; Calcium dynamics;
Implementer(s): Bahl, Armin [bahl at neuro.mpg.de];
Search NeuronDB for information about:  Neocortex U1 L5B pyramidal pyramidal tract GLU cell; I Na,p; I Na,t; I K; I M; I h; I K,Ca; I Calcium; I A, slow;
COMMENT

Deterministic model of kinetics and voltage-dependence of Ih-currents
in layer 5 pyramidal neuron, see Kole et al., 2006. Implemented by
Stefan Hallermann.

Added possibility to shift voltage activiation (vshift) and allowed access to gating variables, Armin Bahl 2009

Predominantly HCN1 / HCN2 

ENDCOMMENT

TITLE Ih-current

UNITS {
	(mA) = (milliamp)
	(mV) = (millivolt)
     (mM) = (milli/liter)

}

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

PARAMETER {
	dt 	   		(ms)
	v 	   		(mV)
        ehd=-47 		(mV) 				       
	gbar=0 (pS/um2)	
	gamma_ih	:not used
	seed		:not used
	vshift = 0
}


NEURON {
	SUFFIX ih
	NONSPECIFIC_CURRENT Iqq
	RANGE Iqq,gbar,vshift,ehd, qtau, qinf, gq
}

STATE {
	qq
}

ASSIGNED {
	Iqq (mA/cm2)
	qtau (ms)
	qinf
	gq	(pS/um2)
	
}

INITIAL {
	qq=alpha(v-vshift)/(beta(v-vshift)+alpha(v-vshift))

	qtau = 1./(alpha(v-vshift) + beta(v-vshift))
	qinf = alpha(v)/(alpha(v-vshift) + beta(v-vshift))
}

BREAKPOINT {
	SOLVE state METHOD cnexp
	
	qtau = 1./(alpha(v-vshift) + beta(v-vshift))
	qinf = alpha(v-vshift)/(alpha(v-vshift) + beta(v-vshift))
	
	gq = gbar*qq
	Iqq = (1e-4)*gq*(v-ehd)
	
}

FUNCTION alpha(v(mV)) {

	alpha = 0.001*6.43*(v+154.9)/(exp((v+154.9)/11.9)-1)
	: parameters are estimated by direct fitting of HH model to
        : activation time constants and voltage activation curve
        : recorded at 34C

}

FUNCTION beta(v(mV)) {
	beta = 0.001*193*exp(v/33.1)			
}

DERIVATIVE state {     : exact when v held constant; integrates over dt step
	qq' = (1-qq)*alpha(v-vshift) - qq*beta(v-vshift)
}