Double and multiple knockout simulations for genome-scale metabolic network reconstructions