Synaptic integration in tuft dendrites of layer 5 pyramidal neurons (Larkum et al. 2009)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:124043
Simulations used in the paper. Voltage responses to current injections in different tuft locations; NMDA and calcium spike generation. Summation of multiple input distribution.
Reference:
1 . Larkum ME, Nevian T, Sandler M, Polsky A, Schiller J (2009) Synaptic integration in tuft dendrites of layer 5 pyramidal neurons: a new unifying principle. Science 325:756-60 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell; Synapse; Dendrite;
Brain Region(s)/Organism:
Cell Type(s): Neocortex V1 L6 pyramidal corticothalamic GLU cell;
Channel(s): I L high threshold; I p,q; I A; I K,leak; I K,Ca; I Sodium;
Gap Junctions:
Receptor(s): GabaA; AMPA; NMDA;
Gene(s):
Transmitter(s): Gaba; Glutamate;
Simulation Environment: NEURON;
Model Concept(s): Dendritic Action Potentials; Active Dendrites; Detailed Neuronal Models; Synaptic Integration;
Implementer(s): Polsky, Alon [alonpol at tx.technion.ac.il];
Search NeuronDB for information about:  Neocortex V1 L6 pyramidal corticothalamic GLU cell; GabaA; AMPA; NMDA; I L high threshold; I p,q; I A; I K,leak; I K,Ca; I Sodium; Gaba; Glutamate;
/
larkumEtAl2009_2
readme.html
ampa.mod
cad2.mod
glutamate.mod *
h.mod *
h2.mod
hh3.mod *
ih.mod
it2.mod *
kap.mod
kca.mod *
kdf.mod
Kdr.mod *
kdr2.mod *
km.mod *
SlowCa.mod *
0.50764
0.55472
070603c2.cll
apic.ses
apical_simulation.hoc
layerV.cll
mosinit.hoc
screenshot1.jpg
screenshot2.jpg
screenshot3.jpg
screenshot4.jpg
screenshot5.jpg
screenshot6.jpg
                            
TITLE I-h channel from Kamilla Angelo, Michael London,Soren R. Christensen, and Michael Hausser

NEURON {
	SUFFIX h
	NONSPECIFIC_CURRENT i
        RANGE ghbar
        GLOBAL qinf,tauh,ratetau
}

UNITS {
	(mA) = (milliamp)
	(mV) = (millivolt)

}

PARAMETER {
	v 		(mV)
      eh  =-40	(mV)        
	ghbar=.0001 	(mho/cm2)
	ratetau=1
}

STATE {
        q
}

ASSIGNED {
	i (mA/cm2)
      qinf      
      tauh
}

INITIAL {
	rate(v)
	q=qinf
}


BREAKPOINT {
	SOLVE states METHOD cnexp
	i = q*ghbar*(v-eh)

}

DERIVATIVE states {  
	rate(v)
      q' =  (qinf - q)/tauh
}

PROCEDURE rate(v (mV)) {
	qinf=1/( 1+exp((90+v)/9.67) )
	:tauh=ratetau*1/0.00062*( exp((v+68)/-22) + exp((v+68)/7.14) )
      tauh=ratetau*1/(0.02*( exp((v+90)/-22) + exp((v+90)/22) ))

}















Loading data, please wait...