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;
/* Integrate the activity trace  - modeling results */
#include <stdio.h>
#include <math.h>
#include <string.h>
#include "nr.h"
#include "intnor.h"

main(int argc, char*argv[])
{
  fl_st fl;
  double dat_ar[Mar][2], ff, xx, nor_ar[Mar], scale_factor;
  int nar, iar;
  char prefix[7]="asked", suffix[4]="ult", file_name[30];

  if (argc >= 3) strcpy(prefix,argv[1]);
  if (argc >= 3) strcpy(suffix,argv[2]);
  if (argc >= 4) sscanf(argv[3], "%lf", &scale_factor);
  printf("scale_factor=%lf\n", scale_factor);

  fl.dat = fopen(strcat(strcat(strcat(strcpy(file_name,prefix),".xx."),
           suffix),".in"),"r");
  fl.out = fopen(strcat(strcpy(file_name,"intn.out."),suffix), "w");
  fl.res = fopen(strcat(strcpy(file_name,"intn.res."),suffix), "w");

  nar = 0;
  while (fscanf(fl.dat, "%lf %lf\n", &ff, &xx) != EOF)
  {
    nar++;
    dat_ar[nar][1] = ff;
    dat_ar[nar][2] = xx;
  }

  for (iar=1; iar<=nar; iar++)
  {
    nor_ar[iar] = dat_ar[iar][2] / dat_ar[3][2];
    fprintf(fl.res,"%lf %lf\n", dat_ar[iar][1], scale_factor * nor_ar[iar]);
  }

  fclose(fl.dat);
  fclose(fl.out);
  fclose(fl.res);
}

Loading data, please wait...