Team:HIT-Harbin/Application

Overview for our application

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.

1 lactobacillus: from lactose to lactic acid

Reaction 1: endocytosis transportation 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: 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.