Dentate Gyrus Feed-forward inhibition (Ferrante et al. 2009)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:124291
In this paper, the model was used to show how that FFI can change a steeply sigmoidal input-output (I/O) curve into a double-sigmoid typical of buffer systems.
Reference:
1 . Ferrante M, Migliore M, Ascoli GA (2009) Feed-forward inhibition as a buffer of the neuronal input-output relation. Proc Natl Acad Sci U S A 106:18004-9 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network; Neuron or other electrically excitable cell;
Brain Region(s)/Organism: Dentate gyrus;
Cell Type(s): Dentate gyrus granule GLU cell; Dentate gyrus MOPP cell;
Channel(s): I Na,t; I A; I K;
Gap Junctions:
Receptor(s): GabaA; AMPA; NMDA;
Gene(s):
Transmitter(s): Gaba; Glutamate;
Simulation Environment: NEURON;
Model Concept(s): Action Potential Initiation; Activity Patterns; Ion Channel Kinetics; Synchronization; Spatio-temporal Activity Patterns; Action Potentials; Noise Sensitivity;
Implementer(s): Migliore, Michele [Michele.Migliore at Yale.edu]; Ferrante, Michele [mferr133 at bu.edu]; Ascoli, Giorgio A [ascoli at gmu.edu];
Search NeuronDB for information about:  Dentate gyrus granule GLU cell; GabaA; AMPA; NMDA; I Na,t; I A; I K; Gaba; Glutamate;
/
FFI
Granule_Fig_1B_right
fna.mod
ichan2.mod *
kdr.mod
nmdanetGR.mod
FIG_1B_right.hoc
fixnseg.hoc *
n500-axon.hoc
                            
TITLE fna.mod   
 
COMMENT
ENDCOMMENT
 
UNITS {
        (mA) = (milliamp)
        (mV) = (millivolt)
	(S) = (siemens)
}
 
NEURON {
        SUFFIX fna
        USEION na READ ena WRITE ina
        RANGE gnabar,  gna, alpha, beta, alphah, betah, sum
        GLOBAL minf, hinf, mtau, htau
}
 
PARAMETER {
        gnabar = .12 (S/cm2)	<0,1e9>
        gl = .0003 (S/cm2)	<0,1e9>
        :celsius (degC)
}
 
STATE {
        m h 
}
 
ASSIGNED {
        ena (mV)
        v (mV)
        :celsius (degC)

	alpha
	beta
	alphah
	betah
	sum
	gna (S/cm2)
        ina (mA/cm2)
        minf hinf 
	mtau (ms) htau (ms) 
}
 
LOCAL mexp, hexp
 
BREAKPOINT {
        SOLVE states METHOD cnexp
        gna = gnabar*m*m*m*h
	ina = gna*(v - ena)
}
 
 
INITIAL {
	rates(v)
	m = minf
	h = hinf
}

DERIVATIVE states {  
        rates(v)
        m' =  (minf-m)/mtau
        h' = (hinf-h)/htau
}
 
LOCAL q10

PROCEDURE rates(v(mV)) {  :Computes rate and other constants at current v.
                      :Call once from HOC to initialize inf at resting v.
        :LOCAL  alphah, betah
        :TABLE minf, mtau, hinf, htau,  DEPEND celsius FROM -100 TO 100 WITH 200

UNITSOFF
        q10=1
	:q10 = 3^((celsius - 6.3)/10)
                :"m" sodium activation system
        alpha = -.3 * vtrap((v+70-25),-5)
        beta =  .3 * vtrap((v+70-53),5)
        sum = alpha + beta
	mtau = 1/(q10*sum)
        minf = alpha/sum

                :"h" sodium inactivation system
        alphah = .23 /exp((v+70-3)/20)
        betah = 3.33 / (1+exp((v+70-55.5)/(-10)) )
        sum = alphah + betah
	htau = 1/(q10*sum)
        hinf = alphah/sum
       
}
 
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...