A network model of the vertebrate retina (Publio et al. 2009)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:124063
In this work, we use a minimal conductance-based model of the ON rod pathways in the vertebrate retina to study the effects of electrical synaptic coupling via gap junctions among rods and among AII amacrine cells on the dynamic range of the retina. The model is also used to study the effects of the maximum conductance of rod hyperpolarization activated current Ih on the dynamic range of the retina, allowing a study of the interrelations between this intrinsic membrane parameter with those two retina connectivity characteristics.
Reference:
1 . Publio R, Oliveira RF, Roque AC (2009) A computational study on the role of gap junctions and rod Ih conductance in the enhancement of the dynamic range of the retina. PLoS One 4:e6970 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism:
Cell Type(s): Retina ganglion GLU cell; Retina photoreceptor cone GLU cell; Retina bipolar GLU cell;
Channel(s):
Gap Junctions: Gap junctions;
Receptor(s):
Gene(s):
Transmitter(s): Glutamate;
Simulation Environment: NEURON;
Model Concept(s): Sensory processing;
Implementer(s): Publio, Rodrigo [publio at oist.jp];
Search NeuronDB for information about:  Retina ganglion GLU cell; Retina photoreceptor cone GLU cell; Retina bipolar GLU cell; Glutamate;
/
PublioEtAl2009
README.html
A2hh_k.mod
A2hh_na.mod
Bip_Ca.mod
Bip_Cad.mod
Bip_ih.mod
Bip_Ka.mod
Bip_Kv.mod
Cone_CPR.mod
Cone_ih.mod
Cone_Kv.mod
Ganglion_hh.mod *
gap.mod
IinjLT.mod
IinjLT_cone.mod
IinjLTDim.mod *
Rod_Ca.mod
Rod_Cad.mod
Rod_Clca.mod
Rod_ih.mod
Rod_Kca.mod
Rod_Kv.mod
Rod_Kx.mod
Rod_leak.mod
syn_bip_gan.mod
syn_rod_bip.mod
A2.tem
Bip.tem
Cone.tem
createcells.hoc
Ganglion.tem
gap.hoc *
init.hoc
mosinit.hoc *
netconnection.hoc
parameters.hoc
Rod.tem
screenshot1.jpg
screenshot2.jpg
session.ses
                            
INDEPENDENT {t FROM 0 TO 1 WITH 1 (ms)}

NEURON {
	SUFFIX IKv
	USEION k READ ek WRITE ik
	RANGE gkvbar
	RANGE m_inf, h_inf
	RANGE tau_m, tau_h
	RANGE m_exp, h_exp

}


UNITS {
	(molar) = (1/liter)
	(mM) = (millimolar)
	(mA) = (milliamp)
	(mV) = (millivolt)

}

PARAMETER {
	gkvbar	= 0.000318	(mho/cm2)
	ek		(mV)
	dt              (ms)
	v               (mV)

}

STATE {
	m h 
}

INITIAL {
: The initial values were determined at a resting value of -66.3232 mV in a single-compartment
: at -60 mV
        m = 0.0345
        h = 0.8594
}

ASSIGNED {
	 ik	(mA/cm2)
       m_inf 
	 h_inf
	 tau_m 
	 tau_h
	 m_exp 
       h_exp

}

BREAKPOINT {
	SOLVE states
	ik = gkvbar * m*m*m *h* (v - ek)

}

PROCEDURE states() {	: exact when v held constant
	evaluate_fct(v)
	m = m + m_exp * (m_inf - m)
	h = h + h_exp * (h_inf - h)

	VERBATIM
	return 0;
	ENDVERBATIM

}

UNITSOFF

PROCEDURE evaluate_fct(v(mV)) { LOCAL am,bm,ah,bh
	
:IA channel

 am = (400) / ((exp(-(v-15)/36)) + 1)
 bm = (exp((-v/13)))
 tau_m = 1 / (am + bm)
 m_inf = am * tau_m
 m_exp = 1 - exp(-dt/tau_m)

 ah = 0.0003 * (exp((-v/7)))
 bh = ((80) / ((exp(-(v+115)/15)) + 1))+0.02
 tau_h = 1 / (ah + bh)
 h_inf = ah * tau_h
 h_exp = 1 - exp(-dt/tau_h)

}

UNITSON

Loading data, please wait...