[2696] | 1 | {Model of hysteresis and irreversible transition in eukaryotic cell} |
---|
| 2 | |
---|
| 3 | #DEFVAR |
---|
| 4 | A = ignore; {protein that activates Cdh1 at finish} |
---|
| 5 | AYp= ignore; {complex} |
---|
| 6 | Yp = ignore; {inactive (phosphorylated) Cdh1-P/APC} |
---|
| 7 | XY = ignore; {complex} |
---|
| 8 | Y = ignore; {active Cdh1/APC complex} |
---|
| 9 | X = ignore; {cyclin/Cdk dimer} |
---|
| 10 | YX = ignore; {complex} |
---|
| 11 | XYX= ignore; {complex} |
---|
| 12 | G = ignore; {complex} |
---|
| 13 | TF = ignore; {complex} |
---|
| 14 | TFp= ignore; {complex} |
---|
| 15 | C1 = ignore; {complex} |
---|
| 16 | C2 = ignore; {complex} |
---|
| 17 | mas= ignore; {mass of the cell} |
---|
| 18 | |
---|
| 19 | |
---|
| 20 | #EQUATIONS |
---|
| 21 | mas = 2mas : (0.01); |
---|
| 22 | 2mas = mas : (0.001); |
---|
| 23 | PROD = X : (0.02); {k1} |
---|
| 24 | X = PROD : (0.05); {k2 prime} |
---|
| 25 | Y + X = YX : (1000); {k2f double prime} |
---|
| 26 | YX = Y + X : (3); {k2r double prime} |
---|
| 27 | YX = Y : (1); {k2 double prime} |
---|
| 28 | X + Y = XY : (1200); {k4f} |
---|
| 29 | XY = X + Y : (35); {k4r} |
---|
| 30 | XY = X + Yp : (80); {k4} |
---|
| 31 | A + Yp = AYp : (4000); {k3f} |
---|
| 32 | AYp = A + Yp : (75); {k3r} |
---|
| 33 | AYp = A + Y : (60); {k3} |
---|
| 34 | X + YX = XYX : (360); {k5f} |
---|
| 35 | XYX = X + YX : (18); {k5r} |
---|
| 36 | XYX = 2X + Yp: (16); {k6} |
---|
| 37 | A = PROD : (0.4); {k6 prime} |
---|
| 38 | AYp = Yp : (0.4); {k6 prime} |
---|
| 39 | G = G + A : (0.00075);{k5 prime} |
---|
| 40 | TF + X = X + TFp: (2.6); {k7f} |
---|
| 41 | TFp = TF : (4.0); {k7r} |
---|
| 42 | TFp + G = C1 : (1.2); {k8f} |
---|
| 43 | C1 = TFp + G : (320); {k8r} |
---|
| 44 | C1 = C1 + A : (0.15);{k9} |
---|
| 45 | TFp + C1= C2 : (65); {k10f} |
---|
| 46 | C2 = C1 + TFp : (0.2); {k10r} |
---|
| 47 | C2 = C2 + A : (3); {k11} |
---|
| 48 | |
---|
| 49 | #LOOKAT X; |
---|
| 50 | #MONITOR X; |
---|
| 51 | |
---|
| 52 | #INITVALUES |
---|
| 53 | CFACTOR = 1; |
---|
| 54 | mas=0.1046E+01; |
---|
| 55 | TF=0.8979E+01; |
---|
| 56 | XY=0.5458E-02; |
---|
| 57 | XYX=0.8918E-02; |
---|
| 58 | YX=0.2989E-02; |
---|
| 59 | X=0.1275E+00; |
---|
| 60 | Y=0.1714E-02; |
---|
| 61 | Yp=0.9605E+00; |
---|
| 62 | AYp=0.2046E-01; |
---|
| 63 | A=0.8577E-03; |
---|
| 64 | C2=0.1692E-02; |
---|
| 65 | C1=0.1769E-04; |
---|
| 66 | TFp=0.1731E-01; |
---|
| 67 | G=0.9983E+00; |
---|
| 68 | |
---|
| 69 | #INLINE F77_INIT |
---|
| 70 | TSTART = 0.d0 |
---|
| 71 | TEND = 1000.d0 |
---|
| 72 | DT = 5.d0 |
---|
| 73 | Volume = 1000.d0 |
---|
| 74 | #ENDINLINE |
---|
| 75 | |
---|
| 76 | #INLINE F77_UTIL |
---|
| 77 | DOUBLE PRECISION FUNCTION CellMass(t) |
---|
| 78 | ! Cell mass evolution in time |
---|
| 79 | DOUBLE PRECISION t |
---|
| 80 | DOUBLE PRECISION m0, mstar, mu, one |
---|
| 81 | PARAMETER (m0 = 1.0d0) |
---|
| 82 | PARAMETER (mstar = 10.0d0) |
---|
| 83 | PARAMETER (mu = 0.01d0) |
---|
| 84 | PARAMETER (one = 1.0d0) |
---|
| 85 | CellMass = mstar/(one-(one-mstar/m0)*EXP(-mu*t)) |
---|
| 86 | END |
---|
| 87 | #ENDINLINE |
---|
| 88 | |
---|
| 89 | #INLINE F90_INIT |
---|
| 90 | TSTART = 0.d0 |
---|
| 91 | TEND = 1000.d0 |
---|
| 92 | DT = 5.d0 |
---|
| 93 | !Volume = 1000.d0 |
---|
| 94 | #ENDINLINE |
---|
| 95 | |
---|
| 96 | #INLINE F90_UTIL |
---|
| 97 | DOUBLE PRECISION FUNCTION CellMass(t) |
---|
| 98 | ! Cell mass evolution in time |
---|
| 99 | DOUBLE PRECISION, INTENT(IN) :: t |
---|
| 100 | DOUBLE PRECISION, PARAMETER:: m0 = 1.0d0 |
---|
| 101 | DOUBLE PRECISION, PARAMETER:: mstar = 10.0d0 |
---|
| 102 | DOUBLE PRECISION, PARAMETER:: mu = 0.01d0 |
---|
| 103 | DOUBLE PRECISION, PARAMETER:: one = 1.0d0 |
---|
| 104 | CellMass = mstar/(one-(one-mstar/m0)*EXP(-mu*t)) |
---|
| 105 | END FUNCTION CellMass |
---|
| 106 | #ENDINLINE |
---|
| 107 | |
---|
| 108 | #INLINE MATLAB_INIT |
---|
| 109 | global TSTART TEND DT TEMP |
---|
| 110 | TSTART = 0; |
---|
| 111 | TEND = 1000; |
---|
| 112 | DT = 5.0; |
---|
| 113 | Volume = 1000.0; |
---|
| 114 | #ENDINLINE |
---|
| 115 | |
---|
| 116 | #INLINE MATLAB_UTIL |
---|
| 117 | function m = CellMass(t) |
---|
| 118 | % Cell mass evolution in time |
---|
| 119 | m = 10/(1+9*exp(-0.01*t)); |
---|
| 120 | return |
---|
| 121 | #ENDINLINE |
---|
| 122 | |
---|
| 123 | #INLINE C_INIT |
---|
| 124 | TSTART = 0.0; |
---|
| 125 | TEND = 1000.0; |
---|
| 126 | DT = 5.0; |
---|
| 127 | Volume = 1000.0; |
---|
| 128 | #ENDINLINE |
---|
| 129 | |
---|
| 130 | #INLINE C_UTIL |
---|
| 131 | double CellMass(double t) |
---|
| 132 | {/* Cell mass evolution in time */ |
---|
| 133 | const double m0 = 1.0; |
---|
| 134 | const double mstar = 10.0; |
---|
| 135 | const double mu = 0.01; |
---|
| 136 | const double one = 1.0; |
---|
| 137 | return mstar/(one-(one-mstar/m0)*exp(-mu*t)); |
---|
| 138 | } |
---|
| 139 | #ENDINLINE |
---|