// 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: EGF R => Ra 
// v2: 2Ra => R2 
// v3: R2 => RP 
// v4: RP => R2 
// v5: RP PLCg => RPLCg 
// v6: RPLCg => RPLCgP 
// v7: RPLCgP => RP PLCgP 
// v8: PLCgP => PLCg 
// v9: RP Grb => RG 
// v10: RG SOS => RGS 
// v11: RGS => GS RP 
// v12: GS => Grb SOS 
// v13: RP Shc => RSh 
// v14: RSh => RShP 
// v15: RShP => RP ShP 
// v16: ShP => Shc 
// v17: RShP Grb => RShG 
// v18: RShG => ShG RP 
// v19: RShG SOS => RShGS 
// v20: RShGS => RP ShGS 
// v21: ShP Grb => ShG 
// v22: ShG SOS => ShGS 
// v23: ShGS => GS ShP 
// v24: RShP GS => RShGS 
// v25: PLCgP => PLCgl 

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

  // variable definitions
  real compartment = 1 L;
  real k1f = .003;
  real k1b = .06;
  real k2f = .01;
  real k2b = .1;
  real k3f = 1;
  real k3b = .01;
  real V4 = 450;
  real K4 = 50;
  real k5f = .06;
  real k5b = .2;
  real k6f = 1;
  real k6b = .05;
  real k7f = .3;
  real k7b = .006;
  real V8 = 1;
  real K8 = 100;
  real k9f = .003;
  real k9b = .05;
  real k10f = .01;
  real k10b = .06;
  real k11f = .03;
  real k11b = .0045;
  real k12f = .0015;
  real k12b = 1E-4;
  real k13f = .09;
  real k13b = .6;
  real k14f = 6;
  real k14b = .06;
  real k15f = .3;
  real k15b = 9E-4;
  real V16 = 1.7;
  real K16 = 340;
  real k17f = .003;
  real k17b = .1;
  real k18f = .3;
  real k18b = 9E-4;
  real k19f = .01;
  real k19b = .0214;
  real k20f = .12;
  real k20b = 2.4E-4;
  real k21f = .003;
  real k21b = .1;
  real k22f = .03;
  real k22b = .064;
  real k23f = .1;
  real k23b = .021;
  real k24f = .009;
  real k24b = .0429;
  real k25f = 1;
  real k25b = .03;
  real ATP = 1;
  real ADP = 1;
  real EGF(time) M;
  real R(time) M;
  real Ra(time) M;
  real R2(time) M;
  real RP(time) M;
  real PLCg(time) M;
  real RPLCg(time) M;
  real RPLCgP(time) M;
  real PLCgP(time) M;
  real Grb(time) M;
  real RG(time) M;
  real SOS(time) M;
  real RGS(time) M;
  real GS(time) M;
  real Shc(time) M;
  real RSh(time) M;
  real RShP(time) M;
  real ShP(time) M;
  real RShG(time) M;
  real ShG(time) M;
  real RShGS(time) M;
  real ShGS(time) M;
  real PLCgl(time) M;
  real v1(time) katal;
  real v2(time) katal;
  real v3(time) katal;
  real v4(time) katal;
  real v5(time) katal;
  real v6(time) katal;
  real v7(time) katal;
  real v8(time) katal;
  real v9(time) katal;
  real v10(time) katal;
  real v11(time) katal;
  real v12(time) katal;
  real v13(time) katal;
  real v14(time) katal;
  real v15(time) katal;
  real v16(time) katal;
  real v17(time) katal;
  real v18(time) katal;
  real v19(time) katal;
  real v20(time) katal;
  real v21(time) katal;
  real v22(time) katal;
  real v23(time) katal;
  real v24(time) katal;
  real v25(time) katal;

  // equations
  when (time=time.min) EGF = 680;
  EGF:time = (-1*v1)/compartment;
  when (time=time.min) R = 100;
  R:time = (-1*v1)/compartment;
  when (time=time.min) Ra = 0;
  Ra:time = (v1 + -2*v2)/compartment;
  when (time=time.min) R2 = 0;
  R2:time = (v2 + -1*v3 + v4)/compartment;
  when (time=time.min) RP = 0;
  RP:time = (v3 + -1*v4 + -1*v5 + v7 + -1*v9 + v11 + -1*v13 + v15 + v18 + v20)/compartment;
  when (time=time.min) PLCg = 105;
  PLCg:time = (-1*v5 + v8)/compartment;
  when (time=time.min) RPLCg = 0;
  RPLCg:time = (v5 + -1*v6)/compartment;
  when (time=time.min) RPLCgP = 0;
  RPLCgP:time = (v6 + -1*v7)/compartment;
  when (time=time.min) PLCgP = 0;
  PLCgP:time = (v7 + -1*v8 + -1*v25)/compartment;
  when (time=time.min) Grb = 85;
  Grb:time = (-1*v9 + v12 + -1*v17 + -1*v21)/compartment;
  when (time=time.min) RG = 0;
  RG:time = (v9 + -1*v10)/compartment;
  when (time=time.min) SOS = 34;
  SOS:time = (-1*v10 + v12 + -1*v19 + -1*v22)/compartment;
  when (time=time.min) RGS = 0;
  RGS:time = (v10 + -1*v11)/compartment;
  when (time=time.min) GS = 0;
  GS:time = (v11 + -1*v12 + v23 + -1*v24)/compartment;
  when (time=time.min) Shc = 150;
  Shc:time = (-1*v13 + v16)/compartment;
  when (time=time.min) RSh = 0;
  RSh:time = (v13 + -1*v14)/compartment;
  when (time=time.min) RShP = 0;
  RShP:time = (v14 + -1*v15 + -1*v17 + -1*v24)/compartment;
  when (time=time.min) ShP = 0;
  ShP:time = (v15 + -1*v16 + -1*v21 + v23)/compartment;
  when (time=time.min) RShG = 0;
  RShG:time = (v17 + -1*v18 + -1*v19)/compartment;
  when (time=time.min) ShG = 0;
  ShG:time = (v18 + v21 + -1*v22)/compartment;
  when (time=time.min) RShGS = 0;
  RShGS:time = (v19 + -1*v20 + v24)/compartment;
  when (time=time.min) ShGS = 0;
  ShGS:time = (v20 + v22 + -1*v23)/compartment;
  when (time=time.min) PLCgl = 0;
  PLCgl:time = (v25)/compartment;
  v1 = k1f*R*EGF-k1b*Ra;
  v2 = k2f*Ra*Ra-k2b*R2;
  v3 = k3f*R2-k3b*RP;
  v4 = V4*RP/(K4+RP);
  v5 = k5f*RP*PLCg-k5b*RPLCg;
  v6 = k6f*RPLCg*ATP-k6b*ADP*RPLCgP;
  v7 = k7f*RPLCgP-k7b*RP*PLCgP;
  v8 = V8*PLCgP/(K8+PLCgP);
  v9 = k9f*RP*Grb-k9b*RG;
  v10 = k10f*RG*SOS-k10b*RGS;
  v11 = k11f*RGS-k11b*RP*GS;
  v12 = k12f*GS-k12b*Grb*SOS;
  v13 = k13f*RP*Shc-k13b*RSh;
  v14 = k14f*RSh-k14b*RShP;
  v15 = k15f*RShP-k15b*ShP*RP;
  v16 = V16*ShP/(K16+ShP);
  v17 = k17f*RShP*Grb-k17b*RShG;
  v18 = k18f*RShG-k18b*RP*ShG;
  v19 = k19f*RShG*SOS-k19b*RShGS;
  v20 = k20f*RShGS-k20b*ShGS*RP;
  v21 = k21f*ShP*Grb-k21b*ShG;
  v22 = k22f*ShG*SOS-k22b*ShGS;
  v23 = k23f*ShGS-k23b*ShP*GS;
  v24 = k24f*RShP*GS-k24b*RShGS;
  v25 = k25f*PLCgP-k25b*PLCgl;

  // variable properties
  compartment.sbmlRole="compartment";
  k1f.sbmlRole="parameter";
  k1b.sbmlRole="parameter";
  k2f.sbmlRole="parameter";
  k2b.sbmlRole="parameter";
  k3f.sbmlRole="parameter";
  k3b.sbmlRole="parameter";
  V4.sbmlRole="parameter";
  K4.sbmlRole="parameter";
  k5f.sbmlRole="parameter";
  k5b.sbmlRole="parameter";
  k6f.sbmlRole="parameter";
  k6b.sbmlRole="parameter";
  k7f.sbmlRole="parameter";
  k7b.sbmlRole="parameter";
  V8.sbmlRole="parameter";
  K8.sbmlRole="parameter";
  k9f.sbmlRole="parameter";
  k9b.sbmlRole="parameter";
  k10f.sbmlRole="parameter";
  k10b.sbmlRole="parameter";
  k11f.sbmlRole="parameter";
  k11b.sbmlRole="parameter";
  k12f.sbmlRole="parameter";
  k12b.sbmlRole="parameter";
  k13f.sbmlRole="parameter";
  k13b.sbmlRole="parameter";
  k14f.sbmlRole="parameter";
  k14b.sbmlRole="parameter";
  k15f.sbmlRole="parameter";
  k15b.sbmlRole="parameter";
  V16.sbmlRole="parameter";
  K16.sbmlRole="parameter";
  k17f.sbmlRole="parameter";
  k17b.sbmlRole="parameter";
  k18f.sbmlRole="parameter";
  k18b.sbmlRole="parameter";
  k19f.sbmlRole="parameter";
  k19b.sbmlRole="parameter";
  k20f.sbmlRole="parameter";
  k20b.sbmlRole="parameter";
  k21f.sbmlRole="parameter";
  k21b.sbmlRole="parameter";
  k22f.sbmlRole="parameter";
  k22b.sbmlRole="parameter";
  k23f.sbmlRole="parameter";
  k23b.sbmlRole="parameter";
  k24f.sbmlRole="parameter";
  k24b.sbmlRole="parameter";
  k25f.sbmlRole="parameter";
  k25b.sbmlRole="parameter";
  ATP.sbmlRole="parameter";
  ADP.sbmlRole="parameter";
  EGF.sbmlRole="species";
  EGF.sbmlCompartment="compartment";
  R.sbmlRole="species";
  R.sbmlCompartment="compartment";
  Ra.sbmlRole="species";
  Ra.sbmlCompartment="compartment";
  R2.sbmlRole="species";
  R2.sbmlCompartment="compartment";
  RP.sbmlRole="species";
  RP.sbmlCompartment="compartment";
  PLCg.sbmlRole="species";
  PLCg.sbmlCompartment="compartment";
  RPLCg.sbmlRole="species";
  RPLCg.sbmlCompartment="compartment";
  RPLCgP.sbmlRole="species";
  RPLCgP.sbmlCompartment="compartment";
  PLCgP.sbmlRole="species";
  PLCgP.sbmlCompartment="compartment";
  Grb.sbmlRole="species";
  Grb.sbmlCompartment="compartment";
  RG.sbmlRole="species";
  RG.sbmlCompartment="compartment";
  SOS.sbmlRole="species";
  SOS.sbmlCompartment="compartment";
  RGS.sbmlRole="species";
  RGS.sbmlCompartment="compartment";
  GS.sbmlRole="species";
  GS.sbmlCompartment="compartment";
  Shc.sbmlRole="species";
  Shc.sbmlCompartment="compartment";
  RSh.sbmlRole="species";
  RSh.sbmlCompartment="compartment";
  RShP.sbmlRole="species";
  RShP.sbmlCompartment="compartment";
  ShP.sbmlRole="species";
  ShP.sbmlCompartment="compartment";
  RShG.sbmlRole="species";
  RShG.sbmlCompartment="compartment";
  ShG.sbmlRole="species";
  ShG.sbmlCompartment="compartment";
  RShGS.sbmlRole="species";
  RShGS.sbmlCompartment="compartment";
  ShGS.sbmlRole="species";
  ShGS.sbmlCompartment="compartment";
  PLCgl.sbmlRole="species";
  PLCgl.sbmlCompartment="compartment";
  v1.sbmlRole="rate";
  v2.sbmlRole="rate";
  v3.sbmlRole="rate";
  v4.sbmlRole="rate";
  v5.sbmlRole="rate";
  v6.sbmlRole="rate";
  v7.sbmlRole="rate";
  v8.sbmlRole="rate";
  v9.sbmlRole="rate";
  v10.sbmlRole="rate";
  v11.sbmlRole="rate";
  v12.sbmlRole="rate";
  v13.sbmlRole="rate";
  v14.sbmlRole="rate";
  v15.sbmlRole="rate";
  v16.sbmlRole="rate";
  v17.sbmlRole="rate";
  v18.sbmlRole="rate";
  v19.sbmlRole="rate";
  v20.sbmlRole="rate";
  v21.sbmlRole="rate";
  v22.sbmlRole="rate";
  v23.sbmlRole="rate";
  v24.sbmlRole="rate";
  v25.sbmlRole="rate";
}