Difference between revisions of "Team:XMU-China/Proof"

Line 61: Line 61:
 
}
 
}
 
@media (min-width:769px){
 
@media (min-width:769px){
−
       .jumbotron h1{width:52%;}
+
       .jumbotron h1{width:58%;}
 
     }
 
     }
 
       @media (max-width: 768px){
 
       @media (max-width: 768px){

Revision as of 07:41, 16 October 2016

Proof Of Concept

Appropriate attribution is half of success.

Experiment

Entire ciruct test

In our experiment, for testing our circuit more conveniently, we use the cfp gene and gfp gene to replace the responder and lysis respectly when constructing our gene circuit. Because the fluorescence is easier to be measured. Also, the strength of the fluorescence can stand for the strength of the corresponding genes' expression.

When AHL around the engineered bacteria is at low concentration, AHL can't be sensed by the Lux sensor part, and the cI gene is expressed at the wild-type levels. While lacI-2 and GFP are repressed since the CI protein is at high concentration. When AHL around the engineered bacteria reach a certain concentration, AHL can be detected by the sensor part. So the promoter LuxpR is activated, CFP and lacI-1 express. So the CFP protein can be tested by fluorescence measurement. And the lacI-1' product-LacR protein can inhibit the expression of cI gene. As a result, lacI-2 and GFP can express. The LacI-2' product can further repress cI gene's expression, and GFP protein will be measured by fluorescence detection to represent the level of lysis gene expression. The time interval of firstly detecting CFP' fluorescence and firstly detecting GFP' fluorescence can represent the delay time which the switch give for responder expression.

Result

1. The expression of our gene circuits

The expression of our gene circuits could be induced by AHL. When AHL is detected by the sensor part, the promoter LuxpR is activated, CFP and lacI-1 express.

The E. coli containing our gene circuits were divided into experimental group and blank control group. We added AHL into experimental group, measuring OD-600 of the E. coli and the fluorescence intensity of CFP and GFP for 900 minutes. The OD-600 showed the growth trend of the E. coli. The fluorescence intensity data and time showed the relationship of the expression of CFP and GFP to time.

Figure 1.1 The OD-600 of E. coli and its variety with time on the AHL of 20ng/μL and 0ng/μL.

The growth curves were similar. It was suggested that they had the same growth trend and the difference of fluorescence intensity was mainly from the expression of circuits.

Figure 1.2 The mean CFP and GFP intensity in E. coli and its variety with time on the AHL of 20ng/μL and 0ng/μL.

The experimental group and control group all were significantly different. The fluorescence intensity between the two groups were compared with independent t-test. It showed that the cyan and green fluorescence intensity in experimental group started to be different from those in control group in the 240th minute and 420th minute separately. The expression time of GFP was delayed to 180mins post CFP.

2. AHL Gradient Induction

We found that the concentration of AHL could have an effect on the expression of our circuits, in order to find out the relation, we decided to make a AHL gradient induction for our gene circuits. As shown in the graph, the fluorescence intensity data and time on the AHL of different concentrations show the the relationship of AHL to the expression of our circuits.

Figure 1.3 The mean CFP intensity in E. coli and its variety with time on the AHL of 1ng/μL, 10ng/μL, 20ng/μL and 0ng/μL.

AHL would influenced the expression of CFP directly. The higher the concentration of AHL was, the more CFP expressed in a range from 1ng/μL to 10ng/μL. In addition, the expression time of CFP was different. T-test showed that CFP in 10ng/μL and 20ng/μL group started to express in the 240th minute and it started to express in the 480th minute in 1ng/μL group.

Figure 1.4 The mean GFP intensity in E. coli and its variety with time on the AHL of 1ng/μL, 10ng/μL, 20ng/μL and 0ng/μL.

The concentration of AHL didn't influence the expression of GFP directly. The GFP in 10ng/μL and 20ng/μL group started to express in the 420th minute and it had no significant difference in expression level. The GFP didn’t express in 1ng/μL group.

Responder test

Model

Abstract: E. coli has been modified in order to produce small interfering RNA molecules (siRNA) to attack the resistance genes that become available for nearby bacteria, when a resistant strain goes through lysis. However, a regulation system should be applied so the productive cells do not become a problem. A circuit has been implemented, this circuit regulates the transcription of said siRNAs among the production of toxins that could lead to death of any nearby cell, including the siRNA producers themselves. Taking in consideration all these elements a mathematical model using differential equations has been designed to approach the number of viable cells within the toxic environment once the circuit has been activated.

Keywords: siRNA, toxin, mRNA, transcription, translation, exponential

Rate of production of toxin mRNA and toxin molecules

The toxins are produced through a gene in the circuit, so the production of toxin depending on the translation of the mRNA produces, which also depends on transcription variables. A differential equation was taken in consideration in order to measure the mRNA rate production.

$$\frac{dr}{dt} = K - D_{r}*r$$

\(K\) represents the transcription rate, which is a constant parameter, in molecules of RNA per second, \(D_{r}\) is the degradation rate expressed in units per seconds as it depends on how much mRNA is actually present on the cell, expressed by the variable “\(r\)” in units of molecules of mRNA produced.

$$r = \frac{K}{D_{r}}+Ce^{-D_{r}t}$$

Afterwards, the original mRNA production equation was obtained through the integration of the differential equation by linear integration.

Establishing a new condition \(r(0)=R_{0}\) then we can complete the system. Substituting we obtain:

$$R_{0} = \frac{K}{D_{r}}+C \rightarrow C = R_{0}-\frac{K}{D_{r}}$$

We finally obtain the equation for mRNA synthesis:

$$r(t)= \frac{K}{D_{r}}+(R_{0}-\frac{K}{D_{r}})e^{-D_{r}t}$$

To simplify the model we assume that each mRNA molecule will be transcribed in order to produce one molecule of protein. This means that each molecule of mRNA produced is equivalent to one molecule of active toxin that in certain concentration can kill bacteria. We assume a fast degradation of mRNA inside the system and almost no significant dispersion of the molecules, so this ratio can be respected.

Productive cells inside the system

The number of E. coli cells within the system can not be considered as a constant, neither as a linear function due to the exponential nature of cellular division. Given normal growth of the cell culture, it is possible to describe the number of E. coli cells as a function of time. Nevertheless, the E. coli concerning the model produces toxins that alter the number of E. coli cells itself, also there is a natural decrease in the number of bacteria caused by cellular death, resulting in the addition of new terms to the traditional equation. An equation considering all the previous factors is shown below.

$$C(t)=C_{0}e^{Gt}(1-D)-M*C$$
$$C(t)=C_{0}e^{Gt}(1-D)-M*C_{0}e^{Gt}(1-D)$$
$$C(t)=C_{0}e^{Gt}(1-D)(1-M)$$

Where:
\(C_{0}\): Initial population [=] cells
\(C\): Total productive cells [=] cells
\(D\): Cell fraction ( \(C_{0}e^{Gt}(1-D)\) ) ​that dies naturally [=] adimensional < 1
\(G\): Rate of cells that are born [=] 1/second
\(M\): Fraction of viable cells ( \(C_{0}e^{Gt}(1-D)\) ) that lyse due to toxins [=] adimensional < 1
\(D\) and \(M\) are not related, therefore they are independent. It is also important to notice that we assumed that the number of toxins equals the number of mRNAs transcribed, so our previous defined function \(r(t)\) can help us to determine the constant \(M\).

\(r(t) < X\): no cell death, \(M = 0\)

\(r(t) \geqslant X\): cell death occurs, \(M=\frac{T r(t)}{(1-D)C_{0}e^{Gt}}\) < 1

where \(X\) is the critical toxin production rate. Also in order to determine the \(M\) term we need other variables.

\(M=\frac{T r(t)}{(1-D)C_{0}e^{Gt}}\), where:

\(r(t)\) = toxins produced in time \(t\) [=] (molecule of toxin)
\(T\) = cells that die by unit of toxin [=] (cells / molecule of toxin)
\(M\) = dead cells / total viable cells [=] adimensional

for \(r(t) < X\): \(C(t) = C_{0}e^{Gt}(1 − D)\) (3)

for \(r(t) \geqslant X\): \(C(t)= C_{0}e^{Gt}(1-D)(1-\frac{T r(t)}{(1-D)C_{0e^{Gt}}})\) (4) then

\(C(t) = C_{0}e^{Gt}(1 − D) − T r(t)\) (5)

The coupled equations yield the following final modelling:

$$C(t)=C_{0}e^{Gt}(1-D)- T[\frac{K}{D_{r}}+(R_{0}-\frac{K}{D_{r}})e^{-D_{r}t}]$$

Discussion

Reference

Name: XMU-China School: Xiamen University


Address: Xiamen University, No. 422, Siming South Road, Xiamen, Fujian, P. R. China 361005

Name: XMU-China School: Xiamen University


Address: Xiamen University, No. 422, Siming South Road, Xiamen, Fujian, P. R. China 361005