Drosophila T4 neuron (Gruntman et al 2018)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:239435
Passive, multi-compartment conductance-based model of a T4 cell. The model reproduces the neuron's response to moving stimuli via integration of spatially offset fast excitatory and slow inhibitory inputs.
Reference:
1 . Gruntman E, Romani S, Reiser MB (2018) Simple integration of fast excitation and offset, delayed inhibition computes directional selectivity in Drosophila. Nat Neurosci 21:250-257 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type:
Brain Region(s)/Organism: Drosophila;
Cell Type(s): Fly lobula plate T4 neuron;
Channel(s):
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Invertebrate;
Implementer(s): Gruntman, Eyal [ gruntmane at janelia.hhmi.org]; Romani, Sandro [ romanis at janelia.hhmi.org];
// compute time courses of conductances and setting the parameters at all synapses

objref condvec_e[num_seclamp_e], condvec_i[num_seclamp_i]


k = 0
count_e = 0
count_i = 0
count_active_e = 0
count_active_i = 0


// defined in flash_start
if (zero_inhib == 1) {
	rel_dur = 0
} else {
	rel_dur = tstop + 1
}

forsec RelDend {


	gpe = a_e * exp( -(mu_e - x_star.x[k])^2 / (2*sigma_e^2) )
	gpi = a_i * exp( -(mu_i - x_star.x[k])^2 / (2*sigma_i^2) )


// activating only the SEClamps that lie within the relevant x_star range (defined in my_start)
    if(selected_section_e.x[k]==1) {

        condvec_e[count_e] = create_cond_vec( tvec, t0, t1, taur_e, taud_e, gpe)
        i_e[count_e].dur1 = 0
    	i_e[count_e].amp1 = 0
        if(x_bound_val.x[k]-1 == ind_x_stim) {

            i_e[count_e].dur1 = tstop + 1
            condvec_e[count_e].play(&i_e[count_e].rs, tvec, 1)
			count_active_e = count_active_e + 1

        }
        count_e = count_e + 1

    }

    if(selected_section_i.x[k]==1) {
        condvec_i[count_i] = create_cond_vec( tvec, t0, t1, taur_i, taud_i, gpi)
        i_i[count_i].dur1 = 0
    	i_i[count_i].amp1 = -70
        if(x_bound_val.x[k]-1 == ind_x_stim) {

            i_i[count_i].dur1 = rel_dur
            condvec_i[count_i].play(&i_i[count_i].rs, tvec, 1)
			count_active_i = count_active_i + 1


        }
        count_i = count_i + 1
    }

	k = k + 1
}

Loading data, please wait...