Team:Valencia Biocampus/Modeling
From 2014.igem.org
(154 intermediate revisions not shown) | |||
Line 17: | Line 17: | ||
<li><a href="#OrthogonalityTab" data-toggle="tab">Orthogonality</a></li> | <li><a href="#OrthogonalityTab" data-toggle="tab">Orthogonality</a></li> | ||
<li><a href="#OpenLicenseTab" data-toggle="tab">Open License</a></li> | <li><a href="#OpenLicenseTab" data-toggle="tab">Open License</a></li> | ||
- | <li><a href="#OtherTab" data-toggle="tab">Other</a></li> | + | <li><a href="#OtherTab" data-toggle="tab">Other projects</a></li> |
<li><a href="#ReferencesTab" data-toggle="tab">References</a></li> | <li><a href="#ReferencesTab" data-toggle="tab">References</a></li> | ||
</ul> | </ul> | ||
Line 26: | Line 26: | ||
<h2>Presentation</h2> | <h2>Presentation</h2> | ||
<p> | <p> | ||
- | + | In this section, we present our models to describe the four legs of the ST$^2$OOL. We have closely followed all the experimental activities, so our models do not pretend to | |
- | + | completely characterize a single dynamical process of a cell, but instead we intend to capture the dependences of our systems on the variables explored | |
+ | by our experiments. | ||
</p> | </p> | ||
<p> | <p> | ||
- | + | <a href="Modeling#StabilityTab"><strong>Stability</strong></a> We have modeled the variation of the activity of cells transformed with Biobricks, when subjected to different temperatures or pH, by considering changes in the Biobrick transcription | |
- | + | and translation rates. | |
- | + | </p> | |
+ | <p> | ||
+ | <a href="Modeling#StandardizationTab"><strong>Standardization</strong></a> We have modeled the expression in different strains by comparing the transcription rates. We have built an app that computes the standardization-stability index and generates the promoter sequence for a given expression rate. | ||
+ | </p> | ||
+ | <p> | ||
+ | <a href="Modeling#OrthogonalityTab"><strong>Orthogonality</strong></a> | ||
+ | We have studied deterministic and stochastic models of the plasmid growth in order to learn about the unspecific interaction between two plasmids in a cell. | ||
+ | </p> | ||
+ | <p> | ||
+ | <b><a href="Modeling#OpenLicenseTab"><strong>Open License</strong></a></b> We have modeled the susceptibility to patent a human creation and define the factors that influence it . | ||
</p> | </p> | ||
- | |||
</div> | </div> | ||
Line 44: | Line 53: | ||
<h2>Stability</h2> | <h2>Stability</h2> | ||
<p> | <p> | ||
- | Living organisms have genetic and biochemical tools to respond to stress. These tools normally up-regulate some stress-specific genes, and, directly or | + | Living organisms have genetic and biochemical tools that allow them to respond to stress. These tools normally up-regulate some stress-specific genes, and, directly or |
- | indirectly, down-regulate the rest of normal genes (most of them housekeeping genes). On the other hand, every | + | indirectly, down-regulate the rest of normal genes (most of them housekeeping genes). On the other hand, every <em>Escherichia coli</em> component works worse under non-optimal conditions |
- | + | (for example, nutrients uptake or ATP production), so that, every <em>E. coli</em> function will be reduced. | |
- | We model the relative expression of biobricks | + | We model the relative expression of biobricks when subjected to variations of temperature, pH and salt concentration. |
</p> | </p> | ||
</html> | </html> | ||
Line 87: | Line 96: | ||
<p> | <p> | ||
The expression level of our heterologous gene can be calculated by equations 1 and 2. Our main goal is to find which of the terms will change most | The expression level of our heterologous gene can be calculated by equations 1 and 2. Our main goal is to find which of the terms will change most | ||
- | significantly with the different stress factors. We | + | significantly with the different stress factors. We are interested in the activity in comparison to a reference one, which we take as the activity of <em>E. coli</em> bacteria grown on the LB medium at <u>38.7 <sup>o</sup>C, pH 7, 1% of NaCl, no UV irradiation and 1 atm of pressure.</u> |
- | + | ||
</p> | </p> | ||
<p> | <p> | ||
Line 99: | Line 107: | ||
</p> | </p> | ||
</html> | </html> | ||
+ | |||
=== Transcription rate === | === Transcription rate === | ||
<html> | <html> | ||
<p> | <p> | ||
- | Among several factors, we first discuss the transcription rate. In the range of temperatures (or pH) that we have studied, the polymerase is | + | Among several factors, we will first discuss the transcription rate. In the range of temperatures (or pH) that we have studied, the polymerase (RNAP) is |
- | sufficiently stable to assume that the elongation rate is not affected and the interaction energies do not change significantly (for example, the promoter | + | sufficiently stable to assume that the elongation rate is not affected and the interaction energies RNAP-promoter do not change significantly (for example, the promoter |
- | specific interaction | + | specific interaction only changes 16%, from 33 <sup>o</sup>C to 44 <sup>o</sup>C). We did not find any significant difference in the nucleotide |
- | compositions of E. coli under stress [2]. Therefore, main variations in the transcription rate are due to the probability of RNAP-promoter bounds | + | compositions of <em>E. coli</em> under stress [2]. Therefore, main variations in the transcription rate are due to the probability of RNAP-promoter bounds |
$$ | $$ | ||
\begin{equation} | \begin{equation} | ||
Line 112: | Line 121: | ||
\end{equation} | \end{equation} | ||
$$ | $$ | ||
- | On the other hand, | + | On the other hand, down-regulation affects the expression level of our Biobrick, because it works as a housekeeping gene (with a normal constitutive promoter for the RNAP). For example, <em>E. coli</em> has a common σ factor for housekeeping genes (σ70), and, another one (σ32) for temperature |
- | + | ||
stress. In the case of temperature stress, the bacteria will synthetize σ32, which will bind RNA polymerase (RNAP) cores that won’t be free for σ70, and so | stress. In the case of temperature stress, the bacteria will synthetize σ32, which will bind RNA polymerase (RNAP) cores that won’t be free for σ70, and so | ||
they will not be able to synthetize mRNA of the normal housekeeping genes. | they will not be able to synthetize mRNA of the normal housekeeping genes. | ||
Line 123: | Line 131: | ||
<p> | <p> | ||
The RNAP is bound to a promoter (unspecific DNA) with $\epsilon_{S}$ ($\epsilon_{NS}$) free energy. There are | The RNAP is bound to a promoter (unspecific DNA) with $\epsilon_{S}$ ($\epsilon_{NS}$) free energy. There are | ||
- | $N_{NS}$ nonspecific sites on the E. coli genome (sequences of 41bp) and R free E$\sigma$70, with $R \ll N_{NS}$ . | + | $N_{NS}$ nonspecific sites on the <em>E. coli</em> genome (sequences of 41bp) and R free E$\sigma$70, with $R \ll N_{NS}$ . |
The combinatory number of unbound states, with energy R$\epsilon_{NS}$, is | The combinatory number of unbound states, with energy R$\epsilon_{NS}$, is | ||
$$ | $$ | ||
Line 168: | Line 176: | ||
</p> | </p> | ||
</html> | </html> | ||
+ | |||
=== Synthesis rate === | === Synthesis rate === | ||
<html> | <html> | ||
Line 173: | Line 182: | ||
Modeling synthesis rate follows the lines discussed in transcription rate, but the working regime is quite different. | Modeling synthesis rate follows the lines discussed in transcription rate, but the working regime is quite different. | ||
The energies of the mRNA-ribosome sequences interactions are much larger than in the RNAP-promoter ones. | The energies of the mRNA-ribosome sequences interactions are much larger than in the RNAP-promoter ones. | ||
- | The average difference of free energies (computed with the RBS calculator for random sequences of E. coli genes) is | + | The average difference of free energies (computed with the RBS calculator for random sequences of <em>E. coli</em> genes) is |
- | $\Delta\epsilon = -20.79 kT$, which | + | $\Delta\epsilon = -20.79 kT$, which means that, in the equation: |
+ | $$ | ||
+ | \begin{equation} | ||
+ | \frac{r_T}{r^0_T} = \frac{Prob_B}{Prob^0_B} = \frac { 1+\frac{R^0}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }} } { \frac {R^0}{R}+\frac{R^0}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }}} , | ||
+ | \end{equation} | ||
+ | $$ | ||
+ | The term $\frac{R^0}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }}>>>\frac {R^0}{R}$. With this, and knowing that there is no effect on the ribosomal protein composition with temperature changes [5], we conclude that the most signifficant factor should be the elongation rate. | ||
</p> | </p> | ||
Therefore, | Therefore, | ||
Line 183: | Line 198: | ||
$$ | $$ | ||
where $r_E$ ($r_I$) is the elongation (initiation) rate and the superscript denotes our reference case. | where $r_E$ ($r_I$) is the elongation (initiation) rate and the superscript denotes our reference case. | ||
- | In E. coli, most of the ribosomes are translating (4 active every one free) [6], thus $\frac{r^0_E}{r^0_I} \simeq\frac{1}{4}$ and | + | In <em>E. coli</em>, most of the ribosomes are translating (4 active every one free) [6], thus $\frac{r^0_E}{r^0_I} \simeq\frac{1}{4}$ and |
$$ | $$ | ||
\begin{equation} | \begin{equation} | ||
Line 191: | Line 206: | ||
</p> | </p> | ||
<p> | <p> | ||
- | We discuss next the dependence of elongation rates with temperature. We use the relation of amino acid concentration with | + | We discuss next the dependence of elongation rates with temperature. We use the relation of amino acid concentration with |
- | temperature [7]. Using data from [8], the | + | temperature [7]. Using data from [8], we know that small variations on the ternary complex does not affects the elongation rate. The slowests steps are: |
- | 1) | + | 1) GTP hydrolysis on ternary complex (0.01 s), 2) dipeptide bond formation (0.03s) and 3) translocation (0.022s). |
$$ | $$ | ||
\begin{equation} | \begin{equation} | ||
Line 199: | Line 214: | ||
\end{equation} | \end{equation} | ||
$$ | $$ | ||
- | where $E^i_a$ is (32, 75.2, | + | where $E^i_a$ is (32, 75.2, 29.260) kJ/mol and $A^{-1}_i$ is (2.5 $\cdot 10^7$, 1.6 $\cdot 10^{14}$, 3.9 $\cdot 10^6$). |
</p> | </p> | ||
</html> | </html> | ||
+ | |||
=== Summary === | === Summary === | ||
<html> | <html> | ||
Line 215: | Line 231: | ||
</html> | </html> | ||
- | === pH | + | === 2) pH === |
<html> | <html> | ||
<p> | <p> | ||
- | + | Contrary to the temperature situations, internal pH hardly varies in <em>E. coli</em> when the external pH does. Therefore, we focus on the intrinsic <em>E. coli</em> element that may be affected by this type of stress, subsequently influencing Biobrick expression. The main factor is the number of free RNAP, which will be reduced due to the synthesis of stress response genes. We estimate the ratio of free RNAP under different conditions by the ratio of their growth rates | |
- | We estimate the ratio of free RNAP under different conditions by the ratio of their growth rates | + | |
$$ | $$ | ||
\begin{equation} | \begin{equation} | ||
Line 233: | Line 248: | ||
</html> | </html> | ||
+ | |||
=== Comparison with wet-lab data === | === Comparison with wet-lab data === | ||
<html> | <html> | ||
<p> | <p> | ||
- | Our model captures the | + | <b><a href="Results#StabilitySection">Our model captures the shape of the observed temperature profile and predicts that additional variations at lower temperatures are independent of the chosen Biobrick</a></b>. We have found that transcription rate decreases with temperature as a consequence of the reduction of free RNAP; and synthesis rate grows above the optimal temperature. Their combination leads to the observed variation of fluorescence with temperature. Unfortunately, this is not the whole story. Our model fails to fit the amount of variation, which is larger experimentally than the predicted one. By comparing our theoretical predictions and experimental results of XL1 Blue strain transformed with Bb1 or Bb3, we find that the ratio between theoretical prediction and experimental results is the same at under-optimal temperatures (in fact, accidentally too close given the error bars). This fact suggests that the temperature dependence not included in our model should come from cell mechanisms; independent from Biobrick parameters (promoter strength, ...). On the contrary, the comparison of theoretical predictions and experimental results at higher temperatures show that our model might require more dynamics which depends on Biobrick parameters. |
- | + | </p> | |
- | fails to fit the amount of variation, which is | + | <p> |
- | + | <b><a href="Results#StabilitySection">Our model captures the shape of the observed pH profile</a></b>. We have found that the lowest fluorescence happens at the pH with lowest growth, while it is mildly dependent on the pH from 6 to 8. | |
- | temperatures | + | |
- | our model should come from cell mechanisms independent | + | |
- | theoretical predictions and experimental results at | + | |
- | + | ||
</p> | </p> | ||
Line 250: | Line 262: | ||
<div class="tab-pane fade" id="StandardizationTab"> | <div class="tab-pane fade" id="StandardizationTab"> | ||
<h2>Standardization</h2> | <h2>Standardization</h2> | ||
- | |||
- | |||
- | |||
<p> | <p> | ||
- | We | + | We will now discuss Biobrick expression in different <em>E. coli</em> strains. The effect of the changes in the number of free RNAP on mRNA synthesis rate depend on the energy of the promoter. The relative transcription rate in different strains can be written as |
$$ | $$ | ||
\begin{equation} | \begin{equation} | ||
Line 261: | Line 270: | ||
\end{equation} | \end{equation} | ||
$$ | $$ | ||
- | where we can see that the dominance of the second term over the first one, controlled by the promoter strength, would lead to similar rate independently on the particular strain chosen. Comparing this formula with the relative transcription rate in stability, we find that promoters whose binding energy makes them more standard also implies more stable promoters. | + | where we can see that the dominance of the second term over the first one, controlled by the promoter strength, would lead to similar rate independently on the particular strain chosen. Comparing this formula with the relative transcription rate in stability, we find that promoters whose binding energy makes them more standard also implies more stable promoters. Strong promoter sequences experience a decrease in synthesis rates when subjected to an increase of promoter energy. We use these results and the energy matrix derived in [3] to create an application that is able to compute promoter sequences for a given relative (to the maximum) expression rate and associated standard and stable rate. The physical parameters behind are the promoter strength (using the data of Anderson’ promoters) and the binding energy and the tool create a sequence of a constitutive promoter with the appropriated trade-off between strength and energy, therefore, with the appropriated trade-off between stability & standardization (ST$^2$)and rate of expression. We define the ST$^2$ index as $$ F_{st2} = (F - 0.012)/0.51 $$ with $$F = \frac{R'}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }}, $$ where $F_{sty}$ has been scaled to run from 1 to 0. We have used a fitting formula for the rate of expression relative to the maximum rate [3] as a function of the binding energy $\Delta\epsilon$, $$ r_T = 2.29 + 0.27 \frac{\Delta\epsilon}{kT}. $$ With this approach, we can infer the ST$^2$ level and the promoter sequence for a given expression rate. We show below our app in javascript. |
- | <p>Move the cursor to a given expression rate (fraction of | + | <div class="col-sm-9 center-block" style="padding:20px;"> |
+ | <p><b>Generate a promoter sequence and compute the ST$^2$ (Stability & Standardization) index for a given expression rate</b></p> | ||
+ | <p>Move the cursor to a given expression rate (as a fraction of its maximum) and press Generate a chain</p> | ||
<p style="text-align:center"> | <p style="text-align:center"> | ||
<input id="sliderApp" type="text" data-slider-min="0" data-slider-max="100" data-slider-step="1" data-slider-value="25" /> | <input id="sliderApp" type="text" data-slider-min="0" data-slider-max="100" data-slider-step="1" data-slider-value="25" /> | ||
<span class="label label-success" id="percent1">25%</span> | <span class="label label-success" id="percent1">25%</span> | ||
- | |||
<input type="hidden" name="textbox1" id="textbox1" value="25"/> | <input type="hidden" name="textbox1" id="textbox1" value="25"/> | ||
<br/><br/> | <br/><br/> | ||
- | < | + | <button type="button" name="button" id="button1" data-loading-text="Generating..." class="btn btn-default">Generate a chain</button> |
</p> | </p> | ||
<br/> | <br/> | ||
- | <p>Promoter sequence (from -41 to -1 nt) is:</p> <p id="output"></p> | + | <p class="text-centered">Promoter sequence (from -41 to -1 nt) is:</p> |
- | <script type="text/javascript" src="https://2014.igem.org/wiki/index.php?title=Template:Team:Valencia_Biocampus/Templates/js/ | + | <p id="output"></p> |
+ | <p class="text-centered">ST$^2$ index: | ||
+ | <span class="label label-danger" id="percent2"></span> | ||
+ | </p> | ||
+ | <script type="text/javascript" src="https://2014.igem.org/wiki/index.php?title=Template:Team:Valencia_Biocampus/Templates/js/NewModelingApp&action=raw&ctype=text/javascript"></script> | ||
<script type="text/javascript"> | <script type="text/javascript"> | ||
- | |||
$("#sliderApp").slider({ | $("#sliderApp").slider({ | ||
formatter: function(value) { | formatter: function(value) { | ||
Line 289: | Line 302: | ||
$("#textbox1").val(slideEvt.value); | $("#textbox1").val(slideEvt.value); | ||
$("#percent1").text(slideEvt.value+'%'); | $("#percent1").text(slideEvt.value+'%'); | ||
- | |||
}); | }); | ||
+ | $('#button1').click(function () { | ||
+ | $("#input").empty(); | ||
+ | $(this).prop("disabled",true); | ||
+ | $(this).text("GENERATING ..."); | ||
+ | |||
+ | CE(); | ||
+ | }); | ||
</script> | </script> | ||
</div> | </div> | ||
</html> | </html> | ||
+ | |||
=== Comparison with wet-lab data === | === Comparison with wet-lab data === | ||
<html> | <html> | ||
<p> | <p> | ||
- | Our wet-lab results show that | + | <b><a href="Results#StandardizationSection">Our wet-lab results show that Biobrick 5 is the most standard one.</a></b>Our model predicts that systems with less modules or factors will be more stable and standard. Biobrick 5 has a tetracycline depending promoter, whose limitating step is the binding of tetracycline to TetR. In this scenario, the intracellular concentration of aTc and the binding to the TetR repressor do not depend on the strain, so there is only one factor in which the strain could have an important effect: the TetR intracellular concentration, which has a constitutive expression. Analyzing the other Biobricks, we found that the strains with less expression of the constitutive promoter are those ones with higher expression when we use the Ptet promoter (but JM109 and BB1), which agree with our prediction. |
Line 313: | Line 333: | ||
<p> | <p> | ||
+ | Plasmid dynamics have been studied exhaustively over these years using mathematical models: from how plasmid growth is regulated to plasmid partitioning after cell division. However, little is known about how two identical plasmids interact with each other. This part of our modelling project aims to shed light upon this issue using previous models of ColE1 plasmids and the properties of two different plasmids with same origin of replication when these are introduced in a cell. We employed both deterministic and stochastic models to try to answer whether two plasmids behave independently, and developed some approximations to highlight important features of the system. | ||
+ | |||
<h3> | <h3> | ||
- | + | Deterministic model | |
- | </h3> | + | </h3> |
- | + | ||
- | + | ||
<br> | <br> | ||
- | + | <p> | |
- | + | This model is taken from [16]. Here, plasmid levels are regulated both by an RNA molecule that inhibits plasmid replication and by the dilution of plasmids as a result of cell growth and division: | |
- | + | ||
- | + | ||
- | + | ||
- | This model is taken from | + | |
- | $$ \frac{d P}{dt} = r \cdot P \cdot R(s) - \frac{\log 2}{\tau} \cdot P $$ | + | $$ \frac{d P}{dt} = r \cdot P \cdot R(s) - \frac{\log 2}{\tau} \cdot P, $$ |
- | $$ \frac{d s}{dt} = \alpha \cdot P - \beta \cdot s $$ | + | $$ \frac{d s}{dt} = \alpha \cdot P - \beta \cdot s, $$ |
- | where | + | where $P$ is the number of plasmids, $s$ is the number of RNA inhibitors, $r$ is the plasmid growth rate, $\tau$ is the cell growth rate, $\alpha$ is the RNA inhibitor synthesis rate |
- | + | , $\beta$ is the RNA inhibitor decay rate and $R(s)$ is the repression effect of RNA inhibitor with | |
- | + | $$ R(s) = \frac{1}{(1+\frac{s}{s_0})^n},$$ | |
- | + | where $s_0$ is the value of RNA inhibitor necessary to repress half of the replication events and $n$ is the cooperativity. Plasmid levels are regulated both by a negative feedback and by cell growth. If we consider several plasmids, we get<br><br> | |
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | $$ | + | $$ \frac{d P_i}{dt} = r_i \cdot P_i \cdot R(s) - \frac{\log 2}{\tau} \cdot P_i,$$ |
- | + | $$ \frac{d s}{dt} = \alpha \cdot \sum \limits_{i=1}^n{P_i} - \beta \cdot s,$$ | |
+ | with a steady state solution<br><br> | ||
- | + | $$ s_{ss} = R^{-1}(\frac{\log 2}{\tau \cdot r_i}),$$ | |
- | + | ||
- | $$ \ | + | $$ \sum \limits_{i=1}^n{P_{i,ss}} = (\frac{ \beta \cdot s_{ss}}{\alpha}).$$ |
- | $$ | + | <br>In a simplified way, with two identical plasmids, we get: |
+ | $$ P_{1,ss} + P_{2,ss} = K.$$ | ||
- | |||
- | |||
- | + | It is noticeable that plasmid levels are not independent. Small perturbations can increase the relative level of one plasmid, favoring its growth. As they are repressed as a function of the total number of plasmids but they grow depending on their number, the system is prone to introduce asymmetries | |
- | $$ \ | + | $$ \frac{d P_1}{d P_2} = \frac{r_1 \cdot P_1}{r_2 \cdot P_2} .$$ |
- | |||
- | |||
- | + | If growth rates are constant, the slope is determined by the initial conditions | |
+ | $$ m = \frac{P_{01}}{P_{02}} .$$ | ||
- | + | <h4>Quasi-steady state approach</h4> | |
+ | To simplify our equations, we assume inhibitor dynamics are faster than plasmid growth, so RNA levels are a direct representation of plasmid levels | ||
- | $$ \frac{d | + | $$ \frac{d P_i}{dt} = r_i \cdot P_i \cdot R(P_i) - \frac{\log 2}{\tau} \cdot P_i .$$ |
- | < | + | <h4>Dependence on initial conditions:</h4> |
- | + | To test our hypothesis, we integrate our differential equations using numerical methods with several initial conditions. As expected, small differences in initial proportions are amplified by the dynamics. | |
- | < | + | |
- | + | ||
- | + | ||
- | + | <div class="thumbnail col-sm-11 center-block"> | |
- | + | <a href="https://static.igem.org/mediawiki/2014/d/d6/VB2014_Det_FinalDistribution.jpg" rel="lightbox"> | |
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | <div | + | |
- | <a href="https://static.igem.org/mediawiki/2014/d/d6/VB2014_Det_FinalDistribution.jpg"> | + | |
<img src="https://static.igem.org/mediawiki/2014/d/d6/VB2014_Det_FinalDistribution.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | <img src="https://static.igem.org/mediawiki/2014/d/d6/VB2014_Det_FinalDistribution.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | ||
</a> | </a> | ||
- | <a href="https://static.igem.org/mediawiki/2014/7/73/VB2014_Det_PhaseMap.jpg"> | + | <a href="https://static.igem.org/mediawiki/2014/7/73/VB2014_Det_PhaseMap.jpg" rel="lightbox"> |
<img src="https://static.igem.org/mediawiki/2014/7/73/VB2014_Det_PhaseMap.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | <img src="https://static.igem.org/mediawiki/2014/7/73/VB2014_Det_PhaseMap.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | ||
</a> | </a> | ||
+ | <div class="caption"><strong>Figure 1. Dependence on initial conditions</strong> A) Steady state distribution of plasmids from different initial conditions. As seen, small differences in initial conditions are amplified. B) Phase space map of plasmid growth. | ||
+ | </div> | ||
</div> | </div> | ||
<br> | <br> | ||
- | |||
- | |||
- | |||
<h3> | <h3> | ||
Stochastic model | Stochastic model | ||
- | </h3 | + | </h3> |
- | + | ||
<br> | <br> | ||
- | First of all, we | + | The deterministic model assumes negligible noise. However, this might not be the case and small variations in plasmid levels could be highly amplified. Hence, we now discuss the importance of noise in plasmid growth. We used the <b>Stochastic Simulation Algorithm (SSA)</b> developed by Gillespie (see [13]). In short, SSA simulates by iteration which processes occur at any time using a random number generator.<br> |
- | + | ||
- | <li> | + | First of all, we used the algorithm to develop a code assuming growth of the plasmids and cell division as independent processes. The steps in the algorithm are:<br> |
- | $$ \alpha_{P_i \rightarrow P_i+1} = | + | |
- | + | <li>Compute the propensity function of each reaction and the overall propensity function: | |
- | $$ | + | $$ \alpha_{P_i \rightarrow P_i+1} = r_i \cdot P_i \cdot R(P_i), $$ |
- | where $\theta_1$ is a random number from a uniform distribution U(0,1). <br><br> | + | $$ \alpha_{0} = \sum \limits_{i=1}^{n}\alpha_{P_i \rightarrow P_{i+1}}. $$ |
- | <li> | + | <li>Compute when the next chemical reaction takes place : |
- | $$ \sum \limits_{i=1}^{j}{\alpha_{P_i \rightarrow P_i+1}} < \theta_2 \leq \sum \limits_{i=1}^{j | + | $$ t = \frac{1}{\alpha_{0}} \cdot \log(\frac{1}{\theta_1}),$$ |
+ | where $\theta_1$ is a random number from a uniform distribution U(0,1) we make use of the random distribution for exponential reactions.<br><br> | ||
+ | <li>Compute which reaction happens: | ||
+ | $$ \frac{1}{\alpha_{0}} \sum \limits_{i=1}^{j-1}{\alpha_{P_i \rightarrow P_i+1}} < \theta_2 \leq \frac{1}{\alpha_{0}} \sum \limits_{i=1}^{j}{\alpha_{P_i \rightarrow P_i+1}},$$ | ||
+ | where $\theta_2$ is a random number given by the uniform distribution U(0,1). | ||
<br> | <br> | ||
<li>Plasmid partition after cell division at t = $\tau$ <br> | <li>Plasmid partition after cell division at t = $\tau$ <br> | ||
- | + | We show below the evolution of a two plasmids system for two different initial conditions: $P_1=P_2=10$ and $P_1 = 10$ and $P_2 = 50$. <br><br> | |
- | + | ||
- | <div | + | <div class="thumbnail col-sm-11 center-block"> |
- | <a href="https://static.igem.org/mediawiki/2014/6/64/VB2014_Sto_Example.jpg"> | + | <a href="https://static.igem.org/mediawiki/2014/6/64/VB2014_Sto_Example.jpg" rel="lightbox"> |
- | <img src="https://static.igem.org/mediawiki/2014/6/64/VB2014_Sto_Example.jpg" width=" | + | <img src="https://static.igem.org/mediawiki/2014/6/64/VB2014_Sto_Example.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> |
</a> | </a> | ||
- | <a href="https://static.igem.org/mediawiki/2014/c/cd/VB2014_Sto_Example2.jpg"> | + | <a href="https://static.igem.org/mediawiki/2014/c/cd/VB2014_Sto_Example2.jpg" rel="lightbox"> |
- | <img src="https://static.igem.org/mediawiki/2014/c/cd/VB2014_Sto_Example2.jpg" width=" | + | <img src="https://static.igem.org/mediawiki/2014/c/cd/VB2014_Sto_Example2.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> |
</a> | </a> | ||
+ | <div class="caption"><strong>Figure 2. Stochastic simulation (solid lines) versus deterministic solution (solid + dots lines) for different initial conditions (left $P_1=P_2=10$, right $P_1 = 10$ and $P_2 = 50$)</strong> This figure depicts the differences between these two modeling approaches. When several cells are taken into account, the dynamics are similar to the deterministic solution. However, the intracellular levels of plasmids can vary from the deterministic solution as shown here. Each plasmid is represented by a color and the total level of plasmids of the deterministic solution is in black. | ||
+ | </div> | ||
</div> | </div> | ||
<br> | <br> | ||
- | As expected, the deterministic model gives a mean version of the dynamics. | + | As expected, the deterministic model gives a mean version of the dynamics. Single cell dynamics shows oscillatory behaviour. The levels of each plasmid can vary but their sum remains unchaged (black line) |
- | |||
<h4> | <h4> | ||
A lineage model | A lineage model | ||
- | </h4> | + | </h4> |
- | Aiming to obtain a better description | + | |
- | As tracking bacteria and their plasmid levels is an important issue, we implemented a tree plot in our code using <a href="http://www.mathworks.co.uk/matlabcentral/fileexchange/35623-tree-data-structure-as-a-matlab-class"><b>a matlab package</b></a> | + | Aiming to obtain a better description, we studied plasmid levels in a lineage from a single cell. <br> |
+ | As tracking bacteria and their plasmid levels is an important issue, we implemented a tree plot in our code using <a href="http://www.mathworks.co.uk/matlabcentral/fileexchange/35623-tree-data-structure-as-a-matlab-class"><b>a matlab package</b></a>. | ||
Plasmids were given to each daughter cell following a gaussian distribution with a standard deviation of $\frac{1}{\sqrt{n}}$ | Plasmids were given to each daughter cell following a gaussian distribution with a standard deviation of $\frac{1}{\sqrt{n}}$ | ||
- | + | ||
- | <div | + | <div class="thumbnail col-sm-11 center-block"> |
- | <a href="https://static.igem.org/mediawiki/2014/6/65/VB2014_Lin_PlasmidExample.jpg"> | + | <a href="https://static.igem.org/mediawiki/2014/6/65/VB2014_Lin_PlasmidExample.jpg" rel="lightbox"> |
<img src="https://static.igem.org/mediawiki/2014/6/65/VB2014_Lin_PlasmidExample.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | <img src="https://static.igem.org/mediawiki/2014/6/65/VB2014_Lin_PlasmidExample.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | ||
</a> | </a> | ||
- | <a href="https://static.igem.org/mediawiki/2014/d/de/VB2014_Lin_TreeExample.jpg"> | + | <a href="https://static.igem.org/mediawiki/2014/d/de/VB2014_Lin_TreeExample.jpg" rel="lightbox"> |
<img src="https://static.igem.org/mediawiki/2014/d/de/VB2014_Lin_TreeExample.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | <img src="https://static.igem.org/mediawiki/2014/d/de/VB2014_Lin_TreeExample.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | ||
</a> | </a> | ||
+ | <div class="caption"><strong>Figure 3. Plasmid dynamics over a lineage (I)</strong> As cells divide, they can suffer from variations in their number of copies of a plasmid, influencing the levels of daughter cells. Left picture shows the evolution in the levels of each plasmid (in different colors) and right picture is a lineage tree to track each cell. | ||
+ | </div> | ||
</div> | </div> | ||
<br> | <br> | ||
- | Also we | + | Also we were able to track the plasmid levels of each bacteria at the time of their division. |
<br> | <br> | ||
- | <div | + | <div class="thumbnail col-sm-11 center-block"> |
- | <a href="https://static.igem.org/mediawiki/2014/4/4e/VB2014_Lin_TreeExampleP1.jpg"> | + | <a href="https://static.igem.org/mediawiki/2014/4/4e/VB2014_Lin_TreeExampleP1.jpg" rel="lightbox"> |
<img src="https://static.igem.org/mediawiki/2014/4/4e/VB2014_Lin_TreeExampleP1.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | <img src="https://static.igem.org/mediawiki/2014/4/4e/VB2014_Lin_TreeExampleP1.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | ||
</a> | </a> | ||
- | <a href="https://static.igem.org/mediawiki/2014/0/08/VB2014_Lin_TreeExampleP2.jpg"> | + | <a href="https://static.igem.org/mediawiki/2014/0/08/VB2014_Lin_TreeExampleP2.jpg" rel="lightbox"> |
<img src="https://static.igem.org/mediawiki/2014/0/08/VB2014_Lin_TreeExampleP2.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | <img src="https://static.igem.org/mediawiki/2014/0/08/VB2014_Lin_TreeExampleP2.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | ||
</a> | </a> | ||
+ | <div class="caption"><strong>Figure 4. Plasmid dynamics over a lineage (II)</strong> Variations of plasmid copy number in each lineage can be track using a lineage tree. Each tree shows the level of its associated plasmid. | ||
+ | </div> | ||
</div> | </div> | ||
<br> | <br> | ||
- | The distribution | + | The distribution at the beginning and at the end of the simulation is: |
- | + | <div class="thumbnail col-sm-11 center-block"> | |
- | <div | + | <a href="https://static.igem.org/mediawiki/2014/e/ea/VB2014_Sto_Dist0-1000.jpg" rel="lightbox"> |
- | <a href="https://static.igem.org/mediawiki/2014/e/ea/VB2014_Sto_Dist0-1000.jpg"> | + | |
<img src="https://static.igem.org/mediawiki/2014/e/ea/VB2014_Sto_Dist0-1000.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | <img src="https://static.igem.org/mediawiki/2014/e/ea/VB2014_Sto_Dist0-1000.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | ||
</a> | </a> | ||
- | <a href="https://static.igem.org/mediawiki/2014/6/64/VB2014_Sto_Dist2000-2500.jpg"> | + | <a href="https://static.igem.org/mediawiki/2014/6/64/VB2014_Sto_Dist2000-2500.jpg" rel="lightbox"> |
<img src="https://static.igem.org/mediawiki/2014/6/64/VB2014_Sto_Dist2000-2500.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | <img src="https://static.igem.org/mediawiki/2014/6/64/VB2014_Sto_Dist2000-2500.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | ||
</a> | </a> | ||
+ | <div class="caption"><strong>Figure 5. Heat maps of time evolution of plasmid distribution in a lineage</strong> Here we show the distributions of the plasmid levels in a lineage for the first 500 simulations (left picture) and the last 500 simulations (right picture), from a single cell with the initial conditions $P_1 = P_2 = 75$. Notice the | ||
+ | broadening of the distribution with the lineage growth. | ||
+ | </div> | ||
</div> | </div> | ||
<br> | <br> | ||
- | + | Comparing with a distribution of orthogonal plasmids, we find that while both are symmetric (respect to the $P_1 = P_2$ line), | |
- | <div | + | non orthogonal plasmids lead to broader heat maps. |
- | <a href="https://static.igem.org/mediawiki/2014/0/07/VB2014_NoDependence.jpg"> | + | |
+ | <div class="thumbnail col-sm-5 center-block"> | ||
+ | <a href="https://static.igem.org/mediawiki/2014/0/07/VB2014_NoDependence.jpg" rel="lightbox"> | ||
<img src="https://static.igem.org/mediawiki/2014/0/07/VB2014_NoDependence.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | <img src="https://static.igem.org/mediawiki/2014/0/07/VB2014_NoDependence.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | ||
</a> | </a> | ||
+ | <div class="caption"><strong>Figure 6. Plasmid distribution in a system with no influence between different plasmids</strong> When each plasmid represses itself, the distribution gets narrower and the time-dependent evolution is lost. | ||
+ | </div> | ||
</div> | </div> | ||
- | |||
<h4> | <h4> | ||
A step beyond: Resistance | A step beyond: Resistance | ||
- | </h4 | + | </h4> |
- | As a proof of concept, we modelled the bacterial time of division as a function of both plasmids, being minimal when both plasmid levels are similar. | + | As a proof of concept, we modelled the bacterial time of division as a function of both plasmids, being minimal when both plasmid levels are similar. For this simple regulation, we obtained that the average level of plasmids in a population is the same for each DNA molecule. |
<br><br> | <br><br> | ||
- | <div | + | <div class="thumbnail col-sm-11 center-block"> |
- | <a href="https://static.igem.org/mediawiki/2014/a/af/VB2014_Res_PlasmidExample.jpg"> | + | <a href="https://static.igem.org/mediawiki/2014/a/af/VB2014_Res_PlasmidExample.jpg" rel="lightbox"> |
<img src="https://static.igem.org/mediawiki/2014/a/af/VB2014_Res_PlasmidExample.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | <img src="https://static.igem.org/mediawiki/2014/a/af/VB2014_Res_PlasmidExample.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | ||
</a> | </a> | ||
- | <a href="https://static.igem.org/mediawiki/2014/b/b4/VB2014_Res_TreeExample.jpg"> | + | <a href="https://static.igem.org/mediawiki/2014/b/b4/VB2014_Res_TreeExample.jpg" rel="lightbox"> |
<img src="https://static.igem.org/mediawiki/2014/b/b4/VB2014_Res_TreeExample.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | <img src="https://static.igem.org/mediawiki/2014/b/b4/VB2014_Res_TreeExample.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | ||
</a> | </a> | ||
+ | <div class="caption"><strong>Figure 7. Plasmid dynamics over a lineage in presence of resistance and plasmid-dependent growth </strong> Introduction of a new feature to the system, a plasmid-dependent growth, cells shows loss of synchrony in growth, favouring cells keeping similar levels of both plasmids. Left picture shows the plasmid copy number over time and right picture shows the lineage tree as a chronogram. | ||
+ | </div> | ||
</div> | </div> | ||
<br> | <br> | ||
- | This way to tune the plasmid levels | + | This way to tune the plasmid levels is unable to make them equal, as depicted in the following picture. Also, there are variations in their growth rates. |
<br><br> | <br><br> | ||
- | <div | + | <div class="thumbnail col-sm-11 center-block"> |
- | <a href="https://static.igem.org/mediawiki/2014/9/91/VB2014_Res_Evol.jpg"> | + | <a href="https://static.igem.org/mediawiki/2014/9/91/VB2014_Res_Evol.jpg" rel="lightbox"> |
<img src="https://static.igem.org/mediawiki/2014/9/91/VB2014_Res_Evol.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | <img src="https://static.igem.org/mediawiki/2014/9/91/VB2014_Res_Evol.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | ||
</a> | </a> | ||
- | <a href="https://static.igem.org/mediawiki/2014/b/b5/VB2014_Res_HistogramTau.jpg"> | + | <a href="https://static.igem.org/mediawiki/2014/b/b5/VB2014_Res_HistogramTau.jpg" rel="lightbox"> |
<img src="https://static.igem.org/mediawiki/2014/b/b5/VB2014_Res_HistogramTau.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | <img src="https://static.igem.org/mediawiki/2014/b/b5/VB2014_Res_HistogramTau.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | ||
</a> | </a> | ||
+ | <div class="caption"> | ||
+ | <strong>Figure 8. Resistance-induced plasmid levels modulation (left) and doubling cell time (right)</strong> Resistance is able to equal the levels of both plasmids, even there was a difference in initial conditions, although there may be a small difference. This is achieved by decreasing the doubling time of cells with similar levels of each plasmid, obtaining a distribution in these times as shown in the right picture. | ||
+ | </div> | ||
</div> | </div> | ||
<br> | <br> | ||
- | |||
- | |||
<h4> | <h4> | ||
Approaches for less computational expensive | Approaches for less computational expensive | ||
- | </h4><br> | + | </h4> |
- | Although these simulations are powerful and helpful and yield informative results, they are computationally expensive. This could become an important pitfall when dealing with a huge number of cells or long times. We | + | <br> |
+ | Although these simulations are powerful and helpful and yield informative results, they are computationally expensive. This could become an important pitfall when dealing with a huge number of cells or long times. We tried to bypass this problem by using an approach if we are working with populations as explained in [16]. We added a degradation process for each plasmid following the equations: <br> | ||
$$ \alpha_{P_i \rightarrow P_i-1} = \frac{log(2)}{\tau} \cdot P_i$$ | $$ \alpha_{P_i \rightarrow P_i-1} = \frac{log(2)}{\tau} \cdot P_i$$ | ||
- | + | This approach can be used to study cell dynamics in a population. With the system of equations, we capture the global dynamics for a cell population. | |
+ | If there is no resistance, we find dynamical examples where a plasmid is lost. | ||
<br><br> | <br><br> | ||
- | <div | + | <div class="thumbnail col-sm-5 center-block"> |
- | <a href="https://static.igem.org/mediawiki/2014/4/4b/VB2014_Approach_Example.jpg"> | + | <a href="https://static.igem.org/mediawiki/2014/4/4b/VB2014_Approach_Example.jpg" rel="lightbox"> |
<img src="https://static.igem.org/mediawiki/2014/4/4b/VB2014_Approach_Example.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | <img src="https://static.igem.org/mediawiki/2014/4/4b/VB2014_Approach_Example.jpg" width="340" height="255" alt="Simple Linear Regression Model" class="thumbnail" style="display: inline; float: none;"/> | ||
</a> | </a> | ||
+ | <div class="caption"> | ||
+ | <strong>Figure 9. Plasmid levels in cell lineages without resistance modulation using a less computational approach.</strong> Simulation of all features are time-consuming and computation expensive. However, simpler approaches can lead to similar results. In this particular example, the drift of a plasmid is shown. | ||
+ | </div> | ||
</div> | </div> | ||
<br> | <br> | ||
- | |||
- | + | ||
- | < | + | </html> |
- | + | ||
- | < | + | === Comparison with wet-lab data === |
- | + | <html> | |
- | + | <p> | |
- | + | The interaction between two plasmids in the same cell is still an uncharted territory. A relevant example is shown in [15], where the authors obtained an emergent behaviour as result of an unexpected relation between the bacterial physiology and plasmid dynamics. Using models in single plasmid dynamics, we switch between deterministic and stochastic models to learn about the unspecific interaction between two plasmids in a cell. <a href="Results#OrthogonalitySection"><b>There is an important variability, which drives the system away from orthogonality. Plasmids with identical origin of replication are not orthogonal</b> <i>per se</i></a>. Although we can introduce different resistances to tune the system, it | |
- | + | is hard to stabilize the system due to the large variability. Plasmids with different mechanisms of replication might solve it but it is not guaranteed. Certainly, new plasmids must be designed to make engineered cells more reliable and predictable. | |
+ | </p> | ||
Line 539: | Line 556: | ||
history: false, | history: false, | ||
highlightOffset: -40 | highlightOffset: -40 | ||
- | }).data("toc- | + | }).data("toc-tocify2"); |
$(".optionName").popover({ trigger: "hover" }); | $(".optionName").popover({ trigger: "hover" }); | ||
Line 549: | Line 566: | ||
<div class="tab-pane fade" id="OpenLicenseTab"> | <div class="tab-pane fade" id="OpenLicenseTab"> | ||
<h2>Open License</h2> | <h2>Open License</h2> | ||
- | + | <p> | |
+ | We model the susceptibility to patent a human creation by | ||
+ | $$ | ||
+ | P.I. = M \cdot N \cdot A ^{log(\frac{A+I}{2})}, | ||
+ | $$ | ||
+ | where P.I. is the patentability index, M is the morality, which is a non-legal term accepted in the E.U., and N refers to novelty. Both factors are described by a discrete variable which can be valued 0 or 1 depending on the results of a questionnaire. A is the inventive step and I is the industrial application which can be valued from 1 to 10 depending of the | ||
+ | results of a test. | ||
+ | </p> | ||
+ | <p> | ||
+ | The Inventive Step is the most relevant variable but is regulated by the industrial application. The highest inventive step can be reduced by almost a factor 2 with the lowest industrial application. On the contrary, the lowest inventive step is not affected by the industrial application. After computing all the discussed factors, the patentability index can run from 0 to 10. The reader can find a detailed discussion of the questionnaires, tests, analysis of the index and references in the <b><a href="HumanPractices"><strong>Open License section</strong></a></b>. | ||
+ | </p> | ||
</div> | </div> | ||
<div class="tab-pane fade" id="OtherTab"> | <div class="tab-pane fade" id="OtherTab"> | ||
<h2>Applications to other projects</h2> | <h2>Applications to other projects</h2> | ||
</html> | </html> | ||
- | |||
<html> | <html> | ||
<p> | <p> | ||
- | + | In our project, we have developed a tool that allows us to predict the stability of <em>E. coli </em>carrying a reporter gene (in particular, the stability when subjected to different temperatures). The variation of activity with temperature was described by: | |
+ | </p> | ||
+ | $$ | ||
+ | \begin{equation} | ||
+ | \frac{A_{c}}{A^0_{c}} = \prod_i \frac{V_i}{V^0_i} , | ||
+ | \end{equation} | ||
+ | $$ | ||
+ | where $V_i$ is a temperature dependent factor, and A<sub>c</sub><sup>0</sup> is the activity in reference conditions (in this case, 37 ºC). | ||
</p> | </p> | ||
<p> | <p> | ||
+ | We found (<a href="#StabilityTab">see Stability</a>): | ||
$$ | $$ | ||
\begin{equation} | \begin{equation} | ||
- | + | \frac{Ac}{Ac^0} = \frac{N_Pr}{N^0_Pr} =\frac{r_T}{r^0_T} \frac{r_P}{r^0_P}, | |
\end{equation} | \end{equation} | ||
$$ | $$ | ||
- | + | with | |
- | + | ||
$$ | $$ | ||
\begin{equation} | \begin{equation} | ||
- | + | \frac{r_T}{r^0_T} = \frac{Prob_B}{Prob^0_B} = \frac { 1+\frac{R^0}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }} } { \frac {R^0}{R}+\frac{R^0}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }}}, | |
\end{equation} | \end{equation} | ||
$$ | $$ | ||
+ | where | ||
+ | $$ | ||
+ | \begin{equation} | ||
+ | \frac{R^0}{R} = \frac {0.027 (T^0)^2 - 2.12 T^0 + 42.2} {0.027 T^2 - 2.12 T + 42.2} , | ||
+ | \end{equation} | ||
+ | $$ | ||
+ | and with $$ | ||
+ | \begin{equation} | ||
+ | \frac{r_P}{r^0_P} = \frac{\frac{5}{4}}{\frac{r^0_E}{r_E}+\frac{1}{4}}, | ||
+ | \end{equation} | ||
+ | $$ | ||
+ | where | ||
+ | $$ | ||
+ | \begin{equation} | ||
+ | r_E = \frac{1}{\frac{1}{A_1}e^{\frac{E^1_a}{RT}}+\frac{1}{A_2}e^{\frac{E^2_a}{RT}}+\frac{1}{A_3}e^{\frac{E^3_a}{RT}}}. | ||
+ | \end{equation} | ||
+ | $$ | ||
+ | In more complex systems, one expects more temperature dependent factors, and varying more significantly. For an enzyme, the activity is: | ||
+ | $$ | ||
+ | \begin{equation} | ||
+ | A_{c} = k e^{-\frac{E_a}{RT } } \frac{N_{Pr}}{1+e^{\frac{\Delta G_F}{RT } } } . | ||
+ | \end{equation} | ||
+ | $$ | ||
+ | where N<sub>Pr</sub> is the number of total proteins, ΔG is the free energy of the folding process, k is the Arrhenius constant of the chemical reaction | ||
+ | and E<sub>a</sub> is the activation energy of the catalyzed reaction. | ||
</p> | </p> | ||
<p> | <p> | ||
- | The | + | A crude estimation can be obtained by assuming: a) all genes are controlled by the same promoter, b) same folding energy for all proteins, and c) same activation energy |
+ | for all reactions. The activation energy for each reaction (mean activation energy is from [1]) is: | ||
+ | $$ | ||
+ | \begin{equation} | ||
+ | \frac{A_{c}}{A^0_{c}} = {(\frac{r_P}{r^0_P} )^n}{({\frac{r_P}{r^0_P}}{\frac{\frac{1}{1+e^{\frac{\Delta G_F}{RT } } }}{\frac{1}{1+e^{\frac{\Delta G_F}{RT^0 } } }}}{\frac{e^{-\frac{E_a}{RT } }}{e^{-\frac{E_a}{RT^0 } }}})^m}, | ||
+ | \end{equation} | ||
+ | $$ | ||
+ | where n is the number of different mRNAs, m is the number of proteins, E$_ a$ is the activation energy of the catalyzed reaction and $\Delta$G is the protein folding energy. Using: | ||
+ | $$ | ||
\begin{equation} | \begin{equation} | ||
\Delta G_F (L,T) = \Delta H(L) + \Delta C_p (L)(T-T_h) - T\Delta S(L) - \Delta C_p (L)T\ln\frac{T}{T_S} \nonumber, | \Delta G_F (L,T) = \Delta H(L) + \Delta C_p (L)(T-T_h) - T\Delta S(L) - \Delta C_p (L)T\ln\frac{T}{T_S} \nonumber, | ||
\end{equation} | \end{equation} | ||
- | + | $$ | |
+ | with | ||
+ | $$ | ||
\begin{eqnarray} | \begin{eqnarray} | ||
\Delta H(L) &=& -5.03 L - 41.6 (kJ/mol)\nonumber, \\ | \Delta H(L) &=& -5.03 L - 41.6 (kJ/mol)\nonumber, \\ | ||
Line 584: | Line 652: | ||
\Delta S(L) &=& -16.8 L -85 (J/mol)\nonumber. | \Delta S(L) &=& -16.8 L -85 (J/mol)\nonumber. | ||
\end{eqnarray} | \end{eqnarray} | ||
+ | $$ | ||
+ | and T<sub>h</sub> = 373.5 K, T<sub>s</sub> = 385 K and L = 300 aa is the <em>E. coli</em> mean protein length is 300 aa, we computed the | ||
+ | relative activity and relative protein number for n mRNAs and m proteins. The results are shown below: | ||
+ | </p> | ||
+ | <div class="thumbnail col-sm-7 center-block"> | ||
+ | <a href="https://static.igem.org/mediawiki/2014/5/59/Vbt_modeling_other_figure1.png" rel="lightbox"> | ||
+ | <img src="https://static.igem.org/mediawiki/2014/5/59/Vbt_modeling_other_figure1.png" alt="Figure 1" /> | ||
+ | </a> | ||
+ | <div class="caption"> | ||
+ | <strong>Figure 1:</strong> | ||
+ | Variation of activity with temperature for a system that has 1, 2, 3, 4 and 5 genes. | ||
+ | </div> | ||
+ | </div> | ||
+ | <div class="thumbnail col-sm-7 center-block"> | ||
+ | <a href="https://static.igem.org/mediawiki/2014/d/d9/Vbt_modeling_other_figure2.png" rel="lightbox"> | ||
+ | <img src="https://static.igem.org/mediawiki/2014/d/d9/Vbt_modeling_other_figure2.png" alt="Figure 2" /> | ||
+ | </a> | ||
+ | <div class="caption"> | ||
+ | <strong>Figure 2:</strong> | ||
+ | Variation of synthetized protein with temperature for a system that has 1, 2, 3, 4 and 5 genes. | ||
+ | </div> | ||
+ | </div> | ||
+ | <div class="thumbnail col-sm-7 center-block"> | ||
+ | <a href="https://static.igem.org/mediawiki/2014/b/b6/Vbt_modeling_other_figure3.png" rel="lightbox"> | ||
+ | <img src="https://static.igem.org/mediawiki/2014/b/b6/Vbt_modeling_other_figure3.png" alt="Figure 3" /> | ||
+ | </a> | ||
+ | <div class="caption"> | ||
+ | <strong>Figure 3:</strong> | ||
+ | Variation of activity with temperature for a system that has 1, 2, 3, 4 and 5 enzymes, but only one mRNA. | ||
+ | </div> | ||
+ | </div> | ||
+ | <h3> | ||
+ | |||
+ | <div class="thumbnail col-sm-7 center-block"> | ||
+ | <a href="https://static.igem.org/mediawiki/2014/b/b2/VBT_graf7.png" rel="lightbox"> | ||
+ | <img src="https://static.igem.org/mediawiki/2014/b/b2/VBT_graf7.png" /> | ||
+ | </a> | ||
+ | <div class="caption"> | ||
+ | <strong>Figure 4:</strong> | ||
+ | Variation of activity with temperature of a system that has 5 proteins, | ||
+ | using 1, 2, 3, 4 and 5 mRNA. | ||
+ | </div> | ||
+ | </div> | ||
+ | |||
+ | </html> | ||
+ | === Summary === | ||
+ | <html> | ||
+ | <p> | ||
+ | - Stability decreases with increasing modules number. In case of many modules, stability is larger when one of the modules is the limiting step. | ||
</p> | </p> | ||
<p> | <p> | ||
- | + | - Stability grows with decreasing number of mRNA at lower temperatures. On the other hand, higher temperatures and lower number of mRNA show higher activity. | |
+ | </p> | ||
+ | |||
+ | </html> | ||
+ | === Application to other iGEM projects === | ||
+ | <html> | ||
+ | <p> | ||
+ | <a href="https://2010.igem.org/Team:Slovenia/ABOUT_US" target="blank">Team Slovenia 2010 (DNA CODING BEYOND TRIPLETS)</a> | ||
+ | developed a system to use a DNA sequence as a scaffold for enzymes, in order to optimize a metabolic pathway (<a href="https://2010.igem.org/Team:Slovenia" target="blank"> TS10 </a>). We apply our model as follows: | ||
+ | </p> | ||
+ | <p> | ||
+ | $$ | ||
\begin{equation} | \begin{equation} | ||
- | A_{c} = | + | \frac{A_{c}}{A^0_{c}} = {(\frac{r_P}{r^0_P} )^n}{({\frac{r_P}{r^0_P}}{\frac{\frac{1}{1+e^{\frac{\Delta G_F}{RT } } }}{\frac{1}{1+e^{\frac{\Delta G_F}{RT^0 } } }}}{\frac{e^{-\frac{E_a}{RT } }}{e^{-\frac{E_a}{RT^0 } }}})^m}{({{\frac{\frac{1}{1+e^{\frac{\Delta G_Z}{RT } } }}{\frac{1}{1+e^{\frac{\Delta G_Z}{RT^0 } } }}}}{{\frac{\frac{1}{1+e^{\frac{\Delta G_I}{RT } } }}{\frac{1}{1+e^{\frac{\Delta G_I}{RT^0 } } }}}})^z}{{\frac{\frac{1}{1+e^{\frac{\Delta G_D}{RT } } }}{\frac{1}{1+e^{\frac{\Delta G_D}{RT^0 } } }}}} |
\end{equation} | \end{equation} | ||
- | where | + | $$ |
+ | where n is the number of required transcripts, m is the number of required translated proteins, z is the number of required zinc-finger—DNA interactions, E$_a$ is the activation energy, $\Delta$G$_F$ is the folding energy of any enzyme required for the pathway, $\Delta$G$_Z$ is the folding energy of every zinc-finger domain [17], $\Delta$G$_I$ is the interaction energy of a normal zinc-finger-DNA interaction [17] and $\Delta$G$_D$ is the DNA melting energy. | ||
+ | </p> | ||
+ | <p> | ||
+ | The goal was to look for a more efficient synthesis of violacein. Using the mean activation energy (50 kJ/mol) for any catalyzed reaction [1] we get | ||
+ | </p> | ||
+ | <div class="thumbnail col-sm-7 center-block"> | ||
+ | <a href="https://static.igem.org/mediawiki/2014/8/8f/Vbt_modeling_other_figure4.png" rel="lightbox"> | ||
+ | <img src="https://static.igem.org/mediawiki/2014/8/8f/Vbt_modeling_other_figure4.png" alt="Figure 4" /> | ||
+ | </a> | ||
+ | <div class="caption"> | ||
+ | <strong>Figure 5:</strong> | ||
+ | Variation of activity with temperature of a metabolic pathway with 1, 2, 3, 4 and 5 enzymes. | ||
+ | </div> | ||
+ | </div> | ||
+ | We conclude that the new modules (zinc-finger and DNA) do not significantly affect the stability. | ||
+ | <p> | ||
+ | <a href="https://2013.igem.org/Team:Heidelberg/Team" target="blank">Team Heidelberg 2013 (PHILOSOPHER’S STONE)</a> | ||
+ | developed a system to synthesize non ribosomal peptides. One of the simplest mechanisms was the synthesis of indigoidine (<a href="https://2013.igem.org/Team:Heidelberg" target="blank">TH13 wiki</a>) : | ||
+ | </p> | ||
+ | <div class="thumbnail col-sm-7 center-block"> | ||
+ | <img src="https://static.igem.org/mediawiki/2014/c/cc/Vbt_modeling_other_schema.png" alt="" /> | ||
+ | </div> | ||
+ | For this synthesis, we need: 1) acquisition of Glu (assuming that the bacteria acquires from the external medium with a specific transporter), 2) transform Glu in cGlu and 3) transform cGlu in Ind (assuming that there ir one lower step). We get: | ||
+ | $$ | ||
+ | \begin{equation} | ||
+ | \frac{A_{c}}{A^0_{c}} = {(\frac{r_P}{r^0_P} )^3}{(\frac{r_P}{r^0_P})^3}{({\frac{\frac{1}{1+e^{\frac{\Delta G_F}{RT } } }}{\frac{1}{1+e^{\frac{\Delta G_F}{RT^0 } } }}}{\frac{e^{-\frac{E_a}{RT } }}{e^{-\frac{E_a}{RT^0 } }}})^2}{\frac{\frac{1}{1+e^{\frac{\Delta G_T}{RT } } }}{\frac{1}{1+e^{\frac{\Delta G_T}{RT^0 } } }}}{\frac{e^{-\frac{E_aT}{RT } }}{e^{-\frac{E_aT}{RT^0 } }}} | ||
+ | \end{equation} | ||
+ | $$ | ||
+ | where E$_a$ is the activation energy of every chemical reaction, E$_{aT}$ is the activation energy of the transport process, $\Delta$G$_F$ is the folding energy of any enzyme required for the pathway and $\Delta$G$_T$ is the folding energy of the transporter. | ||
+ | </p> | ||
+ | <p> | ||
+ | We use the mean activation energy (50 kJ/mol) for any catalyzed reaction [1] and the activation energy for the Lysine transport through membranes derived from [18], | ||
+ | E$_{aT}$=33.244 kJ/mol, we get | ||
+ | </p> | ||
+ | <div class="thumbnail col-sm-7 center-block"> | ||
+ | <a href="https://static.igem.org/mediawiki/2014/5/50/Vbt_modeling_other_figure5.png" rel="lightbox"> | ||
+ | <img src="https://static.igem.org/mediawiki/2014/5/52/VBT_graf8.png" alt="Figure 5" /> | ||
+ | </a> | ||
+ | <div class="caption"> | ||
+ | <strong>Figure 6:</strong> | ||
+ | Variation of the Ind production with temperature compared to the variation of activity of the systems composed by 2 or 3 enzymes. | ||
+ | </div> | ||
+ | </div> | ||
+ | <p> | ||
+ | In this case, we observe deviations to the expectation: the optimum temperature and the maximum value of Ac/Ac$^0$ are lower than in the n,m = 3 case. This implies a significant contribution of the Glu transport module. | ||
</p> | </p> | ||
Line 642: | Line 815: | ||
<li>Paulsson J and Ehrenberg M (2001) <i>Noise in a minimal regulatory network: plasmid copy number control</i> Cambridge Press <a href="http://journals.cambridge.org/action/displayFulltext?type=1&fid=75292&jid=QRB&volumeId=34&issueId=01&aid=75291"><b>link</b></a> | <li>Paulsson J and Ehrenberg M (2001) <i>Noise in a minimal regulatory network: plasmid copy number control</i> Cambridge Press <a href="http://journals.cambridge.org/action/displayFulltext?type=1&fid=75292&jid=QRB&volumeId=34&issueId=01&aid=75291"><b>link</b></a> | ||
+ | |||
+ | <li>Shiro I and Kuldell N (2007) <i>Zinc Finger proteins: from atomic contact to cellular function</i>. Landes Bioscience. <a href="http://www.landesbioscience.com/iu/Iuchi_9781587066429.pdf"><b>link</b></a> | ||
+ | |||
+ | <li>Berlin, R.D (1973) <i>Temperature dependence of nucleoside membrane transport in rabbit alveolar macrophages and polymorphonuclear leukocytes</i>. J. Biol. Chem., 248:4724-4730. <a href="http://www.ncbi.nlm.nih.gov/pubmed/4718740"><b>link</b></a> | ||
</ol> | </ol> |
Latest revision as of 23:35, 17 October 2014
Modeling
Presentation
In this section, we present our models to describe the four legs of the ST$^2$OOL. We have closely followed all the experimental activities, so our models do not pretend to completely characterize a single dynamical process of a cell, but instead we intend to capture the dependences of our systems on the variables explored by our experiments.
Stability We have modeled the variation of the activity of cells transformed with Biobricks, when subjected to different temperatures or pH, by considering changes in the Biobrick transcription and translation rates.
Standardization We have modeled the expression in different strains by comparing the transcription rates. We have built an app that computes the standardization-stability index and generates the promoter sequence for a given expression rate.
Orthogonality We have studied deterministic and stochastic models of the plasmid growth in order to learn about the unspecific interaction between two plasmids in a cell.
Open License We have modeled the susceptibility to patent a human creation and define the factors that influence it .
Stability
Living organisms have genetic and biochemical tools that allow them to respond to stress. These tools normally up-regulate some stress-specific genes, and, directly or indirectly, down-regulate the rest of normal genes (most of them housekeeping genes). On the other hand, every Escherichia coli component works worse under non-optimal conditions (for example, nutrients uptake or ATP production), so that, every E. coli function will be reduced. We model the relative expression of biobricks when subjected to variations of temperature, pH and salt concentration.
Protein production
The mRNA and protein production of a cell can be described by the kinetic equations:
$$ \begin{eqnarray} \frac{d N_{mRNA} }{dt}&=& P \cdot r_T - N_{mRNA} \cdot \lambda_{mRNA} \nonumber ,\\ \frac{d N_{Pr} }{dt}&=& N_{mRNA} \cdot r_P - N_{Pr} \cdot \lambda_{Pr}, \nonumber \end{eqnarray} $$
where $N_{\it{mRNA}} (N_{Pr} )$ is the number of biobrick mRNA (protein) molecules per cell, $\lambda_{mRNA} $ ($\lambda_{Pr}$) is the mRNA (protein) decay rate, $r_T$ ($r_P$) is the transcription (translation) rate, P is the mean number of plasmids per cell.
Assuming stationary messenger and protein content during the cell cycle, $\frac{d N_{mRNA} }{dt}=0$ and $\frac{d N_{Pr} }{dt}=0$, we get
$$ \begin{equation} N_{Pr} = \frac{ r_T \cdot r_P }{ \lambda_{mRNA}\cdot \lambda_{Pr} } P . \end{equation} $$This equation describes the amount of total synthesized protein as a function of calculable rates and the number of free plasmids. It can be easily shown that, after many generations, the mean number of proteins per cell is 1.5 times the protein production during a cell cycle.
1) TEMPERATURE
The expression level of our heterologous gene can be calculated by equations 1 and 2. Our main goal is to find which of the terms will change most significantly with the different stress factors. We are interested in the activity in comparison to a reference one, which we take as the activity of E. coli bacteria grown on the LB medium at 38.7 oC, pH 7, 1% of NaCl, no UV irradiation and 1 atm of pressure.
$$ \begin{equation} \frac{A_{c}}{A^0_{c}} = \prod_i \frac{V_i}{V^0_i} . \end{equation} $$ where $V$ is any of the factors contributing to the activity and the index $0$ refers to the reference case.
Transcription rate
Among several factors, we will first discuss the transcription rate. In the range of temperatures (or pH) that we have studied, the polymerase (RNAP) is sufficiently stable to assume that the elongation rate is not affected and the interaction energies RNAP-promoter do not change significantly (for example, the promoter specific interaction only changes 16%, from 33 oC to 44 oC). We did not find any significant difference in the nucleotide compositions of E. coli under stress [2]. Therefore, main variations in the transcription rate are due to the probability of RNAP-promoter bounds $$ \begin{equation} r_T (T)=\frac{1}{\frac{1}{r_I(\Delta\epsilon)}+\frac{1}{r_E}}Prob_B (\Delta\epsilon, T) \propto Prob_B (T). \nonumber \end{equation} $$ On the other hand, down-regulation affects the expression level of our Biobrick, because it works as a housekeeping gene (with a normal constitutive promoter for the RNAP). For example, E. coli has a common σ factor for housekeeping genes (σ70), and, another one (σ32) for temperature stress. In the case of temperature stress, the bacteria will synthetize σ32, which will bind RNA polymerase (RNAP) cores that won’t be free for σ70, and so they will not be able to synthetize mRNA of the normal housekeeping genes.
Therefore, we just need to find a relation between the amount of σ factors and the expression level, which we derive next with thermodynamic arguments, and also the variation of σ factors with temperature, which can be determined experimentally.
The RNAP is bound to a promoter (unspecific DNA) with $\epsilon_{S}$ ($\epsilon_{NS}$) free energy. There are
$N_{NS}$ nonspecific sites on the E. coli genome (sequences of 41bp) and R free E$\sigma$70, with $R \ll N_{NS}$ .
The combinatory number of unbound states, with energy R$\epsilon_{NS}$, is
$$
\begin{equation}
\frac{N_{NS}!}{R!(N_{NS}-R)!} \simeq \frac{N^R_{NS}}{R!} .
\nonumber
\end{equation}
$$
Similarly, the combinatory number with one bound state, with energy $\epsilon_{S}$+(R-1)$\epsilon_{NS}$, is
$$
\begin{equation}
\frac{N_{NS}!}{(R-1)!(N_{NS}-R-1)!} \simeq \frac{N^{R-1}_{NS}}{(R-1)!} .
\nonumber
\end{equation}
$$
Then, the probability of finding a RNAP bound to a promoter is with the enter
$$
\begin{equation}
Prob_B = \frac {\frac{R}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }}}{1+\frac{R}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }}},
\end{equation}
$$
with $\Delta\epsilon = \epsilon_{S} - \epsilon_{NS}$, calculated with the energy matrix developed in [3] and $k_B$ is
the Bolztmann constant. R $\simeq$ 1000 and $N_{NS} \simeq 10^7$ for E. coli.
Using the measured variation of free E$\sigma$70, in arbitrary units, with temperature, in $^o$C, [4]
$$
\begin{equation}
E\sigma70 = 0.027 T^2 - 2.12 T + 42.2,
\nonumber
\end{equation}
$$
we can determine the relative change of transcription rate
$$
\begin{equation}
\frac{r_T}{r^0_T} = \frac{Prob_B}{Prob^0_B} = \frac { 1+\frac{R^0}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }} } { \frac {R^0}{R}+\frac{R^0}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }}} ,
\end{equation}
$$
with
$$
\begin{equation}
\frac{R^0}{R} = \frac {0.027 (T^0)^2 - 2.12 T^0 + 42.2} {0.027 T^2 - 2.12 T + 42.2} .
\end{equation}
$$
Synthesis rate
Modeling synthesis rate follows the lines discussed in transcription rate, but the working regime is quite different. The energies of the mRNA-ribosome sequences interactions are much larger than in the RNAP-promoter ones. The average difference of free energies (computed with the RBS calculator for random sequences of E. coli genes) is $\Delta\epsilon = -20.79 kT$, which means that, in the equation: $$ \begin{equation} \frac{r_T}{r^0_T} = \frac{Prob_B}{Prob^0_B} = \frac { 1+\frac{R^0}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }} } { \frac {R^0}{R}+\frac{R^0}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }}} , \end{equation} $$ The term $\frac{R^0}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }}>>>\frac {R^0}{R}$. With this, and knowing that there is no effect on the ribosomal protein composition with temperature changes [5], we conclude that the most signifficant factor should be the elongation rate.
Therefore, $$ \begin{equation} \frac{r_P}{r^0_P} = \frac{1+\frac{r^0_E}{r^0_I}}{\frac{r^0_E}{r_E}+\frac{r^0_E}{r^0_I}}, \end{equation} $$ where $r_E$ ($r_I$) is the elongation (initiation) rate and the superscript denotes our reference case. In E. coli, most of the ribosomes are translating (4 active every one free) [6], thus $\frac{r^0_E}{r^0_I} \simeq\frac{1}{4}$ and $$ \begin{equation} \frac{r_P}{r^0_P} = \frac{\frac{5}{4}}{\frac{r^0_E}{r_E}+\frac{1}{4}}. \end{equation} $$We discuss next the dependence of elongation rates with temperature. We use the relation of amino acid concentration with temperature [7]. Using data from [8], we know that small variations on the ternary complex does not affects the elongation rate. The slowests steps are: 1) GTP hydrolysis on ternary complex (0.01 s), 2) dipeptide bond formation (0.03s) and 3) translocation (0.022s). $$ \begin{equation} r_E = \frac{1}{\frac{1}{A_1}e^{\frac{E^1_a}{RT}}+\frac{1}{A_2}e^{\frac{E^2_a}{RT}}+\frac{1}{A_3}e^{\frac{E^3_a}{RT}}}, \end{equation} $$ where $E^i_a$ is (32, 75.2, 29.260) kJ/mol and $A^{-1}_i$ is (2.5 $\cdot 10^7$, 1.6 $\cdot 10^{14}$, 3.9 $\cdot 10^6$).
Summary
Our model includes two main contributions to the changes of the enzymatic activity due to temperature variation: transcription and synthesis rates. In the case of the GFP, the folding free energy term is small and the relative enzymatic activity can be simplified to our master equation for stability $$ \begin{equation} \frac{Ac}{Ac^0} = \frac{r_T}{r^0_T} \frac{r_P}{r^0_P} . \end{equation} $$
2) pH
Contrary to the temperature situations, internal pH hardly varies in E. coli when the external pH does. Therefore, we focus on the intrinsic E. coli element that may be affected by this type of stress, subsequently influencing Biobrick expression. The main factor is the number of free RNAP, which will be reduced due to the synthesis of stress response genes. We estimate the ratio of free RNAP under different conditions by the ratio of their growth rates $$ \begin{equation} \frac{r_T}{r^0_T} = \frac{Prob_B}{Prob^0_B} = \frac { 1+\frac{R^0}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }} } { \frac {R^0}{R}+\frac{R^0}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }}} , \end{equation} $$ with $$ \begin{equation} \frac{R^0}{R} = \frac{growth^0} {growth} . \end{equation} $$
Comparison with wet-lab data
Our model captures the shape of the observed temperature profile and predicts that additional variations at lower temperatures are independent of the chosen Biobrick. We have found that transcription rate decreases with temperature as a consequence of the reduction of free RNAP; and synthesis rate grows above the optimal temperature. Their combination leads to the observed variation of fluorescence with temperature. Unfortunately, this is not the whole story. Our model fails to fit the amount of variation, which is larger experimentally than the predicted one. By comparing our theoretical predictions and experimental results of XL1 Blue strain transformed with Bb1 or Bb3, we find that the ratio between theoretical prediction and experimental results is the same at under-optimal temperatures (in fact, accidentally too close given the error bars). This fact suggests that the temperature dependence not included in our model should come from cell mechanisms; independent from Biobrick parameters (promoter strength, ...). On the contrary, the comparison of theoretical predictions and experimental results at higher temperatures show that our model might require more dynamics which depends on Biobrick parameters.
Our model captures the shape of the observed pH profile. We have found that the lowest fluorescence happens at the pH with lowest growth, while it is mildly dependent on the pH from 6 to 8.
Standardization
We will now discuss Biobrick expression in different E. coli strains. The effect of the changes in the number of free RNAP on mRNA synthesis rate depend on the energy of the promoter. The relative transcription rate in different strains can be written as $$ \begin{equation} \frac{r_T}{r'_T} = \frac{Prob_B}{Prob'_B} = \frac { 1+\frac{R'}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }} } { \frac {R'}{R}+\frac{R'}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }}} , \end{equation} $$ where we can see that the dominance of the second term over the first one, controlled by the promoter strength, would lead to similar rate independently on the particular strain chosen. Comparing this formula with the relative transcription rate in stability, we find that promoters whose binding energy makes them more standard also implies more stable promoters. Strong promoter sequences experience a decrease in synthesis rates when subjected to an increase of promoter energy. We use these results and the energy matrix derived in [3] to create an application that is able to compute promoter sequences for a given relative (to the maximum) expression rate and associated standard and stable rate. The physical parameters behind are the promoter strength (using the data of Anderson’ promoters) and the binding energy and the tool create a sequence of a constitutive promoter with the appropriated trade-off between strength and energy, therefore, with the appropriated trade-off between stability & standardization (ST$^2$)and rate of expression. We define the ST$^2$ index as $$ F_{st2} = (F - 0.012)/0.51 $$ with $$F = \frac{R'}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }}, $$ where $F_{sty}$ has been scaled to run from 1 to 0. We have used a fitting formula for the rate of expression relative to the maximum rate [3] as a function of the binding energy $\Delta\epsilon$, $$ r_T = 2.29 + 0.27 \frac{\Delta\epsilon}{kT}. $$ With this approach, we can infer the ST$^2$ level and the promoter sequence for a given expression rate. We show below our app in javascript.
Generate a promoter sequence and compute the ST$^2$ (Stability & Standardization) index for a given expression rate
Move the cursor to a given expression rate (as a fraction of its maximum) and press Generate a chain
25%
Promoter sequence (from -41 to -1 nt) is:
ST$^2$ index:
Comparison with wet-lab data
Our wet-lab results show that Biobrick 5 is the most standard one.Our model predicts that systems with less modules or factors will be more stable and standard. Biobrick 5 has a tetracycline depending promoter, whose limitating step is the binding of tetracycline to TetR. In this scenario, the intracellular concentration of aTc and the binding to the TetR repressor do not depend on the strain, so there is only one factor in which the strain could have an important effect: the TetR intracellular concentration, which has a constitutive expression. Analyzing the other Biobricks, we found that the strains with less expression of the constitutive promoter are those ones with higher expression when we use the Ptet promoter (but JM109 and BB1), which agree with our prediction.
Orthogonality
Plasmid dynamics have been studied exhaustively over these years using mathematical models: from how plasmid growth is regulated to plasmid partitioning after cell division. However, little is known about how two identical plasmids interact with each other. This part of our modelling project aims to shed light upon this issue using previous models of ColE1 plasmids and the properties of two different plasmids with same origin of replication when these are introduced in a cell. We employed both deterministic and stochastic models to try to answer whether two plasmids behave independently, and developed some approximations to highlight important features of the system.
Deterministic model
This model is taken from [16]. Here, plasmid levels are regulated both by an RNA molecule that inhibits plasmid replication and by the dilution of plasmids as a result of cell growth and division:
$$ \frac{d P}{dt} = r \cdot P \cdot R(s) - \frac{\log 2}{\tau} \cdot P, $$
$$ \frac{d s}{dt} = \alpha \cdot P - \beta \cdot s, $$
where $P$ is the number of plasmids, $s$ is the number of RNA inhibitors, $r$ is the plasmid growth rate, $\tau$ is the cell growth rate, $\alpha$ is the RNA inhibitor synthesis rate
, $\beta$ is the RNA inhibitor decay rate and $R(s)$ is the repression effect of RNA inhibitor with
$$ R(s) = \frac{1}{(1+\frac{s}{s_0})^n},$$
where $s_0$ is the value of RNA inhibitor necessary to repress half of the replication events and $n$ is the cooperativity. Plasmid levels are regulated both by a negative feedback and by cell growth. If we consider several plasmids, we get
$$ \frac{d P_i}{dt} = r_i \cdot P_i \cdot R(s) - \frac{\log 2}{\tau} \cdot P_i,$$
$$ \frac{d s}{dt} = \alpha \cdot \sum \limits_{i=1}^n{P_i} - \beta \cdot s,$$
with a steady state solution
$$ s_{ss} = R^{-1}(\frac{\log 2}{\tau \cdot r_i}),$$
$$ \sum \limits_{i=1}^n{P_{i,ss}} = (\frac{ \beta \cdot s_{ss}}{\alpha}).$$
In a simplified way, with two identical plasmids, we get:
$$ P_{1,ss} + P_{2,ss} = K.$$
It is noticeable that plasmid levels are not independent. Small perturbations can increase the relative level of one plasmid, favoring its growth. As they are repressed as a function of the total number of plasmids but they grow depending on their number, the system is prone to introduce asymmetries
$$ \frac{d P_1}{d P_2} = \frac{r_1 \cdot P_1}{r_2 \cdot P_2} .$$
If growth rates are constant, the slope is determined by the initial conditions
$$ m = \frac{P_{01}}{P_{02}} .$$
Quasi-steady state approach
To simplify our equations, we assume inhibitor dynamics are faster than plasmid growth, so RNA levels are a direct representation of plasmid levels $$ \frac{d P_i}{dt} = r_i \cdot P_i \cdot R(P_i) - \frac{\log 2}{\tau} \cdot P_i .$$Dependence on initial conditions:
To test our hypothesis, we integrate our differential equations using numerical methods with several initial conditions. As expected, small differences in initial proportions are amplified by the dynamics.Stochastic model
The deterministic model assumes negligible noise. However, this might not be the case and small variations in plasmid levels could be highly amplified. Hence, we now discuss the importance of noise in plasmid growth. We used the Stochastic Simulation Algorithm (SSA) developed by Gillespie (see [13]). In short, SSA simulates by iteration which processes occur at any time using a random number generator.
First of all, we used the algorithm to develop a code assuming growth of the plasmids and cell division as independent processes. The steps in the algorithm are:
We show below the evolution of a two plasmids system for two different initial conditions: $P_1=P_2=10$ and $P_1 = 10$ and $P_2 = 50$.
As expected, the deterministic model gives a mean version of the dynamics. Single cell dynamics shows oscillatory behaviour. The levels of each plasmid can vary but their sum remains unchaged (black line)
A lineage model
Aiming to obtain a better description, we studied plasmid levels in a lineage from a single cell.As tracking bacteria and their plasmid levels is an important issue, we implemented a tree plot in our code using a matlab package. Plasmids were given to each daughter cell following a gaussian distribution with a standard deviation of $\frac{1}{\sqrt{n}}$
Also we were able to track the plasmid levels of each bacteria at the time of their division.
The distribution at the beginning and at the end of the simulation is:
Comparing with a distribution of orthogonal plasmids, we find that while both are symmetric (respect to the $P_1 = P_2$ line), non orthogonal plasmids lead to broader heat maps.
A step beyond: Resistance
As a proof of concept, we modelled the bacterial time of division as a function of both plasmids, being minimal when both plasmid levels are similar. For this simple regulation, we obtained that the average level of plasmids in a population is the same for each DNA molecule.This way to tune the plasmid levels is unable to make them equal, as depicted in the following picture. Also, there are variations in their growth rates.
Approaches for less computational expensive
Although these simulations are powerful and helpful and yield informative results, they are computationally expensive. This could become an important pitfall when dealing with a huge number of cells or long times. We tried to bypass this problem by using an approach if we are working with populations as explained in [16]. We added a degradation process for each plasmid following the equations:
$$ \alpha_{P_i \rightarrow P_i-1} = \frac{log(2)}{\tau} \cdot P_i$$ This approach can be used to study cell dynamics in a population. With the system of equations, we capture the global dynamics for a cell population. If there is no resistance, we find dynamical examples where a plasmid is lost.
Comparison with wet-lab data
The interaction between two plasmids in the same cell is still an uncharted territory. A relevant example is shown in [15], where the authors obtained an emergent behaviour as result of an unexpected relation between the bacterial physiology and plasmid dynamics. Using models in single plasmid dynamics, we switch between deterministic and stochastic models to learn about the unspecific interaction between two plasmids in a cell. There is an important variability, which drives the system away from orthogonality. Plasmids with identical origin of replication are not orthogonal per se. Although we can introduce different resistances to tune the system, it is hard to stabilize the system due to the large variability. Plasmids with different mechanisms of replication might solve it but it is not guaranteed. Certainly, new plasmids must be designed to make engineered cells more reliable and predictable.
Open License
We model the susceptibility to patent a human creation by $$ P.I. = M \cdot N \cdot A ^{log(\frac{A+I}{2})}, $$ where P.I. is the patentability index, M is the morality, which is a non-legal term accepted in the E.U., and N refers to novelty. Both factors are described by a discrete variable which can be valued 0 or 1 depending on the results of a questionnaire. A is the inventive step and I is the industrial application which can be valued from 1 to 10 depending of the results of a test.
The Inventive Step is the most relevant variable but is regulated by the industrial application. The highest inventive step can be reduced by almost a factor 2 with the lowest industrial application. On the contrary, the lowest inventive step is not affected by the industrial application. After computing all the discussed factors, the patentability index can run from 0 to 10. The reader can find a detailed discussion of the questionnaires, tests, analysis of the index and references in the Open License section.
Applications to other projects
In our project, we have developed a tool that allows us to predict the stability of E. coli carrying a reporter gene (in particular, the stability when subjected to different temperatures). The variation of activity with temperature was described by:
$$ \begin{equation} \frac{A_{c}}{A^0_{c}} = \prod_i \frac{V_i}{V^0_i} , \end{equation} $$ where $V_i$ is a temperature dependent factor, and Ac0 is the activity in reference conditions (in this case, 37 ºC).We found (see Stability): $$ \begin{equation} \frac{Ac}{Ac^0} = \frac{N_Pr}{N^0_Pr} =\frac{r_T}{r^0_T} \frac{r_P}{r^0_P}, \end{equation} $$ with $$ \begin{equation} \frac{r_T}{r^0_T} = \frac{Prob_B}{Prob^0_B} = \frac { 1+\frac{R^0}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }} } { \frac {R^0}{R}+\frac{R^0}{N_{NS}} e^{-\frac{\Delta\epsilon}{k_BT }}}, \end{equation} $$ where $$ \begin{equation} \frac{R^0}{R} = \frac {0.027 (T^0)^2 - 2.12 T^0 + 42.2} {0.027 T^2 - 2.12 T + 42.2} , \end{equation} $$ and with $$ \begin{equation} \frac{r_P}{r^0_P} = \frac{\frac{5}{4}}{\frac{r^0_E}{r_E}+\frac{1}{4}}, \end{equation} $$ where $$ \begin{equation} r_E = \frac{1}{\frac{1}{A_1}e^{\frac{E^1_a}{RT}}+\frac{1}{A_2}e^{\frac{E^2_a}{RT}}+\frac{1}{A_3}e^{\frac{E^3_a}{RT}}}. \end{equation} $$ In more complex systems, one expects more temperature dependent factors, and varying more significantly. For an enzyme, the activity is: $$ \begin{equation} A_{c} = k e^{-\frac{E_a}{RT } } \frac{N_{Pr}}{1+e^{\frac{\Delta G_F}{RT } } } . \end{equation} $$ where NPr is the number of total proteins, ΔG is the free energy of the folding process, k is the Arrhenius constant of the chemical reaction and Ea is the activation energy of the catalyzed reaction.
A crude estimation can be obtained by assuming: a) all genes are controlled by the same promoter, b) same folding energy for all proteins, and c) same activation energy for all reactions. The activation energy for each reaction (mean activation energy is from [1]) is: $$ \begin{equation} \frac{A_{c}}{A^0_{c}} = {(\frac{r_P}{r^0_P} )^n}{({\frac{r_P}{r^0_P}}{\frac{\frac{1}{1+e^{\frac{\Delta G_F}{RT } } }}{\frac{1}{1+e^{\frac{\Delta G_F}{RT^0 } } }}}{\frac{e^{-\frac{E_a}{RT } }}{e^{-\frac{E_a}{RT^0 } }}})^m}, \end{equation} $$ where n is the number of different mRNAs, m is the number of proteins, E$_ a$ is the activation energy of the catalyzed reaction and $\Delta$G is the protein folding energy. Using: $$ \begin{equation} \Delta G_F (L,T) = \Delta H(L) + \Delta C_p (L)(T-T_h) - T\Delta S(L) - \Delta C_p (L)T\ln\frac{T}{T_S} \nonumber, \end{equation} $$ with $$ \begin{eqnarray} \Delta H(L) &=& -5.03 L - 41.6 (kJ/mol)\nonumber, \\ \Delta C_p(L) &=& -0.062 L + 0.53 (kJ/mol)\nonumber, \\ \Delta S(L) &=& -16.8 L -85 (J/mol)\nonumber. \end{eqnarray} $$ and Th = 373.5 K, Ts = 385 K and L = 300 aa is the E. coli mean protein length is 300 aa, we computed the relative activity and relative protein number for n mRNAs and m proteins. The results are shown below:
Summary
- Stability decreases with increasing modules number. In case of many modules, stability is larger when one of the modules is the limiting step.
- Stability grows with decreasing number of mRNA at lower temperatures. On the other hand, higher temperatures and lower number of mRNA show higher activity.
Application to other iGEM projects
Team Slovenia 2010 (DNA CODING BEYOND TRIPLETS) developed a system to use a DNA sequence as a scaffold for enzymes, in order to optimize a metabolic pathway ( TS10 ). We apply our model as follows:
$$ \begin{equation} \frac{A_{c}}{A^0_{c}} = {(\frac{r_P}{r^0_P} )^n}{({\frac{r_P}{r^0_P}}{\frac{\frac{1}{1+e^{\frac{\Delta G_F}{RT } } }}{\frac{1}{1+e^{\frac{\Delta G_F}{RT^0 } } }}}{\frac{e^{-\frac{E_a}{RT } }}{e^{-\frac{E_a}{RT^0 } }}})^m}{({{\frac{\frac{1}{1+e^{\frac{\Delta G_Z}{RT } } }}{\frac{1}{1+e^{\frac{\Delta G_Z}{RT^0 } } }}}}{{\frac{\frac{1}{1+e^{\frac{\Delta G_I}{RT } } }}{\frac{1}{1+e^{\frac{\Delta G_I}{RT^0 } } }}}})^z}{{\frac{\frac{1}{1+e^{\frac{\Delta G_D}{RT } } }}{\frac{1}{1+e^{\frac{\Delta G_D}{RT^0 } } }}}} \end{equation} $$ where n is the number of required transcripts, m is the number of required translated proteins, z is the number of required zinc-finger—DNA interactions, E$_a$ is the activation energy, $\Delta$G$_F$ is the folding energy of any enzyme required for the pathway, $\Delta$G$_Z$ is the folding energy of every zinc-finger domain [17], $\Delta$G$_I$ is the interaction energy of a normal zinc-finger-DNA interaction [17] and $\Delta$G$_D$ is the DNA melting energy.
The goal was to look for a more efficient synthesis of violacein. Using the mean activation energy (50 kJ/mol) for any catalyzed reaction [1] we get
We conclude that the new modules (zinc-finger and DNA) do not significantly affect the stability.Team Heidelberg 2013 (PHILOSOPHER’S STONE) developed a system to synthesize non ribosomal peptides. One of the simplest mechanisms was the synthesis of indigoidine (TH13 wiki) :
We use the mean activation energy (50 kJ/mol) for any catalyzed reaction [1] and the activation energy for the Lysine transport through membranes derived from [18], E$_{aT}$=33.244 kJ/mol, we get
In this case, we observe deviations to the expectation: the optimum temperature and the maximum value of Ac/Ac$^0$ are lower than in the n,m = 3 case. This implies a significant contribution of the Glu transport module.
References
- Ghosh K., Dill K. (2010). Cellular proteomes have broad distributions of protein stability. Biophys J 99: 3996–4002.
- Soini, J. et al. (2005) Transient increase of ATP as a response to temperature up-shift in Escherichia coli. Microb. Cell Fact. 4, 9.
- Brewster R.C., Jones D.L., Phillips R.(2012) Tuning promoter strength through RNA polymerase binding site design in Escherichia coli. PLoS Comput Biol 8: e1002811.
- Skelly S. et al. (1987) Correlation between the 32-kDa s factor levels and in vitro expression of Escherichia coli heat shock genes. Proc. Natl. Acad. Sci. USA 84:8365–8369.
- Zaritsky, A. (1982). Effects of growth temperature on ribosomes and other physiological properties of Escherichia coli. J. Bacteriol. 151:485-486.
- Bremer, H., Dennis, P. P. (1996) Modulation of chemical composition and other parameters of the cell by growth rate. Neidhardt, et al. eds. Escherichia coli and Salmonella typhimurium: Cellular and Molecular Biology, 2nd ed. chapter 97.
- Piperno J.R., Oxender D.L. (1968) Amino acid transport systems in Escherichia coli K-12. J Biol Chem. 243(22):5914–5920.
- Thompson R.C., Dix D.B., Karim A.M. (1986) The reactions of ribosomes with elongation factor Tu·GTP complexes. Aminoacyl-tRNA-independent reactions in the elongation cycle determine the accuracy of protein synthesis. J Biol Chem, 261:4868–4874.
- Bilgin N. et al. (1992) Kinetic properties of Escherichia coli ribosomes with altered forms of S12. J Mol Biol. 224 (4):1011-27.
- Grigorenko B.L. et al. (2008) Mechanism of the chemical step for the guanosine triphosphate (GTP) hydrolysis catalyzed by elongation factor Tu. Biochim Biophys Acta 1784:1908–1917.
- Gindulyte,A. et al. (2006) The transition state for formation of the peptide bond in the ribosome. Proc. Natl. Acad. Sci. U. S. A., 103:13327–13332.
- Wilson,K. S., Noller, H.F. (1998) Mapping the position of translational elongation factor EF-G in the ribosome by directed hydroxyl radical probing. Cell, 92: 131–139.
- Erban R, Chapman SJ and Maini PK (2007) A PRACTICAL GUIDE TO STOCHASTIC SIMULATIONS OF REACTION-DIFFUSION PROCESSES link
- Klumpp E (2011) Growth-Rate Dependence Reveals Design Principles of Plasmid Copy Number Control PLoS One Vol 6 Issue 5 link
- Marguet P, Tanouchi Y, Splitz E, Smith C and You L (2010) Oscillations by Minimal Bacterial Suicide Circuits Reveal Hidden Facets of Host-Circuit Phisology PLoS One Vol 5 Issue 7 link
- Paulsson J and Ehrenberg M (2001) Noise in a minimal regulatory network: plasmid copy number control Cambridge Press link
- Shiro I and Kuldell N (2007) Zinc Finger proteins: from atomic contact to cellular function. Landes Bioscience. link
- Berlin, R.D (1973) Temperature dependence of nucleoside membrane transport in rabbit alveolar macrophages and polymorphonuclear leukocytes. J. Biol. Chem., 248:4724-4730. link