Team:HIT-Harbin/Application

From 2013.igem.org

(Difference between revisions)
 
(3 intermediate revisions not shown)
Line 11: Line 11:
                 <div class="mainarea">
                 <div class="mainarea">
                 <div class="multi-columns">
                 <div class="multi-columns">
-
<h5 class="hashed"><span>Overview for our application</span></h5>
+
<h5 class="hashed"><span>1 lactobacillus: from lactose to lactic acid</span></h5>
-
<p>In applications, we plan to enhance Bio-electric Interface, the device designed by 2012 Edinburg iGEM team. If the input is biochemical signal molecules and the output become electrons, B-POM can be coupled with Bio-electric Interface, that ephemeral processes in cells will be precisely measurable by simple electronic methods and analyzable by computer. B-POM is also helpful in yoghurt producing, where the control of pH is inaccurate. When B-POM is transplanted into yoghurt-producing bacteria, let hydrogen ions be the input and lacR be the output, and select the proper B-POM parameters, then the pH will be steady around 5.5. In short, our inspired and reliable device is promising in various fields.</p>
+
-
  1 lactobacillus: from lactose to lactic acid
+
<p>Reaction 1: endocytosis transportation of lactose
-
Reaction 1: endocytosis transportation of lactose
+
<p>Table 1 the definition variables
-
Table 1 the definition variables
+
<p>The unit of concentration is mmol/L.
-
Name Definition
+
<p>LS+Lout <=> temp1 : R1(lactose combines with transporting enzyme)
-
LS Concentration of transporting enzyme
+
<p>vf = k1*Lout*LS(the reaction rate is proportional to reactants’ concentrations)
-
Lout Concentration of lactose outside cells
+
<p>vr =k_1*temp1(the reverse reaction rate is proportional to product’s concentration)
-
L Concentration of lactose in cells
+
<p>temp1 => L+LS : R2(lactose is transported into the cell)
-
temp Concentration of lactose combined with transporting enzyme
+
<p>vf = k2*temp1(the reaction rate is proportional to reactants’ concentrations)
-
vf Reaction rate
+
-
vr Reverse reaction rate
+
-
k1 Rate constant of combination of lactose and transporting enzyme
+
-
k_1 Decomposition rate constant
+
-
k2 Transport rate constant of lactose
+
-
The unit of concentration is mmol/L.
+
-
LS+Lout <=> temp1 : R1(lactose combines with transporting enzyme)
+
-
vf = k1*Lout*LS(the reaction rate is proportional to reactants’ concentrations)
+
-
vr =k_1*temp1(the reverse reaction rate is proportional to product’s concentration)
+
-
temp1 => L+LS : R2(lactose is transported into the cell)
+
-
vf = k2*temp1(the reaction rate is proportional to reactants’ concentrations)
+
-
Lout is known as 136mmol/L. We estimate that the time for the overall transportation of lactose into cells is 400s and assume that k1=0.1,k_1=0.05,k2=0.1,LS=5,temp1=0,L=0。Then the simulation result is as below:
+
<p>Lout is known as 136mmol/L. We estimate that the time for the overall transportation of lactose into cells is 400s and assume that k1=0.1,k_1=0.05,k2=0.1,LS=5,temp1=0,L=0。Then the simulation result is as below:
-
+
<img src="https://static.igem.org/mediawiki/2013/1/1d/0001.png" alt=""/>
-
Fig 1. Simulation of lactose transportation
+
<p>Fig 1. Simulation of lactose transportation
-
We can get from the simulation that the transportation is completed when time is 300s.
+
<p>We can get from the simulation that the transportation is completed when time is 300s.
-
Reaction 2: from lactose to lactic acid
+
<p>Reaction 2: from lactose to lactic acid
-
Table 2 the definition variables
 
-
Name Definition Value
 
-
LZ Concentration of galactosidase 5
 
-
L Concentration of lactose in cells 136
 
-
Glu Concentration of glucose 0
 
-
gala Concentration of galactose 0
 
-
PA Concentration of pyruvic acid 0
 
-
LA Concentration of lactic acid 0
 
-
k3 Rate constant of combination of lactose and galactosidase 0.1
 
-
k_3 Decomposition rate constant of lactose-galactosidase complex 0.05
 
-
k4 Decomposition rate constant of lactose 0.1
 
-
k5 Rate constant from glucose to pyruvic acid 0.1
 
