Kv3_X_params = AlphaBetaChannelParams(A_rate=3794.01343804835, A_B=0.0, A_C=1, A_vhalf=-0.035602703223899766, A_vslope=-0.015, B_rate=5159.858275745756, B_B=0.0, B_C=1, B_vhalf=0.07539729677610023, B_vslope=0.015) Kv3_Y_params = TauInfMinChannelParams(T_min=0.005472828692040172, T_vdep=0.02032764942757778, T_vhalf=-0.010800038578430718, T_vslope=0.01, SS_min=0.6, SS_vdep=0.4, SS_vhalf=-0.030800038578430715, SS_vslope=0.01, T_power=1) KvFparam = ChannelSettings(Xpow=4, Ypow=1, Zpow=0, Erev=krev, name='KvF') KvF_X_params = AlphaBetaChannelParams(A_rate=3597.6784536945975, A_B=0.0, A_C=1, A_vhalf=-0.029645387917126724, A_vslope=-0.026, B_rate=3597.6784536945975, B_B=0.0,
Na_s_params= AlphaBetaChannelParams(A_rate = 100,
A_B = 0,
A_C = 1,
A_vhalf = 84e-3,
A_vslope = 10.0e-3,
B_rate = 80,
B_B = 0.0,
B_C = 1,
B_vhalf = -26e-3,
B_vslope = -14e-3)
'''
Na_s_params = TauInfMinChannelParams(T_min=0.5,
T_vdep=2.2,
T_vhalf=-0.032,
T_vslope=0.012,
SS_min=0.15,
SS_vdep=0.85,
SS_vhalf=-0.045,
SS_vslope=0.0054,
T_power=1)
NaFparam = ChannelSettings(Xpow=3, Ypow=1, Zpow=0, Erev=narev, name='NaF')
# persistent Na_channel. Slow down from Dieter Jaeger
'''
NaS_m_params = AlphaBetaChannelParams(A_rate=76e3,
A_B=0,
A_C=1,
A_vhalf=-55.4e-3,
A_vslope=-13.6e-3,
B_rate=70e3,
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
#mtau: Ogata fig 5, no qfactor accounted in mtau, 1.2 will improve spike shape
#activation minf fits Ogata 1990 figure 3C (which is cubed root)
#inactivation hinf fits Ogata 1990 figure 6B
#htau fits the main -50 through -10 slope of Ogata figure 9 (log tau), but a qfact of 2 is already taken into account.
qfactNaF = 2.5
Na_m_params = TauInfMinChannelParams(T_min=0.00011916138208411449,
T_vdep=0.002505368058318507,
T_vhalf=-0.06180983895507333,
T_vslope=0.008,
SS_min=0.0,
SS_vdep=1.0,
SS_vhalf=-0.024809838955073337,
SS_vslope=-0.01,
T_power=2)
Na_h_params = TauInfMinChannelParams(T_min=0.0005699726479054598,
T_vdep=0.002483540949479127,
T_vhalf=-0.03390383342558563,
T_vslope=0.003,
SS_min=0.0,
SS_vdep=1.0,
SS_vhalf=-0.051903833425585626,
SS_vslope=0.006,
T_power=1)
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
# mtau: Ogata fig 5, no qfactor accounted in mtau, 1.2 will improve spike shape
# activation minf fits Ogata 1990 figure 3C (which is cubed root)
# inactivation hinf fits Ogata 1990 figure 6B
# htau fits the main -50 through -10 slope of Ogata figure 9 (log tau), but a qfact of 2 is already taken into account.
qfactNaF = 2.5
Na_m_params = TauInfMinChannelParams(T_min=7.998875749153482e-05,
T_vdep=0.0016817636262595197,
T_vhalf=-0.05308849295680853,
T_vslope=0.008,
SS_min=0.0,
SS_vdep=1.0,
SS_vhalf=-0.01608849295680853,
SS_vslope=-0.01,
T_power=2)
Na_h_params = TauInfMinChannelParams(T_min=0.00027542296648548886,
T_vdep=0.0012001000718321953,
T_vhalf=-0.034845326172449564,
T_vslope=0.003,
SS_min=0.0,
SS_vdep=1.0,
SS_vhalf=-0.05284532617244956,
SS_vslope=0.006,
T_power=1)
Krp_X_params = AlphaBetaChannelParams(A_rate=16,
A_B=0,
A_C=0.0,
A_vhalf=0,
A_vslope=-20e-3,
B_rate=2.4,
B_B=0.0,
B_C=0.0,
B_vhalf=0.0,
B_vslope=40e-3)
Krp_Y_params = TauInfMinChannelParams(T_min=0.287,
T_vdep=4.16,
T_vhalf=-0.042,
T_vslope=0.013,
SS_min=0.13,
SS_vdep=0.87,
SS_vhalf=-0.056,
SS_vslope=0.014)
Krpparam = ChannelSettings(Xpow=2, Ypow=1, Zpow=0, Erev=krev, name='Krp')
SK_Z_params = ZChannelParams(Kd=0.00035,
power=4.6,
tau=0.002,
taumax=0.0037928,
tau_power=4.3,
cahalf=0.002703)
SKparam = ChannelSettings(Xpow=0, Ypow=0, Zpow=1, Erev=krev, name='SKCa')
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
#mtau: Ogata fig 5, no qfactor accounted in mtau, 1.2 will improve spike shape
#activation minf fits Ogata 1990 figure 3C (which is cubed root)
#inactivation hinf fits Ogata 1990 figure 6B
#htau fits the main -50 through -10 slope of Ogata figure 9 (log tau), but a qfact of 2 is already taken into account.
qfactNaF = 2.5
Na_m_params = TauInfMinChannelParams(SS_min=0.0,
SS_vdep=1.0,
SS_vhalf=-25e-3,
SS_vslope=-10e-3,
T_min=0.1e-3 / qfactNaF,
T_vdep=2.1025e-3 / qfactNaF,
T_vhalf=-62e-3,
T_vslope=8e-3,
T_power=2)
Na_h_params = TauInfMinChannelParams(T_min=2 * 0.2754e-3 / qfactNaF,
T_vdep=2 * 1.2e-3 / qfactNaF,
T_vhalf=-42e-3,
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')
Kv3_X_params = AlphaBetaChannelParams(A_rate=6562.550002507868,
A_B=0.0,
A_C=1,
A_vhalf=-0.026020030073901267,
A_vslope=-0.015,
B_rate=8925.0680034107,
B_B=0.0,
B_C=1,
B_vhalf=0.08497996992609873,
B_vslope=0.015)
Kv3_Y_params = TauInfMinChannelParams(T_min=0.005508571871081583,
T_vdep=0.020460409806874448,
T_vhalf=-0.01282016812665924,
T_vslope=0.01,
SS_min=0.6,
SS_vdep=0.4,
SS_vhalf=-0.032820168126659235,
SS_vslope=0.01,
T_power=1)
KvFparam = ChannelSettings(Xpow=4, Ypow=1, Zpow=0, Erev=krev, name='KvF')
KvF_X_params = AlphaBetaChannelParams(A_rate=2929.086037975104,
A_B=0.0,
A_C=1,
A_vhalf=-0.03862764078731916,
A_vslope=-0.026,
B_rate=2929.086037975104,
B_B=0.0,
B_C=1,
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
#mtau: Ogata fig 5, no qfactor accounted in mtau, 1.2 will improve spike shape
#activation minf fits Ogata 1990 figure 3C (which is cubed root)
#inactivation hinf fits Ogata 1990 figure 6B
#htau fits the main -50 through -10 slope of Ogata figure 9 (log tau), but a qfact of 2 is already taken into account.
qfactNaF = 2.5
Na_m_params = TauInfMinChannelParams(SS_min=0.0,
SS_vdep=1.0,
SS_vhalf=-25e-3,
SS_vslope=-10e-3,
T_min=0.1e-3 / qfactNaF,
T_vdep=2.1025e-3 / qfactNaF,
T_vhalf=-62e-3,
T_vslope=8e-3,
T_power=2)
Na_h_params = TauInfMinChannelParams(T_min=2 * 0.2754e-3 / qfactNaF,
T_vdep=2 * 1.2e-3 / qfactNaF,
T_vhalf=-42e-3,
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')
A_B=-0.014e6, #per msec/mV e3 converts to per sec,e3 converts to per V
A_C=-1,
A_vhalf=44e-3,
A_vslope=-2.3e-3,
B_rate=0.0043e3,
B_B=0,
B_C=0,
B_vhalf=44e-3,
B_vslope=34e-3)
Kv13param = ChannelSettings(Xpow=4, Ypow=0, Zpow=0, Erev=krev, name='Kv13')
Ka_X_params = TauInfMinChannelParams(T_min=1e-3,
T_vdep=1e-3,
T_vhalf=0.07,
T_vslope=-0.013,
SS_min=0,
SS_vdep=1,
SS_vhalf=-0.045,
SS_vslope=-0.013)
Ka_Y_params = TauInfMinChannelParams(T_min=0.014,
T_vdep=0,
T_vhalf=37.0e-3,
T_vslope=5.0e-3,
SS_min=0,
SS_vdep=1,
SS_vhalf=-0.077,
SS_vslope=0.008)
Kaparam = ChannelSettings(Xpow=4, Ypow=1, Zpow=0, Erev=carev, name='Ka')
#Dictionary of "standard" channels, to create channels using a loop
Kv3_X_params = AlphaBetaChannelParams(A_rate=3189.410941127774,
A_B=0.0,
A_C=1,
A_vhalf=-0.03776099637932748,
A_vslope=-0.015,
B_rate=4337.598879933773,
B_B=0.0,
B_C=1,
B_vhalf=0.07323900362067252,
B_vslope=0.015)
Kv3_Y_params = TauInfMinChannelParams(T_min=0.008738675971393061,
T_vdep=0.03245793932231708,
T_vhalf=-0.005683681756781953,
T_vslope=0.01,
SS_min=0.6,
SS_vdep=0.4,
SS_vhalf=-0.025683681756781953,
SS_vslope=0.01,
T_power=1)
KvFparam = ChannelSettings(Xpow=4, Ypow=1, Zpow=0, Erev=krev, name='KvF')
KvF_X_params = AlphaBetaChannelParams(A_rate=3261.161091702947,
A_B=0.0,
A_C=1,
A_vhalf=-0.022018415451993303,
A_vslope=-0.026,
B_rate=3261.161091702947,
B_B=0.0,
B_C=1,
Kv3_X_params = AlphaBetaChannelParams(A_rate=4000,
A_B=0,
A_C=1,
A_vhalf=-33e-3,
A_vslope=-15e-3,
B_rate=5440,
B_B=0.0,
B_C=1,
B_vhalf=78e-3,
B_vslope=15e-3)
#Use SSTau Quadratic to better match minimal inactivation of this channel
Kv3_Y_params = TauInfMinChannelParams(SS_min=0.6,
SS_vdep=0.4,
SS_vhalf=-20e-3,
SS_vslope=10e-3,
T_min=0.007,
T_vdep=0.026,
T_vhalf=0e-3,
T_vslope=10e-3)
KvFparam = ChannelSettings(Xpow=4, Ypow=1, Zpow=0, Erev=krev, name='KvF')
KvF_X_params = AlphaBetaChannelParams(A_rate=3e3,
A_B=0,
A_C=1,
A_vhalf=-32e-3,
A_vslope=-26.0e-3,
B_rate=3e3,
B_B=0.0,
B_C=1,
B_vhalf=130e-3,
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
#mtau: Ogata fig 5, no qfactor accounted in mtau, 1.2 will improve spike shape
#activation minf fits Ogata 1990 figure 3C (which is cubed root)
#inactivation hinf fits Ogata 1990 figure 6B
#htau fits the main -50 through -10 slope of Ogata figure 9 (log tau), but a qfact of 2 is already taken into account.
qfactNaF = 2.5
Na_m_params = TauInfMinChannelParams(SS_min=0.0,
SS_vdep=1.0,
SS_vhalf=-25e-3,
SS_vslope=-10e-3,
T_min=0.1e-3 / qfactNaF,
T_vdep=2.1025e-3 / qfactNaF,
T_vhalf=-62e-3,
T_vslope=8e-3,
T_power=2)
Na_h_params = TauInfMinChannelParams(T_min=2 * 0.2754e-3 / qfactNaF,
T_vdep=2 * 1.2e-3 / qfactNaF,
T_vhalf=-42e-3,
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')
Na_s_params= AlphaBetaChannelParams(A_rate = 100,
A_B = 0,
A_C = 1,
A_vhalf = 84e-3,
A_vslope = 10.0e-3,
B_rate = 80,
B_B = 0.0,
B_C = 1,
B_vhalf = -26e-3,
B_vslope = -14e-3)
'''
Na_s_params = TauInfMinChannelParams(T_min=0.5,
T_vdep=2.2,
T_vhalf=-0.032,
T_vslope=0.012,
SS_min=0.15,
SS_vdep=0.85,
SS_vhalf=-0.045,
SS_vslope=0.0054,
T_power=1)
NaFparam = ChannelSettings(Xpow=3, Ypow=1, Zpow=0, Erev=narev, name='NaF')
#persistent Na_channel. Slow down from Dieter Jaeger
'''
NaS_m_params = AlphaBetaChannelParams(A_rate=76e3,
A_B=0,
A_C=1,
A_vhalf=-55.4e-3,
A_vslope=-13.6e-3,
B_rate=70e3,
Na_s_params= AlphaBetaChannelParams(A_rate = 100,
A_B = 0,
A_C = 1,
A_vhalf = 84e-3,
A_vslope = 10.0e-3,
B_rate = 80,
B_B = 0.0,
B_C = 1,
B_vhalf = -26e-3,
B_vslope = -14e-3)
'''
Na_s_params = TauInfMinChannelParams(SS_min=0.15,
SS_vdep=0.85,
SS_vhalf=-0.045,
SS_vslope=0.0054,
T_min=0.01,
T_vdep=2.2,
T_vhalf=-0.032,
T_vslope=0.012,
T_power=2)
NaFparam = ChannelSettings(Xpow=3, Ypow=1, Zpow=0, Erev=narev, name='NaF')
# persistent Na_channel. Slow down from Dieter Jaeger
'''
NaS_m_params = AlphaBetaChannelParams(A_rate=76e3,
A_B=0,
A_C=1,
A_vhalf=-55.4e-3,
A_vslope=-13.6e-3,
B_rate=70e3,
Na_s_params= AlphaBetaChannelParams(A_rate = 100,
A_B = 0,
A_C = 1,
A_vhalf = 84e-3,
A_vslope = 10.0e-3,
B_rate = 80,
B_B = 0.0,
B_C = 1,
B_vhalf = -26e-3,
B_vslope = -14e-3)
'''
Na_s_params = TauInfMinChannelParams(T_min=0.01,
T_vdep=2.2,
T_vhalf=-0.032,
T_vslope=0.012,
SS_min=0.15,
SS_vdep=0.85,
SS_vhalf=-0.038349999999999995,
SS_vslope=0.0054,
T_power=2)
NaFparam = ChannelSettings(Xpow=3, Ypow=1, Zpow=0, Erev=narev, name='NaF')
#persistent Na_channel. Slow down from Dieter Jaeger
'''
NaS_m_params = AlphaBetaChannelParams(A_rate=76e3,
A_B=0,
A_C=1,
A_vhalf=-55.4e-3,
A_vslope=-13.6e-3,
B_rate=70e3,
