Optimization of covalent docking for organophosphates interaction with Anopheles acetylcholinesterase

Organophosphates (OPs) used as potent insecticides for malaria vector control, covalently phosphorylate the catalytic serine residue of Anopheles gambiae AChE (AgAChE) in a reaction that liberates their leaving groups. In the recent 10-year insecticide use assessment, OPs were the most frequently used World Health Organization prequalified insecticides. Molecular modelling programs are best suited to display molecular interactions between ligands and the target proteins. The docking modes that generate ligand poses closer to the binding site show high accuracy in predicting the ligand binding mode. The implicit solvation approach such as molecular mechanics-generalized born surface area (MM-GBSA) is a more reliable method to predict ligand onformations and binding affinities. Apart from covalent docking studies being scarce, current molecular docking programs do not adequately possess the covalent docking reaction algorithm to display the molecular mechanism of OPs at the AgAChE catalytic site. This results into OP docking studies commonly being conducted through noncovalent pannels. The aim of this study was to establish the optimim covalent docking system for OPs through manual customization of Schrödinger's Glide covalent docking reaction algorithm. To achieve this, a newly customized covalent reaction algorithm was assessed on a set of ligands covering aromatic, non-aromatic and hydrophobic OPs and compared to the noncovalent docking results in terms of reliability based on the reported X-ray diffraction molecular interactions and crystal poses. The study established that by virtue of omitting the well-known OP hydrolysis, noncovalent mode suggested molecular interactions that were further from the catalytic triad and could not otherwise occur when the molecule is hydrolyzed as in the customized covalent docking mode. Moreover, the MM-GBSA concurred with the optimized covalent docking in eliminating such inaccurate molecular interactions. Additionally, the covalent docking mode confined the interactions and ligand poses to the catalytic site indicating relatively high accuracy and reliability. This study reports the optimized covalent docking panel that effectively confirmed the molecular mechanisms of OPs, as well as indentifying the corresponding amino acid residues required to stabilize the aromatic, non-aromatic and hydrophobic OPs at the AgAChE catalytic site in line with the reported X-ray diffraction studies. As such, the proposed manual customization of the Schrödinger's Glide covalent docking platform can be used to reliably predict molecular interactions between OPs and AgAChE target.