//New Parametri File - to suit distribution of ion channels in axon/soma/hillock //Last modified 14-March-2008 //Shyam Diwakar M. //Universita' degli Studi di Milano //Updated Na channel percentage to 50% in axon-hillock to maintain propagation velocity //objref panel,morfobox,axonbox,dendbox //A Panel for Channels and someof their controls //panel = new VBox() objref panel panel = new VBox() //new stuff 10 May 2005 Nag = Granule[0].soma.gnabar_GRC_NA Kvg = Granule[0].soma.gkbar_GRC_KV Kmg = Granule[0].soma.gkbar_GRC_KM glL = 5.68e-5//Granule[0].soma.gl_GRC_LKG1 ell = -16.5//Granule[0].soma.gl_GRC_LKG1 el3 = -70 ndend = 4 nsg = 5 naxon = 56 npf=500 ncomp = 1+(4*ndend)+nsg+naxon+npf //for PF NagPf = Nag*10 KvgPf = Kvg*10 RappPf=0.1//((9.76*9.76)/(npf*Granule[0].pf[0].L*Granule[0].pf[0].diam)) Rappaxon = ((9.76*9.76)/(naxon*Granule[0].axon[0].L*Granule[0].axon[0].diam)) Granule[0].soma.gnabar_GRC_NA = 0 Granule[0].soma.gkbar_GRC_KV = 0 KirGmax=0.0009//Granule[0].soma.gkbar_GRC_KIR //Standaradisation reference value - Kir KaGmax=0.0032//Granule[0].soma.gkbar_GRC_KA //Standaradisation reference value -Ka CaGmax=0.00046 //Standaradisation reference value - Ca KCaGmax=0.003//Granule[0].dend_4[3].gkbar_GRC_KCA //Standaradisation reference value - KCa beta=0.6//Granule[0].dend_4[3].beta_GRC_CALC //Standaradisation reference value - removal rate inicon=0.001//Granule[0].dend_4[3].cai0_GRC_CALC //Standaradisation reference value - Initial condition //Mode flags for state/compartment activation Camode1 =0 Camode2 =0 Camode3 =0 Camode4 =1 CamodeS =0 KCamode1 =0 KCamode2 =0 KCamode3 =0 KCamode4 =1 KCamodeS=0 Kirmode1 =1 Kirmode2 =1 Kirmode3 =0 Kirmode4 =0 KirmodeS =0 Kamode1 =0 Kamode2 =0 Kamode3 =0 Kamode4 =0 KamodeS =1 SomaArea=Granule[0].soma.L*Granule[0].soma.diam*PI Dend12Area=Granule[0].dend_1[0].L*Granule[0].dend_1[0].diam*PI Dend34Area=Granule[0].dend_4[0].L*Granule[0].dend_4[0].diam*PI SomascArea=PI*9.76*9.76 RappSomaDend12=SomascArea/(4*Dend12Area) RappSomaDend34=SomascArea/(4*Dend34Area) RappSomaNew=SomascArea/SomaArea RappSomahill=SomascArea/(3.75*PI) RappAH = 3.75/(naxon*Granule[0].axon[0].L*Granule[0].axon[0].diam) gG = Granule[0].soma.ggaba_GRC_LKG2 proc alpKCaM() { alpKCa = ($1==1)+($2==1)+($3==1)+($4==1)+($5==1) } proc alpCaM() { alpCa = ($1==1)+($2==1)+($3==1)+($4==1)+($5==1) } proc alpKaM() { alpKa = ($1==1)+($2==1)+($3==1)+($4==1)+($5==1) } proc alpKirM() { alpKir = ($1==1)+($2==1)+($3==1)+($4==1)+($5==1) } gamma = 0.5 //Percentage of NA/Kv in axon-hillock //For axon and Hill NagH = Nag KvgH = Kvg NagA = Nag KvgA = Kvg KCaD = KCaGmax CaD = CaGmax DendFact=1 //default morphology scaling ratio KaRapp = KaGmax KirRapp = KirGmax CaRapp = CaGmax KCaRapp = KCaGmax //For Dendritic Morphology Scaling proc DendGeomFact(){ //print "Le dimensioni dei dendriti sono (guardando dal soma):" for (i=0;i<4;i=i+1) { Granule[0].dend_1[i].diam=0.75/DendFact Granule[0].dend_2[i].diam=0.75/DendFact Granule[0].dend_3[i].diam=0.75/DendFact Granule[0].dend_4[i].diam=0.75/DendFact Granule[0].dend_1[i].L=5*DendFact Granule[0].dend_2[i].L=5*DendFact Granule[0].dend_3[i].L=2.5*DendFact Granule[0].dend_4[i].L=2.5*DendFact } //print "L1= ",Granule[0].dend_1[0].L," L2= ",Granule[0].dend_2[0].L," L3= ",Granule[0].dend_3[0].L," L4= ",Granule[0].dend_4[0].L //print "D1=D2=D3=D4= ",Granule[0].dend_1[0].diam } //Updating leakage proc glUpdate() { Granule[0].soma.gl_GRC_LKG1 = glL*(RappSomaNew)*(2/3)//-1.8317)//*(1/ncomp) Granule[0].soma.el_GRC_LKG1 = ell for(i=0;i<5;i=i+1) { Granule[0].hillock[i].gl_GRC_LKG1=glL*(RappSomahill)*(1/15)//-11.613)//*(1/5) Granule[0].hillock[i].el_GRC_LKG1=ell } for(i=0;i=1) { for (i=0;i<4;i=i+1) { Granule[0].soma.gkbar_GRC_KA=KaGmax*(1/alpKa)*RappSomaNew*$1 Granule[0].dend_1[i].gkbar_GRC_KA=KaGmax*RappSomaDend12*(1/alpKa)*$2 Granule[0].dend_2[i].gkbar_GRC_KA=KaGmax*RappSomaDend12*(1/alpKa)*$3 Granule[0].dend_3[i].gkbar_GRC_KA=KaGmax*RappSomaDend34*(1/alpKa)*$4 Granule[0].dend_4[i].gkbar_GRC_KA=KaGmax*RappSomaDend34*(1/alpKa)*$5 } } } proc CaU(){ //print "Refresh Ca" alpCaM($1,$2,$3,$4,$5) if(alpCa>=1) { for (i=0;i<4;i=i+1) { Granule[0].soma.gcabar_GRC_CA=CaGmax*(1/alpKa)*RappSomaNew*$1 Granule[0].dend_1[i].gcabar_GRC_CA=CaD*RappSomaDend12*(1/alpCa)*$2 Granule[0].dend_2[i].gcabar_GRC_CA=CaD*RappSomaDend12*(1/alpCa)*$3 Granule[0].dend_3[i].gcabar_GRC_CA=CaD*RappSomaDend34*(1/alpCa)*$4 Granule[0].dend_4[i].gcabar_GRC_CA=CaD*RappSomaDend34*(1/alpCa)*$5 } } } proc KCaU(){ //print "Refresh KCa" //if($1==1) ->addstuff to modify shell d in soma alpKCaM($1,$2,$3,$4,$5) if(alpKCa>=1) { for (i=0;i<4;i=i+1) { Granule[0].soma.gkbar_GRC_KCA=KCaD*(1/alpKCa)*RappSomaNew*$1 Granule[0].dend_1[i].gkbar_GRC_KCA=KCaD*RappSomaDend12*(1/alpKCa)*$2 Granule[0].dend_2[i].gkbar_GRC_KCA=KCaD*RappSomaDend12*(1/alpKCa)*$3 Granule[0].dend_3[i].gkbar_GRC_KCA=KCaD*RappSomaDend34*(1/alpKCa)*$4 Granule[0].dend_4[i].gkbar_GRC_KCA=KCaD*RappSomaDend34*(1/alpKCa)*$5 } } } proc KirU(){ //print "Refresh Kir" alpKirM($1,$2,$3,$4,$5) if(alpKir>=1) { for (i=0;i<4;i=i+1) { Granule[0].soma.gkbar_GRC_KIR=KirGmax*(1/alpKir)*RappSomaNew*$1 Granule[0].dend_1[i].gkbar_GRC_KIR=KirGmax*RappSomaDend12*(1/alpKir)*$2 Granule[0].dend_2[i].gkbar_GRC_KIR=KirGmax*RappSomaDend12*(1/alpKir)*$3 Granule[0].dend_3[i].gkbar_GRC_KIR=KirGmax*RappSomaDend34*(1/alpKir)*$4 Granule[0].dend_4[i].gkbar_GRC_KIR=KirGmax*RappSomaDend34*(1/alpKir)*$5 } } } //New code for Na in axon/hillock -11may05 proc NaAUpdate() { //print "Updating Na in axon" for(i=0;i