// 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
// TC: EmptySet <=> M 
// TL: EmptySet <=> P 
// mRNAD: M <=> EmptySet 
// ProteinD: P <=> EmptySet 
// DBT: P <=> EmptySet 

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

  // variable definitions
  real default = 1 L;
  real CYTOPLASM = 1 L;
  real N_A = 6.0221367E23;
  real EmptySet = 0 M;
  real M(time) M;
  real P(time) M;
  real TC(time) mole/hr;
  real Vm = 1;
  real Pcrit = .1;
  real Keq = 200;
  real TL(time) mole/hr;
  real V = .5;
  real mRNAD(time) mole/hr;
  real D = .1;
  real ProteinD(time) mole/hr;
  real DBT(time) mole/hr;
  real k1 = 10;
  real k2 = .03;
  real J = .05;

  // equations
  when (time=time.min) M = 1/CYTOPLASM;
  M:time = (TC + -1*mRNAD)/CYTOPLASM;
  when (time=time.min) P = 1/CYTOPLASM;
  P:time = (TL + -1*ProteinD + -1*DBT)/CYTOPLASM;
  TC = CYTOPLASM*(Vm/(1+(P*(1-2/(1+(1+8*Keq*P)^.5))/(2*Pcrit))^2));
  TL = V*M*CYTOPLASM;
  mRNAD = D*M*CYTOPLASM;
  ProteinD = D*P*CYTOPLASM;
  DBT = CYTOPLASM*((k1*P*(2/(1+(1+8*Keq*P)^.5))+k2*P)/(J+P));

  // variable properties
  default.sbmlRole="compartment";
  CYTOPLASM.sbmlRole="compartment";
  N_A.sbmlRole="parameter";
  EmptySet.sbmlRole="species";
  EmptySet.sbmlCompartment="default";
  M.sbmlRole="species";
  M.sbmlCompartment="CYTOPLASM";
  P.sbmlRole="species";
  P.sbmlCompartment="CYTOPLASM";
  TC.sbmlRole="rate";
  Vm.sbmlRole="parameter";
  Pcrit.sbmlRole="parameter";
  Keq.sbmlRole="parameter";
  TL.sbmlRole="rate";
  V.sbmlRole="parameter";
  V.sbmlName="Vp";
  mRNAD.sbmlRole="rate";
  D.sbmlRole="parameter";
  D.sbmlName="km";
  ProteinD.sbmlRole="rate";
  DBT.sbmlRole="rate";
  k1.sbmlRole="parameter";
  k1.sbmlName="kp1";
  k2.sbmlRole="parameter";
  k2.sbmlName="kp2";
  J.sbmlRole="parameter";
  J.sbmlName="Jp";
}