-
rm1 Maximum rate of pyruvic acid producing 35
 
-
Km1 Michaelis constant 100
 
-
L+LZ <=> temp2 : R3(lactose combines with galactosidase)
+
<p>L+LZ <=> temp2 : R3(lactose combines with galactosidase)
-
vf = k3*L*LZ(the reaction rate is proportional to reactants’ concentrations)
+
<p>vf = k3*L*LZ(the reaction rate is proportional to reactants’ concentrations)
-
vr = k_3*temp2(the reverse reaction rate is proportional to product’s concentration)
+
<p>vr = k_3*temp2(the reverse reaction rate is proportional to product’s concentration)
-
temp2 => Glu+LZ+gala : R4(lactose decomposes into glucose and galactose)
+
<p>temp2 => Glu+LZ+gala : R4(lactose decomposes into glucose and galactose)
-
vf = k4*temp2(the reaction rate is proportional to reactants’ concentrations)
+
<p>vf = k4*temp2(the reaction rate is proportional to reactants’ concentrations)
-
Glu => 2*PA : R5(glucose decomposes into pyruvic acid, multi-step reaction)
+
<p>Glu => 2*PA : R5(glucose decomposes into pyruvic acid, multi-step reaction)
-
vf = k5*Glu(decomposition rate is proportional to glucose concentration)
+
<p>vf = k5*Glu(decomposition rate is proportional to glucose concentration)
-
PA => LA : R6(from pyruvic acid to lactic acid)
+
<p>PA => LA : R6(from pyruvic acid to lactic acid)
-
vf = rm1*PA/(Km1+PA)(Michealis equation, the reactant is adequate)
+
<p>vf = rm1*PA/(Km1+PA)(Michealis equation, the reactant is adequate)
-
the simulation result is as below:
+
<p>the simulation result is as below:
 
 
   
   
-
Fig 2. Simulation of lactic acid produciton
+
<p>Fig 2. Simulation of lactic acid produciton
-
We can get from the simulation that all lactose becomes lactic acid. The final concentration of lactic acid is 272 mmol/L.  
+
<p>We can get from the simulation that all lactose becomes lactic acid. The final concentration of lactic acid is 272 mmol/L.  
-
Then we consider reaction1 and reaction2 together, the simulation result is as below:
+
<p>Then we consider reaction1 and reaction2 together, the simulation result is as below:
   
   
-
Fig 3. Considering reacton1 and 2 together
+
<p>Fig 3. Considering reacton1 and 2 together
-
We can see that once lactose enters cells, it is consumed, and the transportation of lactose is faster than the production of lactic acid, which ensure that the production of lactic acid is continuing.
+
<p>We can see that once lactose enters cells, it is consumed, and the transportation of lactose is faster than the production of lactic acid, which ensure that the production of lactic acid is continuing.
-
The cells are assumed in normal condition.
+
<p>The cells are assumed in normal condition.
-
Now we are going to analyze the pH switching function f(p)=1-tanh(n*(p-p*)). The curve of f(x)=1-tanh(x) is as below:
+
<p>Now we are going to analyze the pH switching function f(p)=1-tanh(n*(p-p*)). The curve of f(x)=1-tanh(x) is as below:
   
   
-
Fig 4. Curve of f(x)=1-tanh(x)
+
<p>Fig 4. Curve of f(x)=1-tanh(x)
-
When x=1, f(x) is nearly 0, the switch is off. That is, when n=1, the difference between pH in cell and pH threshold should be smaller than 1 that the switch can be on.  
+
<p>When x=1, f(x) is nearly 0, the switch is off. That is, when n=1, the difference between pH in cell and pH threshold should be smaller than 1 that the switch can be on.  
-
The curve of f(x)=1-tanh(30*x) is as below:
+
<p>The curve of f(x)=1-tanh(30*x) is as below:
   
   
-
Fig 5. Curve of f(x)=1-tanh(30*x)
+
<p>Fig 5. Curve of f(x)=1-tanh(30*x)
-
When x=0.1, f(x) is nearly 0, the switch is off. That is, when n=30, the difference between pH in cell and pH threshold should be smaller than 0.1 that the switch can be on.  
+
<p>When x=0.1, f(x) is nearly 0, the switch is off. That is, when n=30, the difference between pH in cell and pH threshold should be smaller than 0.1 that the switch can be on.  
-
Besides, we get the relationship between pH and lactic acid concentration from experiments:  
+
<p>Besides, we get the relationship between pH and lactic acid concentration from experiments:  
-
pH=6.32311-0.02168*[LA]+0.0000426042*[LA]^2(0<[LA]<250)
+
<p>pH=6.32311-0.02168*[LA]+0.0000426042*[LA]^2(0<[LA]<250)
-
its curve is as below:
+
<p>its curve is as below:
   
   
-
Fig 5. Relationship between pH and lactic acid concentration
+
<p>Fig 5. Relationship between pH and lactic acid concentration
 +
 
 +
