A mechanism-based 3D-QSAR approach for classification and prediction of acetylcholinesterase inhibitory potency of organophosphate and carbamate analogs.

Organophosphate (OP) and carbamate esters can inhibit acetylcholinesterase (AChE) by binding covalently to a serine residue in the enzyme active site, and their inhibitory potency depends largely on affinity for the enzyme and the reactivity of the ester. Despite this understanding, there has been no mechanism-based in silico approach for classification and prediction of the inhibitory potency of ether OPs or carbamates. This prompted us to develop a three dimensional prediction framework for OPs, carbamates, and their analogs. Inhibitory structures of a compound that can form the covalent bond were identified through analysis of docked conformations of the compound and its metabolites. Inhibitory potencies of the selected structures were then predicted using a previously developed three dimensional quantitative structure-active relationship. This approach was validated with a large number of structurally diverse OP and carbamate compounds encompassing widely used insecticides and structural analogs including OP flame retardants and thio- and dithiocarbamate pesticides. The modeling revealed that: (1) in addition to classical OP metabolic activation, the toxicity of carbamate compounds can be dependent on biotransformation, (2) OP and carbamate analogs such as OP flame retardants and thiocarbamate herbicides can act as AChEI, (3) hydrogen bonds at the oxyanion hole is critical for AChE inhibition through the covalent bond, and (4) π-π interaction with Trp86 is necessary for strong inhibition of AChE. Our combined computation approach provided detailed understanding of the mechanism of action of OP and carbamate compounds and may be useful for screening a diversity of chemical structures for AChE inhibitory potency.