New methods in computation of reaction fluxes from metabolomics data