Date of this Version
Schroeder & Saha, iScience 23, 100783
Stoichiometric metabolic modeling, particularly genome-scale models (GSMs), is now an indispens- able tool for systems biology. The model reconstruction process typically involves collecting informa- tion from public databases; however, incomplete systems knowledge leaves gaps in any reconstruc- tion. Current tools for addressing gaps use databases of biochemical functionalities to address gaps on a per-metabolite basis and can provide multiple solutions but cannot avoid thermodynami- cally infeasible cycles (TICs), invariably requiring lengthy manual curation. To address these limita- tions, this work introduces an optimization-based multi-step method named OptFill, which performs TIC-avoiding whole-model gapfilling. We applied OptFill to three fictional prokaryotic models of increasing sizes and to a published GSM of Escherichia coli, iJR904. This application resulted in holistic and infeasible cycle-free gapfilling solutions. In addition, OptFill can be adapted to automate inherent TICs identification in any GSM. Overall, OptFill can address critical issues in automated development of high-quality GSMs.