Team:Wageningen UR/project/model

From 2014.igem.org

(Difference between revisions)
 
(81 intermediate revisions not shown)
Line 31: Line 31:
                                     </li>
                                     </li>
                                     <li class="menu-item">
                                     <li class="menu-item">
-
                                         <a href="#performance1">Performance</a>
+
                                         <a href="#performance">Performance</a>
                                         <ul class="three-item">
                                         <ul class="three-item">
                                             <!-- Change the number here to reflect the number of subheaders.-->
                                             <!-- Change the number here to reflect the number of subheaders.-->
Line 51: Line 51:
                                     <li class="menu-item">
                                     <li class="menu-item">
                                         <a href="#parameters1">Parameters</a>
                                         <a href="#parameters1">Parameters</a>
-
                                         <ul class="two-item">
+
                                         <ul class="three-item">
                                             <!-- Change the number here to reflect the number of subheaders.-->
                                             <!-- Change the number here to reflect the number of subheaders.-->
-
                                             <li><a href="#parameters41">Metabolic</a></li>
+
                                             <li><a href="#parameters41">Promoters</a></li>
 +
                                            <!-- if needed only -->
 +
                                            <li><a href="#parameters42">Cost</a></li>
 +
                                            <!-- if needed only -->
                                             <!-- if needed only -->  
                                             <!-- if needed only -->  
-
                                             <li><a href="#parameters42">Dynamic</a></li>
+
                                             <li><a href="#parameters43">Performance</a></li>
-
                                             <!-- if needed only -->                                              
+
                                             <!-- if needed only -->                                                
                                         </ul>
                                         </ul>
                                     </li>
                                     </li>
Line 69: Line 72:
<h2>Modeling Overview</h2>
<h2>Modeling Overview</h2>
<p>
<p>
-
The performance of BananaGuard as biological control agent is defined three factors: its ability to destroy fusarium upon detection, the stability and performance of the genetic circuit and its ability to compete with other micro-organisms in the soil. The ability to destroy fusarium has be tested <a href="https://2014.igem.org/Team:Wageningen_UR/project/fungal_inhibition" class="soft_link">in vitro</a>. This leaves the remaining two questions to be answered. For the performance and stability of the genetic circuit two perspectives were considered. An a priori design of the Kill-switch was modeled aiming to optimize its stability. Including different parts of the genetic circuit the optimized version was subsequently modeled over time. This allowed for the quantification and characterization of the stability and performance of the BananaGuard system.
+
The performance of BananaGuard as biological control agent is defined by <b>three</b> factors: its ability to <b>destroy or inhibit</b> <i>Fusarium oxysporum</i> upon detection, its ability to <b>compete</b> with other micro-organisms in the soil and the <b>stability</b> of the genetic circuit.</p><p>
-
Before this quantification could be done however a parameter was needed. This parameter was provided in the process of answering the competitiveness query. In total three models were made (figure 1):
+
 
-
<br/>
+
If the fungal growth inhibitors cannot destroy or inhibit <i>F. oxysporum</i>, the bananaplants will die. This aspect of BananaGuards performance is being tested <a href="https://2014.igem.org/Team:Wageningen_UR/project/fungal_inhibition" class="soft_link">in vitro</a>. This leaves two  remaining questions to be answered:</p><p> <b>Can the BananaGuard system survive in the presence of competing micro-organisms?</b></p><p> <b>Is the genetic circuit robust to biological fluctuations?</b>
-
</p>
+
 
-
<p>
+
 
-
<b>(I) Promoter Design:</b>
+
</p><p>The performance and stability of the genetic circuit is crucial for BananaGuard to work. To prevent <i>F. oxysporum</i> from functioning, the fungal growth inhibitors in BananaGuard need to be produced when the fungi are present. This requires two aspects of the system to function properly: the toggle-switch within the Kill-switch system needs to be stable and the fusaric acid dependent promoter needs to activate the system in the presence of fusaric acid.
-
Accurate balancing of the different promoters elements in the <a href="https://2014.igem.org/Team:Wageningen_UR/project/kill-switch" class="soft_link"> Kill-switch</a> is required to maintain a bi-stable system. Using statistical mechanics an estimation for the optimal promoter configuration was made by balancing the number of different repressors with number of repressor binding sites and their respective position on the promoter. The obtained results have inspired the design of promoters for the kill-switch.
+
 
-
<br/>
+
</p><p> Optimizing <a href="https://2014.igem.org/Team:Wageningen_UR/project/model#design1" class="soft_link"><b>promoter design</b></a>
-
</p>
+
constitutes the first step of BananaGuard development. Different configurations of promoter elements on the promoters of the Kill-switch were modeled using statistical mechanics. Parameters such as repressor binding strengths and repressor production rates were changed in each simulation. By scoring Kill-switch activation, the promoter element configuration with the most simulations showing a bi-stable switch were considered to be most likely to work. Including the anti-gene transfer inhibition effect on cell growth rates and the fusaric acid dependent promoter, this optimized version of the BananaGuard system was subsequently modeled over a population of dividing cells, assessing both its stability and performance in the form of a stochastic model.
-
<p>
+
 
-
<b>(II) Cost:</b> To ascertain our hosts competitiveness the rate of cell division needs to be investigated on a metabolic level. Cell growth rate is a key factor in the competition between our engineered <i>Pseudomonas putida</i> and the wildtype <i>P. putida</i> found in the rhizosphere of the banana plant roots[1,2]. The integrated genetic circuit uses metabolic resources that would otherwise be dedicated to cellular maintenance or growth. Should this divergence of resources be too large our host will be out competed. A genome scale metabolic model was extended, investigating and outlining the price that has to be paid for introducing our synthetic genes on a metabolic level. The resulting average growth rate has been used in the stochastic model.
+
 
-
</p>
+
</p><p>Running thousands of simulations incorporating <b>noise terms</b>  and <b>promoter leakiness,</b> <a href="https://2014.igem.org/Team:Wageningen_UR/project/model#performance1" class="soft_link">the robustness of the BananaGuard system</a> was quantified. Stochasticity in the system caused by cell division resulted in fractions of fast and slow growing cell populations with either stable or unstable Kill-switches. The same was considered for system performance; upon <b>induction</b> by fusaric acid a portion of the cell population activated their Kill-switch or remained in the same state. In order to make the simulations as accurate as possible a parameter for the maximum growth rate of BananaGuard was required. This parameter was obtained from a metabolic model that described the effects of competition between the BananaGuard system and micro-organisms in the rhizosphere.  
-
<p>
+
 
-
<b>(III) Stability and Performance:</b> A stochastic model was made in order to quantify and characterize the stability and performance of the introduced genetic circuit. The model takes into account growth rates predicted by the metabolic model, cell division, the <a href="https://2014.igem.org/Team:Wageningen_UR/project/gene_transfer" class="soft_link"> anti-gene transfer toxin antitoxin </a> systems effect on cell growth and the kill-switch.  
+
</p><p>To ascertain our hosts competitiveness the rate of cell division needs to be investigated on a metabolic level. Cell growth rate is a key factor in the competition between our engineered <i>Pseudomonas putida</i> and the wildtype <i>P. putida</i> found in the rhizosphere of the banana plant roots [1,2]. The integrated genetic circuit uses metabolic resources that would otherwise be dedicated to cellular maintenance or growth. Should this divergence of resources be too large our host will be outcompeted. A genome scale <a href="https://2014.igem.org/Team:Wageningen_UR/project/model#cost1" class="soft_link">metabolic model</a> was extended, investigating and outlining the cost of introducing our synthetic genes on a metabolic level.
</p>
</p>
<p>
<p>
</p>
</p>
<figure>
<figure>
-
<img src="https://static.igem.org/mediawiki/2014/5/5d/Wageningen_UR_modeling_overview.png" width="90%">
+
<img src="https://static.igem.org/mediawiki/2014/thumb/0/04/Wageningen_UR_model_bob_final_f_version.png/800px-Wageningen_UR_model_bob_final_f_version.png" width="90%">
-
<figcaption>Figure 1. Schematic overview of modeling project indicating overlapping parts</figcaption>
+
<figcaption>Figure 1. Schematic overview of the modeling project indicating overlapping parts.</figcaption>
                         </figure>
                         </figure>
<br/>
<br/>
Line 97: Line 100:
                         <h2 id="introduction1" >Introduction</h2>
                         <h2 id="introduction1" >Introduction</h2>
<p>
<p>
-
Promoter design for a stable kill-switch involves balancing of different regulatory components. Two elements, the number of repressor binding sites and the position of those binding sites with respect to the gene, control two key variables; the chance a repressor binds to the regulatory region and the affinity with which it binds. These variables, in turn, determine the effectiveness of the repressor protein [1]. <br/>
+
Promoter design for a stable <a href="https://2014.igem.org/Team:Wageningen_UR/project/kill-switch" class="soft_link">Kill-switch</a> involves balancing of different regulatory components. Two elements, the number of repressor binding sites and the position of those binding sites with respect to the gene, control two key variables: the chance a repressor binds to the regulatory region and the affinity with which it binds. These variables, in turn, determine the effectiveness of the repressor protein [1,2]. <br/>
-
It has been shown that repressor proteins are more effective when there are more repressor binding sites in front of a gene, confirming that the increased chance of repressor protein binding leads to a significant increase in its effectiveness [2]. With respect to binding affinity, it has been shown for TetR proteins that a repressor binding site positioned in front of the open reading frame is more effective than those positioned at the back [3].<br/>
+
It has been shown that repressor proteins are more effective when there are more repressor binding sites in front of a gene, confirming that the increased chance of repressor protein binding leads to a significant increase in its effectiveness [1,2]. With respect to binding affinity, it has been shown for TetR proteins that a repressor binding site positioned in front of the open reading frame is more effective than those positioned at the back [6].<br/>
-
Because of slight variations in the position of a repressor binding site or variations in experimental conditions, a range of different repressor binding strength values are reported for the same repressor (sometimes varying several orders of magnitude)[4].<br/>
+
Because of slight variations in the position of a repressor binding site or variations in experimental conditions, a range of different repressor binding strength values are reported for the same repressor (sometimes varying several orders of magnitude)[4,5].<br/>
-
This means that case specific designs cannot be modeled accurately but the likelihood of designing a functioning system can be estimated. By taking all the possible repressor binding site combinations and varying the different binding affinities (thereby simulating their position) and production rates, the robustness of the repressor binding site combinations can be determined (figure. 1). Statistical mechanics is an ideal tool because it can incorporate both repressor binding strength and binding sites. Learn more about the statistical mechanics approach or view the results.
+
This means that case specific designs cannot be modeled accurately but it is possible to estimate the likelihood that the designed system will work. By taking all the possible repressor binding site combinations and varying the different binding affinities (thereby simulating their position) and production rates, the robustness of the system can determined given the repressor binding site combinations(figure 2). Statistical mechanics is an ideal tool because it can incorporate both repressor binding strength and binding sites. Continue to learn more about the statistical mechanics approach or view the <a href="https://2014.igem.org/Team:Wageningen_UR/project/model#results1" class="soft_link">results</a>.
</p>
</p>
                         <figure>
                         <figure>
-
                             <img src="https://static.igem.org/mediawiki/2014/7/70/Wagenigen_UR_model_kill-switch.png" width="90%">
+
                             <img src="https://static.igem.org/mediawiki/2014/thumb/f/f1/Wageningen_UR_model_bob2.png/800px-Wageningen_UR_model_bob2.png" width="90%">
-
                             <figcaption><b>Figure 1.</b> The figure shows the kill-switch with arbitrarily chosen repressor binding sites. (I) Altering the position of the TetR repressor binding sites change the effectiveness; a large space of basepairs between the binding site and the LacI will result in weaker repression [8]. (II) Adding a LacI repressor binding site will increase the chance that LacI binds to the promoter region. </figcaption>
+
                             <figcaption><b>Figure 2.</b> The figure shows the kill-switch with arbitrarily chosen repressor binding sites. (I) Altering the position of the TetR repressor binding sites change the effectiveness [6]; a large space of basepairs between the binding site and the LacI will result in weaker repression [5]. (II) Adding a LacI repressor binding site will increase the chance that LacI binds to the promoter region. </figcaption>
                         </figure>
                         </figure>
                         <p></p>
                         <p></p>
-
                             </section>
+
                              
<br/>
<br/>
                
                
Line 114: Line 117:
<h3>Statistical Mechanics</h3>
<h3>Statistical Mechanics</h3>
<p>
<p>
-
Statistical mechanics provides a way to understand how the macroscopic properties of a system arise from the individual components within this system. In effect there are two states to consider in statistical mechanic descriptions. A macrostate specifies a system in terms of quantities that average over the microscopic constituents of the system. For example if you consider a closed container of gas (figure. 2), quantities for the macrostate would be pressure or volume. This means that these macrostate quantities only work if you consider a system with a very large number of particles. Talking about a macrostate of a single molecule makes no sense since the terms are inherently descriptive of large systems. A microstate, however, does describe the properties of a single molecule e.g. kinetic energy and position. The macro- and micro-states come together in the key concept of statistical mechanics. A lot of different microstates can result in the same macrostate, but characterizing the macrostate imposes restraints on the possible microstates[5].  
+
Statistical mechanics provides a way to understand how the macroscopic properties of a system arise from the individual components within this system. In effect there are two states to consider in statistical mechanic descriptions. A macrostate specifies a system in terms of quantities that average over the microscopic constituents of the system. For example if you consider a closed container of gas (figure 3), quantities for the macrostate would be pressure or volume. This means that these macrostate quantities only work if you consider a system with a very large number of particles. Talking about a macrostate of a single molecule makes no sense since the terms are inherently descriptive of large systems. A microstate, however, does describe the properties of a single molecule e.g. kinetic energy and position. The macro- and micro-states come together in the key concept of statistical mechanics. A lot of different microstates can result in the same macrostate, but characterizing the macrostate imposes restraints on the possible microstates [7].  
</p>
</p>
                         <figure>
                         <figure>
                             <img src="https://static.igem.org/mediawiki/2014/6/68/Wageningen_UR_statmech.png" width="90%">
                             <img src="https://static.igem.org/mediawiki/2014/6/68/Wageningen_UR_statmech.png" width="90%">
-
                             <figcaption><b>Figure 2.</b> The image represents the macro and microstates applied in statistical mechanics. The box can be considered a macrostate e.g. volume, temperature, pressure. The molecules with different microstates e.g. kinetic energy, direction. The molecules determine the pressure and temperature but a given temperature provides information on the kinetic energy a molecule is most likely to have.</figcaption>
+
                             <figcaption><b>Figure 3.</b> The image represents the macro and microstates applied in statistical mechanics. The box can be considered a macrostate e.g. volume, temperature, pressure. The molecules with different microstates e.g. kinetic energy, direction. The molecules determine the pressure and temperature but a given temperature provides information on the kinetic energy a molecule is most likely to have.</figcaption>
                         </figure>
                         </figure>
                         <p></p>
                         <p></p>
 +
<br/>
<h3>Modeling Biological Systems using Statistical Mechanics </h3>
<h3>Modeling Biological Systems using Statistical Mechanics </h3>
 +
<p>
<p>
The logic of macrostates and microstates can also be applied to cells, replacing the container of the previous example with the cell and the different particle constituents with RNA polymerases and repressor proteins. The condition and production states of the various proteins can be considered a macrostate.
The logic of macrostates and microstates can also be applied to cells, replacing the container of the previous example with the cell and the different particle constituents with RNA polymerases and repressor proteins. The condition and production states of the various proteins can be considered a macrostate.
-
Modeling with statistical mechanics provides an insight in transcription regulation of a genetic circuit by generating ‘parameter-free’ predictions for the levels of gene expression. Assuming transcription of a gene to be proportional to chance that a RNA polymerase binds to the promoter of interest a statistical mechanics model has two microstate outcomes [1]:
+
Modeling with statistical mechanics provides an insight in transcription regulation of a genetic circuit by generating ‘parameter-free’ predictions for the levels of gene expression. Assuming transcription of a gene to be proportional to the chance that a RNA polymerase binds to the promoter of interest a statistical mechanics model has two microstate outcomes [1,2]:
</p>
</p>
<p>
<p>
-
(I) A state where none of the polymerases are bound to the gene of interest and distributed among non-specific binding sites (non-specific binding).  
+
<b>Non-specific binding:</b> A state where none of the polymerases are bound to the gene of interest and distributed among non-specific binding sites .  
</p>
</p>
<p>
<p>
-
(II) A state where the promoter of the gene of interest is occupied and the remaining polymerases are distributed among the non-specific binding sites (specific binding).  
+
<b>Specific binding:</b> A state where the promoter of the gene of interest is occupied and the remaining polymerases are distributed among the non-specific binding sites.  
-
In order to determine the probability that either state will be realized we need to know the total number of possible micro-sates and the chance that a particular microstate will be achieved. For the non-specific binding by RNA polymerases combinatorics can be used[1].
+
<br/>
 +
In order to determine the probability that either state will be realized we need to know the total number of possible microstates and the chance that a particular microstate will be achieved. For the non-specific binding by RNA polymerases combinatorics can be used [1]. The following equations will describe step for step how statistical mechanics can model biology, starting with combinatorics.
</p>
</p>
</p>  
</p>  
Line 139: Line 145:
</P>
</P>
<P>
<P>
-
Where P is the number of RNA polymerases and N<sub>ns</sub> the number of non-specific binding sites. Because there is a difference in the binding strength between specific and non-specific binding sites a Boltzmann weight needs to be included[1]:
+
Where P is the number of RNA polymerases and N<sub>ns</sub> the number of non-specific binding sites. Because there is a difference in the binding strength between specific and non-specific binding sites a Boltzmann weight needs to be included [1]:
</P>
</P>
<P>
<P>
Line 145: Line 151:
</P>
</P>
<P>
<P>
-
The Boltzmann weight describes difference between strength of specific and non-specific binding. For the non-specific binding of RNA polymerases the Boltzmann weight ε<sub>pd</sub><sup>NS</sup> is going to be very low compared to specific binding compensating for the high number of binding site arrangements. In order to get the total statistical weight the microstate for RNA polymerase bound to the promoter needs to be added[1].
+
To describes difference between strength of specific and non-specific binding a Boltzmann weight is added. For the non-specific binding of RNA polymerases the Boltzmann weight ε<sub>pd</sub><sup>NS</sup> is going to be very low compared to specific binding, compensating for the high number of binding site arrangements. In order to get the total statistical weight the microstate for RNA polymerase bound to the promoter needs to be added [1,2].
</P>
</P>
<P>
<P>
Line 152: Line 158:
</P>
</P>
<P>
<P>
-
ε<sub>pd</sub><sup>S</sup> is the specific binding energy of the RNA polymerase binding to the promoter. To obtain the chance an RNA polymerase is bound we divide the specific binding chance by the total chance[1].
+
ε<sub>pd</sub><sup>S</sup> is the specific binding energy of the RNA polymerase binding to the promoter. To obtain the chance an RNA polymerase is bound we divide the specific binding chance by the total chance [1,2].
</P>
</P>
<P>
<P>
Line 158: Line 164:
</P>
</P>
<P>
<P>
-
By dividing top and bottom by Z(P<sup>-1</sup>)e<sup>(-(ε<sub>pd</sub><sup>S</sup>)/(K<sub>b</sub>T))</sup> the equation takes the form of a Hill function where the assumption was made that the amount of polymerases is a lot smaller than the amount of non-specific binding sites [1].
+
By dividing top and bottom by Z(P<sup>-1</sup>)e<sup>(-(ε<sub>pd</sub><sup>S</sup>)/(K<sub>b</sub>T))</sup> the equation takes the form of a Hill function where the assumption was made that the amount of polymerases is a lot smaller than the amount of non-specific binding sites [1,2].
</P>
</P>
<P>
<P>
Line 165: Line 171:
</P>
</P>
<P>
<P>
-
In this form a fold change of mRNA expression can be correlated to the chance that a polymerase binds. If the situation is thus N<sub>ns</sub>/P e<sup>(-(∆ε<sub>pd</sub><sup>S</sup>)/(K<sub>b</sub>T))</sup> = 0, that the chance of an RNA polymerase binding is equal one, there is no change in expression (i.e. the fold change in expression remains one). This will most likely be the case in a system without interference. Most biological systems however are dependent on regulators, repression is the main method of gene regulation in our designed circuit (fig.1). The effect of repressors can be taken up in previous equation by adding the microstate ‘repressor protein bound’ as a regulatory region to the total statistical weight. The term ‘repressor binding’ states that if a repressor is bound a polymerase cannot bind. The equation assumes a high basal level of expression in the absence of repressors. Note that as equation (1.5) takes the form of a Hill function, the expression fold change is dependent on the number of repressors [1].
+
In this form a fold change of mRNA expression can be correlated to the chance that a polymerase binds (N<sub>ns</sub>/P e<sup>(-(∆ε<sub>pd</sub><sup>S</sup>)/(K<sub>b</sub>T))</sup>). If the term describing RNA polymerase binding equals zero, the the chance that a RNA polymerase binds is equal to one and there is no change in gene expression. This will most likely be the case in a system without interference. Most biological systems, however, are dependent on regulators, repression is the main method of gene regulation in our designed circuit (see <a href="https://2014.igem.org/Team:Wageningen_UR/project/kill-switch" class="soft_link">Kill-switch</a>). The effect of repressors can be taken up in previous equation by adding the microstate ‘repressor protein bound’ as a regulatory region to the total statistical weight (see equation below). The term ‘repressor binding’ states that if a repressor is bound a polymerase cannot bind. The equation assumes a high basal level of expression in the absence of repressors. Note that like the previous equation, it takes the form of a Hill function, the gene expression fold change is dependent on the number of repressors [1,2].
</P>
</P>
<P>
<P>
Line 171: Line 177:
</P>
</P>
<P>
<P>
-
In equation (1.5) R is the number of repressor proteins and ∆ε_d<sup>r</sup> the specific binding energy of the repressor proteins, the same combinatorics rules apply.
+
R is the number of repressor proteins and ∆ε_d<sup>r</sup> the specific binding energy of the repressor proteins, the same combinatorics rules apply repressor binding as for polymerase binding.
-
n is the number of repression regulators should those regulators have the same binding energy. Because they act independently on whether or not an RNA polymerase can bind, the regulatory values are multiplied7.
+
n is the number of repression regulators should those regulators have the same binding energy. Because they act independently on whether or not an RNA polymerase can bind, the regulatory values are multiplied.
-
Equation (1.5) can now be incorporated into the model describing the kill-switch by tying production rates to the fold change of expression (K<sub>prod</sub> x Fold change), i.e. the production rate will decrease proportionally to changes in expression.   
+
The equation can now be incorporated into the model describing the Kill-Switch by tying production rates to the fold change of gene expression (K<sub>prod</sub> x Fold change), i.e. the production rate will decrease proportionally to changes in gene expression.   
</P>
</P>
 +
