Small world networks of Type I and Type II Excitable Neurons (Bogaard et al. 2009)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:121259
Implemented with NEURON 5.9, four model neurons with varying excitability properties affect the spatiotemporal patterning of small world networks of homogeneous and heterogeneous cell population.
Reference:
1 . Bogaard A, Parent J, Zochowski M, Booth V (2009) Interaction of cellular and network mechanisms in spatiotemporal pattern formation in neuronal networks. J Neurosci 29:1677-87 [PubMed]
Citations  Citation Browser
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: Hippocampus;
Cell Type(s):
Channel(s): I Na,t; I A; I K; I h;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Synchronization; Spatio-temporal Activity Patterns; Epilepsy;
Implementer(s):
Search NeuronDB for information about:  I Na,t; I A; I K; I h;
/
bogaard2009
readme.html
h.mod *
h2.mod
kaprox.mod *
kaprox2.mod
kdrca1.mod *
kdrca12.mod
nax.mod
nax2.mod
WMPas.mod
WMPas2.mod
WMPasDend.mod
WMPasDend2.mod
hetrig.hoc
hetrun.hoc
homrig.hoc
homrun.hoc
modelspecs.hoc
mosinit.hoc
screenshothet.jpg
screenshothom.jpg
                            
begintemplate cellA //Cell A is Type I
  public soma, dendrite, nclist
  create soma, dendrite
  objectvar nclist
  
  proc init() {
  
  nclist = new List()
  // Migliore parameters
  gna =  0.3
  gkdr = 0.1
  gka = 0.048
  ghd=0.0001
  vhalfn=13
  celsius = 34.0  
  Vrest = -65
  
  
  // Borg-Graham morphology
  
  connect dendrite(0), soma(0) 	
  
  soma{ 
  /* The soma in the original Working Model is originally a 
  sphere with diam = 35 um. In NEURON it is a cylinder.  
  Also, in the simulator SURF-HIPPO the underlying segment
  circuit is a single-ended scheme, while in NEURON it is a 
  center-tapped scheme. Therefore, the half-axial resistance 
  of the soma must be mached to the half-axial resistance of 
  the first dendritic segment.
  
  ==>  total soma surface = PI*diam*L = 3.85e3 um2       (1)
   
  ==>  half axial resistance of 1st dendritic segment:
          
                       0.5 * 4 * Ra * dendrite.L
           Rhalf = ---------------------------------------  = 2.1 MOhm.
                     PI * dendrite.diam^2 * dendrite.nseg
  
  
                       0.5 * 4 * Ra * soma.L
             ==>   -------------------------- = 2.1 MOhm  (2)
                          PI * soma.diam^2 
  
                   
           Now, substituting L from equation (1) gives 
                        L = 137.5 um
                        diam = 9.0 um
  
  */
  
  	nseg = 1
          L = 135.74
          diam = 9.0246
  	      cm = 0.7
  	insert WMPas			//Borg-Graham passive kinetics
  	insert hd  
          ghdbar_hd=ghd		//Migliore active conductances
  		    vhalfl_hd=-82
  		    ehd_hd=-30
  	insert nax  
          gbar_nax=gna 
  	insert kdr  
          gkdrbar_kdr=gkdr
  		    vhalfn_kdr=vhalfn
    insert kap  
          gkabar_kap = gka
  	ena=55
  	ek=-90
  }
  
  dendrite{
  	nseg = 5	// always put nseg first
  	cm = 0.7
  	L = 1200	// each segment is 240 um long
  	diam = 12
  	insert WMPasDend	//Borg-Graham passive kinetics
  }
  
  forall v=Vrest
  forall Ra = 200
  //forall psection()
  access soma
  }
endtemplate cellA

