Due to its importance in the pharmaceutical industry, ligand dynamic simulations have experienced a great expansion. Using all-atom models and cutting edge hardware, one can perform non-biased ligand migration, active site search and binding studies. In this letter we demonstrate (and validate by PCR mutagenesis) how these techniques, when combined with quantum mechanics, open new possibilities in enzyme engineering. We provide a complete analysis where: 1) biophysical simulations produce ligand diffusion and, 2) biochemical modeling samples the chemical event. Using such broad analysis we engineer a highly stable peroxidase activating the enzyme for new substrate oxidation after rational mutation of two non-conserved surface residues. In particular, we create a new surface-binding site, quantitatively predicting the in vitro change in oxidation rate obtained by mutagenic PCR and achieving a comparable specificity constant to active peroxidases.

This work was supported by the INDOX (KBBE-2013-7-613549 to ATM) European project, and the CTQ2013-48287 (to VG) and BIO2014-56388-R (to FJR-D) projects of the Spanish Ministry of Economy and Competitiveness (MINECO). FJR-D acknowledges a MINECO Ramón&Cajal contract.

Peer Reviewed