Redox cofactor utilization is one of the major barriers to the realization of efficient and cost-competitive cell-free biocatalysis, especially where multiple redox steps are concerned. The design of versatile, cofactor balanced modules for canonical metabolic pathways, such as glycolysis, is one route to overcoming such barriers. Here, we set up a computer-aided design framework to engineer the non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase (GapN) from Streptococcus mutans for enabling an NADH linked efficient cell-free glycolytic pathway with a net zero ATP usage. This rational design approach combines molecular dynamics simulations with a multistate computational design method that allowed us to consider different conformational states encountered along the GapN enzyme catalytic cycle. In particular, the cofactor flip, characteristic of this enzyme family and occurring before product hydrolysis, was taken into account to redesign the cofactor binding pocket for NAD + utilization. While GapN exhibits only trace activity with NAD + , a ∼10,000-fold enhancement of this activity was achieved, corresponding to a recovery of ∼72% of the catalytic efficiency of the wild-type enzyme on NADP + , with a GapN enzyme harboring only 5 mutations.