// 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: s => ACA 
// v2: ACA => p 
// v3: s => PKA 
// v4: PKA => p 
// v5: s => ERK2 
// v6: ERK2 => p 
// v7: s => REGA 
// v8: REGA => p 
// v9: s => cAMP 
// v10: cAMP => p 
// v11: s => excAMP 
// v12: excAMP => p 
// v13: s => CAR1 
// v14: CAR1 => p 

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

  // variable definitions
  real compartment = 1 L;
  real k1 = 2;
  real k2 = .9;
  real k3 = 2.5;
  real k4 = 1.5;
  real k5 = .6;
  real k6 = .8;
  real k7 = 1;
  real k8 = 1.3;
  real k9 = .3;
  real k10 = .8;
  real k11 = .7;
  real k12 = 4.9;
  real k13 = 23;
  real k14 = 4.5;
  real s = 1 M;
  real ACA(time) M;
  real CAR1(time) M;
  real p = 1 M;
  real PKA(time) M;
  real cAMP(time) M;
  real ERK2(time) M;
  real REGA(time) M;
  real excAMP(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;

  // equations
  when (time=time.min) ACA = 3.39;
  ACA:time = (v1 + -1*v2)/compartment;
  when (time=time.min) CAR1 = 2.45;
  CAR1:time = (v13 + -1*v14)/compartment;
  when (time=time.min) PKA = 2.2;
  PKA:time = (v3 + -1*v4)/compartment;
  when (time=time.min) cAMP = 1.38;
  cAMP:time = (v9 + -1*v10)/compartment;
  when (time=time.min) ERK2 = 1.13;
  ERK2:time = (v5 + -1*v6)/compartment;
  when (time=time.min) REGA = 1.24;
  REGA:time = (v7 + -1*v8)/compartment;
  when (time=time.min) excAMP = .48;
  excAMP:time = (v11 + -1*v12)/compartment;
  v1 = k1*CAR1;
  v2 = k2*ACA*PKA;
  v3 = k3*cAMP;
  v4 = k4*PKA;
  v5 = k5*CAR1;
  v6 = k6*PKA*ERK2;
  v7 = k7;
  v8 = k8*ERK2*REGA;
  v9 = k9*ACA;
  v10 = k10*REGA*cAMP;
  v11 = k11*ACA;
  v12 = k12*excAMP;
  v13 = k13*excAMP;
  v14 = k14*CAR1;

  // variable properties
  compartment.sbmlRole="compartment";
  k1.sbmlRole="parameter";
  k2.sbmlRole="parameter";
  k3.sbmlRole="parameter";
  k4.sbmlRole="parameter";
  k5.sbmlRole="parameter";
  k6.sbmlRole="parameter";
  k7.sbmlRole="parameter";
  k8.sbmlRole="parameter";
  k9.sbmlRole="parameter";
  k10.sbmlRole="parameter";
  k11.sbmlRole="parameter";
  k12.sbmlRole="parameter";
  k13.sbmlRole="parameter";
  k14.sbmlRole="parameter";
  s.sbmlRole="species";
  s.sbmlCompartment="compartment";
  ACA.sbmlRole="species";
  ACA.sbmlCompartment="compartment";
  CAR1.sbmlRole="species";
  CAR1.sbmlCompartment="compartment";
  p.sbmlRole="species";
  p.sbmlCompartment="compartment";
  PKA.sbmlRole="species";
  PKA.sbmlCompartment="compartment";
  cAMP.sbmlRole="species";
  cAMP.sbmlCompartment="compartment";
  ERK2.sbmlRole="species";
  ERK2.sbmlCompartment="compartment";
  REGA.sbmlRole="species";
  REGA.sbmlCompartment="compartment";
  excAMP.sbmlRole="species";
  excAMP.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";
}