Pyrimidine and purine nucleoside phosphorylases catalyze the reversible phosphorolytic cleavage and formation of the glycosidic bond of purine and pyrimidine nucleosides, respectively, and are thus, key catalysts for the synthesis of new compounds. The selection of the best combination of enzyme and reaction conditions is not trivial, as each of the two enzymes can perform the reaction in two directions and thus, also competes with the other one for the reaction intermediates. A generic approach to the solution of this problem based on the formulation of a mixed integer dynamic optimization program using MOSAICmodeling is presented.