Team:USTC CHINA/Modeling/KillSwitch

From 2013.igem.org

(Difference between revisions)
 
(10 intermediate revisions not shown)
Line 88: Line 88:
<h1>Introduction</h1>
<h1>Introduction</h1>
<p align="justify">
<p align="justify">
-
To be more user-friendly, 4# circuit contains a reporting system. After melting in water, the spores germinate and express blue pigment protein to report the best using time.</br>
+
To be more user-friendly, 4# circuit contains a reporting system. After melting in water, the spores germinate and express blue pigment protein to report the optimal using time.</br>
-
Meanwhile, 4# circuit could also ensure biosafety. Since other circuit do not have self-killing device, 4# engineered bacterial should kill all the bacterial after using.</p>
+
Meanwhile, 4# circuit could also ensure biosafety. Since other circuits do not contain self-killing device, 4# engineered bacterial would kill all the bacterial after being used.</p>
</div>
</div>
Line 96: Line 96:
<h1>Designing of the suicide system</h1>
<h1>Designing of the suicide system</h1>
<p align="justify">
<p align="justify">
-
We design a circuit of killing switch based on its endogenous genetic system.</br>
+
We designed a circuit of killing switch based on its endogenous genetic system.</br>
-
In B.subtilis, when stationary phase begins, the environmental pressure increases and nutrition becomes limited, forcing B.subtilis to produce spores. Then the community will be divided into two different parts. One is trying to kill the others for enough nutrient, delaying the production of spores and obtaining competitive advantages. Killing is mediated by the exported toxic protein SdpC. SdpI will appear on the membrane surface to avoid itself from being damaged. SdpI could bind free SdpC and autopressor SdpR, removing SdpR’s inhibition against I and R, producing more SdpI to offset SdpC and finally guaranteeing the subgroup alive, thereby delaying the spores production.</br>
+
In <i>B.subtilis</i>, when stationary phase begins, the environmental pressure increases and nutrition becomes limited, forcing <i>B.subtilis</i> to produce spores. Then the community will be divided into two different parts. One is trying to kill the other one for adequate nutrient, delaying the production of spores and obtaining competitive advantages. Killing is mediated by the exported toxic protein SdpC. SdpI will appear on the membrane surface to avoid itself from being damaged. As a membrane protein, SdpI could bind free SdpC and autopressor SdpR, removing SdpR’s inhibition against I and R, producing more SdpI to offset SdpC and finally guaranteeing the subgroup alive, thereby delaying the spores production.</br>
<img src="https://static.igem.org/mediawiki/2013/2/2b/Reporter_3.png" width="500" height="350">
<img src="https://static.igem.org/mediawiki/2013/2/2b/Reporter_3.png" width="500" height="350">
-
<p>We transfer SdpC, which is fused by promoter SdpI/R into high copy plasmids, to damage the balance of the system and kill whole colony. When SdpC appears, SdpI on the membrane will bind free SdpC and adsorb SdpR to cease its inhibition against SdpI P/R, trying to produce more SdpI. At the same time, it will activate the promoter SdpR/I in our circuit and generate more SdpC. The system would fall into an infinite loop, and according to our modeling ,the amount of SdpC increases beyond the ability of SdpI. Thus, the cells with protection mechanism will finally collapse due to too much SdpC. All above forms the killing device.
+
<p>We transfered SdpC, which is fused with promoter SdpI/R into high copy plasmids, to damage the balance of the system and kill whole colony. When SdpC appears, SdpI on the membrane will bind free SdpC and adsorb SdpR to cease its inhibition against SdpI P/R, trying to produce more SdpI. At the same time, it will activate the promoter SdpR/I in our circuit and generate more SdpC. The system would fall into an infinite loop, and according to our modeling, the amount of SdpC can increase beyond the accommodation of SdpI. Thus, the cells with protection mechanism will finally collapse due to too much SdpC. All above forms the killing device.
-
We also designed a test circuit, containing promotor grac and sdpABC only, to check the ability of SdpC.</p>
+
We also designed a test circuit, which contains promotor grac and sdpABC solely, to check the toxicity of SdpC.</p>
<div></br></br>
<div></br></br>
<p align="justify">
<p align="justify">
-
There are both positive and negative feedback loops in this process. On the one hand, SdpI is unable to sequestrate the autorepressor, SdpR, until it captures the toxin, SdpC. The accumulation of SdpC will thus facilitate SdpI to capture more SdpR and thereby relieve the repression of SdpR, stimulating the expression of itself. This is the positive feedback loop which leads to the increasing accumulation of SdpC and finally the death of the bacteria. On the other hand, the removal of SdpR also enhance the expression of SdpI and accelerate the sequestration of SdpC, which forms a negative feedback loop whose effects contradict the positive feedback loop. However, since the copy number of SdpC is much higher, it is believed that the positive loop is strong enough to outweigh the negative one, which guarantees this mechanism will finally leads to collapse instead of equilibrium.
+
There are both positive and negative feedback loops in this process. On the one hand, SdpI is unable to sequestrate the autorepressor, SdpR, until it captures the toxin, SdpC. The accumulation of SdpC will thus facilitate SdpI to capture more SdpR and thereby relieve the repression of SdpR, stimulating the expression of itself. This is the positive feedback loop, which leads to the increasing accumulation of SdpC and the death of the colony. On the other hand, the removal of SdpR also enhances the expression of SdpI and accelerates the sequestration of SdpC, which forms a negative feedback loop whose effects contradict the positive feedback loop. However, since the copy number of SdpC is much higher, we believe that the positive loop is strong enough to outweigh the negative one, which guarantees this circuit will finally leads to collapse instead of equilibrium.
</div>
</div>
<div>
<div>
<h1>The ODE model of singular cells</h1>
<h1>The ODE model of singular cells</h1>
-
There is no denying fact that the essential goal of engineered bacterias who carry this so called “suicide” locus itself is to kill their siblings rather than themselves to ensure the survival of themselves. Surely they can kill their siblings, but can they finally eliminate themselves, as we expect? The trivial experiment protocol and huge uncertainty had put off our experiment, and as expected, we failed to achieve the construction of complete reporter system in our laboratory. Fortunately, we could resort to mathematical models to verify the validity of this locus theoretically.  
+
There is no denying fact that the essential goal of engineered bacterias who carry this so called “suicide” circuit itself is to kill their siblings rather than themselves and ensure the survival of themselves. Surely they can kill their siblings, but can they finally eliminate themselves, as we expect? The trivial experiment protocol and huge uncertainty had put off our experiment, and as expected, we failed to achieve the construction of complete reporter system in our limited laboratory. Fortunately, we could resort to mathematical models to verify the validity of this circuit theoretically.  
There are six independent variables in individual cells, and theoretically if the initial conditions are fixed, all of them will be the univariate functions of time. The following table illustrates the mark and meaning of each variable:</br></br>
There are six independent variables in individual cells, and theoretically if the initial conditions are fixed, all of them will be the univariate functions of time. The following table illustrates the mark and meaning of each variable:</br></br>
<div>To construct reasonable ordinary differential equation (ODE) model to describe and predict the operation of the suicide system, we followed the law of mass action, one basic law of chemistry and biology.</br>
<div>To construct reasonable ordinary differential equation (ODE) model to describe and predict the operation of the suicide system, we followed the law of mass action, one basic law of chemistry and biology.</br>
-
Taken as a statement about kinetics, the law states that the rate of an elementary reaction (a reaction that proceeds through only one transition state, which is one mechanistic step) is proportional to the product of the concentrations of the participating molecules. In modern chemistry this is derived using statistical mechanics. Despite the complicated chemical reactions involved in the process of transcription and translation, it is common and logically sound to view the expression of one particular gene as an elementary reaction and assume the repression effects of the protein itself encodes and the repressor are both linear.</br>
+
Taken as a statement about kinetics, the law states that the rate of an elementary reaction (a reaction that proceeds through only one transition state, which is one mechanistic step) is proportional to the product of the concentrations of the participating molecules. In modern chemistry this is derived using statistical mechanics. Despite the complicated chemical reactions involved in transcription and translation, it is common and logically sound to view the expression of one particular gene as an elementary reaction and assume the repression effects of the protein itself encodes and the repressor are both linear.</br>
According to the law of mass action, we got six independent differential equations of the variables:</br>
According to the law of mass action, we got six independent differential equations of the variables:</br>
<img src="https://static.igem.org/mediawiki/2013/e/e8/For1%281%29.png"></br></br></br>
<img src="https://static.igem.org/mediawiki/2013/e/e8/For1%281%29.png"></br></br></br>
Line 215: Line 215:
<h2>Discussions on the constants</h2>
<h2>Discussions on the constants</h2>
<p align="justify">
<p align="justify">
-
All the constants given above is steady and theoretically measurable when all the conditions are constant. For example, we could measure k<sub>0</sub> by constructing a new engineered bacteria, which contains the gene encoding SdpC and marker gene alone and observing the influence of the concentration of SdpC on its expression. Yet any modification on genome is notoriously time-consuming, which inhibited us from measuring them in person. We also looked up oceans of papers to confer their approximate ranges, but almost all papers are too fragmental to afford any valid information. Therefore, we decided to assume all these constant according to our limited information and make a qualitative analysis instead of quantifiable analysis. All units and dimensions were temporarily ignored. In other words, our model aims at justifying the validity of this suicide mechanism rather than predicting the exact time or any other parameters of the system.
+
All the constants given above is steady and theoretically measurable when all the conditions are constant. For example, we could measure k<sub>0</sub> by constructing a new engineered bacteria, which contains the gene encoding SdpC and marker gene alone, and observing the influence of the concentration of SdpC on its own expression. Yet any modification on genome is notoriously time-consuming, which inhibited us from measuring them in our tiny laboratory. We also looked up oceans of papers to confer their approximate ranges, but almost all papers are too fragmental to provide any valid information. Therefore, we decided to assume all these constant according to our limited information and make a qualitative analysis instead of quantifiable analysis, while in fact the latter one is impossible. All units and dimensions were temporarily ignored. In other words, our model aims at justifying the validity of this suicide mechanism rather than predicting the exact time or any other parameters of the system.
Despite the fact that we have hardly any accurate data on these constants, there are some limitations that we extrapolated from known information before we further explore this model:</br>
Despite the fact that we have hardly any accurate data on these constants, there are some limitations that we extrapolated from known information before we further explore this model:</br>
Line 221: Line 221:
<li>k<sub>0</sub>>>k<sub>4</sub>≈k<sub>7</sub>: k<sub>0</sub>,k<sub>4</sub> and k<sub>7</sub> represent the normal expression rate of SdpC, SdpR and SdpC separately, and the copy number of SdpC is much larger than that of SdpR and SdpI, whereas the value of the latter two is approximately equal;</li>
<li>k<sub>0</sub>>>k<sub>4</sub>≈k<sub>7</sub>: k<sub>0</sub>,k<sub>4</sub> and k<sub>7</sub> represent the normal expression rate of SdpC, SdpR and SdpC separately, and the copy number of SdpC is much larger than that of SdpR and SdpI, whereas the value of the latter two is approximately equal;</li>
<li>k<sub>2</sub>>>k<sub>9</sub>: the existence of free SdpR represses the expression of both SdpI and SdpC, and similarly, since the copy number of SdpC is much higher, we expected the repression effect was stronger accordingly;</li>
<li>k<sub>2</sub>>>k<sub>9</sub>: the existence of free SdpR represses the expression of both SdpI and SdpC, and similarly, since the copy number of SdpC is much higher, we expected the repression effect was stronger accordingly;</li>
-
<li>k<sub>10</sub>>>k<sub>3</sub>,k<sub>6</sub>:it is hard to predict the value of k<sub>3</sub> and k<sub>8</sub>, yet we suppose both of them is much smaller than k<sub>10</sub> because SdpI is a kind of membrane protein inherently, and rarely exists as free protein</li>
+
<li>k<sub>10</sub>>>k<sub>3</sub>,k<sub>6</sub>:it is hard to predict the value of k<sub>3</sub> and k<sub>8</sub>, but we supposed both of them were much smaller than k<sub>10</sub> because SdpI is a membrane protein inherently, and rarely exists as free protein;</li>
-
<li>The primary values of all the six variables are very small or strictly zero. We expect it as the most logical initial status. If the primary value of any variable is relatively large, the suicide mechanism may not run normally</li>
+
<li>The primary values of all the six variables are very small or strictly zero. We expect this as the most logical initial status. If the primary value of any variable is relatively large, the suicide mechanism may not run normally.</li>
</ol>
</ol>
</p>
</p>
Line 228: Line 228:
<h2>Stimulation and discussion</h2>
<h2>Stimulation and discussion</h2>
<p align="justify">
<p align="justify">
-
Simple and rough as the above model is, it does theoretically sound. To test the validity of this model, we first tried to get analytic solution of the ODE set. If this analytic solution exists, we could further investigate the interaction among those variables, and draw some phase planes to get accurate and mathematically perfect description of this model. Unfortunately but expectedly, the existence of analytic solution was negated by MATLAB, and we had to assume groups of values for these constants in advance and analyze the arithmetic solutions instead. These arithmetic solutions not only justified this mechanism is effective enough to commit cell suicide but also indicated some unexpected, or even weird results that beyond our wildest imagination. There are two possibility account for the unexpected results: our model is too rough to include some assignable factor; or there are some implicit but objective limitation inside model, which may be substantiate by later experiments or papers.</br>
+
Simple and rough as the above model is, it does theoretically sound. To test the validity of this model, we first tried to get analytic solution of the ODE set. If this analytic solution exists, we could further investigate the interaction among those variables, and draw some phase planes to get accurate and mathematically perfect descriptions of this model. Unfortunately but expectedly, the existence of analytic solution was negated by MATLAB, and we had to assume groups of values for these constants in advance and analyze the arithmetic solutions instead. These arithmetic solutions not only justified this mechanism is effective enough to commit cell suicide but also indicated some unexpected, or even weird results that beyond our wildest imagination. There are two possibility accounting for the unexpected results: our model is too rough to include some assignable factor; or there are some implicit but objective limitation inside this model, which may be substantiate by later experiments or papers.</br>
When we explored the arithmetic solutions of this ODE set, we received nearly one hundred warnings from MATLAB and for many times our most powerful computer ran out of its 8GB memory, but sometimes we can receive the solution within seconds. We had adjusted our parameters for several times before we got our first solution. Here is the values of parameters for this group, and the graph of arithmetic solutions is also given:</p>
When we explored the arithmetic solutions of this ODE set, we received nearly one hundred warnings from MATLAB and for many times our most powerful computer ran out of its 8GB memory, but sometimes we can receive the solution within seconds. We had adjusted our parameters for several times before we got our first solution. Here is the values of parameters for this group, and the graph of arithmetic solutions is also given:</p>
<div style="margin-top:20px;">
<div style="margin-top:20px;">
Line 275: Line 275:
</div></br></br>
</div></br></br>
<img class="linegraph" src="https://static.igem.org/mediawiki/2013/2/26/Suicide1.png" width="450" height="350" /></br></br>
<img class="linegraph" src="https://static.igem.org/mediawiki/2013/2/26/Suicide1.png" width="450" height="350" /></br></br>
-
<div>At the first glance this graph seemed fine. Initially the concentration of SdpC decreased slightly due to the capturing of SdpI and the repression of float SdpR, but gradually the positive feedback loop works, and C<sub>f</sub> increases rapidly. But when we turned our attention to the curves of other parameters, things seemed not so perfect:</div></br>
+
<div>At the first glance this graph seemed fine. Initially the concentration of SdpC decreased slightly due to the capturing of SdpI and the repression of float SdpR, but gradually the positive feedback loop works, and C<sub>f</sub> increases rapidly. But when we turned our attention to the curves of other parameters, things seemed not so platonic:</div></br>
<div><img class="linegraph" src="https://static.igem.org/mediawiki/2013/c/cf/Suicide2.png" width="450" height="350" /></br></br>
<div><img class="linegraph" src="https://static.igem.org/mediawiki/2013/c/cf/Suicide2.png" width="450" height="350" /></br></br>
-
<div></br></br></br>The curve of I<sub>m</sub>, C<sub>i</sub> and R<sub>i</sub> contradicted our common sense severely. First, I<sub>m</sub>>C<sub>i</sub>>R<sub>i</sub> is expected to be tenable all the time, which precludes the intersects among the three curves; Second, there is no mechanism in this system that could decrease their concentration, and all of them are increasing function; Third and most serious, never will them be negative, as they represent the concentration of real substances.
+
<div></br></br></br>The curve of I<sub>m</sub>, C<sub>i</sub> and R<sub>i</sub> contradicted our common sense severely. First, I<sub>m</sub>>C<sub>i</sub>>R<sub>i</sub> is expected to be tenable all the time, which precludes the intersects among the three curves; Second, there is no mechanism in this system that could decrease their concentration, and all of them are expected to be increasing function; Third and most serious, never will them be negative, as they represent the concentration of real substances.
Then we adjusted the parameters slightly for several times. To eliminate those absurd curves, we reconsidered some assumptions.
Then we adjusted the parameters slightly for several times. To eliminate those absurd curves, we reconsidered some assumptions.
Here we listed another representative group of parament values and relative graph:</div></div></br></br>
Here we listed another representative group of parament values and relative graph:</div></div></br></br>
Line 368: Line 368:
</table></br></br>
</table></br></br>
-
<div style="float:left; margin-left:330px;margin-top:-350px;width:270px;align="justify;">In this group, we gave up one former assumption and set k<sub>2</sub> equal to k<sub>9</sub>. We also gave positive values to I<sub>m</sub>, C<sub>i</sub> and R<sub>i</sub>, which were considered to be zero at first. And by groups of stimulations we realized the value of k<sub>2</sub> does matter, as the derivative of C<sub>f</sub> only increased slightly as k<sub>2</sub> lowers, and the positive values failed to avoid the weird phenomenon in the latter three curves.<br/><br/>
+
<div style="float:left; margin-left:330px;margin-top:-350px;width:270px;align="justify;">In this group, we gave up one former assumption and set k<sub>2</sub> equal to k<sub>9</sub>. We also gave positive values to I<sub>m</sub>, C<sub>i</sub> and R<sub>i</sub>, which were considered zero at first. And by groups of stimulations we realized the value of k<sub>2</sub> does matter, as the derivative of C<sub>f</sub> only increased slightly as k<sub>2</sub> lowers, and the positive values failed to avoid the weird phenomenon in the latter three curves.<br/><br/>
-
We also found that however we adjusted the primary value of I<sub>f</sub> and other parameters, If dropped into approximately zero extremely rapidly at the initial stage and remained balanced, which might account for why the derivatives of the latter curves were abnormally negative. Thus we modified another assumption and increased k<sub>7</sub>. Here is another group of values and corresponding graph:</div></div></br></br>
+
We also found that however we adjusted the primary value of I<sub>f</sub> and other parameters, I<sub>f</sub> dropped into approximately zero extremely rapidly at the initial stage and remained balanced, which might account for why the derivatives of the latter curves were abnormally negative. Thus we modified another assumption and increased k<sub>7</sub>. Here is another group of values and the corresponding graph:</div></div></br></br>
Line 375: Line 375:
<p align="justify">
<p align="justify">
-
Although the derivative of Im is not seriously positive constantly, the three latter curves seemed much more reasonable. Hence, we extrapolated although SdpI and SdpR share the same promoter, the expression of SdpI must much faster than SdpR to ensure successful “suicide.Additionally, the increase of k<sub>7</sub> also represses SdpC, and hence the copy number of SdpC must be larger.
+
Although the derivative of I<sub>m</sub> is not seriously positive constantly, the three latter curves seemed much more reasonable. Hence, we extrapolated that although SdpI and SdpR share the same promoter, the expression of SdpI must much faster than SdpR to ensure successful “suicide”. Additionally, the increase of k<sub>7</sub> also represses SdpC, which requires the copy number of SdpC must be larger.
We kept all other parameters constant and gradually augmented k<sub>0</sub>. The larger k<sub>0</sub>, the more perfect the curve seemed, and here are the values table and graph where k<sub>0</sub> equals 400, 80 times larger than k<sub>4</sub>.</br></p>
We kept all other parameters constant and gradually augmented k<sub>0</sub>. The larger k<sub>0</sub>, the more perfect the curve seemed, and here are the values table and graph where k<sub>0</sub> equals 400, 80 times larger than k<sub>4</sub>.</br></p>
<table border="1" align="center" frame="box">
<table border="1" align="center" frame="box">
Line 467: Line 467:
</table></br></br>
</table></br></br>
<img class="linegraph" src="https://static.igem.org/mediawiki/2013/1/19/Suicide7.png"></br></br>
<img class="linegraph" src="https://static.igem.org/mediawiki/2013/1/19/Suicide7.png"></br></br>
-
Wired but not surprising, there were intersects among the latter three curve, and C<sub>f</sub> decreases continually when it is negative. We continued to try groups of these parameters, and this is the best one where we increased the primary concentration of SdpC and the normal expression rate of SdpI.
+
Wired but not surprising, there were intersects among the latter three curve, and C<sub>f</sub> decreases continually when it is negative. We continued to try groups of these parameters, and this is the best one in which we increased the primary concentration of SdpC and the normal expression rate of SdpI.
<table border="1" align="center" frame="box">
<table border="1" align="center" frame="box">
<tr>
<tr>
Line 514: Line 514:
<img class="linegraph" src="https://static.igem.org/mediawiki/2013/e/eb/Suicide9.png"></br></br>
<img class="linegraph" src="https://static.igem.org/mediawiki/2013/e/eb/Suicide9.png"></br></br>
<div></br></br></br></br></br></br></br></br></br></br></br></br></br>In spite of minimal abnormal phenomenon (C<sub>f</sub> was negative in later stage), this graph roughly testified that in wild bacterial the concentration of float SdpC will drop to nearly zero quickly.
<div></br></br></br></br></br></br></br></br></br></br></br></br></br>In spite of minimal abnormal phenomenon (C<sub>f</sub> was negative in later stage), this graph roughly testified that in wild bacterial the concentration of float SdpC will drop to nearly zero quickly.
-
In sum, the ODE model of singular cells indicate following results:</br>
+
In sum, the ODE model on singular cell indicates following results:</br>
<ol>
<ol>
-
<li>The character of free SdpC is most affected by k<sub>0</sub>, if the copy number of SdpC is large enough, it is theoretically reasonable to commit suicide;</li>
+
<li>The concentration of free SdpC is most affected by k<sub>0</sub>, if the copy number of SdpC is large enough, it is theoretically reasonable to commit suicide;</li>
-
<li>The influence of the value of Imax and k<sub>2</sub> is much limited;</li>
+
<li>The influence of the value of I<sub>max</sub> and k<sub>2</sub> is much limited;</li>
<li>The amount of free SdpI is always near zero;</li>
<li>The amount of free SdpI is always near zero;</li>
<li>SdpC will not increase limitlessly however we transform parameters;</li>
<li>SdpC will not increase limitlessly however we transform parameters;</li>
-
<li>To ensure the success of suicide, it is required k<sub>0</sub>>>k<sub>4</sub>>>k<sub>7</sub>;</li></div>
+
<li>To ensure successful suicide, it is required k<sub>0</sub>>>k<sub>4</sub>>>k<sub>7</sub>;</li></div>
-
<p style="width:300px; margin-left:300px; align="justify"">The last conclusion was our biggest windfall, and we have verified the validity of this suicide mechanism in math. On the one hand, if further experiments proven #4 engineered bacteria will kill both siblings and themselves, it is highly like that the expression rate SdpI is much larger than SdpR even if they share the same promoter; on the other hand, if #4 engineered bacteria are not able to commit suicide, we can try to boost the expression of SdpI to adjust the bacteria.</p>
+
<p style="width:300px; margin-left:300px; align="justify"">The last conclusion was our biggest windfall, and we have verified the validity of this suicide mechanism in math. On the one hand, if further experiments proven #4 engineered bacteria will kill both siblings and themselves, it is highly like that the expression rate SdpI is much larger than SdpR even if they share the same promoter; on the other, if #4 engineered bacteria are not able to commit suicide, we can try to boost the expression of SdpI to adjust the kill switch.</p>
</ol>
</ol>
Line 527: Line 527:
<h1>Discussion on colonies</h1>
<h1>Discussion on colonies</h1>
-
In reality, the engineered bacteria aims at killing its siblings instead of itself, and at first almost all toxin SdpC will be secreted outside the bacteria. Assume the diffusion of toxin among cells comply with diffusion law, that is, the diffusion rate is proportionate with the gradient of concentration. Further assume the death concentration of SdpC is same to all bacteria expect those who contain this locus, the average life expectancy is bacteria will hinge on the rate and distribution of engineered bacteria, and the distribution of life expectancy of bacteria is similar to that of average free path of gas molecules.
+
In reality, the engineered bacteria aims at killing its siblings instead of itself, and at first almost all toxin SdpC will be secreted outside the bacteria. We assume the diffusion of toxin among cells comply with diffusion law, that is, the diffusion rate is proportionate with the gradient of concentration. Further we assume the death concentration of SdpC is same to all bacteria expect those who contain this locus, the average life expectancy is bacteria will hinge on the rate and distribution of engineered bacteria, and the distribution of life expectancy of bacteria is similar to that of average free path of gas molecules.
-
As long the coefficient of diffusion is large enough, any engineered bacterias, no matter how few, is adequate to devastate the whole colony. Alike to the average free path of thin gas, the average suicide time of the whole reporter system is inversely proportional with the square root of the rate of engineer bacteria containing this locus.
+
As long the coefficient of diffusion is large enough, any engineered bacteria, no matter how few, is adequate to devastate the whole colony. Alike to the average free path of thin gas, the average suicide time of the whole reporter system is inversely proportional with the square root of the rate of engineer bacteria containing this circuit.
</div>
</div>
<h1>References</h1>
<h1>References</h1>
<p align="justify">
<p align="justify">
-
Parallel pathways of repression and antirepression governing the transition to stationary phase in Bacillus subtilis
+
Parallel pathways of repression and antirepression governing the transition to stationary phase in <i>Bacillus subtilis</i>
AV Banse, A Chastanet, L Rahn-Lee…,PNAS ,2008 </p>
AV Banse, A Chastanet, L Rahn-Lee…,PNAS ,2008 </p>
</div>
</div>

Latest revision as of 12:59, 28 October 2013

Introduction

To be more user-friendly, 4# circuit contains a reporting system. After melting in water, the spores germinate and express blue pigment protein to report the optimal using time.
Meanwhile, 4# circuit could also ensure biosafety. Since other circuits do not contain self-killing device, 4# engineered bacterial would kill all the bacterial after being used.

Designing of the suicide system

We designed a circuit of killing switch based on its endogenous genetic system.
In B.subtilis, when stationary phase begins, the environmental pressure increases and nutrition becomes limited, forcing B.subtilis to produce spores. Then the community will be divided into two different parts. One is trying to kill the other one for adequate nutrient, delaying the production of spores and obtaining competitive advantages. Killing is mediated by the exported toxic protein SdpC. SdpI will appear on the membrane surface to avoid itself from being damaged. As a membrane protein, SdpI could bind free SdpC and autopressor SdpR, removing SdpR’s inhibition against I and R, producing more SdpI to offset SdpC and finally guaranteeing the subgroup alive, thereby delaying the spores production.

We transfered SdpC, which is fused with promoter SdpI/R into high copy plasmids, to damage the balance of the system and kill whole colony. When SdpC appears, SdpI on the membrane will bind free SdpC and adsorb SdpR to cease its inhibition against SdpI P/R, trying to produce more SdpI. At the same time, it will activate the promoter SdpR/I in our circuit and generate more SdpC. The system would fall into an infinite loop, and according to our modeling, the amount of SdpC can increase beyond the accommodation of SdpI. Thus, the cells with protection mechanism will finally collapse due to too much SdpC. All above forms the killing device. We also designed a test circuit, which contains promotor grac and sdpABC solely, to check the toxicity of SdpC.



