Team:Paris Bettencourt/Modeling

From 2011.igem.org

(Difference between revisions)
(Using diffusion equations)
 
(166 intermediate revisions not shown)
Line 1: Line 1:
{{:Team:Paris_Bettencourt/tpl_test}}
{{:Team:Paris_Bettencourt/tpl_test}}
-
<br>
+
<html>
-
<br>
+
<h1>Modeling</h1>
-
<br>
+
<h2>Modeling in our project</h2>
-
= Modeling introduction =
+
<p>Our modeling was organized around two main questions:
 +
<ul>
 +
<li>Can we <em>explain the transfer</em> through nanotubes?</li>
 +
<li>What will be the <em>behaviour of our constructs</em> and how will it impact our experimental designs?</li>
 +
</ul></p>
-
== What is modeling in synthetic biology? ==
+
<p>Answering those two questions was essential for our project. We needed to know <em>what to expect in order to design our experiments</em> properly and to know what kind of restults we should obtain.</p>
-
Synthetic biology relies heavily on genetic constructs designed to perform a specific task. In order to predict the behaviour of such systems, we need to be able to model it. Knowing in advance how well or how poorly a construct will work can help us during the genetic design phase or when we need to prepare our experiments. We need to describe the future reactions of our system by creating a mathematical model of it. The parameters of the model are then fitted to the experimental results.<br>
 
-
The bible for synthetic biology modeling is Uri Alon's book ''An Introduction to Systems Biology: Design Principles of Biological Circuits''[Ref]. Most of our models are based on his approach to biological circuits. We are now going to explain quickly what is the basic structure behind our simulations.<br>
+
<h2>Investigating nanotube transfer</h2>
-
The gene ''geneX'' is responsible for the production of the corresponding protein ''X''. The promoter ''pX'' controlling the expression ''geneX'' can be:
+
<p>In order to answer the first question, we investigated the <em>physical properties of cell membrane</em> and <em>passive diffusion</em> to comprehend how the transfer could occur. We came up with two different ideas that could explain molecule transfer through the nanotubes, and based our original models on these assumptions, done in Java for passive diffusion and in Matlab for assisted diffusion.These two novel models show that transfer through the nanotubes, whether happening by passive diffusion or the so-called assisted diffusion, is happening too quickly to be accurately measured by fluorescent microscopy. As nanotube transfer is too fast compared to genetic response to allow us to measure its time span correctly, our conclusion was that our designs would not allow us to determine which one of these two processes (passive or assisted diffusion) is dominant during the transfer. Even though this makes it impossible to create a definitive model of molecule transfer through the nanotubes, the information provided by our two alternative models gave us an <em>insight on the time scale of the transfer</em>.</p>
-
* A constitutive promoter, ''geneX'' is activated whatever the conditions
+
<table>
-
* Positively regulated, ''geneX'' is more active when an inducer is present
+
<tr>
-
* Negatively regulated, ''geneX'' is less active when an repressor is present
+
  <td style="width:200px; text-align:center"><a href="https://2011.igem.org/Team:Paris_Bettencourt/Modeling/Diffusion"><img style="width:150px; margin-top:20px;" src="https://static.igem.org/mediawiki/2011/1/1a/Passive-diff-button.png"></a>
 +
  </td>
 +
  <td><b><a href="https://2011.igem.org/Team:Paris_Bettencourt/Modeling/Diffusion">Passive diffusion in nanotubes</a></b> We investigate here the hypothesis of passive diffusion through nanotubes.
 +
  </td>
 +
</tr>
 +
<tr>
 +
  <td style="width:200px; text-align:center"><a href="https://2011.igem.org/Team:Paris_Bettencourt/Modeling/Assisted_diffusion"><img style="width:150px; margin-top:20px;" src="https://static.igem.org/mediawiki/2011/b/b9/Active-diff-button.png"></a>
 +
  </td>
 +
  <td><b><a href="https://2011.igem.org/Team:Paris_Bettencourt/Modeling/Assisted_diffusion">Assisted diffusion</a></b> We propose here a model explaining how diffusion through nanotubes could be "assisted" by the tension differential between cell walls.
 +
  </td>
 +
</tr>
 +
</table>
-
<html><img src='https://static.igem.org/mediawiki/2011/0/03/Constitutive_promoter_general.png' style='width:50%; align:middle'/> <p> Constitutive Promoter</p></html>
 
-
<html><img src='https://static.igem.org/mediawiki/2011/1/16/Positive_promoter_general.png' style='width:50%;'/> <p> Positive autoregulation</p></html>
 
