Modeling in synthetic biology is a crucial tool that helps us to get a comprehensive vision of biological systems and the working principles while providing ways to improve the system. In our case, it significantly led to the design of our project while improving the part choices. In our project we aimed to improve the second generation bioethanol production in which pretreatment process are constructed, from microbial fermentation by increasing the lifespan and productivity of bacteria. In our project, we used E.coli KO11 which is an ethanologenic bacteria thus appropriate for our goal. However, because we designed our gene circuit not just for KO11 but also for other E.coli strains, the modeling was constructed while considering the overall properties of them. While doing so, we used kinetic and visual models to design our gene circuit, calculating the expected behaviours and demonstrating our system.

What we have achieved

1- We constructed our gene circuit with the help of toxicity analysis and enzymatic reaction kinetics.
2- We used microbial growth and fermentation kinetics to simulate the expected behaviors of our system and the effects of our genes with the data obtained from our wet lab team.
3- We improved the understanding of our project by demonstrating the pathways and effects of our genes .

Formation of Our Kinetic Models

Cell Growth Kinetics

We decided to use Monod equation for the cell growth to demonstrate the effects of our genes in which the cell concentration over time depends on the number of cells (X) and the substrate concentration (S) along with several parameters such as Monod Constant (Ks) and maximum specific growth rate (μmax).

Figure 1: Monod Equation

Enzymatic Reaction Kinetics

We decided to use Michaelis and Menten kinetics for the modeling of enzymatic reactions. Thanks to the kinetic model it was possible to estimate the affinity of our genes and the rate of reaction that led us to determine the proper genes for our gene circuit.

Figure 2: Michaelis and Menten Kinetics, where V = reaction rate, S = substrate concentration, Km = half saturation constant and Vmax = maximum reaction rate.

Fermentation Kinetics

Ethanol production rate depends on the cell and substrate concentration. We also considered the mutual inhibition of xylose and glucose fermentation along with the alcohol inhibition. Furthermore, because even when the glucose is consumed completely, the xylose utilization is slower, the phenomenon is shown with a constant (Olsson, Hagerdalt & Zacchi, 1995). Thus, we used the formula developed by Olsson, Hagerdalt & Zacchi (1995) in which the xylose and glucose fermentation was calculated separately and then simultaneously.

G = Glucose concentration (g/L)
K1 = inhibition constant for xylose fermentation
K2 = inhibition constant for nonspecific inhibitors in the condensate in glucose fermentation
K3 = inhibition constant for nonspecific inhibitors in condensate in xylose fermentation
Kms = Monod constant for pentose fermentation (g/L)
Kigs = inhibition constant for pentoses in glucose fermentation (g/L)
Kisg = inhibition constant for glucose in pentose fermentation (g/L)
X = final cell mass concentration
P = product (ethanol) concentration
Rps = ethanol production rate from xylose g/L * h)
Rpg = ethanol production rate from glucose (g/L * h)
Rp = ethanol production rate from both glucose and xylose (g/L *h)
Vmg = maximum ethanol production rate in glucose fermentation(h^-1)
Vms = maximum ethanol production rate in pentose fermentation(h^-1)
S = pentose concentration
Pm = highest amount of ethanol observed

Toxicity Analysis:

Our main issue was determining the genes that would have a positive effect on bacterial fermentation. Thus, we analyzed the toxicity of different inhibitory molecules and substances produced during pretreatment: Furfural, 5-hydroxymethylfurfural (5-HMF) and aliphatic acids; such as formic acid, and levulinic acid.

Toxicity Analysis Results:

Formic acid and Levulinic acid were shown to inhibit the growth significantly but furans were worse (Kim et al., 2013). Hence, we chose to focus on furans. When we discard the lag caused by the inhibitors and focus on the last state, it was possible to calculate what percentage the cell growth was inhibited. Furans (Furfural and HMF) showed approximately 80% inhibition which is the highest number among thus we decided to increase the tolerance of E.coli to furans.

