Role for G Protein G-beta-gamma Isoform Specificity in Synaptic Signal Processing: A Computational Study

Generated code

• The required array sizes are:
  • VARIABLES[64]
  • CONSTANTS[18]
  • BOUND[1]
  • RATES[12]
• The following variables are defined:
  • Variable time(BOUND[0])
  • Variable V(RATES[0] and VARIABLES[0])
  • Variable R(CONSTANTS[0])
  • Variable T(CONSTANTS[1])
  • Variable F(CONSTANTS[2])
  • Variable Cm(CONSTANTS[3])
  • Variable i_app(CONSTANTS[4])
  • Variable i_Na(VARIABLES[12]; recomputed after rate changes)
  • Variable O(RATES[1] and VARIABLES[1])
  • Variable C_G1(RATES[2] and VARIABLES[2])
  • Variable C_G2(RATES[3] and VARIABLES[3])
  • Variable C_G3(RATES[4] and VARIABLES[4])
  • Variable alpha(VARIABLES[13]; recomputed after rate changes)
  • Variable alpha_(VARIABLES[14]; recomputed after rate changes)
  • Variable beta(VARIABLES[15]; recomputed after rate changes)
  • Variable beta_(VARIABLES[16]; recomputed after rate changes)
  • Variable kG_plus(VARIABLES[17]; recomputed after rate changes)
  • Variable a(RATES[5] and VARIABLES[5])
  • Variable ka_plus(CONSTANTS[5])
  • Variable ka_minus(CONSTANTS[6])
  • Variable T(CONSTANTS[7])
  • Variable delta_C1_rxn0(VARIABLES[18]; recomputed after rate changes)
  • Variable delta_C2_rxn0(VARIABLES[19]; recomputed after rate changes)
  • Variable x_infinity(VARIABLES[20]; recomputed after rate changes)
  • Variable alpha_x(VARIABLES[21]; recomputed after rate changes)
  • Variable beta_x(VARIABLES[22]; recomputed after rate changes)
  • Variable i_K(VARIABLES[23]; recomputed after rate changes)
  • Variable n(RATES[6] and VARIABLES[6])
  • Variable alpha_n(VARIABLES[24]; recomputed after rate changes)
  • Variable beta_n(VARIABLES[25]; recomputed after rate changes)
  • Variable i_leak(VARIABLES[26]; recomputed after rate changes)
  • Variable R(RATES[7] and VARIABLES[7])
  • Variable kr_plus(CONSTANTS[8])
  • Variable kr_minus(CONSTANTS[9])
  • Variable Ca(VARIABLES[27]; recomputed after rate changes)
  • Variable Ca_ex(CONSTANTS[10])
  • Variable Ca_open(VARIABLES[28]; recomputed after rate changes)
  • Variable Dc(CONSTANTS[11])
  • Variable r(CONSTANTS[12])
  • Variable sigma(VARIABLES[29]; recomputed after rate changes)
  • Variable i_V(VARIABLES[30]; recomputed after rate changes)
  • Variable g_Ca(CONSTANTS[13])
  • Variable P(CONSTANTS[14])
  • Variable C1(RATES[8] and VARIABLES[8])
  • Variable C2(RATES[9] and VARIABLES[9])
  • Variable C3(RATES[10] and VARIABLES[10])
  • Variable C4(RATES[11] and VARIABLES[11])
  • Variable rate(VARIABLES[31]; recomputed after rate changes)
  • Variable delta_C2_rxn1(VARIABLES[32]; recomputed after rate changes)
  • Variable delta_C3_rxn1(VARIABLES[33]; recomputed after rate changes)
  • Variable rate(VARIABLES[34]; recomputed after rate changes)
  • Variable backrate(VARIABLES[35]; recomputed after rate changes)
  • Variable delta_C3_rxn2(VARIABLES[36]; recomputed after rate changes)
  • Variable delta_C4_rxn2(VARIABLES[37]; recomputed after rate changes)
  • Variable rate(VARIABLES[38]; recomputed after rate changes)
  • Variable backrate(VARIABLES[39]; recomputed after rate changes)
  • Variable delta_C4_rxn3(VARIABLES[40]; recomputed after rate changes)
  • Variable delta_O_rxn3(VARIABLES[41]; recomputed after rate changes)
  • Variable rate(VARIABLES[42]; recomputed after rate changes)
  • Variable backrate(VARIABLES[43]; recomputed after rate changes)
  • Variable delta_C_G1_rxn4(VARIABLES[44]; recomputed after rate changes)
  • Variable delta_C_G2_rxn4(VARIABLES[45]; recomputed after rate changes)
  • Variable rate(VARIABLES[46]; recomputed after rate changes)
  • Variable backrate(VARIABLES[47]; recomputed after rate changes)
  • Variable delta_C_G2_rxn5(VARIABLES[48]; recomputed after rate changes)
  • Variable delta_C_G3_rxn5(VARIABLES[49]; recomputed after rate changes)
  • Variable rate(VARIABLES[50]; recomputed after rate changes)
  • Variable backrate(VARIABLES[51]; recomputed after rate changes)
  • Variable delta_C1_rxn6(VARIABLES[52]; recomputed after rate changes)
  • Variable delta_C_G1_rxn6(VARIABLES[53]; recomputed after rate changes)
  • Variable kG_minus(CONSTANTS[15])
  • Variable rate(VARIABLES[54]; recomputed after rate changes)
  • Variable backrate(VARIABLES[55]; recomputed after rate changes)
  • Variable delta_C2_rxn7(VARIABLES[56]; recomputed after rate changes)
  • Variable delta_C_G2_rxn7(VARIABLES[57]; recomputed after rate changes)
  • Variable kG2_minus(CONSTANTS[16])
  • Variable rate(VARIABLES[58]; recomputed after rate changes)
  • Variable backrate(VARIABLES[59]; recomputed after rate changes)
  • Variable delta_C3_rxn8(VARIABLES[60]; recomputed after rate changes)
  • Variable delta_C_G3_rxn8(VARIABLES[61]; recomputed after rate changes)
  • Variable kG3_minus(CONSTANTS[17])
  • Variable rate(VARIABLES[62]; recomputed after rate changes)
  • Variable backrate(VARIABLES[63]; recomputed after rate changes)