-
<html><img src='https://static.igem.org/mediawiki/2011/0/0c/Negative_promoter_general.png' style='width:50%;'/> <p> Negative autoregulation</p></html>
 
-
The inducer or repressor can be another protein or even the product itself. In the latter case, the gene is auto-regulated, wether positively or negatively. To begin, let's assume ''pX'' is a constitutive promoter for now. We will make another simplification by taking into account the translation step and assuming that the gene "directly" produces protein X.<br>
+
<h2>Predicting the behaviour of our designs</h2>
 +
<p>The second question was a crucial step in the design of our experiments. Taking into account the predicted time scale and other information provided by our two general models, we were able to build <em>models of each of our genetic networks</em>. These models are an improvement of Uri Alon's approach in <i>An Introduction to Systems Biology: Design Principles of Biological Circuits</i> and were done in Matlab. With these models we showed that some designs might work better that the others. For instance, we prioritized the T7 RNA polymerase diffusion and tRNA diffusion systems and decided to concentrate our wet lab experiments and characterizatons in these systems. The ComS system, on the other hand, was less developed because of some disadvantages that our model predicted (high background even without induction, very high activation threshold mainly). Moreover, our models let us <em>evaluate the response time of each of our our constructs</em>. With these estimations, were able to prepare protocols for our <a href ="https://2011.igem.org/Team:Paris_Bettencourt/Experiments/Microscopy">microscopic experiments</a> (by evaluating the characteristic response time of our system, the activation thresholds, etc.).</p>
-
We will know take a look at the parameters involved in modeling this network.
 
-
*''[X]'' is the concentration of protein ''X''
 
-
*<html><img src='https://static.igem.org/mediawiki/2011/a/a8/Alpha.png' style='width:20px;'/></html> is the expression rate of protein ''X'' (molecule.s^-1). It mainly depends on the constitutive promoter.
 
-
*<html><img src='https://static.igem.org/mediawiki/2011/0/05/Delta_dil.png' style='width:24px;'/></html> is the dilution rate, due to cell division (s^-1)
 
-
*<html><img src='https://static.igem.org/mediawiki/2011/6/62/Delta_deg.png' style='width:24px;'/></html> is the degradation rate of protein ''X''
 
-
The equation and solution modeling the behaviour of this system are the following:
 
-
<html><a href='https://2011.igem.org/File:Constitutive.png'><img src='https://static.igem.org/mediawiki/2011/5/5f/Constitutive.png' style='width:60%;'/></a></html>
 
-
Now, let's assume that ''pX'' is auto-regulated, either positively or negatively. We need to introduce a few new parameters.
+
<table>
 +
<tr>
 +
  <td style="width:200px; text-align:center"><a href="https://2011.igem.org/Team:Paris_Bettencourt/what_is_modeling"><img style="width:150px" src="https://static.igem.org/mediawiki/2011/2/2b/Question_mark_button.png"></a>
 +
  </td>
 +
  <td><b><a href="https://2011.igem.org/Team:Paris_Bettencourt/what_is_modeling">The basics about genetic networks modeling</a></b> You can find here an introduction to our methods and the general idea behind most gene network models.
 +
  </td>
 +
</tr>
 +
<tr>
 +
  <td style="width:200px; text-align:center;"><a href="https://2011.igem.org/Team:Paris_Bettencourt/Hypothesis"><img style="width:150px; margin-top:20px;" src="https://static.igem.org/mediawiki/2011/2/21/Hypotheses_button.png"></a>
 +
  </td>
 +
  <td><b><a href="https://2011.igem.org/Team:Paris_Bettencourt/Hypothesis">Our assumptions</a></b> Because of the specificities of our project, we had to adapt the "classic" model to better represent our current situation. On top of that we made and justified a few other hypotheses detailed in this section.
 +
  </td>
 +
</tr>
 +
<tr>
 +
  <td style="width:200px; text-align:center"><a href="https://2011.igem.org/Team:Paris_Bettencourt/Modeling/Designs"><img style="width:150px; margin-top:20px;" src="https://static.igem.org/mediawiki/2011/a/ac/Graph-button.png"></a>
 +
  </td>
 +
  <td><b><a href="https://2011.igem.org/Team:Paris_Bettencourt/Modeling/Designs">Modeling our designs</a></b> Models predicting the behaviour of our designs are detailled in this section.
 +
  </td>
 +
</tr>
 +
</table>
-
*''K'' is the dissociation constant, representing the binding of a inducer/repressor to the promoter
 
