Spiny neuron model with dopamine-induced bistability (Gruber et al 2003)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:39949
These files implement a model of dopaminergic modulation of voltage-gated currents (called kir2 and caL in the original paper). See spinycell.html for details of usage and implementation. For questions about this implementation, contact Ted Carnevale (ted.carnevale@yale.edu)
Reference:
1 . Gruber AJ, Solla SA, Surmeier DJ, Houk JC (2003) Modulation of striatal single units by expected reward: a spiny neuron model displaying dopamine-induced bistability. J Neurophysiol 90:1095-114 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell;
Brain Region(s)/Organism:
Cell Type(s): Neostriatum medium spiny direct pathway GABA cell;
Channel(s): I K,leak; I h; I A, slow;
Gap Junctions:
Receptor(s): D1; Dopaminergic Receptor;
Gene(s): Cav1.3 CACNA1D; D1 DRD1A;
Transmitter(s): Dopamine;
Simulation Environment: NEURON;
Model Concept(s): Intrinsic plasticity;
Implementer(s): Carnevale, Ted [Ted.Carnevale at Yale.edu];
Search NeuronDB for information about:  Neostriatum medium spiny direct pathway GABA cell; D1; Dopaminergic Receptor; I K,leak; I h; I A, slow; Dopamine;
COMMENT
Implements the "spike" mechanism used by Gruber et al. 2003.
Watch this mechanism's spike variable with a NetCon that has 
threshold 0.5, and record that NetCon's spike event times.

NOTES:
  Needs fixed time step integration, 
    won't work with adaptive integration.
  Needs NEURON v. 5.5 or later.
ENDCOMMENT

NEURON {
	POINT_PROCESS GSpike
	RANGE vt, vh, vc, spike, minisi
}

UNITS {
	(mV) = (millivolt)
}

PARAMETER {
	vt = -58 (mV)
	vh = -55 (mV)
	vc = 2.5 (mV)
	minisi = 20 (ms)
	spikewidth = 1 (ms)
}

ASSIGNED {
	v (mV)
	tp (ms)
	spike (1)
}

INITIAL {
	tp = -1e9 (ms)
	spike = 0
}

AFTER SOLVE {
	if (at_time(tp+spikewidth)) {
		spike = 0
	}
	if (v>vt) {
		if (t>tp+minisi*(1+exp(-(v-vh)/vc))) {
			spike = 1
			tp = t
		}
	}
}

Loading data, please wait...