Model Mathematics

Component: environment

Component: membrane

i_stim=stim_amplitudeiftimestim_starttimestim_endtime-stim_start-time-stim_startstim_periodstim_periodstim_duration0otherwise ddtimeV=-i_CaL+i_CaT+i_pCa+i_NaCa+i_Cab+i_Na+i_Nab+i_NaK+i_Kto_f+i_K1+i_Ks+i_Kur+i_Kss+i_Kr+i_ClCa+i_stim

Component: calcium_concentration

ddtimeCai=BiJ_leak+J_xfer-J_up+J_trpn ddtimeCass=BssJ_relVJSRVss-J_xferVmyoVss+i_CaL+i_Cab+i_pCa+i_CaT-2i_NaCaAcapCm2VssF ddtimeCaJSR=BJSRJ_xfer-J_rel ddtimeCaNSR=J_up-J_leakVmyoVNSR-J_trVJSRVNSR Bi=1+CMDN_totKm_CMDNKm_CMDN+Cai2-1 Bss=1+CMDN_totKm_CMDNKm_CMDN+Cass2-1 BJSR=1+CSQN_totKm_CSQNKm_CSQN+CaJSR2-1

Component: calcium_fluxes

J_rel=v1P_O1+P_O2CaJSR-CassP_RyR J_tr=CaNSR-CaJSRtau_tr J_xfer=Cass-Caitau_xfer J_leak=v2CaNSR-Cai J_up=v3Cai2Km_up2+Cai2 J_trpn=k_plus_htrpnCaiHTRPN_tot-HTRPN_Ca+k_plus_ltrpnCaiLTRPN_tot-LTRPN_Ca-k_minus_htrpnHTRPN_Ca+k_minus_ltrpnLTRPN_Ca ddtimeP_RyR=-0.04P_RyR-0.1i_CaLi_CaL_max-V-52648

Component: calcium_buffering

ddtimeLTRPN_Ca=k_plus_ltrpnCaiLTRPN_tot-LTRPN_Ca-k_minus_ltrpnLTRPN_Ca ddtimeHTRPN_Ca=k_plus_htrpnCaiHTRPN_tot-HTRPN_Ca-k_minus_htrpnHTRPN_Ca

Component: ryanodine_receptors

ddtimeP_O1=k_plus_aCassnP_C1+k_minus_bP_O2+k_minus_cP_C2-k_minus_aP_O1+k_plus_bCassmP_O1+k_plus_cP_O1 P_C1=1-P_C2+P_O1+P_O2 ddtimeP_O2=k_plus_bCassmP_O1-k_minus_bP_O2 ddtimeP_C2=k_plus_cP_O1-k_minus_cP_C2

Component: L_type_calcium_current

i_CaL=g_CaLOV-E_CaL ddtimeO=alphaC4+KpcbI1+0.001alphaI2-KpcfO-4betaO+gammaO C1=1-O+C2+C3+C4+I1+I2+I3 ddtimeC2=4alphaC1+2betaC3-betaC2+3alphaC2 ddtimeC3=3alphaC2+3betaC4-2betaC3+2alphaC3 ddtimeC4=2alphaC3+4betaO+0.014KpcbbetaI1-alphagammaC4+0.0024betaI2-KpcfC4+4betaKpcbI3-3betaC4+alphaC4+1gammaKpcfC4 ddtimeI1=gammaO+0.001alphaI3-KpcfI1+0.01alphagammaC4-4betaKpcfI1-KpcbI1 ddtimeI2=0.001KpcfO-alphaI2+KpcbI3+0.002KpcfC4-4betaI2-gammaI2 ddtimeI3=0.001KpcfI1-alphaI3+gammaI2+1gammaKpcfC4-4betaKpcbI3+KpcbI3 alpha=0.4V+12101+0.7-V+40210-0.75-V+2024001+0.12V+1210 beta=0.05-V+1213 gamma=Kpc_maxCassKpc_half+Cass Kpcf=131--V+14.52100

Component: T_type_calcium_current

i_CaT=g_CaTbgV-E_CaT

Component: T_type_calcium_current_b_gate

b_infinity=11+-V+486.1 tau_b=0.1+5.41+V+1006.6 ddtimeb=b_infinity-btau_b

Component: T_type_calcium_current_g_gate

g_infinity=11+V+666.6 tau_g=8+321+V+655 ddtimeg=g_infinity-gtau_g

Component: calcium_pump_current

i_pCa=i_pCa_maxCai2Km_pCa2+Cai2

Component: sodium_calcium_exchange_current

i_NaCa=k_NaCa1K_mNa3+Nao31K_mCa+Cao11+k_sateta-1VFRTetaVFRTNai3Cao-eta-1VFRTNao3Cai

