"
Team:ETH Zurich/modeling/int
From 2014.igem.org
(→Reactions) |
(→Model) |
||
| Line 40: | Line 40: | ||
: Bxb1 + Bxb1 ⇔ DBxb1 | : Bxb1 + Bxb1 ⇔ DBxb1 | ||
: DBxb1 + SI<sub>Bxb1</sub> ⇔ SA<sub>Bxb1</sub> | : DBxb1 + SI<sub>Bxb1</sub> ⇔ SA<sub>Bxb1</sub> | ||
| + | : Bxb1 ⇒ | ||
| + | : DBxb1 ⇒ | ||
<html></b></html></center> | <html></b></html></center> | ||
| Line 46: | Line 48: | ||
: ΦC31 + ΦC31 ⇔ ΦC31 | : ΦC31 + ΦC31 ⇔ ΦC31 | ||
: DΦC31 + SI<sub>ΦC31</sub> ⇔ SA<sub>ΦC31</sub> | : DΦC31 + SI<sub>ΦC31</sub> ⇔ SA<sub>ΦC31</sub> | ||
| + | : ΦC31 ⇒ | ||
| + | : DΦC31 ⇒ | ||
<html></b></html></center> | <html></b></html></center> | ||
| + | |||
| + | === Differential Equations === | ||
| + | |||
| + | Applying mass action kinetic laws, we obtain the following set of differential equations. | ||
| + | |||
| + | $$\frac{d[Bxb1]}{dt}=-2 k_{dimBxb1}[Bxb1]^2+ 2 k_{-dimBxb1}[DBxb1]-d_{Bxb1}[Bxb1]$$ | ||
| + | |||
| + | $$\frac{d[DBxb1]}{dt}=-k_{DNABxb1}[DBxb1][SI_{Bxb1}]+k_{-DNABxb1}[SA_{Bxb1}]+k_{dimBxb1}[Bxb1]^2-k_{-dimBxb1}[DBxb1a]-d_{DBxb1}[DBxb1]$$ | ||
| + | |||
| + | $$\frac{d[SA_{Bxb1}]}{dt}=k_{DNABxb1}[DBxb1][SI_{Bxb1}]-k_{-DNABxb1}[SA_{Bxb1}]$$ | ||
<html></article></html> | <html></article></html> | ||
Revision as of 13:33, 11 October 2014
Integrases
Model
In our design, integrases compute the output of the logic gates. Integrases allow to flip one fragment of DNA. The model we developped is described here.
Chemical Species
- Bxb1: integrase that can fold into two conformations called Bxb1a and Bxb1b. Even if those conformations have the same sequence of amino acids, their tertiary structure is different. As it is experimentally difficult to differentiate these two conformations, we chose to regroup them under the specie Bxb1.
- ΦC31: integrase that can fold into two conformations called ΦC31a and ΦC31b. Similarly to Bxb1, we chose to factorize these two conformations under the specie ΦC31.
- DBxb1: dimer of Bxb1. Bxb1a (respectively Bxb1b) forms a dimer DBxb1a (respectively DBxb1b). To be coherent, dimers re not differentiated depending on their spatial configuration.
- DΦC31: dimer of ΦC31. ΦC31a (respectively ΦC31b) forms a dimer DΦC31a (respectively DΦC31b). The different spatial configuration are not taken into account.
Modelization of DNA-binding sites
Each dimer of integrases can specifically bind to a DNA binding site. As the flipping is irreversible, these DNA binding sites can be three possible states:
- SIIntegraseName: inactive DNA binding site. No dimer is bound to this site, which has never been flipped.
- SAIntegraseName: active DNA binding site. A dimer is to this site.
- SFIntegraseName: flipped DNA binding site. This DNA binding site has been used by a flipping.
Reactions
- For Bxb1
- Bxb1 + Bxb1 ⇔ DBxb1
- DBxb1 + SIBxb1 ⇔ SABxb1
- Bxb1 ⇒
- DBxb1 ⇒
- For ΦC31
- ΦC31 + ΦC31 ⇔ ΦC31
- DΦC31 + SIΦC31 ⇔ SAΦC31
- ΦC31 ⇒
- DΦC31 ⇒
Differential Equations
Applying mass action kinetic laws, we obtain the following set of differential equations.
$$\frac{d[Bxb1]}{dt}=-2 k_{dimBxb1}[Bxb1]^2+ 2 k_{-dimBxb1}[DBxb1]-d_{Bxb1}[Bxb1]$$
$$\frac{d[DBxb1]}{dt}=-k_{DNABxb1}[DBxb1][SI_{Bxb1}]+k_{-DNABxb1}[SA_{Bxb1}]+k_{dimBxb1}[Bxb1]^2-k_{-dimBxb1}[DBxb1a]-d_{DBxb1}[DBxb1]$$
$$\frac{d[SA_{Bxb1}]}{dt}=k_{DNABxb1}[DBxb1][SI_{Bxb1}]-k_{-DNABxb1}[SA_{Bxb1}]$$
Characterization
Data
Assumptions
Parameter fitting
Range of validity of the assumptions
Characterization
Data
Assumptions
Parameter fitting
Range of validity of the assumptions
