fix: getStandardKcat for light ecModels

SysBioChalmers · Jun 6, 2024 · ce31bf7 · ce31bf7
1 parent f5ea969
commit ce31bf7
Showing 1 changed file with 6 additions and 1 deletion.
diff --git a/src/geckomat/gather_kcats/getStandardKcat.m b/src/geckomat/gather_kcats/getStandardKcat.m
@@ -134,7 +134,12 @@
 % Find reactions with GPR but without model.ec entry (for instance due to
 % no protein matching)
 rxnsMissingEnzyme = find(~cellfun(@isempty, model.grRules));
-rxnsMissingEnzyme = find(and(~ismember(model.rxns(rxnsMissingEnzyme),model.ec.rxns), ~contains(model.rxns(rxnsMissingEnzyme),'usage_prot_')));
+if model.ec.geckoLight
+    ecRxnsList = unique(extractAfter(model.ec.rxns,4));
+else
+    ecRxnsList = model.ec.rxns;
+end
+rxnsMissingEnzyme = find(and(~ismember(model.rxns(rxnsMissingEnzyme),ecRxnsList), ~contains(model.rxns(rxnsMissingEnzyme),'usage_prot_')));
 rxnsMissingGPR = [rxnsMissingGPR;rxnsMissingEnzyme];
 
 % Get custom list of reaction IDs to ignore, if existing. First column