Drosophila 3rd instar larval aCC motoneuron (Gunay et al. 2015)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:152028
Single compartmental, ball-and-stick models implemented in XPP and full morphological model in Neuron. Paper has been submitted and correlates anatomical properties with electrophysiological recordings from these hard-to-access neurons. For instance we make predictions about location of the spike initiation zone, channel distributions, and synaptic input parameters.
Reference:
1 . Günay C, Sieling FH, Dharmar L, Lin WH, Wolfram V, Marley R, Baines RA, Prinz AA (2015) Distal spike initiation zone location estimation by morphological simulation of ionic current filtering demonstrated in a novel model of an identified Drosophila motoneuron. PLoS Comput Biol 11:e1004189 [PubMed]
Citations  Citation Browser
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):
Channel(s): I Na,p; I Na,t; I A; I K;
Gap Junctions:
Receptor(s): Cholinergic Receptors;
Gene(s):
Transmitter(s):
Simulation Environment: NEURON; XPP; MATLAB;
Model Concept(s):
Implementer(s): Gunay, Cengiz [cgunay at emory.edu]; Sieling, Fred [fred.sieling at gmail.com]; Prinz, Astrid [astrid.prinz at emory.edu];
Search NeuronDB for information about:  Cholinergic Receptors; I Na,p; I Na,t; I A; I K;
/
Gunay_etal_2014
neuron-model
aCC-L3-neuron.hoc
aCC-L3-neuron+electrode.xml
aCC-L3-neuron-swc.hoc
calc-impedance.hoc
chan-DmKA-Marley.hoc
chan-DmKdr-Marley.hoc
chan-DmNaP-DmNav10.hoc
chan-DmNaT-ODowd.hoc
collapse-neuron-tree.hoc
current-inj-50pA-read-mV_dt_0.025ms.bin
data-axon-tail2-axon-50um-vc-noKdr-long-back-85mV-Na_4_lines_dt_0.025000ms.bin
data-axon-tail2-axon-70um-vc-noKdr-long-back-85mV-Na_4_lines_dt_0.025000ms.bin
data-axon-tail2-axon-70um-vc-noKdr-long-back-85mV-Na-5xNaP_4_lines_dt_0.025000ms.bin
data-axon-tail2-axon-70um-vc-noKdr-long-back-85mV-Na-5xNaT_4_lines_dt_0.025000ms.bin
data-axon-tail2-axon-70um-vc-noKdr-long-back-85mV-passive_4_lines_dt_0.025000ms.bin
data-axon-tail2-chans-axon_11_lines_dt_0.025000ms.bin
data-axon-tail2-chans-axon-last_11_lines_dt_0.025000ms.bin
data-axon-tail2-chans-botdend_11_lines_dt_0.025000ms.bin
data-axon-tail2-chans-ext-axon-70um_11_lines_dt_0.025000ms.bin
data-axon-tail2-chans-in-all_11_lines_dt_0.025000ms.bin
data-i-syn-10syns-20-EPSCs-10x-10ms-VC-60mV_6_lines_dt_0.025000ms.bin
data-i-syn-4dends-50-EPSCs-10x-10ms-VC-60mV_5_lines_dt_0.025000ms.bin
data-i-vclamp-syn-dend-513-180-EPSCs-10x-1ms-saturating_2_lines_dt_0.025000ms.bin
data-syn-dend-357_2_lines_dt_0.025000ms.bin
data-syn-dend-513_2_lines_dt_0.025000ms.bin
data-syn-dend-520_2_lines_dt_0.025000ms.bin
data-syn-dend-685_2_lines_dt_0.025000ms.bin
data-v-syn-10dends-20-EPSCs-10x-10ms-noVC_6_lines_dt_0.025000ms.bin
data-v-syn-4dends-50-EPSCs-10x-10ms-noVC_6_lines_dt_0.025000ms.bin
data-v-syn-dend-513-180-EPSCs-10x-1ms-saturating-noVC_5_lines_dt_0.025000ms.bin
data-v-syn-dend-685-AP_3_lines_dt_0.025000ms.bin
exp-axon-tail2.ses
exp-axon-tail2-chans-axon.ses
exp-axon-tail2-chans-axon-last.ses
exp-axon-tail2-chans-botdend.ses
exp-axon-tail2-chans-ext-axon-50um-onlyNa.ses
exp-axon-tail2-chans-ext-axon-70um.ses
exp-axon-tail2-chans-ext-axon-70um-10alphasynapses.ses
exp-axon-tail2-chans-ext-axon-70um-10x-mimic-sustained.ses
exp-axon-tail2-chans-ext-axon-70um-10x-mimic-sustained-random.ses
exp-axon-tail2-chans-ext-axon-70um-mimic-synapses.ses
exp-axon-tail2-chans-ext-axon-70um-mimic-synapses-sustained-currents.ses
exp-axon-tail2-chans-ext-axon-70um-mimic-synapses-v-change.ses
exp-axon-tail2-chans-ext-axon-70um-onlyNa.ses
exp-axon-tail2-chans-ext-axon-70um-tomasz.ses
exp-axon-tail2-chans-in-all.ses
figures.m
fitfuncs.hoc
graph-i-vc-ext-axon.ses
iclamp-50pA.ses
IClamp-steps.ses
inc-first.ses
lincir-vclamp.hoc
lincir-vclamp.ses
NaP_NaT_data.csv
neuron-CB.ses
neuron-CB+electrode.hoc
neuron-CB-act-electrode-embed-IClamp.ses
neuron-CB-ext-axon.ses
neuron-CB-ext-axon-2pieces.ses
neuron-CB-ext-axon-2pieces-chans-axon.ses
neuron-CB-ext-axon-2pieces-chans-axon-last.ses
neuron-CB-ext-axon-2pieces-chans-botdend.ses
neuron-CB-ext-axon-2pieces-chans-ext-axon-50um-onlyNa.ses
neuron-CB-ext-axon-2pieces-chans-ext-axon-70um.ses
neuron-CB-ext-axon-2pieces-chans-ext-axon-70um-10alphasynapses.ses *
neuron-CB-ext-axon-2pieces-chans-ext-axon-70um-10x-mimic-sustained.ses *
neuron-CB-ext-axon-2pieces-chans-ext-axon-70um-mimic-synapses.ses *
neuron-CB-ext-axon-2pieces-chans-ext-axon-70um-mimic-synapses-v-change.ses *
neuron-CB-ext-axon-2pieces-chans-ext-axon-70um-onlyNa.ses
neuron-CB-ext-axon-2pieces-chans-in-all.ses
neuron-CB-pas-electrode-embed.ses
neuron-CB-pas-electrode-embed-fit-pas.ses
neuron-CB-pas-electrode-embed-fit-pas-VClamp.ses
neuron-CB-pas-electrode-embed-IClamp.ses
neuron-CB-pas-electrode-embed-test-axon-hh-chans.ses
neuron-Import3D-CellBuilder.ses
neuron-NL-CellBuilder.ses
neuron-NL-CellBuilder-pas.ses
neuron-NL-CellBuilder-pas-electrode.ses
neuron-NL-CellBuilder-pas-Na.ses
neuron-PointProcessMgr-ext-axon-2pieces-chans-ext-axon-70um-10alphasynapses.ses
nrn-fit-cap-02_dt_0.025000ms_dy_1e-9nA.bin
shape-plot.ses
SkeletonTree_ORR_aCC_48h1_NL.hoc
soma-vclamp-testbed.ses
stats.hoc
vclamp_-85_to_-25mV.ses
vclamp_soma_-60mV.ses
vclamp_soma_-60mV_syn1234.ses
vclamp_soma_-60mV_syni.ses
vclamp-family.ses
v-graph.ses
v-graph-bigger.ses
v-graph-bigger-axon-2pieces.ses
                            
