MTM_W_sPNaS_sICD.json 3.8 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123
  1. {
  2. "ode": [
  3. "C*(dv/dt) = - I_Na - I_K - I_L - i_p + I_app",
  4. "dm_Na/dt = alpha_m*(1-m_Na)-beta_m*m_Na",
  5. "dh_Na/dt = alpha_h*(1-h_Na)-beta_h*h_Na",
  6. "dn_K/dt = alpha_n*(1-n_K)-beta_n*n_K",
  7. "dNa_i/dt = rho/F*(-parstNa*i_p-I_Na-I_L_Na)*r_fast_vs_ion",
  8. "dK_i/dt = rho/F*(parstK*i_p-I_K-I_L_K)*r_fast_vs_ion",
  9. "dK_o/dt = rho/F*(parstK*i_p-I_K-I_L_K)*-remim*r_fast_vs_ion"
  10. ],
  11. "init_states": {
  12. "m_Na":"0.1",
  13. "h_Na":"0.6",
  14. "n_K":"0.4",
  15. "Na_i":"(10.0*mM)",
  16. "K_i":"(150.0*mM)",
  17. "K_o":"(8.0*mM)",
  18. "v":"(-60.0*mV)"
  19. },
  20. "definitions":
  21. {
  22. "I_Na": "g_Na*(m_Na**3)*h_Na*(v-E_Na)*Q10_gNa**((T-T0)/(10*kelvin))",
  23. "I_K": "g_K*(n_K**4)*(v-E_K)*Q10_gK**((T-T0)/(10*kelvin))",
  24. "I_L": "I_L_K+I_L_Na",
  25. "i_p": "(Piecewise((0,(Na_i<=(Na_sens0/2))),(1,(Na_i>(Na_sens0/2)))))*(I_max*(1/(1+exp(-Na_sens*(Na_i-Na_sens0)))-1/(1+exp(-Na_sens*-Na_sens0/2))))*Q10_gPump**((T-T0)/(10*kelvin))"
  26. },
  27. "functions":
  28. {
  29. "q_m": "Q10_m**((T-T0)/(10*kelvin))",
  30. "q_n": "Q10_n**((T-T0)/(10*kelvin))",
  31. "q_h": "Q10_h**((T-T0)/(10*kelvin))",
  32. "alpha_m": "0.32*(v/mV+alpha_mV)/(1-exp(-(v/mV+alpha_mV)/4))*q_m/ms",
  33. "beta_m": "0.28*(v/mV+beta_mV)/(exp((v/mV+beta_mV)/5)-1)*q_m/ms",
  34. "alpha_h": "0.128*exp(-(v/mV+alpha_hV)/18)*q_h/ms",
  35. "beta_h": "4/(1+exp(-(v/mV+beta_hV)/5))*q_h/ms",
  36. "alpha_n": "0.032*(v/mV+alpha_nV)/(1-exp(-(v/mV+alpha_nV)/5))*q_n/ms",
  37. "beta_n": "0.5*exp(-(v/mV+beta_nV)/40)*q_n/ms",
  38. "m_Na_inf":"alpha_m/(alpha_m+beta_m)",
  39. "h_Na_inf":"alpha_h/(alpha_h+beta_h)",
  40. "n_K_inf":"alpha_n/(alpha_n+beta_n)",
  41. "I_Nainf":"g_Na*(m_Na_inf**3)*h_Na_inf*(v-E_Na)*Q10_gNa**((T-T0)/(10*kelvin))",
  42. "I_Kinf": "g_K*(n_K_inf**4)*(v-E_K)*Q10_gK**((T-T0)/(10*kelvin))",
  43. "I_Linf": "g_L*(v-E_L)*Q10_gL**((T-T0)/(10*kelvin))",
  44. "i_pinf": "(Piecewise((0,Na_i<=(Na_sens0/2)),(1,Na_i>(Na_sens0/2))))*(I_max*(1/(1+exp(-Na_sens*(Na_i-Na_sens0)))-1/(1+exp(-Na_sens*-Na_sens0/2))))*Q10_gPump**((T-T0)/(10*kelvin))",
  45. "E_Na": "R*T/F*log(Na_o/Na_i)",
  46. "E_K":"R*T/F*log(K_o/K_i)",
  47. "E_L":"R*T/F*log((K_o*P_K+Na_o*P_Na)/(K_i*P_K+Na_i*P_Na))",
  48. "I_L_K": "g_L*P_K*(v-E_K)*Q10_gL**((T-T0)/(10*kelvin))",
  49. "I_L_Na": "g_L*P_Na*(v-E_Na)*Q10_gL**((T-T0)/(10*kelvin))"
  50. },
  51. "currents": {
  52. "I_Na": "I_Na",
  53. "I_K": "I_K",
  54. "I_L": "I_L",
  55. "i_p": "i_p"
  56. },
  57. "resting_membrane_pot": {
  58. "E_Na": "E_Na",
  59. "E_K":"E_K",
  60. "E_L":"E_L"
  61. },
  62. "ss_currents": {
  63. "I_Nainf":"I_Nainf",
  64. "I_Kinf": "I_Kinf",
  65. "I_Linf": "I_Linf",
  66. "i_pinf": "i_pinf"
  67. },
  68. "parameters": {
  69. "I_app": "1.0 * uA/cm2",
  70. "C": "1.0 * uF/cm2",
  71. "P_K":"(0.96)",
  72. "P_Na":"(0.04)",
  73. "F":"96484.6*coulomb/mol",
  74. "R":"8314.4*mjoule/(mol*kelvin)",
  75. "T": "293.15 * kelvin",
  76. "rho":"4000/cm",
  77. "remim":"0.2",
  78. "parstNa":"3",
  79. "parstK":"2",
  80. "g_Na":"100.0* msiemens/cm2",
  81. "g_K":"200.0* msiemens/cm2",
  82. "g_L":"0.1* msiemens/cm2",
  83. "alpha_mV":"54.0",
  84. "beta_mV":"27.0",
  85. "alpha_hV":"50.0",
  86. "beta_hV":"27.0",
  87. "alpha_nV":"52.0",
  88. "beta_nV":"57.0",
  89. "T0": "291.15 * kelvin",
  90. "Q10_gL":"1.2",
  91. "Q10_gNa":"1.2",
  92. "Q10_gK":"1.2",
  93. "Q10_gPump":"1.2",
  94. "Q10_n":"2.0",
  95. "Q10_m":"2.0",
  96. "Q10_h":"2.0",
  97. "I_max":"40.0*uA/cm2",
  98. "Na_sens0":"20.0*mM",
  99. "Na_sens":"0.1/mM",
  100. "i_p0":"-2.59284862*uA/cm2",
  101. "Na_o":"140.0*mM",
  102. "r_fast_vs_ion":"0.001*second/ms"
  103. },
  104. "source": "Miles Traub model ++ Na sens Pump, with dynamic concentrations"
  105. }