Computational model of cerebellar tDCS (Zhang et al., 2021)

 Download zip file 
Help downloading and running models
Accession:267189
This archive contains models used in (Zhang et al. 2021) and simulates Purkinje cell, granule cell, and deep cerebellar neuron activities under cerebellar tDCS (transcranial direct current stimulation).
Reference:
1 . Zhang X, Hancock R, Santaniello S (2021) Transcranial Direct Current Stimulation of Cerebellum Alters Spiking Precision in Cerebellar Cortex: A Modeling Study of Cellular Responses. PLoS Comp Biol 17(12):e1009609 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell;
Brain Region(s)/Organism: Cerebellum;
Cell Type(s): Cerebellum Purkinje GABA cell; Cerebellum interneuron granule GLU cell; Cerebellum deep nucleus neuron;
Channel(s):
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s):
Implementer(s):
Search NeuronDB for information about:  Cerebellum Purkinje GABA cell; Cerebellum interneuron granule GLU cell;
% Simulate Purkinje cell activities under climbing fiber synaptic inputs
% and voltage shifts applied to soma or dendritic compartments
% Generate the simulation results used in Fig.7
%
% (c) Xu (Shawn) Zhang, UConn
% xu.3.zhang@uconn.edu
% November 2021
%

clear;clc;

load('dvall.mat');
vshift_all = -1.5:0.25:1.5;

%% Somatic voltage shifts
mkdir('recordings_calcium_soma_vshift');
for n = 1:length(vshift_all)
    disp(n);
    vshift = vshift_all(n);
    
    tmpdv = dv_all*vshift; % Voltage shifts for all compartments
    
    formatSpec = '%.4f \n';
    opt_fileID = fopen('Vshifts_calcium_soma_vshift.txt','w');
    fprintf(opt_fileID,formatSpec,tmpdv);
    fclose(opt_fileID);
    
    formatSpec = '%s \n';
    opt_fileID = fopen('trialname_calcium_soma_vshift.txt','w');
    f_soma = strcat('recordings_calcium_soma_vshift/soma','_',num2str(n),'.txt');
    fprintf(opt_fileID,formatSpec,f_soma);
    fclose(opt_fileID);
    
    [status,cmdout] = system('nrniv mosinit_calcium_soma_vshift.hoc');
end

%% Dendritic voltage shifts
mkdir('recordings_calcium_dend_vshift');
for n = 1:length(vshift_all)
    disp(n);
    vshift = vshift_all(n);
    
    tmpdv = dv_all*vshift; % Voltage shifts for all compartments
    
    formatSpec = '%.4f \n';
    opt_fileID = fopen('Vshifts_calcium_dend_vshift.txt','w');
    fprintf(opt_fileID,formatSpec,tmpdv);
    fclose(opt_fileID);
    
    formatSpec = '%s \n';
    opt_fileID = fopen('trialname_calcium_dend_vshift.txt','w');
    f_soma = strcat('recordings_calcium_dend_vshift/soma','_',num2str(n),'.txt');
    fprintf(opt_fileID,formatSpec,f_soma);
    fclose(opt_fileID);
    
    [status,cmdout] = system('nrniv mosinit_calcium_dend_vshift.hoc');
end

Loading data, please wait...