Changes of ionic concentrations during seizure transitions (Gentiletti et al. 2016)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:222321
"... In order to investigate the respective roles of synaptic interactions and nonsynaptic mechanisms in seizure transitions, we developed a computational model of hippocampal cells, involving the extracellular space, realistic dynamics of Na+, K+, Ca2+ and Cl - ions, glial uptake and extracellular diffusion mechanisms. We show that the network behavior with fixed ionic concentrations may be quite different from the neurons’ behavior when more detailed modeling of ionic dynamics is included. In particular, we show that in the extended model strong discharge of inhibitory interneurons may result in long lasting accumulation of extracellular K+, which sustains the depolarization of the principal cells and causes their pathological discharges. ..."
Reference:
1 . Gentiletti D, Suffczynski P, Gnatkovsky V, de Curtis M (2017) Changes of Ionic Concentrations During Seizure Transitions - A Modeling Study. Int J Neural Syst 27:1750004 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Synapse; Extracellular; Realistic Network;
Brain Region(s)/Organism: Entorhinal cortex;
Cell Type(s):
Channel(s): Na/K pump; KCC2; I Na, leak; I K,leak; I Na,t; I Na,p; I K; I L high threshold; I_AHP; I M; I K,Ca;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s): Ions;
Simulation Environment: NEURON;
Model Concept(s): Epilepsy; Reaction-diffusion; Synchronization; Simplified Models;
Implementer(s): Gentiletti, Damiano [gentiletti.damiano at gmail.com];
Search NeuronDB for information about:  I Na,p; I Na,t; I L high threshold; I K; I K,leak; I M; I K,Ca; Na/K pump; I_AHP; I Na, leak; KCC2; Ions;
/
GentilettiEtAl2016
readme.html
accum.mod
cal.mod
getconc.mod
ikdrd.mod
ikdrs.mod
inad.mod
inas.mod
is.mod
kahp.mod
kc.mod
kcc2.mod
km.mod
leak.mod
LinClamp.mod
nakpump.mod
nap.mod
pkdrd.mod
pkdrs.mod
pnad.mod
pnas.mod
init.hoc
mosinit.hoc *
MySession.ses
screenshot.png
                            
TITLE leak current :passive electrical properties

NEURON {
	SUFFIX leak
	USEION k READ ek WRITE ik
	USEION na READ ena WRITE ina
        NONSPECIFIC_CURRENT ifix
	RANGE gk, ik, gna, ina, gfix, ifix
        RANGE qk, qna
}

UNITS { 
	(mV) = (millivolt)
        (mA) = (milliamp)
        PI		= (pi) (1)
	FARADAY		= 96485.309 (coul/mole)
}

PARAMETER {
	gk	= 1e-5 (mho/cm2)  :potassium leak conductance
	gna	= 1e-5 (mho/cm2)  :sodium leak conductance
	gfix	= 5e-4 (mho/cm2)  :fixed leak conductance
}

ASSIGNED {
	v (mV)
        v_init (mV)
	ik (mA/cm2)
	ek (mV)
	ina (mA/cm2)
	ena (mV)
	ifix (mA/cm2) :fixed leak current
        diam (um)
}

BREAKPOINT {
	ik = gk*(v-ek)
	ina = gna*(v-ena)
	:ifix = gfix*(v-v_init)
	ifix = gfix*(v+70)
        SOLVE integrate METHOD sparse
}

STATE { qk qna }

INITIAL {
	ik = gk*(v-ek)
	ina = gna*(v-ena)
	qk = 0
	qna = 0
}

KINETIC integrate {
	COMPARTMENT diam*diam*PI/4 { qna qk }
	~ qna << ((-ina*diam)*PI*(1e4)/FARADAY )
	~ qk  << ((-ik*diam)*PI*(1e4)/FARADAY )
}

Loading data, please wait...