Layer V PFC pyramidal neuron used to study persistent activity (Sidiropoulou & Poirazi 2012)

 Download zip file 
Help downloading and running models
"... Here, we use a compartmental modeling approach to search for discriminatory features in the properties of incoming stimuli to a PFC pyramidal neuron and/or its response that signal which of these stimuli will result in persistent activity emergence. Furthermore, we use our modeling approach to study cell-type specific differences in persistent activity properties, via implementing a regular spiking (RS) and an intrinsic bursting (IB) model neuron. ... Collectively, our results pinpoint to specific features of the neuronal response to a given stimulus that code for its ability to induce persistent activity and predict differential roles of RS and IB neurons in persistent activity expression. "
1 . Sidiropoulou K, Poirazi P (2012) Predictive features of persistent activity emergence in regular spiking and intrinsic bursting model neurons. PLoS Comput Biol 8:e1002489 [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): Neocortex L5/6 pyramidal GLU cell;
Channel(s): I Na,p; I Na,t; I L high threshold; I A; I K; I K,Ca; I CAN;
Gap Junctions:
Receptor(s): GabaA; GabaB; AMPA; NMDA; IP3;
Transmitter(s): Gaba; Glutamate;
Simulation Environment: NEURON;
Model Concept(s): Activity Patterns; Detailed Neuronal Models;
Implementer(s): Sidiropoulou, Kyriaki [sidirop at];
Search NeuronDB for information about:  Neocortex L5/6 pyramidal GLU cell; GabaA; GabaB; AMPA; NMDA; IP3; I Na,p; I Na,t; I L high threshold; I A; I K; I K,Ca; I CAN; Gaba; Glutamate;
                             TITLE Slow Ca-dependent cation current
                             :   Ca++ dependent nonspecific cation current ICAN
                             :   Differential equations
                             :   Model based on a first order kinetic scheme
                             :       + n cai <->     (alpha,beta)
                             :   Following this model, the activation fct will be half-activated at 
                             :   a concentration of Cai = (beta/alpha)^(1/n) = cac (parameter)
                             :   The mod file is here written for the case n=2 (2 binding sites)
                             :   ---------------------------------------------
                             :   Kinetics based on: Partridge & Swandulla, TINS 11: 69-72, 1988.
                             :   This current has the following properties:
                             :      - inward current (non specific for cations Na, K, Ca, ...)
                             :      - activated by intracellular calcium
                             :      - NOT voltage dependent
                             :   A minimal value for the time constant has been added
                             :   Ref: Destexhe et al., J. Neurophysiology 72: 803-818, 1994.
                             :   See also: ,

:Updated by Kiki Sidiropoulou (2010) so that dADP has slow inactivation kinetics and it is activated after 5 spikes

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

                             NEURON {
                                     SUFFIX ican
                                     USEION n READ en WRITE in VALENCE 1
                                     USEION ca READ cai
				     USEION na WRITE ina
                                     RANGE gbar, m_inf, tau_m, in, mystart
                                     GLOBAL beta, cac, taumin

                             UNITS {
                                     (mA) = (milliamp)
                                     (mV) = (millivolt)
                                     (molar) = (1/liter)
                                     (mM) = (millimolar)

                             PARAMETER {
                                     v               (mV)
                                     celsius = 36    (degC)
                                     en      = -20   (mV)            	: reversal potential
                                     cai     	     (mM)           	: initial [Ca]i
                                     gbar    = 0.0001 (mho/cm2)
                                     beta    = 0.0001 (1/ms) 	 	: backward rate constant
				     cac     = 0.0004 (mM)
				    : middle point of activation fct, for ip3 as somacar, for current injection
                                     taumin  = 0.1   (ms)            	: minimal value of time constant
                		     mystart=50 (ms)             

                             STATE {

                             ASSIGNED {
                                     in      (mA/cm2)
				     ina     (mA/cm2)
                                     tau_m   (ms)

                             BREAKPOINT { 
                                     SOLVE states METHOD cnexp
				if (t>mystart)  {     
				in = gbar * m*m * (v - en)
				ina = 0.7* in}
                             DERIVATIVE states { 

                                     m' = (m_inf - m) / tau_m

                             INITIAL {
                             :  activation kinetics are assumed to be at 22 deg. C
                             :  Q10 is assumed to be 3
                                     tadj = 3.0 ^ ((celsius-22.0)/10)

                                     m = m_inf

                             PROCEDURE evaluate_fct(v(mV),cai(mM)) {  LOCAL alpha2

                                     alpha2 = beta * (cai/cac)^2

                                     tau_m = 1 / (alpha2 + beta) / tadj
                                     m_inf = alpha2 / (alpha2 + beta)

                                     if(tau_m < taumin) { tau_m = taumin }   : min value of time cst

Loading data, please wait...