Cerebellar cortex oscil. robustness from Golgi cell gap jncs (Simoes de Souza and De Schutter 2011)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:139656
" ... Previous one-dimensional network modeling of the cerebellar granular layer has been successfully linked with a range of cerebellar cortex oscillations observed in vivo. However, the recent discovery of gap junctions between Golgi cells (GoCs), which may cause oscillations by themselves, has raised the question of how gap-junction coupling affects GoC and granular-layer oscillations. To investigate this question, we developed a novel two-dimensional computational model of the GoC-granule cell (GC) circuit with and without gap junctions between GoCs. ..."
Reference:
1 . Simões de Souza F, De Schutter E (2011) Robustness effect of gap junctions between Golgi cells on cerebellar cortex oscillations Neural Systems & Circuits 1:7:1-19
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 cell; Cerebellum golgi cell;
Channel(s):
Gap Junctions: Gap junctions;
Receptor(s): GabaA; AMPA; NMDA;
Gene(s): HCN1; HCN2;
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Oscillations; Synchronization; Action Potentials;
Implementer(s): Simoes-de-Souza, Fabio [fabio.souza at ufabc.edu.br];
Search NeuronDB for information about:  Cerebellum interneuron granule cell; GabaA; AMPA; NMDA;
/
network
data
README.txt
gap.mod
Golgi_BK.mod *
Golgi_Ca_HVA.mod *
Golgi_Ca_LVA.mod *
Golgi_CALC.mod *
Golgi_CALC_ca2.mod *
Golgi_hcn1.mod *
Golgi_hcn2.mod *
Golgi_KA.mod *
Golgi_KM.mod *
Golgi_KV.mod *
Golgi_lkg.mod *
Golgi_Na.mod *
Golgi_NaP.mod *
Golgi_NaR.mod *
Golgi_SK2.mod *
GRC_CA.mod *
GRC_CALC.mod *
GRC_KA.mod *
GRC_KCA.mod *
GRC_KIR.mod *
GRC_KM.mod *
GRC_KV.mod *
GRC_LKG1.mod *
GRC_LKG2.mod *
GRC_NA.mod *
K_conc.mod *
Na_conc.mod *
Golgi_ComPanel.hoc *
Golgi_template.hoc
granule_template.hoc
MF_template.hoc
mosinit.hoc
network.hoc
utils.hoc *
                            
// Cerebellum Granular Cell Model

xopen("$(NEURONHOME)/lib/hoc/noload.hoc") 


//DEFINITION OF CELL TEMPLATE

begintemplate grc

public soma
public x,y,z

//object variable will refer to a list that will hold an arbitrary number of NetCon objects.

public nclist
objectvar nclist  

//synapses         
public synAMPA,synNMDA,synGABAa1,synGABAa2 
objectvar synAMPA,synNMDA,synGABAa1,synGABAa2 

//counting spikes
public spiketimes,spikecount
objref spiketimes,spikecount

public voltagem
objref voltagem


tmax=1000   

create soma	//create compartment

public setv

proc setv() {
forall v=$1
}

proc init() {

x=$1
y=$2
z=$3

spiketimes=new Vector()
lastspikecount=0

voltagem=new Vector()

soma {

 pt3dclear()
 pt3dadd(x,y,z,10)	//set position of cell
 pt3dadd(x,y,z+10,10)
 
nseg=1
diam=11.8
L=11.8 
Ra=100
cm=1
celsius = 30

Area = PI*soma.diam*soma.L*1e-8

insert GRC_LKG1 
insert GRC_LKG2
insert GRC_NA 
insert GRC_KV
insert GRC_KA
insert GRC_KIR 
insert GRC_KCA 
insert GRC_CALC
//insert GrC_CaHVA
insert GRC_CA
insert GRC_KM 

gl_GRC_LKG1=5.68e-5	  //2009 model
ggaba_GRC_LKG2 = 3e-05 	  //2009 model
gkbar_GRC_KV= 0.003  	 //2009 model
gkbar_GRC_KA = 0.0032	  //2009 model
gkbar_GRC_KIR = 0.0009  //2009 model
gkbar_GRC_KM= 0.00025   //0.0001  //2009 model
gnabar_GRC_NA= 0.013     //2009 model  
gcabar_GRC_CA= 0.00046    //2009 model
gkbar_GRC_KCA = 0.04   //ajusted Version 3

beta_GRC_CALC = 1.5


usetable_GRC_NA 	= 0			      			
usetable_GRC_KV 	= 0
usetable_GRC_KA 	= 0
usetable_GRC_KIR 	= 0	  
usetable_GRC_KCA 	= 0
usetable_GRC_CA 	= 0
usetable_GRC_KM 	= 0

ena = 87.39
ek = -84.69
//eca = 129.33
ecl = -65

nclist=new List()

//Synapses
synAMPA=new Exp2Syn(0.5) //Adding synapse with rise and decay time
synAMPA.tau1=0.03  // ms rise time 
synAMPA.tau2=0.5  // ms decay time 
synAMPA.e=0     // mV reversal potential 

synNMDA=new Exp2Syn(0.5) //Adding synapse with rise and decay time
synNMDA.tau1=1  // ms rise time
synNMDA.tau2=13.3  // ms decay time 
synNMDA.e=0     // mV reversal potential 

synGABAa1=new Exp2Syn(0.5) //Adding synapse with rise and decay time
synGABAa1.tau1=0.31  // ms rise time 
synGABAa1.tau2=8.8  // ms decay time 
synGABAa1.e=-75     // mV reversal potential 

synGABAa2=new Exp2Syn(0.5) //Adding synapse with rise and decay time
synGABAa2.tau1=0.31  // ms rise time 
synGABAa2.tau2=300  // ms decay time 
synGABAa2.e=-75     // mV reversal potential 

//counting spikes
spikecount=new APCount(0.5)
spikecount.thresh=-20
spikecount.record(spiketimes)

//Saving Vm
voltagem.record(&v(0.5))
}
}

endtemplate grc


/*
 General settings
*/ 
v_init = -70

tstop = 1000 //for IClamp
dt = 0.025

secondorder=0  //method implicit backward euler

celsius = 30 
SogliaSpike=0



Loading data, please wait...