Olfactory Bulb mitral-granule network generates beta oscillations (Osinski & Kay 2016)

 Download zip file 
Help downloading and running models
This model of the dendrodendritic mitral-granule synaptic network generates gamma and beta oscillations as a function of the granule cell excitability, which is represented by the granule cell resting membrane potential.
1 . Osinski BL, Kay LM (2016) Granule cell excitability regulates gamma and beta oscillations in a model of the olfactory bulb dendrodendritic microcircuit. J Neurophysiol 116:522-39 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network; Dendrite;
Brain Region(s)/Organism: Olfactory bulb;
Cell Type(s): Olfactory bulb main interneuron granule MC GABA cell; Olfactory bulb main interneuron granule TC GABA cell; Abstract integrate-and-fire leaky neuron;
Channel(s): I N; I Sodium; I Calcium;
Gap Junctions:
Receptor(s): AMPA; NMDA; Gaba;
Transmitter(s): Gaba; Glutamate;
Simulation Environment: MATLAB;
Model Concept(s): Oscillations; Active Dendrites; Extracellular Fields; Calcium dynamics; Gamma oscillations; Beta oscillations; Olfaction;
Implementer(s): Osinski, Boleslaw [boleszek at uchicago.edu];
Search NeuronDB for information about:  Olfactory bulb main interneuron granule MC GABA cell; Olfactory bulb main interneuron granule TC GABA cell; AMPA; NMDA; Gaba; I N; I Sodium; I Calcium; Gaba; Glutamate;
% Code for running ParamSeep_GCE on the cronusx server

% Boleslaw Osinski (2015)


clear all
input_file = 'OB_params_GCE.txt';
[dt,tsim,numtp,nmit,ngradist,ngraprox,sampf,timevec] ...
    = InitNetwork_GCE(input_file);

%% Run ParamSweep_GCE

% set P1
% P1.line = 44;
% P1.name = 'wGABAGR ';
% P1.val = 0.013:0.001:0.017;

P1.line = 7;
P1.name = 'Wmin ';
P1.val = 0.0128;

% set P2
P2.line = 57;
P2.name = 'Vrest ';
P2.val = 1e-3 * [-75];

% Simulate the model over chosen parameter range
[MCMAT dGCMAT ICMAT param MClfpMAT dGClfpMAT] = ...

%% Calculate LFP frequency and power

FmaxMAT = zeros(length(P1.val),length(P2.val));
maxpwrMAT = zeros(length(P1.val),length(P2.val));

trim = 500; % Ignore first 500 timepoints

L = length(timevec(trim:end-100));  % Length of data that will be analyzed
NFFT = 2^nextpow2(L); % Next power of 2 from L
f = sampf/2*linspace(0,1,NFFT/2+1);
ROI = ceil(8/(f(2)-f(1))):ceil(140/(f(2)-f(1)));
% ROI gives us a window from ~7Hz - 115Hz from which to select peak power.
% This avoids selecting low freuqency components that may dominate the
% spectrum.

% Run FFT on each simulated LFP
for n1 = 1:length(P1.val)
    for n2 = 1:length(P2.val)
        MitLFP = MClfpMAT{n1,n2}.GradistMitGlobal;
        mitFFT = fft(detrend(MitLFP(trim:end-100,1),'constant'),NFFT)/L;

        absmitFFT = 2*abs(mitFFT(1:NFFT/2+1));
        maxpwrMAT(n1,n2) = max(absmitFFT(ROI));

        maxind = find(absmitFFT == maxpwrMAT(n1,n2));
        FmaxMAT(n1,n2) = f(maxind);