Biology Stack Exchange is a question and answer site for biology researchers, academics, and students. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

In my attempt to create a metabolic model of E. coli, I have found a comprehensive list of metabolic reactions and their stoichiometry. The one I am currently using is E. coli model iJO1366 (more here). What I have not yet found is a comprehensive list of kinetic rate constants associated with these reactions. I also require a data source of concentrations of these enzymes and metabolites at steady state in E. coli. Where can I find information on concentrations and rate constants?

share|improve this question
Perhaps the Brenda database might be of some help? – TomD Aug 23 '12 at 11:08
why do you want to compile this? Flux analysis is looking pretty useful and does not require all this data. – shigeta Sep 23 '12 at 4:09
I found SABIO: – Bitwise Mar 21 '13 at 1:02

I would advise you to use the SMBL for iAF1260 E. Coli with a COBRA toolbox on MATLAB. It has all the constants you are looking for optimized for the steady state and have been extensively tested for simulation of single and multiple genes deletion / mutation. Even if you want to stick with your model you will be able to copy the constants from there.

However depending on what you are doing you might not be satisfied by the whole models, since the they are computationally too complicated to simulate and optimize and lead to computationally unstable results, especially when you try to simulate colony collaborative mechanisms.

share|improve this answer
Do you have any literature to support your claim of computational instability? – Coder Mar 21 '13 at 18:02
This is my personal experience. When you use the COBRA toolbox, you use an external module for linear optimization to get your "objective" function. When I computed the base model from iAF1260 E.Coli, it worked well, just as a system of two bacteria living together. But as soon as you started to perform single or double gene deletion analysis, your objective function maximum started to behave erratically. You could make it bit better with it by using a different external linear optimization module, but still no analysis was possible. – Andrei Kucharavy Mar 21 '13 at 18:52
This work was based on this paper and I had an impression that the authors run into pretty much the same limitations: only qualitative analysis is provided with the computational model, whereas they could have gone quantitative with this model and had quantitative experimental results. – Andrei Kucharavy Mar 21 '13 at 18:57

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.