A Model Circuit of Thalamocortical Convergence (Behuret et al. 2013)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:150240
“… Using dynamic-clamp techniques in thalamic slices in vitro, we combined theoretical and experimental approaches to implement a realistic hybrid retino-thalamo-cortical pathway mixing biological cells and simulated circuits. … The study of the impact of the simulated cortical input on the global retinocortical signal transfer efficiency revealed a novel control mechanism resulting from the collective resonance of all thalamic relay neurons. We show here that the transfer efficiency of sensory input transmission depends on three key features: i) the number of thalamocortical cells involved in the many-to-one convergence from thalamus to cortex, ii) the statistics of the corticothalamic synaptic bombardment and iii) the level of correlation imposed between converging thalamic relay cells. In particular, our results demonstrate counterintuitively that the retinocortical signal transfer efficiency increases when the level of correlation across thalamic cells decreases. …”
Reference:
1 . Behuret S, Deleuze C, Gomez L, Fregnac Y and Bal T (2013) Cortically-controlled population stochastic facilitation as a plausible substrate for guiding sensory transfer across the thalamic gateway PLoS Computational Biology 9(12):e1003401
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism: Neocortex; Thalamus; Retina;
Cell Type(s): Thalamus geniculate nucleus/lateral principal neuron; Thalamus reticular nucleus cell; Neocortex U1 L5B pyramidal pyramidal tract GLU cell; Retina ganglion cell; Thalamus lateral geniculate nucleus interneuron;
Channel(s): I Na,t; I T low threshold; I K; I M;
Gap Junctions:
Receptor(s): GabaA; AMPA;
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Synaptic Convergence;
Implementer(s): Behuret, Sebastien [behuret at unic.cnrs-gif.fr];
Search NeuronDB for information about:  Thalamus geniculate nucleus/lateral principal neuron; Thalamus reticular nucleus cell; Retina ganglion cell; Neocortex U1 L5B pyramidal pyramidal tract GLU cell; GabaA; AMPA; I Na,t; I T low threshold; I K; I M;
/
TCconvergenceModel
README.html
cadecay.mod *
ConductancePattern.mod
ConstantCurrent.mod
hh2.mod *
IM.mod
IT.mod
ITGHK.mod
RandomGenerator.mod
RetinalInput.mod
SineWaveCurrent.mod
SynapticNoise.mod
Demo.hoc
DEMO.png
Geometry.hoc
GUI.hoc
mosinit.hoc
Recording.hoc
Run.hoc
screenshot.png
Simulation.hoc
Templates.hoc
                            
//
// Templates.hoc - Cortical pyramidal and thalamocortical relay cell model templates
//
// From Alain Destexhe, original templates can be found in
// http://cns.iaf.cnrs-gif.fr/files/CORTEX.ZIP
// and
// http://cns.iaf.cnrs-gif.fr/files/DENDTC.ZIP
//
// Thalamocortical convergence studies
// Sébastien Béhuret, UNIC/CNRS Paris, 2009
//

begintemplate PyramidalCell
	public Soma
	create Soma

	proc init() { local v_potassium, v_sodium
		v_potassium = -100		// potassium reversal potential 
		v_sodium = 50			// sodium reversal potential 

		create Soma

		Soma {
			Ra = 100		// geometry 
			nseg = 1
			diam = 96
			L = 96			// so that area is about 29000 um2
			cm = 1

			insert pas		// leak current 
			e_pas = -70
			g_pas = 0.0001		// Rin = 34 Meg

			// conversion with McC units: 
			// g(S/cm2) = g(nS)*1e-9/29000e-8
			//	    = g(nS) * 3.45e-6

			insert hh2		// Hodgin-Huxley INa and IK 
			ek = v_potassium
			ena = v_sodium
			vtraub_hh2 = -55	// Resting Vm, BJ was -55
			gnabar_hh2 = 0.05	// McCormick=15 muS, thal was 0.09
			gkbar_hh2 = 0.005	// spike duration of pyr cells

			insert im		// M current 
			taumax_im = 1000
			gkbar_im = 7e-5		// Diego's IM (copyrighted)

		//	insert iahp
		//	alpha_iahp = 1.25e9
		//	beta_iahp = 0.001
		//	gkbar_iahp = 1.7e-5	// McCormick=5 nS

		//	insert cad		// calcium decay
		//	depth_cad = 1		// McCormick= 0.1 um
		//	taur_cad = 5		// McCormick=1 ms !!!
		//	cainf_cad = 2.4e-4	// McCormick=0
		//	kt_cad = 0		// no pump
		}
	}
endtemplate PyramidalCell

begintemplate RelayCell
	public Soma
	create Soma

	proc init() {
  		create Soma

  		Soma {
			nseg = 1
			diam = 76.58
			L = 100
			cm = 0.88
			Ra = 173

			insert pas
			g_pas = 3.79e-5
			e_pas = -76.5

			insert hh2
			ena = 50
			ek = -100
			vtraub_hh2 = -52
			gnabar_hh2 = 0.01
			gkbar_hh2 = 0.01

			insert itGHK
			cai = 2.4e-4
			cao = 2
			eca = 120
			shift_itGHK = -1
			gcabar_itGHK = 0.0002
			qm_itGHK = 2.5
			qh_itGHK = 2.5
			pcabar_itGHK = 8e-5

			insert cad
			depth_cad = 0.1
			kt_cad = 0
			kd_cad = 1e-4
			taur_cad = 5
			cainf_cad = 2.4e-4
		}
	}
endtemplate RelayCell