Parameter estimation for Hodgkin-Huxley based models of cortical neurons (Lepora et al. 2011)

 Download zip file 
Help downloading and running models
Accession:136808
Simulation and fitting of two-compartment (active soma, passive dendrite) for different classes of cortical neurons. The fitting technique indirectly matches neuronal currents derived from somatic membrane potential data rather than fitting the voltage traces directly. The method uses an analytic solution for the somatic ion channel maximal conductances given approximate models of the channel kinetics, membrane dynamics and dendrite. This approach is tested on model-derived data for various cortical neurons.
Reference:
1 . Lepora NF, Overton PG,Gurney K (2012) Efficient fitting of conductance-based model neurons from somatic current clamp. J. Comp. Neuro. 32(1):1-24 [PubMed]
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:
Cell Type(s): Neocortex V1 pyramidal corticothalamic L6 cell; Neocortex V1 pyramidal intratelencephalic L2-5 cell; Neocortex fast spiking (FS) interneuron; Neocortex spiking regular (RS) neuron; Neocortex spiking low threshold (LTS) neuron;
Channel(s): I Na,t; I L high threshold; I T low threshold; I K; I M;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: GENESIS; MATLAB;
Model Concept(s): Parameter Fitting; Simplified Models; Parameter sensitivity;
Implementer(s): Lepora, Nathan [n.lepora at shef.ac.uk];
Search NeuronDB for information about:  Neocortex V1 pyramidal corticothalamic L6 cell; Neocortex V1 pyramidal intratelencephalic L2-5 cell; I Na,t; I L high threshold; I T low threshold; I K; I M;
  
/
demo_modelDB
functions_v0.4
README.txt
fig1.m
fig1_tar.m
                            
Simulation and fitting of two-compartment (active soma, passive dendrite) 
for different classes of cortical neurons.

The fitting technique indirectly matches neuronal currents derived from 
somatic membrane potential data rather than fitting the voltage traces 
directly. The method uses an analytic solution for the somatic ion channel 
maximal conductances given approximate models of the channel kinetics, 
membrane dynamics and dendrite. This approach is tested on model-derived 
data for various cortical neurons. 

The fitting method and models are described in

Lepora N., Overton P. and Gurney K.
Efficient fitting of conductance-based model neurons from somatic current 
clamp. 
Journal Computational Neuroscience 2011

The models are two compartment versions of the single compartment cortical 
neuron models described in (see also modelDB entry)

Pospischil, M., Toledo-Rodriguez, M., Monier, C., Piwkowska, Z., 
Bal, T., Fregnac, Y., Markram, H. and Destexhe, A.
Minimal Hodgkin-Huxley type models for different classes of
cortical and thalamic neurons.
Biological Cybernetics 99: 427-441, 2008.

Intrinsic currents: INa, IKd for action potentials, IM for spike-frequency 
adaptation, ICaL for high-threshold calcium current, ICaT for the 
low-threshold calcium current.  

Usage:
------

Requires GENESIS 2.3 in linux/unix. GENESIS must run from the command line.

Extract zipped folder and navigate to within the directory demo_modelDB.
Run fig1_tar.m in MATLAB to simulate the neurons, for generating the target 
data. Then run fig1.m in MATLAB to run the fitting routines and output jpg 
plots of the fits. The figure is based on fig1 in the paper (with a longer 
time-step to hasten running speed and minimize RAM use).

The properties of the simulated neurons and fit routine parameters can be 
changed by editing these two files.

Note that most of the figures in the paper were generated from compiled 
versions of this software on a small cluster.

Problems:
---------

Tested on three different installations (VMware Ubuntu, and two clusters 
with Scientific Linux).

If you have problems:
1) Can you run GENESIS from the command prompt in the extracted directory 
demo_modelDB?
If not you need the appropriate path (see GENESIS instructions).
2) Do you have enough RAM?
The routines are quite RAM intensive. Try commenting out some of the 
simulations and editing the code to use a larger time-step.



