The following reaction scheme models a situation in which a modified gene produces two mRNAs which stochastically transcribe different reporter proteins, such as GFP and CFP, independently of one another forming the basis of a dual reporter method. The system features feedback from the mRNA given by its production rate being a function
The Gillespie algorithm was implemented in R and run to simulate two situations with the following functions:
$f(x_{1}) = \lambda_{1}$ - Whereby mRNA is transcribed at a constant rate
$f(x_{1}) = \lambda_{1} \frac{K}{K + x_{1}}$ - Whereby mRNA self-repressive as its transcription is inversly proportional to its own concentration
$(x_{1})$ and a rate constant$K$
For both functions, rate constants were set with