There are both positive and negative feedback loops in this process. On the one hand, SdpI is unable to sequestrate the autorepressor, SdpR, until it captures the toxin, SdpC. The accumulation of SdpC will thus facilitate SdpI to capture more SdpR and thereby relieve the repression of SdpR, stimulating the expression of itself. This is the positive feedback loop, which leads to the increasing accumulation of SdpC and the death of the colony. On the other hand, the removal of SdpR also enhances the expression of SdpI and accelerates the sequestration of SdpC, which forms a negative feedback loop whose effects contradict the positive feedback loop. However, since the copy number of SdpC is much higher, we believe that the positive loop is strong enough to outweigh the negative one, which guarantees this circuit will finally leads to collapse instead of equilibrium.

The ODE model of singular cells

There is no denying fact that the essential goal of engineered bacterias who carry this so called “suicide” circuit itself is to kill their siblings rather than themselves and ensure the survival of themselves. Surely they can kill their siblings, but can they finally eliminate themselves, as we expect? The trivial experiment protocol and huge uncertainty had put off our experiment, and as expected, we failed to achieve the construction of complete reporter system in our limited laboratory. Fortunately, we could resort to mathematical models to verify the validity of this circuit theoretically. There are six independent variables in individual cells, and theoretically if the initial conditions are fixed, all of them will be the univariate functions of time. The following table illustrates the mark and meaning of each variable:

