Team:TJU/Modeling


1 Gene Expression Model for Lactate production



Simplification and Assumption

① Simplifications and assumptions for metabolic pathway and gene insertion:supposing that we have adequate glucose and constant pyruvate, the latter used to product lactate can be replenished by glucose in E.coli. Naturally, E.coli can itself produce lactate dehydrogenase (Ldh) through LdhA gene expression. Later, we attempt to introduce plasmids with Ldh gene, which can produce the extra Ldh as an approach to generate even more lactate. However, the excessive Ldh, to a certain degree, will also inhibit Ldh gene expression.

② Simplifications and assumptions for Michaelis-Menten equation: we want to figure out the relations between the speed of lactate production and the amount of Ldh expression. But, there is a precondition in Michaelis-Menten equation, namely, the total amount of enzyme in solution keeps unchanged.

So, we need to modify the existing Michaelis-Menten equation:


   

The above equation reflects a general enzyme catalytic reaction process. With the help of steady state approximation, the equation can transform from first-order reaction to zero-order reaction. Because of the constant substrate concentration as we assumed, the variable in the equation is the variation of ldh amount, which means that the independent variable in the Michaelis-Menten equation turns into the enzyme amount. (variable change from [S] into [E])


   

in which :

(1)


③ Simplification of transmembrane transfer: the original pH in the medium is slightly greater than 7 and the pH can keep constant with the buffer function of KHPO4 together with KH2PO4. Suppose that Shewanella can absorb lactate brought by E.coli timely, thus, pH can remain stable, which is around 7.2 based on our hypothesis. Further, we assume that transmembrane transfer belongs to codiffusion where the concentration difference of H+ and acid groups can drive the whole process. Therefore, the transmembrane transfer is only related to the intra membrane concentration of H+ and acid groups.


Equation Establishment

Firstly, we establish the equation of ldh amount changing with time. The following equations represent the ldh production rate in E.coli without the plasmid insertion and with the plasmid introduction, respectively.


 

(2)


 

(3)


Next, we construct the lactate production rate equation:


  (4)

Finally, the lactate transmembrane transfer equations has been figured out. It shows that lactate will be transferred to the outside at the time it produced. We represent the intra membrane with c and extra membrane with e.


   

The transmembrane transfer rate:


 

(5)


When lactate production is combined with its transfer, we can get the following equation:


 

(6)




Simulation Results



Figure 1. The relations between amount of ldh and lactate with time passed, respectively.


From the simulation results, we can clearly see that the E. coli with Ldh have a significantly increased in number and lactic acid production of lactate dehydrogenase, which is consist with our experimental results.



2 Models for monocultures of Shewanella in MFC



Shewanella Growth Model Based on Carbon Balance

At the very beginning, we construct a fermentation kinetic model based on the Monod Equation to draw a picture of carbon balance between biomass increase and base consumption during S.oneidensis MR-1 growth. Before generation our modeling, we simplify the center process of carbon metabolism in Shewanella via making assumptions as follows:

① Lactate is the preferable carbon source for Shewanella. Lactate can transform into pyruvate by dehydrogenation, which then converts mostly into acetyl coenzyme. Subsequently, acetyl coenzyme can either be oxidized into CO2 and release energy by TCA Circle or transform into acetate.

② When there are too much lactate, the producing rate of acetyl coenzyme will be greater than the consuming rate of that in TCA circle. The excessive acetyl coenzyme will convert into acetate and release to the outside. Inversely, the shortage of lactate leads the utilization of acetate which then maintains the metabolism by acetyl coenzyme formation.

③ Although both acetate and lactate can be used as carbon sources for Shewanella, acetate can restrain its growth. On the one hand, acetate will drop the acidity inside the cell, which, a step forward, can affect the activity of enzyme in vivo resulted in uncompetitive inhibition to cell growth. On the other hand, acetate can itself limit the growth by uncompetitive inhibition. In particular, the process of acetate metabolism and lactate metabolism are competing for acetyl coenzyme, which reduces the lactate utilization and influences cell growth. (Tang & Meadows, 2006)

④ The redox of mediator in the outer membrane is the rate-limiting step of electron transfer chain in Shewanella. TCA circle, an important pathway for the synthesis of energy and biomass , largely depends on the concentration of electron receptor at the end of electron transfer chain. When there is a deficiency of oxidized mediator, TCA circle will be inhibited, thereby, affect the whole process of metabolism as well as biomass increase.


