Coherent and robust modulation of a metabolic network by cytoskeletal organization and dynamics

Generated code

• The required array sizes are:
  • VARIABLES[30]
  • CONSTANTS[72]
  • BOUND[1]
  • RATES[11]
• The following variables are defined:
  • Variable ADP(VARIABLES[11]; recomputed after rate changes)
  • Variable AMP(CONSTANTS[0])
  • Variable Cn(CONSTANTS[1])
  • Variable GTP(CONSTANTS[2])
  • Variable GDP(CONSTANTS[3])
  • Variable H(CONSTANTS[4])
  • Variable NADP(CONSTANTS[5])
  • Variable NADH(CONSTANTS[6])
  • Variable NAD(CONSTANTS[7])
  • Variable CD(VARIABLES[12]; recomputed after rate changes)
  • Variable CT(RATES[0] and VARIABLES[0])
  • Variable kpol(CONSTANTS[8])
  • Variable kf(CONSTANTS[9])
  • Variable kb(CONSTANTS[10])
  • Variable CP(RATES[1] and VARIABLES[1])
  • Variable kdp(CONSTANTS[11])
  • Variable CMTP(CONSTANTS[12])
  • Variable PKp(RATES[2] and VARIABLES[2])
  • Variable kp2(CONSTANTS[13])
  • Variable kp3(CONSTANTS[14])
  • Variable k4(CONSTANTS[15])
  • Variable PKt(VARIABLES[13]; recomputed after rate changes)
  • Variable C_PK(CONSTANTS[16])
  • Variable V_IN(VARIABLES[14]; recomputed after rate changes)
  • Variable Ke_in(CONSTANTS[17])
  • Variable KG_in(CONSTANTS[18])
  • Variable V_IN_max(CONSTANTS[19])
  • Variable V_HK(VARIABLES[15]; recomputed after rate changes)
  • Variable T(VARIABLES[16]; recomputed after rate changes)
  • Variable KpH(CONSTANTS[20])
  • Variable KPEP_r(CONSTANTS[21])
  • Variable KADP_r(CONSTANTS[22])
  • Variable KFDP_r(CONSTANTS[23])
  • Variable cADP(CONSTANTS[24])
  • Variable cFDP(CONSTANTS[25])
  • Variable cPEP(CONSTANTS[26])
  • Variable Lo_PK(CONSTANTS[27])
  • Variable gr_PK(CONSTANTS[28])
  • Variable gt_PK(CONSTANTS[29])
  • Variable n(VARIABLES[17]; recomputed after rate changes)
  • Variable V_PK_max(VARIABLES[18]; recomputed after rate changes)
  • Variable V_PKt_max(CONSTANTS[30])
  • Variable V_PKp_max(CONSTANTS[31])
  • Variable V_TCA(VARIABLES[19]; recomputed after rate changes)
  • Variable KPy_m(CONSTANTS[32])
  • Variable V_TCA_max(CONSTANTS[33])
  • Variable V_ADH(VARIABLES[20]; recomputed after rate changes)
  • Variable KPy__m(CONSTANTS[34])
  • Variable V_ADH_max(CONSTANTS[35])
  • Variable V_ATPase(VARIABLES[21]; recomputed after rate changes)
  • Variable KATP(CONSTANTS[36])
  • Variable time(BOUND[0])
  • Variable G_o(CONSTANTS[37])
  • Variable G(RATES[3] and VARIABLES[3])
  • Variable G6P(RATES[4] and VARIABLES[4])
  • Variable FDP(RATES[5] and VARIABLES[5])
  • Variable G3P(RATES[6] and VARIABLES[6])
  • Variable DPG(RATES[7] and VARIABLES[7])
  • Variable PEP(RATES[8] and VARIABLES[8])
  • Variable Py(RATES[9] and VARIABLES[9])
  • Variable ATP(RATES[10] and VARIABLES[10])
  • Variable PO(CONSTANTS[38])
  • Variable KG_m(CONSTANTS[39])
  • Variable KG_s(CONSTANTS[40])
  • Variable KATP_m(CONSTANTS[41])
  • Variable V_HK_max(CONSTANTS[42])
  • Variable V_PFK(VARIABLES[22]; recomputed after rate changes)
  • Variable KG6P_r(CONSTANTS[43])
  • Variable KATP_r(CONSTANTS[44])
  • Variable KAMP_r(CONSTANTS[45])
  • Variable cATP(CONSTANTS[46])
  • Variable cAMP(CONSTANTS[47])
  • Variable cG6P(CONSTANTS[48])
  • Variable Lo(CONSTANTS[49])
  • Variable gr(CONSTANTS[50])
  • Variable n1(CONSTANTS[51])
  • Variable V_PFK_max(CONSTANTS[52])
  • Variable V_G6PDH(VARIABLES[23]; recomputed after rate changes)
  • Variable TUB(VARIABLES[24]; recomputed after rate changes)
  • Variable KG6P(CONSTANTS[53])
  • Variable KNADP(CONSTANTS[54])
  • Variable KNADP_(CONSTANTS[55])
  • Variable KTUB(CONSTANTS[56])
  • Variable V_G6PDH_max(CONSTANTS[57])
  • Variable V_G6PDH_max_II(CONSTANTS[58])
  • Variable V_ALD(VARIABLES[25]; recomputed after rate changes)
  • Variable KG3P_m(CONSTANTS[59])
  • Variable KFDP_m(CONSTANTS[60])
  • Variable V_ALD_max(CONSTANTS[61])
  • Variable V_ALD_max_r(CONSTANTS[62])
  • Variable V_GAPDH(VARIABLES[26]; recomputed after rate changes)
  • Variable K1(CONSTANTS[63])
  • Variable K2(CONSTANTS[64])
  • Variable K3(CONSTANTS[65])
  • Variable KG3P(CONSTANTS[66])
  • Variable KNAD(CONSTANTS[67])
  • Variable KNADH_i(CONSTANTS[68])
  • Variable V_GAPDH_max(CONSTANTS[69])
  • Variable V_PGK(VARIABLES[27]; recomputed after rate changes)
  • Variable KDPG_m(CONSTANTS[70])
  • Variable V_PGK_max(CONSTANTS[71])
  • Variable V_PK(VARIABLES[28]; recomputed after rate changes)
  • Variable R(VARIABLES[29]; recomputed after rate changes)
