Team:Paris Saclay/simulation

From 2013.igem.org

(Difference between revisions)
(simulation.hsim)
(simulation.hsim)
Line 1: Line 1:
{{Team:Paris_Saclay/incl_debut_generique}}
{{Team:Paris_Saclay/incl_debut_generique}}
='''simulation.hsim'''=
='''simulation.hsim'''=
-
<p>title="FNR regulator simulation";
+
title="FNR regulator simulation scenario 2";
-
<br>
+
 
 +
 
 +
 
geometry = 200:50; // size of E.coli in unit of 10nm 200 : 50=2:0.5 µm
geometry = 200:50; // size of E.coli in unit of 10nm 200 : 50=2:0.5 µm
-
<br>
+
 
 +
 
 +
 
//molecules list
//molecules list
-
<br>
 
-
molecule Pfnr;//fnr promotor<br>
 
-
molecule Ise,ClpXP;//Ise: iron-sulfur cluster assembly proteins : enzyme of dimerisation ; ClpXP: ATP dependent protease that induces degradation for all protein<br>
 
-
molecule QFeFNRa,apoFNR,DFeFNRi;//3 states of FNR: activated, without Fe (inactivated), and inactivated<br>
 
-
molecule GPV,PV;//GPV green protein gene, PV green protein<br>
 
-
molecule GPR,PR;//GPR red protein gene, PR red protein<br>
+
molecule mfnr;
 +
 
 +
molecule FNRi;
 +
 
 +
molecule FNRa;
 +
 
 +
molecule o2;
 +
 
 +
molecule ClpXP;
 +
 
 +
 
 +
 
 +
molecule proteinase1,proteinase2,proteinase3;
 +
 
 +
molecule GPV;
 +
 
 +
molecule PV (0,200,0);
 +
 
 +
molecule GPR;
 +
 
 +
molecule PR (200,0,0);  
 +
 
 +
molecule GPR_binded;
 +
 
 +
molecule FNRa_binded_GPR;
 +
 
 +
molecule FNRa_binded_PV;
 +
 
 +
molecule repression_red_protein;
 +
 
 +
molecule FNRa_binded;
 +
 
 +
 
 +
 
 +
 
 +
 
 +
//size in 10nm unit
 +
 
 +
 
 +
 
 +
size(mfnr)=0.1;
 +
 
 +
 
 +
 
 +
size(ClpXP)=0.5;
 +
 
 +
size(FNRa)=0.6;
 +
 
 +
size(FNRi)=0.3;
 +
 
 +
size(o2)=0.01;
 +
 
 +
size(FNRa_binded)=1.2;
 +
 
 +
size(proteinase1)=0.5;
 +
 
 +
size(proteinase2)=0.5;
 +
 
 +
size(proteinase3)=0.5;
 +
 
 +
 
 +
 
 +
size(GPV)=0.01;
 +
 
 +
size(PV)=1;
 +
 
 +
size(GPR)=0.01;
 +
 
 +
size(PR)=1;
 +
 
 +
size(GPR_binded)=0.02;
 +
 
 +
size(FNRa_binded_GPR)=1.3;
 +
 
 +
size(FNRa_binded_PV)=2.3;
 +
 
 +
size(repression_red_protein)=4;
 +
 
 +
 
 +
 
 +
// speed : probability to move to a distance of 10 nm
 +
 
 +
 
-
molecule QFeFNRa_GPR_binded,QFeFNRa_GPV_binded;//molecule FNR binds on mRNA from GPR and GPV<br>
+
speed(mfnr)=0.1;
-
molecule rep;//complexe of 4 FNRs, prevent production of red protein<br>
 
 +
speed(ClpXP)=0.5;
-
metabolite o2;//declare the o2 as a cytosolic molecule leads to be treated as an homogeneous population<br>
+
speed(FNRa)=0.3;
 +
speed(FNRi)=0.3;
-
//size in 10nm unit<br>
 
-
size(Pfnr)=0.01;<br>
+
speed(GPV)=01;
-
<br>
+
-
size(Ise)=0.5;
+
-
size(ClpXP)=0.5;<br>
+
speed(PV)=0.3;
-
size(QFeFNRa)=2;<br>
+
speed(GPR)=0.1;
-
size(apoFNR)=1;<br>
+
speed(PR)=0.3;
-
size(DFeFNRi)=2;<br>
+
speed(o2)=0.8;
-
size(GPV)=0.01;<br>
+
speed(FNRa_binded)=0.1;
-
size(PV)=1;<br>
+
speed(FNRa_binded_GPR)=0.1;
-
size(GPR)=0.01;<br>
 
-
size(PR)=1;<br>
 
-
size(QFeFNRa_GPR_binded)=2.01;<br>
+
speed(repression_red_protein)=0.01;
-
size(QFeFNRa_GPV_binded)=2.01;<br>
+
speed(GPR_binded)=0.1;
-
size(rep)=4.02;<br>
+
speed(FNRa_binded_PV) = 0.1;
-
//speed <br>
+
maxlinks (GPR) = FNRa_binded(1);
-
speed(Pfnr)=0;<br>
+
maxlinks (GPV) = FNRa_binded(1);
-
speed(Ise)=0.5;<br>
 
