Olfactory bulb network model of gamma oscillations (Bathellier et al. 2006; Lagier et al. 2007)

 Download zip file 
Help downloading and running models
This model implements a network of 100 mitral cells connected with asynchronous inhibitory "synapses" that is meant to reproduce the GABAergic transmission of ensembles of connected granule cells. For appropriate parameters of this special synapse the model generates gamma oscillations with properties very similar to what is observed in olfactory bulb slices (See Bathellier et al. 2006, Lagier et al. 2007). Mitral cells are modeled as single compartment neurons with a small number of different voltage gated channels. Parameters were tuned to reproduce the fast subthreshold oscillation of the membrane potential observed experimentally (see Desmaisons et al. 1999).
1 . Bathellier B, Lagier S, Faure P, Lledo PM (2006) Circuit properties generating gamma oscillations in a network model of the olfactory bulb. J Neurophysiol 95:2678-91 [PubMed]
2 . Lagier S, Panzanelli P, Russo RE, Nissant A, Bathellier B, Sassoè-Pognetto M, Fritschy JM, Lledo PM (2007) GABAergic inhibition at dendrodendritic synapses tunes gamma oscillations in the olfactory bulb. Proc Natl Acad Sci U S A 104:7259-64 [PubMed]
3 . Bathellier B, Lagier S, Faure P, Lledo PM (2006) Corrigendum for Bathellier et al., J Neurophysiol 95 (4) 2678-2691. J Neurophysiol 95:3961-3962
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism: Olfactory bulb;
Cell Type(s): Olfactory bulb main mitral GLU cell;
Channel(s): I Na,p; I Na,t; I A; I K;
Gap Junctions:
Receptor(s): GabaA;
Simulation Environment: C or C++ program;
Model Concept(s): Oscillations; Delay; Olfaction;
Search NeuronDB for information about:  Olfactory bulb main mitral GLU cell; GabaA; I Na,p; I Na,t; I A; I K;

	ChannelRk4.h													JJS 8/29/95
		derived from CONICAL, the Computational Neuroscience Class Library
	A Channel is a type of Current whose voltage V is fixed, but whose
	conductance G varies with time.  A Channel implements an active
	channel; for a passive channel, you can simply use a Current.

	Note that since a Channel needs a voltage to activate the channel,
	it must be either attached to a Compartment, or separately attached
	to a VSink and a VSource (see constructors below).  The former is
	appropriate for a membrane channel; the latter would be appropriate
	for (say) a synapse, whose conductance depends on presynaptic
	voltage, but whose current affects the postsynaptic compartment.

		Current			-- base class
		Stepper			-- base class
		Compartment		-- source of V for channel activation

#ifndef CHANNEL_H
#define CHANNEL_H

#include "CurrentRk4.h"
#include "StepperRk4.h"
#include "CmprtmntRk4.h"

class Channel : public Current, virtual public Stepper
	Channel( Compartment *pTo, real pMaxG=0.1 )					// constructor
	: Current( pTo ), itsComp( pTo ), MaxG( pMaxG ) {}

	Channel( VSink *pTo, VSource *pComp, real pMaxG=0.1 )		// constructor
	: Current( pTo ), itsComp( pComp ), MaxG( pMaxG ) {}
	virtual void Step( const real dt )		// update G
	{ G = 0; }				// (default behavior is passive -- override here)
    virtual void Stepk1( const real dt )		// update Gk1
	{ Gk1 = 0; }	
    virtual void Stepk2( const real dt )		// update Gk2
	{ Gk2 = 0; }	
    virtual void Stepk3( const real dt )		// update Gk3
	{ Gk3 = 0; }	
	virtual void Stepk4( const real dt )		// update Gk4
	{ Gk4 = 0; }
    virtual void Init() {}	
	// public variables:
	real MaxG;					// maximum conductance

	VSource *itsComp;			// compartment whose V affects our G