Pyramidal neuron, fast, regular, and irregular spiking interneurons (Konstantoudaki et al 2014)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:168310
This is a model network of prefrontal cortical microcircuit based primarily on rodent data. It includes 16 pyramidal model neurons, 2 fast spiking interneuron models, 1 regular spiking interneuron model and 1 irregular spiking interneuron model. The goal of the paper was to use this model network to determine the role of specific interneuron subtypes in persistent activity
Reference:
1 . Konstantoudaki X, Papoutsi A, Chalkiadaki K, Poirazi P, Sidiropoulou K (2014) Modulatory effects of inhibition on persistent activity in a cortical microcircuit model. Front Neural Circuits 8:7 [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:
Cell Type(s): Neocortex fast spiking (FS) interneuron; Neocortex spiking regular (RS) neuron; Neocortex spiking low threshold (LTS) neuron; Neocortex spiking irregular interneuron;
Channel(s): I Na,p; I Na,t; I L high threshold; I T low threshold; I A; I K; I h; I_Ks; I_KD;
Gap Junctions:
Receptor(s): GabaA; GabaB; AMPA; NMDA;
Gene(s):
Transmitter(s): Gaba; Glutamate;
Simulation Environment: NEURON;
Model Concept(s): Activity Patterns; Synchronization; Active Dendrites;
Implementer(s): Sidiropoulou, Kyriaki [sidirop at imbb.forth.gr]; Konstantoudaki, Xanthippi [xeniakons at gmail.com];
Search NeuronDB for information about:  GabaA; GabaB; AMPA; NMDA; I Na,p; I Na,t; I L high threshold; I T low threshold; I A; I K; I h; I_Ks; I_KD; Gaba; Glutamate;
/
KonstantoudakiEtAl2014
experiment
data
ampa.mod *
ampain.mod *
cadyn.mod *
cadynin.mod *
cal.mod *
calc.mod *
calcb.mod *
can.mod *
cancr.mod *
canin.mod *
car.mod *
cat.mod *
catcb.mod *
gabaa.mod *
gabaain.mod *
gabab.mod *
h.mod *
hcb.mod *
hin.mod *
ican.mod *
iccb.mod *
iccr.mod *
icin.mod *
iks.mod *
ikscb.mod *
ikscr.mod *
iksin.mod *
kadist.mod *
kadistcr.mod *
kadistin.mod *
kaprox.mod *
kaproxcb.mod *
kaproxin.mod *
kca.mod *
kcain.mod *
kct.mod *
kctin.mod *
kdr.mod *
kdrcb.mod *
kdrcr.mod *
kdrin.mod *
naf.mod *
nafcb.mod *
nafcr.mod *
nafin.mod *
nafx.mod *
nap.mod *
netstim.mod *
NMDA.mod *
NMDAIN.mod *
sinclamp.mod *
cb.hoc
cr.hoc
ExperimentControl.hoc *
final.hoc
incell.hoc
net.hoc
pfc_pc_temp.hoc
run
run_orig
                            
//--------------Set parameters and run this for the net (16 Pcells and 4 Interneurons--> 2PV, 1CB, 1CR)
//--------------Written 8-5-08, modified 08-10-08


load_file("nrngui.hoc")			//load main NEURON library
cvode.active(0)
xopen ("net.hoc")                    //load net structure
strdef syscmd
strdef data_dir
data_dir = "data"
sprint(syscmd, "mkdir -p %s", data_dir)  // make the data directory
system(syscmd)
printf("Opened_pcell")

//--------------------------------------------Set parameters

tstop=5000
steps_per_ms=10
dt=0.1


//-------------------Set synaptic weights


ampaweightpr=0.00008	//Initial stimulus 
ampaweight=0.000065	//according to Wang 2008 
ampaweightin=0.000085	
gabaweight=0.00083	 
gabaweightb=gabaweight*0.35	//ratio gabab to gabaa 0.2	


autogabaweight=0.0073*0.75	//According to Bacci (2003)

nmdaweight=ampaweight*21.5	//according to wang et al, 2008

nmdaweightin=ampaweightin*0.52  //according to wang and gao, Neuropharmacology

ampaweightcb=0.000029			
nmdaweightcb=ampaweightcb*0.86		

ampaweightcr=0.000046			//20pA sta -72 --> 28/05/2010	
nmdaweightcr=ampaweightcr*2.2		

gabaweightcrcb=0.0083*0.25              	//cr-cb connections 

gabaweightcb=0.0006*1.5                  	//cb-pc connections 

gabaweightcr=0.00087*3                 	//cr-pc connections 





//-------------------Set # of synapses

inmaxsyn=120              		//initial stimulation
maxsyn=24	                   	//PC-PC connections
automaxsyn=8                		//autosynapses 
maxsyn1=1                    		//IN-IN connections
maxsyn2=12                		//PC-IN connections 
maxsyn3=15                  		 //IN-PC soma connections
maxsyn4=15				//IN-PC dend1 connections 
maxsyn5=14				//PC-CB connections 
maxsyn6=7				//PC-CR connections 
maxsyn7=2				//CR-CB connections 
maxsyn8=12				//CB-PC connections
maxsyn9=10				//CR-PC connections 

///------------------------------------Pharmacological procedures
proc sadp_soma() {
forsec soma_list{
for(x) {
	fi2=0	//factor for ican
	if(ismembrane("ican")) for(x) { gbar_ican(x)= gbar_ican(x)*fi2 } 
	}}
}
proc sadp_dend() {
forsec dend_list{
for(x) {
	fi2=0	//factor for ican
	if(ismembrane("ican")) for(x) { gbar_ican(x)= gbar_ican(x)*fi2 } 
	}}
}


//--------------------------------------Graphs
xopen("../bash_templates/basic-graphics.hoc") 	
addgraph_2("Pcells[0].soma.v(0.5)",0,tstop, -70, 50)
addgraph_2("INcells[0].soma.v(0.5)",0,tstop, -70, 50)
addgraph_2("CBcells[0].soma.v(0.5)",0,tstop, -70, 50)
addgraph_2("CRcells[0].soma.v(0.5)",0,tstop, -70, 50)



//-------------------------------------------------------Multiple Runs----------------------------------//
strdef syscmd, data_dir, running, tmpstr
//-----Variables
n=int(tstop/dt)
//-----Objects for record data
objref cv
cv=new CVode(0)
objref PCv[nPcells], PCt[nPcells], INv[nINcells], INt[nINcells], CBv[nCBcells], CBt[nCBcells], CRv[nCRcells], CRt[nCRcells]
strdef temp
objref vsoma1, vsoma2, vsoma3, vsoma4, vsoma5, vsoma6, vsoma7, vsoma8, vsoma9, vsoma10, vsoma11, vsoma12, vsoma13, vsoma14, vsoma15, vsoma16, insoma1, insoma2, cbsoma, crsoma, curampa, curnmda, curgabaa, curgabab
objref campa[maxsyn+automaxsyn], campat[maxsyn+automaxsyn], cnmda[maxsyn+automaxsyn], cnmdat[maxsyn+automaxsyn], cgabaa[maxsyn3], cgabaat[maxsyn3], cgabab[maxsyn3], cgababt[maxsyn3], ampac, nmdac, gabaac, gababc	


// to reinstate for loop: comment the runs=0 and the xpanel and
// uncomment the for statement below, the run command and the close
// curly paranethesis at the bottom
//
// for runs = 0, 0 {
runs=0 // comment this line out to reinstate for loop
//--------------------------------------Call procedures
stimulation(runs)
pyramidals(runs)
interneurons(runs)
pyrin(runs)
inpyr(runs)
//inpyrs(runs) DO NOT ACTIVATE!!!
pyrcb(runs)
pyrcr(runs)
crcb(runs)
cbpyr(runs)
crpyr(runs)
noise(runs)
sadp_soma()
print runs

//-------------------------------------Record data

// comment xpanel and uncomment run command below to reinstate for loop:
xpanel("Konstantoudaki et al. 2014")
    xbutton("run","run()")
xpanel()
// run() // remove the comments from these lines to reinstate for loop
// 
//}

Loading data, please wait...