Team:DTU-Denmark/Modelling

From 2010.igem.org

(Difference between revisions)
Line 1: Line 1:
-
__NOTOC__
 
<html>
<html>
<img src="https://static.igem.org/mediawiki/2010/6/6b/Wiki_banner_967px.png"  align="center"> </img>
<img src="https://static.igem.org/mediawiki/2010/6/6b/Wiki_banner_967px.png"  align="center"> </img>
Line 216: Line 215:
</font>
</font>
</table>
</table>
 +
<h1>Modeling approach</h1>
 +
<p align="justify">The simplified version of our <a href="https://2010.igem.org/Team:DTU-Denmark/Switch#final" target="_blank">Bi[o]stable switch</a> illustrated in Figure 1 was used in the modeling. A system must possess two or more stable states for it to function as a switch. The essential components of the switch are the two repressors and two repressible promoters. Each promoter is inhibited by the repressor expressed from the opposing promoter. The aim of our project is to design a robust bistable switch that exhibits bi-stability over a wide-range of parameters and that the two states are tolerant of the fluctuations inherent in gene expression. This will allow the latest induced state to remain constant over cell-generations. The bi-stability will arise as a result of the mutually inhibitory arrangement. As denoted in Table 1, the absence of inducers will leave the cell in its current state, where either promoter 1 expresses repressor 1 or promoter 2 expresses repressor 2. The expression of the repressor will prevent expression from the opposing promoter, preventing the spontaneous change of state.</p>
 +
<p align="justify">The change of state is accomplished by transiently introducing an inducer (or an anti-repressor protein as in our case) of the currently active repressor. Maximal expression from the opposing promoter and establishment of the alternative stable state will ensue.</p>
 +
 +
<h3>Michaelis-Menten Function</h3>
 +
<p align="justify">In our repressor-repressor switch, the simplest switch is one where each repressor binds to its operator by mass action and the rate of transcription is proportional to the amount of unbound operator. A Michaelis-Menten equation can approximate the fraction of bound operator when the molar quantity of the repressor is greater than the operator. The rate of production of repressor protein 1 is given by:</p>
 +
 +
<p align="justify">By expressing K<sub><i>d</i></sub> in terms of repressor :</p>
 +
 +
<p align="justify">And conversely, the rate of production for repressor protein 2 is given by:</p>
 +
 +
<p align="justify">This results in the following repressor-repressor system:</p>
 +
 +
 +
<p align="justify">The nullclines are defined by:</p>
 +
 +
<p align="justify">At equilibrium:</p>
 +
 +
<p align="justify">This results in a linear fractional transformation with a negative second derivative everywhere and is unable to have multiple intersections with the diagonal. As the multiple intersections of the second derivative with the diagonal is a requirement for a switch, the repression described can therefore not support a switch. However, the inclusion of the depletion term and more than one operator site per cell makes a switch possible. The Hill function that takes cooperativity of binding when describing the level of repression is:</p>
 +
 +
<p align="justify">The behavior and conditions for bistability can be elucidated with the following dimensionless model by replacing the Hill equation into the repressor-repressor system:</p>
 +
 +
<p align="justify">The parameters k<sub>1</sub> and k<sub>2</sub> describes the net effect of RNA polymerase binding, open-complex formation, transcription elongation, transcription termination, repressor binding, ribosome binding and translation. The cooperativity, described by n, can be due to the multimerization of the repressor proteins as well as the cooperative binding of repressor multimers to multiple binding sites within the operator. The nullclines of the system are:</p>
 +
 +
<p align="justify">The nullclines intersect at three points, resulting in two steady states and one unstable state. From this it is clear that three intersection points occur due to the sigmoid shape of the graph, which arises when n > 1 showing that cooperativity of binding is necessary. Another key that can be extracted is the importance of similar rates of synthesis of both repressors. If the rates are unbalanced, the nullclines will only intersect once producing only a single stable state. An increase in the rate of repressor synthesis will result in a more robust switch.</p>
 +
<h3>Robustness of the switch</h3>
 +
<p align="justify"></p>
</td>
</td>
<td width="163px" height="100%" valign="top">
<td width="163px" height="100%" valign="top">

Revision as of 19:57, 27 October 2010

Welcome to the DTU iGEM wiki!


UNDER CONSTRUCTION

Introduction

