// 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
// Reaction1: RAF RAFK <=> RAFRAFK 
// Reaction2: RAFRAFK <=> RAF RAFK 
// Reaction3: RAFRAFK <=> RAFK RAFp 
// Reaction4: RAFp RAFPH <=> RAFpRAFPH 
// Reaction5: RAFpRAFPH <=> RAFp RAFPH 
// Reaction6: RAFpRAFPH <=> RAF RAFPH 
// Reaction7: MEK RAFp <=> MEKRAFp 
// Reaction8: MEKRAFp <=> MEK RAFp 
// Reaction9: MEKRAFp <=> MEKp RAFp 
// Reaction10: MEKp MEKPH <=> MEKpMEKPH 
// Reaction11: MEKpMEKPH <=> MEKp MEKPH 
// Reaction12: MEKpMEKPH <=> MEK MEKPH 
// Reaction13: MEKp RAFp <=> MEKpRAFp 
// Reaction14: MEKpRAFp <=> MEKp RAFp 
// Reaction15: MEKpRAFp <=> MEKpp RAFp 
// Reaction16: MEKPH MEKpp <=> MEKppMEKPH 
// Reaction17: MEKppMEKPH <=> MEKPH MEKpp 
// Reaction18: MEKppMEKPH <=> MEKp MEKPH 
// Reaction19: MAPK MEKpp <=> MAPKMEKpp 
// Reaction20: MAPKMEKpp <=> MAPK MEKpp 
// Reaction21: MAPKMEKpp <=> MAPKp MEKpp 
// Reaction22: MAPKp MAPKPH <=> MAPKpMAPKPH 
// Reaction23: MAPKpMAPKPH <=> MAPKp MAPKPH 
// Reaction24: MAPKpMAPKPH <=> MAPK MAPKPH 
// Reaction25: MAPKp MEKpp <=> MAPKpMEKpp 
// Reaction26: MAPKpMEKpp <=> MAPKp MEKpp 
// Reaction27: MAPKpMEKpp <=> MAPKpp MEKpp 
// Reaction28: MAPKPH MAPKpp <=> MAPKppMAPKPH 
// Reaction29: MAPKppMAPKPH <=> MAPKPH MAPKpp 
// Reaction30: MAPKppMAPKPH <=> MAPKp MAPKPH 

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

  // variable definitions
  real Cytoplasm = 1 L;
  real MAPK(time) umol;
  real MAPKMEKpp(time) umol;
  real MAPKp(time) umol;
  real MAPKPH(time) umol;
  real MAPKpMAPKPH(time) umol;
  real MAPKpMEKpp(time) umol;
  real MAPKpp(time) umol;
  real MAPKppMAPKPH(time) umol;
  real MEK(time) umol;
  real MEKp(time) umol;
  real MEKPH(time) umol;
  real MEKpMEKPH(time) umol;
  real MEKpp(time) umol;
  real MEKppMEKPH(time) umol;
  real MEKpRAFp(time) umol;
  real MEKRAFp(time) umol;
  real RAF(time) umol;
  real RAFK(time) umol;
  real RAFp(time) umol;
  real RAFPH(time) umol;
  real RAFpRAFPH(time) umol;
  real RAFRAFK(time) umol;
  real Reaction1(time) microkatal;
  real a1 = 1;
  real Reaction2(time) microkatal;
  real d1 = .4;
  real Reaction3(time) microkatal;
  real k1 = .1;
  real Reaction4(time) microkatal;
  real a2 = .5;
  real Reaction5(time) microkatal;
  real d2 = .5;
  real Reaction6(time) microkatal;
  real k2 = .1;
  real Reaction7(time) microkatal;
  real a3 = 3.3;
  real Reaction8(time) microkatal;
  real d3 = .42;
  real Reaction9(time) microkatal;
  real k3 = .1;
  real Reaction10(time) microkatal;
  real a4 = 10;
  real Reaction11(time) microkatal;
  real d4 = .8;
  real Reaction12(time) microkatal;
  real k4 = .1;
  real Reaction13(time) microkatal;
  real a5 = 3.3;
  real Reaction14(time) microkatal;
  real d5 = .4;
  real Reaction15(time) microkatal;
  real k5 = .1;
  real Reaction16(time) microkatal;
  real a6 = 10;
  real Reaction17(time) microkatal;
  real d6 = .8;
  real Reaction18(time) microkatal;
  real k6 = .1;
  real Reaction19(time) microkatal;
  real a7 = 20;
  real Reaction20(time) microkatal;
  real d7 = .6;
  real Reaction21(time) microkatal;
  real k7 = .1;
  real Reaction22(time) microkatal;
  real a8 = 5;
  real Reaction23(time) microkatal;
  real d8 = .4;
  real Reaction24(time) microkatal;
  real k8 = .1;
  real Reaction25(time) microkatal;
  real a9 = 20;
  real Reaction26(time) microkatal;
  real d9 = .6;
  real Reaction27(time) microkatal;
  real k9 = .1;
  real Reaction28(time) microkatal;
  real a10 = 5;
  real Reaction29(time) microkatal;
  real d10 = .4;
  real Reaction30(time) microkatal;
  real k10 = .1;

  // equations
  when (time=time.min) MAPK = .4;
  MAPK:time = -1*Reaction19 + Reaction20 + Reaction24;
  when (time=time.min) MAPKMEKpp = 0;
  MAPKMEKpp:time = Reaction19 + -1*Reaction20 + -1*Reaction21;
  when (time=time.min) MAPKp = 0;
  MAPKp:time = Reaction21 + -1*Reaction22 + Reaction23 + -1*Reaction25 + Reaction26 + Reaction30;
  when (time=time.min) MAPKPH = .3;
  MAPKPH:time = -1*Reaction22 + Reaction23 + Reaction24 + -1*Reaction28 + Reaction29 + Reaction30;
  when (time=time.min) MAPKpMAPKPH = 0;
  MAPKpMAPKPH:time = Reaction22 + -1*Reaction23 + -1*Reaction24;
  when (time=time.min) MAPKpMEKpp = 0;
  MAPKpMEKpp:time = Reaction25 + -1*Reaction26 + -1*Reaction27;
  when (time=time.min) MAPKpp = 0;
  MAPKpp:time = Reaction27 + -1*Reaction28 + Reaction29;
  when (time=time.min) MAPKppMAPKPH = 0;
  MAPKppMAPKPH:time = Reaction28 + -1*Reaction29 + -1*Reaction30;
  when (time=time.min) MEK = .2;
  MEK:time = -1*Reaction7 + Reaction8 + Reaction12;
  when (time=time.min) MEKp = 0;
  MEKp:time = Reaction9 + -1*Reaction10 + Reaction11 + -1*Reaction13 + Reaction14 + Reaction18;
  when (time=time.min) MEKPH = .2;
  MEKPH:time = -1*Reaction10 + Reaction11 + Reaction12 + -1*Reaction16 + Reaction17 + Reaction18;
  when (time=time.min) MEKpMEKPH = 0;
  MEKpMEKPH:time = Reaction10 + -1*Reaction11 + -1*Reaction12;
  when (time=time.min) MEKpp = 0;
  MEKpp:time = Reaction15 + -1*Reaction16 + Reaction17 + -1*Reaction19 + Reaction20 + Reaction21 + -1*Reaction25 + Reaction26 + Reaction27;
  when (time=time.min) MEKppMEKPH = 0;
  MEKppMEKPH:time = Reaction16 + -1*Reaction17 + -1*Reaction18;
  when (time=time.min) MEKpRAFp = 0;
  MEKpRAFp:time = Reaction13 + -1*Reaction14 + -1*Reaction15;
  when (time=time.min) MEKRAFp = 0;
  MEKRAFp:time = Reaction7 + -1*Reaction8 + -1*Reaction9;
  when (time=time.min) RAF = .3;
  RAF:time = -1*Reaction1 + Reaction2 + Reaction6;
  when (time=time.min) RAFK = .2;
  RAFK:time = -1*Reaction1 + Reaction2 + Reaction3;
  when (time=time.min) RAFp = 0;
  RAFp:time = Reaction3 + -1*Reaction4 + Reaction5 + -1*Reaction7 + Reaction8 + Reaction9 + -1*Reaction13 + Reaction14 + Reaction15;
  when (time=time.min) RAFPH = .3;
  RAFPH:time = -1*Reaction4 + Reaction5 + Reaction6;
  when (time=time.min) RAFpRAFPH = 0;
  RAFpRAFPH:time = Reaction4 + -1*Reaction5 + -1*Reaction6;
  when (time=time.min) RAFRAFK = 0;
  RAFRAFK:time = Reaction1 + -1*Reaction2 + -1*Reaction3;
  Reaction1 = a1*RAF*RAFK;
  Reaction2 = d1*RAFRAFK;
  Reaction3 = k1*RAFRAFK;
  Reaction4 = a2*RAFp*RAFPH;
  Reaction5 = d2*RAFpRAFPH;
  Reaction6 = k2*RAFpRAFPH;
  Reaction7 = a3*MEK*RAFp;
  Reaction8 = d3*MEKRAFp;
  Reaction9 = k3*MEKRAFp;
  Reaction10 = a4*MEKp*MEKPH;
  Reaction11 = d4*MEKpMEKPH;
  Reaction12 = k4*MEKpMEKPH;
  Reaction13 = a5*MEKp*RAFp;
  Reaction14 = d5*MEKpRAFp;
  Reaction15 = k5*MEKpRAFp;
  Reaction16 = a6*MEKPH*MEKpp;
  Reaction17 = d6*MEKppMEKPH;
  Reaction18 = k6*MEKppMEKPH;
  Reaction19 = a7*MAPK*MEKpp;
  Reaction20 = d7*MAPKMEKpp;
  Reaction21 = k7*MAPKMEKpp;
  Reaction22 = a8*MAPKp*MAPKPH;
  Reaction23 = d8*MAPKpMAPKPH;
  Reaction24 = k8*MAPKpMAPKPH;
  Reaction25 = a9*MAPKp*MEKpp;
  Reaction26 = d9*MAPKpMEKpp;
  Reaction27 = k9*MAPKpMEKpp;
  Reaction28 = a10*MAPKPH*MAPKpp;
  Reaction29 = d10*MAPKppMAPKPH;
  Reaction30 = k10*MAPKppMAPKPH;

  // variable properties
  Cytoplasm.sbmlRole="compartment";
  MAPK.sbmlRole="species";
  MAPK.sbmlCompartment="Cytoplasm";
  MAPKMEKpp.sbmlRole="species";
  MAPKMEKpp.sbmlCompartment="Cytoplasm";
  MAPKp.sbmlRole="species";
  MAPKp.sbmlCompartment="Cytoplasm";
  MAPKPH.sbmlRole="species";
  MAPKPH.sbmlCompartment="Cytoplasm";
  MAPKpMAPKPH.sbmlRole="species";
  MAPKpMAPKPH.sbmlCompartment="Cytoplasm";
  MAPKpMEKpp.sbmlRole="species";
  MAPKpMEKpp.sbmlCompartment="Cytoplasm";
  MAPKpp.sbmlRole="species";
  MAPKpp.sbmlCompartment="Cytoplasm";
  MAPKppMAPKPH.sbmlRole="species";
  MAPKppMAPKPH.sbmlCompartment="Cytoplasm";
  MEK.sbmlRole="species";
  MEK.sbmlCompartment="Cytoplasm";
  MEKp.sbmlRole="species";
  MEKp.sbmlCompartment="Cytoplasm";
  MEKPH.sbmlRole="species";
  MEKPH.sbmlCompartment="Cytoplasm";
  MEKpMEKPH.sbmlRole="species";
  MEKpMEKPH.sbmlCompartment="Cytoplasm";
  MEKpp.sbmlRole="species";
  MEKpp.sbmlCompartment="Cytoplasm";
  MEKppMEKPH.sbmlRole="species";
  MEKppMEKPH.sbmlCompartment="Cytoplasm";
  MEKpRAFp.sbmlRole="species";
  MEKpRAFp.sbmlCompartment="Cytoplasm";
  MEKRAFp.sbmlRole="species";
  MEKRAFp.sbmlCompartment="Cytoplasm";
  RAF.sbmlRole="species";
  RAF.sbmlCompartment="Cytoplasm";
  RAFK.sbmlRole="species";
  RAFK.sbmlCompartment="Cytoplasm";
  RAFp.sbmlRole="species";
  RAFp.sbmlCompartment="Cytoplasm";
  RAFPH.sbmlRole="species";
  RAFPH.sbmlCompartment="Cytoplasm";
  RAFpRAFPH.sbmlRole="species";
  RAFpRAFPH.sbmlCompartment="Cytoplasm";
  RAFRAFK.sbmlRole="species";
  RAFRAFK.sbmlCompartment="Cytoplasm";
  Reaction1.sbmlRole="rate";
  a1.sbmlRole="parameter";
  Reaction2.sbmlRole="rate";
  d1.sbmlRole="parameter";
  Reaction3.sbmlRole="rate";
  k1.sbmlRole="parameter";
  Reaction4.sbmlRole="rate";
  a2.sbmlRole="parameter";
  Reaction5.sbmlRole="rate";
  d2.sbmlRole="parameter";
  Reaction6.sbmlRole="rate";
  k2.sbmlRole="parameter";
  Reaction7.sbmlRole="rate";
  a3.sbmlRole="parameter";
  Reaction8.sbmlRole="rate";
  d3.sbmlRole="parameter";
  Reaction9.sbmlRole="rate";
  k3.sbmlRole="parameter";
  Reaction10.sbmlRole="rate";
  a4.sbmlRole="parameter";
  Reaction11.sbmlRole="rate";
  d4.sbmlRole="parameter";
  Reaction12.sbmlRole="rate";
  k4.sbmlRole="parameter";
  Reaction13.sbmlRole="rate";
  a5.sbmlRole="parameter";
  Reaction14.sbmlRole="rate";
  d5.sbmlRole="parameter";
  Reaction15.sbmlRole="rate";
  k5.sbmlRole="parameter";
  Reaction16.sbmlRole="rate";
  a6.sbmlRole="parameter";
  Reaction17.sbmlRole="rate";
  d6.sbmlRole="parameter";
  Reaction18.sbmlRole="rate";
  k6.sbmlRole="parameter";
  Reaction19.sbmlRole="rate";
  a7.sbmlRole="parameter";
  Reaction20.sbmlRole="rate";
  d7.sbmlRole="parameter";
  Reaction21.sbmlRole="rate";
  k7.sbmlRole="parameter";
  Reaction22.sbmlRole="rate";
  a8.sbmlRole="parameter";
  Reaction23.sbmlRole="rate";
  d8.sbmlRole="parameter";
  Reaction24.sbmlRole="rate";
  k8.sbmlRole="parameter";
  Reaction25.sbmlRole="rate";
  a9.sbmlRole="parameter";
  Reaction26.sbmlRole="rate";
  d9.sbmlRole="parameter";
  Reaction27.sbmlRole="rate";
  k9.sbmlRole="parameter";
  Reaction28.sbmlRole="rate";
  a10.sbmlRole="parameter";
  Reaction29.sbmlRole="rate";
  d10.sbmlRole="parameter";
  Reaction30.sbmlRole="rate";
  k10.sbmlRole="parameter";
}