-
*''n'' is the Hill coefficient of the function
 
-
*<html><img src='https://static.igem.org/mediawiki/2011/3/3f/Beta.png' style='width:20px;'/></html> is the maximum production rate of our gene (molecule.s^-1)
 
-
*<html><img src='https://static.igem.org/mediawiki/2011/a/a8/Alpha.png' style='width:20px;'/></html> now represents the basal expression ("leaking") of the gene (molecule.s^-)
 
-
The equations and solutions are now:<br>
+
<!-- PAGE FOOTER -- ITEMS FROM COLUMN ! HAVE BEEN MOVED HERE  -- RDR  -->
-
<html><img src='https://static.igem.org/mediawiki/2011/a/af/Positive_equation_general.png' style='width:40%;'/><a href='https://2011.igem.org/File:Positive.png'><img src='https://static.igem.org/mediawiki/2011/6/6b/Positive.png' style='width:60%;'/></a></html>
+
-
<html><img src='https://static.igem.org/mediawiki/2011/a/ac/Negative_equation_general.png' style='width:40%;'/><a href='https://2011.igem.org/File:Negative.png'><img src='https://static.igem.org/mediawiki/2011/d/d2/Negative.png' style='width:60%;'/></a></html>
+
-
You will note that the regulation is modeled as a Hill function. This type of function is used to represent most of the regulation that takes place in a genetic network.<br>
+
        <div id="footer-wrapper">
 +
        <div id="footer">
 +
              <div id="f-poweredbyico"><a href="http://www.mediawiki.org/"><img src="/wiki/skins/common/images/poweredby_mediawiki_88x31.png" height="31" width="88" alt="Powered by MediaWiki" /></a></div>       <div id="f-copyrightico"><a href="http://creativecommons.org/licenses/by/3.0/"><img src="http://i.creativecommons.org/l/by/3.0/88x31.png" alt="Attribution 3.0 Unported" width="88" height="31" /></a></div>     <ul id="f-list">
 +
-
Of course most genetic networks are more complex than simply auto-regulated nodes. The product of one gene can regulate the activation of another which in turn inhibits a third, etc. By coupling this kind of equations together, we achieved modeling most of our genetic networks pretty easily.<br>
+
  <!-- Recentchanges is not handles well DEBUG -->
 +
    <li id="t-recentchanges"><a href="/Special:RecentChanges"
 +
      title='Recent changes'>Recent changes</a></li>
-
== Our choices for modeling ==
+
    <li id="t-whatlinkshere"><a href="/Special:WhatLinksHere/Team:Paris_Bettencourt/Modeling"
 +
      title="List of all wiki pages that link here [j]" accesskey="j">What links here</a></li>
-
=== Steady state flow from the nanotubes ===
+
                        <li id="t-recentchangeslinked"><a href="/Special:RecentChangesLinked/Team:Paris_Bettencourt/Modeling"
 +
                          title="Recent changes in pages linked from this page [k]" accesskey="k">Related changes</a></li>
-
A lot of our designs are based on a simple ''emitter cell'' - ''nanotubes'' - ''receiver cell'' principle. The basic models however can not take into account the nanotubes transfer. We decided to make separate models for particles transmission through the nanotubes and to assume that the number of signaling protein going through them reaches a '''steady state''' in the receiver cell. It can be reached when the flow through nanotubes into the cell exactly compensate protein degradation and dilution. The reason for such a choice is to simplify the interpretation of data. If we were using a non-steady flow of signaling protein, understanding the behaviour of our reporters would have been much more difficult.<br>
 
-
The way we present our models is the following. We first model both emitting and receiving genes in the same cell to give us a control. We then see what would be the response of the receiver genes alone in a cell for different signaling protein steady states. We can then compare our models to reality and evaluate the influence of nanotubes on the response of the system.<br>
 
-
=== Translation ===
+
                <li id="t-upload"><a href="/Special:Upload"
-
We introduced the translation step in our models. This means each equation described above is transformed into two equations: one describing mRNA production and one describing the translation of this mRNA into a protein. We then need to introduce the protein production rate <html><img src='https://static.igem.org/mediawiki/2011/4/4d/Gamma.png' style='width:20px;' /></html> (s^-1). We also need to introduce a different degradation rate for each product. The first one will be for the degradation of ''mRNA'' associated with ''geneX'' and the second one with protein ''X'' itself. The equations for the auto-repression (other regulation systems are similar) are now:
+
                  title="Upload files [u]" accesskey="u">Upload file</a>
 +
                </li>
 +
                <li id="t-specialpages"><a href="/Special:SpecialPages"
 +
                  title="List of all special pages [q]" accesskey="q">Special pages</a>
