Phase precession through acceleration of local theta rhythm (Castro & Aguiar 2011)

 Download zip file 
Help downloading and running models
Accession:143248
"... Here we present a biophysical spiking model for phase precession in hippocampal CA1 which focuses on the interaction between place cells and local inhibitory interneurons. The model’s functional block is composed of a place cell (PC) connected with a local inhibitory cell (IC) which is modulated by the population theta rhythm. Both cells receive excitatory inputs from the entorhinal cortex (EC). ..."
Reference:
1 . Castro L, Aguiar P (2012) Phase precession through acceleration of local theta rhythm: a biophysical model for the interaction between place cells and local inhibitory neurons. J Comput Neurosci 33:141-50 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism: Hippocampus;
Cell Type(s):
Channel(s):
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: MATLAB;
Model Concept(s): Oscillations; Synchronization; Place cell/field;
Implementer(s):
function new_P = unwindphases( T, P, t0, p0 )
%UNWINDPHASES       unwind a set of cyclic phases
%
%   Description:
%   T is a vector of times
%   P is a vector of phases
%   t0 and p0 are used to set the border line between points ribbons:
%       t0 is the point where the border crosses the phase axis
%       p0 is the point where the border crosses the time axis
%
%  2pi ---------------------...
%     |  ..     \    ..
%   p0|\  :.      \    :.
%     |  \  ::.     \    .::
%     |    \  .::.    \     .:
%     |   .  \  ..      \
%     |    :.  \  :..     \
%   0  ----------------------...
%     0         t0
%
%   Author: Paulo Aguiar, pauloaguiar@fc.up.pt

m = - p0/t0;   % slope

for i = 1:length(T)
    
    % set border equation, of type y = m*x + b
    k = fix( T(i)/t0 );   % border number
    b = (k * 2*pi) + p0;  % border line parameter
    offset = 2.0*pi * ( 1 - k );
    if P(i) < m * T(i) + b      % below border
        new_P(i) = P(i) + offset;
    else                        % above border
        new_P(i) = P(i) + offset - 2.0*pi;
    end
    
end

Loading data, please wait...