Team:Bielefeld-CeBiTec/Results/Modelling

From 2014.igem.org

(Difference between revisions)
(new version)
(done)
 
(7 intermediate revisions not shown)
Line 19: Line 19:
-
<h1>Modelling</h1>
+
<h1>Modeling</h1>
-
 
+
<div class="element">
-
 
+
<div class="element" style="margin:10px; padding:10px">
<div class="element" style="margin:10px; padding:10px">
<div id="text">  
<div id="text">  
-
<br><br>
+
<h6>Summary</h6>
-
<h2>Abstract</h2>
+
-
<div id="text">
+
   <p>
   <p>
-
Different modelling approaches were used to identify enzymatic bottlenecks in the isobutanol production pathways and to predict the formation of the desired product in a given time. First of all we created a map of all metabolic reactions which are part of our project (figure 1). This network not only provides a good overview, it also serves as the basic tool for further considerations. Due to the huge amount of components it does not seem feasible to create a computational model for all reactions at once. Therefore we started our modelling work by carrying out a stochiometric analysis. It reveals that 42 electrons are needed for the production of one isobutanol molecule from CO<sub>2</sub>. As for the isobutanol production pathway, dynamic modeling was carried out, in the course of which bottlenecks could be identified. An increase in expression von ilvD and kivD could increase the isobutanol production after four hours by 100%. Finally we prepared the extension of the existing model to predict the effect of specific carbon dioxide fixing reactions.
+
Different modeling approaches were used to identify enzymatic bottlenecks in the isobutanol production pathway and to predict the formation of the desired product in a given time. First of all we created a map of all metabolic reactions which are part of our project (Figure 1). This network not only provides a good overview, it also serves as the basic tool for further considerations. Due to the huge amount of components it does not seem feasible to create a computational model for all reactions at once. Therefore we started our modeling work by carrying out a stochiometric analysis. It revealed that 42 electrons are needed for the production of one isobutanol molecule from CO<sub>2</sub>. As for the isobutanol production pathway, dynamic modeling was carried out, in the course of which several bottlenecks could be identified. An increase in expression of <i>kivD</i> and <i>adhA</i> could boost the isobutanol production after four hours by 100&nbsp;%. Finally we prepared the extension of the existing model to predict the effect of specific carbon dioxide fixing reactions.
</p>
</p>
</div>
</div>
</div>
</div>
-
</div>
 
-
<br>
 
Line 40: Line 35:
<div class="element" style="margin:10px; padding:10px">
<div class="element" style="margin:10px; padding:10px">
<div id="text">  
<div id="text">  
-
<br><br>
+
<h6>Introduction</h6>
-
<h2>Introduction</h2>
+
   <p>
   <p>
-
Mathematical modelling is crucial to understand complex biological systems (<a href="#schaber2009">Schaber et al., 2009</a>). The analysis of isolated biological components hass been supplemented by a systems biology approach in over ten years (<a href="#chuang2010">Chuang et al., 2010</a>). Mathematical modelling is used to combine biological results (<a href="#kherlopian2008">Kherlopian et al., 2008</a>). Modeling is also a way to achieve results without carrying out experiments in a laboratory. The behaviour of a system can be simulated to get results, which cannot be derived from simply looking at the given system (<a href="#schaber2009">Schaber et al., 2009</a>).
+
Mathematical modeling is crucial to understand complex biological systems (<a href="#schaber2009">Schaber et al., 2009</a>). The analysis of isolated biological components has been supplemented by a systems biology approach in over ten years (<a href="#chuang2010">Chuang et al., 2010</a>). Furthermore mathematical modeling is used to combine biological results (<a href="#kherlopian2008">Kherlopian et al., 2008</a>). Modeling is also a way to achieve results without carrying out experiments in a laboratory. The behaviour of a system can be simulated to get results which cannot be derived from simply looking at the given system (<a href="#schaber2009">Schaber et al., 2009</a>).
-
The most important aim of any modelling approach is the reduction of complexity. As the given biological reality is often diverse and variable, it is important to identify the major rules and principles, which can describe a system.
+
The most important aim of any modeling approach is the reduction of complexity. As the given biological reality is often diverse and variable, it is important to identify the major rules and principles which can describe a system.
</p>
</p>
</div>
</div>
Line 52: Line 46:
<div class="element" style="margin:10px; padding:10px">
<div class="element" style="margin:10px; padding:10px">
<div id="text">  
<div id="text">  
-
<br><br>
 
-
<h2>Our aims</h2>  
+
<h6>Our aims</h6>  
<ul>
<ul>
-
<li>Visualization of the complete meatoblic network </li>
+
<li>Visualization of the complete metabolic network </li>
-
<li>Relating electron input to product output </li>
+
<li>Relating electron input to product formation </li>
<li>Relating electron input to carbon dioxide fixation </li>
<li>Relating electron input to carbon dioxide fixation </li>
<li>The identification of bottlenecks in the isobutanol production pathway </li>
<li>The identification of bottlenecks in the isobutanol production pathway </li>
Line 64: Line 57:
</div>
</div>
</div>
</div>
-
 
+
</div>
-
<br><br>
+
<div class="element" style="margin:10px; padding:10px">
<div class="element" style="margin:10px; padding:10px">
<div id="text">  
<div id="text">  
-
<br><br>
+
<h6>Stoichiometric analysis</h6>
-
<h2>Stoichiometric analysis</h2>
+
<p>
<p>
-
Stoichiometric analysis are useful to get information about the maximal output of a system. In this project the numer of electrons moving into the system limits the amount of product, which could be synhetized in the cells. Therefore the stoichiometric relations of all substances involved in our complex reaction network were calculated (figure 1). The calculation starts with the electrons, which are transportet into the system by mediators. We calculated the resulting production of intracellular molecules based on our map of the metabolic system (figure 1). The results are listed below.
+
Stoichiometric analysis is useful to get information about the maximal output of a system. In this project the number of electrons moving into the system limits the amount of the product which could be synthetized in the cells. Therefore the stoichiometric relations of all substances involved in our complex reaction network were calculated (Figure 1). The calculation starts with the electrons, which are transported into the system by mediators. We calculated the resulting production of intracellular molecules based on our map of the metabolic network (Figure 1). The results are listed below.
</p>
</p>
-
 
-
<br><br>
 
Line 86: Line 75:
</div>
</div>
</center>
</center>
 +
<br>
-
<p>
 
-
The theoretical electron costs of different molecules is listed in table 1. All calculations are based on our pathway map. According to our pathway map there are 42 electrons needed for the production of one molecule isobutanol, if CO<sub>2</sub> is used as sole carbon source. Our calculation does not involve the house keeping metabolism of <i>E. coli</i>, which consumes lots of energy for the survival of the cell. The true number of consumed electrons per produced isobutanol molecule is therefore much higher. The applied electric power can be converted into a number of electrons by the following equation: 1 A = 1 C * s<sup>-1</sup> = 6.2415065 * 10<sup>18</sup> electrons.
 
-
</p>
 
-
<br><br>
 
-
<b>Table1:</b> This table shows the theoretical electron cost of different intracellular molecules. The electron cost is the number of electrons which are needed for the synthesis of this metabolite. The intermediates are substances which are needed for the production of the final product. There costs are included in the final value for each substance.
 
-
<table width="100%" border="1" cellpadding="5" style="background-color:transparent">
+
<div class="element" style="width:400px; float:left">
 +
<font size="2"><b>Table 1:</b> This table shows the theoretical electron cost of different intracellular molecules. The electron cost is the number of electrons which are needed for the synthesis of this metabolite. The intermediates are substances which are needed for the production of the final product. Their costs are included in the final value for each substance.</font>
 +
 
 +
<table width="400px" border="1" cellpadding="5" style="background-color:transparent;float:left">
<tr>
<tr>
<th>Substance</th>
<th>Substance</th>
Line 108: Line 96:
<tr> <td>Isobutanol</td> <td>2x Pyruvate</td> <td>42</td></tr>
<tr> <td>Isobutanol</td> <td>2x Pyruvate</td> <td>42</td></tr>
</table>
</table>
 +
</div>
 +
 +