• No equation needed for Newton-Raphson evaluation.

void SetupFixedConstants(double* CONSTANTS)
{
CONSTANTS[0] = 8314.41;
CONSTANTS[1] = 310;
CONSTANTS[2] = 96485;
CONSTANTS[3] = 1;
CONSTANTS[4] = 40;
CONSTANTS[5] = 200;
CONSTANTS[6] = 0.0015;
CONSTANTS[7] = 1;
CONSTANTS[8] = 0.15;
CONSTANTS[9] = 2.5;
CONSTANTS[10] = 2;
CONSTANTS[11] = 220;
CONSTANTS[12] = 10;
CONSTANTS[13] = 1.2;
CONSTANTS[14] = 6;
CONSTANTS[15] = 0.00025;
CONSTANTS[16] = 0.01;
CONSTANTS[17] = 0.0005;
}

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

void ComputeRates(double* BOUND, double* RATES, double* CONSTANTS, double* VARIABLES)
{
RATES[0] = ((-(VARIABLES[12]+VARIABLES[23]+VARIABLES[26]+CONSTANTS[4]))/CONSTANTS[3]);
RATES[6] = (((VARIABLES[24]*(1.00000-VARIABLES[6]))-(VARIABLES[25]*VARIABLES[6]))) * 0.00100000;
RATES[7] = ((CONSTANTS[8]*VARIABLES[27]*(1.00000-VARIABLES[7]))-(CONSTANTS[9]*VARIABLES[7]));
RATES[8] = (VARIABLES[18]+VARIABLES[52]);
RATES[9] = (VARIABLES[19]+VARIABLES[32]+VARIABLES[56]);
RATES[10] = (VARIABLES[33]+VARIABLES[36]+VARIABLES[60]);
RATES[11] = (VARIABLES[37]+VARIABLES[40]);
RATES[1] = VARIABLES[41];
RATES[2] = (VARIABLES[53]+VARIABLES[44]);
RATES[3] = (VARIABLES[45]+VARIABLES[57]+VARIABLES[48]);
RATES[4] = (VARIABLES[49]+VARIABLES[61]);
RATES[5] = ((CONSTANTS[5]*CONSTANTS[7]*(1.00000-VARIABLES[5]))-(CONSTANTS[6]*VARIABLES[5]));
}

