State and location dependence of action potential metabolic cost (Hallermann et al., 2012)

 Download zip file   Auto-launch 
Help downloading and running models
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:
1 . 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 15:1007-14 [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): Neocortex L5/6 pyramidal GLU cell;
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 [hallermann at medizin.uni-leipzig.de]; Kole, Maarten [m.kole at nin.knaw.nl];
Search NeuronDB for information about:  Neocortex L5/6 pyramidal GLU cell; 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
28_04_10_num19.hoc *
all_28_04_10_num19.ses *
Cell parameters.hoc
charge.hoc *
mosinit.hoc *
                            
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

20140226 ModelDB Administrator: To prevent collision with NEURON's
"charge" the "charge" mechanism has been changed to "charge_" with
changes in the hoc code.  The charge in names of files was left
unchanged.