<p>
 +
<br>
 +
The theoretical electron costs of different molecules is listed in Table 1. All calculations are based on our pathway map. According to our pathway map there are 42 electrons needed for the production of one molecule isobutanol, if CO<sub>2</sub> is used as sole carbon source. Our calculation does not involve the house keeping metabolism of <i>E. coli</i>, which consumes lots of energy for the survival of the cell. The true number of consumed electrons per produced isobutanol molecule is therefore much higher. The applied electric power can be converted into a number of electrons by the following equation: 1 A = 1 C * s<sup>-1</sup> = 6.2415065 * 10<sup>18</sup> electrons.
 +
</p>
 +
<br><br><br><br><br><br><br><br><br><br>
</div>
</div>
</div>
</div>
-
<br><br>
 
<div class="element" style="margin:10px; padding:10px">
<div class="element" style="margin:10px; padding:10px">
<div id="text">
<div id="text">
-
 
+
<h6> Dynamic modeling </h6>
-
<br><br>
+
-
<h2> Dynamic modelling </h2>
+
<p>
<p>
-
After the stoichiometric analysis of the system we designed a dynamic model containing kinetic equations for the metabolite concentrations. It allows the identification of metabolic or enzymatic bottlenecks. This is our major aim. We could even use this information in the next step to modify constructs e.g. exchange an RBS or a promotor sequence. This could improve the different enzyme concentrations for upcoming experiments in the laboratory. Furthermore it was our target to predict the production of isobutanol per substrate in a given time. Our model predicts the isobutanol production in a carbon dioxide fixing cell. To achieve our aims we reduced the complex system shown in figure 1 to the version shown in figure 2. This metabolic network was suitable for our dynamic modelling approach.
+
After the stoichiometric analysis of the system we designed a dynamic model containing kinetic equations for the metabolite concentrations. It allows the identification of metabolic or enzymatic bottlenecks, the major aim of our modeling approach. This information might then be used in the next step to modify constructs e.g. exchange an RBS or a promotor sequence. This could help to balance the different enzyme concentrations in upcoming experiments in the laboratory, it was our goal to predict the production of isobutanol per substrate in a given time. Our model predicts the isobutanol production in a carbon dioxide fixing cell. To achieve our aims we reduced the complex system shown in Figure 1 to the version shown in Figure 2. This metabolic network was suitable for our dynamic modeling approach.
</p>
</p>
-
<br><br>
+
<br>
<center>
<center>
<div class="element" style="width:750px">  
<div class="element" style="width:750px">  
       <a href="https://static.igem.org/mediawiki/2014/b/ba/Bielefeld-CeBiTec_14-10-14_reduced_metabolic_system.jpg" target="_blank"><img src="https://static.igem.org/mediawiki/2014/b/ba/Bielefeld-CeBiTec_14-10-14_reduced_metabolic_system.jpg" width="750px"></a><br>
       <a href="https://static.igem.org/mediawiki/2014/b/ba/Bielefeld-CeBiTec_14-10-14_reduced_metabolic_system.jpg" target="_blank"><img src="https://static.igem.org/mediawiki/2014/b/ba/Bielefeld-CeBiTec_14-10-14_reduced_metabolic_system.jpg" width="750px"></a><br>
       <font size="2" style="text-align:left;">
       <font size="2" style="text-align:left;">
-
<b>Figure 2</b>: Reduced metabolic network of reactions which were selected for modelling. </font>
+
<b>Figure 2</b>: Reduced metabolic network of reactions which were selected for modeling. </font>
</div>
</div>
</center>
</center>
-
<br><br>
+
<br>
-
</div>
+
-
</div>
+
-
 
-
 
-
<div class="element" style="margin:10px; padding:10px">
 
-
<div id="text">
 
-
<h3>Dynamic modelling of the isobutanol production pathway</h3>
 
<p>
<p>
-
The modelling work on the <a href="https://2014.igem.org/Team:Bielefeld-CeBiTec/Project/Isobutanol/GeneticalApproach">isobutanol pathway</a> is based on publications about the isobutanol production pathway <a href="#atsumi2008">(Atsumi et al., 2008</a> and <a href="#atsumi2010"> Atsumi et al., 2010)</a>. We started our work on the isobutanol production pathway by collecting the appropriated kinetic parameters. They were used for the development of a  system of differential equations. As for the choice of the kinetics used, we stick to Michaelis-Menten kinetics. This was published as the best approach, if reaction kinetics are not known (<a href="#breitling2008">Breitling et al., 2008</a>; <a href="#chubukov2014">Chubukov et al., 2014</a>). Additionally kcat and KM values were collected from databases like <a href=“http://www.genome.jp/kegg/“>KEGG</a>, <a href=“http://biocyc.org/“>biocyc</a> and <a href=“http://www.brenda-enzymes.org/“>BRENDA </a>(table 2). Missing values had to be estimated. The starting concentrations for different metabolites were also taken from the literature and from different databases (table 3). <br>
+
The modeling work on the <a href="https://2014.igem.org/Team:Bielefeld-CeBiTec/Project/Isobutanol/GeneticalApproach" target="_blank">isobutanol pathway</a> is based on publications about the isobutanol production pathway <a href="#atsumi2008">(Atsumi et al., 2008</a> and <a href="#atsumi2010"> Atsumi et al., 2010)</a>. We started our work on the isobutanol production pathway by collecting the appropriate kinetic parameters. They were used for the development of a  system of differential equations. Concerning the choice of the kinetics used, we stuck to Michaelis-Menten kinetics as this was published as the best approach if reaction kinetics are not known (<a href="#breitling2008">Breitling et al., 2008</a>; <a href="#chubukov2014">Chubukov et al., 2014</a>). Additionally k<sub>cat</sub> and K<sub>M</sub> values were collected from databases like <a href="http://www.genome.jp/kegg" target="_blank">KEGG</a>, <a href="http://biocyc.org" target="_blank">biocyc</a> and <a href="http://www.brenda-enzymes.org" target="_blank">BRENDA </a>(Table 2). Missing values had to be estimated. The starting concentrations for different metabolites were also taken from the literature and from different databases (Table 3).
-
We decided to use only the forward reactions from pyruvate to isobutanol for different reasons. First it is necessary to get the maximal product concentration, second the reactions are neary irreversibel due to the decarboxylation steps and third there are only a few information on the back reactions available.
+
<br>
 +
We decided to use only the forward reactions from pyruvate to isobutanol for different reasons. First of all, it is necessary to get the maximal product concentration, secondly the reactions are nearly irreversible due to the decarboxylation steps and thirdly there is only few information on the back reactions available.
</p>
</p>
-
<br><br>
+
<div class="element">
-
<b>Table2:</b> This table shows all enzymatic parameters which were used for our dynamic model.
+
<font size="2"><b>Table 2:</b> This table shows all enzymatic parameters which were used for our dynamic model.</font>
Line 155: Line 141:
<tr>
<tr>
<th>Enzyme</th>
<th>Enzyme</th>
-
<th>k<sub>cat</sub></th>
+
<th>k<sub>cat</sub> [-]</th>
<th>K<sub>M</sub> [mM]</th>
<th>K<sub>M</sub> [mM]</th>
<th>Reference</th>
<th>Reference</th>
Line 181: Line 167:
</table>
</table>
-
<br><br>
+
</div>
 +
<br>
-
 
+
<div class="element">
-
<br><br>
+
<font size="2"><b>Table 3:</b> This table shows all metabolite concentrations which were used for our first model. The metabolite concentration was set to zero, if no published value was available.</font>
-
<b>Table3:</b> This table shows all metabolite concentrations which were used for our first model. The metabolite concentration was set to zero, if no published value was available.
+
<table width="100%" border="1" cellpadding="5" style="background-color:transparent">
<table width="100%" border="1" cellpadding="5" style="background-color:transparent">
Line 212: Line 198:
<tr>
<tr>
-
<td>Isobutyraldehyde</td> <td>0.6</td> <td></td>
+
<td>Isobutyraldehyde</td> <td>0.6 (estimated)</td> <td></td>
</tr>
</tr>
Line 219: Line 205:
</tr>
</tr>
</table>
</table>
-
<br><br>
+
</div>
 +
