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;

	Stepmaster.h												JJS 2/21/96
		part of CONICAL, the Computational Neuroscience Class Library

	Any Stepper object must be attached to a Stepmaster.  This serves
	two purposes: (1) it keeps track of which memory index (0 or 1) is
	the current one, via its GetCurIdx() method; (2) it can call the
	Step() method of all its attached Steppers, via StepAll(dt).
	When a Stepmaster object is killed, all the attached Steppers are
	destroyed as well.  If a Stepper is destroyed, the Stepmaster is
	automatically notified but remains intact.  A Stepper can only be
	removed by attaching it to another Stepmaster.
	For ease of use in normal situations, a global gStepmaster is
	automatically defined.  By default, any Steppers created attach
	to gStepmaster.



#ifndef real
#define real double

class Stepper;

// declare a linked list node for keeping track of Steppers
class StepperNode
  friend class Stepmaster;
	StepperNode( StepperNode *pPrev, Stepper& pIt, StepperNode *pNext=0 )
	{ itsStepper = &pIt; if (pPrev) pPrev->itsNext = this; itsNext = pNext; }  
	Stepper	*itsStepper;
	StepperNode	*itsNext;

class Stepmaster
	friend class Stepper;
	Stepmaster( );						// constructor
	~Stepmaster( );						// destructor

	// simulation-building methods
	void Attach( Stepper& );			// add a Stepper to our list

	// simulation-running methods
	void StepAll( const real dt );		// update all attached Steppers
    void doStepk1(const real dt);         //update k1 in all attached steppers 
	void doStepk2(const real dt);         //update k2 in all attached steppers
	void doStepk3(const real dt);         //update k3 in all attached steppers
	void doStepk4(const real dt);         //update k4 in all attached steppers
	void doStepinit(const real dt);       //initialise what need to be...
	int GetCurIdx( ) { return itsCurIdx; }

	void Remove( Stepper& );			// remove a Stepper (dying or defecting)

	// protected variables:
	StepperNode *itsListHead, *itsListTail;	
	int			itsCurIdx;