// 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 MEK => M_MEK_Y 
// reaction_0000002: M_MEK_Y <=> MpY MEK 
// reaction_0000003: MpY MEK => MpY_MEK 
// reaction_0000004: MpY_MEK <=> Mpp MEK 
// reaction_0000005: M MEK => M_MEK_T 
// reaction_0000006: M_MEK_T <=> MpT MEK 
// reaction_0000007: MpT MEK => MpT_MEK 
// reaction_0000008: MpT_MEK <=> Mpp MEK 
// reaction_0000009: Mpp MKP3 => Mpp_MKP3 
// reaction_0000010: Mpp_MKP3 <=> MpT_MKP3_Y 
// reaction_0000011: MpT_MKP3_Y => MpT MKP3 
// reaction_0000012: MpT_MKP3_T <=> M_MKP3_T 
// reaction_0000013: MpT MKP3 => MpT_MKP3_T 
// reaction_0000014: MpY_MKP3 <=> M_MKP3_Y 
// reaction_0000015: M_MKP3_T => M MKP3 
// reaction_0000017: MpY MKP3 => MpY_MKP3 
// reaction_0000019: M_MKP3_Y => M MKP3 

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

  // variable definitions
  real cell = 1 L;
  real k1 = .005;
  real k_1 = 1;
  real k2 = 1.08;
  real k3 = .025;
  real k_3 = 1;
  real k4 = .007;
  real k5 = .05;
  real k_5 = 1;
  real k6 = .008;
  real k7 = .005;
  real k_7 = 1;
  real k8 = .45;
  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 M(time) nM;
  real MpY(time) nM;
  real MpT(time) nM;
  real Mpp(time) nM;
  real MEK(time) nM;
  real MKP3(time) nM;
  real MpY_MEK(time) nM;
  real MpT_MEK(time) nM;
  real M_MEK_Y(time) nM;
  real M_MEK_T(time) nM;
  real Mpp_MKP3(time) nM;
  real MpY_MKP3(time) nM;
  real MpT_MKP3_Y(time) nM;
  real MpT_MKP3_T(time) nM;
  real M_MKP3_T(time) nM;
  real M_MKP3_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_0000012(time) nanokatal;
  real reaction_0000013(time) nanokatal;
  real reaction_0000014(time) nanokatal;
  real reaction_0000015(time) nanokatal;
  real reaction_0000017(time) nanokatal;
  real reaction_0000019(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)/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)/cell;
  when (time=time.min) MEK = 180;
  MEK: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) MKP3 = 100;
  MKP3:time = (-1*reaction_0000009 + reaction_0000011 + -1*reaction_0000013 + reaction_0000015 + -1*reaction_0000017 + reaction_0000019)/cell;
  when (time=time.min) MpY_MEK = 0;
  MpY_MEK:time = (reaction_0000003 + -1*reaction_0000004)/cell;
  when (time=time.min) MpT_MEK = 0;
  MpT_MEK:time = (reaction_0000007 + -1*reaction_0000008)/cell;
  when (time=time.min) M_MEK_Y = 0;
  M_MEK_Y:time = (reaction_0000001 + -1*reaction_0000002)/cell;
  when (time=time.min) M_MEK_T = 0;
  M_MEK_T:time = (reaction_0000005 + -1*reaction_0000006)/cell;
  when (time=time.min) Mpp_MKP3 = 0;
  Mpp_MKP3:time = (reaction_0000009 + -1*reaction_0000010)/cell;
  when (time=time.min) MpY_MKP3 = 0;
  MpY_MKP3:time = (-1*reaction_0000014 + reaction_0000017)/cell;
  when (time=time.min) MpT_MKP3_Y = 0;
  MpT_MKP3_Y:time = (reaction_0000010 + -1*reaction_0000011)/cell;
  when (time=time.min) MpT_MKP3_T = 0;
  MpT_MKP3_T:time = (-1*reaction_0000012 + reaction_0000013)/cell;
  when (time=time.min) M_MKP3_T = 0;
  M_MKP3_T:time = (reaction_0000012 + -1*reaction_0000015)/cell;
  when (time=time.min) M_MKP3_Y = 0;
  M_MKP3_Y:time = (reaction_0000014 + -1*reaction_0000019)/cell;
  reaction_0000001 = cell*(k1*M*MEK-k_1*M_MEK_Y);
  reaction_0000002 = cell*k2*M_MEK_Y;
  reaction_0000003 = cell*(k3*MpY*MEK-k_3*MpY_MEK);
  reaction_0000004 = cell*k4*MpY_MEK;
  reaction_0000005 = cell*(k5*M*MEK-k_5*M_MEK_T);
  reaction_0000006 = cell*k6*M_MEK_T;
  reaction_0000007 = cell*(k7*MpT*MEK-k_7*MpT_MEK);
  reaction_0000008 = cell*k8*MpT_MEK;
  reaction_0000009 = cell*(h1*Mpp*MKP3-h_1*Mpp_MKP3);
  reaction_0000010 = cell*h2*Mpp_MKP3;
  reaction_0000011 = cell*(h3*MpT_MKP3_Y-h_3*MpT*MKP3);
  reaction_0000012 = cell*h5*MpT_MKP3_T;
  reaction_0000013 = cell*(h4*MpT*MKP3-h_4*MpT_MKP3_T);
  reaction_0000014 = cell*h8*MpY_MKP3;
  reaction_0000015 = cell*(h6*M_MKP3_T-h_6*M*MKP3);
  reaction_0000017 = cell*(h7*MpY*MKP3-h_7*MpY_MKP3);
  reaction_0000019 = cell*(h9*M_MKP3_Y-h_9*M*MKP3);

  // 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";
  M.sbmlRole="species";
  M.sbmlCompartment="cell";
  MpY.sbmlRole="species";
  MpY.sbmlCompartment="cell";
  MpT.sbmlRole="species";
  MpT.sbmlCompartment="cell";
  Mpp.sbmlRole="species";
  Mpp.sbmlCompartment="cell";
  MEK.sbmlRole="species";
  MEK.sbmlCompartment="cell";
  MKP3.sbmlRole="species";
  MKP3.sbmlCompartment="cell";
  MpY_MEK.sbmlRole="species";
  MpY_MEK.sbmlCompartment="cell";
  MpT_MEK.sbmlRole="species";
  MpT_MEK.sbmlCompartment="cell";
  M_MEK_Y.sbmlRole="species";
  M_MEK_Y.sbmlCompartment="cell";
  M_MEK_T.sbmlRole="species";
  M_MEK_T.sbmlCompartment="cell";
  Mpp_MKP3.sbmlRole="species";
  Mpp_MKP3.sbmlCompartment="cell";
  MpY_MKP3.sbmlRole="species";
  MpY_MKP3.sbmlCompartment="cell";
  MpT_MKP3_Y.sbmlRole="species";
  MpT_MKP3_Y.sbmlCompartment="cell";
  MpT_MKP3_T.sbmlRole="species";
  MpT_MKP3_T.sbmlCompartment="cell";
  M_MKP3_T.sbmlRole="species";
  M_MKP3_T.sbmlCompartment="cell";
  M_MKP3_Y.sbmlRole="species";
  M_MKP3_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_0000012.sbmlRole="rate";
  reaction_0000013.sbmlRole="rate";
  reaction_0000014.sbmlRole="rate";
  reaction_0000015.sbmlRole="rate";
  reaction_0000017.sbmlRole="rate";
  reaction_0000019.sbmlRole="rate";
}