-
speed(ClpXP)=0.5;<br>
 
-
speed(QFeFNRa)=0.3;<br>
 
-
speed(DFeFNRi)=0.3;<br>
 
-
speed(apoFNR)=0.3;<br>
 
-
speed(GPV)=0;<br>
 
-
speed(PV)=0.3;<br>
 
-
speed(GPR)=0;<br>
 
-
speed(PR)=0.3;<br>
 
-
speed(QFeFNRa_GPR_binded)=0;<br>
 
-
speed(QFeFNRa_GPV_binded)=0;<br>
+
//rules
-
speed(rep)=0;<br>
+
//translation fnr to FNR inactived
 +
mfnr -> mfnr + FNRi [0.001];
-
//reactions<br>
 
 +
//2 FNRi bind
 +
FNRi + FNRi -> FNRa [0.5];
-
//translating Pfnr to FNR<br>
 
-
Pfnr -> Pfnr + apoFNR [0.1];//fnr mRNA generates steagyly apoFNR(inactived)<br>
 
 +
//to many FNRa will reduce FNR production
 +
FNRa + mfnr -> FNRa*mfnr [0.1];
-
//dimerisation <br>
+
FNRa*mfnr -> FNRa + mfnr [0.1];
-
apoFNR + Ise -> QFeFNRa + Ise[0.8];//production of 4Fe-FNR by Ise is steady<br>
 
 +
//deactivation FNR
-
//oxydation<br>
+
FNRa + o2 -> FNRi + FNRi [0.9];
-
QFeFNRa + o2 -> DFeFNRi +o2 [0.8];//aerobically, oxygen inactivates FNR<br>
 
-
DFeFNRi + Ise -> QFeFNRa + Ise [0.5];//but cell continues to reactivate it steadly<br>
 
 +
//degradation FNR
 +
FNRi + ClpXP -> ClpXP [0.05];
-
DFeFNRi + o2 -> apoFNR +o2 [0.6];//return to apoFNR<br>
+
FNRa + proteinase1 -> proteinase1 [0.05];//taux de dissociation est different
-
//degradation FNR<br>
+
//green reporter protein  production
-
apoFNR + ClpXP -> ClpXP [0.05];<br>
+
FNRa + FNRa -> FNRa_binded [0.5];
-
DFeFNRi + ClpXP -> ClpXP [0.05];<br>
+
FNRa_binded -> FNRa + FNRa [0.001];
-
QFeFNRa + ClpXP -> ClpXP [0.05];<br>
 
 +
FNRa_binded + GPV -> FNRa_binded_PV [0.1];
 +
FNRa_binded_PV -> FNRa_binded + GPV[0.001];
 +
FNRa_binded_PV -> PV + FNRa_binded_PV[0.05];
-
//activation/inactivation of green gene<br>
 
-
QFeFNRa + GPV -> QFeFNRa_GPV_binded [0.8];//association<br>
 
-
QFeFNRa_GPV_binded -> QFeFNRa + GPV [0.4];//dissociation<br>
 
-
QFeFNRa_GPV_binded -> PV + QFeFNRa_GPV_binded [0.3];//production of green protein<br>
 
-
//green protein degradation<br>
 
-
PV + ClpXP -> ClpXP [0.05];<br>
+
//red protein repression with 4 FNR activated
-
//repression of red gene<br>
+
FNRa_binded + GPR -> FNRa_binded_GPR [0.8];
-
QFeFNRa + GPR -> QFeFNRa_GPR_binded [0.8];<br>
+
FNRa_binded_GPR + FNRa_binded_GPR -> repression_red_protein [0.5];
-
QFeFNRa_GPR_binded -> QFeFNRa + GPR [0.1];<br>
+
repression_red_protein -> FNRa_binded_GPR + FNRa_binded_GPR [0.00001];
-
QFeFNRa_GPR_binded + QFeFNRa_GPR_binded -> rep [0.5];<br>
+
FNRa_binded_GPR -> FNRa_binded + GPR[0.0001];
-
rep -> QFeFNRa_GPR_binded + QFeFNRa_GPR_binded [0.05];<br>
 
-
//expression the red gene in aerobic condition<br>
 
-
GPR -> GPR + PR [0.3];//production of red protein<br>
+
//production of red protein when aerobic
-
QFeFNRa_GPR_binded -> QFeFNRa_GPR_binded + PR [0.1];<br>
+
GPR + GPR -> GPR_binded[0.05];
 +
GPR_binded -> GPR_binded + PR [0.05];
 +
GPR_binded -> GPR + GPR[0.001];
-
//degradation of red protein<br>
 
-
PR + ClpXP -> ClpXP [0.05];<br>
 
-
//initialtion<br>
+
//degradation of protein
 +
