getReactionsFromEnzyme

 getReactionsFromEnzyme
   Get all reactions that are annotated to a particular enzyme.

 Input:
   ecModel     an ecModel in GECKO 3 format (with ecModel.ec structure)
   proteinId   protein identifier, matching ecModel.ec.enzymes.

 Output:
   rxns        reactions that are associated with this enzyme
   kcat        kcat values of the corresponding reactions
   idx         index of the reactions in ecModel.ec.rxns
   rxnNames    names of the reactions
   grRules     grRules of the reactions

 Usage: [rxns, kcat, idx, rxnNames, grRules] = getReactionsFromEnzyme(ecModel,proteinId)

0001 function [rxns, kcat, idx, rxnNames, grRules] = getReactionsFromEnzyme(ecModel,proteinId)
0002 % getReactionsFromEnzyme
0003 %   Get all reactions that are annotated to a particular enzyme.
0004 %
0005 % Input:
0006 %   ecModel     an ecModel in GECKO 3 format (with ecModel.ec structure)
0007 %   proteinId   protein identifier, matching ecModel.ec.enzymes.
0008 %
0009 % Output:
0010 %   rxns        reactions that are associated with this enzyme
0011 %   kcat        kcat values of the corresponding reactions
0012 %   idx         index of the reactions in ecModel.ec.rxns
0013 %   rxnNames    names of the reactions
0014 %   grRules     grRules of the reactions
0015 %
0016 % Usage: [rxns, kcat, idx, rxnNames, grRules] = getReactionsFromEnzyme(ecModel,proteinId)
0017 
0018 protIdx     = find(strcmpi(ecModel.ec.enzymes,proteinId));
0019 ecRxnIdx    = find(ecModel.ec.rxnEnzMat(:,protIdx));
0020 rxns        = ecModel.ec.rxns(ecRxnIdx);
0021 kcat        = ecModel.ec.kcat(ecRxnIdx);
0022 
0023 idx         = ecRxnIdx;
0024 
0025 [~,rxnIdx]  = ismember(rxns,ecModel.rxns);
0026 rxnNames    = ecModel.rxnNames(rxnIdx);
0027 grRules     =  ecModel.grRules(rxnIdx);
0028 end