Team:Evry/Metabolism model
From 2013.igem.org
Line 175: | Line 175: | ||
</li>lowerbound and upperbound are vector of constraints (minimal and maximal flux values for each reactions)</li> | </li>lowerbound and upperbound are vector of constraints (minimal and maximal flux values for each reactions)</li> | ||
</ul> | </ul> | ||
+ | |||
+ | The values for the lower and upper bounds on the flux of each reactions are either deduced from experiments or put to a very high value when unknown (most of the time) : | ||
+ | |||
+ | <div align="center" style="margin:1%;"> | ||
+ | <img src="https://static.igem.org/mediawiki/2013/3/39/Fluxboundaries.png" alt="FBA boundaries" /> | ||
+ | </div> | ||
+ | |||
+ | Tunning these boundaries allows to represent different experimental conditions, for example reducing an upper bound to a low value may represent a loss of reaction flux due to the scarcity of a certain compound. | ||
</p> | </p> |
Revision as of 00:24, 5 October 2013
Introduction
In this part of our modeling work we focus on genome scale analysis of the enterobactin production pathway. A major concern about our system is its non-orthogonality with the natural metabolic network of the cells : E.Coli already possesses the genes for producing enterobactins. Therefore we wanted to assess the possible interactions between our system and the bacterial metabolism.
Observations
This model stems from the following observations :
- Enterobactin production pathway is a metabolic process;
- Any of the involved metabolites may limit the rate of the reactions.
As can be seen in the figure the pathway is 4 steps long and composed of six different enzymes. Hence there exists 4 possible metabolites which concentration may be limiting :
- Chorismate
- Isochorismate
- 2,3-dihydroxy-2,3-dihydrobenzoate
- 2,3-dihydroxybenzoat
Goals
We highlighted two main kind of interactions between the bacterial and our system :
- The synthetic system interacts with the bacterial metabolism. Leading to scarcity of the metabolites involved in the pathway for the other (possibly essential) metabolic reactions of the cell.
- The other way round, the metabolic reactions could prevent our synthetic system to work as expected by limiting the quantity of the involved metabolites available.
- Is the metabolic model of E.coli able to provide any information about the possible interactions between our system and the metabolism?
- Is the concentration of any metabolite limiting ?
- In the latter case, what is the quantity of this metabolite to add?
Materials and Methods
Model
We used the metabolic model E.Coli iJR904 downloaded from the BiGG model database (http://bigg.ucsd.edu/bigg/home.pl). We chose this model because it contains all the metabolites involved in the enterobactin production pathway.
This model contains 4037 reactions and 625 metabolites but lacks the enterobactin syntase (the last reaction of the enterobactin production pathway). Thus we extended it, adding this reactions and an enterobactin output reaction to be able to consider enterobactin as a component of the model. The inserted reactions are the following :
Name | Formula |
---|---|
ENTSYNTH | 6 ATP + 3 2,3-Dihydroxybenzoate + 3 L-Serine <=> Enterobactin + 6 AMP + 6 Diphosphate |
ENTOUT | => -1 Enterobactin |
Network Reduction
As an analysis of the whole network is too complex, we focused our analysis on the subnetwork presented in Figure 2. We considered the CHORISMIC ACID (hereafter called CHORISMATE) as the first step of the anterobactin production pathway. This is justified by the fact that chorismate can be bought from any compound supplier (for example sigma).
We also found that it would also be possible to by the 2-3-DIHYDROXYBENZOATE compound from the same provider. This compound may be very interesting to test our constructions later as the last precursor of the pathway. For the same reason, and because it does participate in any other reaction in the used model, this compound has no interest in our metabolic modeling.
Flux Balance Analysis
The metabolic model of E.Coli is represented as a stoichiometry matrix S representing the metabolic network. Our modified version has the size 4039 (reactions) * 625 (compounds). The unknown is the flux distribution vector v, a 4039-vector representing the flow of matter (mmol/gDW/h) going through each reactions.
The Flux Balance Analysis method is about finding this flux repartition vector v given an objective function to optimize (usually the growth rate) and a set of constraints of fluxes.
Assumptions of the model
The Two assumptions at the heart of the method are the following :
- steady state: The fluxes are considered to have attained a static equilibrium value and do not change through time.
- No enzyme saturation: The enzymes are supposed to be not saturated, the number of enzymes is always greater than the number of the corresponding reactions happening.
Formalism
The FBA method use a representation of the metabolic reaction network in the form of a stoichiometry matrix S where :
- Each row corresponds to a reaction R_i
- Each column corresponds to a metabolite C_j
- v is the vector of unknown reaction fluxes
- c is a vector of constants defining the objective function
- S is the stoichiometry matrix lowerbound and upperbound are vector of constraints (minimal and maximal flux values for each reactions)
References: