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

 Download zip file 
Help downloading and running models
Accession:58582
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.
Reference:
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;
Channel(s):
Gap Junctions:
Receptor(s): GabaA; GabaB; AMPA;
Gene(s):
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 bgu.ac.il];
Search NeuronDB for information about:  Thalamus geniculate nucleus/lateral principal GLU cell; Thalamus reticular nucleus GABA cell; GabaA; GabaB; AMPA; Gaba; Glutamate;
#!/bin/sh
# Plotting trajectories for all parameter sets, locked to the stimulus
# TC (VPM and POm) and RE populations

pra=tlp
prg=plp

fla=d1
flb=d2

grfl=ptrint.gr

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

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

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 > ptrint.xx.ult.in
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 > ptrint.xx.upt.in
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.in ptrint.xx.ult.t5 ptrint.xx.upt.in 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...