1Physiological states-derived metabolic regulation network determines cellular phenotype but 2 is not considered in current constraint-based models (CBMs). Here, we proposed a novel 3 metabolic flux prediction model, Decrem, which integrated comprehensive regulatory 4 contexts-topologically coupled reactions-derived coactivation regulation and global cell 5 state regulators-derived enzyme kinetics-into canonical CBMs to approximate biologically 6 feasible fluxes. A unique advantage of Decrem is it relies only on the concentrations of 7 identified metabolites. When applied to three model organisms, Escherichia coli, 8Saccharomyces cerevisiae and Bacillus subtilis, Decrem demonstrated high accuracy in 9 metabolic flux prediction across various conditions, particularly for incomplete metabolic 10 models. Growth analysis by Decrem on yeast knockout strains revealed specific overflowed 11 pathway for several low growth rate strains, which are consistent with experiments but not 12 identified by previous methods. Additionally, the Pearson correlation coefficients between 13 experimental and predicted growth rates on 1030 E. coli genome-scale deletion strains were 14 increased to 0.743, compared to the 0.127, 0.103 and 0.281, respectively. This method is 15 expected to accurately simulate dynamic metabolic regulation and phenotype prediction for 16 synthetic biology. 17Keywords: genome-scale growth rate analysis / hierarchical metabolic regulation model / 18 linearized michaelis-menten kinetics / metabolic flux prediction / sparse linear basis 19 decomposition 20