Team:BIOTEC Dresden/Modeling

From 2010.igem.org

(Difference between revisions)
 
(53 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 7: Line 31:
<div id="content_prim">
<div id="content_prim">
   
   
-
                <p>Detection of cell surface cancer markers is a key diagnostic step during cancer therapy as it allows the efficiency of a therapy to be determined. Current antibody-based flow cytometric detection methods are the gold standard. For our 2010 iGEM project, we hope to develop SensorBricks as a reliable and modular system for antigen recognition, signal amplification and quantification. We want to push the boundaries of detection in order to allow earlier diagnosis and thus improve therapeutic prognosis in cancer therapy. Initial steps of SensorBricks will focus on the detection of CD33 and other leukemic markers to increase diagnostic stringency.</p>
+
<h3>Assumptions</h3>
-
                <p>There are three major components in SensorBricks: (i) monoclonal antibodies that bind to an antigen of interest, (ii) a LuxI-Protein A fusion construct which non-specifically binds antibodies and produces the autoinducer N-Acyl homoserine lactone (AHL), and (iii) a Escherichia coli based biosensor which in the presence of AHL strongly amplifies the production of green fluorescence protein (GFP). By coupling signal detection to a genetic circuit, we will be able to amplify the signal in a quantifiable manner, allowing the identification of cancer markers expressed in minute quantities. The presence of CD33+ cells in the blood of a patient will trigger the production of AHL, which in turn, activates the production of GFP in our E. coli biosensor. As such, we should be able to correlate the concentration of CD33 as a function of GFP fluorescence. </p>
+
<p> (1) The cell growth of E.coli DH5α follow logistic growth kinetics with an assigned growth constant µ.<p>
-
                <p>The modularity of SensorBricks is a key feature as it allows the same system and protocols to be used in combination with any antibody of choice to detect any antigen of interest.</p>
+
<p>(2) For simplicity decay of all products follows first-order kinetics. </p>
-
            </div>
+
<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>(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>(7) The enzymatic activity of LuxI within the fusion protein is not decreased by the fusion event.</p>
 +
<h3>Parameters</h3>
-
<div class="visualClear"></div>
+
<table class="tformat" cellspacing="0" cellpadding="0">
 +
<tr class="odd">
 +
<th colspan="3">General Parameters</th>
 +
</tr>
 +
<tr class="even">
 +
<td><b>Growth Rate</b></td>
 +
<td>µ</td>
 +
<td>1.3 h<sup>-1</sup> </td>
 +
</tr>
 +
<tr class="odd">
 +
<td><b>Promoter Activity</b></td>
 +
<td>c</td>
 +
<td>0.01 mM/h</td>
 +
</tr>
 +
</tr>
 +
<tr class="even">
 +
<td><b>Plasmid copy number</b></td>
 +
<td>l</td>
 +
<td>25 </td>
 +
</tr>
 +
<tr class="odd">
 +
<th colspan="3">Degradation Rates</th>
 +
</tr>
 +
<tr class="even">
 +
<td><b>AHL decay rate</b></td>
 +
<td>k<sub>AHL</sub></td>
 +
<td>0.11 h<sup>-1</sup></td>
 +
</tr>
 +
<tr class="odd">
 +
<td><b>LuxR degradation rate</b></td>
 +
<td>k<sub>LuxR</sub></td>
 +
<td>0.11 h<sup>-1</sup></td>
 +
</tr>
 +
</tr>
 +
<tr class="even">
 +
<td><b>GFP degradation rate</b></td>
 +
<td>k<sub>GFP</sub></td>
 +
<td>6.3<sup>-3</sup>min<sup>-1</sup></td>
 +
</tr>
 +
<tr class="odd">
 +
<th colspan="3">Dissociation Constants</th>
 +
</tr>
 +
<tr class="even">
 +
<td><b>AHL-LuxR dissociation constant</b></td>
 +
<td>k<sub>dis</sub></td>
 +
<td>0.09 - 1 µM</td>
 +
</tr>
 +
<tr class="odd">
 +
<th colspan="3">Hill Cooperativity</th>
 +
</tr>
 +
<tr class="even">
 +
<td><b>LuxR cooperativity</b></td>
 +
<td>n<sub>LuxR</sub></td>
 +
<td>2</td>
 +
</tr>
 +
</tr>
 +
<tr class="odd">
 +
<td><b>AHL-LuxR cooperativity</b></td>
 +
<td>n<sub>AHL-LuxR</sub></td>
 +
<td>1</td>
 +
</tr>
 +
<tr class="odd">
 +
<th colspan="3">Enzyme Kinetic Constants</th>
 +
</tr>
 +
<tr class="even">
 +
<td><b>Maximal Velocity LuxI</b></td>
 +
<td>V<sub>max</sub></td>
 +
<td>1.1 mol/min</td>
 +
</tr>
 +
</tr>
 +
<tr class="odd">
 +
<td><b>Km SAM</b></td>
 +
<td>K<sub>m</sub></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>
 +
</table>
 +
 
 +
<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>
 +
 
 +
<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>
 +
 
 +
<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>
 +
<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>
 +
 
 +
 
 +
<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>
 +
 
 +
 
 +
<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>
 +
 
 +
<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>
 +
 
 +
<h3>Results</h3>
 +
 
 +
<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>   
 +
 
 +
<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>        
 +
            </div>
</body>
</body>
</html>
</html>
-
[[Category:BIOTEC_Dresden/Project|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.

Share to Twitter Share to Facebook Share to Orkut Stumble It Email This More...