• • Equation with cmeta ID no cmeta ID in math with cmeta ID 9
  • Equation with cmeta ID no cmeta ID in math with cmeta ID 12
  • Equation with cmeta ID no cmeta ID in math with cmeta ID 14

double NR_minfunc_1(double* CONSTANTS, double* VARIABLES, double* BOUND)
{
  return (((VARIABLES[10]+VARIABLES[11]+CONSTANTS[0])) - (CONSTANTS[1]));
}
double NR_minfunc_2(double* CONSTANTS, double* VARIABLES, double* BOUND)
{
  return (((VARIABLES[12]+VARIABLES[0]+VARIABLES[1])) - (CONSTANTS[12]));
}
double NR_minfunc_3(double* CONSTANTS, double* VARIABLES, double* BOUND)
{
  return (((VARIABLES[13]+VARIABLES[2])) - (CONSTANTS[16]));
}
void SetupFixedConstants(double* CONSTANTS)
{
CONSTANTS[0] = 0.5;
CONSTANTS[1] = 9;
CONSTANTS[2] = 0.95;
CONSTANTS[3] = 0.05;
CONSTANTS[4] = 3.2e-08;
CONSTANTS[5] = 1;
CONSTANTS[6] = 0.01;
CONSTANTS[7] = 1;
CONSTANTS[8] = 10;
CONSTANTS[9] = 3;
CONSTANTS[10] = 2.5;
CONSTANTS[11] = 0.0025;
CONSTANTS[12] = 0.9;
CONSTANTS[13] = 10;
CONSTANTS[14] = 0.05;
CONSTANTS[15] = 0.02;
CONSTANTS[16] = 0.1;
CONSTANTS[17] = 12;
CONSTANTS[18] = 0.001;
CONSTANTS[19] = 10;
CONSTANTS[20] = 9.5e-09;
CONSTANTS[21] = 1;
CONSTANTS[22] = 0.06;
CONSTANTS[23] = 0.025;
CONSTANTS[24] = 1;
CONSTANTS[25] = 0.01;
CONSTANTS[26] = 0.02;
CONSTANTS[27] = 1000;
CONSTANTS[28] = 0.1;
CONSTANTS[29] = 1;
CONSTANTS[30] = 25;
CONSTANTS[31] = 50;
CONSTANTS[32] = 0.329;
CONSTANTS[33] = 10;
CONSTANTS[34] = 0.169;
CONSTANTS[35] = 0.5;
CONSTANTS[36] = 5;
CONSTANTS[37] = 1;
CONSTANTS[38] = 4;
CONSTANTS[39] = 0.11;
CONSTANTS[40] = 0.0062;
CONSTANTS[41] = 0.1;
CONSTANTS[42] = 13;
CONSTANTS[43] = 1;
CONSTANTS[44] = 0.06;
CONSTANTS[45] = 0.025;
CONSTANTS[46] = 1;
CONSTANTS[47] = 0.019;
CONSTANTS[48] = 0.0005;
CONSTANTS[49] = 25000;
CONSTANTS[50] = 10;
CONSTANTS[51] = 2;
CONSTANTS[52] = 30;
CONSTANTS[53] = 0.05;
CONSTANTS[54] = 0.05;
CONSTANTS[55] = 0.05;
CONSTANTS[56] = 0.4;
CONSTANTS[57] = 1.6;
CONSTANTS[58] = 1;
CONSTANTS[59] = 20;
CONSTANTS[60] = 0.5;
CONSTANTS[61] = 2.5;
CONSTANTS[62] = 1;
CONSTANTS[63] = 1.1;
CONSTANTS[64] = 1.5;
CONSTANTS[65] = 2.5;
CONSTANTS[66] = 0.0025;
CONSTANTS[67] = 0.18;
CONSTANTS[68] = 0.0003;
CONSTANTS[69] = 10;
CONSTANTS[70] = 0.002;
CONSTANTS[71] = 3;
}