// Experiment with 2-component axon extention tail, and Na-K channels
// in all compartments

// nrngui stuff and ion chans
load_file("inc-first.ses")

strdef file_name, ses_file
file_name = "ext-axon-70um-mimic-synapses-v-change"

sprint(ses_file, "neuron-CB-ext-axon-2pieces-chans-%s.ses", file_name)

//load CellBuilder
{load_file(ses_file)}


{
xpanel("RunControl", 0)
v_init = -60
xvalue("Init","v_init", 1,"stdinit()", 1, 1 )
xbutton("Init & Run","run()")
xbutton("Stop","stoprun=1")
runStopAt = 5
xvalue("Continue til","runStopAt", 1,"{continuerun(runStopAt) stoprun=1}", 1, 1 )
runStopIn = 1
xvalue("Continue for","runStopIn", 1,"{continuerun(t + runStopIn) stoprun=1}", 1, 1 )
xbutton("Single Step","steprun()")
t = 200
xvalue("t","t", 2 )
tstop = 1200
xvalue("Tstop","tstop", 1,"tstop_changed()", 0, 1 )
dt = 0.025
xvalue("dt","dt", 1,"setdt()", 0, 1 )
steps_per_ms = 40
xvalue("Points plotted/ms","steps_per_ms", 1,"setdt()", 0, 1 )
screen_update_invl = 0.1
xvalue("Scrn update invl","screen_update_invl", 1,"", 0, 1 )
realtime = 1.12
xvalue("Real Time","realtime", 0,"", 0, 1 )
xpanel(893,777)
}