begintemplate cellB //Cell B is Type I
  public soma, dendrite, nclist
  
  create soma, dendrite
  
  objectvar nclist
  
  
  proc init() {
  
  nclist = new List()
  
  // Migliore parameters
  gna =  .15
  gkdr = 0.05
  KMULTP = 0.048
  ghd=0.0005
  celsius = 34.0  
  Vrest = -65
  
  // Borg-Graham morphology
  connect dendrite(0), soma(0) 	
  
  soma{ 
  
  	nseg = 1
        	L = 135.74
        	diam = 9.0246
  	cm = 0.7
  	insert WMPas	//Borg-Graham passive kinetics
  	insert hd 
          ghdbar=ghd	//Migliore active conductances
  	       vhalfl_hd=-82
    insert nax  
          gbar_nax=.5     
    insert kdr 
          gkdrbar_kdr=.5
  	       vhalfn_kdr=13
    insert kap 
          gkabar_kap = KMULTP
  //      insert pas g_pas = 1/RmSoma Ra=RaSoma cm=CmSoma  //Migliore passive kinetics
  }
  
  dendrite{
  	nseg = 5	// always put nseg first
  	cm = 0.7
  	L = 1200	// each segment is 240 um long
  	diam = 12
  	insert WMPasDend	//Borg-Graham passive kinetics
  }
  
  
  forall Ra = 200
  //forall psection()
  
  access soma
  }
endtemplate cellB

begintemplate cellC //Cell C is Type II
  public soma, dendrite, nclist
  create soma, dendrite
  objectvar nclist
  
  proc init() {
  
  nclist = new List()
  // Migliore parameters
  gna =  0.3
  gkdr = 0.08
  gka = 0.048
  ghd=0.00075
  vhalfn=0
  celsius = 34.0  
  Vrest = -65
  
  
  // Borg-Graham morphology
  connect dendrite(0), soma(0) 	
  
  soma{ 
  	nseg = 1
          L = 135.74
          diam = 9.0246
  	cm = 0.7
  	insert WMPas			//Borg-Graham passive kinetics
  	insert hd  
          ghdbar_hd=ghd		//Migliore active conductances
  		    vhalfl_hd=-82
  		    ehd_hd=-30
  	insert nax  
          gbar_nax=gna 
  	insert kdr  
          gkdrbar_kdr=gkdr
  		    vhalfn_kdr=vhalfn
    insert kap
          gkabar_kap = gka
  	ena=55
  	ek=-90
  }
  
  dendrite{
  	nseg = 5	// always put nseg first
  	cm = 0.7
  	L = 1200	// each segment is 240 um long
  	diam = 12
  	insert WMPasDend	//Borg-Graham passive kinetics
  }
  
  forall v=Vrest
  forall Ra = 200
  //forall psection()
  
  access soma
  
  }
endtemplate cellC

begintemplate cellD //Cell D is Type II
  public soma, dendrite, nclist
  
  create soma, dendrite
  
  objectvar nclist
  
  
  proc init() {
  
  nclist = new List()
  
  // Migliore parameters
  gna =  .15
  gkdr = 0.7
  KMULTP = 0.03
  ghd=0.0005
  celsius = 34.0  
  Vrest = -65
  
  // Borg-Graham morphology
  connect dendrite(0), soma(0) 	
  
  soma{ 
  
  	nseg = 1
          L = 135.74
          diam = 9.0246
  	     cm = 0.7
  	insert WMPas2	//Borg-Graham passive kinetics
  	insert hd2
          ghdbar=ghd	//Migliore active conductances
  	      vhalfl_hd2=-82 
          ghdbar_hd2=.0005
    insert nax2  
          gbar_nax2=1.5     
    insert kdr2
          gkdrbar_kdr2=.7 
          vhalfn_kdr2=0
    insert kap2 
          gkabar_kap2 = .03
  }
  
  dendrite{
  	nseg = 5	// always put nseg first
  	cm = 0.7
  	L = 1200	// each segment is 240 um long
  	diam = 12
  	insert WMPasDend	//Borg-Graham passive kinetics
  }
  
  
  forall Ra = 200
  //forall psection()
  
  access soma
  }
endtemplate cellD