To construct reasonable ordinary differential equation (ODE) model to describe and predict the operation of the suicide system, we followed the law of mass action, one basic law of chemistry and biology.
Taken as a statement about kinetics, the law states that the rate of an elementary reaction (a reaction that proceeds through only one transition state, which is one mechanistic step) is proportional to the product of the concentrations of the participating molecules. In modern chemistry this is derived using statistical mechanics. Despite the complicated chemical reactions involved in transcription and translation, it is common and logically sound to view the expression of one particular gene as an elementary reaction and assume the repression effects of the protein itself encodes and the repressor are both linear.
According to the law of mass action, we got six independent differential equations of the variables:


















The following table explain the constants in the above ODE group:

Mark Meaning
Imax Mole number of free SdpI in cytoplasm.
Im Mole number of SdpI in the cell membrane.
Cf Mole number of free SdpC in cytoplasm.
Ci Mole number of SdpC captured by SdpI.
Rf Mole number of free SdpR in cytoplasm.
Ri Mole number of SdpR captured by SdpI


Mark Meaning
Imax The maximal number of SdpI than can be fixed on the cell membrane.
k0 Constant describes the normal expression rate of SdpC
k1 Constant describes the self-repression effects of SdpC
k2 Constant describes the repression of SdpR on the expression of SdpC.
k3 Constant describes the rate of SdpI capturing SdpC
k4 Constant describes the normal expression rate of SdpR
k5 Constant describes the self-repression effects of SdpR
k6 Constant describes the rate of SdpI capturing SdpR
k7 Constant describes the normal expression rate of SdpI
k8 Constant describes the self-repression effects of SdpI
k9 Constant describes the repression of SdpR on the expression of SdpI
k10 Constant describes the rate of SdpI binding to the cell membrane

