ABSTRACT
A series of dialkyl aryl phosphates and dialkyl arylalkyl phosphates were synthesized. Their inhibitory activities were evaluated against acetylcholinesterase (AChE) and butyrylcholinesterase (BChE). The di-n-butyl phosphate series consistently displayed selective inhibition of BChE over AChE. The most potent inhibitors of butyrylcholinesterase were di-n-butyl-3,5-dimethylphenyl phosphate (4b) [KI=1.0±0.4µM] and di-n-butyl 2-naphthyl phosphate (5b) [KI=1.9±0.4µM]. Molecular modeling was used to uncover three subsites within the active site gorge that accommodate the three substituents attached to the phosphate group. Phosphates 4b and 5b were found to bind to these three subsites in analogous fashion with the aromatic groups in both analogs being accommodated by the "lower region," while the lone pairs on the PO oxygen atoms were oriented towards the oxyanion hole. In contrast, di-n-butyl-3,4-dimethylphenyl phosphate (4a) [KI=9±1µM], an isomer of 4b, was found to orient its aromatic group in the "upper left region" subsite as placement of this group in the "lower region" resulted in significant steric hindrance by a ridge-like region in this subsite. Future studies will be designed to exploit these features in an effort to develop inhibitors of higher inhibitory strength against butyrylcholinesterase.
Subject(s)
Butyrylcholinesterase/metabolism , Cholinesterase Inhibitors/chemistry , Cholinesterase Inhibitors/pharmacology , Naphthalenes/chemistry , Naphthalenes/pharmacology , Organophosphates/chemistry , Organophosphates/pharmacology , Organophosphorus Compounds/chemistry , Organophosphorus Compounds/pharmacology , Acetylcholinesterase/drug effects , Acetylcholinesterase/metabolism , Animals , Cattle , Cholinesterase Inhibitors/chemical synthesis , Electrophorus , Horses , Humans , Molecular Docking Simulation , Naphthalenes/chemical synthesis , Organophosphates/chemical synthesis , Organophosphorus Compounds/chemical synthesis , Structure-Activity RelationshipABSTRACT
All-atom molecular dynamics simulations of butyrylcholinesterase (BChE) sans inhibitor and in complex with each of 15 dialkyl phenyl phosphate derivatives were conducted to characterize inhibitor binding modes and strengths. Each system was sampled on the 250 ns timescale in explicit ionic solvent, for a total of over 4 µs of simulation time. A K-means algorithm was used to cluster the resulting structures into distinct binding modes, which were further characterized based on atomic-level contacts between inhibitor chemical groups and active site residues. Comparison of experimentally observed inhibition constants (KI) with the resulting contact tables provides structural explanations for relative binding coefficients and highlights several notable interaction motifs. These include ubiquitous contact between glycines in the oxyanion hole and the inhibitor phosphate group; a sterically driven binding preference for positional isomers that extend aromaticity; a stereochemical binding preference for choline-containing inhibitors, which mimic natural BChE substrates; and the mechanically induced opening of the omega loop region to fully expose the active site gorge in the presence of choline-containing inhibitors. Taken together, these observations can greatly inform future design of BChE inhibitors, and the approach reported herein is generalizable to other enzyme-inhibitor systems and similar complexes that depend on non-covalent molecular recognition.Communicated by Ramaswamy H. Sarma.
Subject(s)
Butyrylcholinesterase/chemistry , Cholinesterase Inhibitors/chemistry , Molecular Docking Simulation , Molecular Dynamics Simulation , Binding Sites , Catalytic Domain , Cholinesterase Inhibitors/pharmacology , Humans , Ligands , Molecular Conformation , Molecular Structure , Protein Binding , Structure-Activity RelationshipABSTRACT
A series of dialkyl phenyl phosphates (DAPPs) were synthesized and evaluated in silico and in vitro for inhibitory activity against acetylcholinesterase and butyrylcholinesterase. Among the compounds examined, several DAPPs were shown to be potent inhibitors of butyrylcholinesterase, while having little activity against acetylcholinesterase. The most potent and selective inhibitors were di-n-butyl phenyl phosphate (K(i)=43 microM), di-n-pentyl phenyl phosphate (K(i)=6 microM), and di-cyclohexyl phenyl phosphate (K(i)=7 microM), the first which was shown to be a competitive inhibitor while the latter two being partial competitive inhibitors. Flexible docking simulations suggested that relative binding affinities generally increased as a function of alkyl chain length, while the strength and nature of inhibitory activity depended on whether the compound bound deeply or midway in the active site gorge, or in the proposed peripheral site.