According to the four assumptions above along with acetate uncompetitive inhibition and competitive inhibition, we decide to use multi-substrates Monod equation(Blanch & Clark, 1997) with two kinds of inhibition to demonstrate the relations between restrictive substrates and biomass increase rate. Considering that the accumulation of NADH will remarkably affect metabolic process under anaerobic conditions, we add a product factor in the Monod equation reflecting metabolites inhibition function. μL and μA refers to the lactate and acetate biomass specific rate, respectively. We quantify the relations between specific rate of lactate/acetate and restrictive base concentration as following:


 

(7)


  (8)

Biomass, lactate and acetate concentration can vary with time and its construction depends on kinetic model which is based on the general kinetics of the fermentation process model. Consumption of lactate and acetate can be both used for biomass increase and endogenous metabolism maintenance. However, considering the rarely energy losses caused by metabolism in vivo, we ignore the carbon consumption brought by metabolism. In this way, lactate can only be used for biomass increase and acetate production while acetate is merely used for biomass increase.


 

(9)


 

(10)


 

(11)



Extracellular Electron Transfer Model

Extracellular Electron Transfer (EET) process can deliver the electrons between electrode and cell membrane with mediating functions of redox active molecules including FMN and riboflavin. It can be subdivided into the oxidation of reduced mediator and reduction of oxidized mediator in electrode.


We establish a model kinetics of electrode process on the basis of Bulter-Volmer equation to illustrate the oxidation of reduced mediator in electrode. Firstly, we make several assumptions just like the following tips:

① Although both of FMN and riboflavin can serve as mediator in electrode redox, we only take riboflavin into consideration for simplification in the model.

② There are two kinds of EET pathway model for riboflavin mediating. One of these called single-electron transfer in which riboflavin acts as cofactor of MtrC protein and can be reduced into semiquinone (Sq) with the standard electrode potential of -260v. Another one is two-electron transfer where the dissociative riboflavin binds to outer membrane c-Cyts as a cofactor in the hydroquinone (Hq) form with shifted potential of -145v. Under conditions of the same electrode and concentration, the exchange current density of single-electron transfer is much high than that of two-electron transfer since enzymatic reaction can greatly decrease its activation energy. Therefore, we assumed that current generation reaction in anode mainly belongs to the single-electron transfer.

③ Oxidation of Sq in electrode is the reversible first order reaction. The net rate for this kind of reaction is made up by the difference between oxidized rate of Sq and reduced rate of riboflavin in electrode.


Accordingly, current density in anode can be calculated by Bulter-Vlomer equation:


 

(12)


Reduction and Oxidation reaction rate constant can be calculated by Arrhenius equation:


 

(13)


 

(14)


As for the reduction of oxidized mediator in outer membrane for Shewanella, we establish a kinetics of enzyme catalyzed reaction on the basis of Michaelis-Menten equation, accordingly.

① Similarly, we also consider the single-electron transfer as the only pathway in vivo.

② In such single-electron transfer process, riboflavin reduction kinetics characteristics follows the general rules of enzymatic reaction.

③ The reduction of riboflavin is related to the reducing state of OM c-Cyts heme center which is also the terminal of EET pathway. Due to the electrons brought by NADH or NADPH from respiration, we argue that the process of riboflavin reduction into Sq is coupled with NADH from respiration process.


Based on the assumptions above, we believe the respiration of OM riboflavin has something to do with its concentration and NADH production rate. According to the kinetics of enzyme catalyzed reaction, we use multi-substrates Michaelis-Menten equation to determine the whole process:


 

(15)


 

(16)


 

(17)



Simulation Results



Figure 2. The growth curve of S.oneidensis in MFC.



Figure 3. The conc. of lactate and acetate varies with time.



Figure 4. The conc. of Sq and Hq varies with time.



Figure 5. The Voltage varies with time.



As can be seen from the simulation results, when the initial inoculation amount Shewanella is 0.1 g/L, it's growth time in anaerobic MFC is approximately 10 hours. The output voltage of the battery is gradually stabilizing after 20 hours, and the voltage drops after 70 hours due to the mass death of Shewanella.Although specific values are different, the simulation results are consistent with the experimental results in the trends of voltage.



Parameters and variables


