// 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
// P0_to_P1: P0 <=> P1 
// T0_to_T1: T0 <=> T1 
// P1_to_P0: P1 <=> P0 
// T1_to_T0: T1 <=> T0 
// P1_to_P2: P1 <=> P2 
// T1_to_T2: T1 <=> T2 
// P2_to_P1: P2 <=> P1 
// T2_to_T1: T2 <=> T1 
// P0_degradation: P0 
// T0_degradation: T0 
// P1_degradation: P1 
// T1_degradation: T1 
// P2_degradation: P2 
// T2_degradation: T2 
// PT_complex_formation: P2 T2 => CC 
// PT_complex_nucleation: CC => Cn 
// PT_complex_degradation: CC 
// PTnucl_complex_degradation: Cn 
// Mp_production: <=> Mp 
// Mt_production: <=> Mt 
// P0_production: <=> P0 
// T0_production: <=> T0 
// Mp_degradation: Mp 
// Mt_degradation: Mt 

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

  // variable definitions
  real Cell = 1 L;
  real compartment_0000002 = 1 L;
  real Pt(time);
  real Tt(time);
  real V_mT = .7;
  real V_dT = 2;
  real P0(time) nM;
  real T0(time) nM;
  real P1(time) nM;
  real T1(time) nM;
  real P2(time) nM;
  real T2(time) nM;
  real CC(time) nM;
  real Cn(time) nM;
  real Mp(time) nM;
  real Mt(time) nM;
  real P0_to_P1(time) nmol/hr;
  real K1_P = 2;
  real V_1P = 8;
  real T0_to_T1(time) nmol/hr;
  real K_1T = 2;
  real V_1T = 8;
  real P1_to_P0(time) nmol/hr;
  real K_2P = 2;
  real V_2P = 1;
  real T1_to_T0(time) nmol/hr;
  real K_2T = 2;
  real V_2T = 1;
  real P1_to_P2(time) nmol/hr;
  real K_3P = 2;
  real V_3P = 8;
  real T1_to_T2(time) nmol/hr;
  real K_3T = 2;
  real V_3T = 8;
  real P2_to_P1(time) nmol/hr;
  real K_4P = 2;
  real V_4P = 1;
  real T2_to_T1(time) nmol/hr;
  real K_4T = 2;
  real V_4T = 1;
  real P0_degradation(time) nmol/hr;
  real k_d = .01;
  real T0_degradation(time) nmol/hr;
  real P1_degradation(time) nmol/hr;
  real T1_degradation(time) nmol/hr;
  real P2_degradation(time) nmol/hr;
  real V_dP = 2;
  real K_dP = .2;
  real T2_degradation(time) nmol/hr;
  real K_dT = .2;
  real PT_complex_formation(time) nmol/hr;
  real k3 = 1.2;
  real k4 = .6;
  real PT_complex_nucleation(time) nmol/hr;
  real k1 = .6;
  real k2 = .2;
  real PT_complex_degradation(time) nmol/hr;
  real k_dC = .01;
  real PTnucl_complex_degradation(time) nmol/hr;
  real k_dN = .01;
  real Mp_production(time) nmol/hr;
  real v_sP = 1;
  real K_IP = 1;
  real n = 4;
  real Mt_production(time) nmol/hr;
  real V_sT = 1;
  real K_IT = 1;
  real P0_production(time) nmol/hr;
  real k_sP = .9;
  real T0_production(time) nmol/hr;
  real k_sT = .9;
  real Mp_degradation(time) nmol/hr;
  real V_mP = .7;
  real K_mP = .2;
  real Mt_degradation(time) nmol/hr;
  real K_mT = .2;

  // equations
  Pt = CC+Cn+P0+P1+P2;
  Tt = CC+Cn+T0+T1+T2;
  when (time=time.min) P0 = 0;
  P0:time = (-1*P0_to_P1 + P1_to_P0 + -1*P0_degradation + P0_production)/Cell;
  when (time=time.min) T0 = 0;
  T0:time = (-1*T0_to_T1 + T1_to_T0 + -1*T0_degradation + T0_production)/Cell;
  when (time=time.min) P1 = 0;
  P1:time = (P0_to_P1 + -1*P1_to_P0 + -1*P1_to_P2 + P2_to_P1 + -1*P1_degradation)/Cell;
  when (time=time.min) T1 = 0;
  T1:time = (T0_to_T1 + -1*T1_to_T0 + -1*T1_to_T2 + T2_to_T1 + -1*T1_degradation)/Cell;
  when (time=time.min) P2 = 0;
  P2:time = (P1_to_P2 + -1*P2_to_P1 + -1*P2_degradation + -1*PT_complex_formation)/Cell;
  when (time=time.min) T2 = 0;
  T2:time = (T1_to_T2 + -1*T2_to_T1 + -1*T2_degradation + -1*PT_complex_formation)/Cell;
  when (time=time.min) CC = 0;
  CC:time = (PT_complex_formation + -1*PT_complex_nucleation + -1*PT_complex_degradation)/Cell;
  when (time=time.min) Cn = 0;
  Cn:time = (PT_complex_nucleation + -1*PTnucl_complex_degradation)/compartment_0000002;
  when (time=time.min) Mp = 0;
  Mp:time = (Mp_production + -1*Mp_degradation)/Cell;
  when (time=time.min) Mt = 0;
  Mt:time = (Mt_production + -1*Mt_degradation)/Cell;
  P0_to_P1 = Cell*V_1P*P0/(K1_P+P0);
  T0_to_T1 = Cell*V_1T*T0/(K_1T+T0);
  P1_to_P0 = Cell*V_2P*P1/(K_2P+P1);
  T1_to_T0 = Cell*V_2T*T1/(K_2T+T1);
  P1_to_P2 = Cell*V_3P*P1/(K_3P+P1);
  T1_to_T2 = Cell*V_3T*T1/(K_3T+T1);
  P2_to_P1 = Cell*V_4P*P2/(K_4P+P2);
  T2_to_T1 = Cell*V_4T*T2/(K_4T+T2);
  P0_degradation = Cell*k_d*P0;
  T0_degradation = Cell*k_d*T0;
  P1_degradation = Cell*k_d*P1;
  T1_degradation = Cell*k_d*T1;
  P2_degradation = Cell*k_d*P2+Cell*V_dP*P2/(K_dP+P2);
  T2_degradation = Cell*k_d*T2+Cell*V_dT*T2/(K_dT+T2);
  PT_complex_formation = Cell*k3*P2*T2-Cell*k4*CC;
  PT_complex_nucleation = Cell*k1*CC-compartment_0000002*k2*Cn;
  PT_complex_degradation = Cell*k_dC*CC;
  PTnucl_complex_degradation = compartment_0000002*k_dN*Cn;
  Mp_production = Cell*v_sP*K_IP^n/(K_IP^n+Cn^n);
  Mt_production = Cell*V_sT*K_IT^n/(K_IT^n+Cn^n);
  P0_production = Cell*k_sP*Mp;
  T0_production = Cell*k_sT*Mt;
  Mp_degradation = Cell*k_d*Mp+Cell*V_mP*Mp/(K_mP+Mp);
  Mt_degradation = Cell*k_d*Mt+Cell*V_mT*Mt/(K_mT+Mt);

  // variable properties
  Cell.sbmlRole="compartment";
  Cell.sbmlName="cytoplasm";
  compartment_0000002.sbmlRole="compartment";
  compartment_0000002.sbmlName="nucleus";
  Pt.sbmlRole="parameter";
  Pt.sbmlName="Total Per";
  Tt.sbmlRole="parameter";
  Tt.sbmlName="Total Tim";
  V_mT.sbmlRole="parameter";
  V_dT.sbmlRole="parameter";
  P0.sbmlRole="species";
  P0.sbmlCompartment="Cell";
  T0.sbmlRole="species";
  T0.sbmlCompartment="Cell";
  P1.sbmlRole="species";
  P1.sbmlCompartment="Cell";
  T1.sbmlRole="species";
  T1.sbmlCompartment="Cell";
  P2.sbmlRole="species";
  P2.sbmlCompartment="Cell";
  T2.sbmlRole="species";
  T2.sbmlCompartment="Cell";
  CC.sbmlRole="species";
  CC.sbmlCompartment="Cell";
  Cn.sbmlRole="species";
  Cn.sbmlCompartment="compartment_0000002";
  Mp.sbmlRole="species";
  Mp.sbmlCompartment="Cell";
  Mt.sbmlRole="species";
  Mt.sbmlCompartment="Cell";
  P0_to_P1.sbmlRole="rate";
  K1_P.sbmlRole="parameter";
  V_1P.sbmlRole="parameter";
  T0_to_T1.sbmlRole="rate";
  K_1T.sbmlRole="parameter";
  V_1T.sbmlRole="parameter";
  P1_to_P0.sbmlRole="rate";
  K_2P.sbmlRole="parameter";
  V_2P.sbmlRole="parameter";
  T1_to_T0.sbmlRole="rate";
  K_2T.sbmlRole="parameter";
  V_2T.sbmlRole="parameter";
  P1_to_P2.sbmlRole="rate";
  K_3P.sbmlRole="parameter";
  V_3P.sbmlRole="parameter";
  T1_to_T2.sbmlRole="rate";
  K_3T.sbmlRole="parameter";
  V_3T.sbmlRole="parameter";
  P2_to_P1.sbmlRole="rate";
  K_4P.sbmlRole="parameter";
  V_4P.sbmlRole="parameter";
  T2_to_T1.sbmlRole="rate";
  K_4T.sbmlRole="parameter";
  V_4T.sbmlRole="parameter";
  P0_degradation.sbmlRole="rate";
  k_d.sbmlRole="parameter";
  T0_degradation.sbmlRole="rate";
  P1_degradation.sbmlRole="rate";
  T1_degradation.sbmlRole="rate";
  P2_degradation.sbmlRole="rate";
  V_dP.sbmlRole="parameter";
  K_dP.sbmlRole="parameter";
  T2_degradation.sbmlRole="rate";
  K_dT.sbmlRole="parameter";
  PT_complex_formation.sbmlRole="rate";
  k3.sbmlRole="parameter";
  k4.sbmlRole="parameter";
  PT_complex_nucleation.sbmlRole="rate";
  k1.sbmlRole="parameter";
  k2.sbmlRole="parameter";
  PT_complex_degradation.sbmlRole="rate";
  k_dC.sbmlRole="parameter";
  PTnucl_complex_degradation.sbmlRole="rate";
  k_dN.sbmlRole="parameter";
  Mp_production.sbmlRole="rate";
  v_sP.sbmlRole="parameter";
  K_IP.sbmlRole="parameter";
  n.sbmlRole="parameter";
  Mt_production.sbmlRole="rate";
  V_sT.sbmlRole="parameter";
  K_IT.sbmlRole="parameter";
  P0_production.sbmlRole="rate";
  k_sP.sbmlRole="parameter";
  T0_production.sbmlRole="rate";
  k_sT.sbmlRole="parameter";
  Mp_degradation.sbmlRole="rate";
  V_mP.sbmlRole="parameter";
  K_mP.sbmlRole="parameter";
  Mt_degradation.sbmlRole="rate";
  K_mT.sbmlRole="parameter";
}