Component: calcium_background_current

i_Cab=g_CabV-E_CaN E_CaN=RT2FlnCaoCai

Component: sodium_concentration

ddtimeNai=-i_Na+i_Nab+3i_NaK+3i_NaCaAcapCmVmyoF

Component: fast_sodium_current

i_Na=g_Nam3hjV-E_Na E_Na=RTFlnNaoNai

Component: fast_sodium_current_m_gate

alpha_m=0.32V+47.131--0.1V+47.13 beta_m=0.08-V11 ddtimem=alpha_m-alpha_m+beta_mm

Component: fast_sodium_current_h_gate

alpha_h=0.13580+V-6.8ifV<-400otherwise beta_h=3.560.079V+3100000.35VifV<-4010.131+V+10.66-11.1otherwise ddtimeh=alpha_h-alpha_h+beta_hh

Component: fast_sodium_current_j_gate

alpha_j=-1271400.2444V+3.474e-5-0.04391VV+37.781+0.311V+79.23ifV<-400otherwise beta_j=0.1212-0.01052V1+-0.1378V+40.14ifV<-400.3-2.535e-7V1+-0.1V+32otherwise ddtimej=alpha_j-alpha_j+beta_jj

Component: sodium_background_current

i_Nab=g_NabV-E_Na

Component: potassium_concentration

ddtimeKi=-i_Kto_f+i_K1+i_Ks+i_Kss+i_Kur+i_Kr-2i_NaKAcapCmVmyoF

Component: fast_transient_outward_potassium_current

i_Kto_f=g_Kto_fato_f6.5ito_fV-E_K E_K=RTFlnKoKi ddtimeato_f=alpha_a1-ato_f-beta_aato_f ddtimeito_f=ito_f_infinity-ito_ftau_ito_f ito_f_infinity=alpha_ialpha_i+beta_i tau_ito_f=10.000152-V+13.570.067083-V+33.57+1+0.00095V+33.570.051335V+33.57+1 alpha_a=0.180640.03577V+30 beta_a=0.3956-0.06237V+30 alpha_i=0.000152-V-3.8115.750.0067083-V+132.0515.75+1 beta_i=0.00095V+132.0515.750.051335V+132.0515.75+1

Component: ultra_rapidly_activating_delayed_rectifier_potassium_current

i_Kur=g_KurauriurV-E_K ddtimeaur=ass-aurtau_aur ddtimeiur=iss-iurtau_iur ass=11+-V+22.57.7 tau_aur=0.493-0.0629V+2.058 iss=11+V+45.25.7 tau_iur=1200-1701+V+45.25.7

Component: non_inactivating_steady_state_potassium_current

i_Kss=g_KssaKssV-E_K ddtimeaKss=ass-aKsstau_Kss tau_Kss=39.3-0.0862V+13.17

Component: time_independent_potassium_current

i_K1=g_K1KoKo+210V-E_K1+0.0896V-E_K

Component: slow_delayed_rectifier_potassium_current

i_Ks=g_KsnKs2V-E_K ddtimenKs=alpha_n1-nKs-beta_nnKs alpha_n=0.00000481333V+26.51--0.128V+26.5 beta_n=0.0000953333-0.038V+26.5

Component: rapid_delayed_rectifier_potassium_current

i_Kr=g_KrO_KV-E_K C_K0=1-C_K1+C_K2+O_K+I_K ddtimeC_K2=kfC_K1+beta_a1O_K-kbC_K2+alpha_a1C_K2 ddtimeC_K1=alpha_a0C_K0+kbC_K2-beta_a0C_K1+kfC_K1 ddtimeO_K=alpha_a1C_K2+beta_iI_K-beta_a1O_K+alpha_iO_K ddtimeI_K=alpha_iO_K-beta_iI_K alpha_a0=0.0223480.01176V beta_a0=0.047002-0.0631V alpha_a1=0.0137330.038198V beta_a1=0.0000689-0.04178V alpha_i=0.0908210.023391V+5 beta_i=0.006497-0.03268V+5

Component: sodium_potassium_pump_current

i_NaK=i_NaK_maxf_NaK11+Km_NaiNai1.5KoKo+Km_Ko f_NaK=11+0.1245-0.1VFRT+0.0365sigma-VFRT sigma=17Nao67300-1

Component: calcium_activated_chloride_current

f_ClCa=CassKm_Cl+Cassn_ClCa i_ClCa=1P_ClCaf_ClCaVF2RTCloVFRT-CliVFRT-1
Source
Derived from workspace Wang, Sobie, 2008 at changeset 4dad3d8d14e2.
Collaboration
To begin collaborating on this work, please use your git client and issue this command:
License
The terms of use/license for this work is unspecified.