In this work, we propose eigenvalue optimization combined with Lyapunov theory concepts to ensure stability of the Embden-Meyerhof-Parnas pathway, the pentosephosphate pathway, the phosphotransferase system and fermentation reactions of Escherichia coli. We address the design of a metabolic network for the maximization of different metabolite production rates. The first case study focuses on serine production, based on a model that consists of 18 differential equations corresponding to dynamic mass balances for extracellular glucose and intracellular metabolites, and thirty kinetic rate expressions. A second case study addresses the design problem to maximize ethanol production, based on a dynamic model that involves mass balances for 25 metabolites and 38 kinetic rate equations. The nonlinear optimization problem including stability constraints has been solved with reduced space Successive Quadratic Programming techniques. Numerical results provide useful insights on the stability properties of the studied kinetic models.