<p>In experiments, the maximum of pH is 6.3, the minimum is 3.6. We assume that the threshold of pH is 5.5.
 +
 
 +
<h5 class="hashed"><span>2 the expression of LacR through pH-induced promoter</span></h5>
-
In experiments, the maximum of pH is 6.3, the minimum is 3.6. We assume that the threshold of pH is 5.5.
+
<p>=> LA
 +
<p>vf = 0
 +
<p>=> mRNA1
 +
<p>vf = v0+v1*f(LA)-dm1*mRNA1
-
2 the expression of LacR through pH-induced promoter
+
<p>=> LR
-
=> LA
+
<p>vf = km*mRNA1-dn1*LR
-
vf = 0
+
-
=> mRNA1
+
-
vf = v0+v1*f(LA)-dm1*mRNA1
+
-
  => LR
+
  <p>=> mRNA2
-
vf = km*mRNA1-dn1*LR
+
<p>vf = vm-dm2*mRNA2-c*a*LR/(1+a*LR)
-
  => mRNA2
+
  <p>=> LS
-
vf = vm-dm2*mRNA2-c*a*LR/(1+a*LR)
+
<p>vf = ks*mRNA2-ds*LS
-
  => LS
+
  <p>=> LZ
-
vf = ks*mRNA2-ds*LS
+
<p>vf = kz*mRNA2-dz*LZ
-
=> LZ
+
<p>Table 3 the value of variables
-
vf = kz*mRNA2-dz*LZ
+
-
Table 3 the value of variables
 
-
NAME VALUE
+
<p>the simulation result is as below:
-
LA(0) 272
+
-
LR(0) 10
+
-
LS(0) 5
+
-
LZ(0) 5
+
-
mRNA1(0) 1
+
-
mRNA2(0) 1
+
-
v0 1
+
-
v1 30
+
-
dm1 0.5
+
-
km 1
+
-
dn1 0.5
+
-
vm 1
+
-
c 1.2
+
-
a 1
+
-
dm2 0.2
+
-
ks 1
+
-
ds 0.5
+
-
kz 1
+
-
dz 0.5
+
-
the simulation result is as below:
+
   
   
-
Fig 6. Simulation of pH induction
+
<p>Fig 6. Simulation of pH induction
-
Form the figure we can see that in 20s (far from 400s) the concentration of transporting enzyme and galactosidase is very small, which means the inhibition resulted from pH induction is transitory for production of lactic acid. In other word, the inhibition can respond to the fast change of pH.
+
<p>Form the figure we can see that in 20s (far from 400s) the concentration of transporting enzyme and galactosidase is very small, which means the inhibition resulted from pH induction is transitory for production of lactic acid. In other word, the inhibition can respond to the fast change of pH.
-
Finally, we put pH induction and lactic acid production together.
+
<p>Finally, we put pH induction and lactic acid production together.
-
When n=1, the simulation is as below:
+
<p>When n=1, the simulation is as below:
   
   
-
Fig 7. Taking pH induction and lactic acid production together (n=1)
+
<p>Fig 7. Taking pH induction and lactic acid production together (n=1)
-
We can see that the concentration of lactic acid doesn’t decline but only slow down (the inhibition). Time to achieve platform phase increases to 800s.
+
<p>We can see that the concentration of lactic acid doesn’t decline but only slow down (the inhibition). Time to achieve platform phase increases to 800s.
-
When n=30, the simulation is as below:
+
<p>When n=30, the simulation is as below:
   
   
-
Fig 8. Taking pH induction and lactic acid production together (n=30)
+
<p>Fig 8. Taking pH induction and lactic acid production together (n=30)
-
The time to achieve platform phase decreases to 100s.
+
<p>The time to achieve platform phase decreases to 100s.

Latest revision as of 04:14, 28 September 2013

HIT-Harbin

Application

