Spatial coupling tunes NMDA receptor responses via Ca2+ diffusion (Iacobucci and Popescu 2019)

 Download zip file 
Help downloading and running models
This code implements a coupled markov model for analysis of positive or negative ion channel coupling from measured unitary currents in patch clamp recordings see our paper: Spatial Coupling Tunes NMDA Receptor Responses via Ca2+ Diffusion Gary J. Iacobucci and Gabriela K. Popescu Journal of Neuroscience 6 November 2019, 39 (45) 8831-8844; DOI:
1 . Iacobucci GJ, Popescu GK (2019) Spatial Coupling Tunes NMDA Receptor Responses via Ca2+ Diffusion. J Neurosci 39:8831-8844 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Synapse; Extracellular;
Brain Region(s)/Organism: Hippocampus;
Cell Type(s):
Channel(s): I Calcium;
Gap Junctions:
Receptor(s): NMDA;
Transmitter(s): Glutamate;
Simulation Environment: MATLAB;
Model Concept(s): Synaptic Plasticity; Development;
Search NeuronDB for information about:  NMDA; I Calcium; Glutamate;
function [Results] = ChannelCouplingAnalysis(coupling_type, method, deadtime)

%Copyright (c) 2016 by Gary Iacobucci
%coupling_type - specify whether you are analyzing channels that are
%positively or negatively coupled. Enter 1 for positive coupling and 2 for
%negative coupling

%method - specify the method of analyzing the transition probability matrix
%(eg. by sequence of state transitions (where all diagonal elements = 0) or
%by time-series (state at each sampled time point)). Enter 1 for time-series
%method or enter 2 for state-sequence method.

%deadtime - enter the deadtime in the form of sample point number. must be
%a whole integer. a value of 1 implies no imposed deadtime. 

%Results - will contain the results of the fitted parameters from the
%'channelcoupling.m' script for each file. Fitted parameters are based on
%the notation from Chung and Kennedy, 1996 and Moreno et al. 2016 eLife
% output.k - coupling coefficient
% output.r - open-to-open probability
% output.z - closed-to-closed probability

%%%%%  HOW TO USE
%prior to running this script, you will need to create an array named 'fileID' containing
%the file names of the files you will analyze. This program is written to
%accept files in the DWT format from QUB software that expresses state
%index and corresponding dwell time in column format. 
%example array format:
% fileID = {'GI-20160810a.dwt';
% 'file1.dwt';
% 'file2.dwt';
% 'file3.dwt';
% 'file4.dwt';
% 'file5.dwt'};
%Example command line input to run script:
% [Results] = ChannelCouplingAnalysis(2, 1, 1);

Results = {};
h = waitbar(0,'Please wait...');
for n = 1:1:numel(fileID)
    filename = fileID{n,1};
    Results{n,1} = output;

Loading data, please wait...