Fly lobular plate VS cell (Borst and Haag 1996, et al. 1997, et al. 1999)

 Download zip file 
Help downloading and running models
In a series of papers the authors conducted experiments to develop understanding and models of fly visual system HS, CS, and VS neurons. This model recreates the VS neurons from those papers with enough success to merit approval by Borst although some discrepancies remain (see readme).
1 . Borst A, Haag J (1996) The intrinsic electrophysiological characteristics of fly lobula plate tangential cells: I. Passive membrane properties. J Comput Neurosci 3:313-36 [PubMed]
2 . Haag J, Theunissen F, Borst A (1997) The intrinsic electrophysiological characteristics of fly lobula plate tangential cells: II. Active membrane properties. J Comput Neurosci 4:349-69 [PubMed]
3 . Haag J, Vermeulen A, Borst A (1999) The intrinsic electrophysiological characteristics of fly lobula plate tangential cells: III. Visual response properties. J Comput Neurosci 7:213-34 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell;
Brain Region(s)/Organism: Drosophila;
Cell Type(s): Fly lobular plate vertical system cell;
Channel(s): I Na,t; I K; I_K,Na;
Gap Junctions:
Simulation Environment: NEURON;
Model Concept(s): Reliability; Vision;
Implementer(s): Carnevale, Ted [Ted.Carnevale at]; Torben-Nielsen, Ben [btorbennielsen at];
Search NeuronDB for information about:  I Na,t; I K; I_K,Na;
Model that tries to mimic the VS model illustrated in the paper by Borst and Haag (1996,1997,1999)

Executing this model will produce a figure similar to fig 10 from Haag 1997.

Created by B. Torben-Nielsen (TENU @ OIST) and Ted Carnevale


/* model specification
for the morphology we took 72Vs4. (

// The model is created with the cell builder loaded through the session. However, for certainty we load the lastest model in the cell builder and re-assign some critical properties

CellBuild[0].continuous = 0

forall {
  // for the sections with vgated channels
  if (ismembrane("emdna")) {
    ena = 100
    ek = -20

// instrumentation, initialization, and run control

objref ccl

dend[504] {
  ccl = new IClamp(0.1666)


proc init() {
  // the ina-dependent k current requires finitialize() to be called twice
  // the first pass initializes ina,
  // and the second pass ensures that ina affects ik properly
  if (MAKE_VREST_UNIFORM == 1) {
    // assume that membrane is either completely passive
    // or that it has ina and ik in addition to pas
    forall if (ismembrane("emdna")) {
      e_pas = (ina + ik + g_pas*v_init)/g_pas
    } else {
      e_pas = v_init
    if ( {
    } else {
  if (MAKE_VREST_UNIFORM == 0) {
  	forall e_pas = v_init

objref g
g = Graph[1] // Graph[1] was created by the ses file loaded above

proc fig10e() { local i
  for i = -6, 6 {
    print "i: ", i, " (nA)"
    ccl.amp = i
    g.exec_menu("Keep Lines") // toggle on & off to keep this trace
    g.exec_menu("Keep Lines")
  ccl.amp = 0 // in case users want to create their own IClamp
    // and explore the model
  print "done"

xpanel("Resting potential")
xstatebutton("Force uniform resting potential",&MAKE_VREST_UNIFORM)
xbutton("Fig. 10e", "fig10e()")
xbutton("Stop", "end()")

Loading data, please wait...