Difference between revisions of "Team:USTC/Model/Performance Analysis/m"

Target of evaluation

---

In the part of system analysis, we have talked about the influence of promoter strengths and plasmid copy numbers to initial steady states and expression level of gene elements after sensing nicotine. The goal of our project is to detect nicotine and transform nicotine to more valuable chemicals, so it is necessary to evaluate the performance of sensing system like threshold and timeliness. Besides, because part of nicotine is degraded to 6-Hydroxynicotine which is out of the pathway of 3-Succinoyl-Pyridine metabolism, we need to calculate the rate of nicotine transformed to the target product. Parameter sweep provides us appropriate parameter intervals for system design, so effectiveness analysis is based on these parameter settings (s2 = 0.5, s1 = 0.5, plamid1 copy number =10, plasmid2 copy number = 10). We mainly focus on the threshold of sensor, time point suitable for fluorescence measurement, the recycle rate and the relationship between nicotine input, promoter1 strength and these three properties.

Threshold of sensing system

---

The usage of sensing system is to detect the existence and concentration of nicotine in the medium, so it has its own detection limit. Besides, in consideration of the complexity of the system, its sensitivity to nicotine input may vary a lot with different concentration of nicotine input. The combination of parameters is another important factor in the discussion of threshold and sensitivity because both low nicotine input and low expression of Vppa leads to inadequate 6-Hydroxynicotine to release the hdnoR repression of GFP signal.

The nicotine input we set is within the interval of [1e-5,1] M and the strength of promoter in is within [0,1]. To get the full vision of GFP change along the time span, the simulation runs for 200000 s. From figure 1 and figure 2, we find that the maximal GFP concentration ranges from 0 to about 2000 nM with the increase of nicotine input. However, the strength of promoter 1 has no significant impact on the GFP maximal concentration. The heat-map shows that when nicotine input is lower than 3.5e-5 M, the GFP concentration is relative low and undetectable.

In experiment, we use microplate reader to detect the fluorescence intensity of GFP. So the specificity and threshold are restricted by the detection limit and specificity of the device in fact. To fit the real condition in experiment, we need to convert GFP concentration (nM) to fluorescence intensity (au). According to the work of 2016 William and Mary iGEM team, we can build up the function between GFP concentration and fluorescence intensity like this:

We define fluorescence intensity more than 100 a.u. as significant result because the data of our interlab shows that the basal fluorescence of medium is about 60 a.u. under the same excitation wavelength light. This requirement determines a detection limit of GFP concentration. The lowest of nicotine concentration to be detected is decided in turn:

To get the value of fluorescence intensity of GFP, we transform the simulation data from GFP concentration (nM) to relative intensity (a.u.). Figure 3 shows the change tendency of GFP fluorescent intensity which is similar to the condition of GFP. To identify the boundary of relative intensity 100 a.u., we draw the contour line of 100 a.u. intensity. From the contour line (Figure 43), we find the when nicotine input is less than 1e-3.6 M, the relative intensity is lower than 100 a.u. and the result is not significant enough according to the criteria mentioned above.

In conclusion, the threshold of our sensing system is above 1e-3.6 M nicotine and the threshold does not strongly shift with the decline of promoter 1 strength.

Time point of signal processing

---

In the experiment, we need to choose the appropriate time to measure the GFP signal to determine the existence and concentration of nicotine in the medium. Besides, to collect the maximal quantity of final product, it is important to select the time point when 3-Succinoyl-Pyridine reaches highest concentration. This session we simulate the change of GFP/NicA2/3-Succinoyl-Pyridine and identify the start time of observable GFP signal and the duration time of GFP signal and NicA2 high expression. At last, according to simulation results, we propose the end time of degradation reaction as the reference to experiment design.

GFP is the fluorescent signal which means the existence of an amount of nicotine in the medium, so we first analyze the process of GFP signal production and attenuation. Figure 7 shows that the variation of time for reaching 100 a.u. is small among different groups and the average time is shorter than 5 min. However, the time for signal complete attenuation (under detection limit) varies a lot which increases greatly with more nicotine input. About the time for reaching maximal GFP intensity, it is interesting that more nicotine input and stronger promoter 1 may lead to longer period for GFP signal maximization, because more Vppa expression means more 6-Hydroxynicotine produced initially to release hdnoR repression quickly. Then we calculate the time interval, which is set to be the time for GFP intensity above 95% GFP maximum, for measuring GFP max intensity. The time interval changes a lot among various nicotine input and the maximal timespan is about 700 min.

The second signal is AHL molecule synthesized by LuxI. The lowest concentration of AHL to activate the transcription of Lux pR is about 0.1 nM, so we calculate the time in which AHL stays in a functional concentration. For nicotine input is lower than 0.1 M, the functional timespan is about 500-1000 min. With higher nicotine input, the state in functional AHL amount can last for over 2500 min.

In the last part, to collect the final product 3-Succinoyl-Pyridine in highest concentration, the time for reaching maximal concentration is calculated. From figure 12, we notice that the range of time is from 0-30 h and the change of time is strongly related to nicotine input but the impact of different promoter 1 strength seems not significant. The possible reason is that the time of GFP signal activation does not vary a lot (0-60min) as described above which means the time for signal induction does not account for a large portion in the total time for nicotine degradation.

The recycle rate of nicotine

---

The final goal of our system is to transform nicotine from tobacco waste to 3-Succinoyl-Pyridine which is a valuable intermediate for chemical synthesis of medicine. Because some portion of nicotine is catalyzed to 6-Hydroxynicotine to operate the sensing system, not all nicotine is transformed to 3-Succinoyl-Pyridine. We need to evaluate the extent of end product production with certain input of nicotine to validate the efficiency of degradation system. To describe the transformation rate, we define a concept called ‘Recycle rate’ calculated as below:

In simulation, we define the nicotine input in the interval of [1e-5, 1] M and promoter1 strength s1 in the interval of [0, 1]. To make sure that 3-Succinoyl-Pyridine reaches the maximal concentration, the simulation time is set to 300000s.

From figure 13, it is transparent that the relation between 1og10 nicotine input and product concentration follows the exponential function. But the amount of final product is not so strong related to promoter1 strength. Only when nicotine input is in high concentration, product amount tends to decline with the increase of s1. It is expectable because stronger promoter1 expresses more Vppa which transformed more nicotine to 6-Hydroxynicotine rather than 3-Succinoyl-Pyridine. About the recycle rate, we notice that when nicotine input is more than 0.01 M, stronger promoter1 leads to the attenuation of recycle rate. The maximal recycle rate is nearly 100% with small s1 which means that our system is efficient in transforming nicotine.