Towards stable kinetics of large metabolic networks: Nonequilibrium potential function approach