Interneuron Specific 3 Interneuron Model (Guet-McCreight et al, 2016)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:223031
In this paper we develop morphologically detailed multi-compartment models of Hippocampal CA1 interneuron specific 3 interneurons using cell current-clamp recordings and dendritic calcium imaging data. In doing so, we developed several variant models, as outlined in the associated README.html file.
Reference:
1 . Guet-McCreight A, Camiré O, Topolnik L, Skinner FK (2016) Using a Semi-Automated Strategy to Develop Multi-Compartment Models That Predict Biophysical Properties of Interneuron-Specific 3 (IS3) Cells in Hippocampus. eNeuro [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: Hippocampus;
Cell Type(s):
Channel(s): I Na,t; I Na,p; I A;
Gap Junctions:
Receptor(s):
Gene(s): Kv3.1 KCNC1; Kv2.1 KCNB1;
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Action Potentials; Active Dendrites; Activity Patterns; Depolarization block; Detailed Neuronal Models; Ion Channel Kinetics;
Implementer(s): Guet-McCreight, Alexandre [alexandre.guet.mccreight at mail.utoronto.ca];
Search NeuronDB for information about:  I Na,p; I Na,t; I A;
/
Guet-McCreightEtAl2016
SD
IKa.mod *
ingauss.mod *
Ksoma.mod *
Nap.mod *
Nasoma.mod *
init.hoc
IS3_SD_TopModel.hoc
                            
COMMENT

Potassium current for the soma
ENDCOMMENT
UNITS {
        (mA) = (milliamp)
        (mV) = (millivolt)
}
 
NEURON {
        SUFFIX Ksoma
        USEION k READ ek WRITE ik
        RANGE gksoma, ik
        GLOBAL ninf, nexp, ntau
}
 
INDEPENDENT {t FROM 0 TO 1 WITH 1 (ms)}
 
PARAMETER {
        v (mV)
        celsius (degC)
        dt (ms)
        gksoma = .0319 (mho/cm2)
        ek (mV)
}
 
STATE {
        n 
}
 
ASSIGNED {
        ik (mA/cm2)
        ninf 
	nexp 
	ntau (ms)
}
 
INITIAL {
	rate(v)
	n = ninf
}

BREAKPOINT {
        SOLVE state METHOD cnexp
	ik = gksoma*n*n*n*n*(v - ek)    
}

DERIVATIVE state {	:exact when v held constant
	rate(v)
	n' = (ninf-n)/ntau
}
UNITSOFF
PROCEDURE rate(v(mV)) {  :Computes rate and other constants at 
		      :current v.
                      :Call once from HOC to initialize inf at resting v.
        LOCAL q10, tinc, alpha, beta
        TABLE ninf, nexp, ntau DEPEND celsius FROM -200 TO 
100 WITH 300
		q10 = 3^((celsius - 24)/10)
		tinc = -dt*q10
		alpha = 0.018*vtrap(-(v-25),25)
		beta = 0.0036*vtrap(v-35,12)
		ntau = 1/(alpha + beta)
		ninf = alpha*ntau
		nexp = 1-exp(tinc/ntau)
}
FUNCTION vtrap(x,y) {	:Traps for 0 in denominator of rate eqns.
		if (fabs(x/y) < 1e-6) {
			vtrap = y*(1 - x/y/2)
		}else{
			vtrap = x/(exp(x/y) - 1)
		}
}
UNITSON

Loading data, please wait...