One of the ideas behind synthetic biology is to engineer artificial systems in biological cells. This can be achieved by using well-characterized and simplified parts that can be found in nature. We will be using parts found in different systems in nature to design a bistable switch, an idea stolen from the world of electronics. A simplified version of our bistable switch is illustrated in Figure 1.

Figure 1: Simplified version of our bistable swicth


This simplified version of our switch, the repressor-repressor switch is analogous SR flip-flop of digital electronics. The SR (set-reset) flip-flop is constructed from a pair of cross-linked NAND or NOR logic gates. An SR flip-flop circuit built with NOR logic gates is illustrated in Figure 2.

Figure 2: SR flip-flop


The circuit, built with the pair of cross-linked NOR logic gates, functions such that output from a gate will only be true if no input is given as demonstrated by the truth table in Table 1. The biological interpretation of this logical operation shows is that the repression by a protein can be prevented either by an inducer (or anti-repressor protein as in our case) (Input 1/2) or by repression of its synthesis by another repressor-protein (Output 2/1). The elements in the circuit do not possess memory in them, but by cross-linking the elements, the circuit is capable of “remembering” and holding its current state. This means that when neither Input 1 nor Input 2 (no inducers) is present, the circuit is bistable with either Output 2 or Output 1 being active (only one of the repressor-proteins are expressed), respectively.

Table 1

Input 1
(Anti-Repressor 1)
Input 2
(Anti-Repressor 2)
Output 1
(Reporter 1)
Output 2
(Reporter 2)
00No changeNo change
0110
1001
11InvalidInvalid

Modeling approach

The simplified version of our Bi[o]stable switch illustrated in Figure 1 was used in the modeling. A system must possess two or more stable states for it to function as a switch. The essential components of the switch are the two repressors and two repressible promoters. Each promoter is inhibited by the repressor expressed from the opposing promoter. The aim of our project is to design a robust bistable switch that exhibits bi-stability over a wide-range of parameters and that the two states are tolerant of the fluctuations inherent in gene expression. This will allow the latest induced state to remain constant over cell-generations. The bi-stability will arise as a result of the mutually inhibitory arrangement. As denoted in Table 1, the absence of inducers will leave the cell in its current state, where either promoter 1 expresses repressor 1 or promoter 2 expresses repressor 2. The expression of the repressor will prevent expression from the opposing promoter, preventing the spontaneous change of state.

The change of state is accomplished by transiently introducing an inducer (or an anti-repressor protein as in our case) of the currently active repressor. Maximal expression from the opposing promoter and establishment of the alternative stable state will ensue.

Michaelis-Menten Function

In our repressor-repressor switch, the simplest switch is one where each repressor binds to its operator by mass action and the rate of transcription is proportional to the amount of unbound operator. A Michaelis-Menten equation can approximate the fraction of bound operator when the molar quantity of the repressor is greater than the operator. The rate of production of repressor protein 1 is given by:

By expressing Kd in terms of repressor :

And conversely, the rate of production for repressor protein 2 is given by:

This results in the following repressor-repressor system:

The nullclines are defined by:

At equilibrium:

This results in a linear fractional transformation with a negative second derivative everywhere and is unable to have multiple intersections with the diagonal. As the multiple intersections of the second derivative with the diagonal is a requirement for a switch, the repression described can therefore not support a switch. However, the inclusion of the depletion term and more than one operator site per cell makes a switch possible. The Hill function that takes cooperativity of binding when describing the level of repression is:

The behavior and conditions for bistability can be elucidated with the following dimensionless model by replacing the Hill equation into the repressor-repressor system:

The parameters k1 and k2 describes the net effect of RNA polymerase binding, open-complex formation, transcription elongation, transcription termination, repressor binding, ribosome binding and translation. The cooperativity, described by n, can be due to the multimerization of the repressor proteins as well as the cooperative binding of repressor multimers to multiple binding sites within the operator. The nullclines of the system are:

The nullclines intersect at three points, resulting in two steady states and one unstable state. From this it is clear that three intersection points occur due to the sigmoid shape of the graph, which arises when n > 1 showing that cooperativity of binding is necessary. Another key that can be extracted is the importance of similar rates of synthesis of both repressors. If the rates are unbalanced, the nullclines will only intersect once producing only a single stable state. An increase in the rate of repressor synthesis will result in a more robust switch.

Robustness of the switch