Furans (HMF & Furfural):

According to the analysis, furans were found to be the most toxic substances when it comes to cell growth (Kim et al., 2013). The most used pretreatment process, dilute acid, gives rise to the formation of furanic aldehydes (Palmqvist and Hahn-Hagerdal, 2000; Larsson et al., 1999; Thomsen et al., 2009; Klinke et al., 2004). They are highly reactive, contributing to the birth of reactive oxygen species (ROS) which damage proteins, nucleic acids and cell organelles (Wierckx et al., 2011). Because of the toxicity provided by furans, cell mass and productivity of fermentation decreases (Almeida et al., 2009; Palmqvist and Hahn-Hagerdal, 2000b; Thomsen et al., 2009). Thus, we examined the pathways of furfural and HMF to find a way to eliminate the setbacks and increase the tolerance.

Reaction Pathways:

It was shown that furfural can be reduced to a less toxic form, furfuryl alcohol, by NAD(P)H dependent oxidoreductases which are transcribed by FucO and YqhD genes (Wierckx et al., 2011). The pathways are shown below:

Furfural Reduction Pathway of YqhD

Enzyme = NADPH Dependent Oxidoreductase:

Furfural + NADPH + Enzyme → NADP+ + Enzyme + Furfuryl Alcohol

HMF + NADPH + Enzyme → NADP+ + Enzyme + 5-Hydroxymethyl-2-furfuryl alcohol

Furfural Reduction Pathway of FucO

Enzyme = NADH dependent Oxidoreductase:

Furfural + NADH + Enzyme → NAD+ + Enzyme + Furfuryl Alcohol

HMF + NADH + Enzyme → NAD+ + Enzyme + 5-Hydroxymethyl-2-furfuryl alcohol

Reaction Kinetics:

The comparison between the reaction rate and Km values of furfural and HMF reductions with different oxidoreductases were done by Michaelis-Menten kinetics and Lineweaver-Burk plot with the data obtained from Miller et al. (2009) and Wang et al. (2011). Matlab’s enzkin function was used to evaluate the results.

Enzyme Reaction Kinetics of FucO

Enzyme Reaction Kinetics of YqhD

Results

It was important to find the appropriate NAD(P)H dependent oxidoreductase that would decrease the harmful effects of furans thus resulting in the improvement on rate of cell mass and bioethanol production (Jarboe et al., 2012; Wang et al., 2011). Therefore we analyzed the reaction kinetics of both FucO (NADH dependent) and YqhD (NADP dependent) with the Michaelis and Menten enzyme kinetics and Lineweaver - Burk plot. At first, we looked through the Km values because they indicate the affinity of enzymes which means that if you have a low Km value then the enzyme is more likely to catalyze the reactions faster and properly (Jarboe et al., 2012). YqhD showed a Km of 5.00 +- 3 mM where Km of FucO was 0.4+- 0.2 mM (Wang et al., 2011). It was shown that that FucO has higher affinity to furfural and is more likely to increase the furfural tolerance (Jarboe et al., 2012 ; Wang et al., 2011). Moreover, YqhD has higher Km for NADPH than most of the key metabolic enzymes such as CysJ (80 μM), which is necessary for sulfate assimilation to form cysteine and methionine; ThrA (90 μM), is important for the formation of threonine; and DapB (17 μM), required for lysine formation (Miller et al.,2009; Jarboe et al., 2012). Therefore, the utilization of YqhD inhibits the growth of the bacteria due to the competition with the important biosynthetic enzymes (Miller et al.,2009). Though, YqhD is found in most of the E.coli strains, due to its lower affinity compared to the FucO, it is possible to eliminate the YqhD gene by the overexpression of FucO (Jarboe et al., 2012). Thus we decided to use the FucO gene coding for L-1,2-propanediol oxidoreductase that is responsible for the furanic compound degradation. Moreover, because furans’ high reactiveness eventually leads to the formation of ROS, the GSH gene producing glutathione synthetase was decided to be utilized in order to decrease the harmful effects and raise tolerance to environmental toxicity.

