//genesis
/*************************** MS Model, Version 9.1 **********************
**************************** globals.g **********************
******************************************************************************/
//Updated 6 November 2014 by Dan Dorman: ELEAK, RA, RM, CM, EREST_ACT, somaLen,gKIRsoma_UI, gKIRdend_UI were updated following optimization to hyperpolarized experimental current injection responses and changing soma morphology to cylindrical
str compartment = "compartment" //set compartment type
//Careful - all of these can be overridden in the .p file
if ({GABAtonic})
float ELEAK = -.18//.29
else
float ELEAK = -.072//-0.029478
end
float PI = 3.1415926
float RA = 1.25//13.7657 Range = 1 - 3
float RM = 1.875//8.97763 Range from KochSegev p175 = 2 - 10
float CM = 0.01//0.035 Range = .005 to .015
float EREST_ACT = -0.086//-.087
float TEMPERATURE = 35
float GHKluge=1e-7 //This param is required because getglobals only has 6 figs to right of decimal
//dividing by GHKluge is reversed in AddCaSpines.g and
//parameters determined by hand tuning to match spike width, AHP shape &, fI curve
float somaLen=11.3137e-6
float prox=26.1e-6// Primary only//42.1e-6 //26.1e-6 //primary and secondary
float mid=37.4e-6// Secondary and Teriary//60.1e-6 //62.1e-6 //tert1
float dist=1000e-6
float gNaFsoma_D1={47500} //50000 0.95x
float gNaFprox_D1={2593} //5700 6000 //includes mid
float gNaFdist_D1={926} //1900 2000
float gNaFsoma_D2={55000} //50000 1.1x
float gNaFprox_D2={2730} //6000 //includes mid
float gNaFdist_D2={975} //2000
float gNaFsoma_D1PD={47500} //50000 0.95x
float gNaFprox_D1PD={2593} //6000 //includes mid
float gNaFdist_D1PD={926} //2000
float gNaFsoma_D2PD={55000} //55000 1.1x
float gNaFprox_D2PD={2730} //6600 //includes mid
float gNaFdist_D2PD={975} //2200
float gNaFsoma_UI={46000} //50000 1.1x
float gNaFprox_UI={2600} //6000 //includes mid
float gNaFdist_UI={0} //2000
/*This set of KAf*/
float gKAfsoma_UI={200} //{300}
float gKAfprox_UI={200} //{550}
float gKAfdist_UI = {200}
//If gKAfdend is too low, get bursting, if too high, no spikes
float gKAfsoma_D1={200} //{300}
float gKAfdend_D1={100} //{550}
//If gKAfdend is too low, get bursting, if too high, no spikes
float gKAfsoma_D2={200} //{300}
float gKAfdend_D2={100} //{550}
//If gKAfdend is too low, get bursting, if too high, no spikes
float gKAfsoma_D1PD={200} //{300}
float gKAfdend_D1PD={100} //{550}
float gKAfsoma_D2PD={200} //{300}
float gKAfdend_D2PD={56} //{550}
float gKAssoma_D1={250}
float gKAsdend_D1={38.93} //includes mid
float gKAssoma_D2={275} //220 1.1x
float gKAsdend_D2={47.02} //22 //includes mid
float gKAssoma_D1PD={250}
float gKAsdend_D1PD={38.93} //includes mid
float gKAssoma_D2PD={275} //220 1.1x
float gKAsdend_D2PD={47.02} //22.2 //includes mid
float gKAssoma_UI={12} //220 1.1x
float gKAsdend_UI={12} //22 //includes mid
//Has effect on number of spikes, but not AHP amplitude
//Need to repeat sims with gKAs blocked - similar results to Shen .. Surmeier?
float gKIRsoma_D1={31.3322} //{11} 1.25x
float gKIRdend_D1={31.3322}
float gKIRsoma_D2={25.0658} //{11}
float gKIRdend_D2={25.0658}
float gKIRsoma_D1PD={31.3322} //{11} 1.25x
float gKIRdend_D1PD={31.3322}
float gKIRsoma_D2PD={25.0658} //{11}
float gKIRdend_D2PD={25.0658}
float gKIRsoma_UI={8.5} //{11}
float gKIRdend_UI={8.5}
float gKrpsoma={9} //10.008 is 14
float gKrpdend={8}
//Has effect on number of spikes, but not AHP amplitude
//Better to lower this and raise gKAs?
float gBKsoma={8}
float gBKdend={8}
//Higher value gave better AHPs.
//Increase in gBKsoma was compensated by lower gKAfdend and gKAssoma
//Perhaps a bit too high - does this prevent KAs from doing its thing?
float gSKsoma={1}
float gSKdend={1}
float gCaL13soma_D1 = {{3e-7}/{GHKluge}} //{{6e-7}/{Cafactor}}
float gCaL13dend_D1 = {{0.5e-8}/{GHKluge}} //{{gCaL13soma}/{Dendfactor}}
float gCaL13soma_D2 = {{0.75}*{3e-7}/{GHKluge}} //{{6e-7}/{Cafactor}} //0.75x
float gCaL13dend_D2 = {{0.75}*{0.5e-8}/{GHKluge}} //{{gCaL13soma}/{Dendfactor}}
float gCaL13soma_D1PD = {{3e-7}/{GHKluge}} //{{6e-7}/{Cafactor}}
float gCaL13dend_D1PD = {{0.5e-8}/{GHKluge}} //{{gCaL13soma}/{Dendfactor}}
float gCaL13soma_D2PD = {{0.75}*{3e-7}/{GHKluge}} //{{6e-7}/{Cafactor}} //0.75x
float gCaL13dend_D2PD = {{0.75}*{0.5e-8}/{GHKluge}} //{{gCaL13soma}/{Dendfactor}}
float gCaL13soma_UI = {{2e-7}/{GHKluge}} //{{6e-7}/{Cafactor}}
float gCaL13dend_UI = {{2e-7}/{GHKluge}} //{{gCaL13soma}/{Dendfactor}}
float gCaT32prox = {{0e-8}/{GHKluge}} //{{10e-8}/{Cafactor}}
float gCaT32dist = {{1e-7}/{GHKluge}}
float gCaT33prox = {{0e-8}/{GHKluge}} //{{10e-8}/{Cafactor}}
float gCaT33dist = {{.5e-8}/{GHKluge}}
float gCaRsoma = {{3e-7}/{GHKluge}} //{{18e-7}/{Cafactor}}
float gCaRdend = {{5e-7}/{GHKluge}} //{{gCaRsoma}/{Dendfactor}}
float gCaNsoma = {{1e-8}/{GHKluge}} //{{10e-7}/{Cafactor}}
float gCaNdend = {{0}/{GHKluge}} //{{gCaNsoma}/{Dendfactor}}
float gCaL12soma_D1 = {{2}*{6e-7}/{GHKluge}} //{{6e-7}/{Cafactor}} //2x
float gCaL12dend_D1 = {{2}*{1e-7}/{GHKluge}} //{{gCaL12soma}/{Dendfactor}}
float gCaL12soma_D2 = {{6e-7}/{GHKluge}} //{{6e-7}/{Cafactor}}
float gCaL12dend_D2 = {{1e-7}/{GHKluge}} //{{gCaL12soma}/{Dendfactor}}
float gCaL12soma_D1PD = {{2}*{6e-7}/{GHKluge}} //{{6e-7}/{Cafactor}} //2x
float gCaL12dend_D1PD = {{2}*{1e-7}/{GHKluge}} //{{gCaL12soma}/{Dendfactor}}
float gCaL12soma_D2PD = {{6e-7}/{GHKluge}} //{{6e-7}/{Cafactor}}
float gCaL12dend_D2PD = {{1e-7}/{GHKluge}} //{{gCaL12soma}/{Dendfactor}}
float gCaL12soma_UI = {{1e-7}/{GHKluge}} //{{6e-7}/{Cafactor}}
float gCaL12dend_UI = {{1e-7}/{GHKluge}} //{{gCaL12soma}/{Dendfactor}}
float gtonicGABAdend={650e-2} //tonic GABA parameters//1300e-2
float x0GABA = 40e-6 //sigmoidal function parameters//20e-6
float kGABA = 0.8e-4
float qfactorKir = 1.2
float qfactorKrp = 3 //3 2
float qfactorNaF = 2.5 //2.5 1.7
float qfactorkAs=2 //3 2
float qfactorkAf=2 //1.5 1.5
float qfactCa = 2 //3 2
//bk and sk qfacts are taken into account in their channel files.