// 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
// Photoreceptor_activation: Pfr <=> Pr 
// Photoreceptor_inactivation: Pr <=> Pfr 
// Transducer_activation: Xi <=> Xa 
// Transducer_inactivation: Xa <=> Xi 
// preS_formation: prepreS <=> preS 
// S_generation: preS <=> S 
// Glucose_sensor_inactivation: Ya Gluc <=> Yi 
// S_formation: <=> S 
// V_formation: <=> V 
// S_degradation: S 
// V_degradation: V 
// Photoreceptor_decay: Pr <=> Pi 

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

  // variable definitions
  real compartment = 1 L;
  real Pfr(time) M;
  real Pr(time) M;
  real Xi(time) M;
  real Xa(time) M;
  real prepreS(time) M;
  real preS(time) M;
  real Ya(time) M;
  real S(time) M;
  real Gluc = 0 M;
  real Yi(time) M;
  real V(time) M;
  real Pi(time) M;
  real Photoreceptor_activation(time) mole/hr;
  real IfrSfrPfr = .1;
  real Photoreceptor_inactivation(time) mole/hr;
  real IrSrPr = 0;
  real Transducer_activation(time) mole/hr;
  real kia = .1;
  real Transducer_inactivation(time) mole/hr;
  real kai = .8;
  real preS_formation(time) mole/hr;
  real kx = .2;
  real S_generation(time) mole/hr;
  real ky = 1;
  real Glucose_sensor_inactivation(time) mole/hr;
  real kG = .1;
  real S_formation(time) mole/hr;
  real alpha1 = 30;
  real V_formation(time) mole/hr;
  real alpha2 = 50;
  real S_degradation(time) mole/hr;
  real kd_s = 1;
  real V_degradation(time) mole/hr;
  real kd_v = 1;
  real Photoreceptor_decay(time) mole/hr;
  real kd = .1;

  // equations
  when (time=time.min) Pfr = 10/compartment;
  Pfr:time = (-1*Photoreceptor_activation + Photoreceptor_inactivation)/compartment;
  when (time=time.min) Pr = 0;
  Pr:time = (Photoreceptor_activation + -1*Photoreceptor_inactivation + -1*Photoreceptor_decay)/compartment;
  when (time=time.min) Xi = 6/compartment;
  Xi:time = (-1*Transducer_activation + Transducer_inactivation)/compartment;
  when (time=time.min) Xa = 0;
  Xa:time = (Transducer_activation + -1*Transducer_inactivation)/compartment;
  when (time=time.min) prepreS = 200/compartment;
  prepreS:time = (-1*preS_formation)/compartment;
  when (time=time.min) preS = 0;
  preS:time = (preS_formation + -1*S_generation)/compartment;
  when (time=time.min) Ya = .9/compartment;
  Ya:time = (-1*Glucose_sensor_inactivation)/compartment;
  when (time=time.min) S = 0;
  S:time = (S_generation + S_formation + -1*S_degradation)/compartment;
  when (time=time.min) Yi = 0;
  Yi:time = (Glucose_sensor_inactivation)/compartment;
  when (time=time.min) V = 30/compartment;
  V:time = (V_formation + -1*V_degradation)/compartment;
  when (time=time.min) Pi = 0;
  Pi:time = (Photoreceptor_decay)/compartment;
  Photoreceptor_activation = compartment*Pfr*IfrSfrPfr;
  Photoreceptor_inactivation = IrSrPr*Pr*compartment;
  Transducer_activation = Xi*kia*Pr*compartment;
  Transducer_inactivation = kai*Xa*compartment;
  preS_formation = prepreS*kx*Xa*compartment;
  S_generation = preS*ky*Ya*compartment;
  Glucose_sensor_inactivation = kG*Ya*Gluc*compartment;
  S_formation = compartment*(alpha1/(1+V^3));
  V_formation = compartment*(alpha2/(1+S^3));
  S_degradation = kd_s*S*compartment;
  V_degradation = compartment*V*kd_v;
  Photoreceptor_decay = compartment*kd*Pr;

  // variable properties
  compartment.sbmlRole="compartment";
  Pfr.sbmlRole="species";
  Pfr.sbmlCompartment="compartment";
  Pr.sbmlRole="species";
  Pr.sbmlCompartment="compartment";
  Xi.sbmlRole="species";
  Xi.sbmlCompartment="compartment";
  Xa.sbmlRole="species";
  Xa.sbmlCompartment="compartment";
  prepreS.sbmlRole="species";
  prepreS.sbmlCompartment="compartment";
  preS.sbmlRole="species";
  preS.sbmlCompartment="compartment";
  Ya.sbmlRole="species";
  Ya.sbmlCompartment="compartment";
  S.sbmlRole="species";
  S.sbmlCompartment="compartment";
  Gluc.sbmlRole="species";
  Gluc.sbmlCompartment="compartment";
  Yi.sbmlRole="species";
  Yi.sbmlCompartment="compartment";
  V.sbmlRole="species";
  V.sbmlCompartment="compartment";
  Pi.sbmlRole="species";
  Pi.sbmlCompartment="compartment";
  Photoreceptor_activation.sbmlRole="rate";
  IfrSfrPfr.sbmlRole="parameter";
  Photoreceptor_inactivation.sbmlRole="rate";
  IrSrPr.sbmlRole="parameter";
  Transducer_activation.sbmlRole="rate";
  kia.sbmlRole="parameter";
  Transducer_inactivation.sbmlRole="rate";
  kai.sbmlRole="parameter";
  preS_formation.sbmlRole="rate";
  kx.sbmlRole="parameter";
  S_generation.sbmlRole="rate";
  ky.sbmlRole="parameter";
  Glucose_sensor_inactivation.sbmlRole="rate";
  kG.sbmlRole="parameter";
  S_formation.sbmlRole="rate";
  alpha1.sbmlRole="parameter";
  V_formation.sbmlRole="rate";
  alpha2.sbmlRole="parameter";
  S_degradation.sbmlRole="rate";
  kd_s.sbmlRole="parameter";
  V_degradation.sbmlRole="rate";
  kd_v.sbmlRole="parameter";
  Photoreceptor_decay.sbmlRole="rate";
  kd.sbmlRole="parameter";
}