Zebrafish Mauthner-cell model (Watanabe et al 2017)

 Download zip file   Auto-launch 
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...