void ComputeVariables(double* BOUND, double* RATES, double* CONSTANTS, double* VARIABLES)
{
VARIABLES[21] = ((0.200000*(VARIABLES[0]+40.0000))/(1.00000-exp(((-(VARIABLES[0]+40.0000))/10.0000))));
VARIABLES[22] = (8.00000*exp((-(VARIABLES[0]+(65.0000/18.0000)))));
VARIABLES[20] = (VARIABLES[21]/(VARIABLES[21]+VARIABLES[22]));
VARIABLES[12] = (120.000*pow(VARIABLES[20],3.00000)*(1.00000-VARIABLES[6])*(VARIABLES[0]-120.000));
VARIABLES[13] = (0.450000*exp((VARIABLES[0]/22.0000)));
VARIABLES[14] = (VARIABLES[13]/8.00000);
VARIABLES[15] = (0.0150000*exp(((-VARIABLES[0])/14.0000)));
VARIABLES[16] = (VARIABLES[15]*8.00000);
VARIABLES[17] = ((3.00000*VARIABLES[5])/(680.000+(320.000*VARIABLES[5])));
VARIABLES[31] = ((4.00000*VARIABLES[13]*VARIABLES[8])+(-(VARIABLES[15]*VARIABLES[9])));
VARIABLES[18] = VARIABLES[31];
VARIABLES[19] = VARIABLES[31];
VARIABLES[23] = (36.0000*pow(VARIABLES[6],4.00000)*(VARIABLES[0]+77.0000));
VARIABLES[24] = ((0.0200000*(VARIABLES[0]+55.0000))/(1.00000-exp(((-(VARIABLES[0]+55.0000))/10.0000))));
VARIABLES[25] = (0.250000*exp(((-(VARIABLES[0]+65.0000))/80.0000)));
VARIABLES[26] = (0.300000*(VARIABLES[0]+54.0000));
VARIABLES[30] = (CONSTANTS[13]*CONSTANTS[14]*((2.00000*CONSTANTS[2]*VARIABLES[0])/(CONSTANTS[0]*CONSTANTS[1]))*(CONSTANTS[10]/(1.00000-exp(((2.00000*CONSTANTS[2]*VARIABLES[0])/(CONSTANTS[0]*CONSTANTS[1]))))));
VARIABLES[29] = (-5.18200*VARIABLES[30]);
VARIABLES[28] = (VARIABLES[29]/(2.00000*CONSTANTS[11]*CONSTANTS[12]*3.14159265358979));
VARIABLES[27] = ((VARIABLES[1]*VARIABLES[28])+0.100000);
VARIABLES[34] = ((3.00000*VARIABLES[13]*VARIABLES[9])+(-(2.00000*VARIABLES[15]*VARIABLES[10])));
VARIABLES[32] = VARIABLES[34];
VARIABLES[35] = ((-(3.00000*VARIABLES[13]*VARIABLES[9]))+(2.00000*VARIABLES[15]*VARIABLES[10]));
VARIABLES[33] = VARIABLES[35];
VARIABLES[38] = ((2.00000*VARIABLES[13]*VARIABLES[10])+(-(3.00000*VARIABLES[15]*VARIABLES[11])));
VARIABLES[36] = VARIABLES[38];
VARIABLES[39] = ((-(2.00000*VARIABLES[13]*VARIABLES[10]))+(3.00000*VARIABLES[15]*VARIABLES[11]));
VARIABLES[37] = VARIABLES[39];
VARIABLES[42] = ((VARIABLES[13]*VARIABLES[11])+(-(4.00000*VARIABLES[15]*VARIABLES[1])));
VARIABLES[40] = VARIABLES[42];
VARIABLES[43] = ((-(VARIABLES[13]*VARIABLES[11]))+(4.00000*VARIABLES[15]*VARIABLES[1]));
VARIABLES[41] = VARIABLES[43];
VARIABLES[46] = ((4.00000*VARIABLES[14]*VARIABLES[2])+(-(VARIABLES[16]*VARIABLES[3])));
VARIABLES[44] = VARIABLES[46];
VARIABLES[47] = ((-(4.00000*VARIABLES[14]*VARIABLES[2]))+(VARIABLES[16]*VARIABLES[3]));
VARIABLES[45] = VARIABLES[47];
VARIABLES[50] = ((3.00000*VARIABLES[14]*VARIABLES[3])+(-(2.00000*VARIABLES[16]*VARIABLES[4])));
VARIABLES[48] = VARIABLES[50];
VARIABLES[51] = ((-(3.00000*VARIABLES[14]*VARIABLES[3]))+(2.00000*VARIABLES[16]*VARIABLES[4]));
VARIABLES[49] = VARIABLES[51];
VARIABLES[54] = ((VARIABLES[17]*VARIABLES[8])+(-(CONSTANTS[15]*VARIABLES[2])));
VARIABLES[52] = VARIABLES[54];
VARIABLES[55] = ((-(VARIABLES[17]*VARIABLES[8]))+(CONSTANTS[15]*VARIABLES[2]));
VARIABLES[53] = VARIABLES[55];
VARIABLES[58] = ((VARIABLES[17]*VARIABLES[9])+(-(CONSTANTS[16]*VARIABLES[3])));
VARIABLES[56] = VARIABLES[58];
VARIABLES[59] = ((-(VARIABLES[17]*VARIABLES[9]))+(CONSTANTS[16]*VARIABLES[3]));
VARIABLES[57] = VARIABLES[59];
VARIABLES[62] = ((VARIABLES[17]*VARIABLES[10])+(-(CONSTANTS[17]*VARIABLES[4])));
VARIABLES[60] = VARIABLES[62];
VARIABLES[63] = ((-(VARIABLES[17]*VARIABLES[10]))+(CONSTANTS[17]*VARIABLES[4]));
VARIABLES[61] = VARIABLES[63];
#ifndef VARIABLES_FOR_RATES_ONLY
#endif
}