// load common funcs
load_file("fitfuncs.hoc")

print_elec_cell()

// calc morph stats
load_file("stats.hoc")

load_file("v-graph-bigger-axon-2pieces.ses")

//load_file("neuron-PointProcessMgr-ext-axon-2pieces-chans-ext-axon-70um-10alphasynapses.ses")


// save state in this file
sprint(state_file_name, "state-neuron-CB+act+elec+ext-axon-2piece-chans-%s.bin", file_name)

// small adjustments

// loop and set random timing for synapses
// use NetStim and NetCon

// load state first because it overrides parameters
strdef state_w_syn
sprint(state_w_syn, "state-neuron-CB+act+elec+ext-axon-2piece-chans-%s-synapse-saved.bin", file_name)

restoreStateFromFile(state_w_syn)

// use ExpSyn object
objref ns1, ns2, ns3, ns4
objref nc1, nc2, nc3, nc4
objref syn1, syn2, syn3, syn4

// Usage: syn(synX, ncX, nsX)
proc set_syn() { $o1 = new ExpSyn(0.5)
    $o1.tau = 5
    $o1.e = 0
    $o2 = new NetStim(.5)
    $o3 = new NetCon($o2, $o1)
    $o2.interval = 5
    $o2.number = 10
    $o2.start = 10
    $o2.noise = 0
    $o3.weight = .00018
    }

// try: (w in uS)
// 685, w= 0.000150 => peak=-9.032 pA, APw=0.0034 (22.66)
// 520, w= 0.000145 => peak=-9.032 pA, APw=0.0033 (22.75)
// 357, w= 0.000175 => peak=-10.67 pA, APw=0.0028 (16.00)
// 513, w= 0.000180 => peak=-11.01 pA, APw=0.0029 (16.11)
dendrite[685] set_syn(syn1, nc1, ns1)
dendrite[520] set_syn(syn2, nc2, ns2)
dendrite[357] set_syn(syn3, nc3, ns3)
dendrite[513] set_syn(syn4, nc4, ns4)

// apply or release voltage clamp from here
load_file("vclamp_soma_-60mV_syn1234.ses")

// customize stimulus params
//nc2.start = 12.5
//nc3.start = 15
//nc4.start = 17.5

wmult = 50
ns1.weight = 0.000150*wmult
ns2.weight = 0.000145*wmult
ns3.weight = 0.000175*wmult
ns4.weight = 0.000180*wmult

// use this to save the data:
// saveLines(2, 1, 2, "data-syn-dend-520")