Emergent Properties of Networks of Biological Signaling Pathways

Generated code

• The required array sizes are:
  • VARIABLES[10]
  • CONSTANTS[30]
  • BOUND[1]
  • RATES[10]
• The following variables are defined:
  • Variable km13(CONSTANTS[0])
  • Variable Vmax13(CONSTANTS[1])
  • Variable GTP_RAS(RATES[0] and VARIABLES[0])
  • Variable GDP_RAS(RATES[1] and VARIABLES[1])
  • Variable kf4(CONSTANTS[2])
  • Variable kb4(CONSTANTS[3])
  • Variable km3(CONSTANTS[4])
  • Variable Vmax3(CONSTANTS[5])
  • Variable GAP(RATES[2] and VARIABLES[2])
  • Variable GAPstar(RATES[3] and VARIABLES[3])
  • Variable Ca4CAM(CONSTANTS[6])
  • Variable Gby(CONSTANTS[7])
  • Variable PKA(CONSTANTS[8])
  • Variable PKC(CONSTANTS[9])
  • Variable kb5(CONSTANTS[10])
  • Variable kf7(CONSTANTS[11])
  • Variable kb7(CONSTANTS[12])
  • Variable kf8(CONSTANTS[13])
  • Variable kb8(CONSTANTS[14])
  • Variable km1(CONSTANTS[15])
  • Variable Vmax1(CONSTANTS[16])
  • Variable km6(CONSTANTS[17])
  • Variable Vmax6(CONSTANTS[18])
  • Variable GEF_inact(RATES[4] and VARIABLES[4])
  • Variable GEF(RATES[5] and VARIABLES[5])
  • Variable Ca4CAM(RATES[6] and VARIABLES[6])
  • Variable Gby(RATES[7] and VARIABLES[7])
  • Variable CaM_GEF(RATES[8] and VARIABLES[8])
  • Variable Gby_GEF(RATES[9] and VARIABLES[9])
  • Variable kf12(CONSTANTS[19])
  • Variable kb12(CONSTANTS[20])
  • Variable km9(CONSTANTS[21])
  • Variable Vmax9(CONSTANTS[22])
  • Variable km10(CONSTANTS[23])
  • Variable Vmax10(CONSTANTS[24])
  • Variable km11(CONSTANTS[25])
  • Variable Vmax11(CONSTANTS[26])
  • Variable t(BOUND[0])
  • Variable kf2(CONSTANTS[27])
  • Variable kb2(CONSTANTS[28])
  • Variable kf5(CONSTANTS[29])
• No equation needed for Newton-Raphson evaluation.

void SetupFixedConstants(double* CONSTANTS)
{
CONSTANTS[0] = 1.0104;
CONSTANTS[1] = 10;
CONSTANTS[2] = 0.1;
CONSTANTS[3] = 0;
CONSTANTS[4] = 3.33333;
CONSTANTS[5] = 4;
CONSTANTS[6] = 1;
CONSTANTS[7] = 1;
CONSTANTS[8] = 1;
CONSTANTS[9] = 1;
CONSTANTS[10] = 1;
CONSTANTS[11] = 1;
CONSTANTS[12] = 0;
CONSTANTS[13] = 1e-05;
CONSTANTS[14] = 1;
CONSTANTS[15] = 7.5;
CONSTANTS[16] = 9;
CONSTANTS[17] = 3.33333;
CONSTANTS[18] = 4;
CONSTANTS[19] = 0.0001;
CONSTANTS[20] = 0;
CONSTANTS[21] = 0.50505;
CONSTANTS[22] = 0.02;
CONSTANTS[23] = 0.50505;
CONSTANTS[24] = 0.02;
CONSTANTS[25] = 0.50505;
CONSTANTS[26] = 0.02;
CONSTANTS[27] = 1;
CONSTANTS[28] = 0;
CONSTANTS[29] = 0.0001;
}

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

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

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