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 = addReaction(copia,'ED1',{'6pgc[c]','h2o[c]','KDPG[c]'},[-1 1 1],true,0,50); %copia = addReaction(copia,'ED2',{'KDPG[c]','pyr[c]','g3p[c]'},[-1 1 1],true,0,50); 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,'ATPM',0,'b'); %copia = changeRxnBounds(copia,'CS',0.3*10*8/7,'u'); %copia = changeRxnBounds(copia,'PFL',1,'b'); %copia = changeRxnBounds(copia,'ME1',0,'b'); %copia = changeRxnBounds(copia,'G6PDH2r',0,'b'); %copia = changeRxnBounds(copia,'PGI',0,'b'); %copia = changeRxnBounds(copia,'PGL',2,'b'); %copia = changeRxnBounds(copia,'THD2',0,'u'); %copia = changeRxnBounds(copia,'NADTRHD',0,'u'); %copia = changeRxnBounds(copia,'EX_ac(e)',0.1,'b'); %copia = changeRxnBounds(copia,'GND',0,'b'); %copia = changeRxnBounds(copia,'EX_etoh(e)',0.1,'b'); %copia = changeRxnBounds(copia,'Ex_cit[c]',0.23,'b'); %copia = changeRxnBounds(copia,'EX_succ(e)',0.15,'b'); %copia = changeRxnBounds(copia,'EX_Pyr(e)',0.01,'b'); copia = changeRxnBounds(copia,'Biomass_Ecoli_core_w_GAM',0,'b'); %copia = changeRxnBounds(copia,'EX_lac_D(e)',0.05,'b'); %copia = changeRxnBounds(copia,'ppdiolDH',5,'u'); %copia = changeRxnBounds(copia,'EX_o2(e)',-2,'b'); %mutante = addReaction(mutante,'transporter',{'glc_D[e]','glc[c]'},[-1 1],true,0,20); %wildtype={'nadp[c]';'nadph[c]'}; %mutant={'nad[c]';'nadh[c]'}; %copia = changeRxnMets(copia,wildtype,mutant,'GND'); %copia = changeRxnMets(copia,wildtype,mutant,'ICDHyr'); %copia = changeRxnMets(copia,wildtype,mutant,'G6PDH2r'); %mutante = addReaction(copia,'NADG6PDH',{'g6p[c]','nad[c]','nadh[c]','6pgl[c]','h[c]'},[-1 -1 1 1 1],true,0,10); %copia = changeObjective(copia,'ATPM'); %copia = changeObjective(copia,'EX_o2(e)'); copia = changeObjective(copia,'EX_13PD[c]'); FBAsolution = optimizeCbModel(copia,'max'); %FBAsolution = optimizeCbModel(copia,'min'); FBAsolution.f %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %Reações que vou "imprimir -- em formato de linha" 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