<br>
 +
 
 +
 
<div id="text">
<div id="text">
<p>
<p>
-
We implemented the system of differential equations (figure 3) in matlab (<a href="https://2014.igem.org/Team:Bielefeld-CeBiTec/Results/Modelling/sourcecode">source code</a>). The predicted changing of the metabolic concentration over the time is shown in figures 4-7. The amount of expressed proteins could differ depending on the distance of the coding sequence downstream of the promotor. The coding sequences for the involved enzymes are located downstream of a common promotor. Therefore we decided to set the enzyme concentration for the first enzyme to 1 and decrease in steps of 0.1 (<a href="https://2014.igem.org/Team:Bielefeld-CeBiTec/Results/Modelling/sourcecode">source code</a>). These different values were tried to identify appropiate concentrations for each enzyme. The results of this dynamic modelling approach could be transfered to the laboratory by using promotors of different strength.
+
We implemented the system of differential equations in Matlab (<a href="https://2014.igem.org/Team:Bielefeld-CeBiTec/Results/Modelling/sourcecode" target="_blank">source code</a>). Two examples are shwon in Figure 3. The predicted changes of the metabolic concentration over the time is shown in Figures 4-7. The amount of expressed proteins could differ depending on the distance of the coding sequence downstream of the promoter. The coding sequences for the involved enzymes are located downstream of a common promoter. Therefore we decided to set the enzyme concentration for the first enzyme to 1 and decrease in steps of 0.1 (<a href="https://2014.igem.org/Team:Bielefeld-CeBiTec/Results/Modelling/sourcecode" target="_blank">source code</a>). These different values were tried to identify appropriate concentrations for each enzyme. The results of this dynamic modeling approach could be transferred to the laboratory by using promoters of different strength.
</p>
</p>
</div>
</div>
-
<br><br>
+
<br>
<center>
<center>
<div class="element" style="width:750px">  
<div class="element" style="width:750px">  
       <a href="https://static.igem.org/mediawiki/2014/7/7e/Bielefeld-CeBiTec_14_10_16_dgls.png" target="_blank"><img src="https://static.igem.org/mediawiki/2014/7/7e/Bielefeld-CeBiTec_14_10_16_dgls.png" width="750px"></a><br>
       <a href="https://static.igem.org/mediawiki/2014/7/7e/Bielefeld-CeBiTec_14_10_16_dgls.png" target="_blank"><img src="https://static.igem.org/mediawiki/2014/7/7e/Bielefeld-CeBiTec_14_10_16_dgls.png" width="750px"></a><br>
       <font size="2" style="text-align:left;">
       <font size="2" style="text-align:left;">
-
<b>Figure 3</b>: Differential equations for the dynamic modelling of the isobutanol production. </font>
+
<b>Figure 3</b>: Examples of differential equations for the dynamic modeling of the isobutanol production. </font>
</div>
</div>
</center>
</center>
-
<br><br>
+
<br>
<center>
<center>
<div class="element" style="width:750px">  
<div class="element" style="width:750px">  
-
       <a href="https://static.igem.org/mediawiki/2014/6/6c/Bielefeld-CeBiTec_14-10-16_without_improvement.jpg"><img src="https://static.igem.org/mediawiki/2014/6/6c/Bielefeld-CeBiTec_14-10-16_without_improvement.jpg" width="750px"></a><br>
+
       <a href="https://static.igem.org/mediawiki/2014/6/6c/Bielefeld-CeBiTec_14-10-16_without_improvement.jpg" target="_blank"><img src="https://static.igem.org/mediawiki/2014/6/6c/Bielefeld-CeBiTec_14-10-16_without_improvement.jpg" width="750px"></a><br>
       <font size="2" style="text-align:left;">
       <font size="2" style="text-align:left;">
<b>Figure 4</b>: Predicted changes in metabolic concentration over time. </font>
<b>Figure 4</b>: Predicted changes in metabolic concentration over time. </font>
</div>
</div>
</center>
</center>
 +
<br>
 +
 +
<br>
<br>
<center>
<center>
<div class="element" style="width:750px">  
<div class="element" style="width:750px">  
-
       <a href="https://static.igem.org/mediawiki/2014/1/1e/Bielefeld-CeBiTec_14-10-16_without_improvement_zoom_in.png"> <img src="https://static.igem.org/mediawiki/2014/1/1e/Bielefeld-CeBiTec_14-10-16_without_improvement_zoom_in.png" width="750px"></a><br>
+
       <a href="https://static.igem.org/mediawiki/2014/1/1e/Bielefeld-CeBiTec_14-10-16_without_improvement_zoom_in.png" target="_blank"> <img src="https://static.igem.org/mediawiki/2014/1/1e/Bielefeld-CeBiTec_14-10-16_without_improvement_zoom_in.png" width="750px"></a><br>
       <font size="2" style="text-align:left;">
       <font size="2" style="text-align:left;">
-
<b>Figure 5</b>: Predicted changes in metabolic concentration over time. The development in the first minutes is shown by zooming into the figure 4. </font>
+
<b>Figure 5</b>: Predicted changes in metabolic concentration over time. The development in the first minutes is shown by zooming into the Figure 4. </font>
</div>
</div>
</center>
</center>
Line 257: Line 249:
<center>
<center>
<div class="element" style="width:750px">  
<div class="element" style="width:750px">  
-
       <a href="https://static.igem.org/mediawiki/2014/8/83/Bielefeld-CeBiTec_14-10-16_improved_enzyme_concentrations.jpg"><img src="https://static.igem.org/mediawiki/2014/8/83/Bielefeld-CeBiTec_14-10-16_improved_enzyme_concentrations.jpg" width="750px"></a><br>
+
       <a href="https://static.igem.org/mediawiki/2014/8/83/Bielefeld-CeBiTec_14-10-16_improved_enzyme_concentrations.jpg" target="_blank"><img src="https://static.igem.org/mediawiki/2014/8/83/Bielefeld-CeBiTec_14-10-16_improved_enzyme_concentrations.jpg" width="750px"></a><br>
       <font size="2" style="text-align:left;">
       <font size="2" style="text-align:left;">
-
<b>Figure 4</b>: Predicted changes in metabolic concentration in an improved system over time. The concentration of IlvD and KivD was increased by factor 3 and 4 respectively. </font>
+
<b>Figure 6</b>: Predicted changes in metabolic concentration in an improved system over time. The concentration of KivD and AdhA was increased by factors 3 and 4, respectively. </font>
</div>
</div>
</center>
</center>
Line 269: Line 261:
" width="750px"></a><br>
" width="750px"></a><br>
       <font size="2" style="text-align:left;">
       <font size="2" style="text-align:left;">
