hong1
The SBML for this model was obtained from the BioModels database (BioModels ID: BIOMD0000000216). Biomodels notes: The figure 2C of the reference publication is reproduced. To reproduce figure 2C, the following initial conditions were used (as per the authors suggestion). This is different from that of the model. M=1.35; CP=0.039; CP2=0.088; TF=0.06; IC=0.44 (TFtot=0.5 & IC=TFtot-TF). JWS Online curation: This model was curated by reproducing the figures as described in the BioModels Notes.
None
None
None
None
None
None
Unit definitions have no effect on the numerical analysis of the model. It remains the responsibility of the modeler to ensure the internal numerical consistency of the model. If units are provided, however, the consistency of the model units will be checked.
| Name | Definition |
|---|
| Id | Name | Spatial dimensions | Size | |
|---|---|---|---|---|
| system | system | 3.0 | 1.0 |
| Id | Name | Initial quantity | Compartment | Fixed | |
|---|---|---|---|---|---|
| CP | Monomer clock proteins (PER) | — | system (system) | ✘ | |
| CP2 | Dimer form of clock proteins (PER/PER) | — | system (system) | ✘ | |
| M | Messenger RNA of the clock proteins (Per mRNA) | — | system (system) | ✘ | |
| TF | Transcription factor (BMAL1/CLK) | — | system (system) | ✘ |
Initial assignments are expressions that are evaluated at time=0. It is not recommended to create initial assignments for all model entities. Restrict the use of initial assignments to cases where a value is expressed in terms of values or sizes of other model entities. Note that it is not permitted to have both an initial assignment and an assignment rule for a single model entity.
| Definition | |
|---|---|
| CP = initialValue_CP | |
| CP2 = initialValue_CP2 | |
| M = initialValue_M | |
| TF = initialValue_TF |
| Id | Name | Objective coefficient | Reaction Equation and Kinetic Law | Flux bounds | |
|---|---|---|---|---|---|
| r1 | r1 | ∅ > M system * function_4_r1 | |||
| r10 | r10 | ∅ > TF + CP2 system * kicd * IC | |||
| r11 | r11 | CP2 > ∅ system * kcp2d * CP2 | |||
| r12 | r12 | TF + CP2 > ∅ system * kica * CP2 * TF | |||
| r13 | r13 | CP2 > ∅ system * function_4_r13 | |||
| r14 | r14 | CP2 > ∅ system * chk2 * CP2 | |||
| r15 | r15 | ∅ > TF system * kcp2d * IC | |||
| r16 | r16 | ∅ > TF system * chk2c * IC | |||
| r17 | r17 | ∅ > TF system * function_4_r17 | |||
| r2 | r2 | ∅ > M system * function_4_r2 | |||
| r3 | r3 | M > ∅ system * kmd * M | |||
| r4 | r4 | ∅ > CP system * function_1 | |||
| r5 | r5 | CP > ∅ system * kcpd * CP | |||
| r6 | r6 | {2.0}CP > CP2 system * ka * pow(CP, 2) | |||
| r7 | r7 | CP2 > {2.0}CP system * kd * CP2 | |||
| r8 | r8 | CP > ∅ system * function_4_r8 | |||
| r9 | r9 | CP > ∅ system * chk2 * CP |
| Id | Value | |
|---|---|---|
| CPtot | 0.0 | |
| Dex | 0.0 | |
| IC | 0.0 | |
| ICtot | 1.0 | |
| IR_pulse_duration | 2.0 | |
| J | 0.3 | |
| Jp | 0.05 | |
| TFtot | 0.5 | |
| chk2 | 0.0 | |
| chk2ToggleOnOff | 1.0 | |
| chk2_value | 0.25 | |
| chk2c | 0.0 | |
| chk2c_value | 0.1 | |
| chk2pulseactive | 1.0 | |
| chk2pulseactive_elim | 1.0 | |
| chk2pulseactive_elim_elim | 0.0 | |
| chk2pulseonoff | 0.0 | |
| function_1 | 0.0 | |
| function_4_r1 | 0.0 | |
| function_4_r13 | 0.0 | |
| function_4_r17 | 0.0 | |
| function_4_r2 | 0.0 | |
| function_4_r8 | 0.0 | |
| initialValue_CP | 0.039 | |
| initialValue_CP2 | 0.088 | |
| initialValue_M | 1.35 | |
| initialValue_TF | 0.06 | |
| ka | 100.0 | |
| kcp2d | 0.0525 | |
| kcpd | 0.525 | |
| kcps | 0.5 | |
| kd | 0.01 | |
| kica | 20.0 | |
| kicd | 0.01 | |
| kmd | 0.1 | |
| kms | 1.0 | |
| kp1 | 10.0 | |
| kp2 | 0.1 | |
| n | 2.0 | |
| trigOFFtime | 0.0 | |
| trigONtime | 4.0 |
| Id | Value | Reaction |
|---|
| Definition | |
|---|---|
| IC = TFtot - TF | |
| function_1 = kcps * M | |
| function_4_r17 = kp2 * IC / (Jp + CP + 2.0 * CP2 + 2.0 * IC) / system | |
| function_4_r2 = kms * pow(TF, n) / (pow(J, n) + pow(TF, n)) / system | |
| chk2pulseonoff = chk2pulseactive - (chk2pulseactive_elim_elim + chk2pulseactive_elim) | |
| function_4_r8 = kp1 * CP / (Jp + CP + 2.0 * CP2 + 2.0 * IC) / system | |
| function_4_r13 = kp2 * CP2 / (Jp + CP + 2.0 * CP2 + 2.0 * IC) / system | |
| function_4_r1 = Dex / system | |
| CPtot = CP + 2.0 * CP2 + 2.0 * IC | |
| trigOFFtime = trigONtime + IR_pulse_duration | |
| chk2c = chk2ToggleOnOff*(chk2pulseonoff * (chk2c_value / IR_pulse_duration)) | |
| chk2 = chk2ToggleOnOff*(chk2pulseonoff * (chk2_value / IR_pulse_duration)) |
| Definition |
|---|
| Definition |
|---|
| Definition |
|---|
| Trigger | Assignments | |
|---|---|---|
| geq(time, trigONtime) | chk2pulseactive_elim = 0 | |
| geq(time, trigOFFtime) | chk2pulseactive_elim_elim = 1 |