Team:Paris/Production modeling

iGEM > Paris >  DryLab >  Vesicle production Model (delay model)

DryLab - Delay model: Improving message quality
 position: absolute; height: 23px; width: 30px; top: 0px; left: 120px; margin-top:10px; padding-top: 7px; background: url(http://2009.igem.org/wiki/images/1/1b/Left_menu_pari.png); z-index:4; }
 * 1) left-side {

height: 25px; width: 400px; position: absolute; top: 0px; left: 130px; margin-top:10px; padding-top: 5px; background: #dadada; z-index:5; }
 * 1) middle-side {

position: absolute; height: 23px; width: 30px; margin-top:10px; padding-top: 7px; top: 0px; left: 520px; background: url(http://2009.igem.org/wiki/images/4/40/Right_menu_paris.png); z-index:4; }   position: absolute; height: 23px; width: 30px; top: 35px; left: 110px; margin-top:10px; padding-top: 7px; background: url(http://2009.igem.org/wiki/images/4/40/Left_menu_paris2.png); z-index:4; }
 * 1) right-side {
 * 1) left-side2 {

height: 25px; width: 430px; position: absolute; top: 35px; left: 120px; margin-top:10px; padding-top: 5px; background: white; z-index:5; }
 * 1) middle-side2 {

position: absolute; height: 23px; width: 30px; margin-top:10px; padding-top: 7px; top: 35px; left: 530px; background: url(http://2009.igem.org/wiki/images/f/f8/Right_menu_paris2.png); z-index:4; }
 * 1) right-side2 {

a.menu_sub { padding-left: 7px; padding-right: 7px; }

a.menu_sub_active { padding-left: 7px; padding-right: 7px; color:#b0310e; font-weight:bold; }  Main |  Delay model|  Vesicle model|  Fec simulation  Introduction |  Transcriptional Cascade|  The feed-forward motif

This part aimed at solving the first problem linked to the creation of messengers :

To improve the quality of signal sent, how can we get a good synchronization between the maximal vesicles production rate and the maximal concentration of proteins to encapsulate

<span/ id="1"><span/ id="2"> Obviously, it is crucial to send vesicles only when the concentration of proteins to encapsulate is at its maximum level in the periplasm. To optimize the quality of our message, we tried to make this peak of vesicles creation corresponding to maximum levels of protein concentration. We assumed that the number of vesicles produced is at its maximum when the concentration of TolRII has reached its peak.

In what follows, we consider two different designs for implementing the delay: using a transcriptional cascade and using feed-forward motifs. We found that the design based on transcriptional cascade is sufficient for obtaining the desired delay[1], whereas the design based on feed-forward motifs allows for more flexibility, at the cost of an increased complexity[2]. These studies are described below.

Implementing the delay by a transcriptional cascade
<span/ id="3"><span/ id="4">

We have decided to construct the following system. The addition of arabinose to the medium leads to the activation of the pBad promoters: two proteins are produced: the protein to encapsulate and LacI-LVA[3][4]. The latter inhibits the expression of TetR and thus releases the inhibition of TolRII, responsible for vesicle formation.



<span/ id="5"> To study this system, we simply used differential equations. We used the method described by U. Alon in his book An Introduction to Systems Biology, with notably Hill functions to describe the regulation of the promoters[5].

The network behavior can be described by a system of differential equations:

<div style="width: 400px; text-align: justify; margin-left:200px;"> These equations describe the evolution of the concentrations of the messenger protein P, LacI repressor, TetR repressor, and the TolRII protein responsible for vesiculation. The gammas are the degradation rates of the respective compounds, while the sigmoidal functions describe the non-linear production rates. The parameters are estimated according to typical values as discussed below

As a first approximation, we assumed that :
 * When the pBad promoter is induced, the concentration of arabinose in the medium is very high and constant during he whole study; as a consequence, we will consider that the creation rate of Protein P and LacI* is constant during the experiment:

<span/ id="6">
 * We considered that all the binding constants are identical[6] and of an average of 40nM which correspond to approximately 40 monomers per cell; we can write that:


 * We chose identical intrinsic promoter activities, &beta;, all equal to 4000 proteins/cell cycle:


 * All the time units were expressed in units of cell cycle (approximately half an hour). As a consequence we chose a dilution rate &gamma; of 1 for protein without special tags. For the LacI protein with a degradation tag (LacI-LVA), the dilution rate is 3 times higher[3]:

For simulation, we used the following initial quantities:
 * 4000 units of TetR (steady state level at inactivation)
 * no protein to be encapsulated
 * no LacI
 * no TolRII

We got the following dynamics using Matlab software, showing that this simple system is capable of generating the desired behavior.



Plot showing that vesiculation happens only once the protein encapsulated has reached its maximal concentration.

Implementing the delay by feed-forward motifs
Using the previous design, we were able to find a way to get the maximal vesicle formation (ie maximal TolRII quantity) when the concentration of proteins to encapsulate is at its maximum. This system satisfies our requirements. However, it is not very flexible. In particular, the concentration of the protein to encapsulate remains always high after induction. This could be detrimental for the cell if this protein is toxic. Consequently, we considered an alternative design allowing for more flexibility, in which the expression of the protein to encapsulate can be transient.

The Feed Forward Motif
<span/ id="7">

The feed forward motif was studied and described by U. Alon in his book "An Introduction to Synthetic Biology" where he describes the different forms of existing feed forward; they can be either coherent or incoherent[7].

In each case, there exist a different network organization to create a feed-forward loop as shown on the picture



The X and Y inputs can be combined either using a logical AND or OR gate thus giving different dynamics properties in each case.

The coherent feed forward loop allow the creation of a delay either on activation or on deactivation depending on the logical function between the two input X and Y. When the output logical gate is an AND function, the delay appears on activation, while when the output gate is an OR function, the delay appears on deactivation. An example of characterized feed-forward loop is the fliA system which was used by the | 2008 Paris team in order to create a oscillator.

In an incoherent feed forward loop, the two inputs of the last gate function have two different logical levels; one molecule (molecule X for example) is an activator while the other molecule is a repressor of the promoter controlling the expression of the output. This way, when the input signal is on, the expression of molecule X starts activating both the promoters controlling the expression of Y and Z ; transcription starts and the amount of Y and Z molecules increases. Nevertheless, when the concentration of Y has reached its repression threshold, repression starts and the concentration of Z molecule decrease to its steady state.



Two elementary incoherent feed forward motifs can be combined in order to create two time delayed bursts as described by U.Alon in its description of the development of Bacillus subtilis spore. Anyway, it seems quite difficult to find a molecule which would be both an activator and a repressor; we decided to try and think about another combination of these two feed forwards so as to get the dynamics we wanted as shown on the following part.

Our design proposition
The two feed-forward loops used are incoherent and organized in a way that each molecule is either an activator or a repressor (and not both). The following scheme gives the functioning of the system.



To study the dynamics of our model, we decided to use a classical approach using Michaelis-Menten equations to model the activation and repression of the various promoters, including a cooperativity of order 2 for the repressor which often the case for 'classical' repressors such as LacI. We obtained the following system of equations:



It is possible to obtain a broad spectrum of different dynamical behaviors by playing with the value of the maximum production rate of Z3, called here &#x03b2;. This maximal production rate impacts deeply both the evolution of the exported protein concentration and the evolution of TolRII concentration. This rate is initially chosen equal to 40 and 400 molecules by cell cycle and we have plotted the different responses for multiples of these initial values.

In each case, the initial conditions are the same :
 * no proteins encapsulated
 * no TolRII
 * no Z1 molecules
 * no Z2 molecules
 * no Z3 molecules
 * a weak amount of Z4 molecules (around a hundred)

Weak maximal Z3 production rate An example of dynamics for a range of weak &#x03b2; values; for each plot, the maximal Z3 production rate is given by: &#x03b2;=k*40

Strong maximal Z3 production rate An example of dynamics for a range of strong &#x03b2; values; for each plot, the maximal Z3 production rate is given by: &#x03b2;=k*400

A more flexible design
Regardless of the steady state concentrations of the molecules, the use of two incoherent feed-forward loops allows for more flexibility in the behaviors of our system. With a simple delay system, we just get one possible dynamic behavior whereas using a variable range of production rates, we can obtain two behaviors:


 * When &#x03b2; is weak (&#x03b2; = 40), the number of vesicles (ie of TolRII) produced slowly increase before returning to zero, whereas the protein to be exported inside the vesicles is created very quickly and remains at a high concentration compared to the TolRII protein. We can therefore send vesicles containing an important amount of interesting proteins and, after a while, stop expressing dangerous proteins for the cell. All of this can be performed with a single input signal.


 * When &#x03b2; is stronger, the behavior is different and we notice a change of concavity of the TolRII curve when the amount of protein exported is at its maximum. During such a change of concavity, the derivative reaches its maximum - in other words, the amount of TolRII produced by time units is at its highest level. Considering the hypothesis according to which the production rate of vesicles is proportional to TolRII production rate, we get a peak in vesicle production when the concentration of the proteins to exported is at its maximum.

Unfortunately, we did not have the time to try to construct this system in the lab and only the first design was implemented. Nevertheless, we thought that this analysis was very useful to get a good insight of the diversity of the possible behaviors of the sender cells.

<div id="paris_content_boxtop"> <div id="paris_content">