"
Team:Paris Bettencourt/what is modeling
From 2011.igem.org
| Line 5: | Line 5: | ||
<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 the preparation of 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.</p> | <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 the preparation of 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.</p> | ||
| - | <p>The bible for synthetic biology modeling is Uri Alon's book < | + | <p>The bible for synthetic biology modeling is Uri Alon's book <i>An Introduction to Systems Biology: Design Principles of Biological Circuits</i> <a href="#references">[1]</a>. 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.</p><br> |
| - | <p>The gene ''geneX'' is responsible for the production of the corresponding protein < | + | <p>The gene ''geneX'' is responsible for the production of the corresponding protein <i>X</i>. The promoter <i>pX</i> controlling the expression <i>geneX</i> can be:</p> |
<ul> | <ul> | ||
| - | <li>A constitutive promoter, < | + | <li>A constitutive promoter, <i>geneX</i> is activated whatever the conditions</li> |
| - | <li>Positively regulated, < | + | <li>Positively regulated, <i>geneX</i> is more active when an inducer is present</li> |
| - | <li>Negatively regulated, < | + | <li>Negatively regulated, <i>geneX</i> is less active when an repressor is present</li> |
</ul> | </ul> | ||
| Line 24: | Line 24: | ||
<br> | <br> | ||
| - | <p>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 < | + | <p>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 <i>pX</i> 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.</p> |
<br> | <br> | ||
<p>We will know take a look at the parameters involved in modeling this network.</p> | <p>We will know take a look at the parameters involved in modeling this network.</p> | ||
<ul> | <ul> | ||
| - | <li>< | + | <li><i>[X]</i> is the concentration of protein <i>X</i></li> |
| - | <li><img src='https://static.igem.org/mediawiki/2011/a/a8/Alpha.png' style='width:20px;'/> is the expression rate of protein < | + | <li><img src='https://static.igem.org/mediawiki/2011/a/a8/Alpha.png' style='width:20px;'/> is the expression rate of protein <i>X</i> (molecule.s<sup>-1</sup>). It mainly depends on the constitutive promoter.</li> |
<li><img src='https://static.igem.org/mediawiki/2011/0/05/Delta_dil.png' style='width:24px;'/> is the dilution rate, due to cell division (s<sup>-1</sup>)</li> | <li><img src='https://static.igem.org/mediawiki/2011/0/05/Delta_dil.png' style='width:24px;'/> is the dilution rate, due to cell division (s<sup>-1</sup>)</li> | ||
| - | <li><img src='https://static.igem.org/mediawiki/2011/6/62/Delta_deg.png' style='width:24px;'/> is the degradation rate of protein < | + | <li><img src='https://static.igem.org/mediawiki/2011/6/62/Delta_deg.png' style='width:24px;'/> is the degradation rate of protein <i>X</i> (s<sup>-1</sup>)</li> |
</ul> | </ul> | ||
<br> | <br> | ||
| Line 42: | Line 42: | ||
<br> | <br> | ||
| - | <p>Now, let's assume that < | + | <p>Now, let's assume that <i>pX</i> is auto-regulated, either positively or negatively. We need to introduce a few new parameters.</p> |
<ul> | <ul> | ||
| - | <li>< | + | <li><i>K</i> is the dissociation constant, representing the binding of a inducer/repressor to the promoter</li> |
| - | <li>< | + | <li><i>n</i> is the Hill coefficient of the function</li> |
<li><img src='https://static.igem.org/mediawiki/2011/3/3f/Beta.png' style='width:20px;'/> is the maximum production rate of our gene (molecule.s<sup>-1</sup>)</li> | <li><img src='https://static.igem.org/mediawiki/2011/3/3f/Beta.png' style='width:20px;'/> is the maximum production rate of our gene (molecule.s<sup>-1</sup>)</li> | ||
<li><img src='https://static.igem.org/mediawiki/2011/a/a8/Alpha.png' style='width:20px;'/> now represents the basal expression ("leaking") of the gene (molecule.s<sup>-1</sup>)</li> | <li><img src='https://static.igem.org/mediawiki/2011/a/a8/Alpha.png' style='width:20px;'/> now represents the basal expression ("leaking") of the gene (molecule.s<sup>-1</sup>)</li> | ||
Revision as of 08:26, 1 September 2011
What is modeling in synthetic biology?
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 the preparation of 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.
The bible for synthetic biology modeling is Uri Alon's book An Introduction to Systems Biology: Design Principles of Biological Circuits [1]. 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.
The gene ''geneX'' is responsible for the production of the corresponding protein X. The promoter pX controlling the expression geneX can be:
- A constitutive promoter, geneX is activated whatever the conditions
- Positively regulated, geneX is more active when an inducer is present
- Negatively regulated, geneX is less active when an repressor is present
Constitutive promoter
Positive autoregulation
Negative autoregulation
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.
We will know take a look at the parameters involved in modeling this network.
- [X] is the concentration of protein X
is the expression rate of protein X (molecule.s-1). It mainly depends on the constitutive promoter.
is the dilution rate, due to cell division (s-1)
is the degradation rate of protein X (s-1)
The equation and solution modeling the behaviour of this system are the following:
Constitutive expression of protein X over time
Now, let's assume that pX is auto-regulated, either positively or negatively. We need to introduce a few new parameters.
- K is the dissociation constant, representing the binding of a inducer/repressor to the promoter
- n is the Hill coefficient of the function
is the maximum production rate of our gene (molecule.s-1)- now represents the basal expression ("leaking") of the gene (molecule.s-1)
The constants used in our examples below are arbitrary and chosen to amplify the behaviours expected. The equations and solutions are now:
Positive autoregulation equation and behaviour
Negative autoregulation equation and behaviour
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.
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.
