cpn  [kholodenko_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:
   MKKK = 90;
   MKKK_P = 10;
   MKK = 280;
   MKK_P = 10;
   MKK_PP = 10;
   MAPK = 280;
   MAPK_P = 10;
   MAPK_PP = 10;

transitions:
continuous:
  a_1
    : [MAPK_PP < 1]
    : [MKKK_P + 1] & [MKKK - 1]
    : sf*V_1*MKKK/((1+(MAPK_PP/K_1)^n)*(K_1_1+MKKK))
    ;
  a_2
    : 
    : [MKKK + 1] & [MKKK_P - 1]
    : V_2*MKKK_P/(K_2+MKKK_P)
    ;
  b_1
    : 
    : [MKKK_P + 1] & [MKK_P + 1] & [MKKK_P - 1] & [MKK - 1]
    : sf*k_3*MKKK_P*MKK/(K_3+MKK)
    ;
  b_2
    : 
    : [MKKK_P + 1] & [MKK_PP + 1] & [MKKK_P - 1] & [MKK_P - 1]
    : sf*k_4*MKKK_P*MKK_P/(K_4+MKK_P)
    ;
  b_4
    : 
    : [MKK + 1] & [MKK_P - 1]
    : V_6*MKK_P/(K_6+MKK_P)
    ;
  b_3
    : 
    : [MKK_P + 1] & [MKK_PP - 1]
    : V_5*MKK_PP/(K_5+MKK_PP)
    ;
  c_1
    : 
    : [MKK_PP + 1] & [MAPK_P + 1] & [MKK_PP - 1] & [MAPK - 1]
    : sf*k_7*MKK_PP*MAPK/(K_7+MAPK)
    ;
  c_2
    : 
    : [MKK_PP + 1] & [MAPK_PP + 1] & [MKK_PP - 1] & [MAPK_P - 1]
    : sf*k_8*MKK_PP*MAPK_P/(K_8+MAPK_P)
    ;
  c_4
    : 
    : [MAPK + 1] & [MAPK_P - 1]
    : V_10*MAPK_P/(K_10+MAPK_P)
    ;
  c_3
    : 
    : [MAPK_P + 1] & [MAPK_PP - 1]
    : V_9*MAPK_PP/(K_9+MAPK_PP)
    ;

}