// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;
unit item = dimensionless;
unit substance = 1 item;

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

// SBML reactions
// EGFBindingReaction: EGF freeEGFReceptor <=> boundEGFReceptor 
// EGFUnbindingReaction: boundEGFReceptor <=> EGF freeEGFReceptor 
// NGFBindingReaction: NGF freeNGFReceptor <=> boundNGFReceptor 
// NGFUnbindingReaction: boundNGFReceptor <=> NGF freeNGFReceptor 
// SosActivationByEGFReaction: SosInactive => SosActive 
// SosActivationByNGFReaction: SosInactive => SosActive 
// SosDeactivationReaction: SosActive => SosInactive 
// RasActivationReaction: RasInactive => RasActive 
// RasDeactivationReaction: RasActive => RasInactive 
// Raf1ByRasActivationReaction: Raf1Inactive => Raf1Active 
// MekbyRaf1ActivationReaction: MekInactive => MekActive 
// MekbyBRafActivationReaction: MekInactive => MekActive 
// ErkActivationReaction: ErkInactive => ErkActive 
// MekDeactivationReaction: MekActive => MekInactive 
// ErkDeactivationReaction: ErkActive => ErkInactive 
// Raf1byPPtaseDeactivationReaction: Raf1Active => Raf1Inactive 
// P90RskActivationReaction: P90RskInactive => P90RskActive 
// PI3KbyEGFRActivationReaction: PI3KInactive => PI3KActive 
// PI3KbyRasActivationReaction: PI3KInactive => PI3KActive 
// AktActivationReaction: AktInactive => AktActive 
// Raf1ByAktDeactivationReaction: Raf1Active => Raf1Inactive 
// C3GActivationReaction: C3GInactive => C3GActive 
// Rap1ActivationReaction: Rap1Inactive => Rap1Active 
// Rap1DeactivationReaction: Rap1Active => Rap1Inactive 
// BRafByRap1ActivationReaction: BRafInactive => BRafActive 
// BRafbyPPtaseDeactivationReaction: BRafActive => BRafInactive 

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

  // variable definitions
  real cell = 1 L;
  real krbEGF = 2.18503E-5;
  real kruEGF = .0121008;
  real krbNGF = 1.38209E-7;
  real kruNGF = .00723811;
  real kEGF = 694.731;
  real KmEGF = 6086070;
  real kNGF = 389.428;
  real KmNGF = 2112.66;
  real kdSos = 1611.97;
  real KmdSos = 896896;
  real kSos = 32.344;
  real KmSos = 35954.3;
  real kRasGap = 1509.36;
  real KmRasGap = 1432410;
  real kRasToRaf1 = .884096;
  real KmRasToRaf1 = 62464.6;
  real kpRaf1 = 185.759;
  real KmpRaf1 = 4768350;
  real kpBRaf = 125.089;
  real KmpBRaf = 157948;
  real kdMek = 2.83243;
  real KmdMek = 518753;
  real kpMekCytoplasmic = 9.85367;
  real KmpMekCytoplasmic = 1007340;
  real kdErk = 8.8912;
  real KmdErk = 3496490;
  real kpP90Rsk = .0213697;
  real KmpP90Rsk = 763523;
  real kPI3K = 10.6737;
  real KmPI3K = 184912;
  real kPI3KRas = .0771067;
  real KmPI3KRas = 272056;
  real kAkt = .0566279;
  real KmAkt = 653951;
  real kdRaf1ByAkt = 15.1212;
  real KmRaf1ByAkt = 119355;
  real kC3GNGF = 146.912;
  real KmC3GNGF = 12876.2;
  real kC3G = 1.40145;
  real KmC3G = 10965.6;
  real kRapGap = 27.265;
  real KmRapGap = 295990;
  real kRap1ToBRaf = 2.20995;
  real KmRap1ToBRaf = 1025460;
  real kdRaf1 = .126329;
  real KmdRaf1 = 1061.71;
  real kdBRaf = 441.287;
  real KmdBRaf = 1.08795E7;
  real EGF(time) substance/L;
  real NGF(time) substance/L;
  real freeEGFReceptor(time) substance/L;
  real boundEGFReceptor(time) substance/L;
  real freeNGFReceptor(time) substance/L;
  real boundNGFReceptor(time) substance/L;
  real SosInactive(time) substance/L;
  real SosActive(time) substance/L;
  real P90RskInactive(time) substance/L;
  real P90RskActive(time) substance/L;
  real RasInactive(time) substance/L;
  real RasActive(time) substance/L;
  real RasGapActive = 1.2E5 substance/L;
  real Raf1Inactive(time) substance/L;
  real Raf1Active(time) substance/L;
  real BRafInactive(time) substance/L;
  real BRafActive(time) substance/L;
  real MekInactive(time) substance/L;
  real MekActive(time) substance/L;
  real ErkInactive(time) substance/L;
  real ErkActive(time) substance/L;
  real PI3KInactive(time) substance/L;
  real PI3KActive(time) substance/L;
  real AktInactive(time) substance/L;
  real AktActive(time) substance/L;
  real C3GInactive(time) substance/L;
  real C3GActive(time) substance/L;
  real Rap1Inactive(time) substance/L;
  real Rap1Active(time) substance/L;
  real RapGapActive = 1.2E5 substance/L;
  real PP2AActive = 1.2E5 substance/L;
  real Raf1PPtase = 1.2E5 substance/L;
  real EGFBindingReaction(time) substance/min;
  real EGFUnbindingReaction(time) substance/min;
  real NGFBindingReaction(time) substance/min;
  real NGFUnbindingReaction(time) substance/min;
  real SosActivationByEGFReaction(time) substance/min;
  real SosActivationByNGFReaction(time) substance/min;
  real SosDeactivationReaction(time) substance/min;
  real RasActivationReaction(time) substance/min;
  real RasDeactivationReaction(time) substance/min;
  real Raf1ByRasActivationReaction(time) substance/min;
  real MekbyRaf1ActivationReaction(time) substance/min;
  real MekbyBRafActivationReaction(time) substance/min;
  real ErkActivationReaction(time) substance/min;
  real MekDeactivationReaction(time) substance/min;
  real ErkDeactivationReaction(time) substance/min;
  real Raf1byPPtaseDeactivationReaction(time) substance/min;
  real P90RskActivationReaction(time) substance/min;
  real PI3KbyEGFRActivationReaction(time) substance/min;
  real PI3KbyRasActivationReaction(time) substance/min;
  real AktActivationReaction(time) substance/min;
  real Raf1ByAktDeactivationReaction(time) substance/min;
  real C3GActivationReaction(time) substance/min;
  real Rap1ActivationReaction(time) substance/min;
  real Rap1DeactivationReaction(time) substance/min;
  real BRafByRap1ActivationReaction(time) substance/min;
  real BRafbyPPtaseDeactivationReaction(time) substance/min;

  // equations
  when (time=time.min) EGF = 1.0002E7;
  EGF:time = (-1*EGFBindingReaction + EGFUnbindingReaction)/cell;
  when (time=time.min) NGF = 456000;
  NGF:time = (-1*NGFBindingReaction + NGFUnbindingReaction)/cell;
  when (time=time.min) freeEGFReceptor = 8E4;
  freeEGFReceptor:time = (-1*EGFBindingReaction + EGFUnbindingReaction)/cell;
  when (time=time.min) boundEGFReceptor = 0;
  boundEGFReceptor:time = (EGFBindingReaction + -1*EGFUnbindingReaction)/cell;
  when (time=time.min) freeNGFReceptor = 1E4;
  freeNGFReceptor:time = (-1*NGFBindingReaction + NGFUnbindingReaction)/cell;
  when (time=time.min) boundNGFReceptor = 0;
  boundNGFReceptor:time = (NGFBindingReaction + -1*NGFUnbindingReaction)/cell;
  when (time=time.min) SosInactive = 1.2E5;
  SosInactive:time = (-1*SosActivationByEGFReaction + -1*SosActivationByNGFReaction + SosDeactivationReaction)/cell;
  when (time=time.min) SosActive = 0;
  SosActive:time = (SosActivationByEGFReaction + SosActivationByNGFReaction + -1*SosDeactivationReaction)/cell;
  when (time=time.min) P90RskInactive = 1.2E5;
  P90RskInactive:time = (-1*P90RskActivationReaction)/cell;
  when (time=time.min) P90RskActive = 0;
  P90RskActive:time = (P90RskActivationReaction)/cell;
  when (time=time.min) RasInactive = 1.2E5;
  RasInactive:time = (-1*RasActivationReaction + RasDeactivationReaction)/cell;
  when (time=time.min) RasActive = 0;
  RasActive:time = (RasActivationReaction + -1*RasDeactivationReaction)/cell;
  when (time=time.min) Raf1Inactive = 1.2E5;
  Raf1Inactive:time = (-1*Raf1ByRasActivationReaction + Raf1byPPtaseDeactivationReaction + Raf1ByAktDeactivationReaction)/cell;
  when (time=time.min) Raf1Active = 0;
  Raf1Active:time = (Raf1ByRasActivationReaction + -1*Raf1byPPtaseDeactivationReaction + -1*Raf1ByAktDeactivationReaction)/cell;
  when (time=time.min) BRafInactive = 1.2E5;
  BRafInactive:time = (-1*BRafByRap1ActivationReaction + BRafbyPPtaseDeactivationReaction)/cell;
  when (time=time.min) BRafActive = 0;
  BRafActive:time = (BRafByRap1ActivationReaction + -1*BRafbyPPtaseDeactivationReaction)/cell;
  when (time=time.min) MekInactive = 6E5;
  MekInactive:time = (-1*MekbyRaf1ActivationReaction + -1*MekbyBRafActivationReaction + MekDeactivationReaction)/cell;
  when (time=time.min) MekActive = 0;
  MekActive:time = (MekbyRaf1ActivationReaction + MekbyBRafActivationReaction + -1*MekDeactivationReaction)/cell;
  when (time=time.min) ErkInactive = 6E5;
  ErkInactive:time = (-1*ErkActivationReaction + ErkDeactivationReaction)/cell;
  when (time=time.min) ErkActive = 0;
  ErkActive:time = (ErkActivationReaction + -1*ErkDeactivationReaction)/cell;
  when (time=time.min) PI3KInactive = 1.2E5;
  PI3KInactive:time = (-1*PI3KbyEGFRActivationReaction + -1*PI3KbyRasActivationReaction)/cell;
  when (time=time.min) PI3KActive = 0;
  PI3KActive:time = (PI3KbyEGFRActivationReaction + PI3KbyRasActivationReaction)/cell;
  when (time=time.min) AktInactive = 1.2E5;
  AktInactive:time = (-1*AktActivationReaction)/cell;
  when (time=time.min) AktActive = 0;
  AktActive:time = (AktActivationReaction)/cell;
  when (time=time.min) C3GInactive = 1.2E5;
  C3GInactive:time = (-1*C3GActivationReaction)/cell;
  when (time=time.min) C3GActive = 0;
  C3GActive:time = (C3GActivationReaction)/cell;
  when (time=time.min) Rap1Inactive = 1.2E5;
  Rap1Inactive:time = (-1*Rap1ActivationReaction + Rap1DeactivationReaction)/cell;
  when (time=time.min) Rap1Active = 0;
  Rap1Active:time = (Rap1ActivationReaction + -1*Rap1DeactivationReaction)/cell;
  EGFBindingReaction = cell*krbEGF*EGF*freeEGFReceptor;
  EGFUnbindingReaction = cell*kruEGF*boundEGFReceptor;
  NGFBindingReaction = krbNGF*NGF*freeNGFReceptor*cell;
  NGFUnbindingReaction = kruNGF*boundNGFReceptor*cell;
  SosActivationByEGFReaction = cell*(kEGF*boundEGFReceptor*SosInactive/(SosInactive+KmEGF));
  SosActivationByNGFReaction = cell*(kNGF*boundNGFReceptor*SosInactive/(SosInactive+KmNGF));
  SosDeactivationReaction = cell*(kdSos*P90RskActive*SosActive/(SosActive+KmdSos));
  RasActivationReaction = cell*(kSos*SosActive*RasInactive/(RasInactive+KmSos));
  RasDeactivationReaction = cell*(kRasGap*RasGapActive*RasActive/(RasActive+KmRasGap));
  Raf1ByRasActivationReaction = cell*(kRasToRaf1*RasActive*Raf1Inactive/(Raf1Inactive+KmRasToRaf1));
  MekbyRaf1ActivationReaction = cell*(kpRaf1*Raf1Active*MekInactive/(MekInactive+KmpRaf1));
  MekbyBRafActivationReaction = cell*(kpBRaf*BRafActive*MekInactive/(MekInactive+KmpBRaf));
  ErkActivationReaction = cell*(kpMekCytoplasmic*MekActive*ErkInactive/(ErkInactive+KmpMekCytoplasmic));
  MekDeactivationReaction = cell*(kdMek*PP2AActive*MekActive/(MekActive+KmdMek));
  ErkDeactivationReaction = cell*(kdErk*PP2AActive*ErkActive/(ErkActive+KmdErk));
  Raf1byPPtaseDeactivationReaction = cell*(kdRaf1*Raf1PPtase*Raf1Active/(Raf1Active+KmdRaf1));
  P90RskActivationReaction = cell*(kpP90Rsk*ErkActive*P90RskInactive/(P90RskInactive+KmpP90Rsk));
  PI3KbyEGFRActivationReaction = cell*(kPI3K*boundEGFReceptor*PI3KInactive/(PI3KInactive+KmPI3K));
  PI3KbyRasActivationReaction = cell*(kPI3KRas*RasActive*PI3KInactive/(PI3KInactive+KmPI3KRas));
  AktActivationReaction = cell*(kAkt*PI3KActive*AktInactive/(AktInactive+KmAkt));
  Raf1ByAktDeactivationReaction = cell*(kdRaf1ByAkt*AktActive*Raf1Active/(Raf1Active+KmRaf1ByAkt));
  C3GActivationReaction = cell*(kC3GNGF*boundNGFReceptor*C3GInactive/(C3GInactive+KmC3GNGF));
  Rap1ActivationReaction = cell*(kC3G*C3GActive*Rap1Inactive/(Rap1Inactive+KmC3G));
  Rap1DeactivationReaction = cell*(kRapGap*RapGapActive*Rap1Active/(Rap1Active+KmRapGap));
  BRafByRap1ActivationReaction = cell*(kRap1ToBRaf*Rap1Active*BRafInactive/(BRafInactive+KmRap1ToBRaf));
  BRafbyPPtaseDeactivationReaction = cell*(kdBRaf*Raf1PPtase*BRafActive/(BRafActive+KmdBRaf));

  // variable properties
  cell.sbmlRole="compartment";
  krbEGF.sbmlRole="parameter";
  kruEGF.sbmlRole="parameter";
  krbNGF.sbmlRole="parameter";
  kruNGF.sbmlRole="parameter";
  kEGF.sbmlRole="parameter";
  KmEGF.sbmlRole="parameter";
  kNGF.sbmlRole="parameter";
  KmNGF.sbmlRole="parameter";
  kdSos.sbmlRole="parameter";
  KmdSos.sbmlRole="parameter";
  kSos.sbmlRole="parameter";
  KmSos.sbmlRole="parameter";
  kRasGap.sbmlRole="parameter";
  KmRasGap.sbmlRole="parameter";
  kRasToRaf1.sbmlRole="parameter";
  KmRasToRaf1.sbmlRole="parameter";
  kpRaf1.sbmlRole="parameter";
  KmpRaf1.sbmlRole="parameter";
  kpBRaf.sbmlRole="parameter";
  KmpBRaf.sbmlRole="parameter";
  kdMek.sbmlRole="parameter";
  KmdMek.sbmlRole="parameter";
  kpMekCytoplasmic.sbmlRole="parameter";
  KmpMekCytoplasmic.sbmlRole="parameter";
  kdErk.sbmlRole="parameter";
  KmdErk.sbmlRole="parameter";
  kpP90Rsk.sbmlRole="parameter";
  KmpP90Rsk.sbmlRole="parameter";
  kPI3K.sbmlRole="parameter";
  KmPI3K.sbmlRole="parameter";
  kPI3KRas.sbmlRole="parameter";
  KmPI3KRas.sbmlRole="parameter";
  kAkt.sbmlRole="parameter";
  KmAkt.sbmlRole="parameter";
  kdRaf1ByAkt.sbmlRole="parameter";
  KmRaf1ByAkt.sbmlRole="parameter";
  kC3GNGF.sbmlRole="parameter";
  KmC3GNGF.sbmlRole="parameter";
  kC3G.sbmlRole="parameter";
  KmC3G.sbmlRole="parameter";
  kRapGap.sbmlRole="parameter";
  KmRapGap.sbmlRole="parameter";
  kRap1ToBRaf.sbmlRole="parameter";
  KmRap1ToBRaf.sbmlRole="parameter";
  kdRaf1.sbmlRole="parameter";
  KmdRaf1.sbmlRole="parameter";
  kdBRaf.sbmlRole="parameter";
  KmdBRaf.sbmlRole="parameter";
  EGF.sbmlRole="species";
  EGF.sbmlCompartment="cell";
  NGF.sbmlRole="species";
  NGF.sbmlCompartment="cell";
  freeEGFReceptor.sbmlRole="species";
  freeEGFReceptor.sbmlCompartment="cell";
  boundEGFReceptor.sbmlRole="species";
  boundEGFReceptor.sbmlCompartment="cell";
  freeNGFReceptor.sbmlRole="species";
  freeNGFReceptor.sbmlCompartment="cell";
  boundNGFReceptor.sbmlRole="species";
  boundNGFReceptor.sbmlCompartment="cell";
  SosInactive.sbmlRole="species";
  SosInactive.sbmlCompartment="cell";
  SosActive.sbmlRole="species";
  SosActive.sbmlCompartment="cell";
  P90RskInactive.sbmlRole="species";
  P90RskInactive.sbmlCompartment="cell";
  P90RskActive.sbmlRole="species";
  P90RskActive.sbmlCompartment="cell";
  RasInactive.sbmlRole="species";
  RasInactive.sbmlCompartment="cell";
  RasActive.sbmlRole="species";
  RasActive.sbmlCompartment="cell";
  RasGapActive.sbmlRole="species";
  RasGapActive.sbmlCompartment="cell";
  Raf1Inactive.sbmlRole="species";
  Raf1Inactive.sbmlCompartment="cell";
  Raf1Active.sbmlRole="species";
  Raf1Active.sbmlCompartment="cell";
  BRafInactive.sbmlRole="species";
  BRafInactive.sbmlCompartment="cell";
  BRafActive.sbmlRole="species";
  BRafActive.sbmlCompartment="cell";
  MekInactive.sbmlRole="species";
  MekInactive.sbmlCompartment="cell";
  MekActive.sbmlRole="species";
  MekActive.sbmlCompartment="cell";
  ErkInactive.sbmlRole="species";
  ErkInactive.sbmlCompartment="cell";
  ErkActive.sbmlRole="species";
  ErkActive.sbmlCompartment="cell";
  PI3KInactive.sbmlRole="species";
  PI3KInactive.sbmlCompartment="cell";
  PI3KActive.sbmlRole="species";
  PI3KActive.sbmlCompartment="cell";
  AktInactive.sbmlRole="species";
  AktInactive.sbmlCompartment="cell";
  AktActive.sbmlRole="species";
  AktActive.sbmlCompartment="cell";
  C3GInactive.sbmlRole="species";
  C3GInactive.sbmlCompartment="cell";
  C3GActive.sbmlRole="species";
  C3GActive.sbmlCompartment="cell";
  Rap1Inactive.sbmlRole="species";
  Rap1Inactive.sbmlCompartment="cell";
  Rap1Active.sbmlRole="species";
  Rap1Active.sbmlCompartment="cell";
  RapGapActive.sbmlRole="species";
  RapGapActive.sbmlCompartment="cell";
  PP2AActive.sbmlRole="species";
  PP2AActive.sbmlCompartment="cell";
  Raf1PPtase.sbmlRole="species";
  Raf1PPtase.sbmlCompartment="cell";
  EGFBindingReaction.sbmlRole="rate";
  EGFUnbindingReaction.sbmlRole="rate";
  NGFBindingReaction.sbmlRole="rate";
  NGFUnbindingReaction.sbmlRole="rate";
  SosActivationByEGFReaction.sbmlRole="rate";
  SosActivationByNGFReaction.sbmlRole="rate";
  SosDeactivationReaction.sbmlRole="rate";
  RasActivationReaction.sbmlRole="rate";
  RasDeactivationReaction.sbmlRole="rate";
  Raf1ByRasActivationReaction.sbmlRole="rate";
  MekbyRaf1ActivationReaction.sbmlRole="rate";
  MekbyBRafActivationReaction.sbmlRole="rate";
  ErkActivationReaction.sbmlRole="rate";
  MekDeactivationReaction.sbmlRole="rate";
  ErkDeactivationReaction.sbmlRole="rate";
  Raf1byPPtaseDeactivationReaction.sbmlRole="rate";
  P90RskActivationReaction.sbmlRole="rate";
  PI3KbyEGFRActivationReaction.sbmlRole="rate";
  PI3KbyRasActivationReaction.sbmlRole="rate";
  AktActivationReaction.sbmlRole="rate";
  Raf1ByAktDeactivationReaction.sbmlRole="rate";
  C3GActivationReaction.sbmlRole="rate";
  Rap1ActivationReaction.sbmlRole="rate";
  Rap1DeactivationReaction.sbmlRole="rate";
  BRafByRap1ActivationReaction.sbmlRole="rate";
  BRafbyPPtaseDeactivationReaction.sbmlRole="rate";
}