Assessing the effect of a series of mutations on the dynamic behavior of phosphite dehydrogenase using molecular docking, molecular dynamics and quantum mechanics/molecular mechanics simulations

Abstract

Discovered in Pseudomonas stutzeri, phosphite dehydrogenase (PTDH) is an enzyme that catalyzes the oxidation of phosphite to phosphate while simultaneously reducing NAD+ to NADH. Despite several investigations into the mechanism of reaction and cofactor regeneration, only a few studies have focused on improving the activity and stability of PTDH. In this study, we combine molecular docking, molecular dynamics (MD) simulation, and Quantum Mechanics/Molecular Mechanics (QM/MM) to identify the impact of 30 mutations on the activity and stability of PTDH. Molecular docking results suggest that E266Q, K76A, K76M, K76R, K76C, and R237K can act on the NAD+ binding site through relatively weak bond development due to their high free binding energy. Moreover, Mulliken population analysis and potential energy barrier indicate that T101A, E175A, E175A/A176R, A176R, and E266Q act on phosphite oxidation. The mutants M53N, M53A, K76R, D79N, D79A, T101A, W134A, W134F Y139F, A146S, E175A, F198I, F198M, E266Q, H292K, S295A, R301K, and R301A were found to act on the structural dynamic of PTDH. The remaining mutants cause the loss of the nitrogen atom of R237 and H292, respectively, inactivating the enzyme. This study provides specific explanations of how mutations affect weak interactions of PTDH. The results should allow researchers to conduct experimental studies to improve PTDH activity and stability. Communicated by Ramaswamy H. Sarma