<br/>
<h3>Model Equations</h3>
<h3>Model Equations</h3>
<P>
<P>
-
Ordinary differential equations describe the genetic circuit, three of which use a statistical mechanics based Hill function. The deg<sub>protein</sub> and Dil<sub>div</sub> terms represent the half-life dilution rate of the proteins. Because the third equation (the CIλ production) relies on an external input and is not dependent on any of the other equations, the ordinary differential equations uses Michaelis-Menten kinetics to simulate a rhamnose input. In every simulation rhamnose is inputted at t = 10 hours [1][10]
+
Ordinary differential equations describe the genetic circuit, three of based on the statistical mechanics approach. The Deg<sub>protein</sub> and Dil<sub>div</sub> terms represent the half-life and dilution rate of the proteins. Because the third equation (the CIλ production) relies on an external input and is not dependent on any of the other equations, the ordinary differential equations uses Michaelis-Menten kinetics to simulate a fusaric acid input. In every simulation fusaric acid is inputted at t = 10 hours [1,2,10].
</P>
</P>
<br\>
<br\>
<P>
<P>
-
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mi mathvariant="bold-italic"> </mi><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">o</mi><msub><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi></mrow></msub></mrow></msub><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mn>1</mn><mo>+</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup></mrow><mrow><mn>1</mn><mo>+</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup><mo>+</mo><mo>(</mo><mfrac><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">V</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi><mo>_</mo><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup><msup><mrow><mo>)</mo></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mo>_</mo><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi></mrow></msup></mrow></mfrac><msub><mrow><mo>-</mo><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic"> </mi></mrow></msub><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic"> </mi><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">D</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">l</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">v</mi></mrow></msub></math>
+
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mi mathvariant="bold-italic"> </mi><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">o</mi><msub><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi></mrow></msub></mrow></msub><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mn>1</mn><mo>+</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup></mrow><mrow><mn>1</mn><mo>+</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup><mo>+</mo><mo>(</mo><mfrac><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">V</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi><mo>_</mo><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup><msup><mrow><mo>)</mo></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mo>_</mo><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi></mrow></msup></mrow></mfrac><msub><mrow><mo>-</mo><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">D</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic"> </mi></mrow></msub><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic"> </mi><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">D</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">l</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">v</mi></mrow></msub><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi></math>
</P>
</P>
<br\>
<br\>
<P>
<P>
-
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mi mathvariant="bold-italic"> </mi><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">o</mi><msub><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow></msub></mrow></msub><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mn>1</mn><mo>+</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup></mrow><mrow><mn>1</mn><mo>+</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup><mi mathvariant="bold-italic"> </mi><mo>+</mo><mi mathvariant="bold-italic"> </mi><mo>(</mo><mfrac><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">V</mi><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup><msup><mrow><mo>)</mo></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi></mrow></msup><mi mathvariant="bold-italic"> </mi><mo>(</mo><mfrac><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">V</mi><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi><mo>_</mo><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup><msup><mrow><mo>)</mo></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi><mo>_</mo><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">r</mi></mrow></msup></mrow></mfrac><msub><mrow><mo>-</mo><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic"> </mi></mrow></msub><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic"> </mi><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">D</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">l</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">v</mi></mrow></msub></math>
+
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mi mathvariant="bold-italic"> </mi><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">o</mi><msub><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi></mrow></msub></mrow></msub><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mn>1</mn><mo>+</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup></mrow><mrow><mn>1</mn><mo>+</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup><mo>+</mo><mo>(</mo><mfrac><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">V</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi><mo>_</mo><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup><msup><mrow><mo>)</mo></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mo>_</mo><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi></mrow></msup></mrow></mfrac><msub><mrow><mo>-</mo><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">D</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic"> </mi></mrow></msub><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic"> </mi><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">D</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">l</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">v</mi></mrow></msub><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi></math>
</P>
</P>
<br\>
<br\>
<P>
<P>
-
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mi mathvariant="bold-italic"> </mi><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">o</mi><msub><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">λ</mi><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi></mrow></msub></mrow></msub><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mo>[</mo><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">H</mi><mi mathvariant="bold-italic">A</mi><mo>]</mo></mrow><mrow><mfenced open="[" close="]" separators="|"><mrow><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">H</mi><mi mathvariant="bold-italic">A</mi></mrow></mfenced><mo>+</mo><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi></mrow></msub></mrow></mfrac><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">λ</mi><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic"> </mi></mrow></msub><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi><mi mathvariant="bold-italic"> </mi><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">D</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">l</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">v</mi></mrow></msub></math>
+
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mi mathvariant="bold-italic"> </mi><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">o</mi><msub><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow></msub></mrow></msub><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mo>[</mo><mi mathvariant="bold-italic">f</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="bold-italic">s</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">d</mi><mo>]</mo></mrow><mrow><mfenced open="[" close="]" separators="|"><mrow><mi mathvariant="bold-italic">f</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="bold-italic">s</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">d</mi></mrow></mfenced><mo>+</mo><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi></mrow></msub></mrow></mfrac><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic"> </mi></mrow></msub><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi><mi mathvariant="bold-italic"> </mi><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">D</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">l</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">v</mi></mrow></msub><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></math>
</P>
</P>
<br\>
<br\>
<P>
<P>
-
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">G</mi><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mi mathvariant="bold-italic"> </mi><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">o</mi><msub><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi></mrow></msub></mrow></msub><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mn>1</mn><mo>+</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup></mrow><mrow><mn>1</mn><mo>+</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup><mo>+</mo><mo>(</mo><mfrac><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">V</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mi mathvariant="bold-italic">G</mi><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">P</mi></mrow></msub></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup><msup><mrow><mo>)</mo></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><msub><mrow><mi mathvariant="bold-italic">t</mi></mrow><mrow><mi mathvariant="bold-italic">G</mi><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">P</mi></mrow></msub></mrow></msup><mo>(</mo><mfrac><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">V</mi><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><msub><mrow><mi mathvariant="bold-italic">λ</mi></mrow><mrow><mi mathvariant="bold-italic">G</mi><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">P</mi></mrow></msub></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup><msup><mrow><mo>)</mo></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi><mo>_</mo><mi mathvariant="bold-italic">G</mi><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">P</mi></mrow></msup></mrow></mfrac><msub><mrow><mo>-</mo><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">G</mi><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">P</mi><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic"> </mi></mrow></msub><mi mathvariant="bold-italic">G</mi><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">P</mi><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">D</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">l</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">v</mi></mrow></msub></math>
+
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mi mathvariant="bold-italic"> </mi><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">o</mi><msub><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi></mrow></msub></mrow></msub><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mn>1</mn><mo>+</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup></mrow><mrow><mn>1</mn><mo>+</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup><mo>+</mo><mo>(</mo><mfrac><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">V</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">n</mi></mrow></msub></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup><msup><mrow><mo>)</mo></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><msub><mrow><mi mathvariant="bold-italic">t</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">n</mi></mrow></msub></mrow></msup><mo>(</mo><mfrac><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">V</mi><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">N</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">s</mi></mrow></msub></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mo>∆</mo><msubsup><mrow><mi mathvariant="bold-italic">ε</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi><mi>t</mi><mi>o</mi><mi>x</mi><mi>i</mi><mi>n</mi></mrow></msubsup></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><mi mathvariant="bold-italic">T</mi></mrow></mfrac></mrow></msup><msup><mrow><mo>)</mo></mrow><mrow><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi><mo>_</mo><mi mathvariant="bold-italic">G</mi><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">P</mi></mrow></msup></mrow></mfrac><msub><mrow><mo>-</mo><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic"> </mi></mrow></msub><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">n</mi><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">D</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">l</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">v</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">n</mi></mrow></msub></math>
-
 
+
</P>
</P>
 +
<h3>Software used<h3>
 +
<p>
 +
All models were made in Python (Python software foundation, Delaware, US) version 2.71 for Mac OS X version 10.5.8. Solving of ODE systems was done using the standard soln function in Python
 +
</p>
-
                         </section>
+
 
 +
                          
<br/>
<br/>
                  
                  
Line 204: Line 215:
                         <h3>System Behaviour and Optimal Repressor Binding Site Configuration</h3>
                         <h3>System Behaviour and Optimal Repressor Binding Site Configuration</h3>
<p>
<p>
-
Before the optimal configuration of repressor binding sites could be determined, the behavior of the system had to be assigned scores representing a functioning or non-functioning toggle switch (see Parameters). The scores range from zero to two, three different behaviors could be identified:
+
Before the optimal configuration of repressor binding sites could be identified, the behaviour of the system had to be assigned scores representing a functioning or non-functioning Kill-switch. Three distinct behaviours could be identified:
</p>
</p>
<br/>
<br/>
<figure>
<figure>
<img src="https://static.igem.org/mediawiki/2014/9/92/Wageningen_UR_modeling_color_blue.png
<img src="https://static.igem.org/mediawiki/2014/9/92/Wageningen_UR_modeling_color_blue.png
-
" width="5%" style="float:left;margin-left:25px; margin-right:15px;"><figcaption style="text-align:left">0: The system does not work, the toggle switch is out of balance and does not function, the system favors either LacI or TetR </figcaption>
+
" width="5%" style="float:left;margin-left:25px; margin-right:15px;"><figcaption style="text-align:left">0: The Kill-switch does not work, the toggle switch is out of balance and does not function. </figcaption>
</figure>
</figure>
<br/>
<br/>
<figure>
<figure>
<img src="https://static.igem.org/mediawiki/2014/c/c8/Wageningen_UR_modeling_color_green.png
<img src="https://static.igem.org/mediawiki/2014/c/c8/Wageningen_UR_modeling_color_green.png
-
" width="5%" style="float:left;margin-left:25px; margin-right:15px;"><figcaption style="text-align:left">1: The system performs less efficiently, though the toggle switch changes state, the GFP promoter is leaky <br/></figcaption>
+
" width="5%" style="float:left;margin-left:25px; margin-right:15px;"><figcaption style="text-align:left">1: The system performs less efficiently, though the toggle switch changes state, the toxin promoter is leaky. <br/></figcaption>
</figure>
</figure>
-
<br/>
 
<br/>
<br/>
<figure>
<figure>
<img src="https://static.igem.org/mediawiki/2014/d/d6/Wageningen_UR_modeling_color_red.png
<img src="https://static.igem.org/mediawiki/2014/d/d6/Wageningen_UR_modeling_color_red.png
-
" width="5%" style="float:left;margin-left:25px; margin-right:15px;"><figcaption style="text-align:left">2: The system performs to design, after a rhamnose input the toggle switch changes state and GFP is produced when CIλ leaves the system</figcaption>
+
" width="5%" style="float:left;margin-left:25px; margin-right:15px;"><figcaption style="text-align:left">2: The system performs to design, after a fusaric acid input the toggle switch changes state and toxin is produced when CIλ leaves the system.</figcaption>
</figure>
</figure>
<br/>
<br/>
 +
<p>
 +
In total eight (A,B,..,H) different configurations of repressor binding sites were tested (Figure 4). Two: F and H, were found to be produce Kill-switch systems robust towards variations in binding strengths and production rates. The result of the best scoring configuration H has been magnified (figure 5). The repressor binding site configuration corresponding to this map has the highest chance of working in an experimental set-up. The x-axis shows the fold change in protein production, the y-axis has increased binding strengths for individual repressors by an order of magnitude each row. The color map indicates that increases in repressor binding strengths can be compensated for by increased production rates of either the target gene or decreased production of the repressor itself. The practical implications for constructing the an optimized version of the Kill-switch were as follows: two TetR repressor binding sites on the LacI promoter, two LacI and one CIλ repressor binding site on the TetR promoter, two TetR and one CIλ repressor binding site on the toxin promoter.
 +
</p>
                         <figure>
                         <figure>
Line 230: Line 243:
                         <figure>
                         <figure>
                             <img src="https://static.igem.org/mediawiki/2014/2/23/Wageningen_UR_model_Bob_Table1.png" width="90%">
                             <img src="https://static.igem.org/mediawiki/2014/2/23/Wageningen_UR_model_Bob_Table1.png" width="90%">
-
                             <figcaption><b>Figure 1.</b> Color maps indicating functioning and non-functioning systems. Each letter represents different repressor binding site configurations. Each small square within the color maps represents a score for a simulation of the system with a unique set of parameters. The colors correspond to a working system (2:red), leaky system (1:green) and a non-functioning system (0:blue). The table underneath shows the configuration of the repressor binding sites. The type of repressor binding site is specified followed by the promoter it is located on. The letters correspond to the colormaps </figcaption>
+
                             <figcaption><b>Figure 4.</b> Color maps indicating functioning and non-functioning systems. Each letter represents different repressor binding site configurations. Each small square within the color maps represents a score for a simulation of the system with a unique set of parameters. The colors correspond to a working system (2:red), leaky system (1:green) and a non-functioning system (0:blue). The table underneath shows the configuration of the repressor binding sites. The type of repressor binding site is specified followed by the promoter it is located on. The letters correspond to the colormaps </figcaption>
                         </figure>
                         </figure>
                         <p></p>
                         <p></p>
<br/>   
<br/>   
-
<P>
+
-
The result of the best scoring configuration H has been magnified (figure 2). The repressor binding site configuration corresponding to this map has the highest chance of working in an experimental set-up. The x-axis shows the production rates in micro-molar per minute, the y-axis has increased binding strengths for individual repressors by an order of magnitude each row. The color map indicates that increases in repressor binding strengths can be compensated for by increased production rates of either the target gene or decreased production of the repressor itself.
+
-
</P>
+
                         <figure>
                         <figure>
                             <img src="https://static.igem.org/mediawiki/2014/d/d8/Wageningen_UR_modelbigcolormap.png" width="75%">
                             <img src="https://static.igem.org/mediawiki/2014/d/d8/Wageningen_UR_modelbigcolormap.png" width="75%">
-
                             <figcaption><b>Figure 1.</b> Color maps indicating functioning and non-functioning systems. Each letter represents different repressor binding site configurations. Each small square within the colour maps represents a score for a simulation of the system with a unique set of parameters. The colours correspond to the previously given description </figcaption>
+
                             <figcaption><b>Figure 5.</b> The best scoring repressor binding site configuration H. Each small square within the colour maps represents a score for a simulation of the system with a unique set of parameters. The x-axis shows the production rates in uM/min. The Y-axis has the repressor binding strengths. Using the lowest value found in literature [11], each row the value for th repressor binding strength increases by an order of magnitude.</figcaption>
                         </figure>
                         </figure>
<br/>
<br/>
<h3>Statistical Mechanics in Practice</h3>
<h3>Statistical Mechanics in Practice</h3>
<P>
<P>
-
Statistical mechanics allows you to describe the system in terms of the number of repressors and repressor binding sites, their binding affinity and the position of the repressor binding on the promoter.
+
There were two main reasons for using statistical mechanics:<br/><br/>
-
In contrast to Michaelis-Menten kinetics, statistical mechanics offers more freedom studying the behavior of differently designed promoters. Should Michaelis-Menten kinetics have been used the dissociation constant of a protein to DNA would need to be experimentally determined for different promoter configurations in order to obtain an indication of system behavior, a ‘parameter dependent’ method of modeling. Statistical mechanics is not bound by experimentation in making a rough predictions. Though the model still needs parameters e.g. number of RNA polymerases in a cell, number of non-specific binding sites, these are accessible without a need for experimentation1.  
+
<b>(I)</b> Statistical mechanics provides a way describe the system in terms of the number of repressors and repressor binding sites, their binding affinity and the position of the repressor binding on the promoter [1,2].<br/><br/>
-
The model predicts the only functioning repressor binding site configurations to be those where LacI and TetR are placed pairwise on the promoter, a result that is confirmed in the literature2.
+
<b>(II)</b> In contrast to Michaelis-Menten kinetics, statistical mechanics offers more freedom studying the behavior of differently designed promoters.Michaelis-Menten requires the dissociation constants of the repressor proteins to be experimentally determined [1,2]. <br/><br/>
-
This means the statistical mechanics approaches can estimate which repressor binding site configuration is most likely to work. The information provided by this model has inspired the design of the promoters for the kill-switch.
+
Based on the results of the model the experimentalists ordered oligonucleotides and assembled six promoters. Building two Kill-switches in parallel the chance of constructing a stable system increases. The assumption that the system is optimized and stable has been incorporated in the model quantifying stability .
-
The next modeling step would be to assess the Cost and Performance of our entire system. This will give insight in the overall capacity and behavior of our biological control agent in the soil.
+
</P>
</P>
-
                            </section>
+
  </section>                        
-
                             </section>
+
                              
<br/>
<br/>
                         <section id="cost1">
                         <section id="cost1">
-
                        <h1 >System Cost</h1>
+
<h1 >System Cost</h1>
<br/>
<br/>
-
                        <h2 id="introduction2">Introduction</h2>
+
<h2 id="introduction2">Introduction</h2>
-
                        <p>
+
<p>
-
                            BananaGuard was designed without regarding the possibility, that our engineered <i>Pseudomonas putida</i> would be out competed by other bacteria and fungi found in the rhizosphere of banana roots[1,2]. However, it is crucial for a successful application that our strain can sustain near the banana roots. This seemed challenging as our system consists of multiple constitutively expressed genes controlling the <a href="https://2014.igem.org/Team:Wageningen_UR/project/kill-switch" class="soft_link">kill switch</a> and <a href="https://2014.igem.org/Team:Wageningen_UR/project/gene_transfer" class="soft_link">toxin-antitoxin</a> systems. Moreover, when the system is switched to the active state by <a href="https://2014.igem.org/Team:Wageningen_UR/project/fungal_sensing#wetlab" class="soft_link">sensing fusaric acid</a>, <i>P. putida</i> produces multiple <a href="https://2014.igem.org/Team:Wageningen_UR/project/fungal_sensing#wetlab" class="soft_link">antifungal enzymes and compounds</a>. The integrated synthetic pathways uses metabolic resources that would otherwise be dedicated to cellular maintenance or growth. Thereby, the synthetic pathways reduce the bacteria’s potential to sustain themselves in the rhizosphere. Because of this reduced potential we wondered to what degree P. putida metabolism was affected by our system and, if need be, how the system could be changed to lessen the impact. To approach this problem a genome-scale metabolic model (GSMM) is used[3]. Our extended GSMM predicts the interaction of our integrated plasmids with genes originating from endogenous pathways. With this prediction of the metabolic stress we get an insight in the metabolic burden caused by the integration of our plasmids and whether the engineered <i>P. putida</i> strain is losing the capability to compete with other rhizosphere-populating microorganisms in the resting state.
+
BananaGuard was designed without knowing whether our engineered <i>Pseudomonas putida</i> would be outcompeted by other bacteria and fungi found in the rhizosphere of banana roots [12, 13]. For a successful application it is crucial that our strain can survive and compete near the banana plant roots. This seemed challenging as our system consists of multiple constitutively expressed genes controlling the <a href="https://2014.igem.org/Team:Wageningen_UR/project/kill-switch" class="soft_link">kill switch</a> and <a href="https://2014.igem.org/Team:Wageningen_UR/project/gene_transfer" class="soft_link">toxin-antitoxin</a> systems. Moreover, when the system is switched to the active state by <a href="https://2014.igem.org/Team:Wageningen_UR/project/fungal_sensing#wetlab" class="soft_link">sensing fusaric acid</a>, BananaGuard produces multiple <a href="https://2014.igem.org/Team:Wageningen_UR/project/fungal_sensing#wetlab" class="soft_link">fungal growth inhibitor enzymes and compounds</a>. The The results display the intrinsic behavior of Kill-switch stability and performance.<br/><br/>integrated synthetic pathways use metabolic resources that would otherwise be dedicated to cellular maintenance or growth. The synthetic pathways therefore reduce the bacteria's potential to sustain themselves in the rhizosphere. Because of this reduced potential we wondered to what degree <i>P. putida</i> metabolism was affected by the introduced synthetic genes and if need be, to what extent the system could be changed to lessen the impact. To approach this problem we used a genome-scale metabolic model (GSMM) [14]. Our GSMM extended by the BananaGuard system predicts the interaction of our integrated plasmids with genes originating from endogenous pathways. With this prediction of the metabolic stress we gain insight in the metabolic burden caused by the integration of our plasmids. This will give an indication about  BananaGuard if it will lose its capability to compete with other rhizosphere-populating micro-organisms in its resting state.
-
                        </P>
+
<br/>
 +
Continue to learn more about a genomic-scale metabolic model or view the <a href="https://2014.igem.org/Team:Wageningen_UR/project/model#results2" class="soft_link">results</a>.
 +
</P>
<br/>
<br/>
-
                        <h2 id="approach2">Theory and experimental design</h2>
+
<h2 id="approach2">Theory and experimental design</h2>
-
                        <p>
+
<p>
-
                            In order to answer our questions with a GSMM we need a model that:
+
In order to answer our questions with a GSMM we need a model that:
-
                        <ul>
+
<ul>
-
                            <li>contains a comprehensive representation of <i>P. putida</i> metabolism</li>
+
<li>contains a comprehensive representation of <i>P. putida</i> metabolism</li>
-
                            <li>is capable of predicting growth rates in various environments</li>
+
<li>is capable of predicting growth rates in various environments</li>
-
                            <li>includes all reactions that are important for the BananaGuard system</li>
+
<li>includes all reactions that are important for the BananaGuard system</li>
-
                        </ul>
+
</ul>
-
                        </p>
+
</p>
-
                        <p>
+
<p>
-
                            For this reason a comparison was made between various GSMM for <i>P. putida</i>.  
+
For this reason a comparison was made between various available GSMMs for <i>P. putida</i>.  
-
                            There are multiple genome-scale metabolic models available for <i>P. Putida</i>. we took a in-depth look in two models that were up-to-date, published and verified: iJP962 and iJN746[4,5].  
+
We took a in-depth look in two models that were up-to-date, published and verified: iJP962 and iJN746[15,16].  
-
                            For the necessary work to model the metabolic stress in <i>P. putida</i>, the best starting model needs to be chosen which has the most potential for accurately modeling the metabolic stress of our system.
+
For the necessary work to model the metabolic stress in <i>P. putida</i>, the best starting model needs to be chosen which has the most potential for accurately modeling the metabolic stress of our system.
-
                        </p>
+
</p>
<figure>
<figure>
-
<figcaption>comparison between iJP962 and iJN746</figcaption>
+
<figcaption>Table 1. Comparison between iJP962 and iJN746</figcaption>
-
                        <ul class="table">
+
<figure>
-
                            <table width="800px">
+
<img src="https://static.igem.org/mediawiki/2014/thumb/d/dd/Wageningen_UR_model_walter_table.png/800px-Wageningen_UR_model_walter_table.png" width="95%">
-
                                <tr>
+
 
-
                                    <th>Genomic Model</th>
+
-
                                    <th>iJP962</th>
+
-
                                    <th>iJN746</th>
+
-
                                </tr>
+
-
                                <tr>
+
-
                                    <td> Reactions </td>
+
-
                                    <td>1070</td>
+
-
                                    <td>957</td>
+
-
                                </tr>
+
-
                                <tr>
+
-
                                    <td> Metabolites </td>
+
-
                                    <td>959</td>
+
-
                                    <td>685</td>
+
-
                                </tr>
+
-
                                <tr>
+
-
                                    <td> Compartments </td>
+
-
                                    <td>1</td>
+
-
                                    <td>2</td>
+
-
                                <tr>
+
-
                                    <td> Focus of the model (based on metabolites and reactions)</td>
+
-
                                    <td> This model has a broad spectrum of reactions in every class, but is lacking the in-depth reactions.</td>
+
-
                                    <td>This model has a focus on the production of interesting compounds for use in biotechnology (aromatic compounds, alkenes, alkanoates). </td>
+
-
                                </tr>
+
-
                            </table>
+
-
                        </ul>
+
</figcaption>
</figcaption>
-
                        <p>
+
<p>
-
                            Taking the above arguments in consideration, model iJP962 is most suitable for modeling our system. For our system in BananaGuard it is not more advantageous to have multiple compartments. Many of the reactions are over the internal membrane and not functional for the core intermediate metabolites. Also iJN746 has a focus on compounds that are in a different class of metabolites then the toxins we chose for our system[5]. Therefore iJP962 is more representative for BananaGuard due to more reactions and more metabolites in a broader spectrum. The primary category of concern for anti-fungal production is the amino acid metabolism and iJP962 has the most reactions in this category[3,4,6]. neither model included all the specific reactions, metabolites and genes that are necessary for modeling the anti-fungal metabolites. These reactions need to be added manually to model iJP962.
+
According to the previously stated requirements,the model iJP962 is most suitable for modeling our system. The biggest difference between these models shown in table 1 is the number of compartments.  But this is not necessarily an advantage. Many of the reactions are over the internal membrane and not functional for the core intermediate metabolites. Added to this, iJN746 has a focus on compounds that are in a different class of metabolites than the fungal growth inhibitors we chose for our system [27]. Therefore iJP962 is more representative due to more reactions and more metabolites in a broader spectrum. The primary category of concern for fungal growth inhibitor production is the amino acid metabolism and iJP962 has the most reactions in this category [14,15,17]. neither model included all the specific reactions, metabolites and genes that are necessary for modeling the fungal growth inhibitor metabolites. These reactions need to be added manually to the GSMM.
-
                        </p>
+
</p>
<br/>
<br/>
-
                        <h3>Extending the model </h3>
+
<h3>Extending the model </h3>
-
                        <p>
+
<p>
-
                            Traditionally, metabolic models only predict metabolic fluxes. Although this is directly applicable to the metabolites that we want to produces, such as DMDS, this does not suffice to answer our questions regarding the impact of the whole system. Therefore, we added the production of all the proteins of interest as if they are metabolites. For this we used the protein sequences to determine the correct stoichiometric ratio of amino acids[7]. The other side-substrates/products are in the same ratio as described for production of proteins for the biomass reaction in the GSMM. In the model <a href="https://2014.igem.org/Team:Wageningen_UR/project/model" class="soft_link">Promoter Design</a> the protein production rate is estimated.
+
Traditionally, metabolic models only predict metabolic fluxes. Although this is directly applicable to the metabolites that we want to produce, such as DMDS, this does not suffice to answer our questions regarding the impact of the whole system. Therefore, we added the production of all the proteins of interest as if they are metabolites. For this we used the protein sequences to determine the correct stoichiometric ratio of amino acids [18]. The other side-substrates/products are in the same ratio as described for production of proteins for the biomass reaction in the GSMM. In the model <a href="https://2014.igem.org/Team:Wageningen_UR/project/model" class="soft_link">Promoter Design</a> the protein production rate is estimated.
-
                        </p>
+
</p>
-
                        <p>
+
<p>
-
                            As our metabolites of interest (DAPG, DMDS, DMTS) were not yet included in the GSMM, we had to deduce their production pathways based on databases and literature information[8,9,10,11](link notebook with all reactions). Furthermore, we took special care that all added pathways are stoichiometrically balanced, which is essential for proper and reliable integration in a GSMM. To get the most realistic results we need to set a biological relevant production rate. This can be done by estimating a production range with in <i>vivo</i> data from literature.
+
As our metabolites of interest (DAPG, DMDS, DMTS) were not yet included in the GSMM, we had to deduce their production pathways based on databases and literature information [19,20,21,22](<a href="https://static.igem.org/mediawiki/2014/7/7d/Wageningen_UR_protocols_Metabolicmodelingpathways.pdf" target="blank" class="soft_link">pathways and metabolites</a>). Furthermore, we took special care that all added pathways are stoichiometrically balanced, which is essential for proper and reliable integration in a GSMM. To get the most realistic results we need to set a biologically relevant production rate. This can be done by estimating a production range with <i>in vivo</i> data from literature.
-
                        </p>
+
</p>
<br/>
<br/>
-
                        <h3> Programs and toolboxes</h3>
+
<h3> Programs and toolboxes</h3>
-
                        <p>
+
<p>
-
                            Matlab[12] was used in combination with the Cobra toolbox[13] for handling the model. The Cobra toolbox is designed for working with genome-scale metabolic models. To solve the flux balance analysis a linear equation solver is needed. Gurobi[14] was chosen because of its compatibility with Matlab and the Cobra toolbox.
+
Matlab[23] was used in combination with the Cobra toolbox[24] for handling the model. The Cobra toolbox is designed for working with genome-scale metabolic models. To solve the flux balance analysis a linear program solver is needed. Gurobi[25] was chosen because of its compatibility with Matlab and the Cobra toolbox.
-
                        </p>
+
</p>
<br/>
<br/>
-
                        <h3> Flux Balance Analysis </h3>
+
<h3> Flux Balance Analysis </h3>
-
                        <p>
+
<p>
-
                            Flux Balance Analysis (FBA) is a framework for investigating the theoretical capabilities of a GSMM and can be solved using a linear programming solver. This is possible because of linear objective function and linear constraints[4,6].   
+
Flux Balance Analysis (FBA) is a framework for investigating the theoretical capabilities of a GSMM and can be solved using a linear programming solver. This is possible because of the linear objective function and linear constraints[4,6].   
-
                            FBA is based on the following assumptions to ensure its linearity:
+
FBA is based on the following assumptions:
-
                        </p>
+
</p>
-
                        <ol>
+
<ol>
-
                            <li>Steady state: The fluxes are considered to have attained a static equilibrium value and do not change through time, this means that accumulation or dissipation is not possible. </li>
+
<li>Steady state: The fluxes are considered to have attained a static equilibrium value and do not change through time, this means that accumulation or dissipation is not possible. </li>
-
                            <li>Thermodynamic constraints: (ir)reversibility of reactions.</li>
+
<li>Thermodynamic constraints: (ir)reversibility of reactions.</li>
-
                            <li>No enzyme saturation: The maximum flux for every reaction is not dependent on the concentration of the enzymes used for that reaction. </li>
+
<li>No enzyme saturation: The maximum flux for every reaction is not dependent on the concentration of the enzymes used for that reaction. </li>
-
                        </ol>
+
</ol>
<br/>
<br/>
-
                        <p>
+
<p>
-
                            The FBA method uses a representation of the metabolic reaction network in the form of a stoichiometric matrix (S) where :  
+
The FBA method uses a representation of the metabolic reaction network in the form of a stoichiometric matrix (S) where :  
-
                        <ul>
+
<ul>
-
                            <li>Each column corresponds to a reaction R<sub>i</sub> </li>
+
<li>Each column corresponds to a reaction R<sub>i</sub> </li>
-
                            <li>Each row corresponds to a metabolite C<sub>j</sub> </li>
+
<li>Each row corresponds to a metabolite C<sub>j</sub> </li>
-
                        </ul>
+
</ul>
-
                        <p>
+
<p>
-
                            The definition of S is :  
+
The definition of S is :  
-
                        </p>
+
</p>
-
                      <math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math" display="block"><msub><mrow><mi>S</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mi>k</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mi>&nbsp;</mi><mi>&nbsp;</mi><mfenced open="{" close="" separators="|"><mrow><mtable><mtr><mtd><mrow><maligngroup/><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>k</mi><mo>&lt;</mo><mi>&nbsp;</mi><mi>&nbsp;</mi><mn>0</mn><mi>&nbsp;</mi><mi>&nbsp;</mi><mo>→</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>j</mi></mrow></msub><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>i</mi><mi>s</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>a</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>s</mi><mi>u</mi><mi>b</mi><mi>s</mi><mi>t</mi><mi>r</mi><mi>a</mi><mi>t</mi><mi>e</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>o</mi><mi>f</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>r</mi><mi>e</mi><mi>a</mi><mi>c</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><msub><mrow><mi>R</mi></mrow><mrow><mi>i</mi></mrow></msub><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>w</mi><mi>i</mi><mi>t</mi><mi>h</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>s</mi><mi>t</mi><mi>o</mi><mi>i</mi><mi>c</mi><mi>h</mi><mi>i</mi><mi>o</mi><mi>m</mi><mi>e</mi><mi>t</mi><mi>r</mi><mi>y</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>k</mi></mrow></mtd></mtr><mtr><mtd><mrow><maligngroup/><mi>k</mi><mo>&gt;</mo><mi>&nbsp;</mi><mi>&nbsp;</mi><mn>0</mn><mi>&nbsp;</mi><mi>&nbsp;</mi><mo>→</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>j</mi></mrow></msub><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>i</mi><mi>s</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>a</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>p</mi><mi>r</mi><mi>o</mi><mi>d</mi><mi>u</mi><mi>c</mi><mi>t</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>o</mi><mi>f</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>r</mi><mi>e</mi><mi>a</mi><mi>c</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><msub><mrow><mi>R</mi></mrow><mrow><mi>i</mi></mrow></msub><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>w</mi><mi>i</mi><mi>t</mi><mi>h</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>s</mi><mi>t</mi><mi>o</mi><mi>i</mi><mi>c</mi><mi>h</mi><mi>i</mi><mi>o</mi><mi>m</mi><mi>e</mi><mi>t</mi><mi>r</mi><mi>y</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>k</mi></mrow></mtd></mtr><mtr><mtd><mrow><maligngroup/><mi>k</mi><mo>=</mo><mn>0</mn><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mo>→</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>j</mi></mrow></msub><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>d</mi><mi>o</mi><mi>e</mi><mi>s</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>n</mi><mi>o</mi><mi>t</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>p</mi><mi>a</mi><mi>r</mi><mi>t</mi><mi>i</mi><mi>c</mi><mi>i</mi><mi>p</mi><mi>a</mi><mi>t</mi><mi>e</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>i</mi><mi>n</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>r</mi><mi>e</mi><mi>a</mi><mi>c</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><msub><mrow><mi>R</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mtd></mtr></mtable></mrow></mfenced></math>
+
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math" display="block"><msub><mrow><mi>S</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mi>k</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mi>&nbsp;</mi><mi>&nbsp;</mi><mfenced open="{" close="" separators="|"><mrow><mtable><mtr><mtd><mrow><maligngroup/><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>k</mi><mo>&lt;</mo><mi>&nbsp;</mi><mi>&nbsp;</mi><mn>0</mn><mi>&nbsp;</mi><mi>&nbsp;</mi><mo>→</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>j</mi></mrow></msub><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>i</mi><mi>s</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>a</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>s</mi><mi>u</mi><mi>b</mi><mi>s</mi><mi>t</mi><mi>r</mi><mi>a</mi><mi>t</mi><mi>e</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>o</mi><mi>f</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>r</mi><mi>e</mi><mi>a</mi><mi>c</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><msub><mrow><mi>R</mi></mrow><mrow><mi>i</mi></mrow></msub><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>w</mi><mi>i</mi><mi>t</mi><mi>h</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>s</mi><mi>t</mi><mi>o</mi><mi>i</mi><mi>c</mi><mi>h</mi><mi>i</mi><mi>o</mi><mi>m</mi><mi>e</mi><mi>t</mi><mi>r</mi><mi>y</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>k</mi></mrow></mtd></mtr><mtr><mtd><mrow><maligngroup/><mi>k</mi><mo>&gt;</mo><mi>&nbsp;</mi><mi>&nbsp;</mi><mn>0</mn><mi>&nbsp;</mi><mi>&nbsp;</mi><mo>→</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>j</mi></mrow></msub><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>i</mi><mi>s</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>a</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>p</mi><mi>r</mi><mi>o</mi><mi>d</mi><mi>u</mi><mi>c</mi><mi>t</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>o</mi><mi>f</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>r</mi><mi>e</mi><mi>a</mi><mi>c</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><msub><mrow><mi>R</mi></mrow><mrow><mi>i</mi></mrow></msub><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>w</mi><mi>i</mi><mi>t</mi><mi>h</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>s</mi><mi>t</mi><mi>o</mi><mi>i</mi><mi>c</mi><mi>h</mi><mi>i</mi><mi>o</mi><mi>m</mi><mi>e</mi><mi>t</mi><mi>r</mi><mi>y</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>k</mi></mrow></mtd></mtr><mtr><mtd><mrow><maligngroup/><mi>k</mi><mo>=</mo><mn>0</mn><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mo>→</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>j</mi></mrow></msub><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>d</mi><mi>o</mi><mi>e</mi><mi>s</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>n</mi><mi>o</mi><mi>t</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>p</mi><mi>a</mi><mi>r</mi><mi>t</mi><mi>i</mi><mi>c</mi><mi>i</mi><mi>p</mi><mi>a</mi><mi>t</mi><mi>e</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>i</mi><mi>n</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><mi>r</mi><mi>e</mi><mi>a</mi><mi>c</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi><mi>&nbsp;</mi><mi>&nbsp;</mi><msub><mrow><mi>R</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mtd></mtr></mtable></mrow></mfenced></math>
-
                        <p>
+
<p>
-
                            The FBA problem is then formulated as a maximization or minimization problem under the determined constraints:  
+
The FBA problem is then formulated as a maximization or minimization problem under the determined constraints:  
-
                        </p>
+
</p>
-
                        <math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math" display="block"><mfenced open="{" close="" separators="|"><mrow><mtable><mtr><mtd><mrow><maligngroup/><mi>M</mi><mi>a</mi><mi>x</mi><mi>i</mi><mi>m</mi><mi>i</mi><mi>z</mi><mi>e</mi><mi>&nbsp;</mi><msup><mrow><mi mathvariant="bold-italic">c</mi></mrow><mrow><mi>T</mi></mrow></msup><mi>v</mi></mrow></mtd></mtr><mtr><mtd><mrow><maligngroup/><mi>S</mi><mi>u</mi><mi>b</mi><mi>j</mi><mi>e</mi><mi>c</mi><mi>t</mi><mi>&nbsp;</mi><mi>t</mi><mi>o</mi><mi>&nbsp;</mi><mi>S</mi><mi>v</mi><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><maligngroup/><mi>L</mi><mi>o</mi><mi>w</mi><mi>e</mi><mi>r</mi><mi>b</mi><mi>o</mi><mi>u</mi><mi>n</mi><mi>d</mi><mo>≤</mo><mi>v</mi><mo>≤</mo><mi>U</mi><mi>p</mi><mi>p</mi><mi>e</mi><mi>r</mi><mi>b</mi><mi>o</mi><mi>u</mi><mi>n</mi><mi>d</mi></mrow></mtd></mtr></mtable></mrow></mfenced></math>
+
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math" display="block"><mfenced open="{" close="" separators="|"><mrow><mtable><mtr><mtd><mrow><maligngroup/><mi>M</mi><mi>a</mi><mi>x</mi><mi>i</mi><mi>m</mi><mi>i</mi><mi>z</mi><mi>e</mi><mi>&nbsp;</mi><msup><mrow><mi mathvariant="bold-italic">c</mi></mrow><mrow><mi>T</mi></mrow></msup><mi>v</mi></mrow></mtd></mtr><mtr><mtd><mrow><maligngroup/><mi>S</mi><mi>u</mi><mi>b</mi><mi>j</mi><mi>e</mi><mi>c</mi><mi>t</mi><mi>&nbsp;</mi><mi>t</mi><mi>o</mi><mi>&nbsp;</mi><mi>S</mi><mi>v</mi><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><maligngroup/><mi>L</mi><mi>o</mi><mi>w</mi><mi>e</mi><mi>r</mi><mi>b</mi><mi>o</mi><mi>u</mi><mi>n</mi><mi>d</mi><mo>≤</mo><mi>v</mi><mo>≤</mo><mi>U</mi><mi>p</mi><mi>p</mi><mi>e</mi><mi>r</mi><mi>b</mi><mi>o</mi><mi>u</mi><mi>n</mi><mi>d</mi></mrow></mtd></mtr></mtable></mrow></mfenced></math>
-
                        <p>
+
<p>
-
                            Where :  
+
Where :  
-
                        <ul>
+
<ul>
-
                            <li><b>v</b> is the vector of to be determined reaction fluxes </li>
+
<li><b>v</b> is the vector of to be determined reaction fluxes </li>
-
                            <li><b>c</b> is a vector of constants defining the objective function </li>
+
<li><b>c</b> is a vector of constants defining the objective function </li>
-
                            <li>S is the stoichiometry matrix </li>
+
<li>S is the stoichiometric matrix </li>
-
                            <li><b>Lowerbound</b> and <b>Upperbound</b> are vectors of constraints indicating, minimal and maximal flux values for each reaction</li>
+
<li><b>Lowerbound</b> and <b>Upperbound</b> are vectors of constraints indicating, minimal and maximal flux values for each reaction</li>
-
                        </ul>
+
</ul>
-
                        </p>
+
</p>
<br/>
<br/>
-
                        <h3> Media composition</h3>
+
<h3> Media composition</h3>
-
                        <p>
+
<p>
-
Within a GSMM the media composition is important for the results. By changing the metabolites and their uptake rates in the media you can dictate what’s available in the model. For most of the nutrients available in the media you need arbitrary bounds as actual uptake rates are unknown. This also enforces the later decision to not include all root components but only the most prevalent. You don’t want that rare compounds have a huge impact on the model predictions. This does imply that the model predictions could be slightly negative: with including the other available compounds the system may work even better.
+
GSMMs consider all metabolic possibilities given the available substrates.  
-
To get the most realistic approximation for BananaGuard in the soil the M9 media composition was used with carbon as limiting nutrient. The carbon source is based on the exudates of banana roots and consists of 34% sucrose and 66% divided over multiple amino acids[15]. We have neglected other secondary metabolites in banana roots exudates because they are exudated in relatively small amounts[16]. 
+
By changing the substrates and their uptake rates in the media you can dictate what’s available in the model. However, for most metabolites the uptake rates are unknown. This can be solved setting these fluxes as unconstrained in the GSMM. This also enforces decision to not include all root components but only the most prevalent. You don’t want that rare compounds have a huge impact on the model predictions. This does imply that the model predictions could be slightly negative: with including the other available compounds the system may work even better.
-
                        </p>
+
</p>
-
                        </p>
+
<p>
 +
To get the most realistic approximation for BananaGuard in the soil the M9 media composition was used with carbon as limiting nutrient. The carbon source is based on the exudates of banana roots and consists of 34% sucrose and 66% divided over multiple amino acids[15, 27].  
 +
</p>
 +
</p>
<br/>
<br/>
-
                      <h2 id="results2">Results</h2>
+
<h2 id="results2">Results</h2>
<p>
<p>
-
<i>P. putida</i> can be found normally in the soil populating the rhizosphere of the roots.
+
<i>P. putida</i> naturally habits the rhizosphere of plant roots.
To check if BananaGuard is still viable in its resting and active state, a comparison is important between the wild type <i>P. putida</i> and BananaGuard. If the growth rate of BananaGuard approximates the growth rate of wild type <i>P. putida</i> in the resting state it is safe to assume that the disadvantage of the introduced synthetic pathway does not have a huge impact on competitiveness with other rhizosphere-populating microorganisms. This means that BananaGuard should be able to populate the rhizosphere of the banana roots.   
To check if BananaGuard is still viable in its resting and active state, a comparison is important between the wild type <i>P. putida</i> and BananaGuard. If the growth rate of BananaGuard approximates the growth rate of wild type <i>P. putida</i> in the resting state it is safe to assume that the disadvantage of the introduced synthetic pathway does not have a huge impact on competitiveness with other rhizosphere-populating microorganisms. This means that BananaGuard should be able to populate the rhizosphere of the banana roots.   
-
For the active state it is not essential anymore to compete with other microorganisms he capability to produce all the anti-fungals with a constant rate around the roots, without overloading a metabolic pathway and causing shortage of intermediates.  
+
For the active state it is not essential to compete with other micro-organisms but the capability to produce all the fungal growth inhibitors with a constant rate around the roots, without overloading a metabolic pathway and causing shortage of intermediates.  
</p>
</p>
<p>
<p>
-
We aim for the best approximation of the environment of the banana roots for this we need the carbon uptake rate. Because these units are not a standard in literature we have to calculate them and for this we need the doubling time of <i>P. putida</i> in the rhizosphere. This doubling time is 3 hours resulting in a growth rate of 0.23 h<sup>-1</sup>[1]. From this we can calculate the carbon uptake rate in the rhizosphere which is around 4.5 mmol gDW<sup>-1</sup>hr<sup>-1</sup>.
+
We aim for the best approximation of the environment of the banana roots. To do this we need to estimate the carbon uptake rate of <i>P. Putida</i>. Because these units are not a standard in literature we have to calculate them. Within this GSMM it is possible to calculate the growth rate with known nutrient uptake rates. This can be reversed, so if we know the the doubling time of <i>P. putida</i> in the rhizosphere the carbon uptake rate can be calculated. This doubling time is 3 hours resulting in a growth rate of 0.23 h<sup>-1</sup> [1]. From this we can calculate the carbon uptake rate in the rhizosphere which is around 4 mmol gDW<sup>-1</sup>hr<sup>-1</sup>.
-
Figure 1 depicts how the growth rate of BananaGuard compares to the growth rate of the wild type <i>P. putida</i> depending on the carbon uptake rate of the organism in the resting state and active state of the kill switch. Glucose was chosen as carbon source reference, because of the efficient degradation of glucose.
+
-
BananaGuard is still able to grow at almost max growth rate (>99%) compared to the wild type in the resting state. We can conclude from this that BananaGuard is not out competed by other rhizosphere-populating microorganisms because of the synthetic pathway.
+
-
When the kill switch is activated by fusaric acid BananaGuard is still able to growth with 50-70% of the maximum growth of the wild type (Figure 1). This indicates that metabolic stress is not a bottleneck for the production of anti-fungals in our fusaric acid activated system.
+
</p>
</p>
<p>
<p>
-
Next we want to know if our conclusions made for glucose is also valid for the environment close to the banana roots. To approach this as realistic possbible the carbon source is changed to have the same ratio of carbohydrates and amino acids as found in the exudates of banana roots[15]. This new carbon composition causes a reduced growth rate because of less efficient degradation of these compounds.
+
Figure 5 depicts how the growth rate of BananaGuard compares to the growth rate of the wild type <i>P. putida</i> depending on the carbon uptake rate of the organism in the resting state and active state of the kill switch. Glucose was chosen as carbon source reference, because of the efficient degradation of glucose.
-
In figure 1 the relative growth rate is shown for the resting and active state based on root exudate composition. The relative growth rate of BananaGuard around the roots has slightly decreased compared to glucose as carbon source but is still above 45% indicating it still has the capacity to grow and maintain the constant anti-fungals production with one condition that BananaGuard needs to stay close to the roots to maintain its carbon uptake rate. From this we can conclude that these changes in the carbon composition do not change our conclusions previously made.
+
BananaGuard is still able to grow at almost max growth rate (>99%) compared to the wild type in the resting state. We can conclude from this that BananaGuard is not outcompeted by other rhizosphere-populating microorganisms because of the synthetic pathway.
 +
When the kill switch is activated by fusaric acid, BananaGuard is still able to growth with 50-65% of the maximum growth of the wild type (Figure 5). This indicates that metabolic stress is not a bottleneck for the production of fungal growth inhibitors in our fusaric acid activated system.
 +
</p>
 +
<p>
 +
Next we want to know if our conclusions made for glucose are also valid for the environment close to the banana roots. To approach this realistic as possible the carbon source is changed to have the same ratio of carbohydrates and amino acids as found in the exudates of banana roots [26]. This new carbon composition causes a reduced growth rate because of less efficient degradation of these compounds.
 +
In figure 5 the relative growth rate is shown for the resting and active state based on root exudate composition. The relative growth rate of BananaGuard around the roots has slightly decreased compared to glucose as carbon source but is still above 45%. This indicates  it still has the capacity to grow and maintain the constant fungal growth inhibitor production with one condition that BananaGuard needs to stay close to the roots to maintain its carbon uptake rate. From this we can conclude that these changes in the carbon composition do not change our previously made conclusions .
</p>
</p>
<figure>
<figure>
<img src="https://static.igem.org/mediawiki/2014/4/4a/Wageningen_UR_modeling_Combined_growth_rates.png" width="80%">
<img src="https://static.igem.org/mediawiki/2014/4/4a/Wageningen_UR_modeling_Combined_growth_rates.png" width="80%">
-
<figcaption>Figure 2: The relative growth rate compared to the wild type <i>P. putida</i> for different carbon uptake rates. The realistic solution is with the banana exudates as carbon source and the other sulution is with glucose as reference. The expected carbon uptake rate of <i>P. putida</i> in the rhizosphere is indicated with transparent red.</figcaption>
+
<figcaption>Figure 5. The relative growth rate compared to the wild type <i>P. putida</i> for different carbon uptake rates. Two different carbon sources banana exudates and and glucose as reference. The expected carbon uptake rate of <i>P. putida</i> in the rhizosphere is indicated with transparent red.</figcaption>
</figure>
</figure>
<br/>
<br/>
-
                        <h3> Oxygen necessity in the soil </h3>
+
<h3> Oxygen necessity in the soil </h3>
-
                        <p>
+
<p>
-
                            <i>P. putida</i> is a strict aerobic bacteria[2]. This means that without oxygen available the bacteria cannot grow. To inhibit the growth of the <i>Fusarium oxysporum</i> BananaGuard must grow in the rhizosphere of the banana roots. If there is not enough oxygen present in the rhizosphere it might not be possible anymore to sustain anti-fungals production or other anaerobic bacteria can outcompete our <i>P. putida</i>. In figure 2 the growth rate in the active state is plotted as a function of the oxygen and carbon availability. If there is little oxygen available the growth rate drops and so will the production of the toxins. However the depletion of oxygen should not be a problem because there is a layer of oxygen present closely to the roots[17]. Because of this small oxygen layer and root exudates the BananaGuard will be contained close to the banana roots. This result is important for BananaGuard because it is designed to protect the banana roots from <i>F. oxysporum</i> invasion.
+
<i>P. putida</i> is a strict aerobic bacteria[13]. This means that the bacteria cannot grow without oxygen available. Another problem could be that it might not be possible anymore to sustain the fungal growth inhibitor production or other anaerobic bacteria can outcompete BananGuard. Figure 6 was made to check the dependancy for oxygen. In figure 6 the growth rate in the active state is plotted as a function of the oxygen and carbon availability. This showed a sharp slope for the oxygen uptake rate and and the growth rate for BananaGuard so this could be a huge bottleneck for BananaGuard. However, after a literature research we concluded that the depletion of oxygen should not be a problem because there is a layer of oxygen present closely to the roots [17].  
-
                        </p>
+
This is advantageous for BananaGuard, because this small oxygen layer and root exudates contains BananaGuard to the banana roots. This result is important for BananaGuard because it is designed to protect the banana roots from <i>F. oxysporum</i> invasion.
-
                        <figure>
+
</p>
-
                            <img src="https://static.igem.org/mediawiki/2014/thumb/6/63/Wageningen_UR_modeling_oxvscarbon.png/800px-Wageningen_UR_modeling_oxvscarbon.png" width="100%">
+
<figure>
-
                            <figcaption>Figure 3 relative growth rate is dependent on the oxygen and carbon uptake rate. Colors indicate the height of the growth rate. Note(1) that this is oxygen uptake rate and not oxygen availability, if the oxygen uptake rate is to high it needs to get rid of the oxygen this causes a decreased growth rate. Note(2) if the wild type growth rate is 0 h<sup>-1</sup>, the relative growth is 0 h<sup>-1</sup> as well.</figcaption>
+
<img src="https://static.igem.org/mediawiki/2014/6/63/Wageningen_UR_modeling_oxvscarbon.png" width="100%">
-
                        </figure>
+
<figcaption>Figure 6. Dependacy of oxygen and carbon uptake rates for the relative growth rate of BananaGuard in the active state. Colors indicate the height of the growth rate. Note(1) that this is oxygen uptake rate and not oxygen availability, if the oxygen uptake rate is to high it needs to get rid of the oxygen this causes a decreased growth rate. Note(2) if the wild type growth rate is 0 h<sup>-1</sup>, the relative growth is 0 h<sup>-1</sup> as well.</figcaption>
 +
</figure>
<br/>
<br/>
-
                        <h3> Parameter analysis results </h3>
+
<h3> Parameter analysis results </h3>
-
                        <p>
+
<p>
-
The parameters used to obtain the above-mentioned results are all based on available literature and published experimental data, except for the pyoverdin production rate. This parameter has been determined from growth experiments with <i>P. putida</i> (Link notebook walter pyoverdin). So every parameter is biological relevant, but because different studies report different values a estimation of a realistic production value has been calculated. This means that the parameter can fluctuate in a small range. To compensate this a parameter analysis has been performed.
+
The parameters used to obtain the above-mentioned results are all based on available literature and published experimental data, except for the pyoverdin production rate. This parameter has been determined from growth experiments with <i>P. putida</i> (<a href=”https://2014.igem.org/Team:Wageningen_UR/notebook/journal" target=”blank” class="soft_link">notebook</a>. So every parameter is biologically relevant, but because different studies report different values an estimation of a realistic production value has been calculated. This means that the parameter can fluctuate in a small range. To compensate this a parameter analysis has been performed.
-
Parameters with a relative small flux value are left out because of the small impact on the results. The other parameters have been checked over a range, dependent on the accuracy of the calculated parameter. The same analysis has been performed for all the fluctuations in the parameter set combined. This analysis shows the possible fluctuations of the modeled system and which parameter has the greatest impact on the system. Figure 3 shows that Pyoverdine and the protein production rate influences the growth rate of BananaGuard the most. From this analysis we can conclude that even in the worst case scenario BananaGuard still can perform its task around the banana roots.  
+
This analysis shows the possible fluctuations of the modeled system and which parameter has the greatest impact on the system. These parameters have been checked over a range, dependent on the accuracy of the calculated parameter. The same analysis has been performed for all the fluctuations in the parameter set combined.
-
                        </p>
+
Figure 7 shows that Pyoverdine and the protein production rate influences the growth rate of BananaGuard the most. From this analysis we can conclude that even in the worst case scenario BananaGuard still performs its task around the banana roots.  
-
                        <figure>
+
</p>
-
                            <img src="https://static.igem.org/mediawiki/2014/thumb/b/ba/Wageningen_UR_modeling_Parameters_Analysis.png/800px-Wageningen_UR_modeling_Parameters_Analysis.png" width="80%">
+
<figure>
-
                            <figcaption>Figure 4 different parameters and their influence on the relative growth rate compared with the WT <i>P. putida</i>. The triangle indicates the the lowest growth rate, the mean of all growth rates and the maximum growth rate of the resting state (red) and active state (green).
+
<img src="https://static.igem.org/mediawiki/2014/b/ba/Wageningen_UR_modeling_Parameters_Analysis.png" width="80%">
-
                                All: all parameters combined, Protein production: [0.1-0.4 μM min<sup>-1</sup>], Non-growth maintenance:[5-9mmol gDW<sup>-1</sup>hr<sup>-1</sup>], Cell density[0.035-0.047 gDW l<sup>-1</sup>], Pyoverdine[100-200 μmol gDW<sup>-1</sup>hr<sup>-1</sup>], DAPG[10-50 μmol gDW<sup>-1</sup>hr<sup>-1</sup>], DMDS[0.08-8 μmol gDW<sup>-1</sup>hr<sup>-1</sup>]
+
<figcaption>Figure 7. Different parameters and their influence on the relative growth rate compared with the WT <i>P. putida</i>. The tetragon indicates the lowest growth rate, the mean of all growth rates and the maximum growth rate of the resting state (red) and active state (green).
-
                            </figcaption>
+
All: all parameters combined, Protein production: [0.1-0.4 μM min<sup>-1</sup>], Non-growth maintenance:[5-9mmol gDW<sup>-1</sup>hr<sup>-1</sup>], Cell density[0.035-0.047 gDW l<sup>-1</sup>], Pyoverdine[100-200 μmol gDW<sup>-1</sup>hr<sup>-1</sup>], DAPG[10-50 μmol gDW<sup>-1</sup>hr<sup>-1</sup>], DMDS[0.08-8 μmol gDW<sup>-1</sup>hr<sup>-1</sup>]
-
                        </figure>
+
</figcaption>
-
                        </p>
+
</figure>
-
                      </section>
+
</p>
 +
</section>
<br/>
<br/>
<br/>
<br/>
-
                      <section id="performance1">
+
 
-
<br/>
+
<h1 id="performance" >System Performance</h1>                    
-
                        <h1>System performance over cell division</h1>
+
                                             
<br/>
<br/>
                         <h2 id="introduction3">Introduction</h2>
                         <h2 id="introduction3">Introduction</h2>
-
<p>In order to gain insight in the performance of our biological control agent in the soil the optimized system was modelled over time. Focusing on the robustness of the toggle-switch in the presence of fluctuating protein concentrations and the activation of the kill-switch upon a fusaric acid influx, the following system characteristics were taken into account:             
+
<p>In order to gain insight in the performance of our biological control agent in the soil the optimized system was modelled over time. Focusing on the robustness of the Kill-switch in the presence of fluctuating protein concentrations and the activation of the Kll-switch upon a fusaric acid influx, the following system characteristics were taken into account:<br/><br/>            
-
As a bi-stable system the kill-switch is subject to both extrinsic noise factors e.g. cell division and intrinsic noise e.g. fluctuations in repressor proteins1. This affects the stability of the toggle-switch within the kill-switch and thereby the system as a whole. Large fluctuations in the toggle-switch repressor protein concentrations will cause the kill-switch to change state intermittently <b>(Noise)</b>.
+
As a bi-stable system the kill-switch is subject to both extrinsic noise factors e.g. cell division and intrinsic noise e.g. fluctuations in repressor proteins. This affects the stability of the toggle-switch within the Kill-switch and thereby the system as a whole. Large fluctuations in the toggle-switch repressor protein concentrations will cause the Kill-switch to change state intermittently <b>(Noise)</b>.
-
The unwanted production of CIλ caused by leaky promoters is another factor contributing to kill-switch instability. A high basal production rate for CIλ will weigh the toggle-switch in favour of the antifungal producing LacI state <b>(Leakiness)</b>.
+
The unwanted production of CIλ caused by leaky promoters is another factor contributing to kill-switch instability. A high basal production rate for CIλ will weigh the toggle-switch in favour of the fungal growth inhibitor producing LacI state <b>(Leakiness)</b>.
-
A state that is vital for our biological control agent to work but only when induced in the presence of Fusarium oxysporum <b>(Induction)</b>.
+
A state that is vital for our biological control agent to work but only when induced in the presence of <i>F. oxysporum</i> oxysporum <b>(Induction)</b>.
-
To that end noise, leakiness and an influx of fusaric acid are the focal points in a stochastic model of our genetic circuit. The model aims to quantify kill-switch instability and characterize the biological control agents’ performance in the presence of Fusarium oxysporum. The model takes into account cell division dividing the cellular contents the moment the cell has doubled in size. The average growth rate provided by the metabolic model and the impact of anti-gene transfer toxin antitoxin system on cell growth rate.
+
To that end noise, leakiness and an influx of fusaric acid are the focal points in a stochastic model of our genetic circuit. The model aims to quantify Kill-switch stability and characterize BananaGuards performance in the presence of <i>F. oxysporum</i>. The model takes into account cell division dividing the cellular contents the moment the cell has doubled in size. The average growth rate provided by the metabolic model and the impact of anti-gene transfer toxin antitoxin system on cell growth rate (figure 8) [34,35,37].
</p>
</p>
<P>
<P>
<figure>
<figure>
                             <img src="https://static.igem.org/mediawiki/2014/0/03/Wageningen_UR_model_bob.png" width="70%">
                             <img src="https://static.igem.org/mediawiki/2014/0/03/Wageningen_UR_model_bob.png" width="70%">
-
                             <figcaption><b>Figure 1.</b> A schematic overview of the components contained within the model: inductuction of the fusaric acid promoter and thereby the system, stability of the kill-switch and impact of the toxin anti-toxin system on cell growth and population dynamics. .</figcaption>
+
                             <figcaption><b>Figure 8.</b> A schematic overview of the components contained within the model: inductuction of the fusaric acid promoter and thereby the system, stability of the kill-switch and impact of the toxin anti-toxin system on cell growth and population dynamics. .</figcaption>
                         </figure>
                         </figure>
</P>
</P>
-
                        </section>                       
+
                                             
<br/>
<br/>
                         <h2 id="approach3" >Theory and experimental design</h2>
                         <h2 id="approach3" >Theory and experimental design</h2>
<h3> Stochastic Modeling and Toxin-Antitoxin Systems </h3>
<h3> Stochastic Modeling and Toxin-Antitoxin Systems </h3>
<P>
<P>
-
The performance of the biological control agent is defined by its ability to destroy fusarium upon detection, its ability to compete with other micro-organisms and the stability of the genetic circuit.
+
The performance of BananaGuard is defined by its ability to destroy <i>F. oxysporum</i> upon detection, its ability to compete with other micro-organisms and the stability of the genetic circuit.
-
The ability to destroy Fusarium has been tested in vitro. The biological control agents’ ability to compete with other soil born micro-organisms has been modeled in the metabolic model. The stability of the genetic switch will be explored in the following model, there are <b3>three</b3> key factors in assessing this aspect of the biological control agent: <b3> Noise, Leakiness and Induction</b3> (figure 1)
+
The ability to destroy <i>F. oxysporum</i> has been tested in vitro. The biological control agents’ ability to compete with other soil born micro-organisms has been modeled in the metabolic model. The stability of the genetic switch will be explored in the following model, there are <b>three</b> key factors in assessing this aspect of the biological control agent: <b> Noise, Leakiness and Induction</b> (figure 9)
</P>
</P>
                         <figure>
                         <figure>
-
                             <img src="https://static.igem.org/mediawiki/2014/a/aa/Wagenigen_UR_modelstoch.png" width="70%">
+
                             <img src="https://static.igem.org/mediawiki/2014/thumb/f/fd/Wageningen_UR_model_bob_Figure55.png/800px-Wageningen_UR_model_bob_Figure55.png" width="100%">
-
                             <figcaption><b>Figure 1.</b> A schematic representation of the reengineered <i>P.putida</i>.<b>(I)</b> Noise can destabilize bistable systems, causing the kill-switch to change state. <b>(II)</b> A leaky promoter can destabalize the toggleswitch by weighing it to the LacI state. <b>(III)</b> the induction of the genetic circuit by fusaric acid determines the performance of <i>P.putida</i> as a biological control agent.</figcaption>
+
                             <figcaption><b>Figure 9.</b> A schematic representation of the reengineered <i>P. putida</i>.<b>(I)</b> Noise can destabilize bistable systems, causing the kill-switch to change state. <b>(II)</b> A leaky promoter can destabalize the toggleswitch by weighing it to the LacI state. <b>(III)</b> the induction of the genetic circuit by fusaric acid determines the performance of BananaGuard as a biological control agent.</figcaption>
                         </figure>
                         </figure>
Line 449: Line 443:
'intrinsic' noise derives from the natural variability in the way that molecules interact. This follows the observed difference in protein expression levels for identically regulated biological systems.
'intrinsic' noise derives from the natural variability in the way that molecules interact. This follows the observed difference in protein expression levels for identically regulated biological systems.
'Extrinsic' noise arises from random fluctuations in the environment or differences in the cellular sate. This comes in part from local pressure and temperature variations or the unequal distribution of cellular organelles after cell division1.
'Extrinsic' noise arises from random fluctuations in the environment or differences in the cellular sate. This comes in part from local pressure and temperature variations or the unequal distribution of cellular organelles after cell division1.
-
Our biological control agent will have to perform in these noisy environments. Recently, more and more experimental evidence in the form of bimodal population distributions indicates that noise plays a very important role in the switching of bistable systems2. In bi-stable systems given a sizeable noise term, prolonged exposure to a single state increases the chance of a state change. These phenomena could be detrimental to the functioning of our biological control agent. If a large part of our re-engineered Pseudomonas putida strain changes from the dormant state to the anti-fungal production state without an external input, the system would not function.
+
Our biological control agent will have to perform in these noisy environments. Recently, more and more experimental evidence in the form of bimodal population distributions indicates that noise plays a very important role in the switching of bistable systems. In bi-stable systems given a sizeable noise term, prolonged exposure to a single state increases the chance of a state change. These phenomena could be detrimental to the functioning of our biological control agent. If a large part of the BananaGuard population changes from the resting state (not induced by fusaric acid) to the fungal growth production state without an external input, the system would not function.
-
The following model therefore aims to quantify the instability of our biological control agent. If a large part of the population is robust to noisy environments we can confirm the functionality of the system.
+
The following model therefore aims to quantify the stability of our biological control agent. If a large part of the population is robust to noisy environments we can confirm the functionality of the system [34,35].
</P>
</P>
<P>
<P>
-
<b>Leakiness (II)</b> is another factor to contributing to the instability of our bi-stable system. Most toggle switch dependent systems do not include a third interaction factor interfering with the stability of one state. The repressed CIλ on the Input-output plasmid of the kill-switch will not be completely repressed. Biologically a repressor does not stay bound, there is a constant turnover of the repressor protein on the promoters3. If in the very short timeframe of promoter non-occupancy a RNA polymerase binds, CIλ will be transcribed. The binding of RNA polymerases occurs intermittently and leads to a basal production rate of CIλ. A high enough basal production rate will weigh the toggle-switch in favour of the anti-fungal LacI state, in effect stopping the biological control agent from functioning. The model therefore aims to investigate the tolerance levels for the leakiness of CIλ and characterize the conditions in which it has the greatest impact.
+
<b>Leakiness (II)</b> is another factor to contributing to the instability of our bi-stable system. Most toggle switch dependent systems do not include a third interaction factor interfering with the stability of one state. The repressed CIλ on the input/output plasmid of the Kill-switch will not be completely repressed. Biologically a repressor does not stay bound, there is a constant turnover of the repressor protein on the promoters. If in the very short timeframe of promoter non-occupancy a RNA polymerase binds, CIλ will be transcribed. The binding of RNA polymerases occurs intermittently and leads to a basal production rate of CIλ. A high enough basal production rate weighs the toggle-switch in favour of the fungal inhibiting LacI state, in effect stopping the biological control agent from functioning. The model therefore aims to investigate the tolerance levels for the leakiness of CIλ and characterize the conditions in which it has the greatest impact [37].
</P>
</P>
<p>
<p>
-
<b>Induction (III)</b> ties into the performance aspect of our biological control. The biological control agent has to perform in the presence of Fusarium oxysporum. It is imperative that the Kill-switch system is activated in upon an influx of fusaric acid. In the soil the influx of fusaric acid will not be continuous but rather interlaced with periods low induction. Most likely leading to an increasing but haggard production of CIλ. The model therefore aims to characterize the optimal conditions for the successful activation of the anti-fungal producing state.
+
<b>Induction (III)</b> ties into the performance aspect of our biological control. The biological control agent has to perform in the presence of <i>F. oxysporum</i>. It is imperative that the Kill-switch system is activated upon an influx of fusaric acid. In the soil the influx of fusaric acid will not be continuous but rather interlaced with periods low induction. Most likely leading to an increasing but haggard production of CIλ. The model therefore aims to characterize the optimal conditions for the successful activation of the fungal growth inhibitor producing state [34].
</P>
</P>
<br/>
<br/>
<h3> Toxin-Antitoxin Systems </h3>
<h3> Toxin-Antitoxin Systems </h3>
<P>
<P>
-
Toxin-antitoxin systems in micro-organisms can act as a survival mechanism for micro-organisms in stressful environments. Operating as an hysteric switch, a small percentage of the population enters a toxin producing state. In this state protein production and cell growth are inhibited. This ‘metabolic inhibition’ makes these persistor cells robust to unfavorable conditions in the environment. Fluctuations in the toxins produced interfere with both the growth rate and protein translation of our biological control agent. This is therefore an important aspect of the stochastic model. Combining the kill-switch with the toxin antitoxin system in a model over cell division will give insight in the effects of varying cell division times on the different destabilizing factors and the induction of the biological control agent by fusaric acid. It will also assess the population dynamics of our biological control agent6.
+
Toxin-antitoxin systems in micro-organisms can act as a survival mechanism for micro-organisms in stressful environments. Operating as an hysteric switch, a small percentage of the population enters a toxin producing state. In this state protein production and cell growth are inhibited causing the cells to become persistor cells. This ‘metabolic inhibition’ makes these persistor cells robust to unfavorable conditions in the environment. Fluctuations in the toxins produced interfere with both the growth rate and protein translation of our biological control agent. This is therefore an important aspect of the stochastic model. Combining the Kill-switch with the toxin antitoxin system in a model over cell division will give insight in the effects of varying cell division times on the different destabilizing factors and the induction of the biological control agent by fusaric acid. It will also outline the population dynamics of our biological control agent. It is important to note that our chassis <i>P. putida</i> already has six native toxin-antitoxin systems. These have been incorporated into the model [38,39].
</P>
</P>
<br/>
<br/>
<h3> Model Equations </h3>
<h3> Model Equations </h3>
<P>
<P>
-
The model consists of eleven differential equations. Equations 1-6 describe the dynamics of the different repressor proteins, in the form of a Hill function. All proteins have an equation describing RNA transcription (equations) and translation (<b>term I</b>), protein degradation (<b>term II</b>) and the stochastic noise term (d = 0.2) (<b>term III</b>). The term below the translation describes the impact of the free proteolytic proteins from the toxin anti-toxin system. Should the level of free toxin increase the value below the fraction will reduce protein translation efficacy.
+
The model consists of eleven differential equations. The first six equations describe the dynamics of the different repressor proteins, in the form of a Hill function. All proteins have an equation describing RNA transcription (equations) and translation (<b>term I</b>), protein degradation (<b>term II</b>) and the stochastic noise term (d = 0.2) (<b>term III</b>). The term below the translation describes the impact of the free proteolytic proteins from the toxin anti-toxin system. Should the level of free toxin increase the value below the fraction will reduce protein translation efficacy [10,38,40,41].
</p>
</p>
 +
<br/>
<P>
<P>
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><msub><mrow><mi mathvariant="bold-italic">I</mi></mrow><mrow><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">P</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">o</mi><msub><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi><mo>_</mo><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub></mrow><mrow><mn>1</mn><mo>+</mo><mo>(</mo><mfrac><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow></msub></mrow></mfrac><msup><mrow><mo>)</mo></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow></msub></mrow></msup></mrow></mfrac><mo>-</mo><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><msub><mrow><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi><mo>_</mo><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><msub><mrow><mi mathvariant="bold-italic">I</mi></mrow><mrow><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub></math>
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><msub><mrow><mi mathvariant="bold-italic">I</mi></mrow><mrow><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">P</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">o</mi><msub><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi><mo>_</mo><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub></mrow><mrow><mn>1</mn><mo>+</mo><mo>(</mo><mfrac><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow></msub></mrow></mfrac><msup><mrow><mo>)</mo></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow></msub></mrow></msup></mrow></mfrac><mo>-</mo><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><msub><mrow><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">I</mi><mo>_</mo><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub><mi mathvariant="bold-italic">L</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">c</mi><msub><mrow><mi mathvariant="bold-italic">I</mi></mrow><mrow><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub></math>
Line 481: Line 476:
</P>
</P>
<P>
<P>
-
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">C</mi><msub><mrow><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow><mrow><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="bold-italic">s</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">P</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">o</mi><msub><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow></msub></mrow><mrow><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="bold-italic">s</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic"> </mi><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="bold-italic">s</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">d</mi></mrow></msub></mrow></mfrac><mo>-</mo><mi mathvariant="bold-italic"> </mi><mo>-</mo><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><msub><mrow><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">λ</mi><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mo>_</mo><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi><mo>_</mo><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></math>
+
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">C</mi><msub><mrow><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow><mrow><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="bold-italic">s</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">P</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">o</mi><msub><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow></msub></mrow><mrow><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="bold-italic">s</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic"> </mi><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="bold-italic">s</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">d</mi></mrow></msub></mrow></mfrac><mo>-</mo><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><msub><mrow><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi><mo>_</mo><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi><mo>_</mo><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></math>
</P>
</P>
<P>
<P>
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">n</mi><msub><mrow><mi mathvariant="bold-italic">s</mi></mrow><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi><mo>_</mo><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub><mi mathvariant="bold-italic">C</mi><msub><mrow><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow><mrow><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub></mrow><mrow><mn>1</mn><mo>+</mo><mfrac><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><msup><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi></mrow></msup></mrow><mrow><mi mathvariant="bold-italic">K</mi><msubsup><mrow><mi mathvariant="bold-italic">T</mi></mrow><mrow><mn>2</mn></mrow><mrow><mi mathvariant="bold-italic">n</mi></mrow></msubsup></mrow></mfrac></mrow></mfrac><mi mathvariant="bold-italic"> </mi><mo>-</mo><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><msub><mrow><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi><mo>_</mo><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi><mo>_</mo><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></math>
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">n</mi><msub><mrow><mi mathvariant="bold-italic">s</mi></mrow><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi><mo>_</mo><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub><mi mathvariant="bold-italic">C</mi><msub><mrow><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow><mrow><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub></mrow><mrow><mn>1</mn><mo>+</mo><mfrac><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><msup><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi></mrow></msup></mrow><mrow><mi mathvariant="bold-italic">K</mi><msubsup><mrow><mi mathvariant="bold-italic">T</mi></mrow><mrow><mn>2</mn></mrow><mrow><mi mathvariant="bold-italic">n</mi></mrow></msubsup></mrow></mfrac></mrow></mfrac><mi mathvariant="bold-italic"> </mi><mo>-</mo><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><msub><mrow><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi><mo>_</mo><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi><mo>_</mo><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></math>
</P>
</P>
 +
<br/>
<P>
<P>
-
The toxin anti-toxin systems can be described by two equations taken from the paper ‘Molecular mechanisms of multiple toxin–antitoxin systems are coordinated to govern the persister phenotype’ by Rick A. Fasani and Michael A. Savageau (2013). They proposed a mathematical description of a hysteric switch, linking it to fluctuating levels in free toxins and antitoxins and the growth rate. If the levels of the free toxin increases the growth rate of the cell decreases. The protein dilution caused by the growth rate was the main factor of protein degradation. This means that if a balanced is tipped the cell changes state and becomes a slow growing persistor cell.<br/><br/>
+
The toxin anti-toxin systems can be described by two equations taken from the paper ‘Molecular mechanisms of multiple toxin–antitoxin systems are coordinated to govern the persister phenotype’ by Rick A. Fasani and Michael A. Savageau (2013). They proposed a mathematical description of a hysteric switch, linking it to fluctuating levels in free toxins and antitoxins and the growth rate. If the levels of the free toxin increases the growth rate of the cell decreases. Protein dilution caused by the growth rate was the main factor of protein degradation. This means that if a balanced is tipped the cell changes state and becomes a slow growing persistor cell.<br/><br/>
-
Equation 7 describes the production of toxin RNA as part of the kill-switch similar to equations 1-6. Equations 8-9 describe the levels of free toxin and antitoxin through synthesis (<b>term I</b>), active degradation (<b>term II</b>), and stochastic fluctuations (<b>term III</b>). The first part in the denominator of the synthesis term describe the impact of protein degradation due to free toxin levels. The second part describes the auto repression of protein synthesis by the various complexes6.
+
The three equations that follow describe the production of toxin RNA as part of the Kill-switch and free toxins and antitoxins(<b>term I</b>), active degradation (<b>term II</b>), and stochastic fluctuations (<b>term III</b>). The equation describing free toxins has a RNA translation term added. The production terms for the free toxins and antitoxins have two components in the denominator. The first term describes the inhibition of protein synthesis caused by the free toxins. The second term describes auto repression of the toxin antitoxin genes by the toxin anti-toxin complexes that are formed [38].
</P>
</P>
<P>
<P>
 +
<br/>
<P>
<P>
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><msub><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">P</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">o</mi><msub><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">n</mi></mrow></msub></mrow><mrow><mn>1</mn><mo>+</mo><mo>(</mo><mfrac><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow></msub></mrow></mfrac><msup><mrow><mo>)</mo></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow></msub></mrow></msup><mo>+</mo><mo>(</mo><mfrac><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow></msub></mrow></mfrac><msup><mrow><mo>)</mo></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow></msub></mrow></msup></mrow></mfrac><mo>-</mo><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><msub><mrow><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><msub><mrow><mi mathvariant="bold-italic">x</mi></mrow><mrow><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub></mrow></msub><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><msub><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub><mi mathvariant="bold-italic"> </mi></math>
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><msub><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">P</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">o</mi><msub><mrow><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">n</mi></mrow></msub></mrow><mrow><mn>1</mn><mo>+</mo><mo>(</mo><mfrac><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow></msub></mrow></mfrac><msup><mrow><mo>)</mo></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">λ</mi></mrow></msub></mrow></msup><mo>+</mo><mo>(</mo><mfrac><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow></msub></mrow></mfrac><msup><mrow><mo>)</mo></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">e</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">R</mi></mrow></msub></mrow></msup></mrow></mfrac><mo>-</mo><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><msub><mrow><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><msub><mrow><mi mathvariant="bold-italic">x</mi></mrow><mrow><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub></mrow></msub><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><msub><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">R</mi><mi mathvariant="bold-italic">N</mi><mi mathvariant="bold-italic">A</mi></mrow></msub><mi mathvariant="bold-italic"> </mi></math>
Line 500: Line 497:
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mi mathvariant="bold-italic">α</mi><mi mathvariant="bold-italic">σ</mi></mrow><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><msup><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi></mrow></msup></mrow><mrow><mi mathvariant="bold-italic">K</mi><msubsup><mrow><mi mathvariant="bold-italic">T</mi></mrow><mrow><mn>2</mn></mrow><mrow><mi mathvariant="bold-italic">n</mi></mrow></msubsup></mrow></mfrac><mo>)</mo><mo>(</mo><mn>1</mn><mo>+</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><msup><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">s</mi></mrow></msup></mrow><mrow><mi mathvariant="bold-italic">K</mi><msubsup><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mn>1</mn></mrow><mrow><mi mathvariant="bold-italic">s</mi></mrow></msubsup></mrow></mfrac><mo>+</mo><mo>(</mo><mn>2</mn><mfrac><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><msup><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">s</mi></mrow></msup><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi></mrow><mrow><mi mathvariant="bold-italic">K</mi><msubsup><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mn>2</mn></mrow><mrow><mi mathvariant="bold-italic">s</mi></mrow></msubsup></mrow></mfrac><msup><mrow><mo>)</mo></mrow><mrow><mi mathvariant="bold-italic">P</mi></mrow></msup><mo>+</mo><mfenced separators="|"><mrow><mfrac><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><msup><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">s</mi></mrow></msup><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><msup><mrow><mi mathvariant="bold-italic">x</mi></mrow><mrow><mi mathvariant="bold-italic">s</mi></mrow></msup></mrow><mrow><mi mathvariant="bold-italic">K</mi><msubsup><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mn>1</mn></mrow><mrow><mi mathvariant="bold-italic">s</mi></mrow></msubsup><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">K</mi><msup><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mi mathvariant="bold-italic">s</mi></mrow></msup></mrow></mfrac></mrow></mfenced><mo>)</mo></mrow></mfrac><mo>-</mo><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><msub><mrow><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">n</mi></mrow></msub><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">n</mi><mo>+</mo><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic"> </mi><msub><mrow><mi mathvariant="bold-italic">ζ</mi></mrow><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">n</mi></mrow></msub></math>
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mi mathvariant="bold-italic">α</mi><mi mathvariant="bold-italic">σ</mi></mrow><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><msup><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi></mrow></msup></mrow><mrow><mi mathvariant="bold-italic">K</mi><msubsup><mrow><mi mathvariant="bold-italic">T</mi></mrow><mrow><mn>2</mn></mrow><mrow><mi mathvariant="bold-italic">n</mi></mrow></msubsup></mrow></mfrac><mo>)</mo><mo>(</mo><mn>1</mn><mo>+</mo><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><msup><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">s</mi></mrow></msup></mrow><mrow><mi mathvariant="bold-italic">K</mi><msubsup><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mn>1</mn></mrow><mrow><mi mathvariant="bold-italic">s</mi></mrow></msubsup></mrow></mfrac><mo>+</mo><mo>(</mo><mn>2</mn><mfrac><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><msup><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">s</mi></mrow></msup><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi></mrow><mrow><mi mathvariant="bold-italic">K</mi><msubsup><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mn>2</mn></mrow><mrow><mi mathvariant="bold-italic">s</mi></mrow></msubsup></mrow></mfrac><msup><mrow><mo>)</mo></mrow><mrow><mi mathvariant="bold-italic">P</mi></mrow></msup><mo>+</mo><mfenced separators="|"><mrow><mfrac><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><msup><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">s</mi></mrow></msup><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><msup><mrow><mi mathvariant="bold-italic">x</mi></mrow><mrow><mi mathvariant="bold-italic">s</mi></mrow></msup></mrow><mrow><mi mathvariant="bold-italic">K</mi><msubsup><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mn>1</mn></mrow><mrow><mi mathvariant="bold-italic">s</mi></mrow></msubsup><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic">K</mi><msup><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mi mathvariant="bold-italic">s</mi></mrow></msup></mrow></mfrac></mrow></mfenced><mo>)</mo></mrow></mfrac><mo>-</mo><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><msub><mrow><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">n</mi></mrow></msub><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">n</mi><mo>+</mo><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic"> </mi><mi mathvariant="bold-italic"> </mi><msub><mrow><mi mathvariant="bold-italic">ζ</mi></mrow><mrow><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">t</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">n</mi></mrow></msub></math>
</P>
</P>
 +
<br/>
<P>
<P>
-
Equation 10 represents cell growth. Cells grow exponentially at a rate of µ_max, the stochastic fluctuation of free toxin levels inhibit maximum growth. The model states that given an initial volume of 1, cell volume will increase according to µV where µ is µ_max divided by the free toxin inhibition term.<br/><br/>
+
dV/dt represents cell growth. Cells grow exponentially at a rate of µ_max, the stochastic fluctuation of free toxin levels inhibit maximum growth. The model states that given an initial volume of 1, cell volume will increase according to µV where µ is µ_max divided by the free toxin inhibition term [38].
-
The moment the cell has doubled in size (volume = 2), the cell contents are divided over the daughter cells following a binomial distribution. The simulation is restarted following one of these daughter cells. Should the kill-switch change state the initial conditions are reset to those at the start of the simulation
+
The moment the cell has doubled in size (volume = 2), the cell contents are divided over the daughter cells following a binomial distribution. The simulation is restarted following one of these daughter cells. Should the kill-switch change state the initial conditions are reset to those at the start of the simulation [42].
</P>
</P>
 +
<br/>
<P>
<P>
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math" display="block"><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">V</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mrow><mi mathvariant="bold-italic">µ</mi></mrow><mrow><mi mathvariant="bold-italic">m</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">x</mi></mrow></msub><mi mathvariant="bold-italic">V</mi></mrow><mrow><mn>1</mn><mo>+</mo><mfrac><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><msup><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi></mrow></msup></mrow><mrow><mi mathvariant="bold-italic">K</mi><msubsup><mrow><mi mathvariant="bold-italic">T</mi></mrow><mrow><mn>1</mn></mrow><mrow><mi mathvariant="bold-italic">n</mi></mrow></msubsup></mrow></mfrac></mrow></mfrac></math>
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math" display="block"><mi mathvariant="bold-italic"> </mi><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">V</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mrow><mi mathvariant="bold-italic">µ</mi></mrow><mrow><mi mathvariant="bold-italic">m</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">x</mi></mrow></msub><mi mathvariant="bold-italic">V</mi></mrow><mrow><mn>1</mn><mo>+</mo><mfrac><mrow><mi mathvariant="bold-italic">T</mi><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">x</mi><mi mathvariant="bold-italic">i</mi><msup><mrow><mi mathvariant="bold-italic">n</mi></mrow><mrow><mi mathvariant="bold-italic">n</mi></mrow></msup></mrow><mrow><mi mathvariant="bold-italic">K</mi><msubsup><mrow><mi mathvariant="bold-italic">T</mi></mrow><mrow><mn>1</mn></mrow><mrow><mi mathvariant="bold-italic">n</mi></mrow></msubsup></mrow></mfrac></mrow></mfrac></math>
</P>
</P>
 +
<br/>
<P>
<P>
-
The last equation describes the concentration of fusaric acid in the cell. Not much is known about fusaric acid dynamics besides its quantities are measured in the soil after infection. The fusaric acid detector experiments have given an indication with respect to the kinetics of the fusaric acid promoter. The required number of fusaric acid molecules for promoter activation is quite low, especially compared to the levels of fusaric acid that have been measured in the soil after infection5. The influx is stochastic and is weighed such that the concentration increases over time. This randomized increasing influx represents the varying influx of fusaric acid into our biological control agent most likely found in nature.
+
The last equation describes the concentration of fusaric acid in the cell. Not much is known about fusaric acid dynamics besides measured quantities in the soil after infection. The fusaric acid detector experiments have given an indication with respect to the kinetics of the fusaric acid promoter. The required number of fusaric acid molecules for promoter activation is quite low, especially compared to the levels of fusaric acid that have been measured in the soil after infection. The influx is stochastic and is weighed such that the concentration increases over time. This randomized increasing influx represents the varying influx of fusaric acid into our BananaGuard population.
</P>
</P>
 +
<br/>
<P>
<P>
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math" display="block"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="bold-italic">s</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">ζ</mi></mrow><mrow><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">f</mi><mi mathvariant="bold-italic">l</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="bold-italic">x</mi></mrow></msub><mo>-</mo><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><msub><mrow><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="bold-italic">s</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">d</mi></mrow></msub><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="bold-italic">s</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">d</mi></math>
<math xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:m="http://schemas.openxmlformats.org/officeDocument/2006/math" display="block"><mfrac><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="bold-italic">s</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">d</mi></mrow><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">t</mi></mrow></mfrac><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">ζ</mi></mrow><mrow><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">n</mi><mi mathvariant="bold-italic">f</mi><mi mathvariant="bold-italic">l</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="bold-italic">x</mi></mrow></msub><mo>-</mo><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">e</mi><msub><mrow><mi mathvariant="bold-italic">g</mi></mrow><mrow><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="bold-italic">s</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">d</mi></mrow></msub><mi mathvariant="bold-italic">F</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="bold-italic">s</mi><mi mathvariant="bold-italic">a</mi><mi mathvariant="bold-italic">r</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">c</mi><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">d</mi></math>
-
                         </section>
+
                          
<br/>
<br/>
 +
 +
</P>
 +
<h3> Model Sensitivity </h3>
 +
<P>
 +
Below a table for the parameter tolerances for the maintenance bi-stability of the kill-switch and functionality in the toxin antitoxin system. A local sensitivity analysis was done in order to find the ranges of the bi-stable region for this parameter set by running simulation in a set state of the Kill-Switch. If the parameters cause instability, the state at the end if the simulation changed. The values for the toxin-antitoxin tolerances were published by Rick A. Fasani and Michael A. Savageau (2013). The indicate the parameter ranges where hysteric switching occurs
 +
</P>
 +
 +
<figure>
 +
<figcaption>Table 2. The parameter tolerances for the maintenance of a bi-stable Kill-switch and hysteric switching in the toxin-antitoxin system .</figcaption>
 +
<img src="https://static.igem.org/mediawiki/2014/7/75/Wageningen_UR_model_bob_table.png" width="90%">
 +
</figure>
 +
 +
<h3>Software used<h3>
 +
<p>
 +
All models were made in Python (Python software foundation, Delaware, US) version 2.71 for Mac OS X version 10.5.8. Solving of ODE systems was done using the standard soln function in Python
 +
</p>
 +
                          
                          
                         <h2 id="results3">Results</h2>
                         <h2 id="results3">Results</h2>
Line 520: Line 538:
<P>
<P>
The results follow the trend of Noise, Leakiness and input. The first set of figures describe the Kill-switch system over cell division without the presence of the toxin anti-toxin system, this means the cell division times are equal for all simulations. The results display the intrinsic behavior of Kill-switch stability and performance.<br/><br/>
The results follow the trend of Noise, Leakiness and input. The first set of figures describe the Kill-switch system over cell division without the presence of the toxin anti-toxin system, this means the cell division times are equal for all simulations. The results display the intrinsic behavior of Kill-switch stability and performance.<br/><br/>
-
The second set of figures follows the model with the toxin-antitoxin system included. The varying growth and protein production rates caused by fluctuations in the free toxin concentration do impact the functionality but the overall system performs in line with the results obtained (Figure 1).
+
The second set of figures follows the model with the toxin-antitoxin system included. The varying growth and protein production rates caused by fluctuations in the free toxin concentration do impact the functionality but the overall system performs in line with the results obtained (Figure 10) [38].
</P>
</P>
                         <figure>
                         <figure>
-
                             <img src="https://static.igem.org/mediawiki/2014/thumb/9/96/Wagenigen_UR_model_bob_A.PNG/800px-Wagenigen_UR_model_bob_A.PNG" width="100%">
+
                             <img src="https://static.igem.org/mediawiki/2014/thumb/9/9d/Wageningen_UR_model_bob_Fig12.png/800px-Wageningen_UR_model_bob_Fig12.png" width="100%">
-
                             <figcaption><b>Figure 1.Plots showing the effect of noise and leakiness on the stability of the kill-switch ans the performance of the killswitch upon induction by fusric acid. Every plot point represent a percentage of unstable (leakyness and noise) or functioning (induction) kill-switches over 10000 simulations.
+
                             <figcaption><b>Figure 10. Plots showing the effect of noise and leakiness on the stability of the Kill-switch and the performance of the Killswitch upon induction by fusric acid. Every plot point represent a percentage of unstable (leakyness and noise) or functioning (induction) kill-switches over 10000 simulations.
<br/><br/>
<br/><br/>
-
<b>(A)</b> The plot shows the impact of a normally distributed noise term over a range of 120 to 240 nM/min for the toggle switch proteins LacI and TetR. 5% of the toggle-switches appears to become unstable given a noise term of 200 nM/min and a 45 minute division time. Should the division time be doubled to 90 minutes less than half a percent becomes unstable.
+
<b>(A)</b> The plot shows the impact of a normally distributed noise term over a range of 120 to 240 nM/min for the toggle switch proteins LacI and TetR. 5% of the Kill-switches becomes unstable given a noise term of 200 nM/min and a 45 minute division time. Should the division time be doubled to 90 minutes less than half a percent becomes unstable.
<br/><br/>
<br/><br/>
-
<b>(B</b> The plot shows The impact of different CIλ protein basal production rates on
+
<b>(B)</b> The plot shows the impact of different CIλ protein basal production rates on
-
toggle-switch stability. The noise term for the toggle-switch proteins LacI and TetR is 200nM. This corresponds to 5% and 0.5% instability due to noise (figure. 1A). The plots show a linear increase in instability for increasing LamdaCI basal production rates given a cell division time of 45 minutes. Proportionally, the increase in instability is a lot higher given a doubling in cell division time umax to 90 minutes. This can be atributed to a lower dilution rate allowing for a build-up of CIλ increasing instability.
+
toggle-switch stability. The noise term for the toggle-switch proteins LacI and TetR is 200nM. This corresponds to 5% and 0.5% instability due to noise (figure. 1A). The plots show a linear increase in instability for increasing CIλ basal production rates given a cell division time of 45 minutes. Proportionally, the increase in instability is a lot higher given a doubling in cell division time umax to 90 minutes. This can be atributed to a lower dilution rate allowing for a build-up of CIλ increasing instability.
<br/><br/>
<br/><br/>
-
<b>C</b> The plot shows the performance levels of the kill-switch upon induction by Fusaric acid. Simulated with four different Fusaric acid inputs and four different growth rates. Testing solely the performance no leakiness was taken into account, a noise term of 200nM was added for the toggle-switch proteins LacI and TetR. The figure show the influx of Fusaric acid is not limiting to the performance of the kill-switch. The repeated doubling of the maximum growth rate doubles performance. Indicating that the build-up and therefore production rate of CIλ is limiting to the functioning of the kill-switch</figcaption>
+
<b>(C)</b> The plot shows the performance levels of the kill-switch upon induction by fusaric acid. Simulated with four different fusaric acid inputs and four different growth rates. Testing solely the performance no leakiness was taken into account, a noise term of 200nM was added for the toggle-switch proteins LacI and TetR. The figure show the influx of fusaric acid is not limiting to the performance of the Kill-switch. The repeated doubling of the maximum growth rate doubles performance. Indicating that the build-up and therefore production rate of CIλ is limiting to the functioning of the Kill-switch.</figcaption>
                         </figure>
                         </figure>
                         <p></p>
                         <p></p>
Line 539: Line 557:
<P>
<P>
-
<b>Noise, A:</b> The significant increase in stability given a longer cell division time coincides with the fact that rapid protein dilution adds to the instability of the toggle switches. Biologically the ratio between the protein levels at steady state has decreased, the chance that the noise term will have an impact therefore increases, resulting in a higher number of toggle switches that have changed state. <br/><br/>
+
<b>Noise, A:</b> The significant increase in stability given a longer cell division time coincides with the fact that rapid protein dilution adds to the instability of the Kill-switches. Biologically the ratio between the protein levels at steady state has decreased, the chance that the noise term will have an impact therefore increases, resulting in a higher number of Kill-switches that have changed state. <br/><br/>
<b>Leakiness, B:</b> The significant increase in instability given a longer cell division and a leaky promoter can be explained in a reversed but similar manner to figure 1A. Lower dilution rates due to increased division times increases the amount of CIλ in the cell. Biologically, a more effective repression of the TetR promoter occurs due higher levels of CIλ. This lowers the steady state levels of the TetR protein decreasing the ratio of the protein levels at steady state, increasing the chance that the noise term will have an impact.<br/><br/>
<b>Leakiness, B:</b> The significant increase in instability given a longer cell division and a leaky promoter can be explained in a reversed but similar manner to figure 1A. Lower dilution rates due to increased division times increases the amount of CIλ in the cell. Biologically, a more effective repression of the TetR promoter occurs due higher levels of CIλ. This lowers the steady state levels of the TetR protein decreasing the ratio of the protein levels at steady state, increasing the chance that the noise term will have an impact.<br/><br/>
<b>Induction, C:</b> The plots correspond production rate and the activation of the CIλ repressor protein to be limiting to the performance of the kill-switch. Similar to the effect of the basal production of CIλ, the increased division times and thereby lowered dilution rates allow CIλ to accumulate. The accumulation gradually lowers the steady state ratios of the toggle switch proteins until the system changes state. The results indicate that for a parameter set aiming to describe low copy plasmid numbers, the performance of our biological control agent genetic circuitry is tied to its growth rate.
<b>Induction, C:</b> The plots correspond production rate and the activation of the CIλ repressor protein to be limiting to the performance of the kill-switch. Similar to the effect of the basal production of CIλ, the increased division times and thereby lowered dilution rates allow CIλ to accumulate. The accumulation gradually lowers the steady state ratios of the toggle switch proteins until the system changes state. The results indicate that for a parameter set aiming to describe low copy plasmid numbers, the performance of our biological control agent genetic circuitry is tied to its growth rate.
</P>
</P>
<P>
<P>
-
Overall the results for the kill-switch modeled over cell division have shown that a fast growth rate affects the ‘basal’ stability of the genetic circuit. High dilution rates increase the chance of instability due to noise. Leakiness and induction in and of the system are more likely to be effective at a slow growth rates. In both cases the build-up of CIλ over time is key in accessing the systems instability and performance.
+
Overall the results for the Kill-switch modeled over cell division have shown that a fast growth rate affects the ‘basal’ stability of the genetic circuit. High dilution rates increase the chance of instability due to noise. Leakiness and induction in and of the system are more likely to be effective at a slow growth rates. In both cases the build-up of CIλ over time is key in accessing the systems instability and performance.
</P>
</P>
<p></p>
<p></p>
Line 559: Line 577:
</p>
</p>
<p>
<p>
-
For both Leakiness and induction, the maximum and the predicted growth rate was simulated. The prior in order to ascertain the general behavior of the system, the latter in order to get a more accurate picture of our biological control agents performance in the soil.
+
For both Leakiness and induction, the maximum and the predicted growth rate was simulated. The prior in order to ascertain the general behavior of the system, the latter in order to get a more accurate picture of our biological control agents performance in the soil (figure 11).
The effect of just noise on the stability of the biological control agent has been neglected because its effect was statistically insignificant.
The effect of just noise on the stability of the biological control agent has been neglected because its effect was statistically insignificant.
</p>
</p>
                         <figure>
                         <figure>
-
                             <img src="https://static.igem.org/mediawiki/2014/thumb/c/c1/Wagenigen_UR_model_bob_resultsAB.png/800px-Wagenigen_UR_model_bob_resultsAB.png" width="100%">
+
                             <img src="https://static.igem.org/mediawiki/2014/1/11/Wageningen_UR_model_bob_Figh.png" width="100%">
-
                             <figcaption><b>Figure 2. The figures show the stability given a leaky promoter . Taking the form of a histogram, the y-axis shows the number of cell divisions that have occurred at a specific time point. The X-axis shows the number of minutes it takes for the cell to divide. The light green bars indicate a functioning system and a specific cell division time, coloured bars indicate non-functioning systems corresponding to a specific basal production rate.
+
                             <figcaption><b>Figure 11. The figures show the stability given a leaky promoter . Taking the form of a histogram, the y-axis shows the number of cell divisions that have occurred at a specific time point. The X-axis shows the number of minutes it takes for the cell to divide. The light green bars indicate a functioning system and a specific cell division time, coloured bars indicate non-functioning systems corresponding to a specific basal production rate.
<br/><br/>
<br/><br/>
-
<b>(A)</b> For a maximum growth rate of 45 minutes, a basal CIλ production of 50 nM/min or higher destabilizes the kill-switch in slow growing cells. The frequency of division times within the population of the cells is spread out. A spread over longer divsion times can be observed after an initial peak at 60 minutes. Protein dilution due to cell division is high, the population dynamics is not impacted. 20000 simulation were run <br/><br/>
+
<b>(A)</b> For a maximum growth rate corresponding to a division time of 45 minutes, a basal CIλ production of 50 nM/min or higher destabilizes the kill-switch in slow growing cells. The frequency of division times within the population of the cells is spread out. A spread over longer divsion times can be observed after an initial peak at 60 minutes. Protein dilution due to cell division is high, the population dynamics is not impacted. A total of 20000 simulation were run. <br/><br/>
-
(B) For a Maximum growth rate of 180 minutes the stability of the kill-switch a basal CIλ production of 50 nM/min or higher destabilizes the kill-switch. The population dynamics are affected. Low protein dilution due to slow growth causes the Kill-switch to leak toxin. Higher basal production rates compensate, increasing the average growth rate but also the instability. A total of 5000 simulation were run.
+
(B) For a maximum growth rate corresponding to a division time of 180 minutes a basal CIλ production of 50 nM/min or higher destabilizes the Kill-switch. The population dynamics are affected. Low protein dilution due to slow growth causes the Kill-switch to leak toxin. Higher basal production rates compensate, increasing the average growth rate but also the instability. A total of 5000 simulation were run.
  </figcaption>
  </figcaption>
                         </figure>
                         </figure>
Line 579: Line 597:
</p>
</p>
<p>
<p>
-
<b>Leakiness, B:</b> Because the population of cells averaged a significantly higher growth rate the percentage of unstable kill-switches increased proportionally too, once again confirming increased instability for lengthy cell divisions cycles. With respect population dynamics, the same logic that is applied to figure 1B applies to figure 1A. the effect however of a leaky CIλ can clearly be seen. A combination of an included toxin on the kill-switch, a promoter subject to leakiness and Lower dilution rates increase the levels of free toxin. Low basal production levels for CIλ increases TetR repression, alleviating toxin repression for which the CIλ cannot compensate if its concentration is to low. Higher levels of basal CIλ production reduce the levels of free toxins speeding up cell growth but decreasing stability. In order to avoid this problem the basal production rate of CIλ needs to be measured and adjusted if necessary, values lower than 50nM for the basal production appear to have no significant effect on the stability of the genetic circuit. For a low basal prodcution of CIλ entails a low basal production of toxin
+
<b>Leakiness, B:</b> Because the population of cells averaged a significantly higher growth rate the percentage of unstable kill-switches increased proportionally too, once again confirming increased instability for lengthy cell divisions cycles. With respect population dynamics, the same logic that is applied to figure 1B applies to figure 1A. the effect however of a leaky CIλ can clearly be seen. A combination of an included toxin on the kill-switch, a promoter subject to leakiness and Lower dilution rates increase the levels of free toxin. Low basal production levels for CIλ increases TetR repression, alleviating toxin repression for which the CIλ cannot compensate if its concentration is to low. Higher levels of basal CIλ production reduce the levels of free toxins speeding up cell growth but decreasing stability. In order to avoid this problem the basal production rate of CIλ needs to be measured and adjusted if necessary, values lower than 50nM for the basal production appear to have no significant effect on the stability of the genetic circuit. For a low basal production of CIλ entails a low basal production of toxin.
  </p>
  </p>
                         <figure>
                         <figure>
                             <img src="https://static.igem.org/mediawiki/2014/0/00/Wageningen_UR_model_bob_C.png" width="100%">
                             <img src="https://static.igem.org/mediawiki/2014/0/00/Wageningen_UR_model_bob_C.png" width="100%">
-
                             <figcaption><b>Figure 3.</b> <b>A</b>The figures show the performance of our biological control agent upon induction. All plots take the form of histograms, the y-axis shows the number of cell divisions that have occurred at a specific time point over 20000 simulations. The X-axis shows the number of minutes it takes for the cell to divide. The green bars indicate a functioning system and a specific cell division time, red bars indicate non-functioning systems given a specific cell division time. Two maximum growth rates uMax were investigated 45 minutes and 180 minutes.<br/><br/>
+
                             <figcaption><b>Figure 12.</b> <b>A</b>The figures show the performance of our biological control agent upon induction. All plots take the form of histograms, the y-axis shows the number of cell divisions that have occurred at a specific time point over 20000 simulations. The X-axis shows the number of minutes it takes for the cell to divide. The green bars indicate a functioning system and a specific cell division time, red bars indicate non-functioning systems given a specific cell division time. Two maximum growth rates uMax were investigated 45 minutes and 180 minutes.<br/><br/>
(A) For a maximum growth rate of 45 minutes 36% of the kill-switches activate, longer division times activate the cells more effectively.<br/><br/> (B) For a maximum growth rate of 180 minutes 98% of the kill-switches activate, longer division times activate the cells more effectively. </figcaption>
(A) For a maximum growth rate of 45 minutes 36% of the kill-switches activate, longer division times activate the cells more effectively.<br/><br/> (B) For a maximum growth rate of 180 minutes 98% of the kill-switches activate, longer division times activate the cells more effectively. </figcaption>
                         </figure>
                         </figure>
Line 590: Line 608:
<b>Induction, A and B:</b> The population dynamics for an induced system varies slightly. The same logic applies with respect to leaky toxins but the time-spans are different. Because the build-up of CIλ occurs faster in an induced system compared to a system with a basal transcription, the time frame in which the toxin promoter leaks is smaller. This results in faster average growth rate compared to a situation with leaky promoters.
<b>Induction, A and B:</b> The population dynamics for an induced system varies slightly. The same logic applies with respect to leaky toxins but the time-spans are different. Because the build-up of CIλ occurs faster in an induced system compared to a system with a basal transcription, the time frame in which the toxin promoter leaks is smaller. This results in faster average growth rate compared to a situation with leaky promoters.
Induction in and of itself works in accordance with the results obtained in figure 1. Slower growth rates activate the system more efficiently if a maximum growth rate 36% of the kill-switches is assumed, given the growth rate obtained from the metabolic model an efficiency of 98% is achieved.</p>
Induction in and of itself works in accordance with the results obtained in figure 1. Slower growth rates activate the system more efficiently if a maximum growth rate 36% of the kill-switches is assumed, given the growth rate obtained from the metabolic model an efficiency of 98% is achieved.</p>
-
<p>
+
 
-
Overall the results provided by the stochastic model gives an initial indication of the system’s efficiency by quantifying its stability and performance. It gives insight into the behaviors the system might display in the soil. Finally it has shown that the biological control agent though potentially slowed down by the toxin anti-toxin system can work.
+
</section>
-
</p>
+
</br>
</br>
-
<h1 id="conclusion">Conclusion</h1>
+
<section id="conclusion">
 +
<h1>Conclusion</h1>
<p>
<p>
-
The performance of BananaGuard as a biological control agent is defined by its ability to destroy fusarium upon detection, its ability to compete with other micro-organisms and the stability of the genetic circuit.
+
The performance of BananaGuard as a biological control agent is defined by its ability to destroy <i>F. oxysporum</i> upon detection, its ability to compete with other micro-organisms and the stability of the genetic circuit.
-
The ability to destroy Fusarium has been tested in vitro. The modeling has therefore been used to predict the answer to the remaining two questions. <br/>
+
The ability to destroy <i>F. oxysporum</i> has been tested in vitro. The modeling has therefore been used to predict the answer to the remaining two questions. <br/><br/>
 +
Can the BananaGuard system survive in the presence of competing micro-organisms?
 +
<br/><br/>
 +
Is the genetic circuit robust to biological fluctuations?
 +
<br/><br/>  
We have demonstrated that different modeling approaches can be taken and integrated with both the lab and one another. The statistical mechanics model, modeling <b>promoter design</b> has shown that a prediction regarding the functionality of the kill-switch can be made with respect to the different constituents on the promoter. The predicted configuration of repressor binding sites has been used in designing a second, optimized version of the kill-switch.
We have demonstrated that different modeling approaches can be taken and integrated with both the lab and one another. The statistical mechanics model, modeling <b>promoter design</b> has shown that a prediction regarding the functionality of the kill-switch can be made with respect to the different constituents on the promoter. The predicted configuration of repressor binding sites has been used in designing a second, optimized version of the kill-switch.
-
The model estimating the <b>cost</b> of the introduced genetic circuit has shown that in its resting state, BananaGuard is capable of competing with other micro-organisms in the soil on a metabolic level. The average growth rate obtained in this model has been used in the final model describing the <b>stability and performance</b> of BananaGuard. The results obtained have characterized several aspects of the kill-switch and the toxin anti-toxin system. First, noise can impact the stability of the kill-switch to a certain degree. Secondly, different basal production levels of CIλ can have different effects on the average cell growth and stability of the kill-switch, a point of attention for final construct of the system. Finally, the kill-switch will perform with 98% efficiency given the slow growth rate in the soil predicted by the metabolic model.
+
The model estimating the <b>cost</b> of the introduced genetic circuit has shown that in its resting state, BananaGuard is capable of competing with other micro-organisms in the soil on a metabolic level. The average growth rate obtained in this model has been used in the final model describing the <b>stability and performance</b> of BananaGuard. The results obtained have characterized several aspects of the kill-switch and the toxin anti-toxin system. First, noise can impact the stability of the kill-switch to a certain degree. Secondly, different basal production levels of CIλ can have different effects on population dynamics, cell growth and the stability of the kill-switch, a point of attention for final construct of the system. Finally, the kill-switch will perform with 98% efficiency given the slow growth rate in the soil predicted by the metabolic model.
</P>
</P>
-
 
</section>
</section>
-
</section>
 
-
 
<br/>
<br/>
 +
<p style="float:right"><b>Continue to <a href="https://2014.igem.org/Team:Wageningen_UR/project/interlab"  class="soft_link"> Interlab Study >></a> </b>
<br/>
<br/>
                             <h1 id="appendix1">Appendix</h1>
                             <h1 id="appendix1">Appendix</h1>
Line 612: Line 632:
<br/>
<br/>
                           <h2>Table of parameters</h1>
                           <h2>Table of parameters</h1>
-
                          <h3 id="parameters41">Metabolic modeling</h2>
 
<br/>
<br/>
 +
                        <h3 id="parameters41">Promoter Design</h3>
 +
<p>
 +
                        <ul class="table">
 +
                            <table width="800px">
 +
                                <tr>
 +
                                    <th> Designation </th>
 +
                                    <th> Value </th>
 +
                                    <th> Description </th>
 +
                                    <th> Reference </th>
 +
                                </tr>
 +
<tr>
 +
                                    <td>K<sub>prod<sub>LacI</sub></sub></td>
 +
                                    <td>1 µM min<sup>-1</sup></td>
 +
                                    <td>LacI production rate</td>
 +
                                    <td>[9]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>K<sub>prod<sub>TetR</sub></sub></td>
 +
                                    <td>1 µM min<sup>-1</sup></td>
 +
                                    <td>TetR production rate</td>
 +
                                    <td>[9]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>K<sub>prod<sub>CIλ</sub></sub></td>
 +
                                    <td>1 µM min<sup>-1</sup></td>
 +
                                    <td>CIλ production rate</td>
 +
                                    <td>[9]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>K<sub>prod<sub>RNAtoxin</sub></sub></td>
 +
                                    <td>0.75 µM min<sup>-1</sup></td>
 +
                                    <td>Toxin production rate</td>
 +
                                    <td>[9]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>Dil</td>
 +
                                    <td>0.023 min<sup>-1</sup> </td>
 +
                                    <td>Degradation due to cell division</td>
 +
                                    <td>[9]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>deg<sub>LacI</sub></td>
 +
                                    <td>0.034 min<sup>-1</sup></td>
 +
                                    <td>Degradation (Half-life LacI)</td>
 +
                                    <td>[3,9]</td>
 +
                                </tr>
 +
                                    <tr>
 +
                                    <td>deg<sub>TetR</sub></td>
 +
                                    <td>0.023 min<sup>-1</sup></td>
 +
                                    <td>Degradation (Half-life TetR)</td>
 +
                                    <td>[4,9,43]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>deg<sub>CIλ</sub></td>
 +
                                    <td>0.0011 min<sup>-1</sup></td>
 +
                                    <td>Degradation (Half-life CIλ)</td>
 +
                                    <td>[4,9,43]</td>
 +
                                </tr>
 +
<tr>
 +
                                    <td>deg<sub>Toxin</sub></td>
 +
                                    <td>0.0025 min<sup>-1</sup></td>
 +
                                    <td>Degradation (Half-life Toxin)</td>
 +
                                    <td>[5]</td>
 +
                                </tr>
 +
                                <tr>
 +
<tr>
 +
                                    <td>A</td>
 +
                                    <td>6.02e<sup>23</sup></td>
 +
                                    <td>Avogadro’s number</td>
 +
                                    <td>[9]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>V</td>
 +
                                    <td>1.5e<sup>-16</sup> dm<sup>-3</sup></td>
 +
                                    <td>Cell volume</td>
 +
                                    <td>[9]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>P</td>
 +
                                    <td>1900</td>
 +
                                    <td>Number of RNA polymerases</td>
 +
                                    <td>[9]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>N<sub>ns</sub></td>
 +
                                    <td>500,000</td>
 +
                                    <td>Non-specific binding sites</td>
 +
                                    <td>[1,2,9]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>K<sub>b</sub></td>
 +
                                    <td>1.38e<sup>-23</sup>m<sup>2</sup>kg s<sup>-</sup></td>
 +
                                    <td>Boltzmann constant</td>
 +
                                    <td>[9]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>T</td>
 +
                                    <td>298</td>
 +
                                    <td>Temperature in Kelvin</td>
 +
                                    <td>-</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>∆ε<sub>pd</sub><sup>S</sup></td>
 +
                                    <td>0.01 nM</td>
 +
                                    <td>Energy of the state, RNA polymerases</td>
 +
                                    <td>[1,2]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>∆ε<sub>d</sub><sup>TetR<sub>LacI</sub></sup></td>
 +
                                    <td>100 pM</td>
 +
                                    <td>Energy of the state, TetR repressor</td>
 +
                                    <td>[4,6,9]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>∆ε<sub>d</sub><sup>TetR<sub>Toxin</sub></sup></td>
 +
                                    <td>100 pM</td>
 +
                                    <td>Energy of the state, TetR repressor</td>
 +
                                    <td>[4,6,9]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>∆ε<sub>d</sub><sup>LacI<sub>TetR</sub></sup></td>
 +
                                    <td>1 pM</td>
 +
                                    <td>Energy of the state, LacI repressor</td>
 +
                                    <td>[3,9,43]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>∆ε<sub>d</sub><sup>CIλ<sub>TetR</sub></sup></td>
 +
                                    <td>1 pM</td>
 +
                                    <td>Energy of the state, CIλ repressor</td>
 +
                                    <td>[4,9]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>∆ε<sub>d</sub><sup>CIλ<sub>Toxin</sub></sup></td>
 +
                                    <td>1 pM</td>
 +
                                    <td>Energy of the state, CIλ repressor</td>
 +
                                    <td>[4,9]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>nTet<sub>LacI</sub></td>
 +
                                    <td>-</td>
 +
                                    <td>Repressor binding sites TetR binding PLacI</td>
 +
                                    <td>-</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>nTet<sub>Toxin</sub></td>
 +
                                    <td>-</td>
 +
                                    <td>Repressor binding sites TetR binding PToxin</td>
 +
                                    <td>-</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>nCIλ<sub>Toxin</sub></td>
 +
                                    <td>-</td>
 +
                                    <td>Repressor binding sites CIλ binding PToxin</td>
 +
                                    <td>-</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>nCIλ<sub>TetR</sub></td>
 +
                                    <td>-</td>
 +
                                    <td>Repressor binding sites CIλ binding PTetR</td>
 +
                                    <td>-</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>nLacI</td>
 +
                                    <td>-</td>
 +
                                    <td>Repressor binding sites Laci binding PTetR</td>
 +
                                    <td>-</td>
 +
                                </tr>
 +
                            </table>
 +
                        </ul>
 +
                        </p>
 +
<br/>
 +
 +
<h3 id="parameters42">System Cost</h2>
                         <p>
                         <p>
                         <ul class="table">
                         <ul class="table">
Line 681: Line 873:
                                     <td>150 μmol gDW<sup>-1</sup>hr<sup>-1</sup> </td>
                                     <td>150 μmol gDW<sup>-1</sup>hr<sup>-1</sup> </td>
                                     <td>Pyoverdine production rate</td>
                                     <td>Pyoverdine production rate</td>
-
                                     <td>[10,11,22]</td>
+
                                     <td>Experiment</td>
                                 </tr>
                                 </tr>
                                 <tr>
                                 <tr>
Line 717: Line 909:
                         </p>
                         </p>
<br/>
<br/>
-
                        <h3 id="parameters42">Dynamic modeling</h2>
+
 
 +
<h3 id="parameters43"> System Performance</h3>
 +
<p>
 +
                        <ul class="table">
 +
                            <table width="800px">
 +
                                <tr>
 +
                                    <th> Designation </th>
 +
                                    <th> Value </th>
 +
                                    <th> Description </th>
 +
                                    <th> Reference </th>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>K<sub>prod RNA<sub>LacI</sub></sub></td>
 +
                                    <td>0.75 µM min<sup>-1</sup></td>
 +
                                    <td>LacI production rate</td>
 +
                                    <td>[9,11]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>K<sub>prod RNA<sub>TetR</sub></sub></td>
 +
                                    <td>0.75 µM min<sup>-1</sup></td>
 +
                                    <td>TetR production rate</td>
 +
                                    <td>[9,11]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>K<sub>prod RNA<sub>CIλ</sub></sub></td>
 +
                                    <td>0.75 µM min<sup>-1</sup></td>
 +
                                    <td>CIλ production rate</td>
 +
                                    <td>[9,11]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>K<sub>prod<sub>RNAtoxin</sub></sub></td>
 +
                                    <td>0.75 µM min<sup>-1</sup></td>
 +
                                    <td>Toxin production rate</td>
 +
                                    <td>[9,11]</td>
 +
                                </tr>
 +
                                      <tr>
 +
                                    <td>deg<sub>ToxinRNA</sub></td>
 +
                                    <td>0.178 min<sup>-1</sup></td>
 +
                                    <td>Degradation (Half-life Toxin RNA)</td>
 +
                                    <td>[5]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>deg<sub>RNA</sub></td>
 +
                                    <td>0.82 min<sup>-1</sup> </td>
 +
                                    <td>Degradation (Half-life RNA)</td>
 +
                                    <td>[45]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>deg<sub>LacI</sub></td>
 +
                                    <td>0.034 min<sup>-1</sup></td>
 +
                                    <td>Degradation (Half-life LacI)</td>
 +
                                    <td>[3,44]</td>
 +
                                </tr>
 +
                                    <tr>
 +
                                    <td>deg<sub>TetR</sub></td>
 +
                                    <td>0.023 min<sup>-1</sup></td>
 +
                                    <td>Degradation (Half-life TetR)</td>
 +
                                    <td>[43]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>deg<sub>CIλ</sub></td>
 +
                                    <td>0.011 min<sup>-1</sup></td>
 +
                                    <td>Degradation (Half-life CIλ)</td>
 +
                                    <td></td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>deg<sub>Toxin</sub></td>
 +
                                    <td>0.000155 min<sup>-1</sup></td>
 +
                                    <td>Degradation (Half-life Toxin)</td>
 +
                                    <td>[38]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>deg<sub>Anti-Toxin</sub></td>
 +
                                    <td>0.0011 min<sup>-1</sup></td>
 +
                                    <td>Degradation (Half-life Anti-Toxin)</td>
 +
                                    <td>[38]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>Trans<sub>LacI</sub></td>
 +
                                    <td>0.12 µM min<sup>-1</sup></td>
 +
                                    <td>Translation Rate LacI</td>
 +
                                    <td>[9,45]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>Trans<sub>TetR</sub></td>
 +
                                    <td>0.08 µM min<sup>-1</sup></td>
 +
                                    <td>Translation Rate TetR</td>
 +
                                    <td>[9,45]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>Trans<sub>CIλ</sub></td>
 +
                                    <td>0.12 µM min<sup>-1</sup></td>
 +
                                    <td>Translation Rate CIλ</td>
 +
                                    <td>[9,45]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>Trans<sub>Zeta</sub></td>
 +
                                    <td>0.12 µM min<sup>-1</sup></td>
 +
                                    <td>Translation Rate Zeta</td>
 +
                                    <td>[9,45]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>K<sub>TetR</sub></td>
 +
                                    <td>4</td>
 +
                                    <td>Dissociation constant (TetR)</td>
 +
                                    <td>[4,6,9,43]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>K<sub>LacI</sub></td>
 +
                                    <td>1.75</td>
 +
                                    <td>Dissociation constant (LacI)</td>
 +
                                    <td>[3,9,11]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>K<sub>CIλ</sub></td>
 +
                                    <td>4</td>
 +
                                    <td>Dissociation constant (CIλ)</td>
 +
                                    <td>[4,43]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>N<sub>TetR</sub></td>
 +
                                    <td>3</td>
 +
                                    <td>Hill constant(TetR)</td>
 +
                                    <td>[1,2,9,10]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>N<sub>LacI</sub></td>
 +
                                    <td>2</td>
 +
                                    <td>Hill constant(LacI)</td>
 +
                                    <td>[3,9,43]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>K<sub>CIλ</sub></td>
 +
                                    <td>2</td>
 +
                                    <td>Hill constant(CIλ)</td>
 +
                                    <td>[4]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>α</td>
 +
                                    <td>1 nM min<sup>-1</sup></td>
 +
                                    <td>Free Toxin-Antitoxin production</td>
 +
                                    <td>[38]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>σ</td>
 +
                                    <td>10</td>
 +
                                    <td>Antitoxin factor</td>
 +
                                    <td>[38]</td>
 +
                                </tr>
 +
                         
 +
                                <tr>
 +
                                    <td>μ<sub>max</sub></td>
 +
                                    <td>0.0153-0.0038</td>
 +
                                    <td>Maximum cell growth rate</td>
 +
                                    <td><b>Metabolic Model</b></td>
 +
              </tr>
 +
<tr>
 +
                                    <td>KP<sub>1</sub></td>
 +
                                    <td>1000 nM</td>
 +
                                    <td>Dissociation constant(Toxin Repression)</td>
 +
                                    <td>[38]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>KP<sub>2</sub></td>
 +
                                    <td>10 nM</td>
 +
                                    <td>Dissociation constant(antiToxin Repression)</td>
 +
                                    <td>[38]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>KT<sub>1</sub></td>
 +
                                    <td>10</td>
 +
                                    <td>Dissociation constant(Growth inhibition)</td>
 +
                                    <td>[38]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>KT<sub>2</sub></td>
 +
                                    <td>100</td>
 +
                                    <td>Dissociation constant(Translation inhibition)</td>
 +
                                    <td>[38]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>n</td>
 +
                                    <td>6</td>
 +
                                    <td>Half of the total number of Toxin-Antitoxin systems</td>
 +
                                    <td>[38]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>p</td>
 +
                                    <td>2</td>
 +
                                    <td>scaling factor</td>
 +
                                    <td>[38]</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>KH</td>
 +
                                    <td>1000 nM</td>
 +
                                    <td>Dissociation constant(Antitoxin)</td>
 +
                                    <td>[38]</td>
 +
                                </tr>
 +
                            </table>
 +
                        </ul>
 +
                        </p>                             
 +
<br/>
                         </section>
                         </section>
Line 725: Line 1,118:
                         <h2>References</h2>
                         <h2>References</h2>
<br/>
<br/>
-
                         <ol>
+
                         <ol class="references">
<h4>Promoter Design</h4>
<h4>Promoter Design</h4>
-
<li>Bintu L, Buchler NE, Garcia HG, Gerland U, Hwa T, Kondev J,Kuhlman T, Phillips R: Transcriptional regulation by the numbers: models. Curr Opin Genet Dev 2005, 15:125-135.</li>
+
<li>Bintu L, Buchler NE, Garcia HG, Gerland U, Hwa T, Kondev J,Kuhlman T, Phillips R (2005): Transcriptional regulation by the numbers: models. Curr Opin Genet Dev, 15:125-135.</li>
-
<li> Bintu L, Buchler NE, Garcia HG, Gerland U, Hwa T, Kondev J,Kuhlman T, Phillips R: Transcriptional regulation by the numbers: application. Curr Opin Genet Dev 2005 (2), 15:125-135.</li>
+
<li> Bintu L, Buchler NE, Garcia HG, Gerland U, Hwa T, Kondev J,Kuhlman T, Phillips R (2005): Transcriptional regulation by the numbers: application. Curr Opin Genet Dev (2), 15:125-135.</li>
  <li>Mitchel Lewis (2005) The Lac repressor. C. R. Biologies 328 (2005) 521–548.</li>
  <li>Mitchel Lewis (2005) The Lac repressor. C. R. Biologies 328 (2005) 521–548.</li>
-
<li>Subhayu Basu1, Yoram Gerchman1, Cynthia H. Collins3, Frances H. Arnold3& Ron Weiss1,2: A synthetic multicellular system for programmed pattern formation. Nature 434 ,28 April 2005, 1130-113.</li>
+
<li>Subhayu Basu1, Yoram Gerchman, Cynthia H. Collins3, Frances H. Arnold& Ron Weiss,2: A synthetic multicellular system for programmed pattern formation. Nature 434 ,28 April 2005, 1130-113.</li>
<li>Falcon C.M and Matthews K.S. (2000) Operator DNA sequence Variation Enhances High Affinity Binding by Hinge Helix Mutants of Lactose Repressor Protein. Biochemistry. 39, 11074-11084.</li>
<li>Falcon C.M and Matthews K.S. (2000) Operator DNA sequence Variation Enhances High Affinity Binding by Hinge Helix Mutants of Lactose Repressor Protein. Biochemistry. 39, 11074-11084.</li>
<li> Christian Berens , Dirk Schnappinger and Wolfgang Hillen ; The Role of the Variable Region in Tet Repressor for Inducibility by Tetracycline, Biological Chemisttry Volume 272, Number 11, Issue of March 14, 1997 pp. 6936-6942.</li>
<li> Christian Berens , Dirk Schnappinger and Wolfgang Hillen ; The Role of the Variable Region in Tet Repressor for Inducibility by Tetracycline, Biological Chemisttry Volume 272, Number 11, Issue of March 14, 1997 pp. 6936-6942.</li>
  <li>ME346A Introduction to Statistical Mechanics – Wei Cai – Stanford University – Win 2011, Homepage; micro.stanford.edu/.</li>
  <li>ME346A Introduction to Statistical Mechanics – Wei Cai – Stanford University – Win 2011, Homepage; micro.stanford.edu/.</li>
<li>Lee T.I, Young R.A: Transcription of eukaryotic protein-coding genes. Annu Rev Genet 2000, 34:77-137. p.109.</li>
<li>Lee T.I, Young R.A: Transcription of eukaryotic protein-coding genes. Annu Rev Genet 2000, 34:77-137. p.109.</li>
-
<li>Milo et al. Nucl. Acids Res. (2010) 38: D750-D753. When using a specific entry from the database it is highly recommended that you also specify the BioNumbers 6 digit ID, e.g. "BNID 100986, Milo et al 2010".</li>
+
<li>Milo et al. Nucl. Acids Res. (2010) 38: D750-D753.</li>
<li>Timothy S. Gardner, Charles R. Cantor1 & James J. Collins; Construction of a genetic toggle switch inEscherichia coli, 20 January 2000, Nature 403, 339-342.</li>
<li>Timothy S. Gardner, Charles R. Cantor1 & James J. Collins; Construction of a genetic toggle switch inEscherichia coli, 20 January 2000, Nature 403, 339-342.</li>
  <li>iGEM, Aberdeen Schotland 2009, homepage: igem.org.</li>
  <li>iGEM, Aberdeen Schotland 2009, homepage: igem.org.</li>
-
 
<br/>
<br/>
 +
<h4>System Cost</h4>
<li>Bennett, R. A., & Lynch, J. M. (1981). Bacterial growth and development in the rhizosphere of gnotobiotic cereal plants. Journal of General Microbiology, 125(1), 95-102.</li>
<li>Bennett, R. A., & Lynch, J. M. (1981). Bacterial growth and development in the rhizosphere of gnotobiotic cereal plants. Journal of General Microbiology, 125(1), 95-102.</li>
<li>Molina, L., Ramos, C., Duque, E., Ronchel, M. C., Garcı́a, J. M., Wyke, L., & Ramos, J. L. (2000). Survival of< i> Pseudomonas putida</i> KT2440 in soil and in the rhizosphere of plants under greenhouse and environmental conditions. Soil Biology and Biochemistry, 32(3), 315-321.</li>
<li>Molina, L., Ramos, C., Duque, E., Ronchel, M. C., Garcı́a, J. M., Wyke, L., & Ramos, J. L. (2000). Survival of< i> Pseudomonas putida</i> KT2440 in soil and in the rhizosphere of plants under greenhouse and environmental conditions. Soil Biology and Biochemistry, 32(3), 315-321.</li>
Line 761: Line 1,154:
<li>BNID 100985, Milo et al 2010.</li>
<li>BNID 100985, Milo et al 2010.</li>
<li>Gross, H., & Loper, J. E. (2009). Genomics of secondary metabolite production by Pseudomonas spp. Natural product reports, 26(11), 1408-1446.</li>
<li>Gross, H., & Loper, J. E. (2009). Genomics of secondary metabolite production by Pseudomonas spp. Natural product reports, 26(11), 1408-1446.</li>
-
<li>Stoker, N. G., Fairweathe, N. F., & Spratt, B. G. (1982). Versatile low-copy-number plasmid vectors for cloning in< i> Escherichia coli</i>. Gene, 18(3), 335-341.</li>
+
<li>Stoker, N. G., Fairweathe, N. F., & Spratt, B. G. (1982). Versatile low-copy-number plasmid vectors for cloning in Escherichia coli. Gene, 18(3), 335-341.</li>
 +
<br/>
 +
<h4>System Design</h4>
 +
<li>Shu C-C, Chatterjee A, Dunny G, Hu W-S, Ramkrishna D (2011) Bistability versus Bimodal Distributions in Gene Regulatory Processes from Population Balance. PLoS Comput Biol 7(8): e1002140.</li>
 +
<li>Volfson, D. & Marciniak, J. et al.(2005) Origins of extrinsic variability in eukaryotic gene expression. Nature 21 December 2005</li>
 +
<li>Martinez Arias, A. & Hayward (2006), P. Filtering transcriptional noise during development: concepts and mechanisms. Nature Rev. Genet. 7, 34–44</li>
 +
<li>Pallavi Penumetcha, Kin Lau, Xiao Zhu, Kelly Davis, Todd T. Eckdahl, A. Malcolm Campbell.  (2010) Improving the Lac System for Synthetic Biology. BIOS. 7-15</li>
 +
<li>R.A. Fasani, M.A. Savageau  Molecular mechanisms of multiple toxin-antitoxin systems are coordinated to govern the persister phenotype (2013) .  Proc. Natl. Acad. Sci. USA, 110, pp. E2528–E2537</li>
 +
<li>Etienne Maisonneuve, Kenn Gerdes (2014), Molecular Mechanisms Underlying Bacterial Persisters. Cell Volume 157, Issue 3,  Pages 539–548</li>
 +
<li>Rué, P., & Garcia-Ojalvo, J. (2013). Modeling gene expression in time and space. Annual review of biophysics, 42, 605-627.</li>
 +
<li>Swain, P. S. (2004). Efficient attenuation of stochasticity in gene expression through post-transcriptional control. Journal of molecular biology, 344(4), 965-976</li>
 +
<li>Gonze, D. (2013). Modeling the effect of cell division on genetic oscillators. Journal of theoretical biology, 325, 22-33.</li>
 +
<li>Braun, D., Basu, S., & Weiss, R. (2005, March). Parameter estimation for two synthetic gene networks: a case study. In Acoustics, Speech, and Signal Processing, 2005. Proceedings.(ICASSP'05). IEEE International Conference on (Vol. 5, pp. v-769). IEEE.</li>
 +
<li>Purcell, O., Grierson, C. S., Bernardo, M., & Savery, N. J. (2012). Temperature dependence of ssrA-tag mediated protein degradation. J. Biol. Eng, 6(10).</li>
 +
<li>Alon, U. (2006). An introduction to systems biology: design principles of biological circuits. CRC press.
 +
Nevozhay, D., Adams, R. M., Murphy, K. F., Josić, K., & Balázsi, G. (2009). Negative autoregulation linearizes the dose–response and suppresses the heterogeneity of gene expression. Proceedings of the National Academy of Sciences, 106(13), 5123-5128.</li>
 +
<li>Braun, D., Basu, S., & Weiss, R. (2005, March). Parameter estimation for two synthetic gene networks: a case study. In Acoustics, Speech, and Signal Processing, 2005. Proceedings.(ICASSP'05). IEEE International Conference on (Vol. 5, pp. v-769). IEEE.</li>
 +
 
</ol>
</ol>
 +
 +
 +
 +
                         </section>
                         </section>
</html>
</html>
{{:Team:Wageningen_UR/templates/footer}}
{{:Team:Wageningen_UR/templates/footer}}

Latest revision as of 03:43, 18 October 2014

Wageningen UR iGEM 2014

Modeling Overview

The performance of BananaGuard as biological control agent is defined by three factors: its ability to destroy or inhibit Fusarium oxysporum upon detection, its ability to compete with other micro-organisms in the soil and the stability of the genetic circuit.

If the fungal growth inhibitors cannot destroy or inhibit F. oxysporum, the bananaplants will die. This aspect of BananaGuards performance is being tested in vitro. This leaves two remaining questions to be answered:

Can the BananaGuard system survive in the presence of competing micro-organisms?

Is the genetic circuit robust to biological fluctuations?

The performance and stability of the genetic circuit is crucial for BananaGuard to work. To prevent F. oxysporum from functioning, the fungal growth inhibitors in BananaGuard need to be produced when the fungi are present. This requires two aspects of the system to function properly: the toggle-switch within the Kill-switch system needs to be stable and the fusaric acid dependent promoter needs to activate the system in the presence of fusaric acid.

Optimizing promoter design constitutes the first step of BananaGuard development. Different configurations of promoter elements on the promoters of the Kill-switch were modeled using statistical mechanics. Parameters such as repressor binding strengths and repressor production rates were changed in each simulation. By scoring Kill-switch activation, the promoter element configuration with the most simulations showing a bi-stable switch were considered to be most likely to work. Including the anti-gene transfer inhibition effect on cell growth rates and the fusaric acid dependent promoter, this optimized version of the BananaGuard system was subsequently modeled over a population of dividing cells, assessing both its stability and performance in the form of a stochastic model.

Running thousands of simulations incorporating noise terms and promoter leakiness, the robustness of the BananaGuard system was quantified. Stochasticity in the system caused by cell division resulted in fractions of fast and slow growing cell populations with either stable or unstable Kill-switches. The same was considered for system performance; upon induction by fusaric acid a portion of the cell population activated their Kill-switch or remained in the same state. In order to make the simulations as accurate as possible a parameter for the maximum growth rate of BananaGuard was required. This parameter was obtained from a metabolic model that described the effects of competition between the BananaGuard system and micro-organisms in the rhizosphere.

To ascertain our hosts competitiveness the rate of cell division needs to be investigated on a metabolic level. Cell growth rate is a key factor in the competition between our engineered Pseudomonas putida and the wildtype P. putida found in the rhizosphere of the banana plant roots [1,2]. The integrated genetic circuit uses metabolic resources that would otherwise be dedicated to cellular maintenance or growth. Should this divergence of resources be too large our host will be outcompeted. A genome scale metabolic model was extended, investigating and outlining the cost of introducing our synthetic genes on a metabolic level.

Figure 1. Schematic overview of the modeling project indicating overlapping parts.


Kill-switch promoter design


Introduction

Promoter design for a stable Kill-switch involves balancing of different regulatory components. Two elements, the number of repressor binding sites and the position of those binding sites with respect to the gene, control two key variables: the chance a repressor binds to the regulatory region and the affinity with which it binds. These variables, in turn, determine the effectiveness of the repressor protein [1,2].
It has been shown that repressor proteins are more effective when there are more repressor binding sites in front of a gene, confirming that the increased chance of repressor protein binding leads to a significant increase in its effectiveness [1,2]. With respect to binding affinity, it has been shown for TetR proteins that a repressor binding site positioned in front of the open reading frame is more effective than those positioned at the back [6].
Because of slight variations in the position of a repressor binding site or variations in experimental conditions, a range of different repressor binding strength values are reported for the same repressor (sometimes varying several orders of magnitude)[4,5].
This means that case specific designs cannot be modeled accurately but it is possible to estimate the likelihood that the designed system will work. By taking all the possible repressor binding site combinations and varying the different binding affinities (thereby simulating their position) and production rates, the robustness of the system can determined given the repressor binding site combinations(figure 2). Statistical mechanics is an ideal tool because it can incorporate both repressor binding strength and binding sites. Continue to learn more about the statistical mechanics approach or view the results.

Figure 2. The figure shows the kill-switch with arbitrarily chosen repressor binding sites. (I) Altering the position of the TetR repressor binding sites change the effectiveness [6]; a large space of basepairs between the binding site and the LacI will result in weaker repression [5]. (II) Adding a LacI repressor binding site will increase the chance that LacI binds to the promoter region.


Theory and experimental design

Statistical Mechanics

Statistical mechanics provides a way to understand how the macroscopic properties of a system arise from the individual components within this system. In effect there are two states to consider in statistical mechanic descriptions. A macrostate specifies a system in terms of quantities that average over the microscopic constituents of the system. For example if you consider a closed container of gas (figure 3), quantities for the macrostate would be pressure or volume. This means that these macrostate quantities only work if you consider a system with a very large number of particles. Talking about a macrostate of a single molecule makes no sense since the terms are inherently descriptive of large systems. A microstate, however, does describe the properties of a single molecule e.g. kinetic energy and position. The macro- and micro-states come together in the key concept of statistical mechanics. A lot of different microstates can result in the same macrostate, but characterizing the macrostate imposes restraints on the possible microstates [7].

Figure 3. The image represents the macro and microstates applied in statistical mechanics. The box can be considered a macrostate e.g. volume, temperature, pressure. The molecules with different microstates e.g. kinetic energy, direction. The molecules determine the pressure and temperature but a given temperature provides information on the kinetic energy a molecule is most likely to have.


Modeling Biological Systems using Statistical Mechanics

The logic of macrostates and microstates can also be applied to cells, replacing the container of the previous example with the cell and the different particle constituents with RNA polymerases and repressor proteins. The condition and production states of the various proteins can be considered a macrostate. Modeling with statistical mechanics provides an insight in transcription regulation of a genetic circuit by generating ‘parameter-free’ predictions for the levels of gene expression. Assuming transcription of a gene to be proportional to the chance that a RNA polymerase binds to the promoter of interest a statistical mechanics model has two microstate outcomes [1,2]:

Non-specific binding: A state where none of the polymerases are bound to the gene of interest and distributed among non-specific binding sites .

Specific binding: A state where the promoter of the gene of interest is occupied and the remaining polymerases are distributed among the non-specific binding sites.
In order to determine the probability that either state will be realized we need to know the total number of possible microstates and the chance that a particular microstate will be achieved. For the non-specific binding by RNA polymerases combinatorics can be used [1]. The following equations will describe step for step how statistical mechanics can model biology, starting with combinatorics.

Nns!P!Nns-P! = Number of non-specific binding possibilities

Where P is the number of RNA polymerases and Nns the number of non-specific binding sites. Because there is a difference in the binding strength between specific and non-specific binding sites a Boltzmann weight needs to be included [1]:

ZP= Nns!P!Nns-P! e-PεpdNSKbT

To describes difference between strength of specific and non-specific binding a Boltzmann weight is added. For the non-specific binding of RNA polymerases the Boltzmann weight εpdNS is going to be very low compared to specific binding, compensating for the high number of binding site arrangements. In order to get the total statistical weight the microstate for RNA polymerase bound to the promoter needs to be added [1,2].

Ztotp =ZP +ZP-1e-εpdSKbT

εpdS is the specific binding energy of the RNA polymerase binding to the promoter. To obtain the chance an RNA polymerase is bound we divide the specific binding chance by the total chance [1,2].

Pbound= Z(P-1)e-εpdSKbTZtotp

By dividing top and bottom by Z(P-1)e(-(εpdS)/(KbT)) the equation takes the form of a Hill function where the assumption was made that the amount of polymerases is a lot smaller than the amount of non-specific binding sites [1,2].

Pbound = 11+ NnsPe-εpdSKbT

In this form a fold change of mRNA expression can be correlated to the chance that a polymerase binds (Nns/P e(-(∆εpdS)/(KbT))). If the term describing RNA polymerase binding equals zero, the the chance that a RNA polymerase binds is equal to one and there is no change in gene expression. This will most likely be the case in a system without interference. Most biological systems, however, are dependent on regulators, repression is the main method of gene regulation in our designed circuit (see Kill-switch). The effect of repressors can be taken up in previous equation by adding the microstate ‘repressor protein bound’ as a regulatory region to the total statistical weight (see equation below). The term ‘repressor binding’ states that if a repressor is bound a polymerase cannot bind. The equation assumes a high basal level of expression in the absence of repressors. Note that like the previous equation, it takes the form of a Hill function, the gene expression fold change is dependent on the number of repressors [1,2].

Fold Change= Pbound(R0)Pbound(R=0)=1+ PNnse-εpdSKbT1+ PNnse-εpdSKbT + (RNnse-εdrKbT)n

R is the number of repressor proteins and ∆ε_dr the specific binding energy of the repressor proteins, the same combinatorics rules apply repressor binding as for polymerase binding. n is the number of repression regulators should those regulators have the same binding energy. Because they act independently on whether or not an RNA polymerase can bind, the regulatory values are multiplied. The equation can now be incorporated into the model describing the Kill-Switch by tying production rates to the fold change of gene expression (Kprod x Fold change), i.e. the production rate will decrease proportionally to changes in gene expression.


Model Equations

Ordinary differential equations describe the genetic circuit, three of based on the statistical mechanics approach. The Degprotein and Dildiv terms represent the half-life and dilution rate of the proteins. Because the third equation (the CIλ production) relies on an external input and is not dependent on any of the other equations, the ordinary differential equations uses Michaelis-Menten kinetics to simulate a fusaric acid input. In every simulation fusaric acid is inputted at t = 10 hours [1,2,10].

dLacIdt= KprodLacI 1+ PNnse-εpdSKbT1+ PNnse-εpdSKbT+(AVTetRNnse-εdTetR_LacKbT)nTet_Lac- DegLacI LacI -DildivLacI

dLacIdt= KprodLacI 1+ PNnse-εpdSKbT1+ PNnse-εpdSKbT+(AVTetRNnse-εdTetR_LacKbT)nTet_Lac- DegLacI LacI -DildivLacI

dCIλdt= KprodCIλ [fusaric acid]fusaric acid+Kd-degCIλ CIλ -DildivCIλ

dToxindt= KprodLacI 1+ PNnse-εpdSKbT1+ PNnse-εpdSKbT+(AVTetRNnse-εdTetRtoxinKbT)nTetToxin(AVCIλNnse-εdCIλtoxinKbT)nCIλ_GFP- degToxin Toxin-Dildivtoxin

Software used

All models were made in Python (Python software foundation, Delaware, US) version 2.71 for Mac OS X version 10.5.8. Solving of ODE systems was done using the standard soln function in Python


Results

System Behaviour and Optimal Repressor Binding Site Configuration

Before the optimal configuration of repressor binding sites could be identified, the behaviour of the system had to be assigned scores representing a functioning or non-functioning Kill-switch. Three distinct behaviours could be identified:


0: The Kill-switch does not work, the toggle switch is out of balance and does not function.

1: The system performs less efficiently, though the toggle switch changes state, the toxin promoter is leaky.

2: The system performs to design, after a fusaric acid input the toggle switch changes state and toxin is produced when CIλ leaves the system.

In total eight (A,B,..,H) different configurations of repressor binding sites were tested (Figure 4). Two: F and H, were found to be produce Kill-switch systems robust towards variations in binding strengths and production rates. The result of the best scoring configuration H has been magnified (figure 5). The repressor binding site configuration corresponding to this map has the highest chance of working in an experimental set-up. The x-axis shows the fold change in protein production, the y-axis has increased binding strengths for individual repressors by an order of magnitude each row. The color map indicates that increases in repressor binding strengths can be compensated for by increased production rates of either the target gene or decreased production of the repressor itself. The practical implications for constructing the an optimized version of the Kill-switch were as follows: two TetR repressor binding sites on the LacI promoter, two LacI and one CIλ repressor binding site on the TetR promoter, two TetR and one CIλ repressor binding site on the toxin promoter.

Figure 4. Color maps indicating functioning and non-functioning systems. Each letter represents different repressor binding site configurations. Each small square within the color maps represents a score for a simulation of the system with a unique set of parameters. The colors correspond to a working system (2:red), leaky system (1:green) and a non-functioning system (0:blue). The table underneath shows the configuration of the repressor binding sites. The type of repressor binding site is specified followed by the promoter it is located on. The letters correspond to the colormaps


Figure 5. The best scoring repressor binding site configuration H. Each small square within the colour maps represents a score for a simulation of the system with a unique set of parameters. The x-axis shows the production rates in uM/min. The Y-axis has the repressor binding strengths. Using the lowest value found in literature [11], each row the value for th repressor binding strength increases by an order of magnitude.

Statistical Mechanics in Practice

There were two main reasons for using statistical mechanics:

(I) Statistical mechanics provides a way describe the system in terms of the number of repressors and repressor binding sites, their binding affinity and the position of the repressor binding on the promoter [1,2].

(II) In contrast to Michaelis-Menten kinetics, statistical mechanics offers more freedom studying the behavior of differently designed promoters.Michaelis-Menten requires the dissociation constants of the repressor proteins to be experimentally determined [1,2].

Based on the results of the model the experimentalists ordered oligonucleotides and assembled six promoters. Building two Kill-switches in parallel the chance of constructing a stable system increases. The assumption that the system is optimized and stable has been incorporated in the model quantifying stability .


System Cost


Introduction

BananaGuard was designed without knowing whether our engineered Pseudomonas putida would be outcompeted by other bacteria and fungi found in the rhizosphere of banana roots [12, 13]. For a successful application it is crucial that our strain can survive and compete near the banana plant roots. This seemed challenging as our system consists of multiple constitutively expressed genes controlling the kill switch and toxin-antitoxin systems. Moreover, when the system is switched to the active state by sensing fusaric acid, BananaGuard produces multiple fungal growth inhibitor enzymes and compounds. The The results display the intrinsic behavior of Kill-switch stability and performance.

integrated synthetic pathways use metabolic resources that would otherwise be dedicated to cellular maintenance or growth. The synthetic pathways therefore reduce the bacteria's potential to sustain themselves in the rhizosphere. Because of this reduced potential we wondered to what degree P. putida metabolism was affected by the introduced synthetic genes and if need be, to what extent the system could be changed to lessen the impact. To approach this problem we used a genome-scale metabolic model (GSMM) [14]. Our GSMM extended by the BananaGuard system predicts the interaction of our integrated plasmids with genes originating from endogenous pathways. With this prediction of the metabolic stress we gain insight in the metabolic burden caused by the integration of our plasmids. This will give an indication about BananaGuard if it will lose its capability to compete with other rhizosphere-populating micro-organisms in its resting state.
Continue to learn more about a genomic-scale metabolic model or view the results.


Theory and experimental design

In order to answer our questions with a GSMM we need a model that:

  • contains a comprehensive representation of P. putida metabolism
  • is capable of predicting growth rates in various environments
  • includes all reactions that are important for the BananaGuard system

For this reason a comparison was made between various available GSMMs for P. putida. We took a in-depth look in two models that were up-to-date, published and verified: iJP962 and iJN746[15,16]. For the necessary work to model the metabolic stress in P. putida, the best starting model needs to be chosen which has the most potential for accurately modeling the metabolic stress of our system.

Table 1. Comparison between iJP962 and iJN746

According to the previously stated requirements,the model iJP962 is most suitable for modeling our system. The biggest difference between these models shown in table 1 is the number of compartments. But this is not necessarily an advantage. Many of the reactions are over the internal membrane and not functional for the core intermediate metabolites. Added to this, iJN746 has a focus on compounds that are in a different class of metabolites than the fungal growth inhibitors we chose for our system [27]. Therefore iJP962 is more representative due to more reactions and more metabolites in a broader spectrum. The primary category of concern for fungal growth inhibitor production is the amino acid metabolism and iJP962 has the most reactions in this category [14,15,17]. neither model included all the specific reactions, metabolites and genes that are necessary for modeling the fungal growth inhibitor metabolites. These reactions need to be added manually to the GSMM.


Extending the model

Traditionally, metabolic models only predict metabolic fluxes. Although this is directly applicable to the metabolites that we want to produce, such as DMDS, this does not suffice to answer our questions regarding the impact of the whole system. Therefore, we added the production of all the proteins of interest as if they are metabolites. For this we used the protein sequences to determine the correct stoichiometric ratio of amino acids [18]. The other side-substrates/products are in the same ratio as described for production of proteins for the biomass reaction in the GSMM. In the model Promoter Design the protein production rate is estimated.

As our metabolites of interest (DAPG, DMDS, DMTS) were not yet included in the GSMM, we had to deduce their production pathways based on databases and literature information [19,20,21,22](pathways and metabolites). Furthermore, we took special care that all added pathways are stoichiometrically balanced, which is essential for proper and reliable integration in a GSMM. To get the most realistic results we need to set a biologically relevant production rate. This can be done by estimating a production range with in vivo data from literature.


Programs and toolboxes

Matlab[23] was used in combination with the Cobra toolbox[24] for handling the model. The Cobra toolbox is designed for working with genome-scale metabolic models. To solve the flux balance analysis a linear program solver is needed. Gurobi[25] was chosen because of its compatibility with Matlab and the Cobra toolbox.


Flux Balance Analysis

Flux Balance Analysis (FBA) is a framework for investigating the theoretical capabilities of a GSMM and can be solved using a linear programming solver. This is possible because of the linear objective function and linear constraints[4,6]. FBA is based on the following assumptions:

  1. Steady state: The fluxes are considered to have attained a static equilibrium value and do not change through time, this means that accumulation or dissipation is not possible.
  2. Thermodynamic constraints: (ir)reversibility of reactions.
  3. No enzyme saturation: The maximum flux for every reaction is not dependent on the concentration of the enzymes used for that reaction.

The FBA method uses a representation of the metabolic reaction network in the form of a stoichiometric matrix (S) where :

  • Each column corresponds to a reaction Ri
  • Each row corresponds to a metabolite Cj

The definition of S is :

Si,j=k  s.t.    k<  0  Cj  is  a  substrate  of  reaction  Ri  with  stoichiometry  kk>  0  Cj  is  a  product  of  reaction  Ri  with  stoichiometry  kk=0    Cj  does  not  participate  in  reaction  Ri

The FBA problem is then formulated as a maximization or minimization problem under the determined constraints:

Maximize cTvSubject to Sv=0LowerboundvUpperbound

Where :

  • v is the vector of to be determined reaction fluxes
  • c is a vector of constants defining the objective function
  • S is the stoichiometric matrix
  • Lowerbound and Upperbound are vectors of constraints indicating, minimal and maximal flux values for each reaction


Media composition

GSMMs consider all metabolic possibilities given the available substrates. By changing the substrates and their uptake rates in the media you can dictate what’s available in the model. However, for most metabolites the uptake rates are unknown. This can be solved setting these fluxes as unconstrained in the GSMM. This also enforces decision to not include all root components but only the most prevalent. You don’t want that rare compounds have a huge impact on the model predictions. This does imply that the model predictions could be slightly negative: with including the other available compounds the system may work even better.

To get the most realistic approximation for BananaGuard in the soil the M9 media composition was used with carbon as limiting nutrient. The carbon source is based on the exudates of banana roots and consists of 34% sucrose and 66% divided over multiple amino acids[15, 27].


Results

P. putida naturally habits the rhizosphere of plant roots. To check if BananaGuard is still viable in its resting and active state, a comparison is important between the wild type P. putida and BananaGuard. If the growth rate of BananaGuard approximates the growth rate of wild type P. putida in the resting state it is safe to assume that the disadvantage of the introduced synthetic pathway does not have a huge impact on competitiveness with other rhizosphere-populating microorganisms. This means that BananaGuard should be able to populate the rhizosphere of the banana roots. For the active state it is not essential to compete with other micro-organisms but the capability to produce all the fungal growth inhibitors with a constant rate around the roots, without overloading a metabolic pathway and causing shortage of intermediates.

We aim for the best approximation of the environment of the banana roots. To do this we need to estimate the carbon uptake rate of P. Putida. Because these units are not a standard in literature we have to calculate them. Within this GSMM it is possible to calculate the growth rate with known nutrient uptake rates. This can be reversed, so if we know the the doubling time of P. putida in the rhizosphere the carbon uptake rate can be calculated. This doubling time is 3 hours resulting in a growth rate of 0.23 h-1 [1]. From this we can calculate the carbon uptake rate in the rhizosphere which is around 4 mmol gDW-1hr-1.

Figure 5 depicts how the growth rate of BananaGuard compares to the growth rate of the wild type P. putida depending on the carbon uptake rate of the organism in the resting state and active state of the kill switch. Glucose was chosen as carbon source reference, because of the efficient degradation of glucose. BananaGuard is still able to grow at almost max growth rate (>99%) compared to the wild type in the resting state. We can conclude from this that BananaGuard is not outcompeted by other rhizosphere-populating microorganisms because of the synthetic pathway. When the kill switch is activated by fusaric acid, BananaGuard is still able to growth with 50-65% of the maximum growth of the wild type (Figure 5). This indicates that metabolic stress is not a bottleneck for the production of fungal growth inhibitors in our fusaric acid activated system.

Next we want to know if our conclusions made for glucose are also valid for the environment close to the banana roots. To approach this realistic as possible the carbon source is changed to have the same ratio of carbohydrates and amino acids as found in the exudates of banana roots [26]. This new carbon composition causes a reduced growth rate because of less efficient degradation of these compounds. In figure 5 the relative growth rate is shown for the resting and active state based on root exudate composition. The relative growth rate of BananaGuard around the roots has slightly decreased compared to glucose as carbon source but is still above 45%. This indicates it still has the capacity to grow and maintain the constant fungal growth inhibitor production with one condition that BananaGuard needs to stay close to the roots to maintain its carbon uptake rate. From this we can conclude that these changes in the carbon composition do not change our previously made conclusions .

Figure 5. The relative growth rate compared to the wild type P. putida for different carbon uptake rates. Two different carbon sources banana exudates and and glucose as reference. The expected carbon uptake rate of P. putida in the rhizosphere is indicated with transparent red.

Oxygen necessity in the soil

P. putida is a strict aerobic bacteria[13]. This means that the bacteria cannot grow without oxygen available. Another problem could be that it might not be possible anymore to sustain the fungal growth inhibitor production or other anaerobic bacteria can outcompete BananGuard. Figure 6 was made to check the dependancy for oxygen. In figure 6 the growth rate in the active state is plotted as a function of the oxygen and carbon availability. This showed a sharp slope for the oxygen uptake rate and and the growth rate for BananaGuard so this could be a huge bottleneck for BananaGuard. However, after a literature research we concluded that the depletion of oxygen should not be a problem because there is a layer of oxygen present closely to the roots [17]. This is advantageous for BananaGuard, because this small oxygen layer and root exudates contains BananaGuard to the banana roots. This result is important for BananaGuard because it is designed to protect the banana roots from F. oxysporum invasion.

Figure 6. Dependacy of oxygen and carbon uptake rates for the relative growth rate of BananaGuard in the active state. Colors indicate the height of the growth rate. Note(1) that this is oxygen uptake rate and not oxygen availability, if the oxygen uptake rate is to high it needs to get rid of the oxygen this causes a decreased growth rate. Note(2) if the wild type growth rate is 0 h-1, the relative growth is 0 h-1 as well.

Parameter analysis results

The parameters used to obtain the above-mentioned results are all based on available literature and published experimental data, except for the pyoverdin production rate. This parameter has been determined from growth experiments with P. putida (notebook. So every parameter is biologically relevant, but because different studies report different values an estimation of a realistic production value has been calculated. This means that the parameter can fluctuate in a small range. To compensate this a parameter analysis has been performed. This analysis shows the possible fluctuations of the modeled system and which parameter has the greatest impact on the system. These parameters have been checked over a range, dependent on the accuracy of the calculated parameter. The same analysis has been performed for all the fluctuations in the parameter set combined. Figure 7 shows that Pyoverdine and the protein production rate influences the growth rate of BananaGuard the most. From this analysis we can conclude that even in the worst case scenario BananaGuard still performs its task around the banana roots.

Figure 7. Different parameters and their influence on the relative growth rate compared with the WT P. putida. The tetragon indicates the lowest growth rate, the mean of all growth rates and the maximum growth rate of the resting state (red) and active state (green). All: all parameters combined, Protein production: [0.1-0.4 μM min-1], Non-growth maintenance:[5-9mmol gDW-1hr-1], Cell density[0.035-0.047 gDW l-1], Pyoverdine[100-200 μmol gDW-1hr-1], DAPG[10-50 μmol gDW-1hr-1], DMDS[0.08-8 μmol gDW-1hr-1]



System Performance


Introduction

In order to gain insight in the performance of our biological control agent in the soil the optimized system was modelled over time. Focusing on the robustness of the Kill-switch in the presence of fluctuating protein concentrations and the activation of the Kll-switch upon a fusaric acid influx, the following system characteristics were taken into account:

As a bi-stable system the kill-switch is subject to both extrinsic noise factors e.g. cell division and intrinsic noise e.g. fluctuations in repressor proteins. This affects the stability of the toggle-switch within the Kill-switch and thereby the system as a whole. Large fluctuations in the toggle-switch repressor protein concentrations will cause the Kill-switch to change state intermittently (Noise). The unwanted production of CIλ caused by leaky promoters is another factor contributing to kill-switch instability. A high basal production rate for CIλ will weigh the toggle-switch in favour of the fungal growth inhibitor producing LacI state (Leakiness). A state that is vital for our biological control agent to work but only when induced in the presence of F. oxysporum oxysporum (Induction). To that end noise, leakiness and an influx of fusaric acid are the focal points in a stochastic model of our genetic circuit. The model aims to quantify Kill-switch stability and characterize BananaGuards performance in the presence of F. oxysporum. The model takes into account cell division dividing the cellular contents the moment the cell has doubled in size. The average growth rate provided by the metabolic model and the impact of anti-gene transfer toxin antitoxin system on cell growth rate (figure 8) [34,35,37].

Figure 8. A schematic overview of the components contained within the model: inductuction of the fusaric acid promoter and thereby the system, stability of the kill-switch and impact of the toxin anti-toxin system on cell growth and population dynamics. .


Theory and experimental design

Stochastic Modeling and Toxin-Antitoxin Systems

The performance of BananaGuard is defined by its ability to destroy F. oxysporum upon detection, its ability to compete with other micro-organisms and the stability of the genetic circuit. The ability to destroy F. oxysporum has been tested in vitro. The biological control agents’ ability to compete with other soil born micro-organisms has been modeled in the metabolic model. The stability of the genetic switch will be explored in the following model, there are three key factors in assessing this aspect of the biological control agent: Noise, Leakiness and Induction (figure 9)

Figure 9. A schematic representation of the reengineered P. putida.(I) Noise can destabilize bistable systems, causing the kill-switch to change state. (II) A leaky promoter can destabalize the toggleswitch by weighing it to the LacI state. (III) the induction of the genetic circuit by fusaric acid determines the performance of BananaGuard as a biological control agent.

Noise (I), or variability of gene expression is of two kinds: 'intrinsic' noise derives from the natural variability in the way that molecules interact. This follows the observed difference in protein expression levels for identically regulated biological systems. 'Extrinsic' noise arises from random fluctuations in the environment or differences in the cellular sate. This comes in part from local pressure and temperature variations or the unequal distribution of cellular organelles after cell division1. Our biological control agent will have to perform in these noisy environments. Recently, more and more experimental evidence in the form of bimodal population distributions indicates that noise plays a very important role in the switching of bistable systems. In bi-stable systems given a sizeable noise term, prolonged exposure to a single state increases the chance of a state change. These phenomena could be detrimental to the functioning of our biological control agent. If a large part of the BananaGuard population changes from the resting state (not induced by fusaric acid) to the fungal growth production state without an external input, the system would not function. The following model therefore aims to quantify the stability of our biological control agent. If a large part of the population is robust to noisy environments we can confirm the functionality of the system [34,35].

Leakiness (II) is another factor to contributing to the instability of our bi-stable system. Most toggle switch dependent systems do not include a third interaction factor interfering with the stability of one state. The repressed CIλ on the input/output plasmid of the Kill-switch will not be completely repressed. Biologically a repressor does not stay bound, there is a constant turnover of the repressor protein on the promoters. If in the very short timeframe of promoter non-occupancy a RNA polymerase binds, CIλ will be transcribed. The binding of RNA polymerases occurs intermittently and leads to a basal production rate of CIλ. A high enough basal production rate weighs the toggle-switch in favour of the fungal inhibiting LacI state, in effect stopping the biological control agent from functioning. The model therefore aims to investigate the tolerance levels for the leakiness of CIλ and characterize the conditions in which it has the greatest impact [37].

Induction (III) ties into the performance aspect of our biological control. The biological control agent has to perform in the presence of F. oxysporum. It is imperative that the Kill-switch system is activated upon an influx of fusaric acid. In the soil the influx of fusaric acid will not be continuous but rather interlaced with periods low induction. Most likely leading to an increasing but haggard production of CIλ. The model therefore aims to characterize the optimal conditions for the successful activation of the fungal growth inhibitor producing state [34].


Toxin-Antitoxin Systems

Toxin-antitoxin systems in micro-organisms can act as a survival mechanism for micro-organisms in stressful environments. Operating as an hysteric switch, a small percentage of the population enters a toxin producing state. In this state protein production and cell growth are inhibited causing the cells to become persistor cells. This ‘metabolic inhibition’ makes these persistor cells robust to unfavorable conditions in the environment. Fluctuations in the toxins produced interfere with both the growth rate and protein translation of our biological control agent. This is therefore an important aspect of the stochastic model. Combining the Kill-switch with the toxin antitoxin system in a model over cell division will give insight in the effects of varying cell division times on the different destabilizing factors and the induction of the biological control agent by fusaric acid. It will also outline the population dynamics of our biological control agent. It is important to note that our chassis P. putida already has six native toxin-antitoxin systems. These have been incorporated into the model [38,39].


Model Equations

The model consists of eleven differential equations. The first six equations describe the dynamics of the different repressor proteins, in the form of a Hill function. All proteins have an equation describing RNA transcription (equations) and translation (term I), protein degradation (term II) and the stochastic noise term (d = 0.2) (term III). The term below the translation describes the impact of the free proteolytic proteins from the toxin anti-toxin system. Should the level of free toxin increase the value below the fraction will reduce protein translation efficacy [10,38,40,41].


dLacIRNAdt= ProdLacI_RNA1+(TetRKTetR)nTetR-degLacI_RNALacIRNA

dLacIdt=TransLacI_RNA LacIRNA 1+ToxinnKT2n-degLacILacI+ ζLacI

dTetRRNAdt= ProdTetR1+(LacIKLacI)nTetR+(λCIKCIλ)nCIλ-degTetR_RNATetRRNA

dTetRdt=TransTetR_RNA TetRRNA 1+ToxinnKT2n-degTetRTetR+ ζTetR

dCIλRNAdt= FusaricAcid ProdCIλFusaricAcid KFusaricAcid- degCIλ_RNA CIλ_RNA

dCIλdt= TransCIλ_RNACIλRNA1+ToxinnKT2n -degCIλ_RNA CIλ_RNA


The toxin anti-toxin systems can be described by two equations taken from the paper ‘Molecular mechanisms of multiple toxin–antitoxin systems are coordinated to govern the persister phenotype’ by Rick A. Fasani and Michael A. Savageau (2013). They proposed a mathematical description of a hysteric switch, linking it to fluctuating levels in free toxins and antitoxins and the growth rate. If the levels of the free toxin increases the growth rate of the cell decreases. Protein dilution caused by the growth rate was the main factor of protein degradation. This means that if a balanced is tipped the cell changes state and becomes a slow growing persistor cell.

The three equations that follow describe the production of toxin RNA as part of the Kill-switch and free toxins and antitoxins(term I), active degradation (term II), and stochastic fluctuations (term III). The equation describing free toxins has a RNA translation term added. The production terms for the free toxins and antitoxins have two components in the denominator. The first term describes the inhibition of protein synthesis caused by the free toxins. The second term describes auto repression of the toxin antitoxin genes by the toxin anti-toxin complexes that are formed [38].


dToxinRNAdt= ProdToxin1+(CIλKCIλ)nCIλ+(TetRKTetR)nTetR-degToxRNAToxinRNA

dToxindt= α(1+ToxinnKT2n)(1+ AntiToxinsKp1s+(2AntitoxinsToxKp2s)P+AntitoxinsToxsKp1s KHs)+ TransToxin_RNAToxin1+ToxinnKT2n-degToxin+ζToxin

dAntiToxindt=ασ(1+ToxinnKT2n)(1+ AntiToxinsKp1s+(2AntitoxinsToxKp2s)P+AntitoxinsToxsKp1s KHs)-degAntiToxinAntiToxin+ ζAntiToxin


dV/dt represents cell growth. Cells grow exponentially at a rate of µ_max, the stochastic fluctuation of free toxin levels inhibit maximum growth. The model states that given an initial volume of 1, cell volume will increase according to µV where µ is µ_max divided by the free toxin inhibition term [38]. The moment the cell has doubled in size (volume = 2), the cell contents are divided over the daughter cells following a binomial distribution. The simulation is restarted following one of these daughter cells. Should the kill-switch change state the initial conditions are reset to those at the start of the simulation [42].


dVdt=µmaxV1+ToxinnKT1n


The last equation describes the concentration of fusaric acid in the cell. Not much is known about fusaric acid dynamics besides measured quantities in the soil after infection. The fusaric acid detector experiments have given an indication with respect to the kinetics of the fusaric acid promoter. The required number of fusaric acid molecules for promoter activation is quite low, especially compared to the levels of fusaric acid that have been measured in the soil after infection. The influx is stochastic and is weighed such that the concentration increases over time. This randomized increasing influx represents the varying influx of fusaric acid into our BananaGuard population.


dFusaricAciddt=ζInflux-degFusaricAcidFusaricAcid

Model Sensitivity

Below a table for the parameter tolerances for the maintenance bi-stability of the kill-switch and functionality in the toxin antitoxin system. A local sensitivity analysis was done in order to find the ranges of the bi-stable region for this parameter set by running simulation in a set state of the Kill-Switch. If the parameters cause instability, the state at the end if the simulation changed. The values for the toxin-antitoxin tolerances were published by Rick A. Fasani and Michael A. Savageau (2013). The indicate the parameter ranges where hysteric switching occurs

Table 2. The parameter tolerances for the maintenance of a bi-stable Kill-switch and hysteric switching in the toxin-antitoxin system .

Software used

All models were made in Python (Python software foundation, Delaware, US) version 2.71 for Mac OS X version 10.5.8. Solving of ODE systems was done using the standard soln function in Python

Results

The results follow the trend of Noise, Leakiness and input. The first set of figures describe the Kill-switch system over cell division without the presence of the toxin anti-toxin system, this means the cell division times are equal for all simulations. The results display the intrinsic behavior of Kill-switch stability and performance.

The second set of figures follows the model with the toxin-antitoxin system included. The varying growth and protein production rates caused by fluctuations in the free toxin concentration do impact the functionality but the overall system performs in line with the results obtained (Figure 10) [38].

Figure 10. Plots showing the effect of noise and leakiness on the stability of the Kill-switch and the performance of the Killswitch upon induction by fusric acid. Every plot point represent a percentage of unstable (leakyness and noise) or functioning (induction) kill-switches over 10000 simulations.

(A) The plot shows the impact of a normally distributed noise term over a range of 120 to 240 nM/min for the toggle switch proteins LacI and TetR. 5% of the Kill-switches becomes unstable given a noise term of 200 nM/min and a 45 minute division time. Should the division time be doubled to 90 minutes less than half a percent becomes unstable.

(B) The plot shows the impact of different CIλ protein basal production rates on toggle-switch stability. The noise term for the toggle-switch proteins LacI and TetR is 200nM. This corresponds to 5% and 0.5% instability due to noise (figure. 1A). The plots show a linear increase in instability for increasing CIλ basal production rates given a cell division time of 45 minutes. Proportionally, the increase in instability is a lot higher given a doubling in cell division time umax to 90 minutes. This can be atributed to a lower dilution rate allowing for a build-up of CIλ increasing instability.

(C) The plot shows the performance levels of the kill-switch upon induction by fusaric acid. Simulated with four different fusaric acid inputs and four different growth rates. Testing solely the performance no leakiness was taken into account, a noise term of 200nM was added for the toggle-switch proteins LacI and TetR. The figure show the influx of fusaric acid is not limiting to the performance of the Kill-switch. The repeated doubling of the maximum growth rate doubles performance. Indicating that the build-up and therefore production rate of CIλ is limiting to the functioning of the Kill-switch.

Noise, A: The significant increase in stability given a longer cell division time coincides with the fact that rapid protein dilution adds to the instability of the Kill-switches. Biologically the ratio between the protein levels at steady state has decreased, the chance that the noise term will have an impact therefore increases, resulting in a higher number of Kill-switches that have changed state.

Leakiness, B: The significant increase in instability given a longer cell division and a leaky promoter can be explained in a reversed but similar manner to figure 1A. Lower dilution rates due to increased division times increases the amount of CIλ in the cell. Biologically, a more effective repression of the TetR promoter occurs due higher levels of CIλ. This lowers the steady state levels of the TetR protein decreasing the ratio of the protein levels at steady state, increasing the chance that the noise term will have an impact.

Induction, C: The plots correspond production rate and the activation of the CIλ repressor protein to be limiting to the performance of the kill-switch. Similar to the effect of the basal production of CIλ, the increased division times and thereby lowered dilution rates allow CIλ to accumulate. The accumulation gradually lowers the steady state ratios of the toggle switch proteins until the system changes state. The results indicate that for a parameter set aiming to describe low copy plasmid numbers, the performance of our biological control agent genetic circuitry is tied to its growth rate.

Overall the results for the Kill-switch modeled over cell division have shown that a fast growth rate affects the ‘basal’ stability of the genetic circuit. High dilution rates increase the chance of instability due to noise. Leakiness and induction in and of the system are more likely to be effective at a slow growth rates. In both cases the build-up of CIλ over time is key in accessing the systems instability and performance.

The stability and performance of the Biological Control Agent over Cell Division

The following results will outline the following:

(I) The impact of the toxin-antitoxin anti-gene gene transfer system on cell growth.

(II) The impact of the toxin antitoxin on promoter leakiness and induction.

For both Leakiness and induction, the maximum and the predicted growth rate was simulated. The prior in order to ascertain the general behavior of the system, the latter in order to get a more accurate picture of our biological control agents performance in the soil (figure 11). The effect of just noise on the stability of the biological control agent has been neglected because its effect was statistically insignificant.

Figure 11. The figures show the stability given a leaky promoter . Taking the form of a histogram, the y-axis shows the number of cell divisions that have occurred at a specific time point. The X-axis shows the number of minutes it takes for the cell to divide. The light green bars indicate a functioning system and a specific cell division time, coloured bars indicate non-functioning systems corresponding to a specific basal production rate.

(A) For a maximum growth rate corresponding to a division time of 45 minutes, a basal CIλ production of 50 nM/min or higher destabilizes the kill-switch in slow growing cells. The frequency of division times within the population of the cells is spread out. A spread over longer divsion times can be observed after an initial peak at 60 minutes. Protein dilution due to cell division is high, the population dynamics is not impacted. A total of 20000 simulation were run.

(B) For a maximum growth rate corresponding to a division time of 180 minutes a basal CIλ production of 50 nM/min or higher destabilizes the Kill-switch. The population dynamics are affected. Low protein dilution due to slow growth causes the Kill-switch to leak toxin. Higher basal production rates compensate, increasing the average growth rate but also the instability. A total of 5000 simulation were run.

Leakiness, A: The population dynamics are effected by free toxins from the toxin-antitoxin system. For a lower maximum growth rate and therefore a lower protein dilution rate the fluctuations in the level of free toxins is attributed to more than just stochasticity. Because a toxin is included on the kill-switch its promoter is subject to leakiness with a low amount of repressor protein present. With lower dilution rates, low amounts of CIλ are able to build up slower over time. The CIλ protein represses both the toxin and the TetR promoter on the toggle switch. This results in lower levels of TetR for which the CIλ cannot compensate over the course of its build-up. A deficiency that translates in a leaky toxin promoter causing a portion of the cells to distribute over longer division times. The effect of leakiness on the stability of the genetic circuit corresponds to the results obtained in figure 1.B. Over the total number simulations only 13% of the kill-switches became unstable given a basal production rate larger than 50nM and a maximum growth rate of 45 minutes. As expected all instances of instability occurred in the slow growing portion of the cell population.

Leakiness, B: Because the population of cells averaged a significantly higher growth rate the percentage of unstable kill-switches increased proportionally too, once again confirming increased instability for lengthy cell divisions cycles. With respect population dynamics, the same logic that is applied to figure 1B applies to figure 1A. the effect however of a leaky CIλ can clearly be seen. A combination of an included toxin on the kill-switch, a promoter subject to leakiness and Lower dilution rates increase the levels of free toxin. Low basal production levels for CIλ increases TetR repression, alleviating toxin repression for which the CIλ cannot compensate if its concentration is to low. Higher levels of basal CIλ production reduce the levels of free toxins speeding up cell growth but decreasing stability. In order to avoid this problem the basal production rate of CIλ needs to be measured and adjusted if necessary, values lower than 50nM for the basal production appear to have no significant effect on the stability of the genetic circuit. For a low basal production of CIλ entails a low basal production of toxin.

Figure 12. AThe figures show the performance of our biological control agent upon induction. All plots take the form of histograms, the y-axis shows the number of cell divisions that have occurred at a specific time point over 20000 simulations. The X-axis shows the number of minutes it takes for the cell to divide. The green bars indicate a functioning system and a specific cell division time, red bars indicate non-functioning systems given a specific cell division time. Two maximum growth rates uMax were investigated 45 minutes and 180 minutes.

(A) For a maximum growth rate of 45 minutes 36% of the kill-switches activate, longer division times activate the cells more effectively.

(B) For a maximum growth rate of 180 minutes 98% of the kill-switches activate, longer division times activate the cells more effectively.

Induction, A and B: The population dynamics for an induced system varies slightly. The same logic applies with respect to leaky toxins but the time-spans are different. Because the build-up of CIλ occurs faster in an induced system compared to a system with a basal transcription, the time frame in which the toxin promoter leaks is smaller. This results in faster average growth rate compared to a situation with leaky promoters. Induction in and of itself works in accordance with the results obtained in figure 1. Slower growth rates activate the system more efficiently if a maximum growth rate 36% of the kill-switches is assumed, given the growth rate obtained from the metabolic model an efficiency of 98% is achieved.


Conclusion

The performance of BananaGuard as a biological control agent is defined by its ability to destroy F. oxysporum upon detection, its ability to compete with other micro-organisms and the stability of the genetic circuit. The ability to destroy F. oxysporum has been tested in vitro. The modeling has therefore been used to predict the answer to the remaining two questions.

Can the BananaGuard system survive in the presence of competing micro-organisms?

Is the genetic circuit robust to biological fluctuations?

We have demonstrated that different modeling approaches can be taken and integrated with both the lab and one another. The statistical mechanics model, modeling promoter design has shown that a prediction regarding the functionality of the kill-switch can be made with respect to the different constituents on the promoter. The predicted configuration of repressor binding sites has been used in designing a second, optimized version of the kill-switch. The model estimating the cost of the introduced genetic circuit has shown that in its resting state, BananaGuard is capable of competing with other micro-organisms in the soil on a metabolic level. The average growth rate obtained in this model has been used in the final model describing the stability and performance of BananaGuard. The results obtained have characterized several aspects of the kill-switch and the toxin anti-toxin system. First, noise can impact the stability of the kill-switch to a certain degree. Secondly, different basal production levels of CIλ can have different effects on population dynamics, cell growth and the stability of the kill-switch, a point of attention for final construct of the system. Finally, the kill-switch will perform with 98% efficiency given the slow growth rate in the soil predicted by the metabolic model.


Continue to Interlab Study >>

Appendix


Table of parameters


Promoter Design

    Designation Value Description Reference
    KprodLacI 1 µM min-1 LacI production rate [9]
    KprodTetR 1 µM min-1 TetR production rate [9]
    KprodCIλ 1 µM min-1 CIλ production rate [9]
    KprodRNAtoxin 0.75 µM min-1 Toxin production rate [9]
    Dil 0.023 min-1 Degradation due to cell division [9]
    degLacI 0.034 min-1 Degradation (Half-life LacI) [3,9]
    degTetR 0.023 min-1 Degradation (Half-life TetR) [4,9,43]
    degCIλ 0.0011 min-1 Degradation (Half-life CIλ) [4,9,43]
    degToxin 0.0025 min-1 Degradation (Half-life Toxin) [5]
    A 6.02e23 Avogadro’s number [9]
    V 1.5e-16 dm-3 Cell volume [9]
    P 1900 Number of RNA polymerases [9]
    Nns 500,000 Non-specific binding sites [1,2,9]
    Kb 1.38e-23m2kg s- Boltzmann constant [9]
    T 298 Temperature in Kelvin -
    ∆εpdS 0.01 nM Energy of the state, RNA polymerases [1,2]
    ∆εdTetRLacI 100 pM Energy of the state, TetR repressor [4,6,9]
    ∆εdTetRToxin 100 pM Energy of the state, TetR repressor [4,6,9]
    ∆εdLacITetR 1 pM Energy of the state, LacI repressor [3,9,43]
    ∆εdCIλTetR 1 pM Energy of the state, CIλ repressor [4,9]
    ∆εdCIλToxin 1 pM Energy of the state, CIλ repressor [4,9]
    nTetLacI - Repressor binding sites TetR binding PLacI -
    nTetToxin - Repressor binding sites TetR binding PToxin -
    nCIλToxin - Repressor binding sites CIλ binding PToxin -
    nCIλTetR - Repressor binding sites CIλ binding PTetR -
    nLacI - Repressor binding sites Laci binding PTetR -


System Cost

    Designation Value Description Reference
    model iJP962 - Genomic-scale metabolic model [6,18]
    carbonsource - Fructose, Glucose, Alanine, Asparagine, Aspartate, L-Leucine, Serine, L-Threonine, L-Proline, D-Glutamate [15]
    glucuptake 4.43 mmol gDW-1hr-1 Carbon uptake rate [1]
    NGAM 3.96 mmol gDW-1hr-1 Non-growth maintenance factor [6]
    KS_rate 0.2 μM min-1 Protein production rate Promoter Design
    TAT_rate 1 nM min-1 Toxin/anti-toxin production rate Promoter Design
    Cell density 0.040 gDW l-1 gDW/OD600=0.1 [20]
    ncells 1e8 cells ml-1 cells/OD600=0.1 [21]
    Chitinase 0.2 μM min-1 Chitinase production rate [7] and Promoter Design
    Pyoverdine 150 μmol gDW-1hr-1 Pyoverdine production rate Experiment
    DAPG 25 μmol gDW-1hr-1 DAPG production rate [8]
    DMDS 0.845 μmol gDW-1hr-1 DMDS production rate [9]
    Nplasmids 10 plasmids cell-1 Number of plasmids per cell (low copy number) [23]
    plasmidTD 0.3 h-1 Time needed for at least a doubling of the number of plasmids per cell [1]
    TDrhizo 3 hr Doubling time P. putida in the rhizosphere [1]


System Performance

    Designation Value Description Reference
    Kprod RNALacI 0.75 µM min-1 LacI production rate [9,11]
    Kprod RNATetR 0.75 µM min-1 TetR production rate [9,11]
    Kprod RNACIλ 0.75 µM min-1 CIλ production rate [9,11]
    KprodRNAtoxin 0.75 µM min-1 Toxin production rate [9,11]
    degToxinRNA 0.178 min-1 Degradation (Half-life Toxin RNA) [5]
    degRNA 0.82 min-1 Degradation (Half-life RNA) [45]
    degLacI 0.034 min-1 Degradation (Half-life LacI) [3,44]
    degTetR 0.023 min-1 Degradation (Half-life TetR) [43]
    degCIλ 0.011 min-1 Degradation (Half-life CIλ)
    degToxin 0.000155 min-1 Degradation (Half-life Toxin) [38]
    degAnti-Toxin 0.0011 min-1 Degradation (Half-life Anti-Toxin) [38]
    TransLacI 0.12 µM min-1 Translation Rate LacI [9,45]
    TransTetR 0.08 µM min-1 Translation Rate TetR [9,45]
    TransCIλ 0.12 µM min-1 Translation Rate CIλ [9,45]
    TransZeta 0.12 µM min-1 Translation Rate Zeta [9,45]
    KTetR 4 Dissociation constant (TetR) [4,6,9,43]
    KLacI 1.75 Dissociation constant (LacI) [3,9,11]
    KCIλ 4 Dissociation constant (CIλ) [4,43]
    NTetR 3 Hill constant(TetR) [1,2,9,10]
    NLacI 2 Hill constant(LacI) [3,9,43]
    KCIλ 2 Hill constant(CIλ) [4]
    α 1 nM min-1 Free Toxin-Antitoxin production [38]
    σ 10 Antitoxin factor [38]
    μmax 0.0153-0.0038 Maximum cell growth rate Metabolic Model
    KP1 1000 nM Dissociation constant(Toxin Repression) [38]
    KP2 10 nM Dissociation constant(antiToxin Repression) [38]
    KT1 10 Dissociation constant(Growth inhibition) [38]
    KT2 100 Dissociation constant(Translation inhibition) [38]
    n 6 Half of the total number of Toxin-Antitoxin systems [38]
    p 2 scaling factor [38]
    KH 1000 nM Dissociation constant(Antitoxin) [38]




References


    Promoter Design

  1. Bintu L, Buchler NE, Garcia HG, Gerland U, Hwa T, Kondev J,Kuhlman T, Phillips R (2005): Transcriptional regulation by the numbers: models. Curr Opin Genet Dev, 15:125-135.
  2. Bintu L, Buchler NE, Garcia HG, Gerland U, Hwa T, Kondev J,Kuhlman T, Phillips R (2005): Transcriptional regulation by the numbers: application. Curr Opin Genet Dev (2), 15:125-135.
  3. Mitchel Lewis (2005) The Lac repressor. C. R. Biologies 328 (2005) 521–548.
  4. Subhayu Basu1, Yoram Gerchman, Cynthia H. Collins3, Frances H. Arnold& Ron Weiss,2: A synthetic multicellular system for programmed pattern formation. Nature 434 ,28 April 2005, 1130-113.
  5. Falcon C.M and Matthews K.S. (2000) Operator DNA sequence Variation Enhances High Affinity Binding by Hinge Helix Mutants of Lactose Repressor Protein. Biochemistry. 39, 11074-11084.
  6. Christian Berens , Dirk Schnappinger and Wolfgang Hillen ; The Role of the Variable Region in Tet Repressor for Inducibility by Tetracycline, Biological Chemisttry Volume 272, Number 11, Issue of March 14, 1997 pp. 6936-6942.
  7. ME346A Introduction to Statistical Mechanics – Wei Cai – Stanford University – Win 2011, Homepage; micro.stanford.edu/.
  8. Lee T.I, Young R.A: Transcription of eukaryotic protein-coding genes. Annu Rev Genet 2000, 34:77-137. p.109.
  9. Milo et al. Nucl. Acids Res. (2010) 38: D750-D753.
  10. Timothy S. Gardner, Charles R. Cantor1 & James J. Collins; Construction of a genetic toggle switch inEscherichia coli, 20 January 2000, Nature 403, 339-342.
  11. iGEM, Aberdeen Schotland 2009, homepage: igem.org.

  12. System Cost

  13. Bennett, R. A., & Lynch, J. M. (1981). Bacterial growth and development in the rhizosphere of gnotobiotic cereal plants. Journal of General Microbiology, 125(1), 95-102.
  14. Molina, L., Ramos, C., Duque, E., Ronchel, M. C., Garcı́a, J. M., Wyke, L., & Ramos, J. L. (2000). Survival of< i> Pseudomonas putida KT2440 in soil and in the rhizosphere of plants under greenhouse and environmental conditions. Soil Biology and Biochemistry, 32(3), 315-321.
  15. Puchałka J, Oberhardt MA, Godinho M, Bielecka A, Regenhardt D, et al. (2008) Genome-Scale Reconstruction and Analysis of the a KT2440 Metabolic Network Facilitates Applications in Biotechnology. PLoS Comput Biol 4(10): e1000210. doi: 10.1371/journal.pcbi.1000210.
  16. Oberhardt, M. A., Puchałka, J., Dos Santos, V. A. M., & Papin, J. A. (2011). Reconciliation of genome-scale metabolic reconstructions for comparative systems analysis. PLoS computational biology, 7(3), e1001116.
  17. Nogales, J., Palsson, B. Ø., & Thiele, I. (2008). A genome-scale metabolic reconstruction of Pseudomonas putida KT2440: iJN746 as a cell factory. BMC systems biology, 2(1), 79.
  18. van Duuren, J. B., Pucha, J., Mars, A. E., Bücker, R., Eggink, G., Wittmann, C., & dos Santos, V. A. (2013). Reconciling in vivo and in silico key biological parameters of Pseudomonas putida KT2440 during growth on glucose under carbon-limited condition. BMC biotechnology, 13(1), 93.
  19. The UniProt Consortium,Activities at the Universal Protein Resource (UniProt),[Pyocin,Q88ID3].
  20. Schnider-Keel, U., Seematter, A., Maurhofer, M., Blumer, C., Duffy, B., Gigot-Bonnefoy, C., ... & Keel, C. (2000). Autoinduction of 2, 4-diacetylphloroglucinol biosynthesis in the biocontrol agent Pseudomonas fluorescensCHA0 and repression by the bacterial metabolites salicylate and pyoluteorin. Journal of Bacteriology, 182(5), 1215-1225.
  21. Arfi, K., Landaud, S., & Bonnarme, P. (2006). Evidence for distinct L-methionine catabolic pathways in the yeast Geotrichum candidum and the bacterium Brevibacterium linens. Applied and environmental microbiology, 72(3), 2155-2162.
  22. Salah El Din, A. L. M., Kyslík, P., Stephan, D., & Abdallah, M. A. (1997). Bacterial iron transport: Structure elucidation by FAB-MS and by 2D NMR (H 1,C 13,N 15) of pyoverdin G4R, a peptidic siderophore produced by a nitrogen-fixing strain of Pseudomonas putida. Tetrahedron, 53(37), 12539-12552.
  23. Boukhalfa, H., Reilly, S. D., Michalczyk, R., Iyer, S., & Neu, M. P. (2006). Iron (III) coordination properties of a pyoverdin siderophore produced by Pseudomonas putida ATCC 33015. Inorganic chemistry, 45(14), 5607-5616.
  24. Matlab
  25. Cobra Toolbox
  26. Gurobi
  27. Buxton, E. W. (1962). Root exudates from banana and their relationship to strains of the Fusarium causing Panama wilt. Annals of Applied Biology, 50(2), 269-282.
  28. Bais, H. P., Weir, T. L., Perry, L. G., Gilroy, S., & Vivanco, J. M. (2006). The role of root exudates in rhizosphere interactions with plants and other organisms. Annu. Rev. Plant Biol., 57, 233-266.
  29. Armstrong, W., & Drew, M. C. (2002). Root growth and metabolism under oxygen deficiency. Plant roots: the hidden half, 3, 729-761.
  30. Oberhardt, M.A., J. Puchalka, V.A.P. Martins dos Santos, and J.A. Papin. 2011. Reconciliation of genome-scale metabolic reconstructions for comparative systems analysis. PLoS Computational Biology, 7(3).
  31. BNID 107924, Milo et al 2010.
  32. BNID 100985, Milo et al 2010.
  33. Gross, H., & Loper, J. E. (2009). Genomics of secondary metabolite production by Pseudomonas spp. Natural product reports, 26(11), 1408-1446.
  34. Stoker, N. G., Fairweathe, N. F., & Spratt, B. G. (1982). Versatile low-copy-number plasmid vectors for cloning in Escherichia coli. Gene, 18(3), 335-341.

  35. System Design

  36. Shu C-C, Chatterjee A, Dunny G, Hu W-S, Ramkrishna D (2011) Bistability versus Bimodal Distributions in Gene Regulatory Processes from Population Balance. PLoS Comput Biol 7(8): e1002140.
  37. Volfson, D. & Marciniak, J. et al.(2005) Origins of extrinsic variability in eukaryotic gene expression. Nature 21 December 2005
  38. Martinez Arias, A. & Hayward (2006), P. Filtering transcriptional noise during development: concepts and mechanisms. Nature Rev. Genet. 7, 34–44
  39. Pallavi Penumetcha, Kin Lau, Xiao Zhu, Kelly Davis, Todd T. Eckdahl, A. Malcolm Campbell. (2010) Improving the Lac System for Synthetic Biology. BIOS. 7-15
  40. R.A. Fasani, M.A. Savageau Molecular mechanisms of multiple toxin-antitoxin systems are coordinated to govern the persister phenotype (2013) . Proc. Natl. Acad. Sci. USA, 110, pp. E2528–E2537
  41. Etienne Maisonneuve, Kenn Gerdes (2014), Molecular Mechanisms Underlying Bacterial Persisters. Cell Volume 157, Issue 3, Pages 539–548
  42. Rué, P., & Garcia-Ojalvo, J. (2013). Modeling gene expression in time and space. Annual review of biophysics, 42, 605-627.
  43. Swain, P. S. (2004). Efficient attenuation of stochasticity in gene expression through post-transcriptional control. Journal of molecular biology, 344(4), 965-976
  44. Gonze, D. (2013). Modeling the effect of cell division on genetic oscillators. Journal of theoretical biology, 325, 22-33.
  45. Braun, D., Basu, S., & Weiss, R. (2005, March). Parameter estimation for two synthetic gene networks: a case study. In Acoustics, Speech, and Signal Processing, 2005. Proceedings.(ICASSP'05). IEEE International Conference on (Vol. 5, pp. v-769). IEEE.
  46. Purcell, O., Grierson, C. S., Bernardo, M., & Savery, N. J. (2012). Temperature dependence of ssrA-tag mediated protein degradation. J. Biol. Eng, 6(10).
  47. Alon, U. (2006). An introduction to systems biology: design principles of biological circuits. CRC press. Nevozhay, D., Adams, R. M., Murphy, K. F., Josić, K., & Balázsi, G. (2009). Negative autoregulation linearizes the dose–response and suppresses the heterogeneity of gene expression. Proceedings of the National Academy of Sciences, 106(13), 5123-5128.
  48. Braun, D., Basu, S., & Weiss, R. (2005, March). Parameter estimation for two synthetic gene networks: a case study. In Acoustics, Speech, and Signal Processing, 2005. Proceedings.(ICASSP'05). IEEE International Conference on (Vol. 5, pp. v-769). IEEE.