void SetupComputedConstants(double* CONSTANTS, double* VARIABLES)
{
VARIABLES[0] = 0;
VARIABLES[1] = 0;
VARIABLES[2] = 0;
VARIABLES[3] = 0;
VARIABLES[4] = 0;
VARIABLES[5] = 0;
VARIABLES[6] = 0;
VARIABLES[7] = 0;
VARIABLES[8] = 0;
VARIABLES[9] = 0;
VARIABLES[10] = 0;
}

void ComputeRates(double* BOUND, double* RATES, double* CONSTANTS, double* VARIABLES)
{
RATES[3] = (VARIABLES[14]-VARIABLES[15]);
RATES[4] = (VARIABLES[15]-(VARIABLES[22]+VARIABLES[23]));
RATES[5] = (VARIABLES[22]-VARIABLES[25]);
RATES[6] = ((2.00000*VARIABLES[25])-VARIABLES[26]);
RATES[7] = (VARIABLES[26]-VARIABLES[27]);
RATES[8] = (VARIABLES[27]-VARIABLES[28]);
RATES[9] = (VARIABLES[28]-(VARIABLES[19]+VARIABLES[20]));
RATES[10] = ((VARIABLES[27]+VARIABLES[28]+(CONSTANTS[38]*VARIABLES[19]))-(VARIABLES[15]+VARIABLES[22]+VARIABLES[21]));
RATES[0] = (-((CONSTANTS[8]*VARIABLES[0]*pow(VARIABLES[1],2.00000))+(CONSTANTS[9]*VARIABLES[12])+(CONSTANTS[10]*VARIABLES[0]*CONSTANTS[3])));
RATES[1] = ((CONSTANTS[8]*VARIABLES[0]*pow(VARIABLES[1],2.00000))-(CONSTANTS[11]*VARIABLES[1]));
RATES[2] = ((0.100000*CONSTANTS[13]*VARIABLES[13]*VARIABLES[1])-((CONSTANTS[14]*VARIABLES[2])+(CONSTANTS[15]*VARIABLES[2]*CONSTANTS[2])));
}

