Endy:Computational modeling of demand/Species and reactions

Transcription species

 * RNAP [P]
 * Promoter [D]
 * RNAP-Promoterclosed [PDc]
 * RNAP-Promoteropen [PDo]
 * RNAP-Promoterinit [PDi]
 * RNAPelong [Pe]
 * mRNA [M]

It should be possible to ignore some of these species as they don't affect the dynamics of the transcription model. TABASCO appears to skip the closed and open RNAP-Promoter species. Lumping the intermediate species together in this way makes some sense as I am unlikely to have access to rate constants for their formation and dissociation. The other advantage to this lumping is that transcription initiation is then analogous to a Michaelis-Menten enzyme reaction.

For a first iteration, I will consider the following species -


 * [P] - stable
 * [D] - stable
 * [PDi] abbreviated to [PD] - stable
 * [Pelong] - stable
 * [M] - unstable

Translation species

 * Ribosome [R]
 * RBS [M]
 * RBS-Ribosomeinit [RM]
 * Ribosomeelong [Relong]
 * Protein [O]

Reactions
Polymerases (P), DNA (D) and Ribosomes (R) are all synthesized at a zeroth order rate. The synthesis rate was set by choosing a desired total number of the particular species per cell. The synthesis was set such that in the abscence of reactions other than dilution due to the growth, the species would reach that desired steady state level. The synthesis rate for species i was - rsyn,i = Css,i*&mu;

Where Css,i is the desired steady state number of species i per cell and &mu; is the rate constant for dilution due to growth.