Practical 3

Linear Programming and Flux Balance Analysis with ScrumPy

In this practical you will try to repeat the results from Fell & Small (1986), which was one of the first attempts to apply LP to metabolic networks. Before staring with the exercises read the documentation of the ScrumPy LP module, here.

1. Download the original paper from this link.

2. Download and and open the adipose model (adipose.spy).

3. Your task is to set up and solve an LP problem where the objective is to minimise glucose consumption, while producing 1 unit of triacylglycerol (TGA).
   1. Start by importing the LP module and create the LP object with the adipose model as argument.
   2. Set the flux of the TGA-synthesizing reaction to 1.0, by entering this information to the appropriate LP method, i.e. we are going to optimise the production of one unit of triglyceride.

      1. Hint 1: the TGA-synthesizing reaction is R41.

      2. Hint 2: the argument to the intended LP method must be a dictionary where keys are reactions and values are flux values.
   3. Since the optimisation direction is minimisation, and this is set by default, you don't need to change the direction.
   4. Next, identify the reaction that should be minimised and enter this as an argument to the suitable method of the LP object (note that the argument must be a list of reaction(s), i.e. your argument should be a list with one reaction).

      • i.Hint: the glucose uptake reaction is named R1.

   5. Solve the LP problem using a suitable method.

   6. To obtain the solution use the LP method GetPrimSol(). This method returns a dictionary object of reactions in the solution as keys and flux values as values, so for convenience assign a name to this solution.

4. What is the objective value, i.e. the flux of the glucose uptake reaction, representing the minimum amount of glucose to make 1 triglyceride? How does this compare to the results in the paper?