-
<html><img src='https://static.igem.org/mediawiki/2011/d/d1/MRNA_equation_general.png' style='width:450px;'/></html>
+
                </li>
 +
                <li><a href='/Special:Preferences'>My preferences</a></li>
 +
            </ul>
 +
        </div> <!-- close footer -->
 +
        <div id='footer'>
 +
    <ul id="f-list">
-
This is justified notably because of our tRNA construction. This design relies directly on the translation process, we therefore need to model it if we want to compare models for different designs.
+
            <li id="t-print"><a href="/wiki/index.php?title=Team:Paris_Bettencourt/Modeling&amp;printable=yes"
 +
              title="Printable version of this page [p]" accesskey="p">Printable version</a>
-
=== Delays for maturation ===
+
            </li>
-
We also introduced a delay for protein production and maturation. Most models can ignore this but our experiments rely heavily on time measurements which means that for proteins with maturation time around 5 min (to compare with the cell division time 40 min) we need to take this into account. We chose to model this simply by adding a delay to the response time as it is shown in the following equations:
+
-
<html><img src='https://static.igem.org/mediawiki/2011/b/bc/Delayed_equation_general.png' style='width:450px;'/></html>
+
            <li id="t-permalink"><a href="/wiki/index.php?title=Team:Paris_Bettencourt/Modeling&amp;oldid=86565"
 +
              title="Permanent link to this revision of the page">Permanent link</a>
 +
            </li>
-
The delay for mRNA production is significantly lower than the one for protein production: <html><img src='https://static.igem.org/mediawiki/2011/2/24/T1.png' style='width:20px;'/></html> << <html><img src='https://static.igem.org/mediawiki/2011/4/49/T2.png' style='width:20px;'/></html>.
 
-
=== No delays for diffusion in receptor cell ===
+
        <li id="privacy"><a href="/2011.igem.org:Privacy_policy" title="2011.igem.org:Privacy policy">Privacy policy</a></li>
-
Most of our designs rely on the following scheme: an emitter cell creates a signaling protein which is received through the nanotubes by a receiver cell which is then activated. However, one could argue that after entering the receiver cell, the signaling proteins spend a significant time diffusing in it before reaching the receptor gene construct. We discussed this at length before finally agreeing not to consider this as a delay.<br>
+
        <li id="disclaimer"><a href="/2011.igem.org:General_disclaimer" title="2011.igem.org:General disclaimer">Disclaimers</a></li>
-
Several models showed us that diffusion, even of a single signaling molecule, happen too fast to add any delay to the response time of the system. We will discuss two of these models here.<br>
+
    </ul>
-
==== Using diffusion equations ====
+
</div> <!-- close footer -->
-
The first step is to study the general principles of diffusion and to apply them to a single molecule. We expected to estimate the order of magnitude for diffusion time of molecules with this model, not to have a precise understanding of the movement of molecules in a cell. Most of our experimental designs rely on time measurement to characterize the nanotubes, it was therefore crucial to see if diffusion time could add a significant delay to the response of receiver cells.<br>
+
</div> <!-- close footer-wrapper -->
-
The principle of this model is quite simple. We use the statistical diffusion equation with a new normalization constant so that it describes the behaviour of one molecule. Rather than obtaining a concentration field, we end up with a distribution of the density of probability to find the molecule at a certain position and a certain time. We did not use any kind of boundary conditions, we therefore only model the "movements" of one molecule floating in an infinite water medium.<br>
+
</div>
-
The equation of diffusion is the following:<br>
+
<div id="scroll_left"><a href="https://2011.igem.org/Team:Paris_Bettencourt/Designs"><img src="https://static.igem.org/mediawiki/2011/0/0a/Arrow-left-big.png" style="width:100%;"></a><a href="https://2011.igem.org/Team:Paris_Bettencourt/Designs">Design overview</a></div>
-
<html><img src='https://static.igem.org/mediawiki/2011/4/4a/Diffusion_equation.png' style='width:230px;'/></html>
+
<div id="scroll_right"><a href="https://2011.igem.org/Team:Paris_Bettencourt/Experiments/List"><img src="https://static.igem.org/mediawiki/2011/e/e0/Arrow-right-big.png" style="width:100%;"></a><a href="https://2011.igem.org/Team:Paris_Bettencourt/Experiments/list">Experiments overview</a></div>
-
 