Parameter Discribtion Number [Unit]
v volume of reactor 140E-3 [L]
μmax,L the maximal specific growth rate taking lactate as substrate 0.52 [h-1]
μmax,A the maximal specific growth rate taking acetate as substrate 0.28 [h-1]
ke endogenous metabolism rate constant 0.08 [h-1]
kc acetate competitive inhibition constant 1.3 [dimensionless]
ku acetate uncompetitive inhibition constant 8.3E-3 [dimensionless]
kR acetate productive inhibition constant 12.6E-3 [mol∙L-1]
ks,L half-saturation constant of Monod equation taking lactate as substrate 12E-3 [mol∙L-1]
ks,A half-saturation constant of Monod equation taking aactate as substrate 12.3E-3 [mol∙L-1]
kMox half-constant of Monod equation taking oxidized mediator as substrate 6E-6 [mol∙L-1]
YX/L apparent biomass coefficient of lactate 19.1 [g dry cell/ mol lactate]
YX/A apparent biomass coefficient of acetate 16.8 [g dry cell/ mol acetate]
n the numbers of electron transfer 1 [mol e-/mol reaction]
kf oxidation rate constant 5.2504E-5 [h-1]
kb reduction rate constant 5.2574E-5 [h-1]
α oxidization transfer coefficient in anode 0.5 [dimensionless]
F faraday's constant 96485 [C∙mol-1]
R ideal gas constant 8.314 [J∙mol-1∙K-1]
T thermodynamic temperature 303.15 [K]
ΔE The potential between riboflavin and the heme left of OM c-Cyts 0.395 [v]
k0 standard rate constant 0.10 [h-1]
vmf the maximum reaction rate constant 6E-5 [h-1]
kiA the dissociation constant of NADH 1E-2 [dimensionless]
kmA Michaelis constant when NADH saturation 1E-2 [dimensionless]
kmB Michaelis constant when riboflavin saturation 1E-2 [dimensionless]
k1 the inhibition constants of riboflavin 5E-4 [dimensionless]
k2 the inhibition constants of NADH 10 [dimensionless]



3 Model for co-culture MFC



Growth Model of E.coli

In order to figure out the kinetics feature of growth, substrate consumption and accumulation of metabolite in lactate producing E.coli, we construct kinetics model of fermentation process based on carbon balance. First of all, the fermentation process of lactate in E.coli under anaerobic conditions has been simplified and we make several assumptions as following:

① Except for carbon sources which serves as restrictive substrates, the growth of E.coli will not be affected by other nutrition in the medium.

E.coli will use carbon sources to maintain its survival, division and proliferation while the accumulation of metabolites is also counted in.

③ The growth model of E.coli is equilibrium increase rather than structural increase, which means the ascending rate of biomass in E.coli increases as the same as the rate for cell weight accumulation.

④ We take lactate as the only one metabolites for accumulation and lactate increase is positively related to its growth rate. Besides, overmuch lactate will inhibit the growth of E.coli.


We further construct the carbon equilibrium calculations for anaerobic lactate fermentation in batches. Therefore, we figure out the following differential equations:


  (18)

  (19)

  (20)

Owing to the assumption of equilibrium increase in E.coli, the relations between specific growth rate of biomass and concentration of restrictive substrates can reflect in the Monod equations that we modeled. However, considering that the accumulation of lactate may inhibit bacteria growth, we introduce an extra inhibition factor for general Monod equation. Due to the positive relations between lactate accumulation and cell growth, the production of metabolites also share the same relation with biomass increase. In this way, Qp can be represented by using the function of μE.


  (21)

  (22)

  (23)


Simulation Results



Figure 6. The growth curve of E.coli.



Figure 7. The conc. of lactate and acetate varies with time.



As can be seen from the above figures, experimental results are very close to the simulation results of Ordinary Differential Equations, thus it is considered a successful model. Based on the growth model of E.coli, We constructed the co-culture MFC model.



Model For Co-culture MFC

We constructed the E.coli growth model and the monocultures Shewanella MFC model into a co-culture MFC model. In this model, E. coli is responsible for the decomposition of glucose to lactic acid, to provide a carbon source for the S.onedensis, S.onedensis responsible for lactate decomposition of carbon dioxide, transporting the electrons to the anode by the redox reaction of riboflavin. The simulation Results is shown as follow:



Figure 8. The Biomass of E.coli and S.onedensis in MFC varies with time.



Figure 9. The Voltage varies with time.



As can be seen in figure 8, when the amount of E. coli is small, due to lack of carbon supply S.onedensis begin to die, the situation changed after the rapid growth of E. coli. According to this result, we reached the following conclusions: In order to maintain the activity of Shewanella at the begining, The primary objective is to ensure the lactate production of E. coli. This conclusion has played an important role in the construction of three-bacterium MFC.



Parameters and variables


Parameter Discribtion Number [Unit]
k lactate inhibition constant 0.085 [L/mol]
kS,G Half-saturation constant taking glucose as substrate 14E-3 [L/mol]
kS,L Half-saturation constant taking lactate as substrate 20E-3 [L/mol]
μEG,max E.coli's maximal specific growth rate taking glucose as substrate 1.3 [h-1]
μE/L,max E.coli's maximal specific growth rate taking lactate as substrate 0.4 [h-1]
kL,C lactate competitive inhibition constant 1 [dimensionless]
kL,U lactate uncompetitive inhibition constant 200E-3 [dimensionless]
mE Endogenous metabolism rate constant of E.coli 0.006 [h-1]
kL,R Lactate productive inhibition constant 12.6E-3 [mol∙L-1]
v volume of reactor 140E-3 [L]
YXE/G apparent biomass coefficient of lactate 8.1 [g dry cell/ mol lactate]
YXE/L apparent biomass coefficient of lactate 5.8 [g dry cell/ mol lactate]
YL/G apparent biomass coefficient of lactate 1.5 [g dry cell/ mol lactate]