Discussions on the constants

All the constants given above is steady and theoretically measurable when all the conditions are constant. For example, we could measure k0 by constructing a new engineered bacteria, which contains the gene encoding SdpC and marker gene alone, and observing the influence of the concentration of SdpC on its own expression. Yet any modification on genome is notoriously time-consuming, which inhibited us from measuring them in our tiny laboratory. We also looked up oceans of papers to confer their approximate ranges, but almost all papers are too fragmental to provide any valid information. Therefore, we decided to assume all these constant according to our limited information and make a qualitative analysis instead of quantifiable analysis, while in fact the latter one is impossible. All units and dimensions were temporarily ignored. In other words, our model aims at justifying the validity of this suicide mechanism rather than predicting the exact time or any other parameters of the system. Despite the fact that we have hardly any accurate data on these constants, there are some limitations that we extrapolated from known information before we further explore this model:

  1. k0>>k4≈k7: k0,k4 and k7 represent the normal expression rate of SdpC, SdpR and SdpC separately, and the copy number of SdpC is much larger than that of SdpR and SdpI, whereas the value of the latter two is approximately equal;
  2. k2>>k9: the existence of free SdpR represses the expression of both SdpI and SdpC, and similarly, since the copy number of SdpC is much higher, we expected the repression effect was stronger accordingly;
  3. k10>>k3,k6:it is hard to predict the value of k3 and k8, but we supposed both of them were much smaller than k10 because SdpI is a membrane protein inherently, and rarely exists as free protein;
  4. The primary values of all the six variables are very small or strictly zero. We expect this as the most logical initial status. If the primary value of any variable is relatively large, the suicide mechanism may not run normally.



