Genome-scale metabolic networks reconstructions from different organisms have become popular in recent years. Genetic engineering is proven to be able to obtain the desirable phenotypes. Optimization algorithms are implemented in previous works to identify the effects of gene knockout on the results. However, the previous works face the problem of falling into local minima. Thus, a hybrid of Bees Algorithm and Flux Balance Analysis (BAFBA) is proposed in this paper to solve the lo...