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;
TITLE inward rectifier potassium (Kir4) channel

COMMENT
This is kinetics of Kir4.1 channels


ENDCOMMENT
NEURON {
    SUFFIX kir4 			
    USEION k READ ko, ek WRITE ik	
    RANGE  ik, gkir, NormK
        
}

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

PARAMETER {
    v 		(mV)
    va1 = -14.83 (mV) 	
	va2 = -105.82 (mV) : 34 (mV)
	va3 = 19.23 (mV)
    gkir = 1.44e-02  (mS/cm2) 
    ek = -70 (mV)
    NormK = 0.81 
	
}

ASSIGNED {
    ik      (mA/cm2)
       
   
	ko      (mM)
}



BREAKPOINT {
        
        ik = (0.001)*gkir * ( v - ek*NormK - va1) *sqrt(((ko)/(1 (mM)))/(1+exp((v-ek*NormK-va2)/va3)))		: calculate ik 
		: printf("v: %g, ko: %g, va2: %g\n", v, ko, va2)
}


Loading data, please wait...