function outModel=removeDrugReactions(inModel)
% outModel=removeDrugReactions(inModel)
%
% Removes reactions from Human-GEM related to drugs
% 
% inModel       The model, typically Human-GEM or Mouse-GEM
% outModel      The resulting model.

drugNames = {'pravastatin';'Gliclazide';'atorvastatin';'fluvastatin';'fluvastain';'fluvstatin';'simvastatin';'cyclosporine';...
             'acetaminophen';'cerivastatin';'Tacrolimus';'ibuprofen';'lovastatin';'Losartan';'nifedipine';...
             'pitavastatin';'rosuvastatin';'Torasemide';'Midazolam'};
drugMets = contains(inModel.metNames,drugNames,'IgnoreCase',true);
drugRxns = (sum(inModel.S(drugMets,:) ~= 0,1) > 0).';
outModel = removeReactions(inModel,inModel.rxns(drugRxns));
end

SysBioChalmers
github-repos/MATLAB/SysBioChalmers-Human-GEM/Human-GEM-0b1bd42adaa2e1d7ac52ee83b989fad8a695759d/code/tINIT/removeDrugReactions.m
Human-GEM
github-repos/MATLAB/SysBioChalmers-Human-GEM
0b1bd42adaa2e1d7ac52ee83b989fad8a695759d