void ComputeVariables(double* BOUND, double* RATES, double* CONSTANTS, double* VARIABLES)
{
VARIABLES[25] = (((CONSTANTS[61]*(VARIABLES[5]/CONSTANTS[60]))-(CONSTANTS[62]*(VARIABLES[6]/CONSTANTS[59])))/(1.00000+(VARIABLES[5]/CONSTANTS[60])+(VARIABLES[6]/CONSTANTS[59])));
NR_MINIMISE(NR_minfunc_1, CONSTANTS, VARIABLES, BOUND, 11);
VARIABLES[26] = (CONSTANTS[69]/(1.00000+(CONSTANTS[66]/VARIABLES[6])+((CONSTANTS[67]/CONSTANTS[7])*(1.00000+(CONSTANTS[0]/CONSTANTS[63])+(VARIABLES[11]/CONSTANTS[64])+(VARIABLES[10]/CONSTANTS[65])))+(((CONSTANTS[66]*CONSTANTS[67])/(VARIABLES[6]*CONSTANTS[7]))*(1.00000+(CONSTANTS[6]/CONSTANTS[68])))+(1.00000+(CONSTANTS[0]/CONSTANTS[63])+(VARIABLES[11]/CONSTANTS[64])+(VARIABLES[10]/CONSTANTS[65]))));
VARIABLES[27] = ((CONSTANTS[71]*VARIABLES[7])/(CONSTANTS[70]+VARIABLES[7]));
VARIABLES[29] = (1.00000+(VARIABLES[8]/CONSTANTS[21])+(VARIABLES[11]/CONSTANTS[22])+(CONSTANTS[28]*(VARIABLES[8]/CONSTANTS[21])*(VARIABLES[11]/CONSTANTS[22])));
VARIABLES[16] = (1.00000+((CONSTANTS[26]*VARIABLES[8])/CONSTANTS[21])+((CONSTANTS[24]*VARIABLES[11])/CONSTANTS[22])+(CONSTANTS[29]*((CONSTANTS[26]*VARIABLES[8])/CONSTANTS[21])*((CONSTANTS[24]*VARIABLES[11])/CONSTANTS[22])));
VARIABLES[18] = (CONSTANTS[30]+((CONSTANTS[31]-CONSTANTS[30])*(VARIABLES[2]/CONSTANTS[16])));
VARIABLES[17] = (4.00000+(VARIABLES[2]/CONSTANTS[16]));
VARIABLES[28] = ((VARIABLES[18]/(1.00000+(CONSTANTS[20]/CONSTANTS[4])))*(((CONSTANTS[28]*(VARIABLES[8]/CONSTANTS[21])*(VARIABLES[11]/CONSTANTS[22])*pow(VARIABLES[29],VARIABLES[17]))+(CONSTANTS[27]*pow(((1.00000+((CONSTANTS[25]*VARIABLES[5])/CONSTANTS[23]))/(1.00000+(VARIABLES[5]/CONSTANTS[23]))),VARIABLES[17])*(VARIABLES[5]/CONSTANTS[23])*CONSTANTS[29]*((CONSTANTS[26]*VARIABLES[8])/CONSTANTS[21])*((CONSTANTS[24]*VARIABLES[11])/CONSTANTS[22])*pow(VARIABLES[16],(1.00000-VARIABLES[17]))))/(pow(VARIABLES[29],VARIABLES[17])+(CONSTANTS[27]*pow(((1.00000+((CONSTANTS[25]*VARIABLES[5])/CONSTANTS[23]))/(1.00000+(VARIABLES[5]/CONSTANTS[23]))),VARIABLES[17])*pow(VARIABLES[16],VARIABLES[17])))));
NR_MINIMISE(NR_minfunc_2, CONSTANTS, VARIABLES, BOUND, 12);
NR_MINIMISE(NR_minfunc_3, CONSTANTS, VARIABLES, BOUND, 13);
VARIABLES[14] = ((CONSTANTS[19]*(CONSTANTS[37]/((CONSTANTS[18]+CONSTANTS[37])*(1.00000+(VARIABLES[4]/CONSTANTS[17])))))-(VARIABLES[3]/((CONSTANTS[18]+VARIABLES[3])*(1.00000+(VARIABLES[4]/CONSTANTS[17])))));
VARIABLES[15] = (CONSTANTS[42]*(1.00000/(1.00000+((CONSTANTS[40]*CONSTANTS[41])/(VARIABLES[3]*VARIABLES[10]))+(CONSTANTS[39]/VARIABLES[3])+(CONSTANTS[41]/VARIABLES[10]))));
VARIABLES[19] = ((CONSTANTS[33]*pow(VARIABLES[9],2.00000))/(pow(CONSTANTS[32],2.00000)+pow(VARIABLES[9],2.00000)));
VARIABLES[20] = ((CONSTANTS[35]*VARIABLES[9])/(CONSTANTS[34]+VARIABLES[9]));
VARIABLES[21] = (CONSTANTS[36]*VARIABLES[10]);
VARIABLES[22] = (CONSTANTS[52]*((CONSTANTS[50]*(VARIABLES[4]/CONSTANTS[43])*(VARIABLES[10]/CONSTANTS[44])*pow((1.00000+(VARIABLES[4]/CONSTANTS[43])+(VARIABLES[10]/CONSTANTS[44])+(CONSTANTS[50]*(VARIABLES[4]/CONSTANTS[43])*(VARIABLES[10]/CONSTANTS[44]))),(CONSTANTS[51]-1.00000)))/(pow((1.00000+(VARIABLES[4]/CONSTANTS[43])+(VARIABLES[10]/CONSTANTS[44])+(CONSTANTS[50]*(VARIABLES[4]/CONSTANTS[43])*(VARIABLES[10]/CONSTANTS[44]))),CONSTANTS[51])+(CONSTANTS[49]*pow(((1.00000+((CONSTANTS[47]*CONSTANTS[0])/CONSTANTS[45]))/(1.00000+(CONSTANTS[0]/CONSTANTS[45]))),CONSTANTS[51])*pow((1.00000+((CONSTANTS[48]*VARIABLES[4])/CONSTANTS[43])+((CONSTANTS[46]*VARIABLES[10])/CONSTANTS[44])+(CONSTANTS[50]*((CONSTANTS[48]*VARIABLES[4])/CONSTANTS[43])*((CONSTANTS[46]*VARIABLES[10])/CONSTANTS[44]))),CONSTANTS[51])))));
VARIABLES[24] = (VARIABLES[0]+VARIABLES[12]);
VARIABLES[23] = ((CONSTANTS[57]/(((CONSTANTS[53]*CONSTANTS[54])/(VARIABLES[4]*CONSTANTS[5]))+(CONSTANTS[53]/VARIABLES[4])+(CONSTANTS[54]/CONSTANTS[5])+1.00000))+(CONSTANTS[58]/(((CONSTANTS[53]*CONSTANTS[55]*CONSTANTS[56])/(VARIABLES[4]*CONSTANTS[5]*VARIABLES[24]))+((CONSTANTS[53]*CONSTANTS[55])/(VARIABLES[4]*CONSTANTS[5]))+((CONSTANTS[55]*CONSTANTS[56])/(CONSTANTS[5]*VARIABLES[24]))+((CONSTANTS[53]*CONSTANTS[56])/(VARIABLES[4]*VARIABLES[24]))+(CONSTANTS[56]/VARIABLES[24])+(CONSTANTS[53]/VARIABLES[4])+(CONSTANTS[55]/CONSTANTS[5])+1.00000)));
#ifndef VARIABLES_FOR_RATES_ONLY
#endif
}