clear
load modelo2
teste = model;

%definição das restrições
teste = changeRxnBounds(teste,'EX_glc(e)',-8.7,'l')
%teste = changeRxnBounds(teste,'THD2',0,'b')
%teste = changeRxnBounds(teste,'EX_ac(e)',0.69*8.4,'l')
%teste = changeRxnBounds(teste,'G6PDH2r',0.27*7.8,'u')
%teste = changeRxnBounds(teste,'ME2',0.07*7.8,'b')
%teste = changeRxnBounds(teste,'ICDHyr',3,'u')
%teste = changeRxnBounds(teste,'PGI',0,'b');
teste = changeRxnBounds(teste,'EX_o2(e)',-11.9,'b');
%teste = changeRxnBounds(teste,'NADTRHD',0,'b');
%wt={'nadp[c]','nadph[c]'};
%mut={'nad[c]','nadh[c]'};
%teste=changeRxnMets(model,wt,mut,'G6PDH2r');
%teste = changeRxnBounds(teste,'EX_co2(e)',(20),'l');

%Definição da função objetivo

teste = changeObjective(teste,'Biomass_Ecoli_core_w_GAM');

%Execução do procedimento de otimização

FBAsolution = optimizeCbModel(teste,'max');
FBAsolution.x;