Table 1. Voltage-dependent conductance parameters (from D’Angelo et al., 2001)

Conductance Vrev

state variables n (mV) a (sec-1) b (sec-1)

gNaF

Activation 3 0.9(V+19)/(1-exp(-(V+19)/10)) 36exp(-0.055(V+44))

Inactivation 1 87.39 0.315exp(-0.3(V+44)) 4.5(1+exp(-(V+11)/5))

gNaR

Activation 1 0.00024-0.015(V-4.5)/((exp(-(V-4.5)/6.8)-1)) 0.14+0.0.047(V+44)/exp((V+44)/0.11)-1))

Inactivation 1 87.39 0.96*exp(-(V+80)/62.5) 0.03*exp((V+83.3)/16.1)

gNaP

Activation 1 87.39 0.091(V+42)/(1-exp(-(V+42)/5)) 0.062(V+42)/(1-exp((V+42)/5))

x=1/(1+exp(-(V+42)/5))

t =5/(a +b )

gKV

Activation 4 -84.69 0.13(V+25)/(1-exp(-(V+25)/10)) 1.69exp(-0.0125(V+35))

gKA

Activation 3 14.67/(1+exp(-(V+9.17)23.32)) 2.98(exp(-(V+18.28)/19.47))

Inactivation 1 -84.69 0.33/(1+exp((V+111.33)/12.84)) 0.31/(1+exp(-(V+49.95)/8.9))

x=1/(1+exp(-(V+46.7)/19.8))

y=1/(1+exp((V+78.8)/8.4))

gKIR

Activation 1 -84.69 0.4exp(-0.041(V+83.94)) 0.51exp(0.028(V+83.94))

gKCa

Activation 1 -84.69 2.5/(1+(1.5e-3*exp(-0.085V)/[Ca])) 1.5/(1+[Ca]/(0.15 *10-3 *(-0.085V)))

gCaHVA

Activation 2 129.33 0.15exp(0.063(V+29.06)) 0.24exp(-0.039(V+18.66))

Inactivation 1 * 0.0039exp(-0.055(V+48)) 0.0039exp(0.012(V+48))

gKslow

Activation 1 -84.69 0.008exp(0.025(V+30)) 0.008exp(-0.05(V+30))

x=1/(1+exp(-(V+30)/6))

This table reports the equations used to calculate a x, b x, a y, b y (see Eqs. 3 and 4) for the membrane conductances used in the model. The number of gating particles (n) and ionic reversal potential(* resting value for Ca2+) used to calculate ionic currents are also shown.