-
<b>Figure 7</b>: Predicted changes in metabolic concentration in an improved system over time. The concentration of IlvD and KivD was increased by factor 3 and 4 respectively. </font>
+
<b>Figure 7</b>: Predicted changes in metabolic concentration in an improved system over time. The concentration of KivD and AdhA was increased by factor 3 and 4 respectively. The development in the first minutes is shown by zooming into the Figure 6. </font>
</div>
</div>
</center>
</center>
Line 275: Line 267:
-
<br><br>
+
<br>
<div id="text">
<div id="text">
<p>
<p>
-
Our modelling results indicated that the concentration of two enzymes are limiting the isobutanol production. They are <a href="https://2014.igem.org/Team:Bielefeld-CeBiTec/Project/Isobutanol/GeneticalApproach">IlvD</a> and <a href="https://2014.igem.org/Team:Bielefeld-CeBiTec/Project/Isobutanol/GeneticalApproach">KivD</a>. An experimentell verification of this effect is the next logical step. This enzymatic bottleneck could be removed by overexpression of the corresponding coding sequences (<a href="https://2014.igem.org/Team:Bielefeld-CeBiTec/Project/Isobutanol/GeneticalApproach"><i>ilvD</i></a> and <a href="https://2014.igem.org/Team:Bielefeld-CeBiTec/Project/Isobutanol/GeneticalApproach"><i>kivD</i></a>). One way to achieve this is the integration of a strong promotor upstream of theire coding sequences. It could be a possibility to improve our isobutanol production.
+
Our modeling results indicated that the concentrations of two enzymes are limiting the isobutanol production <a href="https://2014.igem.org/Team:Bielefeld-CeBiTec/Project/Isobutanol/GeneticalApproach" target="_blank">KivD</a> and <a href="https://2014.igem.org/Team:Bielefeld-CeBiTec/Project/Isobutanol/GeneticalApproach" target="_blank">AdhA</a>. An experimental validation of this effect is the next logical step. This enzymatic bottleneck could be removed by overexpression of the corresponding coding sequences (<a href="https://2014.igem.org/Team:Bielefeld-CeBiTec/Project/Isobutanol/GeneticalApproach" target="_blank"><i>kivD</i></a> and <a href="https://2014.igem.org/Team:Bielefeld-CeBiTec/Project/Isobutanol/GeneticalApproach" target="_blank"><i>adhA</i></a>). One way to achieve this is the integration of a strong promotor upstream of their coding sequences. It could be a possibility to improve our isobutanol production.
-
As shown in figure 4 the isobutanol concentration reaches 2mM after four hours. In the improved system the concentration of isobutanol after four hours is nearly doubled.
+
As shown in Figure 4 the isobutanol concentration reaches 2&nbsp;mM after four hours. In the improved system the concentration of isobutanol after four hours is nearly doubled.
</p>
</p>
</div>
</div>
</div>
</div>
</div>
</div>
 +
<div class="element" style="margin:10px; padding:10px">
<div class="element" style="margin:10px; padding:10px">
<div id="text">
<div id="text">
-
<br><br>
 
-
<h2>Conclusion</h2>
+
<h6>Conclusion</h6>
-
The complete metabolic network of our project is visualized in figure 1. Stoichiometric caculations revealed the number of electrons which is in theory need for the production of a desired substance. There are 42 electrons required for the synthesis of one molecule isobutanol. The fixation of one carbon dioxide molecule would cost seven electrons. We were able to identify a putative enzymatic bottleneck in the isobutanol production pathway. An increased amount of IlvD and KivD could improve the isobutanol production. Additionally the production of isobutanol from pyruvate can be predicted. The next and very important step would be the validation of these prediction by experimental approaches.
+
The complete metabolic network of our project is visualized in Figure 1. Stoichiometric calculations revealed the number of electrons which is in theory need for the production of a desired substance. There are 42 electrons required for the synthesis of one molecule isobutanol. The fixation of one carbon dioxide molecule would cost seven electrons. We were able to identify a putative enzymatic bottleneck in the isobutanol production pathway. An increased amount of KivD and AdhA could improve the isobutanol production. Additionally the production of isobutanol from pyruvate can be predicted. The next and very important step would be the validation of these predictions by experimental approaches.
</div>
</div>
Line 305: Line 297:
<div class="element" style="margin:10px; padding:10px">
<div class="element" style="margin:10px; padding:10px">
<div id="text">
<div id="text">
-
<br><br>
 
-
<h2>Outlook: addition of carbon dioxide fixing reactions</h2>
+
<h6>Outlook: addition of carbon dioxide fixing reactions</h6>
<p>
<p>
-
The next step would be the addition of specific carbon fixing reactions and the pathway leading to pyruvate to the existing dynamic model. Therefore we collected k<sub>cat</sub> and K<sub>M</sub> values for nearly all relevant steps (table 4). They were could be used to formulate additional differential equations, which describe these additional reactions.  
+
The next step would be the addition of specific carbon fixing reactions and the pathway leading to pyruvate in the existing dynamic model. Therefore we collected k<sub>cat</sub> and K<sub>M</sub> values for nearly all relevant steps (Table 4). They can be used to formulate additional differential equations which describe these additional reactions.  
</p>
</p>
-
<br><br>
+
<div class="element">
-
<b>Table4:</b> This table shows all k<sub>cat</sub> and K<sub>M</sub> values of enzymes involved in CO<sub>2</sub>-fixation and the pathway leading to pyruvate.
+
<font size="2"><b>Table 4:</b> This table shows all k<sub>cat</sub> and K<sub>M</sub> values of enzymes involved in CO<sub>2</sub> fixation and the pathway leading to pyruvate.</font>
<table width="100%" border="1" cellpadding="5" style="background-color:transparent">
<table width="100%" border="1" cellpadding="5" style="background-color:transparent">
<tr>
<tr>
<th>Enzyme</th>
<th>Enzyme</th>
-
<th>k<sub>cat</sub></th>
+
<th>k<sub>cat</sub> [-]</th>
<th>K<sub>M</sub> [mM]</th>
<th>K<sub>M</sub> [mM]</th>
<th>Reference</th>
<th>Reference</th>
Line 353: Line 344:
</tr>
</tr>
</table>
</table>
 +
</div>
<br><br>
<br><br>
-
Using all these information the combination of two parts of our project would be possible. The carbon dioxide fixation pathway could be checked for enzymatic bottlenecks. It would be possible to predict the isobutanol production from fixed carbon dioxide.
+
Using these information the combination of two parts of our project would be possible. The carbon dioxide fixation pathway could be checked for enzymatic bottlenecks. It would be possible to predict the isobutanol production from fixed carbon dioxide.
</div>
</div>
Line 369: Line 361:
<div id="text">
<div id="text">
-
<br><br>
+
  <h4 class="member" style="margin-left:20px">References</h4>
-
 
+
-
<div class="element">
+
-
  <div id="text">
+
-
  <h6>References</h6>
+
<ul>
<ul>
-
   
+
 
 +
  <div element="class">