4 Model for P170 promoter



Introduction

P170 is an acid-induced promoter. When rcfB protein is activated by lactic acid, it can bind to P170 and make it sensitive to RNA polymerase. In order to detect and control the concentration of lactic acid, we use P170 and RFP gene, which shows the relationship between lactic acid and fluorescence.



Figure 10. Reaction network of genetic regulation circuit


From the reaction network We can extract follow equation sets From the reaction network:


   
   


Assumptions

1. When genes transcript and translate downward, the mRNA can bind with ribosome and then start translation subsequently, which can insure a constant amount during the whole process. Therefore, we use steady-state approximation on mRNA and mRNA-rib.

2. Because there is some acid buffer in the cell to transfer the produced lactate out of cell, pH will have a rather small influence on enzyme activity. Therefore, we assume that all enzymes activity are not affected by acid concentration, the rate constant keeps unchanged.

3. Concentration of protein will eventually keep constant, which is because expressed protein can have a negative effect on gene expression in return. Thus we assume that the concentration of rcfB protein is constant, and its magnitude is determined by the strength of the promoter upstream.



Equation Establishment

Genetic regulation circuit is shown in Figure 10. When rcfB protein is activated by lactic acid, it can bind to P170 and make it sensitive to RNA polymerase.


  (24)

  (25)

where KLa and K are constants and n is Hill coefficient The acP170 is the activated state, the P170 can be activated in two ways. One is to bind with the lactic acid dimerization and RcfB which can make it easy to be activated. Another one, we assume that a of P170 can be activated directly. Both the binding and unbinding P170 can be promoted while they have different rate constant.

We conclude that:


  (26)

  (27)

According to the equation above, we can consolidate that:


  (28)

Where km1 and km2 are integration parameters.


Next, by using law of mass action and steady state solution, we deduce the equation below.


  (29)

  (30)

Where the concentration of mRNA and mRNA-rib satisfies:


  (31)

  (32)

In the end, the concentration of RFP and the fluorescence can be written as:


  (33)

  (34)

Thus we construct the connection between fluorescence and lactic acid, from which we can know the concentration of lactic acid through the system.



Simulation Results

We set km1 =10 and use promoter J23100, then we get the following result



Figure 11. The result shows how P170,mRNA,mRNA-rib,fluorescence change with the lactic acid.




Fine tuning

From the equation, we can predict that:

1. the integration parameter km1 can characterize the relative level of leakage.

2. the larger the RcfB concentration is, the quicker the fluorescence will be.


At first we have the constitutive promoter fine-tuning. The BIT-China suggests us using different promoters in our construction. The performance of promoter P170 and rCfB play an important role. We choose to use different constitutive promoter J23103 J23109 J23107 and J23100 for the result (Figure 11). According to the figures, we can get the highest rate of change in fluorescence when we use J23100 promoter which is also the most sensitive one among all the promoter testing. Therefore, we choose J23100 promoter since it is convenient to represent lactate concentration by using fluorescence as the measurement. Particularly, when we have a stronger promoter, the amount of rcfB will be more during stable state. As a result gained from equation, the P170 can be activated in a faster speed and consequently, the changing rate of fluorescence will be greater as well.



Figure 12. Fluorescence changes with time when rcfB is constitutively activated by J23100, J23103, J23107 and J23109,respectively.



Next, we will explore the effect brought by leakage which means that P170 can activate and lead to fluorescence expression when lactate doesn’t exist. We have predicted that km1 is the index of leakage. km1 is higher, then we get more leakage shown up. Here we make k_m1 equals to 1, 10, 50 respectively and the result is shown in Figure 12. As predicted, the result is well suited to our calculation, because km1 is a constant unrelated to lactate concentration. It has a remarkable influence on the result when lactate amount is relatively small.
Then we will analyze the constants related to km1, which mainly contains rate k3 that is directly activated by P170 as well as proportion , which can both affect leakage. Our target is reducing the leakage expression, thus, we need to lower those constants as much as possible.



Figure 13. Fluorescence changes with time when k_m1 equals to 1, 10, 50. Difference has been show to illustrate the situation of different km1 under limited lactate concentration.




Parameters and variables




</div>