Zebrafish Mauthner-cell model (Watanabe et al 2017)

 Download zip file 
Help downloading and running models
The NEURON model files encode the channel generator and firing simulator for simulating development and differentiation of the Mauthner cell (M-cell) excitability in zebrafish. The channel generator enables us to generate arbitrary Na+ and K+ channels by changing parameters of a Hodgkin-Huxley model under emulation of two-electrode voltage-clamp recordings in Xenopus oocyte system. The firing simulator simulates current-clamp recordings to generate firing patterns of the model M-cell, which are implemented with arbitrary-generated basic Na+ and K+ conductances and low-threshold K+ channels Kv7.4/KCNQ4 and sole Kv1.1 or Kv1.1 coexpressed with Kvbeta2.
1 . Watanabe T, Shimazaki T, Oda Y (2017) Coordinated Expression of Two Types of Low-Threshold K+ Channels Establishes Unique Single Spiking of Mauthner Cells among Segmentally Homologous Neurons in the Zebrafish Hindbrain. eNeuro [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type:
Brain Region(s)/Organism: Brainstem;
Cell Type(s): Mauthner cell;
Channel(s): I Potassium; I A; I_KLT; I_KHT; I M; I Sodium;
Gap Junctions:
Gene(s): Kv1.1 KCNA1; Kv7.4 KCNQ4; Kvb2 KCNAB2;
Simulation Environment: NEURON;
Model Concept(s): Spike Frequency Adaptation; Bursting; Ion Channel Kinetics;
Implementer(s): Watanabe, Takaki [wtakaki at m.u-tokyo.ac.jp];
Search NeuronDB for information about:  I A; I M; I Sodium; I Potassium; I_KHT; I_KLT;
// Neuron.hoc creates a model Mauthner cell (M-cell) simplified with a single
// cylindrical compartment shown in Materials and Methods. 
// specific resistance 1000 ohm*cm2, specific conductance 0.001 mho/cm2, 0.001*1E9 nS/cm2
// Takaki Watanabe
// wtakaki@m.u-tokyo.ac.jp

create soma
soma {
	cm=1 // 1.0 uF/cm2 = 0.01 pF/um2

proc totalcap(){
somaarea = totcap*1E-6/cm // convert 1pF to 1E-6 uF somaarea cm2
lstd = 1E4*sqrt(somaarea/PI) // 1cm to 1E4 um
soma.L = int(lstd)
soma.diam = int(lstd)


func nstomho() {
	return (1E-9*$1/somaarea) // nS to mho/cm2

proc capcalL(){
somaarea = PI*(1E-4*soma.L)*(1E-4*soma.diam)
totcap = int(1E6*somaarea*cm)

proc capcalD(){
somaarea = PI*(1E-4*soma.L)*(1E-4*soma.diam)
totcap = int(1E6*somaarea*cm)

proc insertoriginal(){	
	insert kht 
	insert na 	
	insert ka 	
	insert leak
	insert kcnq 
	insert kcna
    insert kcnab2
	ek_kht = -90
	erev_leak = -85
proc uninsertall(){	
	uninsert kht 
	uninsert na 
	uninsert ka 
	uninsert leak
    uninsert hh
	uninsert kcnq
	uninsert kcna
	uninsert kcnab2

access soma

Loading data, please wait...