We used molecular docking to generate enzyme-cofactor complexes, MD simulations for the conformational sampling of wild type (WT) and mutant complexes, Density Functional Theory (DFT) calculations for characterization of electronic states, and hybrid QM/MM calculations for exploring the catalytic mechanism. Each of these steps is discussed in detail in the subsequent section.