Stimulation and discussion

Simple and rough as the above model is, it does theoretically sound. To test the validity of this model, we first tried to get analytic solution of the ODE set. If this analytic solution exists, we could further investigate the interaction among those variables, and draw some phase planes to get accurate and mathematically perfect descriptions of this model. Unfortunately but expectedly, the existence of analytic solution was negated by MATLAB, and we had to assume groups of values for these constants in advance and analyze the arithmetic solutions instead. These arithmetic solutions not only justified this mechanism is effective enough to commit cell suicide but also indicated some unexpected, or even weird results that beyond our wildest imagination. There are two possibility accounting for the unexpected results: our model is too rough to include some assignable factor; or there are some implicit but objective limitation inside this model, which may be substantiate by later experiments or papers.
When we explored the arithmetic solutions of this ODE set, we received nearly one hundred warnings from MATLAB and for many times our most powerful computer ran out of its 8GB memory, but sometimes we can receive the solution within seconds. We had adjusted our parameters for several times before we got our first solution. Here is the values of parameters for this group, and the graph of arithmetic solutions is also given:

k0 k1 k2 k3 k4 k5 k6 k8 k8 k9 k10 Imax Cf0 Rf0 If0 Im0 Ci0 Ri0
50 5 5 5 5 5 5 5 5 5 20 500 5 5 1 5 3 2