GSH - Glutathione Synthetase

GSH gene codes for Bifunctional gamma-glutamate-cysteine ligase/glutathione synthetase which is responsible for the mass production of the main antioxidant, glutathione. It is one of the most important antioxidant that works in order to reduce ROSs that are produced during metabolic activities such as furfural reduction and fermentation (Pizzorno, 2014; Forman, 2009). It prevents oxidative stress from building up in the cell metabolism which causes cell damage and eventually, death (Pizzorno, 2014; Forman, 2009). The pathway of glutathione synthetase and glutathione are demonstrated below along with the effects of glutathione on cell mass:

Glutathione Synthesis by Glutathione Synthetase

Glycine + y-glutamyl cysteine + Glutathione Synthetases → Glutathione Synthetase Glutathione

Glutathione Oxidation and Reduction Pathways

Glutathione + ROS → Water + Oxidized Glutathione

Oxidized Glutathione + Glutathione Peroxidases + NADPH → NADP+ + Glutathione

Improving Effects of Glutathione(GSH)

Then, we looked through the effects of glutathione on cell mass in order to observe the improving features of glutathione with the data obtained from Kim & Hahn (2013). In brief even 2mM difference make 30% increase in cell mass which is demonstrated below(Kim & Hahn, 2013).

Reaction Kinetics

Jez, J.M. & Cahoon, R.E (2004) stated that the km values of Glutathione Synthetase(GS) were 39 ± 5 μm for y-glutamyl cysteine and 1510 ± 88 μm for glycine meaning that GS has low affinity to glycine resulting in problematic production of glutathione and low reaction rate. In order to increase the reaction rate and glutathione formation we decided to enhance the gene expression rate by using a strong promoter and RBS. Thus, we looked through the existing promoters and RBS in iGEM database and examined their strengths which were compared by LMU-Munich 2012 iGEM team and the Alverno_Ca team.

Promoter Analysis

RBS Analysis

Decision

Thanks to the modeling of different promoter and RBS strengths, we decided to use the promoter J23100 because of it being the most powerful promoter among the Anderson promoter family and B0034 as RBS because of it being one of the most strongest RBS and the encouraging comments. Then, we decided to construct three different modules to improve the lifespan and ethanol production of E.coli KO11. Our first module consists of a strong promoter (J23100), an RBS(B0034), FucO gene only and a double terminator (B0015). Our second module consists of a strong promoter (J23100), an RBS(B0034), GSH gene only and a double terminator (B0015). Our last module is named Bio-E because it contains both genes with the same promoter, RBS and double terminator. The modules are tested to observe the lifespan against un-engineered E.coli KO11 then, simulated and compared by Monod equation with the data obtained from our wet lab team. Moreover, we simulated the expected ethanol production results by fermentation kinetics with the estimated parameters.

Growth and Fermentation

Because fermentation mainly depends on the cell growth of bacteria, there are no striking differences between cell growth and ethanol production (Luong, 1985). We demonstrated the kinetics of cell growth and fermentation with the alcohol and furans inhibition. We tested four modeles: Un-engineered E.coli strain KO11, with the overexpression of GSH, with the overexpression of FucO and with the overexpression of both GSH and FucO. The equations and results are shown below:

Growth

The growth rate of our bacteria is important to determine the effects of our genes because we aimed to increase the toxicity tolerance of E.coli KO11 and consequently leading to a better cell growth rate and lifespan. The specific growth rate and the substrate concentration data were used to define the maximum specific growth rate and Monod Constant.

In order to obtain the parameters, Monod growth model is linearized by inverting and factoring out umax and we get the Lineweaver-Burk plot (Rorke & Kana, 2017 ).

