assignGPRA Assign each reaction a Gene-protein-reaction association model = assignGPRA(model,gpraFile) INPUTS model COBRA model structure gpraFile SimPheny GPRA file OUTPUT model COBRA model with gene-protein-reaction assoction 10/4/06 Markus Herrgard
0001 function model = assignGPRA(model,gpraFile) 0002 %assignGPRA Assign each reaction a Gene-protein-reaction association 0003 % 0004 % model = assignGPRA(model,gpraFile) 0005 % 0006 %INPUTS 0007 % model COBRA model structure 0008 % gpraFile SimPheny GPRA file 0009 % 0010 %OUTPUT 0011 % model COBRA model with gene-protein-reaction assoction 0012 % 0013 % 10/4/06 Markus Herrgard 0014 0015 [rxnList,genes,rxnGeneMat,subSystems] = parseSimPhenyGPRA(gpraFile); 0016 0017 [tmp,gpraRxnInd,modelRxnInd] = intersect(rxnList,model.rxns); 0018 0019 model.genes = genes; 0020 0021 model.rxnGeneMat = sparse(length(model.rxns),length(genes)); 0022 model.rxnGeneMat(modelRxnInd,:) = rxnGeneMat(gpraRxnInd,:);