The bigg_e_coli_core model is a very simplified metabolic network model of Escherichia coli K12.
This toy model is almost identical to the one in
et al. (2010). Note that we can only recommend reading this excellent primer as a complement to our hands-on tutorial.
- Reset your session.
- Pick the bigg_e_coli_core E. coli model from the model and pathway repository.
- What are the total number of reactions, chemical compounds, and compartments in the model?
- This is automatically displayed in the sandbox: #reac=97, #chem=56 and #comp=3. Indeed a very small model.
- What is the total number of external (or boundary) reactions?
- What are the total numbers of reversible or irreversible reactions?
- How many reactions occur in a single compartment only? What kind of reactions are the remaining
reactions? What is their total number?
- The answers can be found in the BC (Basic Classification) summary
- Restrict the display to external reactions in the reac_list view
by typing BOUNDARY in the search field.
What can you tell about the growth conditions, i.e. about the "growth medium" of
E. coli K12 as modelled by bigg_e_coli_core? What does it mean that some
external reactions are unconstrained/unidirectional/bidirectional?
- D-glucose is the sole carbon source that can be consumed with an absolute uptake rate of 10.0
(the negative sign indicates that this boundary reaction proceeds from right to left).
- Phosphate, ammonium, water, proton, oxygen, and carbon dioxide are freely available.
They can be exchanged in both consumed and secreted as no flux constraints are placed on these reactions.
- The remaining organic compounds (e.g. acetate) occurring in the external compartment can only be exported and thus correspond to possibly secreted metabolites.
The bigg_e_coli_core model is a simplified metabolic network.
- What can you tell about its growth reaction? Is it biologically relevant?
- What can you tell about its sulfur metabolism?
- Identify the growth reaction in the reac_list view by, for example, typing BIOMASS in the search field... What a simplified growth reaction!! Only two amino acids are required for growth: glutamate and glutamine.
- There are no sulfur compounds in the set of boundary reactions - the sulfur metabolism is completely ignored although sulfur atoms occur in some compounds, e.g. CoA!
Run the GCR (groups of coupled reactions) and FBA (flux balance analysis) on bigg_e_coli_core model.
- Does the model grow?
- Can you identify any reaction correlation groups? How do they relate to the flux distribution?
- Can you identify any reaction correlation groups that cannot carry a flux in steady state? If so, why?
The idea is to create a new model for the anaerobic growth of E. coli K12:
- Determine the Reac ID of the boundary reaction of O2
- Export the bigg_e_coli_core model as an Excel spreadsheet (or in any other format, but Excel is rather convenient for performing a quick change)
- Find the row with the correct Reac ID in the enzymes sheet and set the flux boundaries of O2 to [0;0]. What does this mean?
- Save the model on the user side using a new filename.
- Upload the new model using a new and informative model identifier, e.g. anaerobic_textbook
Inspect the original and modified models using the reac_list view. Can you identify any differences between the original and the modified model? If so, which ones?
Run FBA on the anaerobic bigg_e_coli_core model.
- Can you identify any differences in the external flux distributions between the aerobic and anaerobic models?
- Can you establish the overall growth stoichiometries for both models (hint: examine at the flux distribution of the boundary reactions). Are they compatible with E. coli K12 biology?
- Can you identify any internal reaction(s) of the two models with altered directionality?
Run the reaction knockout analysis (RKO) on both the aerobic and the anaerobic model.
- How many reactions are essential in each model (aerobic vs anaerobic growth)?
- Are there reactions that are essential in only one model?
- Are there reactions leading to a beneficial, i.e. producing more biomass, knockout? If so, why?
- The following cytoplasmic reactions are essential for anaerobic growth:
|1 ADP + 1 H(+) + 1 D-fructofuranose 1,6-bisphosphate(4-)|| <-- ||1 ATP + 1 D-fructofuranose 6-phosphate|
|1 phosphoenolpyruvate + 1 CO(2) + 1 H2O|| --> ||1 Oxaloacetate + 1 H(+) + 1 hydrogenphosphate|
|1 glyceraldehyde 3-phosphate(2-) + 1 glycerone phosphate||<==>||1 D-fructofuranose 1,6-bisphosphate(4-)|
|1 D-glucopyranose 6-phosphate||<==>||1 D-fructofuranose 6-phosphate|
|1 glyceraldehyde 3-phosphate(2-)||<==>||1 glycerone phosphate|
However, there is no reaction that is essential for aerobic growth!
A more realistic model
- Add the model iAF1260 to your workspace and compare it to the bigg_e_coli_core model with respect to the growth reaction and the boundary conditions.
- Rerun the previous analyses and manipulations to get a feeling for a more realistic and complete model.