<li id="atsumi2008">
<li id="atsumi2008">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Atsumi, Shota, Taizo Hanai, und James C. Liao. „Non-Fermentative Pathways for Synthesis of Branched-Chain Higher Alcohols as Biofuels“. <a href="http://www.nature.com/nature/journal/v451/n7174/full/nature06450.html">Nature</a> 451, no. 7174 (2008): 86–89.
+
Atsumi, Shota, Taizo Hanai, und James C. Liao. „Non-Fermentative Pathways for Synthesis of Branched-Chain Higher Alcohols as Biofuels“. <a href="http://www.nature.com/nature/journal/v451/n7174/full/nature06450.html" target="_blank">Nature</a> 451, no. 7174 (2008): 86–89.
</div>
</div>
</div>
</div>
Line 387: Line 376:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Atsumi, Shota, Tung-Yun Wu, Eva-Maria Eckl, Sarah D. Hawkins, Thomas Buelter, und James C. Liao. „Engineering the isobutanol biosynthetic pathway in Escherichia coli by comparison of three aldehyde reductase/alcohol dehydrogenase genes“. <a href="http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2802489/">Applied Microbiology and Biotechnology</a> 85, no. 3 (2010): 651–57.
+
Atsumi, Shota, Tung-Yun Wu, Eva-Maria Eckl, Sarah D. Hawkins, Thomas Buelter, und James C. Liao. „Engineering the isobutanol biosynthetic pathway in Escherichia coli by comparison of three aldehyde reductase/alcohol dehydrogenase genes“. <a href="http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2802489" target="_blank">Applied Microbiology and Biotechnology</a> 85, no. 3 (2010): 651–57.
</div>
</div>
</div>
</div>
Line 397: Line 386:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Yang, Y. T., G. N. Bennett, und K. Y. San. „The Effects of Feed and Intracellular Pyruvate Levels on the Redistribution of Metabolic Fluxes in Escherichia Coli“. <a href="http://bioc.rice.edu/bennetlab/references/ref40-sdarticle.pdf">Metabolic Engineering</a> 3, no. 2 (2001): 115–23.
+
Yang, Y. T., G. N. Bennett, und K. Y. San. „The Effects of Feed and Intracellular Pyruvate Levels on the Redistribution of Metabolic Fluxes in Escherichia Coli“. <a href="http://bioc.rice.edu/bennetlab/references/ref40-sdarticle.pdf" target="_blank">Metabolic Engineering</a> 3, no. 2 (2001): 115–23.
</div>
</div>
</div>
</div>
Line 406: Line 395:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Wadano, Akira, Keisuke Nishikawa, Tomohiro Hirahashi, Ryohei Satoh, und Toshio Iwaki. „Reaction Mechanism of Phosphoribulokinase from a Cyanobacterium, Synechococcus PCC7942“. <a href="http://link.springer.com/article/10.1023%2FA%3A1005979801741">Photosynthesis Research</a> 56, no. 1 (1998): 27–33.
+
Wadano, Akira, Keisuke Nishikawa, Tomohiro Hirahashi, Ryohei Satoh, und Toshio Iwaki. „Reaction Mechanism of Phosphoribulokinase from a Cyanobacterium, Synechococcus PCC7942“. <a href="http://link.springer.com/article/10.1023%2FA%3A1005979801741" target="_blank">Photosynthesis Research</a> 56, no. 1 (1998): 27–33.
</div>
</div>
</div>
</div>
Line 415: Line 404:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Kobayashi, Daisuke, Masahiro Tamoi, Toshio Iwaki, Shigeru Shigeoka, und Akira Wadano. „Molecular Characterization and Redox Regulation of Phosphoribulokinase from the Cyanobacterium Synechococcus Sp. PCC 7942“. <a href="http://www.ncbi.nlm.nih.gov/pubmed/12668773/">Plant & Cell Physiology</a> 44, no. 3 (2003): 269–76.
+
Kobayashi, Daisuke, Masahiro Tamoi, Toshio Iwaki, Shigeru Shigeoka, und Akira Wadano. „Molecular Characterization and Redox Regulation of Phosphoribulokinase from the Cyanobacterium Synechococcus Sp. PCC 7942“. <a href="http://www.ncbi.nlm.nih.gov/pubmed/12668773" target="_blank">Plant & Cell Physiology</a> 44, no. 3 (2003): 269–76.
</div>
</div>
</div>
</div>
Line 424: Line 413:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Lan, Yun, und Keith A. Mott. „Determination of Apparent Km Values for Ribulose 1,5-Bisphosphate Carboxylase/Oxygenase (Rubisco) Activase Using the Spectrophotometric Assay of Rubisco Activity“. <a href="http://www.plantphysiol.org/content/95/2/604.abstract">Plant Physiology</a> 95, Nr. 2 (1991): 604–9.
+
Lan, Yun, und Keith A. Mott. „Determination of Apparent Km Values for Ribulose 1,5-Bisphosphate Carboxylase/Oxygenase (Rubisco) Activase Using the Spectrophotometric Assay of Rubisco Activity“. <a href="http://www.plantphysiol.org/content/95/2/604.abstract" target="_blank">Plant Physiology</a> 95, Nr. 2 (1991): 604–9.
</div>
</div>
</div>
</div>
Line 434: Line 423:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Sage, Rowan F. „Variation in the K(cat) of Rubisco in C(3) and C(4) Plants and Some Implications for Photosynthetic Performance at High and Low Temperature“. <a href="http://www.ncbi.nlm.nih.gov/pubmed/11886880">Journal of Experimental Botany</a> 53, no. 369 (2002): 609–20.
+
Sage, Rowan F. „Variation in the K(cat) of Rubisco in C(3) and C(4) Plants and Some Implications for Photosynthetic Performance at High and Low Temperature“. <a href="http://www.ncbi.nlm.nih.gov/pubmed/11886880" target="_blank">Journal of Experimental Botany</a> 53, no. 369 (2002): 609–20.
</div>
</div>
</div>
</div>
Line 444: Line 433:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Fifis, T., und R. K. Scopes. „Purification of 3-Phosphoglycerate Kinase from Diverse Sources by Affinity Elution Chromatography“. <a href="http://www.ncbi.nlm.nih.gov/pubmed/?term=fifis+1978">The Biochemical Journal</a> 175, no. 1 (1978): 311–19.
+
Fifis, T., und R. K. Scopes. „Purification of 3-Phosphoglycerate Kinase from Diverse Sources by Affinity Elution Chromatography“. <a href="http://www.ncbi.nlm.nih.gov/pubmed/?term=fifis+1978" target="_blank">The Biochemical Journal</a> 175, no. 1 (1978): 311–19.
</div>
</div>
</div>
</div>
Line 454: Line 443:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Zhao G, Halbur T, Pankratz DC. Colonization of oropharynx and nasal cavity of CDCD pigs by a nontoxigenic strain of Pasteurella multocida type D. <a href="http://www.aasv.org/shap/abstracts/abstract.php?v3n3p113">J Swine Health Prod</a> 1995;3(3):113-115.
+
Zhao G, Halbur T, Pankratz DC. Colonization of oropharynx and nasal cavity of CDCD pigs by a nontoxigenic strain of Pasteurella multocida type D. <a href="http://www.aasv.org/shap/abstracts/abstract.php?v3n3p113" target="_blank">J Swine Health Prod</a> 1995;3(3):113-115.
</div>
</div>
</div>
</div>
Line 463: Line 452:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Fraser, H. I., M. Kvaratskhelia, und M. F. White. „The Two Analogous Phosphoglycerate Mutases of Escherichia Coli“. <a href="http://www.ncbi.nlm.nih.gov/pubmed/10437801">FEBS Letters</a> 455, no. 3 (1999): 344–48.
+
Fraser, H. I., M. Kvaratskhelia, und M. F. White. „The Two Analogous Phosphoglycerate Mutases of Escherichia Coli“. <a href="http://www.ncbi.nlm.nih.gov/pubmed/10437801" target="_blank">FEBS Letters</a> 455, no. 3 (1999): 344–48.
</div>
</div>
</div>
</div>
Line 472: Line 461:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
White, Malcolm F., und Linda A. Fothergill-Gilmore. „Development of a Mutagenesis, Expression and Purification System for Yeast Phosphoglycerate Mutase“. <a href="http://onlinelibrary.wiley.com/doi/10.1111/j.1432-1033.1992.tb17099.x/pdf">European Journal of Biochemistry</a> 207, no. 2 (1992): 709–14.
+
White, Malcolm F., und Linda A. Fothergill-Gilmore. „Development of a Mutagenesis, Expression and Purification System for Yeast Phosphoglycerate Mutase“. <a href="http://onlinelibrary.wiley.com/doi/10.1111/j.1432-1033.1992.tb17099.x/pdf" target="_blank">European Journal of Biochemistry</a> 207, no. 2 (1992): 709–14.
</div>
</div>
</div>
</div>
Line 481: Line 470:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Spring, Thomas G., und Finn Wold. „The Purification and Characterization of Escherichia Coli Enolase“. <a href="http://www.jbc.org/content/246/22/6797.long">Journal of Biological Chemistry</a> 246, no. 22 (1971): 6797–6802.
+
Spring, Thomas G., und Finn Wold. „The Purification and Characterization of Escherichia Coli Enolase“. <a href="http://www.jbc.org/content/246/22/6797.long" target="_blank">Journal of Biological Chemistry</a> 246, no. 22 (1971): 6797–6802.
</div>
</div>
</div>
</div>
Line 490: Line 479:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Albe, Kathy R., Margaret H. Butler, und Barbara E. Wright. „Cellular concentrations of enzymes and their substrates“. <a href="http://www.sciencedirect.com/science/article/pii/S0022519305802668">Journal of Theoretical Biology</a> 143, no. 2 (1990): 163–95.
+
Albe, Kathy R., Margaret H. Butler, und Barbara E. Wright. „Cellular concentrations of enzymes and their substrates“. <a href="http://www.sciencedirect.com/science/article/pii/S0022519305802668" target="_blank">Journal of Theoretical Biology</a> 143, no. 2 (1990): 163–95.
</div>
</div>
</div>
</div>
Line 499: Line 488:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Oria-Hernández, Jesús, Nallely Cabrera, Ruy Pérez-Montfort, und Leticia Ramírez-Silva. „Pyruvate Kinase Revisited: The Activating Effect of K+“. <a href="http://www.ncbi.nlm.nih.gov/pubmed/16147999">The Journal of Biological Chemistry</a> 280, no. 45 (2005): 37924–29.
+
Oria-Hernández, Jesús, Nallely Cabrera, Ruy Pérez-Montfort, und Leticia Ramírez-Silva. „Pyruvate Kinase Revisited: The Activating Effect of K+“. <a href="http://www.ncbi.nlm.nih.gov/pubmed/16147999" target="_blank">The Journal of Biological Chemistry</a> 280, no. 45 (2005): 37924–29.
</div>
</div>
</div>
</div>
Line 508: Line 497:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Breitling, Rainer, David Gilbert, Monika Heiner, und Richard Orton. „A Structured Approach for the Engineering of Biochemical Network Models, Illustrated for Signalling Pathways“. <a href="http://bib.oxfordjournals.org/content/9/5/404.short">Briefings in Bioinformatics</a> 9, no. 5 (2008): 404–21.
+
Breitling, Rainer, David Gilbert, Monika Heiner, und Richard Orton. „A Structured Approach for the Engineering of Biochemical Network Models, Illustrated for Signalling Pathways“. <a href="http://bib.oxfordjournals.org/content/9/5/404.short" target="_blank">Briefings in Bioinformatics</a> 9, no. 5 (2008): 404–21.
</div>
</div>
</div>
</div>
Line 517: Line 506:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Chubukov, Victor, Luca Gerosa, Karl Kochanowski, und Uwe Sauer. „Coordination of Microbial Metabolism“. <a href="http://www.nature.com/nrmicro/journal/v12/n5/abs/nrmicro3238.html">Nature Reviews Microbiology</a> 12, no. 5 (Mai 2014): 327–40.
+
Chubukov, Victor, Luca Gerosa, Karl Kochanowski, und Uwe Sauer. „Coordination of Microbial Metabolism“. <a href="http://www.nature.com/nrmicro/journal/v12/n5/abs/nrmicro3238.html" target="_blank">Nature Reviews Microbiology</a> 12, no. 5 (Mai 2014): 327–40.
</div>
</div>
</div>
</div>
Line 526: Line 515:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Schaber, J., W. Liebermeister, und E. Klipp. „Nested Uncertainties in Biochemical Models“. <a href="http://www.ncbi.nlm.nih.gov/pubmed/19154080">IET Systems Biology</a> 3, no. 1 (2009): 1–9.
+
Schaber, J., W. Liebermeister, und E. Klipp. „Nested Uncertainties in Biochemical Models“. <a href="http://www.ncbi.nlm.nih.gov/pubmed/19154080" target="_blank">IET Systems Biology</a> 3, no. 1 (2009): 1–9.
</div>
</div>
</div>
</div>
Line 537: Line 526:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Chuang, Han-Yu, Matan Hofree, und Trey Ideker. „A Decade of Systems Biology“. <a href="http://www.ncbi.nlm.nih.gov/pubmed/20604711">Annual Review of Cell and Developmental Biology</a> 26 (2010): 721–44.
+
Chuang, Han-Yu, Matan Hofree, und Trey Ideker. „A Decade of Systems Biology“. <a href="http://www.ncbi.nlm.nih.gov/pubmed/20604711" target="_blank">Annual Review of Cell and Developmental Biology</a> 26 (2010): 721–44.
</div>
</div>
</div>
</div>
Line 546: Line 535:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Kherlopian, Armen R., Ting Song, Qi Duan, Mathew A. Neimark, Ming J. Po, John K. Gohagan, und Andrew F. Laine. „A Review of Imaging Techniques for Systems Biology“. <a href="http://www.biomedcentral.com/1752-0509/2/74">BMC Systems Biology</a> 2, no. 1 (12. August 2008): 74.
+
Kherlopian, Armen R., Ting Song, Qi Duan, Mathew A. Neimark, Ming J. Po, John K. Gohagan, und Andrew F. Laine. „A Review of Imaging Techniques for Systems Biology“. <a href="http://www.biomedcentral.com/1752-0509/2/74" target="_blank">BMC Systems Biology</a> 2, no. 1 (12. August 2008): 74.
</div>
</div>
</div>
</div>
Line 555: Line 544:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Tyagi, R., Lee, Y.-T., Guddat, L.W., and Duggleby, R.G. (2005). Probing the mechanism of the bifunctional enzyme ketol-acid reductoisomerase by site-directed mutagenesis of the active site. <a href="http://onlinelibrary.wiley.com/doi/10.1111/j.1742-4658.2004.04506.x/pdf">FEBS J.</a> 272, 593–602.
+
Tyagi, R., Lee, Y.-T., Guddat, L.W., and Duggleby, R.G. (2005). Probing the mechanism of the bifunctional enzyme ketol-acid reductoisomerase by site-directed mutagenesis of the active site. <a href="http://onlinelibrary.wiley.com/doi/10.1111/j.1742-4658.2004.04506.x/pdf" target="_blank">FEBS J.</a> 272, 593–602.
</div>
</div>
</div>
</div>
Line 564: Line 553:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Atsumi, S., Li, Z., and Liao, J.C. (2009). Acetolactate Synthase from Bacillus subtilis Serves as a 2-Ketoisovalerate Decarboxylase for Isobutanol Biosynthesis in Escherichia coli. <a href="http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2753059/">Appl Environ Microbiol</a> 75, 6306–6311.
+
Atsumi, S., Li, Z., and Liao, J.C. (2009). Acetolactate Synthase from Bacillus subtilis Serves as a 2-Ketoisovalerate Decarboxylase for Isobutanol Biosynthesis in Escherichia coli. <a href="http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2753059" target="_blank">Appl Environ Microbiol</a> 75, 6306–6311.
</div>
</div>
</div>
</div>
Line 572: Line 561:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Flint, D.H., Emptage, M.H., Finnegan, M.G., Fu, W., and Johnson, M.K. (1993). The role and properties of the iron-sulfur cluster in Escherichia coli dihydroxy-acid dehydratase. <a href="http://www.ncbi.nlm.nih.gov/pubmed/8325851">J. Biol. Chem.</a> 268, 14732–14742.
+
Flint, D.H., Emptage, M.H., Finnegan, M.G., Fu, W., and Johnson, M.K. (1993). The role and properties of the iron-sulfur cluster in Escherichia coli dihydroxy-acid dehydratase. <a href="http://www.ncbi.nlm.nih.gov/pubmed/8325851" target="_blank">J. Biol. Chem.</a> 268, 14732–14742.
</div>
</div>
</div>
</div>
Line 581: Line 570:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Gocke, D., Nguyen, C.L., Pohl, M., Stillger, T., Walter, L., and Müller, M. (2007). Branched-Chain Keto Acid Decarboxylase from Lactococcus lactis (KdcA), a Valuable Thiamine Diphosphate-Dependent Enzyme for Asymmetric CC Bond Formation. <a href="http://onlinelibrary.wiley.com/doi/10.1002/adsc.200700057/pdf">Adv. Synth. Catal.</a> 349, 1425–1435.
+
Gocke, D., Nguyen, C.L., Pohl, M., Stillger, T., Walter, L., and Müller, M. (2007). Branched-Chain Keto Acid Decarboxylase from Lactococcus lactis (KdcA), a Valuable Thiamine Diphosphate-Dependent Enzyme for Asymmetric CC Bond Formation. <a href="http://onlinelibrary.wiley.com/doi/10.1002/adsc.200700057/pdf" target="_blank">Adv. Synth. Catal.</a> 349, 1425–1435.
</div>
</div>
</div>
</div>
Line 590: Line 579:
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
<div class="element" style="margin_10px 10px 10px 10px; padding:10px 10px 10px 10px">
  <div id="text">
  <div id="text">
