Coupling constrained-based flux sampling and clustering to tackle cancer metabolic heterogeneity