//============================================================================== // CellML file: D:\ICC\ICC_model.cellml // CellML model: ICC_model // Date and time: 4/21/2008 at 1:36:00 at PM //------------------------------------------------------------------------------ // Conversion from CellML 1.0 to C was done using COR (0.9.31.781) // Copyright 2002-2008 Dr Alan Garny // http://COR.physiol.ox.ac.uk/ - COR@physiol.ox.ac.uk //------------------------------------------------------------------------------ // http://www.CellML.org/ //============================================================================== #include "ICC_model.h" //------------------------------------------------------------------------------ #include //------------------------------------------------------------------------------ // State variables //------------------------------------------------------------------------------ double Y[_NB_OF_STATE_VARIABLES_]; double dY[_NB_OF_STATE_VARIABLES_]; // 0: ICC_Membrane___Ca_i (millimolar) // 1: ICC_Membrane___Vm (voltage_units) // 2: PU_unit___ADP_i (millimolar) // 3: PU_unit___ADP_m (millimolar) // 4: PU_unit___Ca_ER (millimolar) // 5: PU_unit___Ca_PU (millimolar) // 6: PU_unit___Ca_m (millimolar) // 7: PU_unit___NADH_m (millimolar) // 8: PU_unit___deltaPsi (voltage_units) // 9: PU_unit___h (dimensionless) // 10: d_CaCl___d_CaCl (dimensionless) // 11: d_ERG___d_ERG (dimensionless) // 12: d_Ltype___d_Ltype (dimensionless) // 13: d_NSCC___d_NSCC (dimensionless) // 14: d_Na___d_Na (dimensionless) // 15: d_VDDR___d_VDDR (dimensionless) // 16: d_kv11___d_kv11 (dimensionless) // 17: f_Ltype___f_Ltype (dimensionless) // 18: f_Na___f_Na (dimensionless) // 19: f_VDDR___f_VDDR (dimensionless) // 20: f_ca_Ltype___f_ca_Ltype (dimensionless) // 21: f_kv11___f_kv11 (dimensionless) //------------------------------------------------------------------------------ // Constants //------------------------------------------------------------------------------ double Environment___Ca_o; // millimolar double Environment___Cl_o; // millimolar double Environment___F; // F_units double Environment___K_o; // millimolar double Environment___Na_o; // millimolar double Environment___Q10Ca; // dimensionless double Environment___Q10K; // dimensionless double Environment___Q10Na; // dimensionless double Environment___R; // R_units double Environment___T; // Temperature_units double Environment___T_exp; // Temperature_units double ICC_Membrane___Cl_i; // millimolar double ICC_Membrane___Cm; // capacitance_units double ICC_Membrane___K_i; // millimolar double ICC_Membrane___Na_i; // millimolar double ICC_Membrane___P_cyto; // dimensionless double ICC_Membrane___Vol; // volume_units double ICC_Membrane___fc; // dimensionless double I_BK___G_max_BK; // conductance_units double I_CaCl___G_max_CaCl; // conductance_units double I_ERG___G_max_ERG; // conductance_units double I_Ltype___G_max_Ltype; // conductance_units double I_NSCC___G_max_NSCC; // conductance_units double I_NSCC___NaPerm_o_Kperm; // dimensionless double I_Na___G_max_Na; // conductance_units double I_VDDR___G_max_VDDR; // conductance_units double I_bk___G_max_bk; // conductance_units double I_kv11___G_max_kv11; // conductance_units double J_PMCA___J_max_PMCA; // millimolar_per_second double PU_unit___Cmito; // millifarads double PU_unit___Glc; // millimolar double PU_unit___IP3; // millimolar double PU_unit___J_ERleak; // rate_constants_units double PU_unit___J_hyd_max; // millimolar_per_second double PU_unit___J_max_leak; // rate_constants_units double PU_unit___J_red_basal; // millimolar_per_second double PU_unit___Jmax_ANT; // millimolar_per_second double PU_unit___Jmax_IP3; // rate_constants_units double PU_unit___Jmax_NaCa; // millimolar_per_second double PU_unit___Jmax_serca; // millimolar_per_second double PU_unit___Jmax_uni; // rate_constants_units double PU_unit___KCa_PDH; // millimolar double PU_unit___K_Ca; // millimolar double PU_unit___K_F1; // millimolar double PU_unit___K_Glc; // millimolar double PU_unit___K_Na; // millimolar double PU_unit___K_act; // millimolar double PU_unit___K_hyd; // rate_constants_units double PU_unit___K_res; // dimensionless double PU_unit___K_trans; // millimolar double PU_unit___L; // dimensionless double PU_unit___P_ER; // dimensionless double PU_unit___P_PU; // dimensionless double PU_unit___P_mito; // dimensionless double PU_unit___Pi_m; // millimolar double PU_unit___b; // dimensionless double PU_unit___beta1; // per_millimolar double PU_unit___beta2; // per_millimolar double PU_unit___beta3; // per_millimolar double PU_unit___beta4; // per_millimolar double PU_unit___beta5; // per_millimolar double PU_unit___beta6; // per_millimolar double PU_unit___beta7; // per_millimolar double PU_unit___beta_max; // rate_constants_units double PU_unit___conc; // millimolar double PU_unit___d_ACT; // millimolar double PU_unit___d_INH; // millimolar double PU_unit___d_IP3; // millimolar double PU_unit___deltaPsi_B; // voltage_units double PU_unit___deltaPsi_star; // voltage_units double PU_unit___deltapH; // dimensionless double PU_unit___fe; // dimensionless double PU_unit___fm; // dimensionless double PU_unit___frac; // dimensionless double PU_unit___g; // dimensionless double PU_unit___g_H; // millimolar_per_second_per_millivolt double PU_unit___k_serca; // millimolar double PU_unit___n; // dimensionless double PU_unit___na; // dimensionless double PU_unit___nhyd; // dimensionless double PU_unit___p1; // dimensionless double PU_unit___p2; // dimensionless double PU_unit___p3; // dimensionless double PU_unit___pa; // rate_constants_units double PU_unit___pb; // rate_constants_units double PU_unit___pc1; // rate_constants_units double PU_unit___pc2; // rate_constants_units double PU_unit___r1; // dimensionless double PU_unit___r2; // dimensionless double PU_unit___r3; // dimensionless double PU_unit___ra; // rate_constants_units double PU_unit___rb; // rate_constants_units double PU_unit___rc1; // rate_constants_units double PU_unit___rc2; // rate_constants_units double PU_unit___rho_F1; // millimolar double PU_unit___rho_res; // millimolar double PU_unit___tauh; // time_units double PU_unit___total_ANP_i; // millimolar double PU_unit___total_ANP_m; // millimolar double PU_unit___total_NAD_m; // millimolar double PU_unit___u1; // dimensionless double PU_unit___u2; // dimensionless double d_CaCl___tau_d_CaCl; // time_units double d_NSCC___tau_d_NSCC; // time_units //------------------------------------------------------------------------------ // Computed variables //------------------------------------------------------------------------------ double Environment___FoRT; // Inverse_Voltage_units double Environment___RToF; // voltage_units double Environment___T_correction_BK; // conductance_units double Environment___T_correction_Ca; // dimensionless double Environment___T_correction_K; // dimensionless double Environment___T_correction_Na; // dimensionless double ICC_Membrane___V_cyto; // volume_units double I_BK___E_K; // voltage_units double I_BK___I_BK; // current_units double I_CaCl___E_Cl; // voltage_units double I_CaCl___I_CaCl; // current_units double I_ERG___E_K; // voltage_units double I_ERG___I_ERG; // current_units double I_Ltype___E_Ca; // voltage_units double I_Ltype___I_Ltype; // current_units double I_NSCC___E_NSCC; // voltage_units double I_NSCC___I_NSCC; // current_units double I_Na___E_Na; // voltage_units double I_Na___I_Na; // current_units double I_VDDR___E_Ca; // voltage_units double I_VDDR___I_VDDR; // current_units double I_bk___E_K; // voltage_units double I_bk___I_bk; // current_units double I_kv11___E_K; // voltage_units double I_kv11___I_kv11; // current_units double J_PMCA___J_PMCA; // millimolar_per_second double PU_unit___ADP3_i; // millimolar double PU_unit___ADP3_m; // millimolar double PU_unit___ADP_ifree; // millimolar double PU_unit___ADP_mfree; // millimolar double PU_unit___ATP4_i; // millimolar double PU_unit___ATP4_m; // millimolar double PU_unit___ATP_i; // millimolar double PU_unit___ATP_m; // millimolar double PU_unit___A_F1; // voltage_units double PU_unit___A_res; // voltage_units double PU_unit___J_ANT; // millimolar_per_second double PU_unit___J_ERout; // millimolar_per_second double PU_unit___J_HF1; // millimolar_per_second double PU_unit___J_Hleak; // millimolar_per_second double PU_unit___J_Hres; // millimolar_per_second double PU_unit___J_NaCa; // millimolar_per_second double PU_unit___J_SERCA; // millimolar_per_second double PU_unit___J_glyTotal; // millimolar_per_second double PU_unit___J_hyd; // millimolar_per_second double PU_unit___J_hydSS; // millimolar_per_second double PU_unit___J_leak; // millimolar_per_second double PU_unit___J_o; // millimolar_per_second double PU_unit___J_pF1; // millimolar_per_second double PU_unit___J_pGly; // millimolar_per_second double PU_unit___J_pTCA; // millimolar_per_second double PU_unit___J_red; // millimolar_per_second double PU_unit___J_uni; // millimolar_per_second double PU_unit___MWC; // millimolar double PU_unit___MgADP_i; // millimolar double PU_unit___NAD_m; // millimolar double PU_unit___PMF; // voltage_units double PU_unit___V_ER; // volume_units double PU_unit___V_MITO; // volume_units double PU_unit___V_PU; // volume_units double PU_unit___f_PDHa; // dimensionless double d_BK___d_BK; // dimensionless double d_CaCl___d_inf_CaCl; // dimensionless double d_ERG___d_inf_ERG; // dimensionless double d_ERG___tau_d_ERG; // time_units double d_Ltype___d_inf_Ltype; // dimensionless double d_Ltype___tau_d_Ltype; // time_units double d_NSCC___d_inf_NSCC; // dimensionless double d_Na___d_inf_Na; // dimensionless double d_Na___tau_d_Na; // time_units double d_VDDR___d_inf_VDDR; // dimensionless double d_VDDR___tau_d_VDDR; // time_units double d_kv11___d_inf_kv11; // dimensionless double d_kv11___tau_d_kv11; // time_units double f_Ltype___f_inf_Ltype; // dimensionless double f_Ltype___tau_f_Ltype; // time_units double f_Na___f_inf_Na; // dimensionless double f_Na___tau_f_Na; // time_units double f_VDDR___f_inf_VDDR; // dimensionless double f_VDDR___tau_f_VDDR; // time_units double f_ca_Ltype___f_ca_inf_Ltype; // dimensionless double f_ca_Ltype___tau_f_ca_Ltype; // time_units double f_kv11___f_inf_kv11; // dimensionless double f_kv11___tau_f_kv11; // time_units //------------------------------------------------------------------------------ // Initialisation //------------------------------------------------------------------------------ void init() { //--------------------------------------------------------------------------- // State variables //--------------------------------------------------------------------------- Y[0] = 0.00000993087; // ICC_Membrane___Ca_i (millimolar) Y[1] = -67.0; // ICC_Membrane___Vm (voltage_units) Y[2] = 0.0077282; // PU_unit___ADP_i (millimolar) Y[3] = 2.60093454; // PU_unit___ADP_m (millimolar) Y[4] = 0.007299; // PU_unit___Ca_ER (millimolar) Y[5] = 0.0000902; // PU_unit___Ca_PU (millimolar) Y[6] = 0.000136; // PU_unit___Ca_m (millimolar) Y[7] = 0.101476; // PU_unit___NADH_m (millimolar) Y[8] = 164.000044; // PU_unit___deltaPsi (voltage_units) Y[9] = 0.9397; // PU_unit___h (dimensionless) Y[10] = 0.0; // d_CaCl___d_CaCl (dimensionless) Y[11] = 0.0; // d_ERG___d_ERG (dimensionless) Y[12] = 0.0; // d_Ltype___d_Ltype (dimensionless) Y[13] = 0.0; // d_NSCC___d_NSCC (dimensionless) Y[14] = 0.0; // d_Na___d_Na (dimensionless) Y[15] = 0.0; // d_VDDR___d_VDDR (dimensionless) Y[16] = 0.0; // d_kv11___d_kv11 (dimensionless) Y[17] = 1.0; // f_Ltype___f_Ltype (dimensionless) Y[18] = 1.0; // f_Na___f_Na (dimensionless) Y[19] = 1.0; // f_VDDR___f_VDDR (dimensionless) Y[20] = 1.0; // f_ca_Ltype___f_ca_Ltype (dimensionless) Y[21] = 1.0; // f_kv11___f_kv11 (dimensionless) //--------------------------------------------------------------------------- // Constants //--------------------------------------------------------------------------- Environment___Ca_o = 2.5; // millimolar Environment___Cl_o = 134.0; // millimolar Environment___F = 96.4846; // F_units Environment___K_o = 7.0; // millimolar Environment___Na_o = 137.0; // millimolar Environment___Q10Ca = 2.1; // dimensionless Environment___Q10K = 1.5; // dimensionless Environment___Q10Na = 2.45; // dimensionless Environment___R = 8.3144; // R_units Environment___T = 310.0; // Temperature_units Environment___T_exp = 297.0; // Temperature_units ICC_Membrane___Cl_i = 88.0; // millimolar ICC_Membrane___Cm = 0.025; // capacitance_units ICC_Membrane___K_i = 120.0; // millimolar ICC_Membrane___Na_i = 30.0; // millimolar ICC_Membrane___P_cyto = 0.7; // dimensionless ICC_Membrane___Vol = 1.0e-12; // volume_units ICC_Membrane___fc = 0.01; // dimensionless I_BK___G_max_BK = 23.0; // conductance_units I_CaCl___G_max_CaCl = 10.1; // conductance_units I_ERG___G_max_ERG = 2.5; // conductance_units I_Ltype___G_max_Ltype = 2.0; // conductance_units I_NSCC___G_max_NSCC = 12.15; // conductance_units I_NSCC___NaPerm_o_Kperm = 1.056075; // dimensionless I_Na___G_max_Na = 20.0; // conductance_units I_VDDR___G_max_VDDR = 3.0; // conductance_units I_bk___G_max_bk = 0.15; // conductance_units I_kv11___G_max_kv11 = 6.3; // conductance_units J_PMCA___J_max_PMCA = 0.088464; // millimolar_per_second PU_unit___Cmito = 0.006995; // millifarads PU_unit___Glc = 1.0; // millimolar PU_unit___IP3 = 0.0006; // millimolar PU_unit___J_ERleak = 1.666667; // rate_constants_units PU_unit___J_hyd_max = 0.037625; // millimolar_per_second PU_unit___J_max_leak = 0.01; // rate_constants_units PU_unit___J_red_basal = 0.3333; // millimolar_per_second PU_unit___Jmax_ANT = 15.0; // millimolar_per_second PU_unit___Jmax_IP3 = 50000.0; // rate_constants_units PU_unit___Jmax_NaCa = 0.05; // millimolar_per_second PU_unit___Jmax_serca = 1.8333; // millimolar_per_second PU_unit___Jmax_uni = 5000.0; // rate_constants_units PU_unit___KCa_PDH = 0.00005; // millimolar PU_unit___K_Ca = 0.003; // millimolar PU_unit___K_F1 = 1.71e9; // millimolar PU_unit___K_Glc = 8.7; // millimolar PU_unit___K_Na = 9.4; // millimolar PU_unit___K_act = 0.00038; // millimolar PU_unit___K_hyd = 0.05125; // rate_constants_units PU_unit___K_res = 1.35e18; // dimensionless PU_unit___K_trans = 0.006; // millimolar PU_unit___L = 50.0; // dimensionless PU_unit___P_ER = 0.1; // dimensionless PU_unit___P_PU = 0.001; // dimensionless PU_unit___P_mito = 0.12871; // dimensionless PU_unit___Pi_m = 20.0; // millimolar PU_unit___b = 0.5; // dimensionless PU_unit___beta1 = 1.66; // per_millimolar PU_unit___beta2 = 0.0249; // per_millimolar PU_unit___beta3 = 4.0; // per_millimolar PU_unit___beta4 = 2.83; // per_millimolar PU_unit___beta5 = 1.3; // per_millimolar PU_unit___beta6 = 2.66; // per_millimolar PU_unit___beta7 = 0.16; // per_millimolar PU_unit___beta_max = 2.055; // rate_constants_units PU_unit___conc = 0.001; // millimolar PU_unit___d_ACT = 0.001; // millimolar PU_unit___d_INH = 0.0014; // millimolar PU_unit___d_IP3 = 0.00025; // millimolar PU_unit___deltaPsi_B = 50.0; // voltage_units PU_unit___deltaPsi_star = 91.0; // voltage_units PU_unit___deltapH = -0.4; // dimensionless PU_unit___fe = 0.01; // dimensionless PU_unit___fm = 0.0003; // dimensionless PU_unit___frac = 0.5; // dimensionless PU_unit___g = 0.85; // dimensionless PU_unit___g_H = 0.0033333; // millimolar_per_second_per_millivolt PU_unit___k_serca = 0.00042; // millimolar PU_unit___n = 2.0; // dimensionless PU_unit___na = 2.8; // dimensionless PU_unit___nhyd = 2.7; // dimensionless PU_unit___p1 = 1.346e-8; // dimensionless PU_unit___p2 = 7.739e-7; // dimensionless PU_unit___p3 = 6.65e-15; // dimensionless PU_unit___pa = 1.656e-5; // rate_constants_units PU_unit___pb = 3.373e-7; // rate_constants_units PU_unit___pc1 = 9.651e-14; // rate_constants_units PU_unit___pc2 = 4.845e-19; // rate_constants_units PU_unit___r1 = 2.077e-18; // dimensionless PU_unit___r2 = 1.728e-9; // dimensionless PU_unit___r3 = 1.059e-26; // dimensionless PU_unit___ra = 6.394e-10; // rate_constants_units PU_unit___rb = 1.762e-13; // rate_constants_units PU_unit___rc1 = 2.656e-19; // rate_constants_units PU_unit___rc2 = 8.632e-27; // rate_constants_units PU_unit___rho_F1 = 0.7; // millimolar PU_unit___rho_res = 0.4; // millimolar PU_unit___tauh = 4.0; // time_units PU_unit___total_ANP_i = 2.0; // millimolar PU_unit___total_ANP_m = 12.0; // millimolar PU_unit___total_NAD_m = 8.0; // millimolar PU_unit___u1 = 15.0; // dimensionless PU_unit___u2 = 1.1; // dimensionless d_CaCl___tau_d_CaCl = 0.03; // time_units d_NSCC___tau_d_NSCC = 0.35; // time_units } //------------------------------------------------------------------------------ // Computation //------------------------------------------------------------------------------ void compute(double t) { // t: time (time_units) Environment___FoRT = Environment___F/(Environment___R*Environment___T); Environment___RToF = Environment___R*Environment___T/Environment___F; Environment___T_correction_Ca = pow(Environment___Q10Ca, (Environment___T-Environment___T_exp)/10.0); Environment___T_correction_K = pow(Environment___Q10K, (Environment___T-Environment___T_exp)/10.0); Environment___T_correction_Na = pow(Environment___Q10Na, (Environment___T-Environment___T_exp)/10.0); Environment___T_correction_BK = 1.1*(Environment___T-Environment___T_exp); ICC_Membrane___V_cyto = ICC_Membrane___Vol*ICC_Membrane___P_cyto; I_Na___E_Na = Environment___RToF*log(Environment___Na_o/ICC_Membrane___Na_i); I_Na___I_Na = I_Na___G_max_Na*Y[18]*Y[14]*(Y[1]-I_Na___E_Na); I_Ltype___E_Ca = 0.5*Environment___RToF*log(Environment___Ca_o/Y[0]); I_Ltype___I_Ltype = I_Ltype___G_max_Ltype*Y[17]*Y[12]*Y[20]*(Y[1]-I_Ltype___E_Ca); I_VDDR___E_Ca = 0.5*Environment___RToF*log(Environment___Ca_o/Y[0]); I_VDDR___I_VDDR = I_VDDR___G_max_VDDR*Y[19]*Y[15]*(Y[1]-I_VDDR___E_Ca); I_kv11___E_K = Environment___RToF*log(Environment___K_o/ICC_Membrane___K_i); I_kv11___I_kv11 = I_kv11___G_max_kv11*Y[21]*Y[16]*(Y[1]-I_kv11___E_K); I_ERG___E_K = Environment___RToF*log(Environment___K_o/ICC_Membrane___K_i); I_ERG___I_ERG = I_ERG___G_max_ERG*Y[11]*(Y[1]-I_ERG___E_K); d_BK___d_BK = 1.0/(1.0+exp(Y[1]/-17.0-2.0*log(Y[0]/0.001))); I_BK___E_K = Environment___RToF*log(Environment___K_o/ICC_Membrane___K_i); I_BK___I_BK = (I_BK___G_max_BK+Environment___T_correction_BK)*d_BK___d_BK*(Y[1]-I_BK___E_K); I_CaCl___E_Cl = Environment___RToF*log(ICC_Membrane___Cl_i/Environment___Cl_o); I_CaCl___I_CaCl = I_CaCl___G_max_CaCl*Y[10]*(Y[1]-I_CaCl___E_Cl); I_NSCC___E_NSCC = Environment___RToF*log((Environment___K_o+Environment___Na_o*I_NSCC___NaPerm_o_Kperm)/(ICC_Membrane___K_i+ICC_Membrane___Na_i*I_NSCC___NaPerm_o_Kperm)); I_NSCC___I_NSCC = I_NSCC___G_max_NSCC*Y[13]*(Y[1]-I_NSCC___E_NSCC); I_bk___E_K = Environment___RToF*log(Environment___K_o/ICC_Membrane___K_i); I_bk___I_bk = I_bk___G_max_bk*(Y[1]-I_bk___E_K); J_PMCA___J_PMCA = J_PMCA___J_max_PMCA*1.0/(1.0+0.000298/Y[0]); dY[1] = -1.0*1.0/ICC_Membrane___Cm*(I_Na___I_Na+I_Ltype___I_Ltype+I_VDDR___I_VDDR+I_kv11___I_kv11+I_ERG___I_ERG+I_BK___I_BK+I_CaCl___I_CaCl+I_NSCC___I_NSCC+I_bk___I_bk+J_PMCA___J_PMCA*2.0*1000000.0*1000000.0*Environment___F*ICC_Membrane___V_cyto); PU_unit___J_leak = PU_unit___J_max_leak*(Y[5]-Y[0]); dY[0] = ICC_Membrane___fc*((-1.0*I_Ltype___I_Ltype+-1.0*I_VDDR___I_VDDR)/(2.0*1000000.0*1000000.0*Environment___F*ICC_Membrane___V_cyto)+PU_unit___J_leak+-1.0*J_PMCA___J_PMCA); PU_unit___V_MITO = ICC_Membrane___Vol*PU_unit___P_mito; PU_unit___V_PU = ICC_Membrane___Vol*PU_unit___P_PU; PU_unit___V_ER = ICC_Membrane___Vol*PU_unit___P_ER; PU_unit___J_ERout = (PU_unit___Jmax_IP3*pow(PU_unit___IP3/(PU_unit___IP3+PU_unit___d_IP3), 3.0)*pow(Y[5]/(Y[5]+PU_unit___d_ACT), 3.0)*pow(Y[9], 3.0)+PU_unit___J_ERleak)*(Y[4]-Y[5]); PU_unit___J_SERCA = PU_unit___Jmax_serca*pow(Y[5], 2.0)/(pow(PU_unit___k_serca, 2.0)+pow(Y[5], 2.0)); PU_unit___MWC = PU_unit___conc*Y[5]/PU_unit___K_trans*pow(1.0+Y[5]/PU_unit___K_trans, 3.0)/(pow(1.0+Y[5]/PU_unit___K_trans, 4.0)+PU_unit___L/pow(1.0+Y[5]/PU_unit___K_act, PU_unit___na)); PU_unit___J_uni = PU_unit___Jmax_uni*(PU_unit___MWC-Y[6]*exp(-2.0*Environment___FoRT*(Y[8]-PU_unit___deltaPsi_star)))*2.0*Environment___FoRT*(Y[8]-PU_unit___deltaPsi_star)/(1.0-exp(-2.0*Environment___FoRT*(Y[8]-PU_unit___deltaPsi_star))); PU_unit___J_NaCa = PU_unit___Jmax_NaCa*exp(PU_unit___b*Environment___FoRT*(Y[8]-PU_unit___deltaPsi_star))/((1.0+pow(PU_unit___K_Na/ICC_Membrane___Na_i, PU_unit___n))*(1.0+PU_unit___K_Ca/Y[6])); PU_unit___NAD_m = PU_unit___total_NAD_m-Y[7]; PU_unit___A_res = Environment___RToF*log(PU_unit___K_res*sqrt(Y[7])/sqrt(PU_unit___NAD_m)); PU_unit___J_o = PU_unit___rho_res*0.5*((PU_unit___ra*pow(10.0, 6.0*PU_unit___deltapH)+PU_unit___rc1*exp(6.0*PU_unit___deltaPsi_B*Environment___FoRT))*exp(PU_unit___A_res*Environment___FoRT)+-1.0*PU_unit___ra*exp(PU_unit___g*6.0*Environment___FoRT*Y[8])+PU_unit___rc2*exp(Environment___FoRT*PU_unit___A_res)*exp(Environment___FoRT*Y[8]*6.0*PU_unit___g))/((1.0+PU_unit___r1*exp(Environment___FoRT*PU_unit___A_res))*exp(Environment___FoRT*PU_unit___deltaPsi_B*6.0)+(PU_unit___r2+PU_unit___r3*exp(Environment___FoRT*PU_unit___A_res))*exp(Environment___FoRT*Y[8]*PU_unit___g*6.0)); PU_unit___J_Hres = PU_unit___rho_res*3.966*(PU_unit___ra*pow(10.0, 6.0*PU_unit___deltapH)*exp(Environment___FoRT*PU_unit___A_res)+PU_unit___rb*pow(10.0, 6.0*PU_unit___deltapH)+-1.0*(PU_unit___ra+PU_unit___rb)*exp(PU_unit___g*Environment___FoRT*Y[8]*6.0))/((1.0+PU_unit___r1*exp(Environment___FoRT*PU_unit___A_res))*exp(6.0*Environment___FoRT*PU_unit___deltaPsi_B)+(PU_unit___r2+PU_unit___r3*exp(Environment___FoRT*PU_unit___A_res))*exp(PU_unit___g*6.0*Environment___FoRT*Y[8])); PU_unit___ATP_i = PU_unit___total_ANP_i-Y[2]; PU_unit___J_glyTotal = PU_unit___beta_max*(1.0+PU_unit___beta1*PU_unit___Glc)*PU_unit___beta2*PU_unit___Glc*PU_unit___ATP_i/(1.0+PU_unit___beta3*PU_unit___ATP_i+(1.0+PU_unit___beta4*PU_unit___ATP_i)*PU_unit___beta5*PU_unit___Glc+(1.0+PU_unit___beta6*PU_unit___ATP_i)*PU_unit___beta7*PU_unit___Glc); PU_unit___f_PDHa = 1.0/(1.0+PU_unit___u2*(1.0+PU_unit___u1/pow(1.0+Y[6]/PU_unit___KCa_PDH, 2.0))); PU_unit___J_red = PU_unit___J_red_basal+6.3944*PU_unit___f_PDHa*PU_unit___J_glyTotal; PU_unit___J_pTCA = PU_unit___J_red_basal/3.0+0.84*PU_unit___f_PDHa*PU_unit___J_glyTotal; PU_unit___ATP_m = PU_unit___total_ANP_m-Y[3]; PU_unit___ADP_mfree = 0.8*Y[3]; PU_unit___A_F1 = Environment___RToF*log(PU_unit___K_F1*PU_unit___ATP_m/(PU_unit___ADP_mfree*PU_unit___Pi_m)); PU_unit___J_pF1 = PU_unit___rho_F1*-1.0*((PU_unit___pa*pow(10.0, 3.0*PU_unit___deltapH)+PU_unit___pc1*exp(3.0*Environment___FoRT*PU_unit___deltaPsi_B))*exp(Environment___FoRT*PU_unit___A_F1)+-1.0*PU_unit___pa*exp(3.0*Environment___FoRT*Y[8])+PU_unit___pc2*exp(Environment___FoRT*PU_unit___A_F1)*exp(3.0*Environment___FoRT*Y[8]))/((1.0+PU_unit___p1*exp(Environment___FoRT*PU_unit___A_F1))*exp(3.0*Environment___FoRT*PU_unit___deltaPsi_B)+(PU_unit___p2+PU_unit___p3*exp(Environment___FoRT*PU_unit___A_F1))*exp(3.0*Environment___FoRT*Y[8])); PU_unit___J_HF1 = -1.0*PU_unit___rho_F1*3.0*(PU_unit___pa*pow(10.0, 3.0*PU_unit___deltapH)*exp(Environment___FoRT*PU_unit___A_F1)+PU_unit___pb*pow(10.0, 3.0*PU_unit___deltapH)+-1.0*(PU_unit___pa+PU_unit___pb)*exp(3.0*Environment___FoRT*Y[8]))/((1.0+PU_unit___p1*exp(Environment___FoRT*PU_unit___A_F1))*exp(3.0*Environment___FoRT*PU_unit___deltaPsi_B)+(PU_unit___p2+PU_unit___p3*exp(Environment___FoRT*PU_unit___A_F1))*exp(3.0*Environment___FoRT*Y[8])); PU_unit___ATP4_i = 0.05*PU_unit___ATP_i; PU_unit___ADP3_m = 0.45*PU_unit___ADP_mfree; PU_unit___ADP_ifree = 0.3*Y[2]; PU_unit___ADP3_i = 0.45*PU_unit___ADP_ifree; PU_unit___ATP4_m = 0.05*PU_unit___ATP_m; PU_unit___J_ANT = PU_unit___Jmax_ANT*(1.0-PU_unit___ATP4_i*PU_unit___ADP3_m/(PU_unit___ADP3_i*PU_unit___ATP4_m)*exp(-1.0*Environment___FoRT*Y[8]))/((1.0+PU_unit___ATP4_i/PU_unit___ADP3_i*exp(-1.0*PU_unit___frac*Environment___FoRT*Y[8]))*(1.0+PU_unit___ADP3_m/PU_unit___ATP4_m)); PU_unit___PMF = Y[8]-2.303*Environment___RToF*PU_unit___deltapH; PU_unit___J_Hleak = PU_unit___g_H*PU_unit___PMF; PU_unit___J_pGly = 0.15*PU_unit___J_glyTotal; PU_unit___J_hydSS = PU_unit___J_hyd_max/(1.0+pow(PU_unit___K_Glc/PU_unit___Glc, PU_unit___nhyd)); PU_unit___J_hyd = PU_unit___K_hyd*PU_unit___ATP_i+PU_unit___J_hydSS; dY[7] = PU_unit___J_red-PU_unit___J_o; dY[3] = PU_unit___J_ANT+-1.0*PU_unit___J_pTCA+-1.0*PU_unit___J_pF1; dY[2] = -1.0*PU_unit___J_ANT*PU_unit___V_MITO/ICC_Membrane___V_cyto+PU_unit___J_hyd+-1.0*PU_unit___J_pGly; PU_unit___MgADP_i = 0.55*PU_unit___ADP_ifree; dY[5] = ICC_Membrane___fc*((PU_unit___J_NaCa-PU_unit___J_uni)*PU_unit___V_MITO/PU_unit___V_PU+(PU_unit___J_ERout-PU_unit___J_SERCA)*PU_unit___V_ER/PU_unit___V_PU+-1.0*PU_unit___J_leak*ICC_Membrane___V_cyto/PU_unit___V_PU); dY[6] = PU_unit___fm*(PU_unit___J_uni-PU_unit___J_NaCa); dY[4] = PU_unit___fe*(PU_unit___J_SERCA-PU_unit___J_ERout); dY[8] = -1.0*Environment___F*PU_unit___V_MITO*1000000.0*1.0/PU_unit___Cmito*(PU_unit___J_Hleak+-1.0*PU_unit___J_Hres+PU_unit___J_ANT+PU_unit___J_HF1+2.0*PU_unit___J_uni); dY[9] = 1.0*(PU_unit___d_INH-Y[9]*(Y[5]+PU_unit___d_INH))/PU_unit___tauh; d_CaCl___d_inf_CaCl = 1.0/(1.0+pow(0.00014/Y[0], 3.0)); dY[10] = (d_CaCl___d_inf_CaCl-Y[10])/d_CaCl___tau_d_CaCl; d_ERG___d_inf_ERG = 0.2+0.8/(1.0+exp((Y[1]+20.0)/-1.8)); d_ERG___tau_d_ERG = Environment___T_correction_K*0.003; dY[11] = (d_ERG___d_inf_ERG-Y[11])/d_ERG___tau_d_ERG; d_Ltype___d_inf_Ltype = 1.0/(1.0+exp((Y[1]+17.0)/-4.3)); d_Ltype___tau_d_Ltype = Environment___T_correction_Ca*0.001; dY[12] = (d_Ltype___d_inf_Ltype-Y[12])/d_Ltype___tau_d_Ltype; d_NSCC___d_inf_NSCC = 1.0/(1.0+pow(0.0000745/Y[5], -85.0)); dY[13] = (d_NSCC___d_inf_NSCC-Y[13])/d_NSCC___tau_d_NSCC; d_Na___d_inf_Na = 1.0/(1.0+exp((Y[1]+47.0)/-4.8)); d_Na___tau_d_Na = Environment___T_correction_Na*0.003; dY[14] = (d_Na___d_inf_Na-Y[14])/d_Na___tau_d_Na; d_VDDR___d_inf_VDDR = 1.0/(1.0+exp((Y[1]+26.0)/-6.0)); d_VDDR___tau_d_VDDR = Environment___T_correction_Ca*0.006; dY[15] = (d_VDDR___d_inf_VDDR-Y[15])/d_VDDR___tau_d_VDDR; d_kv11___d_inf_kv11 = 1.0/(1.0+exp((Y[1]+25.0)/-7.7)); d_kv11___tau_d_kv11 = Environment___T_correction_K*0.005; dY[16] = (d_kv11___d_inf_kv11-Y[16])/d_kv11___tau_d_kv11; f_Ltype___f_inf_Ltype = 1.0/(1.0+exp((Y[1]+43.0)/8.9)); f_Ltype___tau_f_Ltype = Environment___T_correction_Ca*0.086; dY[17] = (f_Ltype___f_inf_Ltype-Y[17])/f_Ltype___tau_f_Ltype; f_Na___f_inf_Na = 1.0/(1.0+exp((Y[1]+78.0)/7.0)); f_Na___tau_f_Na = Environment___T_correction_Na*0.0016; dY[18] = (f_Na___f_inf_Na-Y[18])/f_Na___tau_f_Na; f_VDDR___f_inf_VDDR = 1.0/(1.0+exp((Y[1]+66.0)/6.0)); f_VDDR___tau_f_VDDR = Environment___T_correction_Ca*0.04; dY[19] = (f_VDDR___f_inf_VDDR-Y[19])/f_VDDR___tau_f_VDDR; f_ca_Ltype___f_ca_inf_Ltype = 1.0-1.0/(1.0+exp((Y[0]-0.0001-0.000214)/-0.0000131)); f_ca_Ltype___tau_f_ca_Ltype = Environment___T_correction_Ca*0.002; dY[20] = (f_ca_Ltype___f_ca_inf_Ltype-Y[20])/f_ca_Ltype___tau_f_ca_Ltype; f_kv11___f_inf_kv11 = 0.5+0.5/(1.0+exp((Y[1]+44.8)/4.4)); f_kv11___tau_f_kv11 = Environment___T_correction_K*0.005; dY[21] = (f_kv11___f_inf_kv11-Y[21])/f_kv11___tau_f_kv11; } //============================================================================== // End of file //==============================================================================