`removeReaction!` etc. doesn't update all fields in a `CoreModel`

LCSB-BioCore / COBREXA.jl

Constraint-Based Reconstruction and EXascale Analysis

https://lcsb-biocore.github.io/COBREXA.jl/

Apache License 2.0

42 stars 8 forks source link

`removeReaction!` etc. doesn't update all fields in a `CoreModel` #768

Open htpusa opened 1 year ago

htpusa commented 1 year ago

using COBREXA
model = load_model(CoreModel,"test/downloaded/e_coli_core.xml")
removeReaction!(model,1)
n_reactions(model) # 94
length(model.grrs) # 95

Only some of the fields are updated so, for example, the gene associations end up not matching anymore.

exaexa commented 1 year ago

ah this. We added grrs later so this is an avoidable but kinda expected bug. Thanks for report