Reaction scheme where the products are created from the reactants and the change of a product quantity is proportional to the product of reactant activities. The reaction scheme does not include any reverse process that creates the reactants from the products. The change of a product quantity is proportional to the quantity of one reactant.
Two-gene model of the Ostreococcus circadian clock
This is a model of the circadian clock of Ostreococcus tauri, with a negative feedback loop between TOC1 and CCA1 (a.k.a. LHY) and multiple light inputs. It was used and described in Troein et al., Plant Journal (2011).
This version of the model has eliminated one light input by setting the light accumulator acc to 1, corresponding to Figure 2F of Dixon et al., New Phytologist (2014). The model has been entrained in LD 12:12 and will be given one LD cycle before going into constant light. The length of the day before LL is defined in the function called light.
The model incorporates luciferase reporters, and in this SBML model the four different versions of the model for transcriptional and translational reporter lines (pTOC1::LUC, pCCA1::LUC, TOC1-LUC and CCA1-LUC) are all accessible by setting one of the rep_X parameters to 1 and the others to 0. You can also set all four to 0 to only simulate the non-reporter core of the system.
Input to the system should be provided by modifying the "light" function. Implementations of an LD cycles and transition to LL are provided as examples, but the model was also used with more complicated light regimes that vary between data sets and are not convenient to express directly in SBML.
The functions "ox_cca1" and "ox_toc1" can be altered to add overexpression of CCA1 and TOC1. Setting either to x gives additional, constitutive transcription at x times the maximal (and typically not realizable) transcription rate of the native gene. The overexpression mutant fits in Figure 7 of Troein et al. (2011) used ox_cca1 = 0.115 and oc_toc1 = 0.0584, respectively.
The functions "copies_toc1" and "copies_cca1" are normally 1 but can be lowered to simulate knockdown experiments. The functions "transcription", "translation" and "proteasome" can be modified to simulate the effects of altering the overall rate of transcription, translation and protein degradation.
The parameters were fitted specifically to data from transgenic reporter lines TOC8, pTOC3, LHY7 and pLHY7 (Corellou et al., Plant Cell 2009). Parameters that begin with "effcopies" describe the effective number of copies of CCA1 or TOC1 in the respective translational fusion lines, with anything above 1 due to the fusion proteins.
For the model fitting, the initial values were fitted to the data in the various time courses. The initial values given here correspond to the limit cycle of the system in LD 12:12. The system converges to the limit cycle in just a few days under most light conditions, so these initial values are biologically meaningful.
The species cca1luc_c and cca1luc_n have been merged into cca1luc (which corresponds to the observable luminescence signal), because Copasi refused to run the system otherwise. For TOC1-LUC, the predicted output signal is the sum of toc1luc_1 and toc1luc_2.