-
Werther, T., Spinka, M., Tittmann, K., Schütz, A., Golbik, R., Mrestani-Klaus, C., Hübner, G., and König, S. (2008). Amino Acids Allosterically Regulate the Thiamine Diphosphate-dependent α-Keto Acid Decarboxylase from Mycobacterium tuberculosis. <a href="http://www.jbc.org/content/283/9/5344.long">J. Biol. Chem.</a> 283, 5344–5354.
+
Werther, T., Spinka, M., Tittmann, K., Schütz, A., Golbik, R., Mrestani-Klaus, C., Hübner, G., and König, S. (2008). Amino Acids Allosterically Regulate the Thiamine Diphosphate-dependent α-Keto Acid Decarboxylase from Mycobacterium tuberculosis. <a href="http://www.jbc.org/content/283/9/5344.long" target="_blank">J. Biol. Chem.</a> 283, 5344–5354.
</div>
</div>
</div>
</div>
Line 598: Line 587:
-
 
+
</div>
-
 
+
</ul>
</ul>
-
</div>
 
-
</div>
 
</div>
</div>
</html>
</html>

Latest revision as of 03:04, 18 October 2014


Modeling

Summary

Different modeling approaches were used to identify enzymatic bottlenecks in the isobutanol production pathway and to predict the formation of the desired product in a given time. First of all we created a map of all metabolic reactions which are part of our project (Figure 1). This network not only provides a good overview, it also serves as the basic tool for further considerations. Due to the huge amount of components it does not seem feasible to create a computational model for all reactions at once. Therefore we started our modeling work by carrying out a stochiometric analysis. It revealed that 42 electrons are needed for the production of one isobutanol molecule from CO2. As for the isobutanol production pathway, dynamic modeling was carried out, in the course of which several bottlenecks could be identified. An increase in expression of kivD and adhA could boost the isobutanol production after four hours by 100 %. Finally we prepared the extension of the existing model to predict the effect of specific carbon dioxide fixing reactions.

