Team:HIT-Harbin/Modeling

From 2013.igem.org

Revision as of 06:03, 14 August 2013 by Sema1688 (Talk | contribs)
Home Team Official Team Profile Project Parts Submitted to the Registry Modeling Notebook Safety Attributions

Before experiment, a mathematical model is developed to describe our system and make a prediction about our following experiment. Our model is based on the following assumption.

(1) Transcription follows Hill equation, other reactions follows mass action principle and saturation kinetics.

(2) Parameter values are partly found on current publications[1,2] and partly guess in a reasonable way.

(3) The amount of protein binding to promoter is negligible because the amount of protein binding to promoter is trace compared to the total amount of the protein

Model of AND gate

Figure 1 : A Schematic of the Biological Andgate

As described before, HrpS binds to HrpR forming a complex HrpRS. HrpRS then binds to promoter hrpL triggering transcription of the downstream gene. We use two inducible promoters--the isopropylthiogalactoside (IPTG)-inducible Plac and the arabinose-inducible PBAD as the AND gate inputs. A quantitative dynamic model based on protein and protein-promoter interactions for the hrp gene regulatory is described below:

1 Transcription of hrpR and hrpS mRNA

A transcription function named Hill function[3] was used to explain how IPTG, Arabinose and HrpRS influence the promoters. The Hill function and its curve are

Figure 2 : Example of regulation functions (a)Hill function;(b)Step function;(c)logoid function

θj is the threshold for the regulatory influence of the xj on a target gene, m is the steepness parameter. h+ explain the positive regulation on the target gene, while h- is the negative regulation function.

IPTG binds to promoter Plac and trigger the transcription of hrpR_mRNA at a maximum rate constant vmr,

while hrpR_mRNA degrade at a rate constant Drnar, so the dynamic change for hrpR_mRNA is

Similarly, we get the reaction for PBAD, they are

2 Translation of HrpR and HrpS protein

To simplify our model, we made some assumption

(1) the level of translation product is assumed to be proportional to mRNA levels

(2) There is no delay in synthesis of either component or delay because of protein transportation

HrpR and HrpS are synthesized at the rate constants Kr and Ks respectively, and degrade at the rate constant Dr and Ds.

The net forming rate of HrpR and HrpS are

3 Interaction between protein HrpRS and promoter hrpL

HrpR and HrpS binding rate constant is define as Krs, and their dissociation rate constant is K_rs. Complex HrpRS degrade at a rate constant Drs

The net forming rate of HrpRS is

Activating for the promoter hrpL is like Plac and PBAD

Promoter hrpL trigger the transcription of our target protein. Ko and Do denote the synthetic and degradation rate constant for output protein.

The whole model for the hrp AND gate are summarize below