Furthermore, instead of just a point estimate of the fit, we wanted study the predictive posterior distribution of the model. By them we calculated the model fit for a randomly selected subset of the chain and calculate the predictive envelope of the model. The grey areas in the plot correspond to 50%, 90%, 95%, and 99% posterior regions with the data obtained from Olsson, Hagerdalt & Zacchi (1995) for un-engineered E.coli KO11.

After determining the parameters, we used Monod equation with the fit data to demonstrate the cell growth.

Because our experiments were incomplete and we were not able to measure the change in substrate concentration, only the specific growth rate of our bacteria were calculated with the data obtained from our wet lab team: KO11, KO11 with FucO, KO11 with GSH, KO11 with GSH and FucO. They measured the change in cell concentration with OD600 in 10mM furfural. Then the change in the specific growth rate were compared in percent. The results were used to redesign the Monod model that are shown below.

Result

Insertion of FucO gene increased the cell growth rate by 10 %. Insertion GSH gene increased the cell growth rate by 18% percent. Insertion of GSH and FucO gene increased the growth rate by 120 %. Thus, the best improvement was observed with the dual expression of FucO and GSH gene and demonstrated by Monod equation.

Fermentation

Ethanol production rate depends on the cell mass, sugar concentration and the inhibitor substances (Olsson, Hagerdalt & Zacchi, 1995). In our project, the improvement of ethanol production was decided to be accomplished by increasing the tolerance of bacteria to inhibitors. Thus, the kinetic module developed by Olsson, Hagerdalt & Zacchi (1995) was used to emphasize the expected ethanol production rate of normal and improved E.coli KO11. Because we couldn’t construct a fermentation experiment, we used the parameters’ value from Olsson, Hagerdalt & Zacchi (1995) which were obtained from the utilization of KO11 in toxic substances. The Monod equation which was redesigned by our experimental data, was integrated into the fermentation kinetic model. The kinetic model were constructed by Matlab and shown below.

Result:

According to our kinetic models, the best ethanol yield was shown with the dual expression of FucO and GSH gene and demonstrated by fermentation model while portraying the expected results.

References

