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

 Download zip file 
Help downloading and running models
Accession:91387
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).
Reference:
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;
Gene(s):
Transmitter(s):
Simulation Environment: C or C++ program;
Model Concept(s): Oscillations; Delay; Olfaction;
Implementer(s):
Search NeuronDB for information about:  Olfactory bulb main mitral GLU cell; GabaA; I Na,p; I Na,t; I A; I K; 
/**************************************************************************

	ChanStd.h													JJS 9/08/95
	
		part of CONICAL, the Computational Neuroscience Class Library
	
	Objects of the ChanStd class are active channels whose gating
	variables (M and H) are updated by one of three standard equation
	forms.  All three formulas use the same parameters (A, B, and V0)
	but these have very different meanings from form to form.  The
	equations used are defined in StdForms.h.
	
	Requires:
		ChanTAU			-- base class
		StdForms.h		-- standard formula definitions
		
**************************************************************************/

#ifndef CHANTM_H
#define CHANTM_H

#include "ChanTAU.h"
#include "StdForms.h"

// one channel needs four equations, denoted by:
enum StdEqType {
	kStdTauM,
	kStdInfM,
	kStdTauH,
	kStdInfH,
kQtyStdEqTypes };

// each equation has three parameters, denoted by:
enum {
	kStdParaA,
	kStdParaB,
	kStdParaV0,
	kStdParaC,
kQtyStdParas };

class ChanTM : public ChanTAU
{
  public:

	// constructors
	ChanTM( Compartment *pTo, real pMaxG=0.1, real pMexp=1, real pHexp=1 )
	: ChanTAU( pTo, pMaxG, pMexp, pHexp ) {}
	ChanTM( VSink *pTo, VSource *pComp, real pMaxG=0.1, real pMexp=1, real pHexp=1 )
	: ChanTAU( pTo, pComp, pMaxG, pMexp, pHexp ) {}
	
	// function defining method (a setter)
	virtual void SetFunc( StdEqType pType, StdForm pForm, real pA, real pB, real pV0, real pC )
	{	itsForm[pType] = pForm;
		itsParam[pType][kStdParaA] = pA;
		itsParam[pType][kStdParaB] = pB;
		itsParam[pType][kStdParaV0]= pV0;
		itsParam[pType][kStdParaC]= pC;
	}
		
	// channel activation/inactivation functions
	virtual real TauM( const real V ) const
	{	return StdFormula( itsForm[kStdTauM], V,
			itsParam[kStdTauM][kStdParaA], 
			itsParam[kStdTauM][kStdParaB],
			itsParam[kStdTauM][kStdParaV0],
            itsParam[kStdTauM][kStdParaC] ); }
		
	virtual real InfM( const real V) const
	{	return StdFormula( itsForm[kStdInfM], V,
			itsParam[kStdInfM][kStdParaA], 
			itsParam[kStdInfM][kStdParaB],
			itsParam[kStdInfM][kStdParaV0],
            itsParam[kStdInfM][kStdParaC] ); }

	virtual real TauH( const real V ) const
	{	return StdFormula( itsForm[kStdTauH], V,
			itsParam[kStdTauH][kStdParaA], 
			itsParam[kStdTauH][kStdParaB],
			itsParam[kStdTauH][kStdParaV0],
            itsParam[kStdTauH][kStdParaC] ); }

	virtual real InfH( const real V ) const
	{	return StdFormula( itsForm[kStdInfH], V,
			itsParam[kStdInfH][kStdParaA], 
			itsParam[kStdInfH][kStdParaB],
			itsParam[kStdInfH][kStdParaV0],
            itsParam[kStdInfH][kStdParaC] ); }

	// public variables which define the functions:

	StdForm	itsForm[kQtyStdEqTypes];	// format (kExpForm, kSigForm, kLinForm)
	real	itsParam[kQtyStdEqTypes][kQtyStdParas];
	
};

#endif