-
Notifications
You must be signed in to change notification settings - Fork 49
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
bug: full ecHumanGEM has multiple usage reactions for the same proteins #353
Comments
The problem is that |
@wshao1 How do you generate adapterLocation = fullfile(findGECKOroot,'tutorials','light_ecModel','HumanGEMAdapter.m');
adapter = ModelAdapterManager.setDefault(adapterLocation);
model = loadConventionalGEM();
[ecModel, noUniprot] = makeEcModel(model,false);
find(strcmp(ecModel.rxns,'usage_prot_Q0WX57'))
ans =
0×1 empty double column vector |
The only differences between my adapter and the one provided in the
It is unclear to me why this simplification is required, and I could not find the |
I can replicate this with Human-GEM v1.16. The problem is with reaction MAR09490 whose gene association contains multiple isoenzymes that are all annotated to the same UniProt ID. As a work-around, the duplicated
The above code should work, but I haven't tested the resulting model with In the meanwhile, we'll work on changes to Regarding |
Thanks @edkerk, removing the duplicated reactions allows The resulting model shows zero growth rate however.
I have also tried extracting an ecHT29 using the HT29-GEM from the
Integrating proteomics data and then flexibilizing enzyme concentrations results in the following error, suggesting getSubsetEcModel may have returned a non-functional model.
I will investigate this further. |
The context-specific ecModel failing to simulate is likely the result of the large ecModel and the context-specific non-ec GEM that are used as input by This is now clarified in the function documentation, and |
Description of the bug:
I have an ecHumanGEM created using the full_ecModel protocol, which I want to reduce to a context-specific GEM using
getSubsetEcModel
. However, this returns an error due to the ecHumanGEM having multiple protein usage reactions with the same ID (i.e. 'usage_prot_Q0WX57'). This results in the following traceback:Reproducing these results:
System information
I hereby confirm that:
main
branch of the repository.The text was updated successfully, but these errors were encountered: