// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;

// SBML property definitions
property sbmlRole=string;
property sbmlName=string;
property sbmlCompartment=string;

// SBML reactions
// reaction_0000001: M MAPKK => M_MAPKK_Y 
// reaction_0000002: M_MAPKK_Y <=> MpY MAPKK 
// reaction_0000003: MpY MAPKK => MpY_MAPKK 
// reaction_0000004: MpY_MAPKK <=> Mpp MAPKK 
// reaction_0000005: M MAPKK => M_MAPKK_T 
// reaction_0000006: M_MAPKK_T <=> MpT MAPKK 
// reaction_0000007: MpT MAPKK => MpT_MAPKK 
// reaction_0000008: MpT_MAPKK <=> Mpp MAPKK 
// reaction_0000009: Mpp MKP => Mpp_MKP_Y 
// reaction_0000010: Mpp_MKP_Y <=> MpT_MKP_Y 
// reaction_0000011: MpT_MKP_Y => MpT MKP 
// reaction_0000013: MpT MKP => MpT_MKP_T 
// reaction_0000012: MpT_MKP_T <=> M_MKP_T 
// reaction_0000015: M_MKP_T => M MKP 
// reaction_0000017: MpY MKP => MpY_MKP_Y 
// reaction_0000014: MpY_MKP_Y <=> M_MKP_Y 
// reaction_0000019: M_MKP_Y => M MKP 
// reaction_0000020: Mpp MKP => Mpp_MKP_T 
// reaction_0000021: Mpp_MKP_T <=> MpY_MKP_T 
// reaction_0000022: MpY_MKP_T => MpY MKP 

math main {
  realDomain time second;
  time.min=0;
  extern time.max;
  extern time.delta;

  // variable definitions
  real cell = 1 L;
  real k1 = .02;
  real k_1 = 1;
  real k2 = .01;
  real k3 = .032;
  real k_3 = 1;
  real k4 = 15;
  real k5 = .02;
  real k_5 = 1;
  real k6 = .01;
  real k7 = .032;
  real k_7 = 1;
  real k8 = 15;
  real h1 = .045;
  real h_1 = 1;
  real h2 = .092;
  real h3 = 1;
  real h_3 = .01;
  real h4 = .01;
  real h_4 = 1;
  real h5 = .5;
  real h6 = .086;
  real h_6 = .0011;
  real h7 = .01;
  real h_7 = 1;
  real h8 = .47;
  real h9 = .14;
  real h_9 = .0018;
  real h10 = .045;
  real h_10 = 1;
  real h11 = .092;
  real h12 = 1;
  real h_12 = .01;
  real M(time) nM;
  real MpY(time) nM;
  real MpT(time) nM;
  real Mpp(time) nM;
  real MAPKK(time) nM;
  real MKP(time) nM;
  real MpY_MAPKK(time) nM;
  real MpT_MAPKK(time) nM;
  real M_MAPKK_Y(time) nM;
  real M_MAPKK_T(time) nM;
  real Mpp_MKP_Y(time) nM;
  real Mpp_MKP_T(time) nM;
  real MpY_MKP_Y(time) nM;
  real MpY_MKP_T(time) nM;
  real MpT_MKP_Y(time) nM;
  real MpT_MKP_T(time) nM;
  real M_MKP_T(time) nM;
  real M_MKP_Y(time) nM;
  real reaction_0000001(time) nanokatal;
  real reaction_0000002(time) nanokatal;
  real reaction_0000003(time) nanokatal;
  real reaction_0000004(time) nanokatal;
  real reaction_0000005(time) nanokatal;
  real reaction_0000006(time) nanokatal;
  real reaction_0000007(time) nanokatal;
  real reaction_0000008(time) nanokatal;
  real reaction_0000009(time) nanokatal;
  real reaction_0000010(time) nanokatal;
  real reaction_0000011(time) nanokatal;
  real reaction_0000013(time) nanokatal;
  real reaction_0000012(time) nanokatal;
  real reaction_0000015(time) nanokatal;
  real reaction_0000017(time) nanokatal;
  real reaction_0000014(time) nanokatal;
  real reaction_0000019(time) nanokatal;
  real reaction_0000020(time) nanokatal;
  real reaction_0000021(time) nanokatal;
  real reaction_0000022(time) nanokatal;

  // equations
  when (time=time.min) M = 800;
  M:time = (-1*reaction_0000001 + -1*reaction_0000005 + reaction_0000015 + reaction_0000019)/cell;
  when (time=time.min) MpY = 0;
  MpY:time = (reaction_0000002 + -1*reaction_0000003 + -1*reaction_0000017 + reaction_0000022)/cell;
  when (time=time.min) MpT = 0;
  MpT:time = (reaction_0000006 + -1*reaction_0000007 + reaction_0000011 + -1*reaction_0000013)/cell;
  when (time=time.min) Mpp = 0;
  Mpp:time = (reaction_0000004 + reaction_0000008 + -1*reaction_0000009 + -1*reaction_0000020)/cell;
  when (time=time.min) MAPKK = 180;
  MAPKK:time = (-1*reaction_0000001 + reaction_0000002 + -1*reaction_0000003 + reaction_0000004 + -1*reaction_0000005 + reaction_0000006 + -1*reaction_0000007 + reaction_0000008)/cell;
  when (time=time.min) MKP = 100;
  MKP:time = (-1*reaction_0000009 + reaction_0000011 + -1*reaction_0000013 + reaction_0000015 + -1*reaction_0000017 + reaction_0000019 + -1*reaction_0000020 + reaction_0000022)/cell;
  when (time=time.min) MpY_MAPKK = 0;
  MpY_MAPKK:time = (reaction_0000003 + -1*reaction_0000004)/cell;
  when (time=time.min) MpT_MAPKK = 0;
  MpT_MAPKK:time = (reaction_0000007 + -1*reaction_0000008)/cell;
  when (time=time.min) M_MAPKK_Y = 0;
  M_MAPKK_Y:time = (reaction_0000001 + -1*reaction_0000002)/cell;
  when (time=time.min) M_MAPKK_T = 0;
  M_MAPKK_T:time = (reaction_0000005 + -1*reaction_0000006)/cell;
  when (time=time.min) Mpp_MKP_Y = 0;
  Mpp_MKP_Y:time = (reaction_0000009 + -1*reaction_0000010)/cell;
  when (time=time.min) Mpp_MKP_T = 0;
  Mpp_MKP_T:time = (reaction_0000020 + -1*reaction_0000021)/cell;
  when (time=time.min) MpY_MKP_Y = 0;
  MpY_MKP_Y:time = (reaction_0000017 + -1*reaction_0000014)/cell;
  when (time=time.min) MpY_MKP_T = 0;
  MpY_MKP_T:time = (reaction_0000021 + -1*reaction_0000022)/cell;
  when (time=time.min) MpT_MKP_Y = 0;
  MpT_MKP_Y:time = (reaction_0000010 + -1*reaction_0000011)/cell;
  when (time=time.min) MpT_MKP_T = 0;
  MpT_MKP_T:time = (reaction_0000013 + -1*reaction_0000012)/cell;
  when (time=time.min) M_MKP_T = 0;
  M_MKP_T:time = (reaction_0000012 + -1*reaction_0000015)/cell;
  when (time=time.min) M_MKP_Y = 0;
  M_MKP_Y:time = (reaction_0000014 + -1*reaction_0000019)/cell;
  reaction_0000001 = cell*(k1*M*MAPKK-k_1*M_MAPKK_Y);
  reaction_0000002 = cell*k2*M_MAPKK_Y;
  reaction_0000003 = cell*(k3*MpY*MAPKK-k_3*MpY_MAPKK);
  reaction_0000004 = cell*k4*MpY_MAPKK;
  reaction_0000005 = cell*(k5*M*MAPKK-k_5*M_MAPKK_T);
  reaction_0000006 = cell*k6*M_MAPKK_T;
  reaction_0000007 = cell*(k7*MpT*MAPKK-k_7*MpT_MAPKK);
  reaction_0000008 = cell*k8*MpT_MAPKK;
  reaction_0000009 = cell*(h1*Mpp*MKP-h_1*Mpp_MKP_Y);
  reaction_0000010 = cell*h2*Mpp_MKP_Y;
  reaction_0000011 = cell*(h3*MpT_MKP_Y-h_3*MpT*MKP);
  reaction_0000013 = cell*(h4*MpT*MKP-h_4*MpT_MKP_T);
  reaction_0000012 = cell*h5*MpT_MKP_T;
  reaction_0000015 = cell*(h6*M_MKP_T-h_6*M*MKP);
  reaction_0000017 = cell*(h7*MpY*MKP-h_7*MpY_MKP_Y);
  reaction_0000014 = cell*h8*MpY_MKP_Y;
  reaction_0000019 = cell*(h9*M_MKP_Y-h_9*M*MKP);
  reaction_0000020 = cell*(h10*Mpp*MKP-h_10*Mpp_MKP_T);
  reaction_0000021 = cell*h11*Mpp_MKP_T;
  reaction_0000022 = cell*(h12*MpY_MKP_T-h_12*MpY*MKP);

  // variable properties
  cell.sbmlRole="compartment";
  k1.sbmlRole="parameter";
  k_1.sbmlRole="parameter";
  k2.sbmlRole="parameter";
  k3.sbmlRole="parameter";
  k_3.sbmlRole="parameter";
  k4.sbmlRole="parameter";
  k5.sbmlRole="parameter";
  k_5.sbmlRole="parameter";
  k6.sbmlRole="parameter";
  k7.sbmlRole="parameter";
  k_7.sbmlRole="parameter";
  k8.sbmlRole="parameter";
  h1.sbmlRole="parameter";
  h_1.sbmlRole="parameter";
  h2.sbmlRole="parameter";
  h3.sbmlRole="parameter";
  h_3.sbmlRole="parameter";
  h4.sbmlRole="parameter";
  h_4.sbmlRole="parameter";
  h5.sbmlRole="parameter";
  h6.sbmlRole="parameter";
  h_6.sbmlRole="parameter";
  h7.sbmlRole="parameter";
  h_7.sbmlRole="parameter";
  h8.sbmlRole="parameter";
  h9.sbmlRole="parameter";
  h_9.sbmlRole="parameter";
  h10.sbmlRole="parameter";
  h_10.sbmlRole="parameter";
  h11.sbmlRole="parameter";
  h12.sbmlRole="parameter";
  h_12.sbmlRole="parameter";
  M.sbmlRole="species";
  M.sbmlCompartment="cell";
  MpY.sbmlRole="species";
  MpY.sbmlCompartment="cell";
  MpT.sbmlRole="species";
  MpT.sbmlCompartment="cell";
  Mpp.sbmlRole="species";
  Mpp.sbmlCompartment="cell";
  MAPKK.sbmlRole="species";
  MAPKK.sbmlCompartment="cell";
  MKP.sbmlRole="species";
  MKP.sbmlCompartment="cell";
  MpY_MAPKK.sbmlRole="species";
  MpY_MAPKK.sbmlCompartment="cell";
  MpT_MAPKK.sbmlRole="species";
  MpT_MAPKK.sbmlCompartment="cell";
  M_MAPKK_Y.sbmlRole="species";
  M_MAPKK_Y.sbmlCompartment="cell";
  M_MAPKK_T.sbmlRole="species";
  M_MAPKK_T.sbmlCompartment="cell";
  Mpp_MKP_Y.sbmlRole="species";
  Mpp_MKP_Y.sbmlCompartment="cell";
  Mpp_MKP_T.sbmlRole="species";
  Mpp_MKP_T.sbmlCompartment="cell";
  MpY_MKP_Y.sbmlRole="species";
  MpY_MKP_Y.sbmlCompartment="cell";
  MpY_MKP_T.sbmlRole="species";
  MpY_MKP_T.sbmlCompartment="cell";
  MpT_MKP_Y.sbmlRole="species";
  MpT_MKP_Y.sbmlCompartment="cell";
  MpT_MKP_T.sbmlRole="species";
  MpT_MKP_T.sbmlCompartment="cell";
  M_MKP_T.sbmlRole="species";
  M_MKP_T.sbmlCompartment="cell";
  M_MKP_Y.sbmlRole="species";
  M_MKP_Y.sbmlCompartment="cell";
  reaction_0000001.sbmlRole="rate";
  reaction_0000002.sbmlRole="rate";
  reaction_0000003.sbmlRole="rate";
  reaction_0000004.sbmlRole="rate";
  reaction_0000005.sbmlRole="rate";
  reaction_0000006.sbmlRole="rate";
  reaction_0000007.sbmlRole="rate";
  reaction_0000008.sbmlRole="rate";
  reaction_0000009.sbmlRole="rate";
  reaction_0000010.sbmlRole="rate";
  reaction_0000011.sbmlRole="rate";
  reaction_0000013.sbmlRole="rate";
  reaction_0000012.sbmlRole="rate";
  reaction_0000015.sbmlRole="rate";
  reaction_0000017.sbmlRole="rate";
  reaction_0000014.sbmlRole="rate";
  reaction_0000019.sbmlRole="rate";
  reaction_0000020.sbmlRole="rate";
  reaction_0000021.sbmlRole="rate";
  reaction_0000022.sbmlRole="rate";
}