Cerebellar granular layer (Maex and De Schutter 1998)

 Download zip file 
Help downloading and running models
Circuit model of the granular layer representing a one-dimensional array of single-compartmental granule cells (grcs) and Golgi cells (Gocs). This paper examines the effects of feedback inhibition (grc -> Goc -> grc) versus feedforward inhibition (mossy fibre -> Goc -> grc) on synchronization and oscillatory behaviour.
1 . Maex R, De Schutter E (1998) Synchronization of golgi and granule cell firing in a detailed network model of the cerebellar granule cell layer. J Neurophysiol 80:2521-37 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism: Cerebellum;
Cell Type(s): Cerebellum interneuron granule GLU cell; Cerebellum golgi cell;
Channel(s): I Na,t; I A; I h; I K,Ca; I L high threshold; I_KD;
Gap Junctions:
Receptor(s): AMPA; GabaA; NMDA;
Transmitter(s): Gaba; Glutamate;
Simulation Environment: GENESIS;
Model Concept(s): Synchronization; Oscillations;
Implementer(s): Maex, Reinoud [reinoud at bbf.uia.ac.be];
Search NeuronDB for information about:  Cerebellum interneuron granule GLU cell; GabaA; AMPA; NMDA; I Na,t; I L high threshold; I A; I h; I K,Ca; I_KD; Gaba; Glutamate;
// genesis

//float EK=-0.085
float Gbar = 1.0

float offset = 0.01

int tab_ydivs = {tab_xdivs}
float tab_ymin = {CCaI}
float tab_ymax = {Ca_tab_max}
float Temp = 22
float temperature = 5

/* non-inactivating BK-type Ca-dependent K current
** Moczydlowski and Latorre 1983, J Gen Physiol 82, 511-542.
** Uses original parameters. 
** Includes Paul Smolen's bug fix. 
/* scaled for units: V, sec, mM */
function make_Moczyd_KC
	int i, j
	float ginf, itau, c, dc, cmin, cmax
	float x, dx, y, dy
	float a, b
	float ZFbyRT = 23210/(273.15 + (Temp))

	if (!({exists Moczyd_KC}))
		create tab2Dchannel Moczyd_KC
		setfield Moczyd_KC Ek {EK} Gbar {Gbar}  \
		    Xindex {VOLT_C1_INDEX} Xpower 1 Ypower 0 Zpower 0

		call Moczyd_KC TABCREATE X {tab_xdivs} {tab_xmin}  \
		    {tab_xmax} {tab_ydivs} {tab_ymin} {tab_ymax}
echo diag Moczyd_KC 1
	dx = ({tab_xmax} - {tab_xmin})/{tab_xdivs}
	dy = (tab_ymax - tab_ymin)/tab_ydivs
	x = {tab_xmin} - {offset}
	for (i = 0; i <= ({tab_xdivs}); i = i + 1)
		y = tab_ymin
		for (j = 0; j <= (tab_ydivs); j = j + 1)
			/* \
			     Must check that the following are scaled correctly!!  \
			a = 2.5e3/(1 + ((1.5e-3*{exp {-0.085e3*x}})/y))
			b = 1.5e3/(1 + (y/(150e-6*{exp {-0.077e3*x}})))
			itau = a + b
			ginf = a/itau
			setfield Moczyd_KC X_A->table[{i}][{j}] {temperature * a}
			setfield Moczyd_KC X_B->table[{i}][{j}] {temperature * itau}
			y = y + dy
		x = x + dx
    setfield Moczyd_KC X_A->calc_mode {LIN_INTERP}
    setfield Moczyd_KC X_B->calc_mode {LIN_INTERP}
echo diag Moczyd_KC 3
create neutral comp
setfield comp x -0.10
setfield comp y 0.10
addmsg comp Moczyd_KC VOLTAGE x
addmsg comp Moczyd_KC CONCEN y
echo diag Moczyd_KC 4

Loading data, please wait...