Ensemble-based, high-throughput virtual screening of potential inhibitor targeting putative farnesol dehydrogenase of Metisa plana (Lepidoptera: Psychidae)

Abstract

Metisa plana (Lepidoptera: Psychidae) bagworm is a leaf-eater caterpillar ubiquitously found as a damaging pest in oil palm plantations, specifically in Malaysia. Various strategies have been implemented, including the usage of chemical insecticides. However, the main challenges include the development of insecticide resistance and its detrimental effects on the environment and non-target organisms. Therefore, a biorational insecticide is introduced by targeting the juvenile hormone (JH) biosynthetic pathway, which is mainly present in the insect and vital for the insect's growth, diapause, metamorphosis, and adult reproduction. This study aimed to investigate the potential inhibitor for the rate-limiting enzyme involved in the JH pathway known as farnesol dehydrogenase. A 255 amino acids sequence encoded for the putative M. plana farnesol dehydrogenase (MpFolDH) open reading frame had been identified and isolated. The three-dimensional structure of MpFolDH was predicted to have seven β- sheets with α-helices at both sides, showing typical characteristics for classical short-chain dehydrogenase and associated with oxidoreductase activity. Then, the ensemble-based virtual screening was conducted based on the ZINC20 database, in which 43 768 compounds that fulfilled pesticide-likeness criteria were screened by site-specific molecular docking. After a short molecular dynamics simulation (5 ns) was conducted towards 102 compounds, only the top 10 compounds based on their most favourable binding energy were selected for a more extended simulation (100 ns). Based on the protein-ligand stability, protein compactness, residues rigidity, binding interaction, binding energy throughout the 100 ns simulation, and physicochemical analysis, ZINC000408743205 was selected as a potential inhibitor for this enzyme. Amino acids decomposition analysis indicates Ile18, Ala95, Val198 and Val202 were the critical contributor residues for MpFolDH-inhibitors(s) complex.