The development and application of a novel in vivo selection method enabled the testing of more than >106 variants of formate dehydrogenase (FDH) to switch its cofactor specificity in one single round of protein directed evolution, using multiple mutational sites selected based on previous studies and initial computational modeling. Due to this new approach, key