Model Mathematics

Component: Environment

Component: Fixed_parameters

H_e=1⁢10-pH_e W_x=0.9⁢W_m W_i=0.1⁢W_m H_i=H_e K_i=K_e k_dHPi=1⁢10-6.75 k_dHatp=1⁢10-6.48 k_dHadp=1⁢10-6.29

Component: Adjustable_parameters

Component: Proton_motive_force

dG_H=F⁢dPsi+RT⁢ln⁡H_iH_x

Component: Dehydrogenase_flux

J_DH=x_DH⁢r⁢NAD_x-NADH_x⁢1+Pi_xk_Pi11+Pi_xk_Pi2

Component: Electron_flux_complex_I

dG_C1op=dG_C1o-RT⁢ln⁡H_x1 e -7-RT⁢ln⁡QQH2 J_C1=x_C1⁢ⅇ-dG_C1op+4⁢dG_HRT⁢NADH_x-NAD_x

Component: Electron_flux_complex_III

dG_C3op=dG_C3o+2⁢RT⁢ln⁡H_x1 e -7-RT⁢ln⁡QH2Q J_C3=x_C3⁢1+Pi_xk_Pi31+Pi_xk_Pi4⁢ⅇ-dG_C3op+4⁢dG_H-2⁢F⁢dPsi2⁢RT⁢Cox-Cred

Component: Electron_flux_complex_IV

dG_C4op=dG_C4o-2⁢RT⁢ln⁡H_x1 e -7-RT2⁢ln⁡O21 J_C4=x_C4⁢11+k_O2O2⁢CredCtot⁢ⅇ-dG_C4op+2⁢dG_H2⁢RT⁢Cred-Cox⁢ⅇF⁢dPsiRT

Component: ATP_synthesis_flux

J_F1=x_F1⁢ⅇ-dG_F1o-n_A⁢dG_HRT⁢K_DDK_DT⁢ADP_mx⁢Pi_x-ATP_mx⁢1

Component: ANT_flux

Psi_x=-0.65⁢dPsi Psi_i=0.35⁢dPsi J_ANT=x_ANT⁢ADP_fiADP_fi+ATP_fi⁢ⅇ-F⁢Psi_iRT-ADP_fxADP_fx+ATP_fx⁢ⅇ-F⁢Psi_xRT⁢ADP_fik_mADP+ADP_fiifADP_fi>mincond∨ATP_fi>mincond0otherwise

Component: MgATPx_binding_flux

ATP_fx=ATP_x-ATP_mx J_MgATPx=x_MgA⁢ATP_fx⁢Mg_x-K_DT⁢ATP_mx

Component: MgADPx_binding_flux

ADP_fx=ADP_x-ADP_mx J_MgADPx=x_MgA⁢ADP_fx⁢Mg_x-K_DD⁢ADP_mx

Component: MgATPi_binding_flux

ATP_fi=ATP_i-ATP_mi J_MgATPi=x_MgA⁢ATP_fi⁢Mg_i-K_DT⁢ATP_mi

Component: MgADPi_binding_flux

ADP_fi=ADP_i-ADP_mi J_MgADPi=x_MgA⁢ADP_fi⁢Mg_i-K_DD⁢ADP_mi

Component: ATP_substrate_flux

J_ATP=gamma⁢p_A⁢ATP_e-ATP_i

Component: ADP_substrate_flux

J_ADP=gamma⁢p_A⁢ADP_e-ADP_i

Component: AMP_substrate_flux

J_AMP=gamma⁢p_A⁢AMP_e-AMP_i

Component: Pi_substrate_flux

J_Pi2=gamma⁢x_Pi2⁢Pi_e-Pi_i

Component: Phosphate_hydrogen_cotransporter_flux

H2PIi=Pi_i⁢H_iH_i+k_dHPi H2PIx=Pi_x⁢H_xH_x+k_dHPi J_Pi1=x_Pi1⁢H_x⁢H2PIi-H_i⁢H2PIxH2PIi+k_PiH

Component: Adenylate_kinase_flux

J_AKi=x_AK⁢K_AK⁢ADP_i⁢ADP_i-AMP_i⁢ATP_i

Component: Hydrogen_leak_flux

J_Hle=x_Hle⁢dPsi⁢H_i⁢ⅇF⁢dPsiRT-H_xⅇF⁢dPsiRT-1

Component: Passive_potassium_flux

J_K=x_K⁢dPsi⁢K_i⁢ⅇF⁢dPsiRT-K_xⅇF⁢dPsiRT-1

Component: Potassium_hydrogen_flux

J_KH=x_KH⁢K_i⁢H_x-K_x⁢H_i

Component: NAD_x_concentration

NAD_x=NADtot-NADH_x

Component: Q_concentration

Q=Qtot-QH2

Component: Cox_concentration

Cox=Ctot-Cred

Component: ADP_binding

ADP_fe=ADP_e-ADP_me ADP_me=K_DD+ADP_e+Mg_tot-K_DD+ADP_e+Mg_tot2-4⁢Mg_tot⁢ADP_e2

Component: Mg_binding

Mg_e=Mg_tot-ADP_me Mg_i=Mg_e

Component: dH_x_dt

ddtH_x=x_buff⁢H_x⁢J_DH-5⁢J_C1-2⁢J_C3-4⁢J_C4+n_A-1⁢J_F1+2⁢J_Pi1+J_Hle-J_KHW_x

Component: dK_x_dt

ddtK_x=J_KH+J_KW_x

Component: dMg_x_dt

ddtMg_x=-J_MgATPx-J_MgADPxW_x

Component: dNADH_x_dt

ddtNADH_x=J_DH-J_C1W_x

Component: dQH2_dt

ddtQH2=J_C1-J_C3W_x

Component: dCred_dt

ddtCred=2⁢J_C3-2⁢J_C4W_i

Component: dATP_x_dt

ddtATP_x=J_F1-J_ANTW_x

Component: dADP_x_dt

ddtADP_x=-J_F1+J_ANTW_x

Component: dATP_mx_dt

ddtATP_mx=J_MgATPxW_x

Component: dADP_mx_dt

ddtADP_mx=J_MgADPxW_x

Component: dPi_x_dt

ddtPi_x=-J_F1+J_Pi1W_x

Component: dATP_i_dt

ddtATP_i=J_ATP+J_ANT+J_AKiW_i

Component: dADP_i_dt

ddtADP_i=J_ADP-J_ANT-2⁢J_AKiW_i

Component: dAMP_i_dt

ddtAMP_i=J_AMP+J_AKiW_i

Component: dATP_mi_dt

ddtATP_mi=J_MgATPiW_i

Component: dADP_mi_dt

ddtADP_mi=J_MgADPiW_i

Component: dPi_i_dt

ddtPi_i=-J_Pi1+J_Pi2W_i

Component: dPsi_dt

ddtdPsi=4⁢J_C1+2⁢J_C3+4⁢J_C4-n_A⁢J_F1-J_ANT-J_Hle-J_KC_im

Component: dO2_dt

ddtO2=0