5-neuron-model of neocortex for producing realistic extracellular AP shapes (Van Dijck et al. 2012)

 Download zip file 
Help downloading and running models
Accession:226812
This is a 5-neuron model of neocortex, containing one tufted layer-5 pyramidal cell, two non-tufted pyramidal cells, and two inhibitory interneurons. It was used to reproduce extracellular spike shapes in a study comparing algorithms for spike sorting and electrode selection. The neuron models are adapted from Dyhrfjeld-Johnsen et al. (2005).
Reference:
1 . Van Dijck G, Seidl K, Paul O, Ruther P, Van Hulle MM, Maex R (2012) Enhancing the yield of high-density electrode arrays through automated electrode selection. Int J Neural Syst 22:1-19 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Extracellular; Neuron or other electrically excitable cell; Realistic Network;
Brain Region(s)/Organism:
Cell Type(s): Neocortex U1 L5B pyramidal pyramidal tract GLU cell; Neocortex U1 L2/6 pyramidal intratelencephalic GLU cell;
Channel(s):
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: GENESIS;
Model Concept(s):
Implementer(s): Maex, Reinoud [reinoud at bbf.uia.ac.be];
Search NeuronDB for information about:  Neocortex U1 L5B pyramidal pyramidal tract GLU cell; Neocortex U1 L2/6 pyramidal intratelencephalic GLU cell;
/
Five-neuron-neocortex
L5P37C-onlybasal
channels
README
Axon_chans.g *
Axon_chans_tab.g
Axon_comps.g
DiffRm.g *
DS1_141099_rot2_sc_defmesh_axon_onlybasal.p
electrodes_fixbug.g
electrodes_try.g *
Excitatory_fibres.g
Fibres.g *
Firing_rate_modulation.g *
Firing_rate_profile.g *
Golgi_grc_multicomp_ax.g
Gran_chan_tab.g *
Gran_comp_soma_dend_axon.g
Gran_const.g *
Gran_synchan.g *
Harsch-Robinson_modulation.g *
Hgradient.g *
Inhibitory_fibres.g *
L5P_ascout.g
L5P_ascout_exp.g *
L5P_chans.g *
L5P_chans_tab.g
L5P_chans_tab_Temp.g *
L5P_chans_Temp.g *
L5P_comps.g *
L5P_comps+axon+syn.g
L5P_const.g *
L5P_const+axon+syn.g *
L5P_graph.g
L5P_history.g
L5P_onlybasal_make.g
L5P_synchan.g
L5P37C_onlybasal.g
nsynapses.g *
                            
// genesis

// Creates for each compartment element of L5P an associated randomspike element
// located at the same position, and connected to the L5P cell using volumeconnect.
// The randomspike elements are collected in /Excitatory_fibres.
// The connections may be laid horizontally, using only the z-coordinate as a criterion,
// or radially, allowing the fibre to connect to all compartments at about the same
// position from the soma.

// Inaddition (from 6Oct07) creates for each randomspike element a diffamp from which 
// it will receive a RATE msg.

// Moreover, both randomspike and diffamp populations are divided into 3 equal parts:
// feedforward afferents, intracolumn feedbak afferents and intercolumn feedback afferents.

 
int  index


include Fibres.g

make_fibres /Inhibitory_fibres


 
// Connecting fibres to L5P *****

ce /Inhibitory_fibres


// First connection scheme: each fiber runs horizontally and is allowed to make a synapse
// on each compartment in (approximately) the same horizontal plane


echo Connecting inhibitory feedforward fibres 

      volumeconnect FF/FF_fibre[] \
                    /L5P/p#[][TYPE=compartment]/GABA \
                    -relative \
                    -sourcemask box -1 -1 -1 1 1 1 \
                    -destmask ellipsoid 0 0 0 1 1 2e-6 \ 
                    -probability 0.025 // 0.05 // 1.0 // {P_parallel_fiber_to_Golgi_cell_synapse}

      volumeweight FF/FF_fibre[] \
                   /L5P/p#[][TYPE=compartment]/GABA  -fixed \
                   {weight_GABA_synapse} -uniform {weight_distribution}

      volumedelay   FF/FF_fibre[] \
                    /L5P/p#[][TYPE=compartment]/GABA \
                    -fixed 0.001 // 08 // 0.001 // 0.002
// -radial {E_fibre_conduction_velocity} -uniform {delay_distribution}


echo Connecting inhibitory intracolumnar feedback fibres 

      volumeconnect FBintra/FBintra_fibre[] \
                    /L5P/p#[][TYPE=compartment]/GABA \
                    -relative \
                    -sourcemask box -1 -1 -1 1 1 1 \
                    -destmask ellipsoid 0 0 0 1 1 2e-6 \ 
                    -probability 0.05 // 1.0 // {P_parallel_fiber_to_Golgi_cell_synapse}

      volumeweight FBintra/FBintra_fibre[] \
                   /L5P/p#[][TYPE=compartment]/GABA  -fixed \
                   {weight_GABA_synapse} -uniform {weight_distribution}

      volumedelay   FBintra/FBintra_fibre[] \
                    /L5P/p#[][TYPE=compartment]/GABA \
                    -fixed 0.02