Introduction

Mathematical modeling is crucial to understand complex biological systems (Schaber et al., 2009). The analysis of isolated biological components has been supplemented by a systems biology approach in over ten years (Chuang et al., 2010). Furthermore mathematical modeling is used to combine biological results (Kherlopian et al., 2008). Modeling is also a way to achieve results without carrying out experiments in a laboratory. The behaviour of a system can be simulated to get results which cannot be derived from simply looking at the given system (Schaber et al., 2009). The most important aim of any modeling approach is the reduction of complexity. As the given biological reality is often diverse and variable, it is important to identify the major rules and principles which can describe a system.

Our aims
  • Visualization of the complete metabolic network
  • Relating electron input to product formation
  • Relating electron input to carbon dioxide fixation
  • The identification of bottlenecks in the isobutanol production pathway
  • Prediction of maximal isobutanol production over the time
Stoichiometric analysis

Stoichiometric analysis is useful to get information about the maximal output of a system. In this project the number of electrons moving into the system limits the amount of the product which could be synthetized in the cells. Therefore the stoichiometric relations of all substances involved in our complex reaction network were calculated (Figure 1). The calculation starts with the electrons, which are transported into the system by mediators. We calculated the resulting production of intracellular molecules based on our map of the metabolic network (Figure 1). The results are listed below.


Figure 1: Complete metabolic network of reactions involved in our project.

Table 1: This table shows the theoretical electron cost of different intracellular molecules. The electron cost is the number of electrons which are needed for the synthesis of this metabolite. The intermediates are substances which are needed for the production of the final product. Their costs are included in the final value for each substance.
Substance Intermediates Number of electrons
FADH2 - 2
NADH+H+ - 2
ATP - 1
Triosephosphate 9 ATP + 6 NADH 21
CO2 fixation - 7
Pyruvate Triosephosphate 19
Isobutanol 2x Pyruvate 42


The theoretical electron costs of different molecules is listed in Table 1. All calculations are based on our pathway map. According to our pathway map there are 42 electrons needed for the production of one molecule isobutanol, if CO2 is used as sole carbon source. Our calculation does not involve the house keeping metabolism of E. coli, which consumes lots of energy for the survival of the cell. The true number of consumed electrons per produced isobutanol molecule is therefore much higher. The applied electric power can be converted into a number of electrons by the following equation: 1 A = 1 C * s-1 = 6.2415065 * 1018 electrons.











Dynamic modeling

After the stoichiometric analysis of the system we designed a dynamic model containing kinetic equations for the metabolite concentrations. It allows the identification of metabolic or enzymatic bottlenecks, the major aim of our modeling approach. This information might then be used in the next step to modify constructs e.g. exchange an RBS or a promotor sequence. This could help to balance the different enzyme concentrations in upcoming experiments in the laboratory, it was our goal to predict the production of isobutanol per substrate in a given time. Our model predicts the isobutanol production in a carbon dioxide fixing cell. To achieve our aims we reduced the complex system shown in Figure 1 to the version shown in Figure 2. This metabolic network was suitable for our dynamic modeling approach.



Figure 2: Reduced metabolic network of reactions which were selected for modeling.

The modeling work on the isobutanol pathway is based on publications about the isobutanol production pathway (Atsumi et al., 2008 and Atsumi et al., 2010). We started our work on the isobutanol production pathway by collecting the appropriate kinetic parameters. They were used for the development of a system of differential equations. Concerning the choice of the kinetics used, we stuck to Michaelis-Menten kinetics as this was published as the best approach if reaction kinetics are not known (Breitling et al., 2008; Chubukov et al., 2014). Additionally kcat and KM values were collected from databases like KEGG, biocyc and BRENDA (Table 2). Missing values had to be estimated. The starting concentrations for different metabolites were also taken from the literature and from different databases (Table 3).
We decided to use only the forward reactions from pyruvate to isobutanol for different reasons. First of all, it is necessary to get the maximal product concentration, secondly the reactions are nearly irreversible due to the decarboxylation steps and thirdly there is only few information on the back reactions available.

Table 2: This table shows all enzymatic parameters which were used for our dynamic model.
Enzyme kcat [-] KM [mM] Reference
AlsS 121 13.6 Atsumi et al., 2009
IlvC 2.2 0.25 Tyagi et al., 2005
IlvD 10 (estimated) 1.5 Flint et al., 1993
KivD 20 (estimated) 5 (estimated) Werther et al., 2008; Gorcke et al., 2007
AdhA 6.6 9.1 Atsumi et al., 2010

Table 3: This table shows all metabolite concentrations which were used for our first model. The metabolite concentration was set to zero, if no published value was available.
Metabolite Concentration [mM] Reference
Pyruvate 10 Yang et al., 2000
2-Acetolactate 0 -
2,3-Dihydroxyisovalerate 0 -
2-Ketoisovalerate 0 -
Isobutyraldehyde 0.6 (estimated)
Isobutanol variable -

We implemented the system of differential equations in Matlab (source code). Two examples are shwon in Figure 3. The predicted changes of the metabolic concentration over the time is shown in Figures 4-7. The amount of expressed proteins could differ depending on the distance of the coding sequence downstream of the promoter. The coding sequences for the involved enzymes are located downstream of a common promoter. Therefore we decided to set the enzyme concentration for the first enzyme to 1 and decrease in steps of 0.1 (source code). These different values were tried to identify appropriate concentrations for each enzyme. The results of this dynamic modeling approach could be transferred to the laboratory by using promoters of different strength.



Figure 3: Examples of differential equations for the dynamic modeling of the isobutanol production.


Figure 4: Predicted changes in metabolic concentration over time.



Figure 5: Predicted changes in metabolic concentration over time. The development in the first minutes is shown by zooming into the Figure 4.


Figure 6: Predicted changes in metabolic concentration in an improved system over time. The concentration of KivD and AdhA was increased by factors 3 and 4, respectively.


Figure 7: Predicted changes in metabolic concentration in an improved system over time. The concentration of KivD and AdhA was increased by factor 3 and 4 respectively. The development in the first minutes is shown by zooming into the Figure 6.

Our modeling results indicated that the concentrations of two enzymes are limiting the isobutanol production KivD and AdhA. An experimental validation of this effect is the next logical step. This enzymatic bottleneck could be removed by overexpression of the corresponding coding sequences (kivD and adhA). One way to achieve this is the integration of a strong promotor upstream of their coding sequences. It could be a possibility to improve our isobutanol production. As shown in Figure 4 the isobutanol concentration reaches 2 mM after four hours. In the improved system the concentration of isobutanol after four hours is nearly doubled.

