Our original model was based on the assumption that the T7 RNA polymerase is constitutively expressed and thus that IPTG addition instantaneously “turns on” the expression of the mRNA and FAP (reactions 3-12 in Figure 1). In this model, reactions 3 - 7 involve mRNA. The T7 RNAP polymerase is not consumed in producing mRNA and mRNA is not consumed in producing the FAP.

We assumed that the fluorescence signal is proportional to the concentration of the [mRNA-DFHBI]. The proportionality constant K includes factors specific to the fluorimeter settings and K was one of the variable parameters in our fitting algorithm.

The dosage curve indicated that the dye binding reactions (6, 7 and 10, 11) are not rate-limiting and allowed us to measure the equilibrium constant K_D for these reactions. K_D is equal to the ratio of the “off” rate (k_d) to “on” rate (k_a).

The degradation of [dye-mRNA] and [dye-FAP] complexes affects only the biological component of the complex and not the dye. We assumed that the two complexes degrade with the same rate constant, α and β, as their mRNA and FAP unbound counterparts, respectively.

The model described above could not lead to good simulations of the time dependence of the fluorescence of the [mRNA-DFHBI] complex. The apparent exponential dependence of the [mRNA-DFHBI] on time could be rationalized if one includes in the model two additional reaction that account for the fact that in the regulatory mechanism the T7 RNA polymerase in the BL21(DE3) cell strain, there is a LacO site upstream of the promoter for the polymerase. We accounted for this process by including reactions 1 and 2 in the Figure 1. Reaction 1 represents the induction of the T7 RNAP by adding a fixed amount of IPTG and reaction 2 represents the degradation of T7 RNAP.

The following system of differential equations represent the chemical reactions in figure 1.

K_D's were taken from the original papers Szent-Gyorgi et al and Paige et al.

The initial concentrations of dyes and IPTG are known. The [IPTG] value is constant. k₁ should be determined so that the steady state for T7RNAP is accurate. We implemented these equations in MATLAB using an ODE function and a stochastic simulation approach.

For low dye concentration, the fluorescence increases until all the dye is bound. One can take advantage of this property in a titration experiment to determine the maximum concentration of RNA or protein produced by the cell.

All dye bound condition:

Excess dye condition:

Limited Dye Added:

Reaction 1: IPTG->T7RNAP+IPTG
Reaction 2: T7RNAP->0
Reaction 3: T7RNAP->mRNA+T7RNAP
Reaction 4: mRNA->0
Reaction 5: mRNA->FAP+mRNA
Reaction 6: mRNA+DFHBI->[mRNA-DFHBI]
Reaction 7: [mRNA-DFHBI]->mRNA+DFHBI
Reaction 8: [mRNA-DFHBI]->0+DFHBI
Reaction 9: [mRNA-DFHBI]->FAP+[mRNA-DFHBI]
Reaction 10: FAP+MG->[FAP-MG]
Reaction 11: [FAP-MG]->FAP+MG
Reaction 12: FAP->0
Reaction 13: [FAP-MG]->0+MG