In Silico Design and Validation of a Novel HPPD-Inhibiting Herbicide Candidate Based on Benzofuran and Arylthioacetic...

Inhibition of 4-hydroxyphenylpyruvate dioxygenase (HPPD) is a well-established strategy for weed control, yet the emergence of resistance underscores the need for more potent inhibitors. In this study, datasets of benzofuran analogues and arylthioacetic acid-derived triketones were analyzed using multivariate image analysis of quantitative structure-activity relationships (MIA-QSARs) to guide the identification of promising candidates. Molecular docking was conducted on HPPD, including Co(II) substitution to explore metal-ligand interactions, and molecular dynamics simulations, parametrized with MCPB.py, evaluated the stability and interaction patterns of the complexes. Fourteen candidates were proposed, six of which exhibited higher predicted activity and improved performance relative to mesotrione. Among them, candidate P1 emerged as the most promising, reproducing key interactions of mesotrione while displaying lower binding energy and stable convergence during dynamics. These results demonstrate P1 as a novel HPPD inhibitor and highlight the utility of combining MIA-QSAR, docking, and MCPB.py-parametrized dynamics in rational metalloenzyme inhibitor design.