Ion concentration dynamics as a mechanism for neuronal bursting (Barreto & Cressman 2011)

 Download zip file 
Help downloading and running models
Accession:142630
"We describe a simple conductance-based model neuron that includes intra and extracellular ion concentration dynamics and show that this model exhibits periodic bursting. The bursting arises as the fast-spiking behavior of the neuron is modulated by the slow oscillatory behavior in the ion concentration variables and vice versa. By separating these time scales and studying the bifurcation structure of the neuron, we catalog several qualitatively different bursting profiles that are strikingly similar to those seen in experimental preparations. Our work suggests that ion concentration dynamics may play an important role in modulating neuronal excitability in real biological systems."
Reference:
1 . Barreto E, Cressman JR (2011) Ion concentration dynamics as a mechanism for neuronal bursting. J Biol Phys 37:361-73 [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: Generic;
Cell Type(s): Hodgkin-Huxley neuron;
Channel(s):
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: C or C++ program; XPP;
Model Concept(s): Bursting; Oscillations; Simplified Models; Depolarization block;
Implementer(s): Barreto, Ernest ;
#
# This is the one-cell model with fixed [K]o AND [Na]i used in
#
# E. Barreto and J.R. Cressman, "Ion Concentration Dynamics as a Mechanism for Neuronal Bursting",
# Journal of Biological Physics 37, 361-373 (2011).
#
# Link to the paper: http://www.springerlink.com/content/v52215p195159211/
# Author-generated version available at: http://arxiv.org/abs/1012.3124
#
# The variables are:
#	V=y[1]=V, the membrane voltage
#	n=y[2]=n, gating variable
#	h=y[3]=h, gating variable
#

par E_cl=-81.93864549
par Cm=1.0, g_na=100.0, g_naL=0.0175, g_k=40.0, g_kL=0.05
par g_ahp=0.01, g_clL=0.05, g_ca=0.1, phi=3.0

par Ko=4.0, Nai=18.0, I=0.0

########

alpha_n=0.01*(V+34.0)/(1.0-exp(-0.1*(V+34.0)))
beta_n=0.125*exp(-(V+44.0)/80.0)
alpha_m=0.1*(V+30.0)/(1.0-exp(-0.1*(V+30.0)))
beta_m=4.0*exp(-(V+55.0)/18.0)
alpha_h=0.07*exp(-(V+44.0)/20.0)
beta_h=1.0/(1.0+exp(-0.1*(V+14.0)))
m_inf=alpha_m/(alpha_m+beta_m)

Kin=158.0-Nai
Naout=144.0-7.0*(Nai-18.0)
E_k=26.64*log((Ko/Kin))
E_na=26.64*log((Naout/Nai))

Ina=g_na*(m_inf*m_inf*m_inf)*h*(V-E_na)+g_naL*(V-E_na)
Ik=(g_k*n*n*n*n)*(V-E_k)+g_kL*(V-E_k)
Icl=g_clL*(V-E_cl)

# differential equations

V'=(1.0/Cm)*(-Ina-Ik-Icl+I)
n'=phi*(alpha_n*(1-n)-beta_n*n)
h'=phi*(alpha_h*(1-h)-beta_h*h)

####

init V=-50,n=0.08553,h=0.96859
@ TOTAL=1000,BOUND=10000
done


Loading data, please wait...