Rational design based on bioactive conformation analysis of pyrimidinylbenzoates as acetohydroxyacid synthase inhibitors by integrating molecular docking, CoMFA, CoMSIA, and DFT calculations

Pyrimidinylthiobenzoates constitute an important kind of herbicides targeting acetohydroxyacid synthase (AHAS, EC 2.2.1.6), which catalyze the first common step in branched-chain amino acid biosynthesis. Due to the symmetry of 4,6-dimethoxypyrimidyl, there are two kinds of conformation of pyrimidinylthiobenzoates: one's phenyl is left-extending (named conformation-L); the other's phenyl is right-extending (named conformation-R). On the basis of the assumption that 3D quantitative structure - activity relationship (QSAR) models derived from the bioactive conformation should give the best result, a strategy of density-functional-theory-based 3D-QSAR was proposed to identify the bioactive conformation of pyrimidinylthiobenzoates by integrating the techniques of molecular docking, comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), and density functional theory calculation. The combination of three criteria of q(2), r(2), and r(pred)(2) obtained from CoMFA and CoMSIA analyses clearly indicated that conformation-R rather than conformation-L might be the bioactive conformation for pyrimidinylthiobenzoates. A further comparison between the two binding modes indicated that pyrimidinylthiobenzoates and sulfonylureas have very similar binding sites, such as Trp586, Arg380, and Pro192. However, Lys251 formed H bonds with sulfonylureas rather than pyrimidinylthiobenzoates. In addition, the orientation of phenyl groups of the two classes of compounds in the binding pocket were revealed to be opposite, which explained why the mutation of Pro192 displayed different sensitivity to sulfonylureas and pyrimidinylthiobenzoates. On the basis of the understanding of interactions between pyrimidinyl-thiobenzoates and AHAS, we designed and synthesized six 8-(4,6-dimethoxypyrimidin-2-yloxy)-4-methylphthalazin-1-one derivatives according to the 3D-QSAR models. The excellent correlation between the tested K-i values against wildtype A. thaliana acetohydroxyacid synthase and the predicted IC50 values demonstrated the high reliability of the established 3D-QSAR models. To our knowledge, this is the first report highlighting the binding mode of herbicidal pyrimidinylthiobenzoates, which consisted of the reported results of herbicide resistance.