Team:BIOTEC Dresden/Modeling
From 2010.igem.org
(41 intermediate revisions not shown) | |||
Line 2: | Line 2: | ||
<html> | <html> | ||
+ | <head> | ||
+ | <style type="text/css"> | ||
+ | .lin { display: none; } | ||
+ | .den, .nom { display: block; width:100%; text-align:center } | ||
+ | .nom { text-decoration: underline; } | ||
+ | .dwn{ bottom:-5px; left:-25px; position:relative; } | ||
+ | </style> | ||
+ | <script type="text/javascript" src="http://endo.bio.tu-dresden.de/igem/scripthost/fancybox/jquery.easing-1.3.pack.js"></script> | ||
+ | <script type="text/javascript" src="http://endo.bio.tu-dresden.de/igem/scripthost/fancybox/jquery.fancybox-1.3.2.pack.js"></script> | ||
+ | <link rel="stylesheet" href="http://endo.bio.tu-dresden.de/igem/scripthost/fancybox/jquery.fancybox-1.3.2.css" type="text/css" media="screen" /> | ||
+ | <script type="text/javascript"> | ||
+ | $(document).ready(function() { | ||
+ | $("a[rel='lightbox']").fancybox({ | ||
+ | 'transitionIn' : 'elastic', | ||
+ | 'transitionOut' : 'elastic', | ||
+ | 'speedIn' : 600, | ||
+ | 'speedOut' : 200, | ||
+ | 'overlayShow' : false, | ||
+ | 'titleShow': false | ||
+ | }); | ||
+ | }); | ||
+ | </script> | ||
+ | |||
+ | </head> | ||
Line 12: | Line 36: | ||
<p> (3) All promoters have the same activity resulting in a constant transcription rate. </p> | <p> (3) All promoters have the same activity resulting in a constant transcription rate. </p> | ||
<p>(4 ) The rate of translation is considered to be equal yet the degradation rate of proteins differs. </p> | <p>(4 ) The rate of translation is considered to be equal yet the degradation rate of proteins differs. </p> | ||
- | <p>(5 )AHL transport is only restricted by diffusion through the cell membrane following the concentration gradient, AHL<sub>out</sub> > AHL<sub>in<sub>.<p> | + | <p>(5 )AHL transport is only restricted by diffusion through the cell membrane following the concentration gradient, AHL<sub>out</sub> > AHL<sub>in</sub>.<p> |
<p>(6) Binding and dissociation constants of AHL and luxR are based on mass-action kinetics.</p> | <p>(6) Binding and dissociation constants of AHL and luxR are based on mass-action kinetics.</p> | ||
<p>(7) The enzymatic activity of LuxI within the fusion protein is not decreased by the fusion event.</p> | <p>(7) The enzymatic activity of LuxI within the fusion protein is not decreased by the fusion event.</p> | ||
Line 23: | Line 47: | ||
</tr> | </tr> | ||
<tr class="even"> | <tr class="even"> | ||
- | <td><b> | + | <td><b>Growth Rate</b></td> |
<td>µ</td> | <td>µ</td> | ||
<td>1.3 h<sup>-1</sup> </td> | <td>1.3 h<sup>-1</sup> </td> | ||
Line 37: | Line 61: | ||
<td>l</td> | <td>l</td> | ||
<td>25 </td> | <td>25 </td> | ||
- | |||
</tr> | </tr> | ||
<tr class="odd"> | <tr class="odd"> | ||
Line 93: | Line 116: | ||
<td>K<sub>m</sub></td> | <td>K<sub>m</sub></td> | ||
<td>130 µM</td> | <td>130 µM</td> | ||
+ | </tr> | ||
+ | <tr class="odd"> | ||
+ | <th colspan="3">Diffusion Coefficients</th> | ||
+ | </tr> | ||
+ | <tr class="even"> | ||
+ | <td><b>AHL diffusion coefficient</b></td> | ||
+ | <td>k<sub>dif</sub></td> | ||
+ | <td>0.001 mm<sup>2</sup>/min</td> | ||
</tr> | </tr> | ||
</table> | </table> | ||
Line 98: | Line 129: | ||
<h3>Simplified Model</h3> | <h3>Simplified Model</h3> | ||
<p>In a first simple model a single cell and its kinetics in output will be analyzed. The major input in this simplified scenario is AHL produced from SAM catalyzed by LuxI. The direct output is GFP or any other fluorescent reporter like RFP and YFP. In the gene design it can be seen that GFP expression is coupled to its induction upon binding of LuxR-AHL complex to the luxR promoter. </p> | <p>In a first simple model a single cell and its kinetics in output will be analyzed. The major input in this simplified scenario is AHL produced from SAM catalyzed by LuxI. The direct output is GFP or any other fluorescent reporter like RFP and YFP. In the gene design it can be seen that GFP expression is coupled to its induction upon binding of LuxR-AHL complex to the luxR promoter. </p> | ||
- | + | ||
+ | <a alt="Complete Gene Circuit of SensorBricks" rel="lightbox" href="https://static.igem.org/mediawiki/2010/0/07/BiotecDresden_Fusion_Construct.png"> <img id="Fusion_Construct" class="border thumb" src="https://static.igem.org/mediawiki/2010/0/07/BiotecDresden_Fusion_Construct.png" style="margin-left: 30%;"></a> | ||
<p>In the extracellular environment the enzyme LuxI is present in form of the fusion construct, proteinA-LuxI. The conversion of S-adenosyl methionine (SAM) to acylated homoserine lactone (AHL) catalyzed by LuxI follows enzyme kinetics by Michaelis-Menten:</p> | <p>In the extracellular environment the enzyme LuxI is present in form of the fusion construct, proteinA-LuxI. The conversion of S-adenosyl methionine (SAM) to acylated homoserine lactone (AHL) catalyzed by LuxI follows enzyme kinetics by Michaelis-Menten:</p> | ||
- | <span class="nom"><span class=" | + | <span class="nom"><span class="dwn">v = </span><span><i>V<sub>max</sub> x SAM<sub>ex</sub></i></span><span class="lin">/</span></span> |
- | class=" | + | <span class="den"><span class="lin">(</span><i>K<sub>m</sub></i> + <i>SAM<sub>ex</sub></i><span class="lin">)</span></span></span> |
- | <span class=" | + | |
<p>Next, AHL has to be transported across the membrane of the cell via diffusion. The inside concentration is therefore determined by the diffusion rate across the membrane (k<sub>dif</sub>), the initial concentration (AHL<sub>ex</sub>) as well as the decay constant of AHL inside the cell (k<sub>AHL</sub>). </p> | <p>Next, AHL has to be transported across the membrane of the cell via diffusion. The inside concentration is therefore determined by the diffusion rate across the membrane (k<sub>dif</sub>), the initial concentration (AHL<sub>ex</sub>) as well as the decay constant of AHL inside the cell (k<sub>AHL</sub>). </p> | ||
- | + | ||
- | <p>Furthermore AHL as well as constitutively expressed LuxR form a complex which can bind to the luxR promoter. The amount of AHL-LuxR complex binding the luxR promoter is dependent on the concentration of AHL, the concentration of LuxR and its degradation (k<sub>LuxR</sub> | + | |
- | + | <span class="nom" style="width: 90%;"><i>d[AHL]</i></span> | |
+ | <span class="den" style="width: 90%;"><i>dt</i></span><span style="position: relative; left: 48%; bottom: 33px;">= k<sub>diff</sub>(AHL<sub>ex</sub>- AHL<sub>in</sub>) -(k<sub>dec</sub> * AHL<sub>in</sub>)</span> | ||
+ | |||
+ | <p>Furthermore AHL as well as constitutively expressed LuxR form a complex which can bind to the luxR promoter. The amount of AHL-LuxR complex binding the luxR promoter is dependent on the concentration of AHL, the concentration of LuxR and its degradation (k<sub>LuxR</sub>) the binding and dissociation constant of AHL and LuxR (k<sub>dis</sub>=. </p> | ||
+ | <span style="position: relative; left: 42%;">[LuxR-AHL]=k<sub>dis</sub> x [AHL] x [LuxR]</span> | ||
<p>A signal can be detected using different reporters like fluorescent proteins, pigments or substances that can be detected. We focused on the detection using fluorescent proteins like GFP which could be detected using absorbance measurements. In this easy model the kinetics of GFP as output signal will be displayed. </p> | <p>A signal can be detected using different reporters like fluorescent proteins, pigments or substances that can be detected. We focused on the detection using fluorescent proteins like GFP which could be detected using absorbance measurements. In this easy model the kinetics of GFP as output signal will be displayed. </p> | ||
- | + | ||
+ | <span class="nom"><i>G(t)</i></span> | ||
+ | <span class="den"><i>dt</i></span><span style="position: relative; left: 55%; bottom: 33px;">= x(t)-γG(t)</span> | ||
<p>The formula displays GFP expression over time G(t). The term X(t) refers to the mRNA production and degradation as well as fluorescent read out obtained in RFU, composed of kt x M, where kt is an unknown constant and M is mRNA production remaining constant for all proteins expressed in our model. The Term γ stands for dilution rate which can also be represented by the growth rate µ of the E.coli Dh5α. The half life of GFP was experimentally determined to be around 24 hours. Therefore the value of γ is about 0.029 h<sup>-1</sup> and can be neglected as it is so much smaller than dilution rate. </p> | <p>The formula displays GFP expression over time G(t). The term X(t) refers to the mRNA production and degradation as well as fluorescent read out obtained in RFU, composed of kt x M, where kt is an unknown constant and M is mRNA production remaining constant for all proteins expressed in our model. The Term γ stands for dilution rate which can also be represented by the growth rate µ of the E.coli Dh5α. The half life of GFP was experimentally determined to be around 24 hours. Therefore the value of γ is about 0.029 h<sup>-1</sup> and can be neglected as it is so much smaller than dilution rate. </p> | ||
<h3>Results</h3> | <h3>Results</h3> | ||
- | <p>For the simulation first the model of a growing culture with direct AHL input and GFP output was analyzed. | + | |
+ | <p>For the simulation first the model of a growing culture with direct AHL input and GFP output was analyzed.</p> | ||
+ | |||
+ | <a href="https://static.igem.org/mediawiki/2010/b/bb/BiotecDresden_Simple_GFP_reporter.png" rel="lightbox" title="Simple Model"><img class="border thumb" src="https://static.igem.org/mediawiki/2010/b/bb/BiotecDresden_Simple_GFP_reporter.png" class="border thumb" style="display: inline; margin-left: 11%; height: 250px;"></a> | ||
+ | |||
+ | <a href="https://static.igem.org/mediawiki/2010/8/89/BiotecDresden_AHL_diffusion.jpg" rel="lightbox" title="Simple Model"><img class="border thumb" src="https://static.igem.org/mediawiki/2010/8/89/BiotecDresden_AHL_diffusion.jpg" class="border thumb" style="margin-left: 1%; height: 250px;"></a> | ||
+ | <p>The three graphs shown below demonstrate the behavior of the gene construct upon AHL induction. The higher the concentration of AHL used for the induction the higher is also the fluorescence signal detected from GFP. AHL diffuses basically without any constriction into the bacteria clearly depicted by the exponentially following decay of the graphs (Figure 1). Correspondingly the endpoint fluorescent signals of GFP vary depending on the AHL input which can be seen in Figure 2. For high AHL concentrations the endpoint GFP signal is elevated compared to lower AHL concentrations. </p> | ||
+ | |||
+ | |||
+ | <a href="https://static.igem.org/mediawiki/2010/e/e1/BiotecDresden_GFP_output.jpg" rel="lightbox" title="Simple Model"><img class="border thumb" src="https://static.igem.org/mediawiki/2010/e/e1/BiotecDresden_GFP_output.jpg" class="border thumb" style="display: inline; margin-left: 11%; height: 250px;"></a> | ||
+ | |||
+ | <a href="https://static.igem.org/mediawiki/2010/1/15/LuxR-AHL_complex.jpg" rel="lightbox" title="Simple Model"><img class="border thumb" src="https://static.igem.org/mediawiki/2010/1/15/LuxR-AHL_complex.jpg" class="border thumb" style="margin-left: 1%; height: 250px;"></a> | ||
+ | |||
<p>In Figure 3 the LuxR-AHL complex formation is displayed which also increases in endpoint value for increased AHL concentrations caused by the availability of binding partner. </p> | <p>In Figure 3 the LuxR-AHL complex formation is displayed which also increases in endpoint value for increased AHL concentrations caused by the availability of binding partner. </p> | ||
Line 121: | Line 170: | ||
<h3>Discussion</h3> | <h3>Discussion</h3> | ||
<p>Even this simplified model shows that the output GFP signal can be used to quantify the amount of AHL present in the sample. This will allow us in the future to accurately determine the amount of cancer cells present in the patient’s blood sample. Using amplification constructs like introducing LuxI downstream of the Lux pR together with a reporter protein even more AHL will be produced. The set-up can be seen in Figure 4. </p> | <p>Even this simplified model shows that the output GFP signal can be used to quantify the amount of AHL present in the sample. This will allow us in the future to accurately determine the amount of cancer cells present in the patient’s blood sample. Using amplification constructs like introducing LuxI downstream of the Lux pR together with a reporter protein even more AHL will be produced. The set-up can be seen in Figure 4. </p> | ||
- | |||
</div> | </div> | ||
- | |||
- | |||
- | |||
- | |||
- | |||
</body> | </body> | ||
</html> | </html> | ||
- | [[Category:BIOTEC_Dresden/Project| | + | [[Category:BIOTEC_Dresden/Project|Modeling]] |
{{Biotec_Dresden/Bottom}} | {{Biotec_Dresden/Bottom}} |
Latest revision as of 02:41, 28 October 2010
Assumptions
(1) The cell growth of E.coli DH5α follow logistic growth kinetics with an assigned growth constant µ.
(2) For simplicity decay of all products follows first-order kinetics.
(3) All promoters have the same activity resulting in a constant transcription rate.
(4 ) The rate of translation is considered to be equal yet the degradation rate of proteins differs.
(5 )AHL transport is only restricted by diffusion through the cell membrane following the concentration gradient, AHLout > AHLin.
(6) Binding and dissociation constants of AHL and luxR are based on mass-action kinetics.
(7) The enzymatic activity of LuxI within the fusion protein is not decreased by the fusion event.
Parameters
General Parameters | ||
---|---|---|
Growth Rate | µ | 1.3 h-1 |
Promoter Activity | c | 0.01 mM/h |
Plasmid copy number | l | 25 |
Degradation Rates | ||
AHL decay rate | kAHL | 0.11 h-1 |
LuxR degradation rate | kLuxR | 0.11 h-1 |
GFP degradation rate | kGFP | 6.3-3min-1 |
Dissociation Constants | ||
AHL-LuxR dissociation constant | kdis | 0.09 - 1 µM |
Hill Cooperativity | ||
LuxR cooperativity | nLuxR | 2 |
AHL-LuxR cooperativity | nAHL-LuxR | 1 |
Enzyme Kinetic Constants | ||
Maximal Velocity LuxI | Vmax | 1.1 mol/min |
Km SAM | Km | 130 µM |
Diffusion Coefficients | ||
AHL diffusion coefficient | kdif | 0.001 mm2/min |
Simplified Model
In a first simple model a single cell and its kinetics in output will be analyzed. The major input in this simplified scenario is AHL produced from SAM catalyzed by LuxI. The direct output is GFP or any other fluorescent reporter like RFP and YFP. In the gene design it can be seen that GFP expression is coupled to its induction upon binding of LuxR-AHL complex to the luxR promoter.
In the extracellular environment the enzyme LuxI is present in form of the fusion construct, proteinA-LuxI. The conversion of S-adenosyl methionine (SAM) to acylated homoserine lactone (AHL) catalyzed by LuxI follows enzyme kinetics by Michaelis-Menten:
v = Vmax x SAMex/ (Km + SAMex)Next, AHL has to be transported across the membrane of the cell via diffusion. The inside concentration is therefore determined by the diffusion rate across the membrane (kdif), the initial concentration (AHLex) as well as the decay constant of AHL inside the cell (kAHL).
d[AHL] dt= kdiff(AHLex- AHLin) -(kdec * AHLin)Furthermore AHL as well as constitutively expressed LuxR form a complex which can bind to the luxR promoter. The amount of AHL-LuxR complex binding the luxR promoter is dependent on the concentration of AHL, the concentration of LuxR and its degradation (kLuxR) the binding and dissociation constant of AHL and LuxR (kdis=.
[LuxR-AHL]=kdis x [AHL] x [LuxR]A signal can be detected using different reporters like fluorescent proteins, pigments or substances that can be detected. We focused on the detection using fluorescent proteins like GFP which could be detected using absorbance measurements. In this easy model the kinetics of GFP as output signal will be displayed.
G(t) dt= x(t)-γG(t)The formula displays GFP expression over time G(t). The term X(t) refers to the mRNA production and degradation as well as fluorescent read out obtained in RFU, composed of kt x M, where kt is an unknown constant and M is mRNA production remaining constant for all proteins expressed in our model. The Term γ stands for dilution rate which can also be represented by the growth rate µ of the E.coli Dh5α. The half life of GFP was experimentally determined to be around 24 hours. Therefore the value of γ is about 0.029 h-1 and can be neglected as it is so much smaller than dilution rate.
Results
For the simulation first the model of a growing culture with direct AHL input and GFP output was analyzed.
The three graphs shown below demonstrate the behavior of the gene construct upon AHL induction. The higher the concentration of AHL used for the induction the higher is also the fluorescence signal detected from GFP. AHL diffuses basically without any constriction into the bacteria clearly depicted by the exponentially following decay of the graphs (Figure 1). Correspondingly the endpoint fluorescent signals of GFP vary depending on the AHL input which can be seen in Figure 2. For high AHL concentrations the endpoint GFP signal is elevated compared to lower AHL concentrations.
In Figure 3 the LuxR-AHL complex formation is displayed which also increases in endpoint value for increased AHL concentrations caused by the availability of binding partner.
Discussion
Even this simplified model shows that the output GFP signal can be used to quantify the amount of AHL present in the sample. This will allow us in the future to accurately determine the amount of cancer cells present in the patient’s blood sample. Using amplification constructs like introducing LuxI downstream of the Lux pR together with a reporter protein even more AHL will be produced. The set-up can be seen in Figure 4.