Modeling

Potassium sensor

Phosphate sensor

Model's Assumption

The effect of the endogenous Kdp system of E.coli was neglected.

In our engineered E.coli, it contains the inserted plasmid of PKdpF plus green fluorescence protein generator (BBa_E0240) in pSB3K3 backbone as well as the endogenous Kdp operon, and both operons have the same promoter, PKdpF. As a matter of fact, titration by the endogenous kdp operon of the transcription inducer, phosphorylated KdpE which binds to PKdpF was expected initially. And the titration of phosphorylated KdpE is anticipated to lower the expression of GFP. However, when the DNA copy number of both endogenous and the inserted operons were determined, it was found that the number of endogenous kdp operon is 10 times smaller than that of the inserted one:

This implied that the number of endogenous promoter PKdpF is ten times smaller than that of the inserted operon and accounts only 8.33% of the total number of promoter PKdpF in the engineered E.coli. Therefore, the titration effect of phosphorylated KdpE becomes insignificant. As a result, the effect of endogenous Kdp system was negligible.

Level of KdpD, KdpE and KdpF were assumed to be constant

In accordance to [Kremling A. 04], in the potassium ions concentration range which we were studying, from 0mM to 0.02 mM, the fluctuation of the concentration KdpF as well as KdpD and KdpE was only within 10 µM and 3 µM respectively. Due to the small fluctuation range compared to the gene expression of GFP reporter, it was reasonable to assume the concentration of KdpD, KdpE and KdpF to be constant in the model.

It was assumed that the initial concentration of mRNA for GFP, immature GFP and mature GFP equal to zero.

It was assumed that all reactions below were in steady state such that:

Equations of the Model:

Phosphorylation of KdpD:

Phosphyl-group Transfer:

Binding of KdpE to promoter PkdpF:

Transcription:

Translation:

Green Fluorescent Protein maturation:

Data point fitting to the model:

The experimental data from FACS were used to fit in the prediction model; then we adapt the unit conversion from [Caitlin C, Jeniffer B 13] to convert GFP per cell fluorescence intensity to concentration of GFP per cell:

Prediction model:

Parameters and Variables list

Parameters:	Values:
k₁: Rate constant of forward auto-phosphorylation	0.23 1/µM hr
k_-1:Rate constant of backward auto-phosphorylation	0.0000051 1/µM hr
k₂: Rate constant of forward phosphyl-group transfer	0.00227 1/µM hr
k_-2: Rate constant of backward phosphyl-group transfer	0.00087 1/µM hr
K_b: Dissociation constant (KdpE-DNA)	0.0532 1/µM
K_h: Dissociation constant (KdpF-KdpE dephosphorylation system)	0.0532 µM (adjustable)
k⁰_h: Reference k_h value	0.00023 1/µM hr (adjustable)
C⁰_K+: Reference concentration of potassium ions	0.02mM
k_tr: transcription rate constant	10600 hr^-1
k_z: mRNA degradation rate of kdpFABC and kdpDE coding region	21.74 hr^-1
k_d: degradation rate of KdpDE	0.2 hr^-1
k_d,F: degradation of KdpF	4.8 hr^-1
k_tlD: translation rate constant of KdpD	160 hr^-1
k_tlE: translation rate constant of KdpE	8100 hr^-1
k_tlF: translation rate constant of KdpF	5.4 hr^-1
u: growth rate of the cell (in K minimal medium, 2^nd order)	0.275466667 hr^-1 (exp. obtained)
Ψ: basal activity of the promoter	0.0017
k_deg,GFP,mRNA: GFP mRNA degradation rate (1^st order)	0.0048s^-1 =17.28 /hr
γ_I: Immature GFP degradation rate	0.7 hr^-1
γ_m: Mature GFP degradation rate	0.021 hr^-1
ρ: Immature GFP synthesis rate (trans rate, 1^st order)	1440 protein/hr/mRNA (0.4protein/s/mRNA)
w: GFP maturation rate (1^st order)	6.48 hr^-1

Table 1. Parameters list 1.

Parameters	Values
ADP concentration	0.305 mM
ATP concentration	3.05 mM
Initial concentration of mRNA (GFP)	0/cell
Concentration of KdpF	10 µM
Concentration of KdpE	2.6 µM
Concentration of KdpD	2.6 µM
Concentration of inserted plasmid	0.0183 µM
Initial concentration of immature and mature GFP	0 /cell

Table 2. Parameters list 2.

Variables	Denotation	Unit
	Concentration of KdpD	µM
	Concentration of phosphorylated KdpD	µM
	Concentration of phosphorylated KdpE	µM
	Concentration of ADP molecules	mM
	Concentration of ATP molecules	mM
	Concentration of KdpF	µM
	Rate constant of dephosphorylation of KdpE

	Concentration of potassium ions	mM
	Concentration of KdpE bound promoter P_kdpF	µM
	Concentration of Unbound promoter P_kdpF	µM
	Total concentration of promoter P_kdpF	µM
	Concentration of mRNA molecules for green fluorescence protein
	Concentration of immature green fluorescent protein	µM
	Concentration of mature green fluorescent protein	µM

Table 3. Variables list.

Reference:

Heermann R, Zigann K, Gayer S, Rodriguez-Fernandez M, Banga JR, et al. (2014) Dynamics of an Interactive Network Composed of a Bacterial Two- Component System, a Transporter and K+ as Mediator. PLoS ONE 9(2): e89671. doi:10.1371/journal.pone.0089671

Brewster RC, Jones DL, Phillips R (2012) Tuning Promoter Strength through RNA Polymerase Binding Site Design in Escherichia coli. PLoS Comput Biol 8(12): e1002811. doi:10.1371/journal.pcbi.1002811

Modeling, Simulation and Identification of the Dynamics of K Uptake in E. coli. (2014). Universitätsbibliothek der TU München.

Kelly, Jason et al. “Measuring the activity of BioBrick promoters using an in vivo reference standard.” Journal of Biological Engineering 3.1 (2009): 4.

J. Gayer, Stefan. "Modeling, Simulation and Identification of the Dynamics of K Uptake in E. Coli." Technische Universität München Fachgebiet Für Systembiotechnologie (2013). Print.

Kremling, A., Heermann, R., Centler, F., & Gilles, E. (2004). Analysis of two-component signal transduction by mathematical modeling using the KdpD/KdpE system of Escherichia coli.

Conboy, C., & Braff, J. (2013, May 29). Molecules of Equivalent GFP. Retrieved from http://openwetware.org/wiki/MEG