Lepora NF, Overton PG,Gurney K (2012) Efficient fitting of conductance-based model neurons from somatic current clamp. J. Comp. Neuro. 32(1):1-24[PubMed]

References and models cited by this paper

References and models that cite this paper

Achard P, De Schutter E (2006) Complex parameter landscape for a complex neuron model. PLoS Comput Biol 2:e94 [PubMed]

Booth V, Rinzel J (1995) A minimal, compartmental model for a dendritic origin of bistability of motoneuron firing patterns. J Comput Neurosci 2:299-312 [Journal] [PubMed]

Bower JM, Beeman D (1998) The Book Of Genesis: Exploring Realistic Neural Models With The General Neural Simulation System

Bush K, Knight J, Anderson C (2005) Optimizing conductance parameters of cortical neural models via electrotonic partitions. Neural Netw 18:488-96 [PubMed]

Connors BW, Gutnick MJ (1990) Intrinsic firing patterns of diverse neocortical neurons. Trends Neurosci 13:99-104 [PubMed]

Davis GW (2006) Homeostatic control of neural activity: from phenomenology to molecular design. Annu Rev Neurosci 29:307-23 [PubMed]

Druckmann S, Banitt Y, Gidon A, Schurmann F, Markram H, Segev I (2007) A novel multiple objective optimization framework for constraining conductance-based neuron models by experimental data. Front Neurosci 1:7-18 [PubMed]

Druckmann S, Berger TK, Hill S, Schurmann F, Markram H, Segev I (2008) Evaluating automated parameter constraining procedures of neuron models by experimental and surrogate data. Biol Cybern 99:371-9 [PubMed]

Edelman GM, Gally JA (2001) Degeneracy and complexity in biological systems. Proc Natl Acad Sci U S A 98:13763-8 [PubMed]

Gerstner W, Naud R (2009) Neuroscience. How good are neuron models? Science 326:379-80 [PubMed]

Goldman MS, Golowasch J, Marder E, Abbott LF (2001) Global structure, robustness, and modulation of neuronal models. J Neurosci 21:5229-38 [Journal] [PubMed]

   Global structure, robustness, and modulation of neuronal models (Goldman et al. 2001) [Model]

Golowasch J, Goldman MS, Abbott LF, Marder E (2002) Failure of averaging in the construction of a conductance-based neuron model. J Neurophysiol 87:1129-31 [Journal] [PubMed]

Gunay C, Edgerton JR, Jaeger D (2008) Channel density distributions explain spiking variability in the globus pallidus: a combined physiology and computer simulation database approach. J Neurosci 28:7476-91 [Journal] [PubMed]

   [7 reconstructed morphologies on NeuroMorpho.Org]
   Globus pallidus multi-compartmental model neuron with realistic morphology (Gunay et al. 2008) [Model]

Herz AV, Gollisch T, Machens CK, Jaeger D (2006) Modeling single-neuron dynamics and computations: a balance of detail and abstraction. Science 314:80-5 [PubMed]

Hines M (1984) Efficient computation of branched nerve equations. Int J Biomed Comput 15:69-76 [PubMed]

Hobbs KH, Hooper SL (2008) Using complicated, wide dynamic range driving to develop models of single neurons in single recording sessions. J Neurophysiol 99:1871-83 [PubMed]

Hodgkin AL, Huxley AF (1952) A quantitative description of membrane current and its application to conduction and excitation in nerve. J Physiol 117:500-44 [Journal] [PubMed]

   Squid axon (Hodgkin, Huxley 1952) (LabAXON) [Model]
   Squid axon (Hodgkin, Huxley 1952) (NEURON) [Model]
   Squid axon (Hodgkin, Huxley 1952) (SNNAP) [Model]
   Squid axon (Hodgkin, Huxley 1952) used in (Chen et al 2010) (R language) [Model]
   Squid axon (Hodgkin, Huxley 1952) (SBML, XPP, other) [Model]

