NADPH regeneration capacity is attracting growing research attention due to its important role in resisting oxidative stress. Besides, NADPH availability has been regarded as a limiting factor in production of industrially valuable compounds. The central carbon metabolism carries the carbon skeleton flux supporting the operation of NADPH-regenerating enzyme and offers flexibility in coping with NADPH demand for varied intracellular environment. To acquire an insightful understanding of its NADPH regeneration capacity, the elementary mode method was employed to compute all elementary flux modes (EFMs) of a network representative of central carbon metabolism. Based on the metabolic flux distributions of these modes, a cluster analysis of EFMs with high NADPH regeneration rate was conducted using the self-organizing map clustering algorithm. The clustering results were used to study the relationship between the flux of total NADPH regeneration and the flux in each NADPH producing enzyme. The results identified several reaction combinations supporting high NADPH regeneration, which are proven to be feasible in cells via thermodynamic analysis and coincident with a great deal of previous experimental report. Meanwhile, the reaction combinations showed some common characteristics: there were one or two decarboxylation oxidation reactions in the combinations that produced NADPH and the combination constitution included certain gluconeogenesis pathways. These findings suggested cyclization pathways as a powerful way for NADPH regeneration capacity of bacterial central carbon metabolism.