Model Mathematics

Component: environment

Component: sarcomere_geometry

sovr_ze=len_thick2iflen_thick2<SL2SL2otherwise sovr_cle=SL2-SL-len_thinifSL2-SL-len_thin>len_hbare2len_hbare2otherwise len_sovr=sovr_ze-sovr_cle SOVFThick=len_sovr2len_thick-len_hbare SOVFThin=len_sovrlen_thin

Component: Ca_binding_to_troponin_to_thin_filament_regulation

konT=kdHCam+H_consmkdHCam+HmkonQkonTmpC-31010 H=10e610-pH H_cons=10e610-7.15 koffLT=koffLkoffmodQkoffTmpC-31010 koffHT=koffHkoffmodQkoffTmpC-31010 kn_pT=kn_ppermtotQkn_pTmpC-31010 kp_nT=kp_ninprmtQkp_nTmpC-31010 permtot=|11+perm50Tropregnperm| inprmt=1permtotif1permtot<100100otherwise ddtimeTRPNCaL=dTRPNCaL dTRPNCaL=konTCai1-TRPNCaL-koffLTTRPNCaL ddtimeTRPNCaH=dTRPNCaH dTRPNCaH=konTCai1-TRPNCaH-koffHTTRPNCaH Tropreg=1-SOVFThinTRPNCaL+SOVFThinTRPNCaH

Component: thin_filament_regulation_and_crossbridge_cycling_rates

fappT=fappxbmodspQfappTmpC-31010 gapslmd=1+1-SOVFThickgslmod gappT=gappgapslmdxbmodspQgappTmpC-31010 hfmd=-xXBprer|xXBprer|hfmdcxXBprerx_02 hbmd=1 hfT=hfhfmdxbmodspQhfTmpC-31010 hbT=hbhbmdxbmodspQhbTmpC-31010 gxbmd=sigmapx_0-xXBpostrx_02ifxXBpostr<x_0sigmanxXBpostr-x_0x_02otherwise gxbT=gxbgxbmdxbmodspQgxbTmpC-31010

Component: regulation_and_crossbridge_cycling_state_equations

ddtimeN_NoXB=kp_nTP_NoXB-kn_pTN_NoXB ddtimeP_NoXB=kn_pTN_NoXB-kp_nTP_NoXB N=1-P+XBprer+XBpostr ddtimeP=kn_pTN-kp_nT+alpha1_plusP ddtimeXBprer=dXBprer dXBprer=alpha1_plusP+alpha2_minusXBpostr-alpha1_minus+alpha2_plusXBprer ddtimeXBpostr=dXBpostr dXBpostr=alpha3_minusP+alpha2_plusXBprer-alpha2_minus+alpha3_plusXBpostr alpha1_plus=fappT alpha2_plus=hfT alpha3_plus=MgATPgxbTMgATP_conskdADP+MgADP_conskdADP+MgADP alpha1_minus=xPigappTxPi_cons alpha2_minus=HhbTH_conskdADP+MgADP_consMgADP_consMgADPkdADP+MgADP alpha3_minus=fxbT fxbT=kdADPfappThfTgxbTMgATP_consgappTxPi_conshbTH_conskMgATP

Component: mean_strain_of_strongly_bound_states

ddtimexXBprer=dxXBprer dxXBprer=dSL2+xPsidutyprer-alpha1_plusxXBprer+alpha2_minusxXBpostr-x_0+xXBprer ddtimexXBpostr=dxXBpostr dxXBpostr=dSL2+xPsidutypostralpha2_plusxXBprer+x_0-xXBpostr dutyprer=alpha3_minusalpha2_minus+alpha3_plusalpha1_plus+alpha2_minusalpha1_plusalpha1_plusalpha2_plus+alpha3_minusalpha1_minus+alpha3_minusalpha2_plus+alpha3_minusalpha2_minus+alpha3_plusalpha1_plus+alpha2_minusalpha1_plus+alpha2_plusalpha3_plus+alpha3_minusalpha1_minus+alpha3_plusalpha1_minus dutypostr=alpha1_plusalpha2_plus+alpha3_minusalpha1_minus+alpha3_minusalpha2_plusalpha1_plusalpha2_plus+alpha3_minusalpha1_minus+alpha3_minusalpha2_plus+alpha3_minusalpha2_minus+alpha3_plusalpha1_plus+alpha2_minusalpha1_plus+alpha2_plusalpha3_plus+alpha3_minusalpha1_minus+alpha3_plusalpha1_minus

Component: normalised_active_and_passive_force

fxb=kdADPfapphfgxbMgATP_consgappxPi_conshbH_conskMgATP SSXBprer=fxbhb+gxbfapp+hbTfapphfgxb+hbgapp+gxbgapp+fxbhb+gxbfapp+hbfapp+fxbgapp+fapphf+fxbhb SSXBpostr=fapphf+fxbgapp+fxbhbhfgxb+hbgapp+gxbgapp+fxbhb+gxbfapp+hbfapp+fxbgapp+fapphf+fxbhb F_active=SOVFThickxXBpostrXBpostr+xXBprerXBprerx_0SSXBpostr ppforce_t=SL-SLrest|SL-SLrest|PCon_tPExp_t|SL-SLrest|-1 ppforce_c=PCon_cPExp_c|SL-SL_c|-1ifSL>SL_c0otherwise ppforce=ppforce_t+ppforce_c preload=|SLset-SLrest|SLset-SLrestPCon_tPExp_t|SLset-SLrest|-1 afterload=KSESLset-SLifSEon=10otherwise ddtimeSL=dSL dSL=0 ddtimeintf=preload+afterload-ppforce+F_active

Component: calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding

FrSBXB=XBpostr+XBprerSSXBpostr+SSXBprer dFrSBXB=dXBpostr+dXBprerSSXBpostr+SSXBprer dsovr_ze=-0.5dSLifSL<len_thick0otherwise dsovr_cle=-0.5dSLif2len_thin-SL>len_hbare0otherwise dlen_sovr=dsovr_ze-dsovr_cle dSOVFThin=dlen_sovrlen_thin dSOVFThick=2dlen_sovrlen_thick-len_hbare TropTot=Trop_conc1-SOVFThinTRPNCaL+SOVFThinFrSBXBTRPNCaH+1-FrSBXBTRPNCaL dTropTot=Trop_conc-dSOVFThinTRPNCaL+1-SOVFThindTRPNCaL+dSOVFThinFrSBXBTRPNCaH+1-FrSBXBTRPNCaL+SOVFThindFrSBXBTRPNCaH+FrSBXBdTRPNCaH+1-FrSBXBdTRPNCaL-dFrSBXBTRPNCaL

Component: model_parameters