Huys QJ, Ahrens MB, Paninski L (2006) Efficient estimation of detailed single-neuron models. J Neurophysiol 96:872-90 [Journal] [PubMed]

   Efficient estimation of detailed single-neuron models (Huys et al. 2006) [Model]

Huys QJ, Paninski L (2009) Smoothing of, and parameter estimation from, noisy biophysical recordings. PLoS Comput Biol 5:e1000379-80 [PubMed]

Izhikevich EM (2003) Simple model of spiking neurons. IEEE Trans Neural Netw 14:1569-72 [Journal] [PubMed]

   Artificial neuron model (Izhikevich 2003, 2004, 2007) [Model]

Keren N, Peled N, Korngreen A (2005) Constraining compartmental models using multiple voltage recordings and genetic algorithms. J Neurophysiol 94:3730-42 [Journal] [PubMed]

Kita T, Kita H, Kitai ST (1984) Passive electrical membrane properties of rat neostriatal neurons in an in vitro slice preparation. Brain Res 300:129-39 [PubMed]

Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, Devon K, Dewar K, Doyle M, Fit (2001) Initial sequencing and analysis of the human genome. Nature 409:860-921

Lemasson G, Maex R (2001) Introduction to equation solving and parameter fitting. Computational neuroscience: Realistic modelling for experimentalists 125:125

Lepora NF, Overton P, Gurney K (2009) Efficient current-based optimization techniques for parameter estimation in multi-compartment neuronal models Bmc Neuroscience 10(S1):347

Lindau M, Neher E (1988) Patch-clamp techniques for time-resolved capacitance measurements in single cells. Pflugers Arch 411:137-46 [PubMed]

Llano I, Marty A, Armstrong CM, Konnerth A (1991) Synaptic- and agonist-induced excitatory currents of Purkinje cells in rat cerebellar slices. J Physiol 434:183-213 [PubMed]

London M, Hausser M (2005) Dendritic computation. Annu Rev Neurosci 28:503-32 [PubMed]

Mahon S, Deniau JM, Charpier S, Delord B (2000) Role of a striatal slowly inactivating potassium current in short-term facilitation of corticostriatal inputs: a computer simulation study. Learn Mem 7:357-62 [Journal] [PubMed]

   Striatal Output Neuron (Mahon, Deniau, Charpier, Delord 2000) [Model]

Mainen ZF, Sejnowski TJ (1996) Influence of dendritic structure on firing pattern in model neocortical neurons. Nature 382:363-6 [Journal] [PubMed]

   [2 reconstructed morphologies on NeuroMorpho.Org]
   Pyramidal Neuron Deep, Superficial; Aspiny, Stellate (Mainen and Sejnowski 1996) [Model]

Marder E, Goaillard JM (2006) Variability, compensation and homeostasis in neuron and network function. Nat Rev Neurosci 7:563-74 [PubMed]

Marder E, Prinz AA (2002) Modeling stability in neuron and network function: the role of activity in homeostasis. Bioessays 24:1145-54 [PubMed]

Markram H (2006) The blue brain project. Nat Rev Neurosci 7:153-60 [Journal] [PubMed]

   [241 reconstructed morphologies on NeuroMorpho.Org]

Mennerick S, Que J, Benz A, Zorumski CF (1995) Passive and synaptic properties of hippocampal neurons grown in microcultures and in mass cultures. J Neurophysiol 73:320-32 [Journal] [PubMed]

Moore EH (1920) On the reciprocal of the general algebraic matrix Bulletin Of The American Mathematical Society 26:394-395

Morse TM, Davison AP, Hines ML (2001) Parameter space reduction in neuron model optimization through minimization of residual voltage clamp current Program No. 606.5. 2001 Abstract Viewer-Itinerary Planner

Penrose R (1955) A generalized inverse for matrices Proc Cambridge Philosophy Soci 51:406-413

Pinsky PF, Rinzel J (1994) Intrinsic and network rhythmogenesis in a reduced Traub model for CA3 neurons. J Comput Neurosci 1:39-60 [Journal] [PubMed]

   CA3 pyramidal cell: rhythmogenesis in a reduced Traub model (Pinsky, Rinzel 1994) [Model]

