Difference between revisions of "Team:Toulouse/Modeling"

As said before, the aim of our project is to create a biological system able to produce two molecules: butyric acid and formic acid.
To achieve this, we need to modify the existing balance between the metabolic pathways present in E. Coli. Indeed, we want to optimize butyrate and formate productions in our bacterium by adjusting environmental conditions in order to obtain the desired concentrations of the associated acids. The following metabolic network represents all of the known metabolites and metabolic pathways for E. coli K12 MG1655 (best known model) as of today. It was obtained from the KEGG database. Our first step was to identify the pathways in which our molecules take part, in order to have a clear understanding of their role and effect.

Formate is naturally produced by E.coli but at a level that is quite low. Our project requires that Apicoli produces more. Hence we had to optimize its biosynthesis by studying the genes coding for the enzymes involved in the pathway. We decided to focus our efforts on the Pyruvate Formate Lyase (PFL), the enzyme that causes degradation of pyruvate, thus yielding formate.

The subnetwork presented below was obtained from MetExplore platform[2] and presents all reactions from the KEGG and ByoCyc databases involved in the production or consumption of formate. This map will help us predict the likely consequences of a PFL-induced formate overproduction in Apicoli. In fact, formate is harmful to our bacterium and is normally metabolized to other products. We thus have to find the balance between producing enough formate to kill the varroa without killing Apicoli.

Contrary to formate, butyrate is not naturally produced by E. Coli. E.coli possesses an enzyme called (Butyryl-coA transferase (or Acetate-coA transferase) that yields butyrate, but this reaction cannot happen spontaneously in the organism due to the lack of Butanoyl-coA, its substrate. Indeed, study of the biosynthesis pathway shows that the enzymes responsible for Butanoyl-coA production (EC.2.1.3.19 phosphate butyryltransferase , EC.1.3.1.44 trans-2-enoyl-CoA reductase (NAD+) , etc.) cannot be found in our strain. Hence, in order to obtain butyrate, we chose to introduce a complete production pathway relying on genes coming from different organisms in Apicoli (see Attract).

To go further in the development of our project, we decided to use a method called Flux Balance Analysis (FBA) and Flux Variability Analysis (FVA). It is based on the model EC_iJO1366 [1]. This is the most recent model concerning E.coli K12 MG1655. It is a stoichiometric model defining all metabolic ways known to this day in this particular strain. It has been adapted to better suit our own bacterial system. Here you will find the associated XML file. The modifications we made are such that the described strain is now capable of producing butyrate.
Our modeling aims at determining the maximum butyrate and formate quantities our system would be able to produce, depending on two initial conditions: oxygen and glucose flux. We set them as follows:

It is interesting to note that FBA provides the produced and consumed metabolites in a quantitative way.
The setpoint for oxygen flux is rather low in order to simulate microaerobic conditions. The glucose flux setpoint was chosen according to the results of our tests with the Biosilta kit (see “Preliminary part”).
This kit was used to ensure a stable availability of glucose over time for our bacteria. Indeed the medium contains enzymes capable of catabolizing starch, thus gradually releasing glucose in the culture.
For known initial quantities of starch and enzymes, we are able to deduce the glucose release flux. Thus, we chose the appropriate enzyme concentration and polymer quantity in order to have a glucose rate of 0.3998 mmol.gDW^-1.h^-1.

The model provides results for the metabolites flux in mmol.gDW^-1.h^-1 but it is difficult to get an idea of the actual quantity this represents so we will convert it to mmol/L. To do this we chose a length of time of 13 hours for the day and 7 hours for the night, since our solution will primarily be deployed during summer. This means that butyric acid production time is estimated to 13 hours and formic acid production time to 7 hours.
Concerning biomass it is more complicated since the bacterium grows all the time. Thus in order to ensure a minimum production, biomass concentration (X) at the beginning is defined to X = 0.56 gDW.L^-1.

$$ \textrm{Real unit} (mmol.L^{-1})= \textrm{Model unit} (mmol.gDW^{-1}\cdot h^{-1})\times X \times time $$ $$ P_A = H_A\cdot C_{A,eq}, \textrm{Henry's law} $$ $$ C_{A,eq} = \frac{19,964}{0,019} = 1,019mol.L^{-1}$$