A two-layer biophysical olfactory bulb model of cholinergic neuromodulation (Li and Cleland 2013)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:149739
This is a two-layer biophysical olfactory bulb (OB) network model to study cholinergic neuromodulation. Simulations show that nicotinic receptor activation sharpens mitral cell receptive field, while muscarinic receptor activation enhances network synchrony and gamma oscillations. This general model suggests that the roles of nicotinic and muscarinic receptors in OB are both distinct and complementary to one another, together regulating the effects of ascending cholinergic inputs on olfactory bulb transformations.
Reference:
1 . Li G, Cleland TA (2013) A two-layer biophysical model of cholinergic neuromodulation in olfactory bulb. J Neurosci 33:3037-58 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism:
Cell Type(s): Olfactory bulb main mitral GLU cell; Olfactory bulb main interneuron periglomerular GABA cell; Olfactory bulb main interneuron granule MC GABA cell;
Channel(s): I Na,p; I L high threshold; I T low threshold; I A; I M; I h; I K,Ca; I CAN; I Sodium; I Calcium; I Potassium; I_Ks; I Cl, leak; I Ca,p;
Gap Junctions:
Receptor(s): Nicotinic; GabaA; Muscarinic; AMPA; NMDA;
Gene(s):
Transmitter(s): Acetylcholine;
Simulation Environment: NEURON; MATLAB;
Model Concept(s): Sensory processing; Sensory coding; Neuromodulation; Olfaction;
Implementer(s): Li, Guoshi [guoshi_li at med.unc.edu];
Search NeuronDB for information about:  Olfactory bulb main mitral GLU cell; Olfactory bulb main interneuron periglomerular GABA cell; Olfactory bulb main interneuron granule MC GABA cell; Nicotinic; GabaA; Muscarinic; AMPA; NMDA; I Na,p; I L high threshold; I T low threshold; I A; I M; I h; I K,Ca; I CAN; I Sodium; I Calcium; I Potassium; I_Ks; I Cl, leak; I Ca,p; Acetylcholine;
/
ACh_ModelDB
data
Input
Readme.txt
cadecay.mod *
cadecay2.mod *
Caint.mod *
Can.mod *
CaPN.mod *
CaT.mod *
GradeAMPA.mod *
GradeGABA.mod *
GradNMDA.mod *
hpg.mod *
kAmt.mod *
KCa.mod *
KDRmt.mod *
kfasttab.mod *
kM.mod *
KS.mod *
kslowtab.mod *
LCa.mod *
nafast.mod *
NaP.mod *
Naxn.mod *
Nicotin.mod *
nmdanet.mod *
OdorInput.mod *
Background.hoc
Connect.hoc
GC_def.hoc
GC_save.hoc *
GC_Stim.hoc
Input.hoc
MC_def.hoc
MC_save.hoc
MC_Stim.hoc
mod_func.c
mosinit.hoc
nrniv.exe.stackdump
OB.hoc
Parameter.hoc
PG_def.hoc
PG_save.hoc *
PG_Stim.hoc
SaveData.hoc
tabchannels.dat *
tabchannels.hoc
                            
NEURON CODE FOR IMPLEMENTING SIMULATIONS OF CHOLINERGIC
NEUROMODULATION IN OLFACTORY BULB

				Guoshi Li
				Computational Physiology Lab
				Department of Psychology
				Cornell University
				Ithaca, NY 14853

A full description of the model may be found in:

Li G and Cleland TA(2013) A two-layer biophysical model of cholinergic
neuromodulation in olfactory bulb.  Journal of Neuroscience
33:3037-3058.

For questions, please email: gl275@cornell.edu

The OB network model contains 25 mitral cells (MCs), 25 periglomerular
cells (PGs) and 100 granule cell (GCs)

The package contains five folders:

"celldata" folder to store data from single cell simulations
"data" folder to store data from network simulations
"SP" folder to store the timing of random background spikes to MCs
"Connection" folder to store the connectivity information between MCs
     and GCs
"Input" folder contains two pre-generated data files read by the
     program to set the pre-odor and steady-state odor values

For network simulation, run mosinit.hoc
For MC single-cell simulation, run MC_Stim.hoc
For GC single-cell simulation, run GC_Stim.hoc
For PG single-cell simulation, run PG_Stim.hoc

The simulation step is set to 0.001 ms for results presented in the
Journal of Neuroscience Paper (Li and Cleland 2013)
Simulation time of the full model is about 3.5 hours for one run in a
workstation
Small time step is used to ensure accuracy of results
To reproduce identical results shown in the paper, make sure a time
step of 0.001 ms is used
Larger step could be used for testing purpose

Some key parameters:
NTCE = 0: Full model with EPL simulation
NTCE = 1: Glomerular model without EPL (i.e., no granule cells)
NICOTIN  = 0: nAChRs inactive
NICOTIN  = 1: nAChRs activated
MUSCARIN = 0: mAChRs inactive
MUSCARIN = 1: mAChRs activated

The default setting is the full model under the Control case (nAChRs
and mAChRs inactive)

The data saved in the "data" folder after simulation is analyzed using
the following custom Matalb scripts

PlotV.m: plot cell voltages
PlotG.m: plot GABAa conductances
Raster.m: generate raster plots of spikes
LFP.m: frequency analysis of the sLFP
Phase.m: generate phase distribution plot and raster plot of spike phases

Changelog
---------
2022-05: Updated MOD files to compile with the latest neuron releases where
         ion variables used as STATE can not be declared as GLOBAL.