Basal ganglia-corticothalamic (BGCT) network (Chen et al., 2014)

 Download zip file 
Help downloading and running models
Accession:152113
We developed a biophysical model of the basal ganglia-corticothalamic network in this work. "... We demonstrate that the typical absence seizure activities can be controlled and modulated by the direct GABAergic projections from the substantia nigra pars reticulata (SNr) to either the thalamic reticular nucleus (TRN) or the specific relay nuclei (SRN) of thalamus, through different biophysical mechanisms. ... results highlight the bidirectional functional roles of basal ganglia in controlling and modulating absence seizures, and might provide novel insights into the therapeutic treatments of this brain disorder."
References:
1 . Chen M, Guo D, Wang T, Jing W, Xia Y, Xu P, Luo C, Valdes-Sosa PA, Yao D (2014) Bidirectional control of absence seizures by the basal ganglia: a computational evidence. PLoS Comput Biol 10:e1003495 [PubMed]
2 . Chen M, Guo D, Li M, Ma T, Wu S, Ma J, Cui Y, Xia Y, Xu P, Yao D (2015) Critical Roles of the Direct GABAergic Pallido-cortical Pathway in Controlling Absence Seizures. PLoS Comput Biol 11:e1004539 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neural mass;
Brain Region(s)/Organism: Neocortex; Thalamus; Basal ganglia;
Cell Type(s):
Channel(s):
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: XPP; MATLAB;
Model Concept(s): Epilepsy;
Implementer(s): Guo, Daqing [dqguo at uestc.edu.cn];
len1=31;
trails=30; 

delay=0.05;
v_sr=linspace(0.4,1.76,len1)*1.0e-3;
v_p1xi=3.0e-4;
open1=1;
open2=1;
KK=1;

matlabpool local 3


Pmax1=zeros(len1,trails);
Pmax2=zeros(len1,trails);
Pmin1=zeros(len1,trails);
Pmin2=zeros(len1,trails);
FR=zeros(len1,trails);
State=zeros(len1,trails);
FD=zeros(len1,trails);

for i=1:len1
    tic
    vv_sr=v_sr(i); 
    pmax1=[];
    pmax2=[];
    pmin1=[];
    pmin2=[];
    state=[];
    fd=[];
    parfor j=1:trails
        [p1,p2,p3,p4,p5,p6]=BGCT_subfun2(delay,vv_sr,v_p1xi,open1,open2,KK);
        pmax1=[pmax1,p1];
        pmax2=[pmax2,p2];
        pmin1=[pmin1,p3];
        pmin2=[pmin2,p4];
        state=[state,p5];
        fd=[fd,p6];
    end
    
    Pmax1(i,:)=pmax1;
    Pmax2(i,:)=pmax2;
    Pmin1(i,:)=pmin1;
    Pmin2(i,:)=pmin2;
    State(i,:)=state;
    FD(i,:)=fd;
    toc
end

matlabpool close

v_sr=v_sr*1000;
plot(v_sr,log10(Pmax1),'.r',v_sr,log10(Pmax2),'.b',v_sr,log10(Pmin2),'.m',v_sr,log10(Pmin1),'.k')

Loading data, please wait...