This pathway consists of two signal inputs, three possible steady states: a strong signal at Pulse1 will cause a Gene of interest1 expression, a strong signal at Pulse2 will result in a Gene of interest 2 expression, and when both have strong inputs, both Not expressed.After we have designed this pathway, naturally, 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 the program.
-
Brief Parameter Table of DDEs Model
more details -
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 he 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 is suppressed 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 quickly.
-
Test of filtering performance
We do 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 CI and CI2 in the cell. The balanced system made of these two components works as the capacitance in the circuit to adjust signal fluctuations. When CI is increased due to signal fluctuations, CI2 synthesis rate increases to consume more CI, 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 will analyze 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 adjust them one by one and observe the change of the steady-state and the response time of the model
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)
Similar to Kmcro, the response times 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 this to control expression levels.
3.KmcII(apparent association constant for CII binding with pRE)
From the figure, 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. By providing another competitive binding substance to "change" KmcII to implement control G2 by the switch.
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.
