Emergent Properties of Networks of Biological Signaling Pathways

Generated code

• The required array sizes are:
  • VARIABLES[12]
  • CONSTANTS[27]
  • BOUND[1]
  • RATES[12]
• The following variables are defined:
  • Variable t(BOUND[0])
  • Variable kf1(CONSTANTS[0])
  • Variable kb1(CONSTANTS[1])
  • Variable kf3(CONSTANTS[2])
  • Variable km8(CONSTANTS[3])
  • Variable km12(CONSTANTS[4])
  • Variable Vmax2(CONSTANTS[5])
  • Variable Vmax4(CONSTANTS[6])
  • Variable Vmax6(CONSTANTS[7])
  • Variable Vmax8(CONSTANTS[8])
  • Variable Vmax12(CONSTANTS[9])
  • Variable kb3(CONSTANTS[10])
  • Variable kf5(CONSTANTS[11])
  • Variable kb5(CONSTANTS[12])
  • Variable kf7(CONSTANTS[13])
  • Variable kb7(CONSTANTS[14])
  • Variable kf10(CONSTANTS[15])
  • Variable kb10(CONSTANTS[16])
  • Variable kf11(CONSTANTS[17])
  • Variable DAG(RATES[0] and VARIABLES[0])
  • Variable PIP2CaPLA2(RATES[1] and VARIABLES[1])
  • Variable PIP2PLA2(RATES[2] and VARIABLES[2])
  • Variable DAGCaPLA2(RATES[3] and VARIABLES[3])
  • Variable CaPLA2_star(RATES[4] and VARIABLES[4])
  • Variable kf13(CONSTANTS[18])
  • Variable kb13(CONSTANTS[19])
  • Variable km2(CONSTANTS[20])
  • Variable km4(CONSTANTS[21])
  • Variable km6(CONSTANTS[22])
  • Variable AA(RATES[5] and VARIABLES[5])
  • Variable APC(RATES[6] and VARIABLES[6])
  • Variable kb11(CONSTANTS[23])
  • Variable km9(CONSTANTS[24])
  • Variable Vmax9(CONSTANTS[25])
  • Variable PIP2_star(RATES[7] and VARIABLES[7])
  • Variable PLA2_cyt(RATES[8] and VARIABLES[8])
  • Variable MAPK(CONSTANTS[26])
  • Variable PLA2_star(RATES[9] and VARIABLES[9])
  • Variable Ca(RATES[10] and VARIABLES[10])
  • Variable CaPLA2(RATES[11] and VARIABLES[11])
• No equation needed for Newton-Raphson evaluation.

void SetupFixedConstants(double* CONSTANTS)
{
CONSTANTS[0] = 2e-09;
CONSTANTS[1] = 0.5;
CONSTANTS[2] = 1.6667e-07;
CONSTANTS[3] = 20;
CONSTANTS[4] = 20;
CONSTANTS[5] = 11.04;
CONSTANTS[6] = 5.4;
CONSTANTS[7] = 36;
CONSTANTS[8] = 60;
CONSTANTS[9] = 120;
CONSTANTS[10] = 0.1;
CONSTANTS[11] = 2e-08;
CONSTANTS[12] = 0.1;
CONSTANTS[13] = 5e-09;
CONSTANTS[14] = 4;
CONSTANTS[15] = 0.17;
CONSTANTS[16] = 0;
CONSTANTS[17] = 1e-05;
CONSTANTS[18] = 0.4;
CONSTANTS[19] = 0;
CONSTANTS[20] = 20;
CONSTANTS[21] = 20;
CONSTANTS[22] = 20;
CONSTANTS[23] = 0.1;
CONSTANTS[24] = 25.6417;
CONSTANTS[25] = 20;
CONSTANTS[26] = 0.1;
}

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

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

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