doi:10.1016/S0361-9230(99)00256-7
Copyright © 2000 Elsevier Science Inc. All rights reserved.
Articles
A reduced compartmental model of the mitral cell for use in network models of the olfactory bulb
Andrew P. Davison
,
, a, Jianfeng Fenga and David Browna
a Laboratory of Computational Neuroscience, The Babraham Institute, Babraham, Cambridge, UK
Received 6 August 1999;
revised 3 November 1999.
Available online 10 March 2000.
References and further reading may be available for this article. To view references and further reading you must
purchase this article.
Abstract
We have developed two-, three- and four-compartment models of a mammalian olfactory bulb mitral cell as a reduction of a complex 286-compartment model [1]. A minimum of three compartments, representing soma, secondary (basal) dendrites and the glomerular tuft of the primary dendrite, is required to adequately reproduce the behaviour of the full model over a broad range of firing rates. Adding a fourth compartment to represent the shaft of the primary dendrite gives a substantial improvement. The reduced models exhibit behaviours in common with the full model which were not used in fitting the model parameters. The reduced models run 75 or more times faster than the full model, making their use in large, realistic network models of the olfactory bulb practical.
Author Keywords: Olfaction; Single-neuron models; Simplified models; Biological neural networks
FIG. 1. Response of the full model to a total constant current of 0.4 nA to the glomerulus, divided evenly among the compartments (
4.3 pA per compartment). (A) Somatic membrane potential trace and (B) Membrane potential trace at base of glomerular tuft, for varying K in the tuft. These figures show the double somatic spike and extended dendritic calcium spike seen with low K, that are removed by increasing the value of the parameter. Increasing Kfurther delays the somatic spike and abolishes the dendritic calcium spike.
FIG. 2. Schematic of the four-compartment model showing ion channels and points of stimulation. Currents are sodium (INa), slow potassium delayed rectifier (IK), fast potassium delayed rectifier (IKfast), potassium anomalous rectifier (IKA), calcium-dependent potassium (IKCa) and L-type calcium (ILCa). The conductances between soma and secondary dendrite, soma and primary dendrite, and primary dendrite and glomerulus are gsd, gsp and gpg, respectively.
FIG. 3. Somatic membrane potential trace for full and reduced models with Ifull = 0.4 μA cm−2 after fitting to spike shape. Both three- and four-compartment models give good fits to the full model. The two-compartment model has a markedly raised reset potential.
FIG. 4. Somatic membrane potential trace for full and reduced models with Ifull = 0.8 μA cm−2 after fitting to spike shape with Ifull = 0.4 μA cm−2. The miniature spike on the down-slope of the action potential in the three-compartment model is due to the delay between spiking in the soma and in the secondary dendrite.
FIG. 5. Comparing fit-to-shape to fit-to-time for the four-compartment model. (A) Firing frequency, (B) latency of first spike. Current was injected in the soma compartment. The four compartment model with fit-to-shape parameters shows almost zero deviation from the full model with current injection 0.4 μA cm−2, the level at which it was fitted, but large differences with higher or lower inputs. In contrast, the four-compartment model with fit-to-time parameters fits reasonably well over a wide range of input levels.
FIG. 6. Somatic membrane potential trace for full and four-compartment models with Ifull = 0.2 μA cm−2 and 1.6 μA cm−2 after fitting to spike timing. (A) Current injection to soma. (B) Current injection to glomerulus. The models fit well in terms of first spike latency and firing frequency, but there are discrepancies in the detailed shape of the membrane potential traces. This illustrates one of the problems with very-reduced compartmental models.
FIG. 7. Comparing reduced models with two, three and four compartments, with fit-to-time parameters: (A) Firing frequency as a function of current for stimulation of soma; (B) as A for stimulation of glomerulus; (C) Latency of first spike as a function of current for stimulation of soma; (D) as C for stimulation of glomerulus. The three- and four-compartment models fit the full model closely for a wide range of input levels and for both somatic and glomerular inputs. The two-compartment model gives a tolerable fit for somatic input and for low glomerular input, but deviates considerably for high glomerular input.
Table 1. Best fit parameter values
