clear
load modelo2;
copia = model;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

copia = addReaction(copia,'GLYCDH',{'glyc[c]','nad[c]','dha[c]','nadh[c]','h[c]'},[-1 -1 1 1 1],true,-20,20);
copia = addReaction(copia,'DHAK',{'dha[c]','atp[c]','dhap[c]','adp[c]','h[c]'},[-1 -1 1 1 1],true,-20,20);
copia = addReaction(copia,'glycdeshidratase',{'glyc[c]','h2o[c]','3HP[c]'},[-1 1 1],true,0,20);
copia = addReaction(copia,'ppdiolDH',{'3HP[c]','nadh[c]','h[c]','13PD[c]','nad[c]'},[-1 -1 -1 1 1],true,0,20);

copia = addExchangeRxn(copia,{'glyc[c]'},-20,20);
copia = addExchangeRxn(copia,{'13PD[c]'},0,20);

copia = changeRxnBounds(copia,'EX_glc(e)',0,'b');
copia = changeRxnBounds(copia,'EX_glyc[c]',-10,'l');
copia = changeRxnBounds(copia,'Biomass_Ecoli_core_w_GAM',0,'b');


copia = changeObjective(copia,'EX_13PD[c]');

for i=0:20
copia = changeRxnBounds(copia,'EX_o2(e)',-1*0.5*i,'b');
FBAsolution = optimizeCbModel(copia,'max');
rendimentos(i+1)=FBAsolution.x(findRxnIDs(copia,'ppdiolDH'));
oxigenio(i+1)=-1*FBAsolution.x(findRxnIDs(copia,'EX_o2(e)'));
end

plot(oxigenio,rendimentos)
FBAsolution.f

%%%%%%%%%%%%%%%%%%%%

reac = {'ICL', 'MALS', 'ME1', 'ME2', 'PPC', 'PPCK', 'CS', 'ICDHyr', 'AKGDH', 'SUCDi', 'MDH', 'GLCpts', 'PGI', 'PFK', 'FBA', 'TPI', 'GAPD', 'ENO', 'PYK', 'PDH', 'PPS', 'FBP', 'ATPS4r', 'NADTRHD', 'THD2', 'G6PDH2r', 'GND', 'RPE', 'RPI', 'TALA', 'TKT1', 'TKT2', 'PFL', 'ALCD2x', 'PTAr', 'EX_lac_D(e)', 'Biomass_Ecoli_core_w_GAM', 'GLYCDH', 'ppdiolDH', 'ATPM', 'CYTBD','EX_o2(e)'};
rxnID = findRxnIDs(copia,reac);

%Salvando os fluxos
for i=1:length(rxnID)
    
    fluxos(i)=FBAsolution.x(rxnID(i));
    
end

% CAPTURA DA IMAGEM
%img = imread('SemEDingles.jpg');
img = imread('MetabEcolisemED1.jpg');
hFig = figure;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% SUPRIME APARICÃO DE UM WARNING
s = warning('off', 'Images:initSize:adjustingMag');
imshow(img, 'InitialMagnification',100, 'Border','tight')
warning(s);
%image(img);

% CONVERSAO DOS VALORES NUMERICOS DOS FLUXOS EN CARACTERES PARA SEREM IMPRESSOS

for i = 1:length(reac)
    
    carteles{i} = num2str(sprintf('%g',round(fluxos(i)*100)/100)); 

end

% COMANDOS PARA CONSEGUIR A IMPRESSAO DOS FLUXOS ACIMA DA FIGURA
for i=1:length(reac)
text(valorx(i),valory(i),carteles(i));
end

% TROCANDO AS PROPRIEDADES DOS TEXTOS QUE SAEM NA FIGURA

figureHandle = gcf;

set(findall(figureHandle,'type','text'),'Color','black','fontSize',8.5,'fontWeight','bold');


% Apagando TODAS AS VARIAVEIS TEMPORAIS QUE NAO INTERESA CONSERVAR
clear ans reac rxnID figureHandle carteles fluxos hFig i img s FBAsolution copia valores