At the first glance this graph seemed fine. Initially the concentration of SdpC decreased slightly due to the capturing of SdpI and the repression of float SdpR, but gradually the positive feedback loop works, and Cf increases rapidly. But when we turned our attention to the curves of other parameters, things seemed not so platonic:






The curve of Im, Ci and Ri contradicted our common sense severely. First, Im>Ci>Ri is expected to be tenable all the time, which precludes the intersects among the three curves; Second, there is no mechanism in this system that could decrease their concentration, and all of them are expected to be increasing function; Third and most serious, never will them be negative, as they represent the concentration of real substances. Then we adjusted the parameters slightly for several times. To eliminate those absurd curves, we reconsidered some assumptions. Here we listed another representative group of parament values and relative graph:


k0 k1 k2 k3 k4 k5 k6 k8 k8 k9 k10 Imax Cf0 Rf0 If0 Im0 Ci0 Ri0
100 5 5 5 5 5 5 20 5 5 20 500 5 5 1 5 3 2




k0 k1 k2 k3 k4 k5 k6 k8 k8 k9 k10 Imax Cf0 Rf0 If0 Im0 Ci0 Ri0
100 5 5 5 5 5 5 20 5 5 20 500 5 5 1 5 3 2


In this group, we gave up one former assumption and set k2 equal to k9. We also gave positive values to Im, Ci and Ri, which were considered zero at first. And by groups of stimulations we realized the value of k2 does matter, as the derivative of Cf only increased slightly as k2 lowers, and the positive values failed to avoid the weird phenomenon in the latter three curves.

