Modelling

Introduction

Mathematical modelling is essentiell to understand complex biological systems (Klipp et al., 2009). The analysis of isolated biological components is supplemented by a systems biology approach since ten years ago (Chuang et al., 2010). Mathematical modelling is used to combine biological results (Kherlopian et al., 2008). Modelling is also a way to gain 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 (Klipp et al., 2009). The most important aim of any modelling approach is the reduction of complexity. The given biological reality is often very divers and variable. Therefor it is important to identify the major rules and principles which can describe a system.

Results

First of all it was our aim to predict the production of isobutanol. Our model should give information about the optimal concentration of each enzyme in the isobutanol production pathway. The next aim was the prediction of isobutanol production in a carbon dioxide fixing cell. The complete system is shown in figure 1. This complex network of reactions was reduced to the system shown in fig. 2. This reduced version was used for modelling.

[komplette Grafik] Fig.1: Complete metabolic network of reactions which describes a part of our project. [reduzierte Grafik = modellierter Teil des Stoffwechsels] Fig.2: Reduced metabolic network of reactions which were selected for modelling.

We started our modelling work by reading publications about the isobutanol production pathway (Atsumi et al., 2008 and Atsumi et al., 2008). The first modelling approach was a system of differential equations using Michealis-Menten kinetics. This was published as the best approach if reaction kinetics are not known (REFERENZ EINFÜGEN). All needed Vmax and KM values were colleted from the literature and from databases like KEGG, biocyc and BRENDA(table 1).

Table1: This table shows all enzymatic parameters which were used for our first model.

The starting concentrations for different metabolites were also taken from the literature and from different databases (table 2).

Table2: 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.

Enzyme	Vmax	KM	Reference
AlsS
IlvC
IlvD
KivD
AdhA

We implemented the system of differential equations in matlab (link to source code) and created first results (fig. 3).

[Plot für Vmax und Km] Fig.3: Predicted changes in metabolic concentration over time.

To improve our prediction we decided to switch from Vmax and Km to kcat and the amount of the different enzymes (table 3). The amount of expressed proteins could differ depending on the distance of the coding sequence downstream of the promotor. Different values can be used to simulate the usage of promotors of different strength. This approach also allows the modelling of different growth states. The growth is represented by an increase in the amount of enzyme.

Table3: This table shows all kcat values which were used for modelling of the isobutanol production pathway. kcat was set to 10, if no published value was available.

Metabolite	Concentration	Reference
Pyruvate
2-Acetolactate
2,3-Dihydroxyisovalerate
2-Ketoisovalerate
Isobutyraldehyde
Isobutanol

The modelling results indicated that the concentration of IlvD is limiting the isobutanol production. This bottle neck could be removed by overexpression of the corresponding coding sequence. This could be achived by using a strong promotor and RBS in front of this coding sequence.

The next model improvement was the addition of some of the carbon fixing reactions and the pathway leading to pyruvate. We used kcat values for all relevant steps (fig.2 and table 4).

Table4: This table shows all kcat values of enzymes involved in CO2-fixation and the pathway leading to pyruvate.

Enzyme	kcat	Reference
AlsS
IlvC
IlvD
KivD
AdhA

Stoichiometric analysis

We calculated the stoichiometric relations of all substances involved in our complex reaction network (fig. x). Starting with the electrons which are transported into the system by mediators we calculated the resulting production of all other molecules. The results are shown in fig. XXX and listed below.

[Grafik mit Zahlen an Interaktionen in Pathway]

Fig.X: stoichiometric relations of all substances

In theory there are XXXX 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 its survival. The number of consumed electrons per produced isobutanol molecule is therefor 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 * 10^18 electrons.

References

Atsumi, Shota, Taizo Hanai, und James C. Liao. „Non-Fermentative Pathways for Synthesis of Branched-Chain Higher Alcohols as Biofuels“. Nature 451, Nr. 7174 (3. Januar 2008): 86–89. doi:10.1038/nature06450. http://www.nature.com/nature/journal/v451/n7174/full/nature06450.html
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, Nr. 3 (Januar 2010): 651–57. doi:10.1007/s00253-009-2085-6. http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2802489/

Enzyme	kcat	Reference
AlsS
IlvC
IlvD
KivD
AdhA

Team:Bielefeld-CeBiTec/Results/Modelling/erster/test/123

From 2014.igem.org