Team:Tsinghua-A/Modeling/GILLESPIE/part2

From 2012.igem.org

(Difference between revisions)
Line 22: Line 22:
<div id="topTitle">
<div id="topTitle">
-
<h2 id="titleText" style="color:rgb(137,202,154);">Tsinghua-A::Modeling::<span style="color: rgb(12,117,34);">Part2</span></h2>
+
<h2 id="titleText" style="color:rgb(137,202,154);">Tsinghua-A::Modeling::GILLESPIE::<span style="color: rgb(12,117,34);">Part2</span></h2>
</div>
</div>
<div id="topWrapper">
<div id="topWrapper">

Revision as of 11:39, 26 September 2012

Tsinghua-A::Modeling::GILLESPIE::Part2

Our model

In our model , we use Gillespie algorithm to simulate the process of inversion. Here we use Markov chain to describe the process of inversion.
We assume that the state ‘A’ , ‘B’ , ‘C’ , ‘D’ , is the original state , the state after inversion , the Intermediate state when A forms holiday junction, the Intermediate state when B forms holiday junction respectively.

(VAN DUYNE G D, 2001)

(VAN DUYNE G D, 2001)

The following are the reactions which are related to the process of inversion.

Then we get the following Markov chain.

Here,are the probabilities of the reaction in an Infinitesimal time. According to the property of the Markov chain , we can get:

In the model,is proportional to the concentration of Cre-Loxp andare much larger than Here we see 1000 genes which are in state ‘A’ as a whole and we get the following results.
The concentration of Cre-Loxp :

The number of the genes which are in the original state and in the state after inversion:

Return