Team:HUST-China/Model/model-pro

Modeling

Prokaryote

This pathway consists of two signal inputs, three possible steady states: a strong signal at Pulse1 will cause a Gene of interest 1 expression, a strong signal at Pulse 2 will result in a Gene of interest 2 expression, and when both have strong inputs, both do NOT express. After designing this pathway, we want to figure out some properties, such as switching-response time, the stability of the entire system, etc. And these will be able to lead us to improve our circuit.

Modeling

In order to simulate the reaction of this prokaryotic pathway in engineering bacteria, we establish a delay differential equation system based on Michaelis-Menten equation and chemical reaction kinetics to simulate the operation of this pathway, and then translate it into mathematical language and into program.

  • Brief Parameter Table of DDEs Model

    more details
  • ParameterDescription ValueSource
    trc pulse1rate of pulse1 transcription -estimated
    trc pulse2rate of pulse2 transcription-estimated
    copynumcopy number of plasmid50estimated
    trl CⅠrate of CⅠ translation 0.05estimated
    trl CⅡrate of CⅡ translation 0.2estimated
    trl CⅢrate of CⅢ translation 0.3estimated
    trl mCrorate of Cro translation 0.6estimated
    trl Gene_of_interset_1rate of Gene_of_interset_1 translation 1estimated
    trl Gene_of_interset_2rate of Gene_of_interset_2 translation1estimated
    deg mCⅠrate of mCⅠ degradation 0.5estimated
    deg mCⅡrate of mCⅡ degradation 3.7estimated
    deg mCⅢrate of mCⅢ degradation 2.5estimated
    deg mCrorate of mCro degradation 3.65estimated
    deg mGene_of_interset_1rate of mGene_of_interset_1 degradation 5.7estimated
    deg mGene_of_interset_2rate of mGene_of_interset_2 degradation0.5estimated
    deg CⅠrate of CⅠ degradation 0.03estimated
    deg CⅡrate of CⅡ degradation 0.74estimated
    deg CⅢrate of CⅢ degradation 0.5estimated
    deg Crorate of Cro degradation 0.5estimated
    deg Gene_of_interset_1rate of Gene_of_interset_1 degradation 1estimated
    deg Gene_of_interset_1rate of Gene_of_interset_2 degradation1estimated
    Vmax mCⅡmaximum transcription rate of promoter pRE0.2estimated
    Vmax mCⅢmaximum transcription rate of promoter pR0.15estimated
    Vmax mCromaximum transcription rate of promoter pRE0.2estimated
    Vmax mgene_of_ineterest_1maximum transcription rate of promoter pR5estimated
    Vmax mgene_of_ineterest_2maximum transcription rate of promoter pRM5estimated
    KmCⅡapparent association constant for CⅠ binding with pRE0.3estimated
    KmCⅠ2apparent association constant for CⅠ2 binding with pR0.5estimated
    KmCroapparent association constant for Cro binding with pRM0.025estimated
    K1CI dimerization reaction constant0.05Burz et al.(1994); Shea and Ackers(1985)
    K2CI dimerization reaction constant0.6estimated
    VcⅡmaxmaxinum rate of Ftsh break down CⅡ2.4estimated
    VcⅢmaxmaxinum rate of Ftsh break down CⅢ2.5estimated
    KCⅡMichaelis-Menten constant of FtsH breaks down CⅡ1.5estimated
    KmCⅢMichaelis-Menten constant of FtsH breaks down CⅢ1estimated
    τ1time for CⅠ2 binding with promoter pR1estimated
    τ2time for CⅡ activating promoter pRE1estimated
    τ3time for CⅡ activating promoter pRE1estimated
    τ4time for Cro binding with promoter pRM1estimated
    τ5time for CⅠ2 binding with promoter pR(Gene)1estimated

  • Formulary

    more details

Analysis

Then we can solve some specific problems with the help of these equations and mathematical tools, and get some basic properties about the pathway and to evaluate its working state:

  • Test of switch function

    After designing a switch-enabled path, the primary task is to determine whether it can work properly, and secondly, to understand how the signal input conditions will lead to switching.

    We adjust the pulse1 and pulse2 input signal strength ratio and statistics of product expression at this time to observe the work of the state of this pathway.

    We can see that our switching circuit faithfully fulfills his duty to express only one protein under a specific signal input. While the other one is suppressed to be at a low level of expression. Besides, from the figure we can clearly figure out when the ratio between the two signals is about 10, we can realize the conversion.

  • Test of filtering performance

    We did the following test to check out our filtering performance:

    When we keep ratio between the G1 / G2 signal about 1:10, we input a sinusoidal signal with Gaussian white noise to G2:

    The ratio of this noise signal is 1. And if our filter can handle it, it shows that our filter can work as a weakening or even eliminating part against most of the biological noise.

    We can see that the output does not contain too much noise in the waveform, which is still a typical protein generation curve. And we have the output signal analysis. We find that the output noise signal ratio increased to 75 (without filtering the ratio is 1), the majority of the noise has been removed.

    By tracking the variables in the delay differential equation, we can explain the principle of our filter from another point of view:

    We have the following reaction in our design:When there are still Cl and Cl2 in the cell. The balanced system made of these two components works as the capacitance in the circuit to adjust signal fluctuations. When Cl is increased due to signal fluctuations, Cl2 synthesis rate increases to consume more Cl, vice versa. It also explains why such a system of G2 filtering effect performances better than a simply attenuation as the G1. Since this balanced system can play a role as a capacitor in the general circuit, we can thus imitate the basic circuit in electronic circuits (such as LC shocks, etc.) to design the biological circuit to realize different functions.

  • Parameter Sensitivity Analysis / Robustness Evaluation

    Next, we analyzed the sensitivity of our pathway parameters and take a look at two important parameters in our pathway, the robustness of the response time and the expression level:

    Firstly, we obtained the parameters that are more sensitive to our model by pre-analysis, and then we adjusted them one by one and observed the change of the steady-state and the response time of the model.Another reason why we chose them was that we could use a variety of ways to “change” them, like providing another competitive binding substance or modify the promoter binding sites.

    1. KmCI2(apparent association constant for CI2 binding with pR)


    We can see that when the value of KmCI2 increases, the response time of G1 decreases but the expression level rises to a stable value, and G2 shows considerable stability, which gives us a revelation: when we improve the circuit, we can use the promoter of larger inhibition constant to reduce the response time while keeping the expression level unchanged.


    2. Kmcro(apparent association constant for Cro binding with pRM)


    Similarly, the response time of G1 and G2 showed fairly good robustness, but the expression of G2 was sensitive to the value of Kmcro . With the increase of Kmcro, the expression of G2 also increased. But we can also use its sensitivity to control expression levels.

    3.KmcII(apparent association constant for CII binding with pRE)

    From the figure above, we can see that the response time and expression level of G1 are relatively stable, while G2 is more sensitive to this parameter. However, considering that KmcII value is very small, G2 is almost non-expressed and the leakage of gene expression cannot be avoided, the fluctuation is due to G2 leakage expression. It can be considered like that G2 response time is stable. Another interesting point is that the sudden increase in G2 expression can be considered as a switch, which provides a way for us to design more multivariable switches later. We can turn G2 "off" by providing another competitive binding substance to "change" the value of Kmcll.In addition, we can also modify the promoter binding sites to achieve our goal.

Summary

In general, the designed filter and switch functions are achieved. And in the case of little external disturbance our system can be relatively stable. This prokaryotic pathway also exhibits a high level of plasticity, which provides us with many ways for further application and future improvement.