Organophosphorus compound esterase profiles as predictors of therapeutic and toxic effects

Certain organophosphorus compounds (OPCs) inhibit various serine esterases (EOHs) via phosphorylation of their active site serines. We focused on 4 EOHs of particular toxicological interest: acetylcholinesterase (AChE: acute neurotoxicity; cognition enhancement), butyrylcholinesterase (BChE: inhibition of drug metabolism and/or stoichiometric scavenging of EOH inhibitors; cognition enhancement), carboxylesterase (CaE: inhibition of drug metabolism and/or stoichiometric scavenging of EOH inhibitors), and neuropathy target esterase (NTE: delayed neurotoxicity, OPIDN). The relative degree of inhibition of these EOHs constitutes the “esterase profile” of an OPC and serves as a major determinant of its net physiological effects. Thus, understanding and controlling the esterase profile of OPC activity and selectivity toward these 4 target enzymes is a significant undertaking. In the present study, we analyzed the inhibitor properties of 52 OPCs against the 4 EOHs, along with pairwise and multitarget selectivities between them, using 2 QSAR approaches: Hansch modeling and Molecular Field Topology Analysis (MFTA). The general formula of the OPCs was (RO)2P(O)X, where R = alkyl, X = - SCH(Hal)COOEt (Hal = Cl, Br), -SCHCl2, -SCH2Br, -OCH(CF3)R1 (R1 = C6H5, CF3, COOEt, COOMe). The Hansch model showed that increasing neuropathic potential correlated with rising R hydrophobicity; moreover, OPC binding to scavenger EOHs (BChE and CaE) had different effects on potential acute and delayed neurotoxicity. Predicted protective roles of BChE and CaE against acute toxicity were enhanced with increasing hydrophobicity, but projected protection against OPIDN was decreased. Next, Molecular Field Topology Analysis (MFTA) models were built, considering atomic descriptors, e.g., effective charge, van der Waals radius of environment, and group lipophilicity. Activity/selectivity maps confirmed predictions from Hansch models and revealed other structural factors affecting activity and selectivity. Virtual screening based on multitarget selectivity MFTA models was used to design libraries of OPCs with favorable esterase profiles for potential application as selective inhibitors of CaE without untoward side effects.

Highlights► Esterase profiles of organophosphorus compounds were analyzed by two QSAR models. ► Molecular Field Topology Analysis (MFTA) and Hansch analysis give consistent results. ► Major factors influencing activity and selectivity were identified. ► Selective carboxylesterase inhibitors were designed based on MFTA models.