Metabolic analysis

  • For this task you have to install a Python package for constraint based modeling of metabolism, called cobrapy. Load the EColi reconstructed model. Use this documentation link to get the metabolic reactions and the genes associated to them.

  • Read the reactions into a network having the metabolites as nodes, and find out the most central metabolites both in terms of degree and in terms of betweeness centrality. Find all genes associated to reactions containing them. Perform GO enrichment on the set of genes looking in particular at cellular components. Are the central metabolites active in particular locations inside the cell?

Extra points:

  • Run a flux balance analysis and put the flux values into the network as well. How are the results changing?

In [ ]: