-
Notifications
You must be signed in to change notification settings - Fork 0
/
example2.m
37 lines (36 loc) · 1.06 KB
/
example2.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
function [] = example2()
% load the ecoli core model.
load('ecoli_core_model.mat')
%
% GridProd is employed.
% The target metabolite is "succ[c]".
% The glucose reaction is 'EX_glc(e)'.
% The oxygen reaction is 'EX_o2(e)'.
% The biomass objective fucntion reaction is
% 'Biomass_Ecoli_core_w_GAM'.
%
% The options are specified as follows.
% the glucose uptake ratio is 8,
% the oxygen uptake ratio is 18.5,
% the minimum growth ratio is 0.01.
% P is 10.
%
[targetProduction, minFlux, maxFlux, blockedRxns, usedRxns, biomass]=...
GridProd(model,{'succ[c]'},'EX_glc(e)','EX_o2(e)','Biomass_Ecoli_core_w_GAM',...
'GUR',8,'OUR',18.5,'minGrowth',0.01,'P',10);
minFlux
maxFlux
%
% 9.9904 is obtained as the minimum and maximum production rates for succinate.
%
save('example2.mat');
paintGrid('results/analyzeResult_69.mat')
%
% 'paintGrid' shows the heatmap that shows the levels of succinate
% production for each grid.
% Note that 69 is ID of succinate in the ecoli core model.
save('example2.mat');
%
% Jul. 20, 2017, Takeyuki TAMURA
%
end