// 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
// v1: alpha 
// v2: Ste2 => Ste2a 
// v3: Ste2a => Ste2 
// v4: Ste2a => p 
// v5: Ste2 
// v6: Gabc => GaGTP Gbc 
// v7: GaGTP => GaGDP 
// v8: GaGTP => GaGDP 
// v9: GaGDP Gbc => Gabc 
// v10: Gbc complexC => complexD 
// v11: complexD => Gbc complexC 
// v12: Ste11 Ste5 => complexA 
// v13: complexA => Ste11 Ste5 
// v14: Fus3 Ste7 => complexB 
// v15: complexB => Fus3 Ste7 
// v16: complexA complexB => complexC 
// v17: complexC => Fus3 Ste11 Ste7 Ste5 
// v18: complexD Ste20 => complexE 
// v19: complexE => complexD Ste20 
// v20: complexE => complexF 
// v21: complexE => Gbc Ste7 Ste5 Fus3 Ste20 Ste11 
// v22: complexF => complexG 
// v23: complexF => Gbc Ste7 Ste5 Fus3 Ste20 Ste11 
// v24: complexG => complexH 
// v25: complexG => Gbc Ste7 Ste5 Fus3 Ste20 Ste11 
// v26: complexH => complexI 
// v27: complexH => Gbc Ste7 Ste5 Fus3 Ste20 Ste11 
// v28: complexI => complexL Fus3PP 
// v29: complexL Fus3 => complexK 
// v30: complexK => complexL Fus3 
// v31: complexK => complexI 
// v32: complexL => Gbc Ste7 Ste5 Ste20 Ste11 
// v33: Fus3PP => Fus3 
// v34: Ste12 Fus3PP => Ste12a 
// v35: Ste12a => Ste12 Fus3PP 
// v36: Bar1 => Bar1a 
// v37: Bar1a => Bar1 
// v38: Bar1a => Bar1aex 
// v39: Far1 => Far1PP 
// v40: Far1PP => Far1 
// v41: Far1 => Far1U 
// v42: Gbc Far1PP => complexM 
// v43: complexM => Gbc Far1PP 
// v44: complexN => Cdc28 Far1PP 
// v45: Cdc28 Far1PP => complexN 
// v46: p => Sst2 
// v47: Sst2 => p 

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

  // variable definitions
  real Extracellular = 1 L;
  real compartment = 1 L;
  real alpha(time) nM;
  real Ste2(time) nM;
  real Ste2a(time) nM;
  real Gabc(time) nM;
  real GaGTP(time) nM;
  real Gbc(time) nM;
  real GaGDP(time) nM;
  real complexC(time) nM;
  real complexD(time) nM;
  real Ste5(time) nM;
  real Ste11(time) nM;
  real complexA(time) nM;
  real Ste7(time) nM;
  real Fus3(time) nM;
  real complexB(time) nM;
  real Ste20(time) nM;
  real complexE(time) nM;
  real complexF(time) nM;
  real complexG(time) nM;
  real complexH(time) nM;
  real complexI(time) nM;
  real complexL(time) nM;
  real Fus3PP(time) nM;
  real complexK(time) nM;
  real Ste12(time) nM;
  real Ste12a(time) nM;
  real Bar1(time) nM;
  real Bar1a(time) nM;
  real Bar1aex(time) nM;
  real Far1(time) nM;
  real Far1PP(time) nM;
  real Far1U(time) nM;
  real complexM(time) nM;
  real complexN(time) nM;
  real Cdc28(time) nM;
  real Sst2(time) nM;
  real p = 0 nM;
  real v1(time) nmol/min;
  real k1 = .03 1/nmol*L/min;
  real v2(time) nmol/min;
  real k2 = .0012 1/nmol*L/min;
  real v3(time) nmol/min;
  real k3 = .6 1/min;
  real v4(time) nmol/min;
  real k4 = .24 1/min;
  real v5(time) nmol/min;
  real k5 = .024 1/min;
  real v6(time) nmol/min;
  real k6 = .0036 1/nmol*L/min;
  real v7(time) nmol/min;
  real k7 = .24 1/min;
  real v8(time) nmol/min;
  real k8 = .033 1/nmol*L/min;
  real v9(time) nmol/min;
  real k9 = 2E3 1/nmol*L/min;
  real v10(time) nmol/min;
  real k10 = .1 1/nmol*L/min;
  real v11(time) nmol/min;
  real k11 = 5 1/min;
  real v12(time) nmol/min;
  real k12 = 1 1/nmol*L/min;
  real v13(time) nmol/min;
  real k13 = 3 1/min;
  real v14(time) nmol/min;
  real k14 = 1 1/nmol*L/min;
  real v15(time) nmol/min;
  real k15 = 3 1/min;
  real v16(time) nmol/min;
  real k16 = 3 1/nmol*L/min;
  real v17(time) nmol/min;
  real k17 = 100 1/min;
  real v18(time) nmol/min;
  real k18 = 5 1/nmol*L/min;
  real v19(time) nmol/min;
  real k19 = 1 1/min;
  real v20(time) nmol/min;
  real k20 = 10 1/min;
  real v21(time) nmol/min;
  real k21 = 5 1/min;
  real v22(time) nmol/min;
  real k22 = 47 1/min;
  real v23(time) nmol/min;
  real k23 = 5 1/min;
  real v24(time) nmol/min;
  real k24 = 345 1/min;
  real v25(time) nmol/min;
  real k25 = 5 1/min;
  real v26(time) nmol/min;
  real k26 = 50 1/min;
  real v27(time) nmol/min;
  real k27 = 5 1/min;
  real v28(time) nmol/min;
  real k28 = 140 1/min;
  real v29(time) nmol/min;
  real k29 = 10 1/nmol*L/min;
  real v30(time) nmol/min;
  real k30 = 1 1/min;
  real v31(time) nmol/min;
  real k31 = 250 1/min;
  real v32(time) nmol/min;
  real k32 = 5 1/min;
  real v33(time) nmol/min;
  real k33 = 50 1/min;
  real v34(time) nmol/min;
  real k34 = 18 1/nmol*L/min;
  real v35(time) nmol/min;
  real k35 = 10 1/min;
  real v36(time) nmol/min;
  real k36 = .1 1/nmol*L/min;
  real v37(time) nmol/min;
  real k37 = .1 1/min;
  real v38(time) nmol/min;
  real k38 = .01 1/min;
  real v39(time) nmol/min;
  real k39 = 18 1/min;
  real v40(time) nmol/min;
  real k40 = 1 1/min;
  real v41(time) nmol/min;
  real k41 = .02 1/nmol*L/min;
  real v42(time) nmol/min;
  real k42 = .1 1/nmol*L/min;
  real v43(time) nmol/min;
  real k43 = .01 1/min;
  real v44(time) nmol/min;
  real k44 = .01 1/min;
  real v45(time) nmol/min;
  real k45 = .1 1/nmol*L/min;
  real v46(time) nmol/min;
  real k46 = 200 nmol/L/min;
  real v47(time) nmol/min;
  real k47 = 1 1/min;

  // equations
  when (time=time.min) alpha = 100;
  alpha:time = (-1*v1)/Extracellular;
  when (time=time.min) Ste2 = 1.6666667E3;
  Ste2:time = (-1*v2 + v3 + -1*v5)/compartment;
  when (time=time.min) Ste2a = 0;
  Ste2a:time = (v2 + -1*v3 + -1*v4)/compartment;
  when (time=time.min) Gabc = 1.6666667E3;
  Gabc:time = (-1*v6 + v9)/compartment;
  when (time=time.min) GaGTP = 0;
  GaGTP:time = (v6 + -1*v7 + -1*v8)/compartment;
  when (time=time.min) Gbc = 0;
  Gbc:time = (v6 + -1*v9 + -1*v10 + v11 + v21 + v23 + v25 + v27 + v32 + -1*v42 + v43)/compartment;
  when (time=time.min) GaGDP = 0;
  GaGDP:time = (v7 + v8 + -1*v9)/compartment;
  when (time=time.min) complexC = 2.3572494E2;
  complexC:time = (-1*v10 + v11 + v16 + -1*v17)/compartment;
  when (time=time.min) complexD = 0;
  complexD:time = (v10 + -1*v11 + -1*v18 + v19)/compartment;
  when (time=time.min) Ste5 = 1.5833177E2;
  Ste5:time = (-1*v12 + v13 + v17 + v21 + v23 + v25 + v27 + v32)/compartment;
  when (time=time.min) Ste11 = 1.5833177E2;
  Ste11:time = (-1*v12 + v13 + v17 + v21 + v23 + v25 + v27 + v32)/compartment;
  when (time=time.min) complexA = 1.059433E2;
  complexA:time = (v12 + -1*v13 + -1*v16)/compartment;
  when (time=time.min) Ste7 = 36.39970164;
  Ste7:time = (-1*v14 + v15 + v17 + v21 + v23 + v25 + v27 + v32)/compartment;
  when (time=time.min) Fus3 = 6.863997E2;
  Fus3:time = (-1*v14 + v15 + v17 + v21 + v23 + v25 + v27 + -1*v29 + v30 + v33)/compartment;
  when (time=time.min) complexB = 77.87536257;
  complexB:time = (v14 + -1*v15 + -1*v16)/compartment;
  when (time=time.min) Ste20 = 1E3;
  Ste20:time = (-1*v18 + v19 + v21 + v23 + v25 + v27 + v32)/compartment;
  when (time=time.min) complexE = 0;
  complexE:time = (v18 + -1*v19 + -1*v20 + -1*v21)/compartment;
  when (time=time.min) complexF = 0;
  complexF:time = (v20 + -1*v22 + -1*v23)/compartment;
  when (time=time.min) complexG = 0;
  complexG:time = (v22 + -1*v24 + -1*v25)/compartment;
  when (time=time.min) complexH = 0;
  complexH:time = (v24 + -1*v26 + -1*v27)/compartment;
  when (time=time.min) complexI = 0;
  complexI:time = (v26 + -1*v28 + v31)/compartment;
  when (time=time.min) complexL = 0;
  complexL:time = (v28 + -1*v29 + v30 + -1*v32)/compartment;
  when (time=time.min) Fus3PP = 0;
  Fus3PP:time = (v28 + -1*v33 + -1*v34 + v35)/compartment;
  when (time=time.min) complexK = 0;
  complexK:time = (v29 + -1*v30 + -1*v31)/compartment;
  when (time=time.min) Ste12 = 200;
  Ste12:time = (-1*v34 + v35)/compartment;
  when (time=time.min) Ste12a = 0;
  Ste12a:time = (v34 + -1*v35)/compartment;
  when (time=time.min) Bar1 = 200;
  Bar1:time = (-1*v36 + v37)/compartment;
  when (time=time.min) Bar1a = 0;
  Bar1a:time = (v36 + -1*v37 + -1*v38)/compartment;
  when (time=time.min) Bar1aex = 0;
  Bar1aex:time = (v38)/Extracellular;
  when (time=time.min) Far1 = 500;
  Far1:time = (-1*v39 + v40 + -1*v41)/compartment;
  when (time=time.min) Far1PP = 0;
  Far1PP:time = (v39 + -1*v40 + -1*v42 + v43 + v44 + -1*v45)/compartment;
  when (time=time.min) Far1U = 0;
  Far1U:time = (v41)/compartment;
  when (time=time.min) complexM = 0;
  complexM:time = (v42 + -1*v43)/compartment;
  when (time=time.min) complexN = 0;
  complexN:time = (-1*v44 + v45)/compartment;
  when (time=time.min) Cdc28 = 300;
  Cdc28:time = (v44 + -1*v45)/compartment;
  when (time=time.min) Sst2 = 0;
  Sst2:time = (v46 + -1*v47)/compartment;
  v1 = Extracellular*alpha*Bar1aex*k1;
  v2 = compartment*Ste2*alpha*k2;
  v3 = compartment*Ste2a*k3;
  v4 = compartment*Ste2a*k4;
  v5 = compartment*Ste2*k5;
  v6 = compartment*Ste2a*Gabc*k6;
  v7 = compartment*GaGTP*k7;
  v8 = compartment*GaGTP*Sst2*k8;
  v9 = compartment*GaGDP*Gbc*k9;
  v10 = compartment*Gbc*complexC*k10;
  v11 = compartment*complexD*k11;
  v12 = compartment*Ste5*Ste11*k12;
  v13 = compartment*complexA*k13;
  v14 = compartment*Ste7*Fus3*k14;
  v15 = compartment*complexB*k15;
  v16 = compartment*complexA*complexB*k16;
  v17 = compartment*complexC*k17;
  v18 = compartment*complexD*Ste20*k18;
  v19 = compartment*complexE*k19;
  v20 = compartment*complexE*k20;
  v21 = compartment*complexE*k21;
  v22 = compartment*complexF*k22;
  v23 = compartment*complexF*k23;
  v24 = compartment*complexG*k24;
  v25 = compartment*complexG*k25;
  v26 = compartment*complexH*k26;
  v27 = compartment*complexH*k27;
  v28 = compartment*complexI*k28;
  v29 = compartment*complexL*Fus3*k29;
  v30 = compartment*complexK*k30;
  v31 = compartment*complexK*k31;
  v32 = compartment*complexL*k32;
  v33 = compartment*Fus3PP*k33;
  v34 = compartment*Ste12*Fus3PP*k34;
  v35 = compartment*Ste12a*k35;
  v36 = compartment*Ste12a*Bar1*k36;
  v37 = compartment*Bar1a*k37;
  v38 = compartment*Bar1a*k38;
  v39 = compartment*(Far1*Fus3PP*Fus3PP/(100*100+Fus3PP*Fus3PP))*k39;
  v40 = compartment*Far1PP*k40;
  v41 = compartment*Far1*Cdc28*k41;
  v42 = compartment*Gbc*Far1PP*k42;
  v43 = compartment*complexM*k43;
  v44 = compartment*complexN*k44;
  v45 = compartment*Far1PP*Cdc28*k45;
  v46 = compartment*(Fus3PP^2/(4^2+Fus3PP^2))*k46;
  v47 = compartment*Sst2*k47;

  // variable properties
  Extracellular.sbmlRole="compartment";
  compartment.sbmlRole="compartment";
  compartment.sbmlName="cell";
  alpha.sbmlRole="species";
  alpha.sbmlCompartment="Extracellular";
  Ste2.sbmlRole="species";
  Ste2.sbmlCompartment="compartment";
  Ste2a.sbmlRole="species";
  Ste2a.sbmlCompartment="compartment";
  Gabc.sbmlRole="species";
  Gabc.sbmlCompartment="compartment";
  GaGTP.sbmlRole="species";
  GaGTP.sbmlCompartment="compartment";
  Gbc.sbmlRole="species";
  Gbc.sbmlCompartment="compartment";
  GaGDP.sbmlRole="species";
  GaGDP.sbmlCompartment="compartment";
  complexC.sbmlRole="species";
  complexC.sbmlCompartment="compartment";
  complexD.sbmlRole="species";
  complexD.sbmlCompartment="compartment";
  Ste5.sbmlRole="species";
  Ste5.sbmlCompartment="compartment";
  Ste11.sbmlRole="species";
  Ste11.sbmlCompartment="compartment";
  complexA.sbmlRole="species";
  complexA.sbmlCompartment="compartment";
  Ste7.sbmlRole="species";
  Ste7.sbmlCompartment="compartment";
  Fus3.sbmlRole="species";
  Fus3.sbmlCompartment="compartment";
  complexB.sbmlRole="species";
  complexB.sbmlCompartment="compartment";
  Ste20.sbmlRole="species";
  Ste20.sbmlCompartment="compartment";
  complexE.sbmlRole="species";
  complexE.sbmlCompartment="compartment";
  complexF.sbmlRole="species";
  complexF.sbmlCompartment="compartment";
  complexG.sbmlRole="species";
  complexG.sbmlCompartment="compartment";
  complexH.sbmlRole="species";
  complexH.sbmlCompartment="compartment";
  complexI.sbmlRole="species";
  complexI.sbmlCompartment="compartment";
  complexL.sbmlRole="species";
  complexL.sbmlCompartment="compartment";
  Fus3PP.sbmlRole="species";
  Fus3PP.sbmlCompartment="compartment";
  complexK.sbmlRole="species";
  complexK.sbmlCompartment="compartment";
  Ste12.sbmlRole="species";
  Ste12.sbmlCompartment="compartment";
  Ste12a.sbmlRole="species";
  Ste12a.sbmlCompartment="compartment";
  Bar1.sbmlRole="species";
  Bar1.sbmlCompartment="compartment";
  Bar1a.sbmlRole="species";
  Bar1a.sbmlCompartment="compartment";
  Bar1aex.sbmlRole="species";
  Bar1aex.sbmlCompartment="Extracellular";
  Far1.sbmlRole="species";
  Far1.sbmlCompartment="compartment";
  Far1PP.sbmlRole="species";
  Far1PP.sbmlCompartment="compartment";
  Far1U.sbmlRole="species";
  Far1U.sbmlCompartment="compartment";
  complexM.sbmlRole="species";
  complexM.sbmlCompartment="compartment";
  complexN.sbmlRole="species";
  complexN.sbmlCompartment="compartment";
  Cdc28.sbmlRole="species";
  Cdc28.sbmlCompartment="compartment";
  Sst2.sbmlRole="species";
  Sst2.sbmlCompartment="compartment";
  p.sbmlRole="species";
  p.sbmlCompartment="compartment";
  v1.sbmlRole="rate";
  k1.sbmlRole="parameter";
  v2.sbmlRole="rate";
  k2.sbmlRole="parameter";
  v3.sbmlRole="rate";
  k3.sbmlRole="parameter";
  v4.sbmlRole="rate";
  k4.sbmlRole="parameter";
  v5.sbmlRole="rate";
  k5.sbmlRole="parameter";
  v6.sbmlRole="rate";
  k6.sbmlRole="parameter";
  v7.sbmlRole="rate";
  k7.sbmlRole="parameter";
  v8.sbmlRole="rate";
  k8.sbmlRole="parameter";
  v9.sbmlRole="rate";
  k9.sbmlRole="parameter";
  v10.sbmlRole="rate";
  k10.sbmlRole="parameter";
  v11.sbmlRole="rate";
  k11.sbmlRole="parameter";
  v12.sbmlRole="rate";
  k12.sbmlRole="parameter";
  v13.sbmlRole="rate";
  k13.sbmlRole="parameter";
  v14.sbmlRole="rate";
  k14.sbmlRole="parameter";
  v15.sbmlRole="rate";
  k15.sbmlRole="parameter";
  v16.sbmlRole="rate";
  k16.sbmlRole="parameter";
  v17.sbmlRole="rate";
  k17.sbmlRole="parameter";
  v18.sbmlRole="rate";
  k18.sbmlRole="parameter";
  v19.sbmlRole="rate";
  k19.sbmlRole="parameter";
  v20.sbmlRole="rate";
  k20.sbmlRole="parameter";
  v21.sbmlRole="rate";
  k21.sbmlRole="parameter";
  v22.sbmlRole="rate";
  k22.sbmlRole="parameter";
  v23.sbmlRole="rate";
  k23.sbmlRole="parameter";
  v24.sbmlRole="rate";
  k24.sbmlRole="parameter";
  v25.sbmlRole="rate";
  k25.sbmlRole="parameter";
  v26.sbmlRole="rate";
  k26.sbmlRole="parameter";
  v27.sbmlRole="rate";
  k27.sbmlRole="parameter";
  v28.sbmlRole="rate";
  k28.sbmlRole="parameter";
  v29.sbmlRole="rate";
  k29.sbmlRole="parameter";
  v30.sbmlRole="rate";
  k30.sbmlRole="parameter";
  v31.sbmlRole="rate";
  k31.sbmlRole="parameter";
  v32.sbmlRole="rate";
  k32.sbmlRole="parameter";
  v33.sbmlRole="rate";
  k33.sbmlRole="parameter";
  v34.sbmlRole="rate";
  k34.sbmlRole="parameter";
  v35.sbmlRole="rate";
  k35.sbmlRole="parameter";
  v36.sbmlRole="rate";
  k36.sbmlRole="parameter";
  v37.sbmlRole="rate";
  k37.sbmlRole="parameter";
  v38.sbmlRole="rate";
  k38.sbmlRole="parameter";
  v39.sbmlRole="rate";
  k39.sbmlRole="parameter";
  v40.sbmlRole="rate";
  k40.sbmlRole="parameter";
  v41.sbmlRole="rate";
  k41.sbmlRole="parameter";
  v42.sbmlRole="rate";
  k42.sbmlRole="parameter";
  v43.sbmlRole="rate";
  k43.sbmlRole="parameter";
  v44.sbmlRole="rate";
  k44.sbmlRole="parameter";
  v45.sbmlRole="rate";
  k45.sbmlRole="parameter";
  v46.sbmlRole="rate";
  k46.sbmlRole="parameter";
  v47.sbmlRole="rate";
  k47.sbmlRole="parameter";
}