Hodgkin-Huxley model of persistent activity in prefrontal cortex neurons (Winograd et al. 2008)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:113997
The paper demonstrate a form of graded persistent activity activated by hyperpolarization. This phenomenon is modeled based on a slow calcium regulation of Ih, similar to that introduced earlier for thalamic neurons (see Destexhe et al., J Neurophysiol. 1996). The only difference is that the calcium signal is here provided by the high-threshold calcium current (instead of the low-threshold calcium current in thalamic neurons).
Reference:
1 . Winograd M, Destexhe A, Sanchez-Vives MV (2008) Hyperpolarization-activated graded persistent activity in the prefrontal cortex. Proc Natl Acad Sci U S A 105:7298-303 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell; Channel/Receptor;
Brain Region(s)/Organism: Prefrontal cortex (PFC);
Cell Type(s):
Channel(s): I Na,t; I L high threshold; I K; I M; I h;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Activity Patterns; Ion Channel Kinetics;
Implementer(s): Destexhe, Alain [Destexhe at iaf.cnrs-gif.fr];
Search NeuronDB for information about:  I Na,t; I L high threshold; I K; I M; I h;
/
HAGPA
README.html
APCounter2.mod *
Cadynamics.mod *
HH2.mod *
ICaL.mod *
Ih.mod *
IKM.mod *
ipulse3.mod *
pasi.mod *
demo_HPGA_non-saturating.hoc
demo_HPGA_non-saturating_noIh.hoc
demo_HPGA_saturating.hoc
init.hoc
mosinit.hoc *
screenshot.jpg
                            
COMMENT
Hines' counter for action potentials
after a run, n indicates the number of action potentials that crossed
some threshold value (thresh)
AD: store the action potential times in vector "times"; 
calculate ISI and instantaneous rate
ENDCOMMENT

INDEPENDENT {t FROM 0 TO 1 WITH 1 (ms)}

DEFINE MAXAP 20000

NEURON {
	POINT_PROCESS APCounter2
	RANGE n, thresh, times, isi, rate
}

UNITS {
	(mV) = (millivolt)
}

PARAMETER {
	n
	thresh = 0 (mV)
}

ASSIGNED {
	firing
	times[MAXAP]	(ms)
	isi		(ms)
	rate		(1/s)
}

INITIAL {
	n = 0
	firing = 0
	isi = 0
	rate = 0
	check()
}

BREAKPOINT {
	SOLVE check
}

PROCEDURE check() {
	if (v >= thresh && !firing) {
		firing = 1
		times[n] = t
		n = n + 1
		if(n>1) {
			isi = t-times[n-2]
			rate = 1000/isi
		}
	}
	if (firing && v < thresh) {
		firing = 0
	}
}

Loading data, please wait...