# B_B = 0.0,
# B_C = 0.0,
# B_vhalf = 0.0,
# B_vslope = 33.5e-3)
# KDr_Y_params = []
Krpparam = ChannelSettings(Xpow=2, Ypow=1, Zpow=0, Erev=krev, name='Krp')
#Act tuned to fit Nisenbaum 1996 fig6C (minf^2) and fig 8C (mtau)
qfactKrp = 3 #Used by RE
Krp_X_params = AlphaBetaChannelParams(A_rate=48,
A_B=0,
A_C=0.0,
A_vhalf=0,
A_vslope=-0.02,
B_rate=7.199999999999999,
B_B=0.0,
B_C=0.0,
B_vhalf=0.0,
B_vslope=0.04)
# tuned to fit Nisenbaum 1996 fig 9D (hinf, 87% inactivating) and 9B (htau)
Krp_Y_params = TauInfMinChannelParams(T_min=0.8639999999999999,
T_vdep=12.48,
T_vhalf=-0.042,
T_vslope=0.013000000000000001,
SS_min=0.13,
SS_vdep=0.87,
SS_vhalf=-0.056,
SS_vslope=0.015,
T_power=1)
B_B = 0.0,
B_C = 0.0,
B_vhalf = 0.0,
B_vslope = 33.5e-3)
KDr_Y_params = []
'''
Krpparam = ChannelSettings(Xpow=2, Ypow=1, Zpow=0, Erev=krev, name='Krp')
# Act tuned to fit Nisenbaum 1996 fig6C (minf^2) and fig 8C (mtau)
qfactKrp = 3 # Used by RE
Krp_X_params = AlphaBetaChannelParams(A_rate=93.17594708227949,
A_B=0.0,
A_C=0.0,
A_vhalf=-0.0023501473737281243,
A_vslope=-0.02,
B_rate=13.976392062341922,
B_B=0.0,
B_C=0.0,
B_vhalf=-0.0023501473737281243,
B_vslope=0.04)
# tuned to fit Nisenbaum 1996 fig 9D (hinf, 87% inactivating) and 9B (htau)
Krp_Y_params = TauInfMinChannelParams(T_min=1.6771670474810307,
T_vdep=24.22574624139267,
T_vhalf=-0.044350147373728124,
T_vslope=0.013000000000000001,
SS_min=0.13,
SS_vdep=0.87,
SS_vhalf=-0.05835014737372812,
SS_vslope=0.015,
T_power=1)
T_vslope=3e-3,
SS_min=0.0,
SS_vdep=1.0,
SS_vhalf=-60e-3,
SS_vslope=6e-3)
NaFparam = ChannelSettings(Xpow=3, Ypow=1, Zpow=0, Erev=narev, name='NaF')
#This is from Migliore.
KDrparam = ChannelSettings(Xpow=1, Ypow=0, Zpow=0, Erev=krev, name='KDr')
KDr_X_params = AlphaBetaChannelParams(A_rate=28.2,
A_B=0,
A_C=0.0,
A_vhalf=0,
A_vslope=-12.5e-3,
B_rate=6.78,
B_B=0.0,
B_C=0.0,
B_vhalf=0.0,
B_vslope=33.5e-3)
KDr_Y_params = []
Krpparam = ChannelSettings(Xpow=2, Ypow=1, Zpow=0, Erev=krev, name='Krp')
#Act tuned to fit Nisenbaum 1996 fig6C (minf^2) and fig 8C (mtau)
qfactKrp = 3 #Used by RE
Krp_X_params = AlphaBetaChannelParams(A_rate=16 * qfactKrp,
A_B=0,
A_C=0.0,
A_vhalf=0,
ZpowCDI = 2
VMIN = -120e-3
VMAX = 50e-3
VDIVS = 3401 #0.5 mV steps
#units for calcium concentration: mM
CAMIN = 0.01e-3 #10 nM
CAMAX = 40e-3 #40 uM, might want to go up to 100 uM with spines
CADIVS = 4001 #10 nM steps
kDr_X_params = AlphaBetaChannelParams(A_rate=6.5,
A_B=0,
A_C=0,
A_vhalf=0,
A_vslope=-12.5e-3,
B_rate=24,
B_B=0,
B_C=0,
B_vhalf=0,
B_vslope=33.5e-3)
kDr_Y_params = []
kDrparams = ChannelSettings(Xpow=1, Ypow=0, Zpow=0, Erev=-0.09, name='kDr')
KAdist_m_params = AlphaBetaChannelParams(A_rate=0.7,
A_B=0,
A_C=0,
A_vhalf=0,
A_vslope=-50e-3,
B_rate=0.15,
T_vslope = 3e-3,
SS_min = 0.0,
SS_vdep = 1.0,
SS_vhalf = -60e-3,
SS_vslope = 6e-3)
NaFparam = ChannelSettings(Xpow=3, Ypow=1, Zpow=0, Erev=narev, name='NaF')
#This is from Migliore.
KDrparam = ChannelSettings(Xpow=1, Ypow=0, Zpow=0, Erev=krev, name='KDr')
KDr_X_params = AlphaBetaChannelParams(A_rate = 28.2,
A_B = 0,
A_C = 0.0,
A_vhalf = 0,
A_vslope = -12.5e-3,
B_rate = 6.78,
B_B = 0.0,
B_C = 0.0,
B_vhalf = 0.0,
B_vslope = 33.5e-3)
KDr_Y_params = []
Krpparam = ChannelSettings(Xpow=2, Ypow=1, Zpow=0, Erev=krev, name='Krp')
#Act tuned to fit Nisenbaum 1996 fig6C (minf^2) and fig 8C (mtau)
qfactKrp=3 #Used by RE
Krp_X_params = AlphaBetaChannelParams(A_rate = 16*qfactKrp,
A_B = 0,
A_C = 0.0,
A_vhalf = 0,
VMAX = 50e-3 # units(Volts) # Maximum allowed membrane potential.
VDIVS = 3000 # No units # Range of VMIN to VMAX is equally split in VDIVS sections. # 0.5 mV steps
# contains all gating parameters and reversal potentials
# Gate equations have the form:
# AlphaBetaChannelParams (specify forward and backward transition rates):
# alpha(v) or beta(v) = (rate + B * v) / (C + exp((v + vhalf) / vslope))
# descrition: alpha(v)/beta(v) are the implicit Action potential driving equations.
# alpha(v)/beta(v) equations adjustes the steady states[SS](alpha(v)/(alpha(v) + beta(v))
# and time constants (1/(alpha(v) + beta(v))).
# m,h,n are HH-channel gating variables
Na_m_params = AlphaBetaChannelParams(A_rate=-1e5 * (-25e-3 - EREST_ACT),
A_B=-1e5,
A_C=-1.0,
A_vhalf=-25e-3 - EREST_ACT,
A_vslope=-10e-3,
B_rate=4e3,
B_B=0.0,
B_C=0.0,
B_vhalf=0.0 - EREST_ACT,
B_vslope=18e-3)
Na_h_params = AlphaBetaChannelParams(A_rate=70,
A_B=0.0,
A_C=0.0,
A_vhalf=0.0 - EREST_ACT,
A_vslope=0.02,
B_rate=1000.0,
B_B=0.0,
B_C=1.0,
B_vhalf=-30e-3 - EREST_ACT,
B_vslope=-0.01)
# tau(v) = taumin + tauVdep / (1 + exp((v + tauVhalf) / tauVslope))
# or if tau_power=2: tau(v) = taumin + tauVdep / (1 + exp((v + tauVhalf) / tauVslope))* 1 / (1 + exp((v + tauVhalf) / -tauVslope))
#
# where v is membrane potential in volts, vhalf and vslope have units of volts
# C, min and max are dimensionless; and C should be either +1, 0 or -1
# Rate has units of per sec, and B has units of per sec per volt
# taumin and tauVdep have units of per sec
qfactNaF = 1.0
#These values were too fast - change rate from 35e3 to 16e3
Na_m_params = AlphaBetaChannelParams(A_rate=37792.709096491206,
A_B=0.0,
A_C=1,
A_vhalf=0.03533635901010369,
A_vslope=-0.005,
B_rate=37792.709096491206,
B_B=0.0,
B_C=1,
B_vhalf=0.03533635901010369,
B_vslope=0.005)
Na_h_params = AlphaBetaChannelParams(A_rate=5654.606181605965,
A_B=0.0,
A_C=1,
A_vhalf=0.07689738878303974,
A_vslope=0.009,
B_rate=5654.606181605965,
B_B=0.0,
B_C=1,
B_vhalf=0.03689738878303973,
B_vslope=-0.005)
# tau(v) = taumin + tauVdep / (1 + exp((v + tauVhalf) / tauVslope))
# or if tau_power=2: tau(v) = taumin + tauVdep / (1 + exp((v + tauVhalf) / tauVslope))* 1 / (1 + exp((v + tauVhalf) / -tauVslope))
#
# where v is membrane potential in volts, vhalf and vslope have units of volts
# C, min and max are dimensionless; and C should be either +1, 0 or -1
# Rate has units of per sec, and B has units of per sec per volt
# taumin and tauVdep have units of per sec
qfactNaF = 1.0
#These values were too fast - change rate from 35e3 to 16e3
Na_m_params = AlphaBetaChannelParams(A_rate=24160.0,
A_B=0,
A_C=1,
A_vhalf=0.0401,
A_vslope=-0.005,
B_rate=24160.0,
B_B=0.0,
B_C=1,
B_vhalf=0.0401,
B_vslope=0.005)
Na_h_params = AlphaBetaChannelParams(A_rate=6040.0,
A_B=0,
A_C=1,
A_vhalf=0.0761,
A_vslope=0.009,
B_rate=6040.0,
B_B=0.0,
B_C=1,
B_vhalf=0.0361,
B_vslope=-0.005)
CADIVS = 5999 #10 nM steps
#contains all gating parameters and reversal potentials
# Gate equations have the form:
# AlphaBetaChannelParams (specify forward and backward transition rates):
# alpha(v) or beta(v) = (rate + B * v) / (C + exp((v + vhalf) / vslope))
#Fast sodium - same as GP, but no slow gating
qfactNaF = 1.0
#These values were too fast - change rate from 35e3 to 12e3
Na_m_params = AlphaBetaChannelParams(A_rate=8e3,
A_B=0,
A_C=1,
A_vhalf=36e-3,
A_vslope=-5e-3,
B_rate=8e3,
B_B=0.0,
B_C=1,
B_vhalf=36e-3,
B_vslope=5e-3)
Na_h_params = AlphaBetaChannelParams(A_rate=4500,
A_B=0,
A_C=1,
A_vhalf=70e-3,
A_vslope=9e-3,
B_rate=3000,
B_B=0.0,
B_C=1,
B_vhalf=32e-3,
B_vslope=-5e-3)
VMAX = 50e-3
VDIVS = 3401 #0.5 mV steps
#units for calcium concentration: mM
CAMIN = 0.01e-3 #10 nM
CAMAX = 40e-3 #40 uM, might want to go up to 100 uM with spines
CADIVS = 4001 #10 nM steps
###### All channels from Kotaleski & Blackwell FSI model
qfactNaF = 1
NaF_X_params = AlphaBetaChannelParams(A_rate=3020e3,
A_B=-40.0e6,
A_C=-1,
A_vhalf=-75.5e-3,
A_vslope=-13.5e-3,
B_rate=1.2262e3,
B_B=0.0,
B_C=0.0,
B_vhalf=0.0,
B_vslope=42.248e-3)
NaF_Y_params = AlphaBetaChannelParams(A_rate=0.0035e3,
A_B=0.0,
A_C=0.0,
A_vhalf=0.0,
A_vslope=24.186e-3,
B_rate=-0.017e6 * 51.25e-3,
B_B=-0.017e6,
B_C=-1.0,
B_vhalf=51.25e-3,
B_vslope=-5.2e-3)
CAMAX = 60e-3 #40 uM, might want to go up to 100 uM with spines
CADIVS = 5999 #10 nM steps
#contains all gating parameters and reversal potentials
# Gate equations have the form:
# AlphaBetaChannelParams (specify forward and backward transition rates):
# alpha(v) or beta(v) = (rate + B * v) / (C + exp((v + vhalf) / vslope))
qfactNaF = 1.0
#These values were too fast - change rate from 35e3 to 16e3
Na_m_params = AlphaBetaChannelParams(A_rate=16e3,
A_B=0,
A_C=1,
A_vhalf=36e-3,
A_vslope=-5e-3,
B_rate=16e3,
B_B=0.0,
B_C=1,
B_vhalf=36e-3,
B_vslope=5e-3)
Na_h_params = AlphaBetaChannelParams(A_rate=4.0e3,
A_B=0,
A_C=1,
A_vhalf=72e-3,
A_vslope=9e-3,
B_rate=4e3,
B_B=0.0,
B_C=1,
B_vhalf=32e-3,
B_vslope=-5e-3)
# tau(v) = taumin + tauVdep / (1 + exp((v + tauVhalf) / tauVslope))
# or if tau_power=2: tau(v) = taumin + tauVdep / (1 + exp((v + tauVhalf) / tauVslope))* 1 / (1 + exp((v + tauVhalf) / -tauVslope))
#
# where v is membrane potential in volts, vhalf and vslope have units of volts
# C, min and max are dimensionless; and C should be either +1, 0 or -1
# Rate has units of per sec, and B has units of per sec per volt
# taumin and tauVdep have units of per sec
qfactNaF = 1.0
#These values were too fast - change rate from 35e3 to 16e3
Na_m_params = AlphaBetaChannelParams(A_rate=8929.205626365298,
A_B=0.0,
A_C=1,
A_vhalf=0.029425351443600374,
A_vslope=-0.005,
B_rate=8929.205626365298,
B_B=0.0,
B_C=1,
B_vhalf=0.029425351443600374,
B_vslope=0.005)
Na_h_params = AlphaBetaChannelParams(A_rate=2668.1033154895335,
A_B=0.0,
A_C=1,
A_vhalf=0.06443186709257158,
A_vslope=0.009,
B_rate=2668.1033154895335,
B_B=0.0,
B_C=1,
B_vhalf=0.02443186709257159,
B_vslope=-0.005)
# inf(v) = min + max / (1 + exp((v + vhalf) / vslope))
# tau(v) = taumin + tauVdep / (1 + exp((v + tauVhalf) / tauVslope))
# or if tau_power=2: tau(v) = taumin + tauVdep / (1 + exp((v + tauVhalf) / tauVslope))* 1 / (1 + exp((v + tauVhalf) / -tauVslope))
#
# where v is membrane potential in volts, vhalf and vslope have units of volts
# C, min and max are dimensionless; and C should be either +1, 0 or -1
# Rate has units of per sec, and B has units of per sec per volt
# taumin and tauVdep have units of per sec
qfactNaF = 2.0
NaF_X_params = AlphaBetaChannelParams(A_rate=-120.0e6 * -71.5e-3 * qfactNaF,
A_B=-120.0e6 * qfactNaF,
A_C=-1,
A_vhalf=-52.5e-3,
A_vslope=-9.5e-3,
B_rate=.36786e3 * qfactNaF,
B_B=0.0,
B_C=0.0,
B_vhalf=0.0,
B_vslope=32.248e-3)
NaF_Y_params = AlphaBetaChannelParams(A_rate=2000,
A_B=0.0,
A_C=1,
A_vhalf=0.067,
A_vslope=0.009,
B_rate=2000,
B_B=0.0,
B_C=1,
B_vhalf=0.027,
B_vslope=-0.015)
# tau(v) = taumin + tauVdep / (1 + exp((v + tauVhalf) / tauVslope)) # or if tau_power=2: tau(v) = taumin + tauVdep / (1 + exp((v + tauVhalf) / tauVslope))* 1 / (1 + exp((v + tauVhalf) / -tauVslope)) # # where v is membrane potential in volts, vhalf and vslope have units of volts # C, min and max are dimensionless; and C should be either +1, 0 or -1 # Rate has units of per sec, and B has units of per sec per volt # taumin and tauVdep have units of per sec qfactNaF = 1.0 #These values were too fast - change rate from 35e3 to 16e3 Na_m_params = AlphaBetaChannelParams(A_rate=48869.52576705926, A_B=0.0, A_C=1, A_vhalf=0.03256589042652369, A_vslope=-0.005, B_rate=48869.52576705926, B_B=0.0, B_C=1, B_vhalf=0.03256589042652369, B_vslope=0.005) Na_h_params = AlphaBetaChannelParams(A_rate=6810.764574104278, A_B=0.0, A_C=1, A_vhalf=0.07745302320358723, A_vslope=0.009, B_rate=6810.764574104278, B_B=0.0, B_C=1, B_vhalf=0.037453023203587225, B_vslope=-0.005)
# tau(v) = taumin + tauVdep / (1 + exp((v + tauVhalf) / tauVslope))
# or if tau_power=2: tau(v) = taumin + tauVdep / (1 + exp((v + tauVhalf) / tauVslope))* 1 / (1 + exp((v + tauVhalf) / -tauVslope))
#
# where v is membrane potential in volts, vhalf and vslope have units of volts
# C, min and max are dimensionless; and C should be either +1, 0 or -1
# Rate has units of per sec, and B has units of per sec per volt
# taumin and tauVdep have units of per sec
qfactNaF = 1.0
#These values were too fast - change rate from 35e3 to 16e3
Na_m_params = AlphaBetaChannelParams(A_rate=74805.13121562618,
A_B=0.0,
A_C=1,
A_vhalf=0.03185275093333961,
A_vslope=-0.005,
B_rate=74805.13121562618,
B_B=0.0,
B_C=1,
B_vhalf=0.03185275093333961,
B_vslope=0.005)
Na_h_params = AlphaBetaChannelParams(A_rate=8892.34887056553,
A_B=0.0,
A_C=1,
A_vhalf=0.07005373487255172,
A_vslope=0.009,
B_rate=8892.34887056553,
B_B=0.0,
B_C=1,
B_vhalf=0.030053734872551714,
B_vslope=-0.005)
# tau(v) = taumin + tauVdep / (1 + exp((v + tauVhalf) / tauVslope))
# or if tau_power=2: tau(v) = taumin + tauVdep / (1 + exp((v + tauVhalf) / tauVslope))* 1 / (1 + exp((v + tauVhalf) / -tauVslope))
#
# where v is membrane potential in volts, vhalf and vslope have units of volts
# C, min and max are dimensionless; and C should be either +1, 0 or -1
# Rate has units of per sec, and B has units of per sec per volt
# taumin and tauVdep have units of per sec
qfactNaF = 1.0
# These values were too fast - change rate from 35e3 to 16e3
Na_m_params = AlphaBetaChannelParams(A_rate=17191.515248231462,
A_B=0.0,
A_C=1,
A_vhalf=0.032543629783614365,
A_vslope=-0.005,
B_rate=17191.515248231462,
B_B=0.0,
B_C=1,
B_vhalf=0.032543629783614365,
B_vslope=0.005)
Na_h_params = AlphaBetaChannelParams(A_rate=2000.0011080946113,
A_B=0.0,
A_C=1,
A_vhalf=0.06964827061745116,
A_vslope=0.009,
B_rate=2000.0011080946113,
B_B=0.0,
B_C=1,
B_vhalf=0.029648270617451164,
B_vslope=-0.005)
# tau(v) = taumin + tauVdep / (1 + exp((v + tauVhalf) / tauVslope))
# or if tau_power=2: tau(v) = taumin + tauVdep / (1 + exp((v + tauVhalf) / tauVslope))* 1 / (1 + exp((v + tauVhalf) / -tauVslope))
#
# where v is membrane potential in volts, vhalf and vslope have units of volts
# C, min and max are dimensionless; and C should be either +1, 0 or -1
# Rate has units of per sec, and B has units of per sec per volt
# taumin and tauVdep have units of per sec
qfactNaF = 1.0
#These values were too fast - change rate from 35e3 to 16e3
Na_m_params = AlphaBetaChannelParams(A_rate=8929.205626365298,
A_B=0.0,
A_C=1,
A_vhalf=0.029425351443600374,
A_vslope=-0.005,
B_rate=8929.205626365298,
B_B=0.0,
B_C=1,
B_vhalf=0.029425351443600374,
B_vslope=0.005)
Na_h_params = AlphaBetaChannelParams(A_rate=2668.1033154895335,
A_B=0.0,
A_C=1,
A_vhalf=0.06443186709257158,
A_vslope=0.009,
B_rate=2668.1033154895335,
B_B=0.0,
B_C=1,
B_vhalf=0.02443186709257159,
B_vslope=-0.005)