Quantitative Structure Activity Relationship Studies and Molecular Dynamics Simulations of 2-(Aryloxyacetyl)cyclohexane-1,3-Diones Derivatives as 4-Hydroxyphenylpyruvate Dioxygenase Inhibitors

4-Hydroxyphenylpyruvate dioxygenase (HPPD) is a significant enzyme in the biosynthesis of plastoquinone and tocopherol. Moreover, it is also a potential target to develop new herbicide. The technology of computer-aided drug design (CADD) is a useful tool in the efficient discovery of new HPPD inhibitors. Forty-three compounds with known activities were used to generate comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models based on common framework and molecular docking. The structural contribution to the activity was determined, which provided further information for the design of novel inhibitors. Molecular docking was used to explain the changes in activity caused by the binding mode between ligand and protein. The molecular dynamics (MD) results indicated that the electrostatic energy was the major driving force for ligand–protein interaction and the Phe403 made the greatest contribution to the binding. The present work has provided useful information for the rational design of novel HPPD inhibitors with improved activity.