1 lactobacillus: from lactose to lactic acid

Reaction 1: endocytosis transportation of lactose

Table 1 the definition variables

The unit of concentration is mmol/L.

LS+Lout <=> temp1 : R1(lactose combines with transporting enzyme)

vf = k1*Lout*LS(the reaction rate is proportional to reactants’ concentrations)

vr =k_1*temp1(the reverse reaction rate is proportional to product’s concentration)

temp1 => L+LS : R2(lactose is transported into the cell)

vf = k2*temp1(the reaction rate is proportional to reactants’ concentrations)

Lout is known as 136mmol/L. We estimate that the time for the overall transportation of lactose into cells is 400s and assume that k1=0.1,k_1=0.05,k2=0.1,LS=5,temp1=0,L=0。Then the simulation result is as below:

Fig 1. Simulation of lactose transportation

We can get from the simulation that the transportation is completed when time is 300s.

Reaction 2: from lactose to lactic acid

L+LZ <=> temp2 : R3(lactose combines with galactosidase)

vf = k3*L*LZ(the reaction rate is proportional to reactants’ concentrations)

vr = k_3*temp2(the reverse reaction rate is proportional to product’s concentration)

temp2 => Glu+LZ+gala : R4(lactose decomposes into glucose and galactose)

vf = k4*temp2(the reaction rate is proportional to reactants’ concentrations)

Glu => 2*PA : R5(glucose decomposes into pyruvic acid, multi-step reaction)

vf = k5*Glu(decomposition rate is proportional to glucose concentration)

PA => LA : R6(from pyruvic acid to lactic acid)

vf = rm1*PA/(Km1+PA)(Michealis equation, the reactant is adequate)

the simulation result is as below:  

Fig 2. Simulation of lactic acid produciton

We can get from the simulation that all lactose becomes lactic acid. The final concentration of lactic acid is 272 mmol/L.

Then we consider reaction1 and reaction2 together, the simulation result is as below:

Fig 3. Considering reacton1 and 2 together

We can see that once lactose enters cells, it is consumed, and the transportation of lactose is faster than the production of lactic acid, which ensure that the production of lactic acid is continuing.

The cells are assumed in normal condition.

Now we are going to analyze the pH switching function f(p)=1-tanh(n*(p-p*)). The curve of f(x)=1-tanh(x) is as below:

Fig 4. Curve of f(x)=1-tanh(x)

When x=1, f(x) is nearly 0, the switch is off. That is, when n=1, the difference between pH in cell and pH threshold should be smaller than 1 that the switch can be on.

The curve of f(x)=1-tanh(30*x) is as below:

Fig 5. Curve of f(x)=1-tanh(30*x)

When x=0.1, f(x) is nearly 0, the switch is off. That is, when n=30, the difference between pH in cell and pH threshold should be smaller than 0.1 that the switch can be on.

Besides, we get the relationship between pH and lactic acid concentration from experiments:

pH=6.32311-0.02168*[LA]+0.0000426042*[LA]^2(0<[LA]<250)

its curve is as below:

Fig 5. Relationship between pH and lactic acid concentration

In experiments, the maximum of pH is 6.3, the minimum is 3.6. We assume that the threshold of pH is 5.5.

2 the expression of LacR through pH-induced promoter

=> LA

vf = 0

=> mRNA1

vf = v0+v1*f(LA)-dm1*mRNA1

=> LR

vf = km*mRNA1-dn1*LR

=> mRNA2

vf = vm-dm2*mRNA2-c*a*LR/(1+a*LR)

=> LS

vf = ks*mRNA2-ds*LS

=> LZ

vf = kz*mRNA2-dz*LZ

Table 3 the value of variables

the simulation result is as below:

Fig 6. Simulation of pH induction

Form the figure we can see that in 20s (far from 400s) the concentration of transporting enzyme and galactosidase is very small, which means the inhibition resulted from pH induction is transitory for production of lactic acid. In other word, the inhibition can respond to the fast change of pH.

Finally, we put pH induction and lactic acid production together.

When n=1, the simulation is as below:

Fig 7. Taking pH induction and lactic acid production together (n=1)

We can see that the concentration of lactic acid doesn’t decline but only slow down (the inhibition). Time to achieve platform phase increases to 800s.

When n=30, the simulation is as below:

Fig 8. Taking pH induction and lactic acid production together (n=30)

The time to achieve platform phase decreases to 100s.