cpn  [kholodenko_M_MAPk_cascade1]
{
constants:
parameter:
  double V_1 = 2.5;
  double V_2 = 0.25;
  double V_5 = 0.75;
  double V_6 = 0.75;
  double V_9 = 0.5;
  double V_10 = 0.5;
  double K_1 = 9;
  double K_1_1 = 10;
  double K_2 = 8;
  double K_3 = 15;
  double n = 1;
  double K_4 = 15;
  double K_5 = 15;
  double K_6 = 15;
  double K_7 = 15;
  double K_8 = 15;
  double K_9 = 15;
  double K_10 = 15;
  double k_3 = 0.025;
  double k_4 = 0.025;
  double k_7 = 0.025;
  double k_8 = 0.025;
  double sf = 1;

places:
continuous:
   M_MKKK = 90;
   M_MKKK_P = 10;
   M_MKK = 280;
   M_MKK_P = 10;
   M_MKK_PP = 10;
   M_MAPK = 280;
   M_MAPK_P = 10;
   M_MAPK_PP = 10;

transitions:
continuous:
  R_a_1
    : [M_MAPK_PP < 1]
    : [M_MKKK_P + 1] & [M_MKKK - 1]
    : sf*V_1*M_MKKK/((1+(M_MAPK_PP/K_1)^n)*(K_1_1+M_MKKK))
    ;
  R_a_2
    : 
    : [M_MKKK + 1] & [M_MKKK_P - 1]
    : V_2*M_MKKK_P/(K_2+M_MKKK_P)
    ;
  R_b_1
    : 
    : [M_MKKK_P + 1] & [M_MKK_P + 1] & [M_MKKK_P - 1] & [M_MKK - 1]
    : sf*k_3*M_MKKK_P*M_MKK/(K_3+M_MKK)
    ;
  R_b_2
    : 
    : [M_MKKK_P + 1] & [M_MKK_PP + 1] & [M_MKKK_P - 1] & [M_MKK_P - 1]
    : sf*k_4*M_MKKK_P*M_MKK_P/(K_4+M_MKK_P)
    ;
  R_b_4
    : 
    : [M_MKK + 1] & [M_MKK_P - 1]
    : V_6*M_MKK_P/(K_6+M_MKK_P)
    ;
  R_b_3
    : 
    : [M_MKK_P + 1] & [M_MKK_PP - 1]
    : V_5*M_MKK_PP/(K_5+M_MKK_PP)
    ;
  R_c_1
    : 
    : [M_MKK_PP + 1] & [M_MAPK_P + 1] & [M_MKK_PP - 1] & [M_MAPK - 1]
    : sf*k_7*M_MKK_PP*M_MAPK/(K_7+M_MAPK)
    ;
  R_c_2
    : 
    : [M_MKK_PP + 1] & [M_MAPK_PP + 1] & [M_MKK_PP - 1] & [M_MAPK_P - 1]
    : sf*k_8*M_MKK_PP*M_MAPK_P/(K_8+M_MAPK_P)
    ;
  R_c_4
    : 
    : [M_MAPK + 1] & [M_MAPK_P - 1]
    : V_10*M_MAPK_P/(K_10+M_MAPK_P)
    ;
  R_c_3
    : 
    : [M_MAPK_P + 1] & [M_MAPK_PP - 1]
    : V_9*M_MAPK_PP/(K_9+M_MAPK_PP)
    ;
}