We also found that however we adjusted the primary value of If and other parameters, If dropped into approximately zero extremely rapidly at the initial stage and remained balanced, which might account for why the derivatives of the latter curves were abnormally negative. Thus we modified another assumption and increased k7. Here is another group of values and the corresponding graph:




Although the derivative of Im is not seriously positive constantly, the three latter curves seemed much more reasonable. Hence, we extrapolated that although SdpI and SdpR share the same promoter, the expression of SdpI must much faster than SdpR to ensure successful “suicide”. Additionally, the increase of k7 also represses SdpC, which requires the copy number of SdpC must be larger. We kept all other parameters constant and gradually augmented k0. The larger k0, the more perfect the curve seemed, and here are the values table and graph where k0 equals 400, 80 times larger than k4.

k0 k1 k2 k3 k4 k5 k6 k8 k8 k9 k10 Imax Cf0 Rf0 If0 Im0 Ci0 Ri0
400 5 5 5 5 5 5 20 5 5 20 500 5 5 1 5 3 2




Take the curve of Cf and Rf separately, the curves seemed more perfect:















In wild bacteria who are unable to produced SdpC, naturally k0 equals zero. We expected Cf would decreased gradually and finally approximate zero, and here are the corresponding table and graph:


k0 k1 k2 k3 k4 k5 k6 k8 k8 k9 k10 Imax Cf0 Rf0 If0 Im0 Ci0 Ri0
0 5 5 5 5 5 5 30 5 5 20 500 8 5 1 5 3 2




Wired but not surprising, there were intersects among the latter three curve, and Cf decreases continually when it is negative. We continued to try groups of these parameters, and this is the best one in which we increased the primary concentration of SdpC and the normal expression rate of SdpI.
k0 k1 k2 k3 k4 k5 k6 k8 k8 k9 k10 Imax Cf0 Rf0 If0 Im0 Ci0 Ri0
0 5 5 5 5 5 5 5 5 5 20 500 8 5 1 5 3 2




The curve of Cf and Rf alone:














In spite of minimal abnormal phenomenon (Cf was negative in later stage), this graph roughly testified that in wild bacterial the concentration of float SdpC will drop to nearly zero quickly. In sum, the ODE model on singular cell indicates following results:
  1. The concentration of free SdpC is most affected by k0, if the copy number of SdpC is large enough, it is theoretically reasonable to commit suicide;
  2. The influence of the value of Imax and k2 is much limited;
  3. The amount of free SdpI is always near zero;
  4. SdpC will not increase limitlessly however we transform parameters;
  5. To ensure successful suicide, it is required k0>>k4>>k7;

The last conclusion was our biggest windfall, and we have verified the validity of this suicide mechanism in math. On the one hand, if further experiments proven #4 engineered bacteria will kill both siblings and themselves, it is highly like that the expression rate SdpI is much larger than SdpR even if they share the same promoter; on the other, if #4 engineered bacteria are not able to commit suicide, we can try to boost the expression of SdpI to adjust the kill switch.

Discussion on colonies

In reality, the engineered bacteria aims at killing its siblings instead of itself, and at first almost all toxin SdpC will be secreted outside the bacteria. We assume the diffusion of toxin among cells comply with diffusion law, that is, the diffusion rate is proportionate with the gradient of concentration. Further we assume the death concentration of SdpC is same to all bacteria expect those who contain this locus, the average life expectancy is bacteria will hinge on the rate and distribution of engineered bacteria, and the distribution of life expectancy of bacteria is similar to that of average free path of gas molecules. As long the coefficient of diffusion is large enough, any engineered bacteria, no matter how few, is adequate to devastate the whole colony. Alike to the average free path of thin gas, the average suicide time of the whole reporter system is inversely proportional with the square root of the rate of engineer bacteria containing this circuit.

References

Parallel pathways of repression and antirepression governing the transition to stationary phase in Bacillus subtilis AV Banse, A Chastanet, L Rahn-Lee…,PNAS ,2008