+
</html>
-
Where ''c'' is the concentration of particles in the cell, function of <html><img src='https://static.igem.org/mediawiki/2011/0/06/X_vector.png' style='width:20px;'/></html> (position) and ''t'' (time). ''D'' is the diffusion coefficient.<br>
+
-
 
+
-
We decided to use as initial conditions a Dirac function centered on the origin of space. The solution for such an equation is:<br>
+
-
 
+
-
<html><img src='https://static.igem.org/mediawiki/2011/0/04/Diffusion_solution_general.png' style='width:230px;'/></html>
+
-
 
+
-
 
+
-
 
+
-
This solution shows that through time, the density of probability "diffuses" in all directions equally. So we have take advantage of the spherical symetry of the problem.
+
-
 
+
-
<html><a href='https://static.igem.org/mediawiki/2011/3/3b/Diffusion_3D.png'><img src='https://static.igem.org/mediawiki/2011/3/3b/Diffusion_3D.png' style='width:100%;'/></a></html>
+
-
<html><a href='https://2011.igem.org/File:Diffusion_1D.png'><img src='https://static.igem.org/mediawiki/2011/0/0d/Diffusion_1D.png' style='width:100%;'/></a></html>
+
-
(Note that this serie of graph is here only as an example. The units are arbitrary.)
+
-
 
+
-
 
+
-
What is of interest for us is the instant when we can consider that the probability of presence of the molecule is "the same" for every point in the cell. The cell here is an arbitrary boundary which has no direct impact on the model. In our approach, the cell is a sphere centered on the origin and with a diameter of 1 micrometer (roughly Subtilis length).  We chose to consider that the probability is "the same" when the lowest probability within the cell is 95% of the highest probability in the cell (i.e. in the center).
+
-
 
+
-
<html><a href='https://2011.igem.org/File:95_percent_1D.png'><img src='https://static.igem.org/mediawiki/2011/f/f7/95_percent_1D.png' style='width:60%;'/></a></html>
+
-
 
+
-
We can then solve analytically the equation. Let us call (t0) the time when the probability is the same (at 95%) within the cell boundaries. This mean that for concentration for ''x=cell radius'' is 0.95 times concentration for ''x=0''.
+
-
 
+
-
<html><a href='https://2011.igem.org/File:95_percent_diffusion_time_equation.png'><img src='https://static.igem.org/mediawiki/2011/2/2b/95_percent_diffusion_time_equation.png' style='width:100%;'/></a></html>
+
-
 
+
-
We therefore see that the particle diffusing could be anywhere in the cell a mere XX s after entering the cell. This is too small to be mesured acurately through any biological reporter.
+
-
 
+
-
We therefore see that the particle diffusing could be anywhere in the cell a mere XX s after entering the cell. This is too small to be mesured accurately through any biological reporter.
+
-
 
+
-
==== Using a stochastic model ====
+
-
 
+
-
This model is very similar to the one used in [Ribosome kinetics and aa-tRNA competition determine rate and fidelity of peptide synthesis. Fluitt A, Pienaar E, Viljoen H. Comput Biol Chem. 2007 Oct;31(5-6):335-46. Epub 2007 Aug 15.]. We consider that the particle diffusing in the cell is a random walker.<br>
+
-
 
+
-
We want to see how long it takes for a particle of a given size to diffuse to any point of a cell. We use the following parameters:
+
-
*''V'' volume of the cell (10^-18 m^3)
+
-
*<html><img src='https://static.igem.org/mediawiki/2011/d/d2/Lambda.png' style='width:20px;'/></html> characteristic size of the particle (m)
+
-
*''D'' diffusion coefficient of the particle (m^2.s^-1)
+
-
 