Conclusion
The complete metabolic network of our project is visualized in Figure 1. Stoichiometric calculations revealed the number of electrons which is in theory need for the production of a desired substance. There are 42 electrons required for the synthesis of one molecule isobutanol. The fixation of one carbon dioxide molecule would cost seven electrons. We were able to identify a putative enzymatic bottleneck in the isobutanol production pathway. An increased amount of KivD and AdhA could improve the isobutanol production. Additionally the production of isobutanol from pyruvate can be predicted. The next and very important step would be the validation of these predictions by experimental approaches.
Outlook: addition of carbon dioxide fixing reactions

The next step would be the addition of specific carbon fixing reactions and the pathway leading to pyruvate in the existing dynamic model. Therefore we collected kcat and KM values for nearly all relevant steps (Table 4). They can be used to formulate additional differential equations which describe these additional reactions.

Table 4: This table shows all kcat and KM values of enzymes involved in CO2 fixation and the pathway leading to pyruvate.
Enzyme kcat [-] KM [mM] Reference
PrkA 72.6 0.09 Wadano et al., 1998,Kobayashi et al., 2003
RubisCO 20 (estimated) 0.02 (estimated) Lan and Mott, 1991,Sage, 2002
Pgk 480 1 (estimated) Fifis and Scopes., 1978
GapA - 0.5 Zhao et al., 1995
GpmA 490 (in S.cerevisiae) 0.15 Fraser et al., 1999,White and Fothergil-Gilmore, 1992
Eno 17600 0.1 Spring and Wold, 1972, Albe et al., 1990
PykF 3.2 0.3 (estimated) Oria-Hernandez et al., 2005


Using these information the combination of two parts of our project would be possible. The carbon dioxide fixation pathway could be checked for enzymatic bottlenecks. It would be possible to predict the isobutanol production from fixed carbon dioxide.

References

  • Atsumi, Shota, Taizo Hanai, und James C. Liao. „Non-Fermentative Pathways for Synthesis of Branched-Chain Higher Alcohols as Biofuels“. Nature 451, no. 7174 (2008): 86–89.
  • Atsumi, Shota, Tung-Yun Wu, Eva-Maria Eckl, Sarah D. Hawkins, Thomas Buelter, und James C. Liao. „Engineering the isobutanol biosynthetic pathway in Escherichia coli by comparison of three aldehyde reductase/alcohol dehydrogenase genes“. Applied Microbiology and Biotechnology 85, no. 3 (2010): 651–57.
  • Yang, Y. T., G. N. Bennett, und K. Y. San. „The Effects of Feed and Intracellular Pyruvate Levels on the Redistribution of Metabolic Fluxes in Escherichia Coli“. Metabolic Engineering 3, no. 2 (2001): 115–23.
  • Wadano, Akira, Keisuke Nishikawa, Tomohiro Hirahashi, Ryohei Satoh, und Toshio Iwaki. „Reaction Mechanism of Phosphoribulokinase from a Cyanobacterium, Synechococcus PCC7942“. Photosynthesis Research 56, no. 1 (1998): 27–33.
  • Kobayashi, Daisuke, Masahiro Tamoi, Toshio Iwaki, Shigeru Shigeoka, und Akira Wadano. „Molecular Characterization and Redox Regulation of Phosphoribulokinase from the Cyanobacterium Synechococcus Sp. PCC 7942“. Plant & Cell Physiology 44, no. 3 (2003): 269–76.
  • Lan, Yun, und Keith A. Mott. „Determination of Apparent Km Values for Ribulose 1,5-Bisphosphate Carboxylase/Oxygenase (Rubisco) Activase Using the Spectrophotometric Assay of Rubisco Activity“. Plant Physiology 95, Nr. 2 (1991): 604–9.
  • Sage, Rowan F. „Variation in the K(cat) of Rubisco in C(3) and C(4) Plants and Some Implications for Photosynthetic Performance at High and Low Temperature“. Journal of Experimental Botany 53, no. 369 (2002): 609–20.
  • Fifis, T., und R. K. Scopes. „Purification of 3-Phosphoglycerate Kinase from Diverse Sources by Affinity Elution Chromatography“. The Biochemical Journal 175, no. 1 (1978): 311–19.
  • Zhao G, Halbur T, Pankratz DC. Colonization of oropharynx and nasal cavity of CDCD pigs by a nontoxigenic strain of Pasteurella multocida type D. J Swine Health Prod 1995;3(3):113-115.
  • Fraser, H. I., M. Kvaratskhelia, und M. F. White. „The Two Analogous Phosphoglycerate Mutases of Escherichia Coli“. FEBS Letters 455, no. 3 (1999): 344–48.
  • White, Malcolm F., und Linda A. Fothergill-Gilmore. „Development of a Mutagenesis, Expression and Purification System for Yeast Phosphoglycerate Mutase“. European Journal of Biochemistry 207, no. 2 (1992): 709–14.
  • Spring, Thomas G., und Finn Wold. „The Purification and Characterization of Escherichia Coli Enolase“. Journal of Biological Chemistry 246, no. 22 (1971): 6797–6802.
  • Albe, Kathy R., Margaret H. Butler, und Barbara E. Wright. „Cellular concentrations of enzymes and their substrates“. Journal of Theoretical Biology 143, no. 2 (1990): 163–95.
  • Oria-Hernández, Jesús, Nallely Cabrera, Ruy Pérez-Montfort, und Leticia Ramírez-Silva. „Pyruvate Kinase Revisited: The Activating Effect of K+“. The Journal of Biological Chemistry 280, no. 45 (2005): 37924–29.
  • Breitling, Rainer, David Gilbert, Monika Heiner, und Richard Orton. „A Structured Approach for the Engineering of Biochemical Network Models, Illustrated for Signalling Pathways“. Briefings in Bioinformatics 9, no. 5 (2008): 404–21.
  • Chubukov, Victor, Luca Gerosa, Karl Kochanowski, und Uwe Sauer. „Coordination of Microbial Metabolism“. Nature Reviews Microbiology 12, no. 5 (Mai 2014): 327–40.
  • Schaber, J., W. Liebermeister, und E. Klipp. „Nested Uncertainties in Biochemical Models“. IET Systems Biology 3, no. 1 (2009): 1–9.
  • Chuang, Han-Yu, Matan Hofree, und Trey Ideker. „A Decade of Systems Biology“. Annual Review of Cell and Developmental Biology 26 (2010): 721–44.
  • Kherlopian, Armen R., Ting Song, Qi Duan, Mathew A. Neimark, Ming J. Po, John K. Gohagan, und Andrew F. Laine. „A Review of Imaging Techniques for Systems Biology“. BMC Systems Biology 2, no. 1 (12. August 2008): 74.
  • Tyagi, R., Lee, Y.-T., Guddat, L.W., and Duggleby, R.G. (2005). Probing the mechanism of the bifunctional enzyme ketol-acid reductoisomerase by site-directed mutagenesis of the active site. FEBS J. 272, 593–602.
  • Atsumi, S., Li, Z., and Liao, J.C. (2009). Acetolactate Synthase from Bacillus subtilis Serves as a 2-Ketoisovalerate Decarboxylase for Isobutanol Biosynthesis in Escherichia coli. Appl Environ Microbiol 75, 6306–6311.
  • Flint, D.H., Emptage, M.H., Finnegan, M.G., Fu, W., and Johnson, M.K. (1993). The role and properties of the iron-sulfur cluster in Escherichia coli dihydroxy-acid dehydratase. J. Biol. Chem. 268, 14732–14742.
  • Gocke, D., Nguyen, C.L., Pohl, M., Stillger, T., Walter, L., and Müller, M. (2007). Branched-Chain Keto Acid Decarboxylase from Lactococcus lactis (KdcA), a Valuable Thiamine Diphosphate-Dependent Enzyme for Asymmetric CC Bond Formation. Adv. Synth. Catal. 349, 1425–1435.
  • Werther, T., Spinka, M., Tittmann, K., Schütz, A., Golbik, R., Mrestani-Klaus, C., Hübner, G., and König, S. (2008). Amino Acids Allosterically Regulate the Thiamine Diphosphate-dependent α-Keto Acid Decarboxylase from Mycobacterium tuberculosis. J. Biol. Chem. 283, 5344–5354.