Exploring stereochemical specificity of phosphotriesterase by MM-PBSA and MM-GBSA calculation and steered molecular...

Wild-type phosphotriesterase (PTE) prefers the S_P-enantiomers over the corresponding R_P-enantiomers by factors ranging from 10 to 90. To satisfy the binding modes of the PTE of S_P- and R_P-enantiomers, all-atom molecular dynamics simulations were carried out on two paraoxon S_P and R_P derivatives, namely, S_p-1 and R_p-1. Molecular mechanics Poisson-Boltzmann surface area and molecular mechanics generalized Born surface area (MM-PBSA and MM-GBSA) calculations indicated that His230 in S_p-1-PTE had a closer interaction with the substrate than that in R_p-1-PTE and that such interaction increased the catalytic efficiency of PTE for S_p-1. The steered molecular dynamics simulation indicated that, compared with S_p-1, R_p-1 in the unbinding (binding) may hinder some residue displacement, thus requiring more effort to escape the binding pocket of PTE. In addition, Trp131, Phe306, and Tyr309 are deemed important residues for the S_p-1 unbinding pathway via PTE, whereas Tyr309 alone is considered an important residue for the R_p-1 unbinding pathway. These results demonstrate the possibility of dramatically altering the stereoselectivity and overall reactivity of the native enzyme toward chiral substrates by modifying specific residues located within the active site of PTE.