//                    -radial {E_fibre_conduction_velocity} -uniform {delay_distribution}

/*
echo Connecting inhibitory intercolumnar feedback fibres

      volumeconnect FBinter/FBinter_fibre[] \
                    /L5P/p#[][TYPE=compartment]/GABA \
                    -relative \
                    -sourcemask box -1 -1 -1 1 1 1 \
                    -destmask ellipsoid 0 0 0 1 1 2e-6 \ 
                    -probability 0.05 // 1.0 // {P_parallel_fiber_to_Golgi_cell_synapse}

      volumeweight FBinter/FBinter_fibre[] \
                   /L5P/p#[][TYPE=compartment]/GABA  -fixed \
                   {weight_GABA_synapse} -uniform {weight_distribution}

      volumedelay   FBinter/FBinter_fibre[] \
                    /L5P/p#[][TYPE=compartment]/GABA \
                    -radial {E_fibre_conduction_velocity} -uniform {delay_distribution}
*/

/*

// Second connection scheme: each fiber is allowed to make a synapse on compartments
// located from the soma (the origin) at the same distance as the compartment to which the fibre
// has been associated (remember that a fibre for each compartment was created).
int i
float r, x1, y1, z1

      ce {FF}

      for (i = 0; i < 686; i = {i} + 1)
          x1 = {getfield FF_fibre[{i}] x}
          y1 = {getfield FF_fibre[{i}] y}
          z1 = {getfield FF_fibre[{i}] z}
          r = {sqrt {x1*x1 + y1*y1 + z1*z1}}
//          echo {r}

          volumeconnect FF_fibre[{i}] \
                    /L5P/p#[][TYPE=compartment]/GABA \
                    -sourcemask box -1 -1 -1 1 1 1 \
                    -destmask ellipsoid 0 0 0 {r + 1e-6} {r + 1e-6} {r + 1e-6} \ 
                    -desthole ellipsoid 0 0 0 {r - 1e-6} {r - 1e-6} {r - 1e-6} \ 
                    -probability 0.1 // 0.05
      end

      volumeweight FF_fibre[] \
                   /L5P/p#[][TYPE=compartment]/GABA  -fixed \
                   {weight_AMPA_synapse} -uniform {weight_distribution}

      volumedelay   FF_fibre[] \
                    /L5P/p#[][TYPE=compartment]/GABA \
                    -radial {E_fibre_conduction_velocity} -uniform {delay_distribution}


      ce {FBintra}

      for (i = 0; i < 686; i = {i} + 1)
          x1 = {getfield FBintra_fibre[{i}] x}
          y1 = {getfield FBintra_fibre[{i}] y}
          z1 = {getfield FBintra_fibre[{i}] z}
          r = {sqrt {x1*x1 + y1*y1 + z1*z1}}
//          echo {r}
  
          volumeconnect FBintra_fibre[{i}] \
                    /L5P/p#[][TYPE=compartment]/GABA \
                    -sourcemask box -1 -1 -1 1 1 1 \
                    -destmask ellipsoid 0 0 0 {r + 1e-6} {r + 1e-6} {r + 1e-6} \ 
                    -desthole ellipsoid 0 0 0 {r - 1e-6} {r - 1e-6} {r - 1e-6} \ 
                    -probability 0.1 // 0.05
      end

      volumeweight FBintra_fibre[] \
                   /L5P/p#[][TYPE=compartment]/GABA  -fixed \
                   {weight_AMPA_synapse} -uniform {weight_distribution}

      volumedelay   FBintra_fibre[] \
                    /L5P/p#[][TYPE=compartment]/GABA \
                    -radial {E_fibre_conduction_velocity} -uniform {delay_distribution}




      ce {FBinter}

      for (i = 0; i < 686; i = {i} + 1)
          x1 = {getfield FBinter_fibre[{i}] x}
          y1 = {getfield FBinter_fibre[{i}] y}
          z1 = {getfield FBinter_fibre[{i}] z}
          r = {sqrt {x1*x1 + y1*y1 + z1*z1}}
//          echo {r}
  
          volumeconnect FBinter_fibre[{i}] \
                    /L5P/p#[][TYPE=compartment]/GABA \
                    -sourcemask box -1 -1 -1 1 1 1 \
                    -destmask ellipsoid 0 0 0 {r + 1e-6} {r + 1e-6} {r + 1e-6} \ 
                    -desthole ellipsoid 0 0 0 {r - 1e-6} {r - 1e-6} {r - 1e-6} \ 
                    -probability 0.1 // 0.05
      end

      volumeweight FBinter_fibre[] \
                   /L5P/p#[][TYPE=compartment]/GABA  -fixed \
                   {weight_AMPA_synapse} -uniform {weight_distribution}

      volumedelay   FBinter_fibre[] \
                    /L5P/p#[][TYPE=compartment]/GABA \
                    -radial {E_fibre_conduction_velocity} -uniform {delay_distribution}




*/