Disentangling astroglial physiology with a realistic cell model in silico (Savtchenko et al 2018)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:243508
"Electrically non-excitable astroglia take up neurotransmitters, buffer extracellular K+ and generate Ca2+ signals that release molecular regulators of neural circuitry. The underlying machinery remains enigmatic, mainly because the sponge-like astrocyte morphology has been difficult to access experimentally or explore theoretically. Here, we systematically incorporate multi-scale, tri-dimensional astroglial architecture into a realistic multi-compartmental cell model, which we constrain by empirical tests and integrate into the NEURON computational biophysical environment. This approach is implemented as a flexible astrocyte-model builder ASTRO. As a proof-of-concept, we explore an in silico astrocyte to evaluate basic cell physiology features inaccessible experimentally. ..."
Reference:
1 . Savtchenko LP, Bard L, Jensen TP, Reynolds JP, Kraev I, Medvedev N, Stewart MG, Henneberger C, Rusakov DA (2018) Disentangling astroglial physiology with a realistic cell model in silico. Nat Commun 9:3554 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Glia;
Brain Region(s)/Organism: Hippocampus;
Cell Type(s): Astrocyte;
Channel(s): I Calcium; I Potassium; Kir;
Gap Junctions: Gap junctions;
Receptor(s):
Gene(s):
Transmitter(s): Glutamate;
Simulation Environment: NEURON; MATLAB; C or C++ program;
Model Concept(s): Calcium waves; Calcium dynamics; Potassium buffering; Volume transmission; Membrane Properties;
Implementer(s): Savtchenko, Leonid P [leonid.savtchenko at ucl.ac.uk];
Search NeuronDB for information about:  I Calcium; I Potassium; Kir; Glutamate;
COMMENT
Longitudinal diffusion of potassium (no buffering)
Savtchenko et al., 2018
ENDCOMMENT

NEURON {
	SUFFIX kdifl
	USEION k READ ik, ki WRITE ki
	RANGE Dk, ki0, iextra
}

UNITS {
	
	(mM) = (milli/liter)
	(um) = (micron)
	FARADAY = (faraday) (coulomb)
	PI = (pi) (1)
	
}

INITIAL {
	
	ki = ki0
	
	ka = ki
}
PARAMETER {
    ki0 = 110 (mM)
	Dk = 0.6 (micron2/ms)
	iextra = 0 (milliamp/cm2)
	
}

ASSIGNED {
	ik (milliamp/cm2)
	
	diam (um)
	ki       (mM)
}

STATE {
	ka (mM)
}

BREAKPOINT {
	SOLVE conc METHOD sparse
}

KINETIC conc {
	COMPARTMENT PI*diam*diam/4 {ka}
	LONGITUDINAL_DIFFUSION Dk*diam*diam {ka}
	: LONGITUDINAL_DIFFUSION Dk {ka}
	~ ka << (-(ik-iextra)/(FARADAY)*PI*diam*(1e4))
	ki = ka
}

Loading data, please wait...