ModelDB is moving. Check out our new site at The corresponding page is

Coding of stimulus frequency by latency in thalamic networks (Golomb et al 2005)

 Download zip file 
Help downloading and running models
The paper presents models of the rat vibrissa processing system including the posterior medial (POm) thalamus, ventroposterior medial (VPm) thalamus, and GABAB- mediated feedback inhibition from the reticular thalamic (Rt) nucleus. A clear match between the experimentally measured spike-rates and the numerically calculated rates for the full model occurs when VPm thalamus receives stronger brainstem input and weaker GABAB-mediated inhibition than POm thalamus.
1 . Golomb D, Ahissar E, Kleinfeld D (2006) Coding of stimulus frequency by latency in thalamic networks through the interplay of GABAB-mediated feedback and stimulus shape. J Neurophysiol 95:1735-50 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism: Thalamus;
Cell Type(s): Thalamus geniculate nucleus/lateral principal GLU cell; Thalamus reticular nucleus GABA cell;
Gap Junctions:
Receptor(s): GabaA; GabaB; AMPA;
Transmitter(s): Gaba; Glutamate;
Simulation Environment: C or C++ program;
Model Concept(s): Simplified Models; Rate-coding model neurons;
Implementer(s): Golomb, David [golomb at];
Search NeuronDB for information about:  Thalamus geniculate nucleus/lateral principal GLU cell; Thalamus reticular nucleus GABA cell; GabaA; GabaB; AMPA; Gaba; Glutamate;
# Plotting trajectories for all parameter sets, locked to the stimulus
# TC (VPM and POm) and RE populations



tlp.ex $fla
plp.ex $fla
stm.ex $fla

cat > ptrint.xx.fsed <<EOF
s/INPUT ff 2.0 5.0 8.0 11.0/INPUT ff parmin=1.0 parmax=11.5 npar=21/

sed -f ptrint.xx.fsed $pra.n.$fla > $pra.n.$flb

tlp.ex $flb
plp.ex $flb

awk '{if ($2 == "1") print $1, 1000.0 * $7}' $prg.res.$flb >
awk '{if ($2 == "1") print $1, 1000.0 * $6}' $prg.res.$flb > ptrint.xx.ult.t5
awk '{if ($2 == "2") print $1, 1000.0 * $7}' $prg.res.$flb >
awk '{if ($2 == "2") print $1, 1000.0 * $6}' $prg.res.$flb > ptrint.xx.upt.t5

awk '{if ($1 != "#") print $1, $2}' delay_exp > delay.exp.vpm.xx
awk '{if ($1 != "#") print $1, $3}' delay_exp > delay.exp.pom.xx

awk '{print $1, $2 * 40.0}' intg.res.a3 > intg.res.a3.xx
awk '{print $1, $2 * 40.0}' intg.res.a4 > intg.res.a4.xx

intnor.ex ptrint ult 40.0
intnor.ex ptrint upt 40.0

xmgrace \
        -graph 0 plp.xx.$fla.1.ult plp.xx.$fla.2.ult \
                 plp.xx.$fla.3.ult plp.xx.$fla.4.ult \
        -graph 1 plp.xx.$fla.1.upt plp.xx.$fla.2.upt \
                 plp.xx.$fla.3.upt plp.xx.$fla.4.upt \
        -graph 2 stm.flt.$fla \
        -graph 3 stm.fpt.$fla \
        -graph 4 intn.res.ult intg.res.a3.xx \
                 ptrint.xx.ult.t5 delay.exp.vpm.xx \
        -graph 5 intn.res.upt intg.res.a4.xx \
                 ptrint.xx.upt.t5 delay.exp.pom.xx \
        -graph 6 plp.xx.$fla.1.ura plp.xx.$fla.2.ura \
                 plp.xx.$fla.3.ura plp.xx.$fla.4.ura \
        -graph 7 plp.xx.$fla.1.urb plp.xx.$fla.2.urb \
                 plp.xx.$fla.3.urb plp.xx.$fla.4.urb \
        -hdevice EPS -p $grfl -printfile ptrint.$fla.eps

/bin/rm $pra.out.$fla $pra.tmp.$fla $pra.col.$fla
/bin/rm $pra.out.$flb $pra.tmp.$flb $pra.col.$flb
/bin/rm plp.xx.$fla.*.*
/bin/rm stm.flt.$fla stm.fpt.$fla stm.ffr.$fla
/bin/rm ptrint.xx.fsed
/bin/rm plp.xx.$flb.*.*
/bin/rm ptrint.xx.ult.t5 ptrint.xx.upt.t5
/bin/rm delay.exp.vpm.xx delay.exp.pom.xx
/bin/rm intg.res.a3.xx  intg.res.a4.xx
/bin/rm intn.out.ult intn.res.ult intn.out.upt intn.res.upt

Loading data, please wait...