+
-
We divide the cytoplasm volume ''V'' into <html><img src='https://static.igem.org/mediawiki/2011/7/72/Occupation_sites_general.png' style='width:80px;'/></html> occupation sites for the walker. The characteristic time <html><img src='https://static.igem.org/mediawiki/2011/8/84/Tau.png' style='width:20px;'/></html>  associated with the transition from one site to another is: <html><img src='https://static.igem.org/mediawiki/2011/8/81/Characteristic_time_general.png' style='width:90px;'/></html> [http://www.ece.gatech.edu/research/labs/bwn/papers/2011/c6.pdf]
+
-
 
+
-
If we have ''R'' walkers of this type, the probability that a molecule arrives at a given occupation site during the time interval <html><img src='https://static.igem.org/mediawiki/2011/8/84/Tau.png' style='width:20px;'/></html> is: <html><img src='https://static.igem.org/mediawiki/2011/b/bb/Probability_general.png' style='width:80px;'/></html>. In our case we study only one molecule (worst case scenario) so let us assume R=1.
+
-
 
+
-
The average time that elpases before the arrival of a particle is:<html><img src='https://static.igem.org/mediawiki/2011/9/9d/Characteristic_time_dev_general.png' style='width:190px;'/></html>.
+
-
 
+
-
You will find below a table of time before arrival for different molecules.
+
-
 
+
-
INSERT TABLE.
+
-
 
+
-
=== Parameters ===
+
-
Finally, we made some assumptions for certain parameters that are used in most our models. These assumptions are discussed at length in the Parameters section.
+
-
 
+
-
ADD LINK
+
-
 
+
-
== Direct observation ==
+
-
 
+
-
== Characterization ==
+
-
[[File:Parameters.png|thumb|center|upright=3.0|Relevant parameters for modeling]]
+
-
[[File:T7_equations.png|thumb|center|upright=3.0|Allosteric equations for modeling]]
+
-
 
+
-
=== T7 system ===
+
-
 
+
-
[[Team:Paris_Bettencourt/Modeling/T7_diffusion|The T7 diffusion model.]]
+
-
 
+
-
=== tRNA_amber system ===
+
-
[[Team:Paris_Bettencourt/Modeling/tRNA_diffusion|The tRNA diffusion model.]]
+
-
 
+
-
== Distribution of tRNA_amber in mRNA populations ==
+
-
[[Team:Paris_Bettencourt/Modeling/tRNA_distribution|The tRNA distribution model.]]
+
-
 
+
-
== Master/Slave ==
+
-
 
+
-
== Bi-directional communication ==
+
-
 
+
-
 
+
-
 
+
-
<br><br><br><br><br><br><br><br>
+

Latest revision as of 03:14, 29 October 2011

Team IGEM Paris 2011

Modeling

Modeling in our project

Our modeling was organized around two main questions:

  • Can we explain the transfer through nanotubes?
  • What will be the behaviour of our constructs and how will it impact our experimental designs?

Answering those two questions was essential for our project. We needed to know what to expect in order to design our experiments properly and to know what kind of restults we should obtain.

Investigating nanotube transfer

In order to answer the first question, we investigated the physical properties of cell membrane and passive diffusion to comprehend how the transfer could occur. We came up with two different ideas that could explain molecule transfer through the nanotubes, and based our original models on these assumptions, done in Java for passive diffusion and in Matlab for assisted diffusion.These two novel models show that transfer through the nanotubes, whether happening by passive diffusion or the so-called assisted diffusion, is happening too quickly to be accurately measured by fluorescent microscopy. As nanotube transfer is too fast compared to genetic response to allow us to measure its time span correctly, our conclusion was that our designs would not allow us to determine which one of these two processes (passive or assisted diffusion) is dominant during the transfer. Even though this makes it impossible to create a definitive model of molecule transfer through the nanotubes, the information provided by our two alternative models gave us an insight on the time scale of the transfer.

Passive diffusion in nanotubes We investigate here the hypothesis of passive diffusion through nanotubes.
Assisted diffusion We propose here a model explaining how diffusion through nanotubes could be "assisted" by the tension differential between cell walls.

Predicting the behaviour of our designs

The second question was a crucial step in the design of our experiments. Taking into account the predicted time scale and other information provided by our two general models, we were able to build models of each of our genetic networks. These models are an improvement of Uri Alon's approach in An Introduction to Systems Biology: Design Principles of Biological Circuits and were done in Matlab. With these models we showed that some designs might work better that the others. For instance, we prioritized the T7 RNA polymerase diffusion and tRNA diffusion systems and decided to concentrate our wet lab experiments and characterizatons in these systems. The ComS system, on the other hand, was less developed because of some disadvantages that our model predicted (high background even without induction, very high activation threshold mainly). Moreover, our models let us evaluate the response time of each of our our constructs. With these estimations, were able to prepare protocols for our microscopic experiments (by evaluating the characteristic response time of our system, the activation thresholds, etc.).

The basics about genetic networks modeling You can find here an introduction to our methods and the general idea behind most gene network models.
Our assumptions Because of the specificities of our project, we had to adapt the "classic" model to better represent our current situation. On top of that we made and justified a few other hypotheses detailed in this section.
Modeling our designs Models predicting the behaviour of our designs are detailled in this section.