• Almeida, J.R.M., Modig, T., Petersson, A., Hahn-Hagerdal, B., Liden, G., Gorwa-Grauslund, M.F. (2007). Increased tolerance and conversion of inhibitors in lignocellulosic hydrolysates by Saccharomyces cerevisiae. J Chem Tech Biotechnol 82(4):340–349. https://doi.org/10.1002/jctb.1676
• Han, K., & Levenspiel, O. (1988). Extended monod kinetics for substrate, product, and cell inhibition. Biotechnology and Bioengineering, 32(4), 430–447. doi: 10.1002/bit.260320404 Heer, D., Sauer, U. (2008). Identification of furfural as a key toxin in lignocellulosic hydrolysates and evolution of a tolerant yeast strain. Microbial Biotechnol 1(6):497–506. doi: 10.1111/j.1751-7915.2008.00050.x
• Klinke, H.B., Thomsen, A.B., Ahring, B.K. (2004). Inhibition of ethanol-producing yeast and bacteria by degradation products produced during pre-treatment of biomass. Appl Microbiol Biotechnol 66(1):10–26. doi: 10.1007/s00253-004-1642-2
• Larsson, S., Palmqvist, E., Hahn-Hagerdal, B., Tengborg, C., Stenberg, K., Zacchi, G., Nilvebrant N.O. (1999). The generation of fermentation inhibitors during dilute acid hydrolysis of softwood. Enzyme Microb Technol 24(3–4):151–159. https://doi.org/10.1016/S0141-0229(98)00101-X
• Miller, E. N., Jarboe, L. R., Yomano, L. P., York, S. W., Shanmugam, K. T., & Ingram, L. O. (2009). Silencing of NADPH-Dependent Oxidoreductase Genes (yqhD and dkgA) in Furfural-Resistant Ethanologenic Escherichia coli. Applied and Environmental Microbiology, 75(13), 4315–4323. http://doi.org/10.1128/AEM.00567-09
• Palmqvist, E., Hahn-Hagerdal, B. (2000a) Fermentation of lignocellulosic hydrolysates. I: inhibition and detoxification. Bioresource Technol 74:17–24. https://doi.org/10.1016/S0960-8524(99)00160-1
• Palmqvist, E., Hahn-Hagerdal, B. (2000b) Fermentation of lignocellulosic hydrolysates. II: inhibitors and mechanisms of inhibition. Bioresource Technol 74:25–33. https://doi.org/10.1016/S0960-8524(99)00161-3
• Thomsen, M.H., Thygesen, A., Thomsen, A.B. (2009) Identification and characterization of fermentation inhibitors formed during hydrothermal treatment and following SSF of wheat straw. Appl Microbiol Biotechnol 83(3):447–455. doi: 10.1007/s00253-009-1867-1
• Zaldivar, J., Martinez, A., Ingram, L.O. (1999).Effect of selected aldehydes on the growth and fermentation of ethanologenic Escherichia coli. Biotechnol Bioeng 65(1):24–33. https://doi.org/10.1002/(SICI)1097-0290(19991005)65:1<24::AID-BIT4>3.0.CO;2-2
• Luong, J.H.(1985).Kinetics of Ethanol Inhibition in Alcohol Fermentation. Biotechnol Bioeng 27, 280-285. doi:10.1002/bit.260270311
• Rorke, D.C.S., Kana E.B.G.(2017).Kinetics of Bioethanol Production from Waste Sorghum Leaves Using Saccharomyces cerevisiae BY473. Fermentation 3(2), 19. doi: https://doi.org/10.3390/fermentation3020019
• Wierckx, N., Koopman, F., Ruijssenaars, H. J., & de Winde, J. H. (2011). Microbial degradation of furanic compounds: biochemistry, genetics, and impact. Applied Microbiology and Biotechnology, 92(6), 1095–1105. http://doi.org/10.1007/s00253-011-3632-5
• Wang, X., Miller, E. N., Yomano, L. P., Zhang, X., Shanmugam, K. T., & Ingram, L. O. (2011). Increased Furfural Tolerance Due to Overexpression of NADH-Dependent Oxidoreductase FucO in Escherichia coli Strains Engineered for the Production of Ethanol and Lactate. Applied and Environmental Microbiology, 77(15), 5132–5140. http://doi.org/10.1128/AEM.05008-11
• Kim, S. K., Park, D.H., Song, S.H., Wee, Y.J., Jeong, G.T. (2013) Effect of fermentation inhibitors in the presence and absence of activated charcoal on the growth of Saccharomyces cerevisiae, Bioprocess Biosyst Eng. 36(6), 659-666. doi: 10.1007/s00449-013-0888-4
• Jez, J.M., Cahoon, R.E. (2004). Kinetic mechanism of glutathione synthetase from Arabidopsis thaliana. J Biol Chem. 279(41), 42726-42731. doi: 10.1074/jbc.M407961200
• Jarboe, L. R., Liu, P., Kautharapu, K. B., & Ingram, L. O. (2012). Optimization of enzyme parameters for fermentative production of biorenewable fuels and chemicals. Computational and Structural Biotechnology Journal, 3, e201210005. http://doi.org/10.5936/csbj.201210005
• Kim, D., & Hahn, J.,S. (2013). Roles of the Yap1 Transcription Factor and Antioxidants in Saccharomyces cerevisiae’s Tolerance to Furfural and 5-Hydroxymethylfurfural, Which Function as Thiol-Reactive Electrophiles Generating Oxidative Stress. Applied and Environmental Microbiology, 79(16), 5069–5077. http://doi.org/10.1128/AEM.00643-13
• Olsson, L., Hagerdalt, B. & Zacchi, G. (1995). Kinetics of ethanol production by recombinant Escherichia coli KO11. Biotechnol Bioeng, 20;45(4) ,356-65. doi: 10.1002/bit.260450410