Pospischil M, Toledo-Rodriguez M, Monier C, Piwkowska Z, Bal T, Fregnac Y, Markram, Destexhe (2008) Minimal Hodgkin-Huxley type models for different classes of cortical and thalamic neurons. Biol Cybern 99:427-41 [Journal] [PubMed]

   Hodgkin-Huxley models of different classes of cortical neurons (Pospischil et al. 2008) [Model]

Prinz AA, Billimoria CP, Marder E (2003) Alternative to hand-tuning conductance-based models: construction and analysis of databases of model neurons. J Neurophysiol 90:3998-4015 [Journal] [PubMed]

RALL W (1960) Membrane potential transients and membrane time constant of motoneurons. Exp Neurol 2:503-32 [PubMed]

Rall W (1969) Time constants and electrotonic length of membrane cylinders and neurons. Biophys J 9:1483-508 [PubMed]

Schulz DJ, Goaillard JM, Marder E (2006) Variable channel expression in identified single and electrically coupled neurons in different animals. Nat Neurosci 9:356-62 [PubMed]

Swensen AM, Bean BP (2005) Robustness of burst firing in dissociated purkinje neurons with acute or long-term reductions in sodium conductance. J Neurosci 25:3509-20 [PubMed]

Tabak J, Murphey CR, Moore LE (2001) Parameter estimation methods for single neuron models. J Comput Neurosci 9:215-36 [Journal] [PubMed]

Taylor AL, Goaillard JM, Marder E (2009) How multiple conductances determine electrophysiological properties in a multicompartment model. J Neurosci 29:5573-86 [PubMed]

Taylor AL, Hickey TJ, Prinz AA, Marder E (2006) Structure and visualization of high-dimensional conductance spaces. J Neurophysiol 96:891-905 [PubMed]

Tien JH, Guckenheimer J (2008) Parameter estimation for bursting neural models. J Comput Neurosci 24:358-73

Tobin AE, Calabrese RL (2006) Endogenous and half-center bursting in morphologically inspired models of leech heart interneurons. J Neurophysiol 96:2089-106 [Journal] [PubMed]

Turrigiano GG, Nelson SB (2004) Homeostatic plasticity in the developing nervous system. Nat Rev Neurosci 5:97-107 [PubMed]

Van Geit W, Achard P, De Schutter E (2007) Neurofitter: a parameter tuning package for a wide range of electrophysiological neuron models. Front Neuroinformatics 1:1-42 [PubMed]

Van Geit W, De Schutter E, Achard P (2008) Automated neuron model optimization techniques: a review. Biol Cybern 99:241-51

Vanier MC, Bower JM (1999) A comparative survey of automated parameter-search methods for compartmental neural models. J Comput Neurosci 7:149-71 [Journal] [PubMed]

Venter JC, Adams MD, Myers EW, Li PW, Mural RJ, Sutton GG, Smith HO, Yandell M, Evans CA, Hol (2001) The sequence of the human genome. Science 291:1304-51 [PubMed]

Weaver CM, Wearne SL (2006) The role of action potential shape and parameter constraints in optimization of compartment models Neurocomputing 69:1053-1057 [Journal]

   AP shape and parameter constraints in optimization of compartment models (Weaver and Wearne 2006) [Model]

Weaver CM, Wearne SL (2008) Neuronal Firing Sensitivity to Morphologic and Active Membrane Parameters. PLoS Comput Biol 4:e11 [Journal] [PubMed]

Wood R, Gurney KN (2002) A new parametric search technique for biophysical models. a study in the rat neostriatal projection neuron 32nd Annual Meeting of the Society for Neuroscience Abstracts

Wood R, Gurney KN, Wilson CJ (2004) A novel parameter optimization technique for compartmental models applied to a model of a striatal mediumspiny neuron Neurocomputing 60:1109-1116

(58 refs)