Enhancing Paraoxon Binding to Organophosphorus Hydrolase Active Site

Organophosphorus (OP) compounds are among the most toxic of chemical substances and widely used as insecticides, pesticides, and chemical warfare agents. The most important enzyme inhibited by OP compounds is acetylcholinesterase (AChe). Inactivation of AChe function results in the accumulation of neurotransmitter, leading to death due to serious respiratory disorders. Organophosphorus hydrolase (OPH), also called phosphotriesterase, is a homo-dimeric metalloenzyme that can hydrolyze various OP agents in the circulatory system, resulting in products that are generally of reduced toxicity. The best OPH substrate found to date is the insecticide diethyl p-nitrophenyl phosphate (paraoxon). Most structural and kinetic studies assume that the binding orientation of paraoxon is identical to that of diethyl 4-methylbenzylphosphonate, which is the only substrate analog co-crystallized with OPH. In the current work, we used a combined docking and molecular dynamics (MD) approach to predict the likely binding mode of paraoxon in the OPH active site. We identified a potential binding mode of paraoxon that does not match the binding mode of diethyl 4-methylbenzylphosphonate. Then, we used the predicted binding mode to run MD simulations on the wild type (WT) OPH complexed with paraoxon, and OPH mutants complexed with paraoxon. Additionally, we identified 3 hot-spot residues (D253, H254, and I255) involved in the stability of the OPH active site. To further assess these predictions, we then experimentally assayed single and double mutants involving these residues (D253E, H254S, I255S, D253E-H254R and D253E-I255G) for hydrolytic activity against paraoxon. Computational structural analysis of protein-substrate dynamics shows different hydrogen bonding profiles for mutants involving D253 (D253E, D253E-H254R, and D253E-I255G) compared to WT OPH. Additionally, the binding free energy calculations and the experimental kinetics (particularly, k_cat and K_M) of the reactions between each OPH mutant and paraoxon show that mutated forms D253E, D253E-H254R, and D253E-I255G exhibit enhanced activity over WT OPH. Interestingly, our experimental results show that the activity of the double mutant D253E-H254R increased by 19-fold compared to WT OPH.