PR + proteinase2 -> proteinase2 [0.1];
 +
PV + proteinase3 -> proteinase3 [0.1];
-
init(100,Pfnr);<br>
 
-
init(100,Ise);<br>
 
-
init(100,ClpXP);<br>
+
//initialtion
-
init(100,QFeFNRa);<br>
 
-
init(100,DFeFNRi);<br>
 
-
init(100,apoFNR);<br>
+
init(10,mfnr);
-
init(100,GPV);<br>
 
-
init(10,PV);<br>
 
-
init(100,GPR);<br>
+
init(100,ClpXP);
-
init(10,PR);<br>
+
init(100,GPV);
-
init(10,o2);<br>
+
init(100,GPR);
-
init(10,QFeFNRa_GPR_binded);<br>
+
init(1000,o2);
-
init(10,QFeFNRa_GPV_binded);<br>
+
init(100,proteinase1);
-
init(10,rep);<br>
+
init(100,proteinase2);
-
</p>
+
 +
init(100,proteinase3);
{{Team:Paris_Saclay/incl_fin}}
{{Team:Paris_Saclay/incl_fin}}

Revision as of 17:18, 3 October 2013

simulation.hsim

title="FNR regulator simulation scenario 2";


geometry = 200:50; // size of E.coli in unit of 10nm 200 : 50=2:0.5 µm


//molecules list


molecule mfnr;

molecule FNRi;

molecule FNRa;

molecule o2;

molecule ClpXP;


molecule proteinase1,proteinase2,proteinase3;

molecule GPV;

molecule PV (0,200,0);

molecule GPR;

molecule PR (200,0,0);

molecule GPR_binded;

molecule FNRa_binded_GPR;

molecule FNRa_binded_PV;

molecule repression_red_protein;

molecule FNRa_binded;



//size in 10nm unit


size(mfnr)=0.1;


size(ClpXP)=0.5;

size(FNRa)=0.6;

size(FNRi)=0.3;

size(o2)=0.01;

size(FNRa_binded)=1.2;

size(proteinase1)=0.5;

size(proteinase2)=0.5;

size(proteinase3)=0.5;


size(GPV)=0.01;

size(PV)=1;

size(GPR)=0.01;

size(PR)=1;

size(GPR_binded)=0.02;

size(FNRa_binded_GPR)=1.3;

size(FNRa_binded_PV)=2.3;

size(repression_red_protein)=4;


// speed : probability to move to a distance of 10 nm


speed(mfnr)=0.1;


speed(ClpXP)=0.5;

speed(FNRa)=0.3;

speed(FNRi)=0.3;


speed(GPV)=01;

speed(PV)=0.3;

speed(GPR)=0.1;

speed(PR)=0.3;

speed(o2)=0.8;

speed(FNRa_binded)=0.1;

speed(FNRa_binded_GPR)=0.1;


speed(repression_red_protein)=0.01;

speed(GPR_binded)=0.1;

speed(FNRa_binded_PV) = 0.1;


maxlinks (GPR) = FNRa_binded(1);

maxlinks (GPV) = FNRa_binded(1);






//rules

//translation fnr to FNR inactived

mfnr -> mfnr + FNRi [0.001];


//2 FNRi bind

FNRi + FNRi -> FNRa [0.5];


//to many FNRa will reduce FNR production

FNRa + mfnr -> FNRa*mfnr [0.1];

FNRa*mfnr -> FNRa + mfnr [0.1];


//deactivation FNR

FNRa + o2 -> FNRi + FNRi [0.9];


//degradation FNR

FNRi + ClpXP -> ClpXP [0.05];

FNRa + proteinase1 -> proteinase1 [0.05];//taux de dissociation est different


//green reporter protein production

FNRa + FNRa -> FNRa_binded [0.5];

FNRa_binded -> FNRa + FNRa [0.001];


FNRa_binded + GPV -> FNRa_binded_PV [0.1];

FNRa_binded_PV -> FNRa_binded + GPV[0.001];

FNRa_binded_PV -> PV + FNRa_binded_PV[0.05];




//red protein repression with 4 FNR activated


FNRa_binded + GPR -> FNRa_binded_GPR [0.8];

FNRa_binded_GPR + FNRa_binded_GPR -> repression_red_protein [0.5];

repression_red_protein -> FNRa_binded_GPR + FNRa_binded_GPR [0.00001];

FNRa_binded_GPR -> FNRa_binded + GPR[0.0001];



//production of red protein when aerobic

GPR + GPR -> GPR_binded[0.05];

GPR_binded -> GPR_binded + PR [0.05];

GPR_binded -> GPR + GPR[0.001];



//degradation of protein

PR + proteinase2 -> proteinase2 [0.1];

PV + proteinase3 -> proteinase3 [0.1];


//initialtion


init(10,mfnr);


init(100,ClpXP);

init(100,GPV);

init(100,GPR);

init(1000,o2);

init(100,proteinase1);

init(100,proteinase2);

init(100,proteinase3);