Computational investigation of the molecular conformation-dependent binding mode of (E)-β-farnesene analogs with a heterocycle to aphid odorant-binding proteins.

Odorant-binding proteins (OBPs) play an important role as ligand-transfer filters in olfactory recognition in insects. (E)-β-farnesene (EBF) is the main component of the aphid alarm pheromone and could keep aphids away from crops to prevent damage. Computational insight into the molecular binding mode of EBF analogs containing a heterocycle based on the structure of Megoura viciae OBP 3 (MvicOBP3) was obtained by molecular docking and molecular dynamics simulations. The results showed that high affinity EBF analogs substituted with an aromatic ring present a unique binding conformation in the surface cavity of MvicOBP3. A long EBF chain was located inside the cavity and was surrounded by many hydrophobic residues, while the substituted aromatic ring was exposed to the outside due to limitations from the formation of multiple hydrogen bonds. However, the low activity EBF analogs displayed an exactly inverted binding pose, with EBF loaded on the external side of the protein cavity. The affinity of the recently synthesized EBF analogs containing a triazine ring was evaluated in silico based on the binding modes described above and in vitro through fluorescence competitive binding assay reported later. Compound N1 not only showed a similar binding conformation to that of the high affinity analogs but was also found to have a much higher docking score and binding affinity than the other analogs. In addition, the docking score results correlated well with the predicted logP values for these EBF analogs, suggesting highly hydrophobic interactions between the protein and ligand. These studies provide an in silico screening model for the binding affinity of EBF analogs in order to guide their rational design based on aphid OBPs.