Cyclohexanone monooxygenases (CHMO) consume molecular oxygen and NADPH to catalyze the valuable oxidation of cyclic ketones. However, CHMO usage is restricted by poor stability and stringent specificity for NADPH. Efforts to engineer CHMO have been limited by the sensitivity of the enzyme to perturbations in conformational dynamics and long-range interactions that cannot be predicted. We demonstrate an aerobic, high-throughput growth selection platform in Escherichia coli for oxygenase evolution based on NADH redox balance. We applied this NADH-dependent selection to alter the cofactor specificity of CHMO to accept NADH, a less expensive cofactor than NADPH. We first identified the variant CHMO DTNP (S208D-K326T-K349N-L143P) with a ∼1200-fold relative cofactor specificity switch from NADPH to NADH compared to the wild type through semirational design. Molecular modeling suggests CHMO DTNP activity is driven by cooperative fine-tuning of cofactor contacts. Additional evolution of CHMO DTNP through random mutagenesis yielded the variant CHMO DTNPY with a ∼2900-fold relative specificity switch compared to the wild type afforded by an additional distal mutation, H163Y. These results highlight the difficulty in engineering functionally innovative variants from static models and rational designs, and the need for high throughput selection methods. Our introduced tools for oxygenase engineering accelerate the advancements of characteristics essential for industrial feasibility.