SenseLab Home ModelDB Home

State and location dependence of action potential metabolic cost (Hallermann et al., 2012)
Accession: 144526
With this model of a layer 5 pyramidal neuron the state and location dependence of the ATP usage and the metabolic efficiency of action potentials can be analyzed. Model parameters were constrained by direct subcellular recordings at dendritic, somatic and axonal compartments.
Reference: Hallermann S, de Kock CP, Stuart GJ, Kole MH (2012) State and location dependence of action potential metabolic cost in cortical pyramidal neurons Nat Neurosci., epub before print [PubMed]
Citations  Citation Browser
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:  Neocortex;
Cell Type(s):  Neocortical pyramidal neuron: deep;  
Channel(s):  I Na,t; I K;  
Gap Junctions:  
Receptor(s):  
Gene(s):  
Transmitter(s):  
Simulation Environment:  Neuron;
Model Concept(s):  Action Potential Initiation; Action Potentials;
Implementer(s):  Hallermann, Stefan [s.hallermann at eni-g.de]; Kole, Maarten [m.kole at nin.knaw.nl];
Search NeuronDB for information about:  Neocortical pyramidal neuron: deep; I Na,t; I K;
\
HallermannEtAl2012
readme.txt
Cad.mod
CaH.mod
CaT.mod
charge.mod
h.mod
Kca.mod
Kv.mod
Kv1_axonal.mod
Kv7.mod
na8st.mod
nax8st.mod
charge.hoc
mosinit.hoc
Cell parameters.hoc
28_04_10_num19.hoc
all_28_04_10_num19.ses
                            
This model reproduces Figure 6a and 6b and Supplementary Figure 5b of
the article

"State and location dependence of action potential metabolic cost in
cortical pyramidal neurons". Nature Neuroscience (2012)
doi:10.1038/nn.3132

When the folder is compiled using the mknrndll script the simulation
can be started with "mosinit.hoc".

The newly developed model of the axonal K+ currents (Kv1_axonal.mod)
was based on step evoked K+ currents including activation,
deactivation and inactivation protocols as well as steady state
curves. To describe the activation, eight Hodgkin-Huxley gates were
used (n^8). To describe the biphasic inactivation and recovery from
inactivation, two populations of channels with different
Hodgkin-Huxley inactivation gates were used (n^8*h1^+n^8*h2).

To analyze the AP efficiency and the ATP usage on a subcellular level,
the Na+/K+ charge overlap, the excess Na+ influx and the ATP usage
were calculated as RANGE variables (charge.mod).

Executing "ch()" in the terminal shell (see charge.hoc) will print
average ATP usage per cellular compartment as shown below (some of the
information is exported into text files).

-------------------------------------------------------------------------

The Na influx was analyzed from 517.000000 to 535.000000 ms
Required ATP molecules to pump Na ions out:
Total: 1.998503e+08
Soma: 1.416016e+07
Basal dendrite: 4.609883e+07
Apical dendrite: 9.149446e+07
Axon Initial Segment: 2.498228e+07
Axon collaterals: 2.095502e+07
Myelin: 7.227968e+05
Node: 1.436733e+06

-------------------------------------------------------------------------
-------------------------------------------------------------------------

Required ATP molecules to pump Na ions out per um^2:
Total: 2.818287e+03
Soma: 6.810007e+03
Basal dendrite: 2.531221e+03
Apical dendrite: 2.213969e+03
Axon Initial Segment: 5.531797e+04
Axon collaterals: 3.932075e+03
Myelin: 2.081271e+02
Node: 3.517895e+04
-------------------------------------------------------------------------
Total area (um^2): 7.091197e+04
Soma area (um^2): 2.079316e+03
-------------------------------------------------------------------------

For questions please contact s.hallermann@eni-g.de or
m.kole@nin.knaw.nl


June 12th 2012, Stefan Hallermann


ModelDB Home  SenseLab Home   Help
Questions, comments, problems? Email the ModelDB Administrator
How to cite ModelDB
This site is Copyright 2012 Shepherd Lab, Yale University