|
Size: 2136
Comment:
|
Size: 2057
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 8: | Line 8: |
| 1. Download and open the adipose model ([[http://mudsharkstatic.brookes.ac.uk/Download/adipose.spy|adipose.spy]]). | 2. Download and open the adipose model ([[http://mudsharkstatic.brookes.ac.uk/Download/adipose.spy|adipose.spy]]). |
| Line 10: | Line 10: |
| 1. 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). | 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). |
| Line 13: | Line 13: |
| a. 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. | b. Set the flux of the TGA-synthesising 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. [[AccliPhot/WorkshopOne/prac3/hint_1 | Hint 1]] |
| Line 15: | Line 16: |
| i. Hint 1: the TGA-synthesizing reaction is ''R41''. | [[AccliPhot/WorkshopOne/prac3/hint_2 | Hint 2]] |
| Line 17: | Line 18: |
| i. Hint 2: the argument to the intended LP method must be a dictionary where keys are reactions and values are flux values. | |
| Line 19: | Line 19: |
| a. Since the optimisation direction is minimisation, and this is set by default, you don't need to change the direction. | c. Since the optimisation direction is minimisation, and this is set by default, you don't need to change the direction. |
| Line 21: | Line 21: |
| a. 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). | d. 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). |
| Line 25: | Line 25: |
| a. Solve the LP problem using a suitable method. | e. Solve the LP problem using a suitable method. |
| Line 27: | Line 27: |
| a. 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. | f. 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. |
| Line 29: | Line 29: |
| 1. 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? | 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? |
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.
Download the original paper from this link.
Download and open the adipose model (adipose.spy).
- 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).
- Start by importing the LP module and create the LP object with the adipose model as argument. b. Set the flux of the TGA-synthesising 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. c. Since the optimisation direction is minimisation, and this is set by default, you don't need to change the direction. d. 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.
f. 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.
- Start by importing the LP module and create the LP object with the adipose model as argument. b. Set the flux of the TGA-synthesising 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. c. Since the optimisation direction is minimisation, and this is set by default, you don't need to change the direction. d. 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).
- 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?