Emergent Properties of Networks of Biological Signaling Pathways

Generated code

• The required array sizes are:
  • VARIABLES[11]
  • CONSTANTS[30]
  • BOUND[1]
  • RATES[11]
• The following variables are defined:
  • Variable km6(CONSTANTS[0])
  • Variable km7(CONSTANTS[1])
  • Variable km8(CONSTANTS[2])
  • Variable km9(CONSTANTS[3])
  • Variable Vmax6(CONSTANTS[4])
  • Variable Vmax7(CONSTANTS[5])
  • Variable Vmax8(CONSTANTS[6])
  • Variable Vmax9(CONSTANTS[7])
  • Variable PP2A(CONSTANTS[8])
  • Variable MAPKK(RATES[0] and VARIABLES[0])
  • Variable Vmax3(CONSTANTS[9])
  • Variable Vmax4(CONSTANTS[10])
  • Variable PP2A(CONSTANTS[11])
  • Variable Raf(RATES[1] and VARIABLES[1])
  • Variable Raf_star(RATES[2] and VARIABLES[2])
  • Variable Raf_star_star(RATES[3] and VARIABLES[3])
  • Variable GTPRas(RATES[4] and VARIABLES[4])
  • Variable GTPRasRaf_star(RATES[5] and VARIABLES[5])
  • Variable MAPKK_star(RATES[6] and VARIABLES[6])
  • Variable MAPKK_star_star(RATES[7] and VARIABLES[7])
  • Variable km10(CONSTANTS[12])
  • Variable km11(CONSTANTS[13])
  • Variable km12(CONSTANTS[14])
  • Variable km13(CONSTANTS[15])
  • Variable Vmax10(CONSTANTS[16])
  • Variable Vmax11(CONSTANTS[17])
  • Variable Vmax12(CONSTANTS[18])
  • Variable Vmax13(CONSTANTS[19])
  • Variable MAPK(RATES[8] and VARIABLES[8])
  • Variable MAPK_tyr(RATES[9] and VARIABLES[9])
  • Variable MAPK_star(RATES[10] and VARIABLES[10])
  • Variable MKP1(CONSTANTS[20])
  • Variable PKC(CONSTANTS[21])
  • Variable t(BOUND[0])
  • Variable kf5(CONSTANTS[22])
  • Variable kb5(CONSTANTS[23])
  • Variable km1(CONSTANTS[24])
  • Variable km2(CONSTANTS[25])
  • Variable km3(CONSTANTS[26])
  • Variable km4(CONSTANTS[27])
  • Variable Vmax1(CONSTANTS[28])
  • Variable Vmax2(CONSTANTS[29])
• No equation needed for Newton-Raphson evaluation.

void SetupFixedConstants(double* CONSTANTS)
{
CONSTANTS[0] = 0.159092;
CONSTANTS[1] = 0.159092;
CONSTANTS[2] = 15.6565;
CONSTANTS[3] = 15.6565;
CONSTANTS[4] = 0.105;
CONSTANTS[5] = 0.105;
CONSTANTS[6] = 6;
CONSTANTS[7] = 6;
CONSTANTS[8] = 0.224;
CONSTANTS[9] = 6;
CONSTANTS[10] = 6;
CONSTANTS[11] = 0.224;
CONSTANTS[12] = 0.0462967;
CONSTANTS[13] = 0.0462967;
CONSTANTS[14] = 0.0666667;
CONSTANTS[15] = 0.0666667;
CONSTANTS[16] = 0.15;
CONSTANTS[17] = 0.15;
CONSTANTS[18] = 1;
CONSTANTS[19] = 1;
CONSTANTS[20] = 0.032;
CONSTANTS[21] = 0.1;
CONSTANTS[22] = 4e-05;
CONSTANTS[23] = 0.5;
CONSTANTS[24] = 66.6667;
CONSTANTS[25] = 25.6417;
CONSTANTS[26] = 15.6565;
CONSTANTS[27] = 15.6565;
CONSTANTS[28] = 4;
CONSTANTS[29] = 10;
}

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

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

void ComputeVariables(double* BOUND, double* RATES, double* CONSTANTS, double* VARIABLES)
{
#ifndef VARIABLES_FOR_RATES_ONLY
#endif
}