To refine the draft model, the gapFind [56] program in the COBRA Toolbox [22] was used to identify metabolic gaps in draft model and literature data were used to fill these gaps. The metabolites in each reaction were characterized based on their chemical formulae and neutral charges, which were obtained using CHEBI [57] and PubChem [58].