A Combined Computational Methodology for the Discovery of Hit Compounds with Putative Insect Repellency Properties

2D finger-printing and molecular docking studies identified potent mosquito repellents targeting odorant binding protein 1

Personal protection measures against the mosquitoes like the use of repellents constitute valuable tools in the effort to prevent the transmission of vector-borne diseases. Therefore, the discovery of novel repellent molecules which will be effective at lower concentrations and provide a longer duration of protection remains an urgent need. Mosquito Odorant-Binding Proteins (OBPs) involved in the initial steps of the olfactory signal transduction cascade have been recognized not only as passive carriers of odors and pheromones but also as the first molecular filter to discriminate semiochemicals, hence serving as molecular targets for the design of novel pest control agents. Among the three-dimensional structures of mosquito OBPs solved in the last decades, the OBP1 complexes with known repellents have been widely used as reference structures in docking analysis and molecular dynamics simulation studies for the structure-based discovery of new molecules with repellent activity. Herein, ten compounds known to be active against mosquitoes and/or displaying a binding affinity for Anopheles gambiae AgamOBP1 were used as queries in an in silico screening of over 96 million chemical samples in order to detect molecules with structural similarity. Further filtering of the acquired hits on the basis of toxicity, vapor pressure, and commercial availability resulted in 120 unique molecules that were subjected to molecular docking studies against OBP1. For seventeen potential OBP1-binders, the free energy of binding (FEB) and mode of interaction with the protein were further estimated by molecular docking simulations leading to the selection of eight molecules exhibiting the highest similarity with their parental compounds and favorable energy values. The in vitro determination of their binding affinity to AgamOBP1 and the evaluation of their repellent activity against female Aedes albopictus mosquitoes revealed that our combined ligand similarity screening and OBP1 structure-based molecular docking successfully detected three molecules with enhanced repellent properties. A novel DEET-like repellent with lower volatility (8.55 × 10^-4 mmHg) but a higher binding affinity for OBP1 than DEET (1.35 × 10^-3 mmHg). A highly active repellent molecule that is predicted to bind to the secondary Icaridin (sIC)-binding site of OBP1 with higher affinity than to the DEET-site and, therefore, represents a new scaffold to be exploited for the discovery of binders targeting multiple OBP sites. Finally, a third potent repellent exhibiting a high degree of volatility was found to be a strong DEET-site binder of OBP1 that could be used in slow-release formulations.

The structure of AgamOBP5 in complex with the natural insect repellents Carvacrol and Thymol: Crystallographic, fluorescence and thermodynamic binding studies

Among several proteins participating in the olfactory perception process of insects, Odorant Binding Proteins (OBPs) are today considered valid targets for the discovery of compounds that interfere with their host-detection behavior. The 3D structures of Anopheles gambiae mosquito AgamOBP1 in complex with the known synthetic repellents DEET and Icaridin have provided valuable information on the structural characteristics that govern their selective binding. However, no structure of a plant-derived repellent bound to an OBP has been available until now. Herein, we present the novel three-dimensional crystal structures of AgamOBP5 in complex with two natural phenolic monoterpenoid repellents, Carvacrol and Thymol, and the MPD molecule. Structural analysis revealed that both monoterpenoids occupy a binding site (Site-1) by adopting two alternative conformations. An additional Carvacrol was also bound to a secondary site (Site-2) near the central cavity entrance. A protein-ligand hydrogen-bond network supplemented by van der Waals interactions spans the entire binding cavity, bridging α4, α6, and α3 helices and stabilizing the overall structure. Fluorescence competition and Differential Scanning Calorimetry experiments verified the presence of two binding sites and the stabilization effect on AgamOBP5. While Carvacrol and Thymol bind to Site-1 with equal affinity in the submicromolar range, they exhibit a significantly lower and distinct binding capacity for Site-2 with Kd's of ~7 μΜ and ~18 μΜ, respectively. Finally, a comparison of AgamOBP5 complexes with the AgamOBP4-Indole structure revealed that variations of ligand-interacting aminoacids such as A109T, I72M, A112L, and A105T cause two structurally similar and homologous proteins to display different binding specificities.

Characterization in Effective Stimulation on the Magnitude, Gating, Frequency Dependence, and Hysteresis of INa Exerted by Picaridin (or Icaridin), a Known Insect Repellent

Picaridin (icaridin), a member of the piperidine chemical family, is a broad-spectrum arthropod repellent. Its actions have been largely thought to be due to its interaction with odorant receptor proteins. However, to our knowledge, to what extent the presence of picaridin can modify the magnitude, gating, and/or the strength of voltage-dependent hysteresis (Hys(V)) of plasmalemmal ionic currents, such as, voltage-gated Na+ current [INa], has not been entirely explored. In GH3 pituitary tumor cells, we demonstrated that with exposure to picaridin the transient (INa(T)) and late (INa(L)) components of voltage-gated Na+ current (INa) were differentially stimulated with effective EC50's of 32.7 and 2.8 μM, respectively. Upon cell exposure to it, the steady-state current versus voltage relationship INa(T) was shifted to more hyperpolarized potentials. Moreover, its presence caused a rightward shift in the midpoint for the steady-state inactivate curve of the current. The cumulative inhibition of INa(T) induced during repetitive stimuli became retarded during its exposure. The recovery time course from the INa block elicited, following the conditioning pulse stimulation, was satisfactorily fitted by two exponential processes. Moreover, the fast and slow time constants of recovery from the INa block by the same conditioning protocol were noticeably increased in the presence of picaridin. However, the fraction in fast or slow component of recovery time course was, respectively, increased or decreased with an increase in picaridin concentrations. The Hys(V)'s strength of persistent INa (INa(P)), responding to triangular ramp voltage, was also enhanced during cell exposure to picaridin. The magnitude of resurgent INa (INa(R)) was raised in its presence. Picaritin-induced increases of INa(P) or INa(R) intrinsically in GH3 cells could be attenuated by further addition of ranolazine. The predictions of molecular docking also disclosed that there are possible interactions of the picaridin molecule with the hNaV1.7 channel. Taken literally, the stimulation of INa exerted by the exposure to picaridin is expected to exert impacts on the functional activities residing in electrically excitable cells.