Difference between revisions of "Team:METU HS Ankara/Model"
| Line 16: | Line 16: | ||
</header> | </header> | ||
| − | + | ||
| − | + | ||
| − | + | ||
<script src="https://2018.igem.org/common/MathJax-2.5-latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"/> | <script src="https://2018.igem.org/common/MathJax-2.5-latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"/> | ||
Revision as of 10:51, 16 October 2018
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 modelling 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.
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).
We decided to use Michaelis and Menten kinetics for 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.
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 K1 (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)
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.
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.
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.
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:
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.
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 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:
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 use a strong promoter and RBS that would increase the gene expression. 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.
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.
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:
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 ).
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 ethanol production rate of normal and improved E.coli KO11.
An empirical kinetic model developed by Olsson, Hagerdalt & Zacchi (1995) was used to emphasize the ethanol production rate of normal and improved E.coli KO11. The model that uses the Monod kinetics including substrate and product inhibition, is demonstrated below.
Modelling
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
Enzymatic Reaction Kinetics
Fermentation Kinetics
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:
Figure 3:
The data were obtained from Kim et al. (2013). According to Kim et al. (2013), inhibitor substances were put individually into the medium in 5g/L and relative cell growths were calculated. In the presence of furfural and HMF, relative growth was lower than in the presence of formic acid, and levulinic acid (Kim et al., 2013).
Toxicity Analysis Results
Furans (HMF & Furfural)
Reaction Pathways
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:
Enzyme Reaction Kinetics of FucO
Enzyme Reaction Kinetics of YqhD
Results
GSH - Glutathione Synthetase
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
Reaction Kinetics
Promoter Analysis
RBS Analysis
Decision
Growth and Fermentation
Growth
Fermentation