Molecular layer interneurons in cerebellum encode valence in associative learning (Ma et al 2020)

 Download zip file 
Help downloading and running models
Accession:266578
We used two-photon microscopy to study the role of ensembles of cerebellar molecular layer interneurons (MLIs) in a go-no go task where mice obtain a sugar water reward. In order to begin understanding the circuit basis of our findings in changes in lick behavior with chemogenetics in the go-no go associative learning olfactory discrimination task we generated a simple computational model of MLI interaction with PCs.
Reference:
1 . Ma M, Futia GL, De Souza FM, Ozbay BN, Llano I, Gibson EA, Restrepo D (2020) Molecular layer interneurons in the cerebellum encode for valence in associative learning Nature Communications, accepted
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: Cerebellum; Mouse;
Cell Type(s): Cerebellum Purkinje GABA cell; Cerebellum interneuron stellate GABA cell;
Channel(s):
Gap Junctions:
Receptor(s): AMPA; GabaA;
Gene(s):
Transmitter(s): Gaba; Glutamate;
Simulation Environment: NEURON;
Model Concept(s): Action Potentials; Detailed Neuronal Models;
Implementer(s): Simoes-de-Souza, Fabio [fabio.souza at ufabc.edu.br];
Search NeuronDB for information about:  Cerebellum Purkinje GABA cell; Cerebellum interneuron stellate GABA cell; GabaA; AMPA; Gaba; Glutamate;
/
MaEtAl2020
README.html
bkpkj.mod *
cad.mod *
cadiff.mod *
cae.mod *
cap2.mod *
captain.mod *
cat.mod *
cha.mod *
erg.mod *
gkca.mod *
Golgi_Ca_LVA.mod *
Golgi_KA.mod *
Golgi_KV.mod
Golgi_Na.mod *
hpkj.mod *
k23.mod *
ka.mod *
kc3.mod *
kd.mod *
kdyn.mod *
khh.mod *
km.mod *
kpkj.mod *
kpkj2.mod *
kpkjslow.mod *
kv1.mod *
leak.mod *
lkpkj.mod *
myexchanger.mod *
myexchangersoma.mod *
mypump.mod *
mypumpsoma.mod *
nadifl.mod *
narsg.mod *
newnew.mod *
pump.mod *
10480.tmp
2_compartment_template.hoc
full.ses *
lesbos.ses *
mosinit.hoc
mosinit_PC_SC_SminusCNO.hoc
mosinit_PC_SC_SminusSaline.hoc
mosinit_PC_SC_SplusCNO.hoc
mosinit_PC_SC_SplusSaline.hoc
PC_alx.swc
PF_template.hoc
Plot_results.m
SC_morphology.hoc
SC_template.hoc
SC_withoutaxon.swc
screenshot.png
                            
: FORREST MD (2014) Two Compartment Model of the Cerebellar Purkinje Neuron

COMMENT

Ih current	 - hyperpolarization-activated nonspecific Na and K channel
		 - contributes to the resting membrane potential
		 - controls the afterhyperpolarization
Reference:

1.	Maccaferri, G. and McBain, C.J. The hyperpolarization-activated current
	(Ih) and its contribution to pacemaker activity in rat CA1 hippocampal
	stratum oriens-alveus interneurons, J. Physiol. 497.1:119-130,
	1996.

		V1/2 = -84.1 mV
		   k = 10.2
		reversal potential = -32.9 +/- 1.1 mV

at -70 mV, currents were fitted by a single exponetial of t = 2.8+/- 0.76 s
at -120 mV, two exponentials were required, t1 = 186.3+/-33.6 ms 
t2 = 1.04+/-0.16 s


2.	Maccaferri, G. et al. Properties of the
	Hyperpoarization-activated current in rat hippocampal CA1 Pyramidal
	cells. J. Neurophysiol. Vol. 69 No. 6:2129-2136, 1993.

		V1/2 = -97.9 mV
		   k = 13.4
		reversal potential = -18.3 mV

3.	Pape, H.C.  Queer current and pacemaker: The
	hyperpolarization-activated cation current in neurons, Annu. Rev. 
	Physiol. 58:299-327, 1996.

		single channel conductance is around 1 pS
		average channel density is below 0.5 um-2
		0.5 pS/um2 = 0.00005 mho/cm2 = 0.05 umho/cm2		
4.	Magee, J.C. Dendritic Hyperpolarization-Activated Currents Modify
	the Integrative Properties of Hippocampal CA1 Pyramidal Neurons, J.
	Neurosci., 18(19):7613-7624, 1998

Deals with Ih in CA1 pyramidal cells.  Finds that conductance density
increases with distance from the soma.

soma g = 0.0013846 mho/cm2
dendrite g (300-350 um away) = 0.0125 mho/cm2
see Table 1 in th paper

ENDCOMMENT

 UNITS {
        (mA) = (milliamp)
        (mV) = (millivolt)
}
 
NEURON {
        SUFFIX cha
        USEION h READ eh WRITE ih VALENCE 1
        RANGE gkhbar,ih
        GLOBAL rinf, rexp, tau_r
}
 
INDEPENDENT {t FROM 0 TO 1 WITH 1 (ms)}
 
PARAMETER {
        v (mV)
        p = 5 (degC)
        dt (ms)
        gkhbar = 0.001385       : 0.001385 (mho/cm2)			
        eh = -32.9 (mV)
}
 
STATE {
        r
}
 
ASSIGNED {
        ih (mA/cm2)
	rinf rexp
	tau_r
}
 
BREAKPOINT {
        SOLVE deriv METHOD derivimplicit
        ih = gkhbar*r*(v - eh)
}
 
INITIAL {
	rates(v)
	r = rinf
}

DERIVATIVE deriv { :Computes state variable h at current v and dt.
	rates(v)
	r' = (rinf - r)/tau_r
}

PROCEDURE rates(v) {  :Computes rate and other constants at current v.
                      :Call once from HOC to initialize inf at resting v.
        TABLE rinf, rexp, tau_r DEPEND dt, p FROM -200
TO 100 WITH 300
	rinf = 1/(1 + exp((v+84.1)/10.2))
	rexp = 1 - exp(-dt/(tau_r))
	tau_r = 100 + 1/(exp(-17.9-0.116*v)+exp(-1.84+0.09*v))
}
 
UNITSON



Loading data, please wait...