sciences In Silico Study and Excito-Repellent Activity of Vitex negundo L. Essential Oil against Anopheles gambiae

: (1) Background: Essential oil from Vitex negundo is known to have repellent and insecticidal properties toward the Anopheles gambiae and this is linked to its monoterpene and sesquiterpene content. In this work, an effort is made to delineate the constitution of V. negundo essential oil (VNEO) and their interaction with odorant-binding proteins (OBPs) of A. gambiae and hence access its repellent efﬁciency as cost-effective and safer malaria vector control alternatives. (2) Methods: Anopheles species authentication was performed by genomic DNA analysis and was subjected to behavioral analysis. GC-MS proﬁling was used to identify individual components of VNEO. Anopheles OBPs were obtained from the RCSB protein data bank and used for docking studies. Determination of ligand efﬁciency metrics and QSAR studies were performed using Hyper Chem Professional 8.0.3, and molecular dynamics simulations were performed using the Desmond module. (3) Results: GC-MS analysis of VNEO showed 28 compounds (monoterpenes, 80.16%; sesquiterpenes, 7.63%; and unknown constituents, 10.88%). The ligand efﬁciency metrics of all four ligands against the OBP 7 were within acceptable ranges. β -selinene ( − 12.2 kcal/mol), β -caryophellene ( − 9.5 kcal/mol), sulcatone ( − 10.9 kcal/mol), and α -ylangene ( − 9.3 kcal/mol) showed the strongest binding afﬁnities for the target proteins. The most stable hydrophobic interactions were observed between β -selinene (Phe111 and Phe120), Sulcatone (Phe54 and Phe120), and α -ylangene (Phe111), while only sulcatone (Tyr49) presented H-bond interactions in the simulated environment. (4) Conclusions: Sulcatone and β -caryophyllene presented the best log p values, 6.45 and 5.20, respectively. These lead phytocompounds can be used in their purest as repellent supplement or as a natural anti-mosquito agent in product formulations.


Introduction
In hot and humid regions of the world, malaria is a huge health burden. It is an acute fever sickness caused by Plasmodium parasites that are transmitted to humans by the bites of infected mosquitoes [1]. Seven mosquito species from the genus Anopheles, teins [17]. Odorant-binding proteins (OBPs) are water-soluble, small extracellular proteins (15)(16)(17) found in the sensillum lymph of the insect sensilla [18]. They are thought to bind the odorant molecule and transport it to the receptors on chemosensory neurons via the aqueous sensillar lymph. OBPs, which are often expressed outside the chemosensory organs, are responsible exclusively as carriers of chemicals that, once solubilized, were transported to the olfactory receptors [19,20]. While it is undeniable that OBPs can bind odorants, their mechanism of interaction with odorant compounds is less clear [21]. OBP 1, 7, and 4 are members of the odorant-binding protein family and orthologs involved in host recognition signals and repellent alarm pheromones in Aedes aegypti [22]. OBPs are known targets for plant-based anti-mosquito insecticides [23]. One of those plants is Vitex negundo L. (Verbenaceae), which has been established to have repellent, insecticidal, and mosquitocidal properties linked strongly to the monoterpenes and sesquiterpenes in its content [24,25]. Progress in the identification of new targets (particularly in mosquitoes) for insecticides has been hampered over the past decade by a lack of comprehensive data on OBP mechanisms of action [26]. However, due to recent developments, availability of more current data, and scientific proofs, the structural data of key proteins implicated in OBPs action pathways may now be accessed and used.
Computer-aided molecular design is a rational strategy that combines a variety of theoretical and computational approaches routinely used in current drug research [27]. It has evolved into a powerful technique for lead screening, lead optimization, and the development of novel drug molecules, including repellant compounds [28,29]. When compared to traditional experimental approaches, the aim is to speed up the process while simultaneously lowering cost and allowing for the capacity to deal with big, diversified databases of ligands. The goal of this study is to evaluate the potential of anopheles mosquito OBPs as targets for V. negundo-based anti-mosquito repellent molecules in silico and in vitro using X-ray crystallographic structural data of mosquito OBPs in the complex.

Collection and Identification of V. negundo Leaves
In August 2021, the leaves of V. negundo were harvested in Benue State, Nigeria. The leaves were identified, and voucher specimen number NARICT/BSH/1610 was deposited at the National Research Institute for Chemical Technology (NARICT) in Zaria, Nigeria.

Leave Handling and Extraction
Fresh leaves of V. negundo were collected, washed, and extracted separately within 12 h of collection. Extraction was performed according to the procedure of Okoli et al. [30] with a 25 kg capacity fabricated essential oil distillation system (EDS). The plant was extracted for 45 min, and the distillates were recovered and separated into essential oil and hydrosol using a 2 L separatory funnel. The essential oil was then dried over anhydrous Na 2 SO 4 before being stored for further analysis.

GC-MS Profiling of the Essential Oils
GC-MS analysis of V. negundo essential oil was performed using a Varian CP3800 gas chromatograph with an HP5 capillary column (30 mm 0.25 mm, thickness 0.25 m), a flow rate of 1.2 mL/min, and gas carrier ion nitrogen grains. Varian Saturn 2000 Trap Detector Temperatures in the oven were gradually increased from 50 to 280 • C at a rate of 3 • C/min. Analysis conditions: injection temperature and transmission line temperatures of 220 and 240 • C, respectively. Injection volume: 0.2 L of 10% hexane solution, 1:30 ratio. The co-production of essential oils with solutions comprising C8-C22 alkanes of the same class results in a precise maintenance index for all molecules. Individual components are identified using storage indicators, which are then compared to previously described chemicals in the literature [31,32]. In addition, a computer library linked to GC

Preparation of Targets
Uniprot (http://www.uniprot.org/) (accessed on 11 May 2022) was used to compile all of the protein information. OBP 1 (PDB ID 3N7H), OBP 7 (PDB ID 3R1O), OBP 4 (PDB ID 3Q8I), and OBP 5 (PDB ID 3Q8I) are four A. gambiae OBP targets that were retrieved from the RCSB protein data bank (http://www.rcsb.org) (accessed on 11 May 2022) (PDB ID 2ERB). To view the protein and eliminate all water molecules from the structure, we used the PYMOL (version 1.7.4.5 program, Madison, WI, USA). The Swiss PDB (viewer v4.1.0, Basel, Switzerland) was used to repair void atomic spaces and crystallographic disturbances using energy reduction. The ".pdb" format was used to preserve the improved protein structure.

Ligands Preparation
In this study, all phytocompounds identified from the GC-MS analysis in the essential oils were retrieved from the PubChem database (www.pubchem.ncbi.nlm.nih.gov) (accessed on 21 March 2022) and saved in the Structure Data Format (SDF). Further, N,Ndiethyl-3-methylbenzamide (DEET) (PubChem CID 4284) was selected as a positive control because it is widely used as a chemical repellent against a variety of insects [33,34], strong electrophysiological responses [35], and interference with the olfactory receptor neurons/receptors resulting in avoidance behavior [36].

QSAR Studies
The physico-chemical properties of the FLP in correlation to biological activity were studied using QSAR modeling on the Hyper Chem (Professional 8.0.3 program, New Delh, India). The structure of each of the FLP was optimized with a semi-empirical PM3 method using a (MM+) force field, while a Fletcher-Reeves conjugate gradient algorithm was used for the energy minimization. A range of QSAR parameters, including hydration energy, free energy, refractivity, total energy, surface area, polarizability, volume, dipole moment, RMS gradient, mass, and partition coefficient (log P), were calculated.

Molecular Docking and Binding Site Prediction
In silico molecular docking with the aid of AutoDock vina (version 4; The Scripps Research Institute, La Jolla, CA, USA) was used to find the suitable binding orientations and conformations of the ligand with the targeted proteins. This is necessary for calculating the ligand s preferred orientations with the highest binding affinities for the protein s active sites, which are connected with structural pockets and cavities. The drug s binding to specific amino acid residues was displayed in the BIOVIA Discovery Studio visualizer (v16.1.0.15350, Paris, France) after molecular docking.

Molecular Dynamics Simulation
The stability of the lead phytochemicals from the docking analysis complexed with the odorant-binding protein 7 was observed over 100 ns molecular dynamics simulation MDS studies using the Desmond module of Schrodinger 2019-4 in the explicit solvent system with an OPLS3 force field.
A 10 Å buffer region of orthorhombic periodic boundary conditions was used with a solvated system of TIP3P molecules. The system was subsequently neutralized by counter ions (Na+) in an ensemble (NPT) of Nose-Hoover thermostats with a barostat to maintain a consistent temperature (310 K) and constant pressure (1 bar) for all the systems [37]. The hybrid energy minimization strategy, which was followed by conjugate gradient techniques, was created with the steepest descent of over 1000 steps. For energy reduction, the Broyden-Fletcher-Goldfarb-Shanno (LBFGS) method with a convergence threshold gradient of 1 Kcal/mol/Å was also used.
Appl. Sci. 2022, 12, 7500 5 of 22 The smooth particle mesh Ewald technique was used to calculate electrostatic interactions within a cut-off radius of 9 Å for short-range van der Waals and Coulomb interactions. The reference system propagator algorithms (RESPA) integration for multiple time steps was used in the dynamics study of bonded (2 fs), near (2 fs), and far-bonded (6 fs) interactions. The data were collected for every 100 ps, and the obtained trajectory was analyzed with Maestro graphical interphase. Various structural and thermodynamic parameters were computed from the trajectory files of the individual systems.
2.9. Rearing, Identification, Genomic DNA Extraction, and PCR Amplification of Anopheles Species Collection of mosquito larvae was performed from a 0.12 m × 2.5 m deep temporary pool with grass vegetation in Kaduna State, Nigeria, and identification of larvae was performed at the sites of breeding, along with morphological classification. The larvae were then transported to the Biological Sciences laboratory insectary of Kaduna State University, Nigeria. They were then sorted into the Anophelinae subfamilies using the Coetzee, [38] identification key under a compound microscope. After a 2 h acclimatization period, the immature larvae were fed a low-fat flour-baked product and then transferred to breeding chambers. There, they were reared to adulthood in a separate square wooden chamber for three weeks at 25 • C, 65% relative humidity, and a regulated light-dark (14/10 h) cycle.
Using a basic Olympus light microscope, the emerged adults were identified morphologically using taxonomic characters such as palps, proboscis, wing venation, markings, or tuffs on legs or abdomen as provided by the dichotomous keys used by Coetzee [38]. The adults in the cages were fed a 10% sucrose solution after eclosion from their pupal cases and allowed to rest and grow for 2 to 3 days. Only newly emerged adult female A. gambiae were manually aspirated into a 200 mL perforated plastic container and allowed to rest for one hour before being exposed to the essential oil mixture and pure components ( Figure 1).

Mosquito Behavioral Study
The mosquito behavioral responses to the V. negundo essential oil were conducted at an optimal repellency concentration of 0.48% v/v, determined in a previous study by Okoli et al. [30].; while the pure constituents (analytical purity, Sigma-Aldrich (Johannesburg, South Africa) were selected based on their binding energies. A 1.5 g cotton ball was im- Genomic DNA (gDNA) was extracted from 20 adult Anopheles mosquitoes using the Quick-DNATM Miniprep Plus Kit (D4069) product by ZYMO Research Company according to the manufacturer s protocol, as described by Okoli et al. [30]. Adult mosquitoes from the Aedes gambiae (s.l) complex were tested using PCR and genomic DNA assays for species identification, molecular form identification, and molecular analysis. Ahmadu Bello University, Zaria, Nigeria, conducted the molecular analyses.

Mosquito Behavioral Study
The mosquito behavioral responses to the V. negundo essential oil were conducted at an optimal repellency concentration of 0.48% v/v, determined in a previous study by Okoli et al. [30]; while the pure constituents (analytical purity, Sigma-Aldrich (Johannesburg, South Africa) were selected based on their binding energies. A 1.5 g cotton ball was immersed into the solution of essential oil and selected compounds in petrolatum, which was then mounted on the cotton ball holder of the test arm of the olfactometer [39]. Pure petrolatum was used as a negative control, and the air stream within the chamber was controlled at a rate of 180 mL/min using an electric air pump.
Behavioral responses were monitored and recorded for 1 h in a glass olfactometer containing precisely 150 female A. gambiae mosquitoes [40]. The olfactometer was air cleaned with a stream of hot air (>60 • C) after each investigation, and the cotton ball was removed, and the holder cleaned. The olfactory test was carried out three times. The repellent rate was calculated using the equation below (1) [41]: % Repelled = 100 − mean number of mosquitoes selecting essential oil/pure compound 100 − mean number of mosquitoes not selecting essential oil/pure compound 100 (1) The Probit analysis model in IBM SPSS v.25 statistical software (Chicago, IL, USA) was used to estimate the 50% mosquito repelled rate.

PCR Confirmation of Anopheles Gambiae s.s
The PCR amplicons examined under the trans-illuminator UV light revealed positive bandwidths of 390 bp for A. gambiae and 315 bp for A. arabiensis (Plate S1). Positive bandwidths of 390 bp were seen in the PCR amplicons when seen under the trans-illuminator UV light. After conditioning the PCR, samples 1-7 of randomly selected A. gambiae s.l displayed DNA bandwidths of 475 bp, confirming the species to be A. gambiae s.s (Plate S2). Table 1 shows the results of the V. negundo essential oils GC-MS analysis, and the chromatograms are presented in Figure S1. The essential oil contained 28 known compounds. Monoterpene and sesquiterpene contents are 80.16% and 7.63%, respectively, while other unknown constitutes account for about 10.88%. The observed compounds are in consonance with the reports of Gill et al. [42] and Huang et al. [43].

Molecular Docking
The binding energies of the GC-MS-identified phytocompounds from essential oils of V. negundo to four A. gambiae OBPs proteins are reported in Table S1. Four phytocompounds that demonstrated a multiplicity of binding properties to the four proteins with varying degrees of interaction within the active pockets of the proteins were selected as the lead phy-tocompounds. In Table S1, β-caryophellene (−9.5 kcal/mol), β-selinene (−12.2 kcal/mol), sulcatone (−10.9 kcal/mol), and α-ylangene (−9.3 kcal/mol) demonstrated the lowest binding energy for OBP7 compared to OBP1, OBP4, and OBP. β-selinene had the highest binding energies (−11.1, −11.2, and −12.2 Kcla/mol) to three of the OBPs (OBP 1, OBP 4, and OBP, respectively), thereby exhibiting a multiplicity of binding tendencies. β-selinene, in part or in association with the four lead phytochemicals, may be involved in the functional blocking of the olfactory receptor co-receptors, which underlines the activities reported. The four lead phytocompounds from the docking studies were selected for further computational investigation.

Amino Acid Interaction of Lead Phytocompounds with Selected A. gambiae OBP
The four lead phytocompounds demonstrated the highest binding energy to the OBP 7 (3R1O) than the other 3 OBPs; hence, it was selected for further interactive analysis with Appl. Sci. 2022, 12, 7500 8 of 22 the lead phytocompounds. The interaction of the FLP with amino acid residues of the binding pocket of OBP 7 (3R1O) is presented in Table 2. Table 2. Amino acid interaction of lead phytocompounds with A. gambiae OBP 7 (3R1O).

S/No
Compounds Binding Energy (Kcal) Interacting Amino Acid Residues It was observed that all FLP interacted with amino acid residues Phe54, Phe111, Phe120, and Pro41 majorly through hydrophobic interactions ( Figure 2). Previous experimental binding studies have established that the binding grooves of the A. gambiae OBPs are elongated hydrophobic groves that are capable of accommodating both linear and polycyclic compounds. The FLP interacted with binding site residues of known inhibitory compounds, including the AZO and palmitic acid binding sites. The properties of the interacting residues were responsible for the multiple hydrophobic interactions. It was also observed that the FLP exhibited a similar orientation in the hydrophobic pocket of the selected A. gambiae OBP. Though the pocket is elongated, the FLP is accommodated in the same region of the binding pocket ( Figure 3). Previous experimental binding studies have established that the binding grooves of the A. gambiae OBPs are elongated hydrophobic groves that are capable of accommodating both linear and polycyclic compounds. The FLP interacted with binding site residues of known inhibitory compounds, including the AZO and palmitic acid binding sites. The properties of the interacting residues were responsible for the multiple hydrophobic interactions. It was also observed that the FLP exhibited a similar orientation in the hydrophobic pocket of the selected A. gambiae OBP. Though the pocket is elongated, the FLP is accommodated in the same region of the binding pocket ( Figure 3).  Table 3 shows the ligand efficiency metrics of the selected ligands. Ligand efficiency (LE), ligand lipophilic efficiency (LLE), and fit quality (FQ) are expected to have threshold values of 0.3, 3, and 0.8 for a molecule to be classified as a potential hit quantitatively [48]. During lead discovery, the ligand efficiency lipophilic price (LELP) is estimated to be between −10 and 10 [49]. The ligand efficiency metrics of the four ligands against the OBP are within the criteria, qualifying them as a possible odorant-binding protein repellent compound.

QSAR Studies
The physico-chemical properties of the FLP in correlation to biological activity as computed from their QSAR parameters are presented in Table 4.   Table 3 shows the ligand efficiency metrics of the selected ligands. Ligand efficiency (LE), ligand lipophilic efficiency (LLE), and fit quality (FQ) are expected to have threshold values of 0.3, 3, and 0.8 for a molecule to be classified as a potential hit quantitatively [48]. During lead discovery, the ligand efficiency lipophilic price (LELP) is estimated to be between −10 and 10 [49]. The ligand efficiency metrics of the four ligands against the OBP are within the criteria, qualifying them as a possible odorant-binding protein repellent compound.

QSAR Studies
The physico-chemical properties of the FLP in correlation to biological activity as computed from their QSAR parameters are presented in Table 4. Sulcatone and beta-selinene presented the best QSAR parameters. For instance, both compounds exhibited the highest partition coefficient log p-value, which is important in evaluating the rate of permeability into cell membranes [50].

Frontier Molecular Orbitals (FMO) Studies
The chemical reactivity and the most likely reactive sites of the FLP were studied using FMO analysis (Figure 4). Sulcatone and beta-selinene presented the best QSAR parameters. For instance, bot compounds exhibited the highest partition coefficient log p-value, which is important i evaluating the rate of permeability into cell membranes [50].

Frontier Molecular Orbitals (FMO) Studies
The chemical reactivity and the most likely reactive sites of the FLP were studie using FMO analysis (Figure 4). The calculated energies, EHOMO and ELUMO, of compounds help to explain th global reactivity descriptors such as chemical potential, chemical hardness, and electro philicity of the FLP were elucidated from the computed EHOMO and ELUMO, while th negative values obtained for the EHOMO and ELUMO were confirmatory of the stabilit of the FLP. It was observed that β-caryophyllene had the least band energy gap (EHOMO The calculated energies, EHOMO and ELUMO, of compounds help to explain the global reactivity descriptors such as chemical potential, chemical hardness, and electrophilicity of the FLP were elucidated from the computed EHOMO and ELUMO, while the negative values obtained for the EHOMO and ELUMO were confirmatory of the stability of the FLP. It was observed that β-caryophyllene had the least band energy gap (EHOMO-ELUMO) compared to other phytocompounds. Therefore, β-caryophyllene demonstrated the highest chemical reactivity. The parameters that measure the chemical reactivity of the FLP, such as chemical potential (µ), electronegativity (χ), global hardness (η), global electrophilicity index (ω), and global softness (S), were calculated (Table S2). Sulcatone had the least chemical softness (S) hence, predicted to have higher stability than the other FLP. The electrophilicity (χ) is another important parameter that measures the ability of a compound to accept electron(s) from its surrounding. B-caryophyllene showed less electrophilicity than other compounds.

Interpretation of Molecular Dynamic Analysis
A full all-atom 100 ns molecular dynamics simulation was performed on the complexes of the lead phytochemicals (β-selinene, β-caryophellene, sulcatone, and α-ylangene) with the odorant-binding protein 7. The obtained MD simulation trajectories for the four systems were further analyzed using different structural and thermodynamic parameters.

Ligand Properties
The lead phytochemicals were analyzed relative to the reference conformation in the protein. Various structural parameters were plotted for each of the systems, including the root mean square deviation (RMSD), radius of gyration (rGyr), intramolecular hydrogen bonds (intraHB), molecular surface area (MolSA) within 1.4 Å probe radius, solvent accessible surface area (SASA), and polar surface area (PSA) (Figures S2-S5). After equilibration at the beginning of the simulation, a stable RMSD was observed for the β-caryophellene and sulcatone systems, while minimal RMSD fluctuation (<3.0 Å) was observed for the βselinene and α-ylangene systems, thus indicative of no immense dynamical alteration [51]. The rGyr profiles for the four systems were stable throughout the simulation period, suggesting no alterations in the chemical structure of the ligands [52]. Besides sulcatone, the three other phytochemical systems displayed no intramolecular hydrogen bond (intraHB) plots. The MolSA and SASA for the four phytochemicals in the referenced systems were in acceptable ranges, indicating that they were imposingly stable in the complexes during the simulation run.

Protein Secondary Structure
The secondary structure elements (SSE) of the odorant-binding protein 7 upon the binding of the lead phytochemicals such as alpha-helices and beta-strands were analyzed during the simulation period. Figure 5 shows the SSE distribution by residue index throughout the protein structure. The total SSE (%) for the protein upon the binding of β-selinene, β-caryophellene, sulcatone, and α-ylangene is 55.38, 54.51, 54.12, and 47.83, while the alpha-helices (%) and beta-strands (%) were, respectively, 53.35 and 2.02; 54.36 and 0.14; 53.09 and 1.03; 47.74 and 0.09 for β-selinene, β-caryophellene, sulcatone, and α-ylangene bound systems ( Figure 5A-D). The analysis of the four systems shows that binding of the α-ylangene to the protein minimally reduced the SSE of the odorant protein, hence the conformational integrity of the protein ( Figure 5D).

Protein Root Mean Square Fluctuation
The RMSF residual index is valuable for depicting local changes along the protein chain during the simulation. Figure 6A   The RMSF residual index is valuable for depicting local changes along the protein chain during the simulation. Figure 6A-D shows the RMSF plots of odorant-binding protein 7 complexed with the lead phytochemicals. The peaks indicate areas of the protein with the highest fluctuations during the course of the simulation. Besides the fluctuation at the terminals of the plots, peaks were observed between the 40 and 80 ns time frames for all four systems. This corresponds to the unstructured region of the proteins. The structured regions (alpha-helical and betastrand) presented less fluctuation than the unstructured region (usually the loop region) of the protein. The high peaks did not necessarily reflect a compromise in the integrity of the protein structure but motion due to the absence of a structured unit at that time frame.

Protein-Ligand Contacts
The stacked bar interaction plots that are categorized by type (hydrogen bonds, hydrophobic, ionic, and water bridges) of interactions were used to study the protein-lead phytochemical interactions or contact throughout the simulation period (Figure 7). Using the distance of 2.5 Å between the donor and acceptor atoms as the geometric criteria for H-bond, only the sulcatone bound system presented several H-bond interactions in the simulated environment. Among other H-bonds, the highest H-bond interaction fraction with about 92% contact time of the whole simulation time was between sulcatone and residue Tyr49. The peaks indicate areas of the protein with the highest fluctuations during the course of the simulation. Besides the fluctuation at the terminals of the plots, peaks were observed between the 40 and 80 ns time frames for all four systems. This corresponds to the unstructured region of the proteins. The structured regions (alpha-helical and beta-strand) presented less fluctuation than the unstructured region (usually the loop region) of the protein. The high peaks did not necessarily reflect a compromise in the integrity of the protein structure but motion due to the absence of a structured unit at that time frame.

Protein-Ligand Contacts
The stacked bar interaction plots that are categorized by type (hydrogen bonds, hydrophobic, ionic, and water bridges) of interactions were used to study the protein-lead phytochemical interactions or contact throughout the simulation period (Figure 7). Using the distance of 2.5 Å between the donor and acceptor atoms as the geometric criteria for H-bond, only the sulcatone bound system presented several H-bond interactions in the simulated environment. Among other H-bonds, the highest H-bond interaction fraction with about 92% contact time of the whole simulation time was between sulcatone and residue Tyr49.
Sulcatone also made several salt bridges with various residues of the odorant protein. The presented hydrogen-bonded protein-ligand interactions that were mediated by water molecules in the plot were selected using the geometric criteria for a protein-water or waterligand H-bond of 2.8 Å. Figure 7B shows that Glu116, Lys119, and Arg9 made the highest salt bridges with the ligand groups of sulcatone. Several hydrophobic contacts, including pi-cation, pi-pi, and other non-specific interactions, were observed with different residues of the odorant proteins. Among the several hydrophobic interactions, Phe111 and Phe120 of the odorant protein maintained the highest hydrophobic contact of more than 45% of the simulation time with β-selinene. Sulcatone maintained the highest hydrophobic contact with Phe54 and Phe120 for more than 40% of the simulation period. β-caryophellene presented several contacts of less than 20% of the simulation with the odorant protein, while the hydrophobic contact between α-ylangene and residue Phe111 was maintained for more than 40% of the simulation time. An in-depth analysis plot of the various types of interactions that occurred over the simulation period of 100 ns was further presented in a timeline representation in Figure 8. Sulcatone also made several salt bridges with various residues of the odorant protein. The presented hydrogen-bonded protein-ligand interactions that were mediated by water molecules in the plot were selected using the geometric criteria for a protein-water or water-ligand H-bond of 2.8 Å. Figure 7B shows that Glu116, Lys119, and Arg9 made the highest salt bridges with the ligand groups of sulcatone. Several hydrophobic contacts, including pi-cation, pi-pi, and other non-specific interactions, were observed with different residues of the odorant proteins. Among the several hydrophobic interactions, Phe111 and Phe120 of the odorant protein maintained the highest hydrophobic contact of more than 45% of the simulation time with β-selinene. Sulcatone maintained the highest hydrophobic contact with Phe54 and Phe120 for more than 40% of the simulation period. βcaryophellene presented several contacts of less than 20% of the simulation with the odorant protein, while the hydrophobic contact between α-ylangene and residue Phe111 was maintained for more than 40% of the simulation time. An in-depth analysis plot of the various types of interactions that occurred over the simulation period of 100 ns was further presented in a timeline representation in Figure 8. The top panel shows the detailed contacts the odorant protein made with the lead phytochemicals over the course of the trajectory, while the bottom panel shows the residue in contact with the ligand in each of the trajectory frames.

Root Mean Square Deviation Analysis
The root mean square deviation (RMSD) measures the alteration in the displacement of atoms for a particular frame with respect to a reference frame for all the trajectories. Figure 9A-D shows the RMSD of the odorant-binding protein 7 in complex with the lead phytochemicals. The red line (left Y-axis) in the plots shows the RMSD evolution of the protein after all protein frames are first aligned on the reference frame backbone. The ligand RMSD (right Y-axis) shows how stable the lead phytochemicals are with respect to the protein. The protein-ligand complex was first aligned on the reference protein backbone. Figure 9A shows the plots of the RMSD of the β-selinene complex system.

Root Mean Square Deviation Analysis
The root mean square deviation (RMSD) measures the alteration in the displacement of atoms for a particular frame with respect to a reference frame for all the trajectories. Figure 9A-D shows the RMSD of the odorant-binding protein 7 in complex with the lead phytochemicals. The red line (left Y-axis) in the plots shows the RMSD evolution of the protein after all protein frames are first aligned on the reference frame backbone. The ligand RMSD (right Y-axis) shows how stable the lead phytochemicals are with respect to the protein. The protein-ligand complex was first aligned on the reference protein backbone. Figure 9A shows the plots of the RMSD of the β-selinene complex system.
The protein RMSD was equilibrated before the first 10 ns, while the ligand RMSD was equilibrated after the first 20 ns. Both systems were stabilized with minimal fluctuations and converged at 100 ns. Both the protein and ligand RMSD for the sulcatone complexed system were equilibrated before 10 ns of the simulation with minimal fluctuation (less than 2 Å) after equilibration. The protein RMSD of the β-caryophellene bound system was equilibrated after 20 ns with a sharp peak after equilibration, while the ligand RMSD The protein RMSD was equilibrated before the first 10 ns, while the ligand RMSD was equilibrated after the first 20 ns. Both systems were stabilized with minimal fluctuations and converged at 100 ns. Both the protein and ligand RMSD for the sulcatone complexed system were equilibrated before 10 ns of the simulation with minimal fluctuation (less than 2 Å) after equilibration. The protein RMSD of the β-caryophellene bound system was equilibrated after 20 ns with a sharp peak after equilibration, while the ligand RMSD was equilibrated at the beginning of the simulation. Both systems experienced minimal fluctuations that were within the acceptable range of <3 Å, and convergence was also reached toward 100 ns. The protein RMSD plots for the α-ylangene bound system followed a different pattern from the ligand RMSD. Both systems experienced minimal fluctuations that were within the acceptable range of <3 Å, and convergence was also reached toward 100 ns. The protein RMSD plots for the α-ylangene bound system followed a different pattern from the ligand RMSD. Both systems were equilibrated before 20 ns but experienced several deep and peak fluctuations throughout the simulation. The deviance in the RMSD plots further explains the reduction in the secondary structure elements (alpha-helices and beta-strands) ( Figure 5D). The binding of α-ylangene reduced the stability of the protein [53]. In general, the odorant protein complexed with the lead phytochemicals except α-ylangene demonstrated a high degree of structural stability and compactness [54]; hence, the systems can be adapted for other experimental models.
fluctuations that were within the acceptable range of <3 Å, and convergence was also reached toward 100 ns. The protein RMSD plots for the α-ylangene bound system followed a different pattern from the ligand RMSD. Both systems experienced minimal fluctuations that were within the acceptable range of <3 Å, and convergence was also reached toward 100 ns. The protein RMSD plots for the α-ylangene bound system followed a different pattern from the ligand RMSD. Both systems were equilibrated before 20 ns but experienced several deep and peak fluctuations throughout the simulation. The deviance in the RMSD plots further explains the reduction in the secondary structure elements (alpha-helices and beta-strands) ( Figure 5D). The binding of α-ylangene reduced the stability of the protein [53]. In general, the odorant protein complexed with the lead phytochemicals except α-ylangene demonstrated a high degree of structural stability and compactness [54]; hence, the systems can be adapted for other experimental models.  Figure 10 indicates five ligands identified in the V.negundo essential oil, and these results are reported in Table S1.  3.14. Response of A. gambiae to the Essential Oil and Selected Compounds Figure 10 indicates five ligands identified in the V.negundo essential oil, and these results are reported in Table S1. The observed repellency of the essential oil compared to the chemical-based repellent underlines the potency of the oil, which is in tandem with the ethno-application of the V. negundo as a mosquito repellent as reported by Hazarika et al. [55] and Zaidan et al. [56].

Response of A. gambiae to the Essential Oil and Selected Compounds
As the exposure time increases, the repellency activity of the compounds increases to a concentration where there are no observable changes in activity. The essential oil showed a significant increase in the percentage of mosquitoes repelled at 49.9 min compared to the compounds; however, no significant difference (p > 0.05) was observed in the percentage of mosquitoes repelled between the pure compounds.
As the exposure time increases, the repellency activity of the compounds increases to a concentration where there are no observable changes in activity. The essential oil showed a significant increase in the percentage of mosquitoes repelled at 49.9 min compared to the compounds; however, no significant difference (p > 0.05) was observed in the percentage of mosquitoes repelled between the pure compounds.
The contact time of the tested samples varied depending on their diffusivity ( Figures  3-5); in the first 10 min of exposure to V. negundo oil and DEET, they repelled about 76% and 85% of the A. gambiae (Figure 11a,b, while 24%, 18%, 14%, and 19% of the A. gambiae responded to the presence of β-selinene, β-caryophellene, sulcatone, and α-ylangene, respectively (Figure 11c- In 60 min, sulcatone, α-ylangene, β-caryophellene, and β-selinene, repelled 75%, 77%, 85%, and 83%, while the oil and DEET repelled 100% of the A. gambiae. There is a significant difference in the responses of the A. gambiae to the oil and DEET compared to the pure compounds (Figure 11c-f). The optimal repellencies of the oil and DEET were recorded at 29.8 min and 35.9 min. Pure compounds showed an increase in the amount of mosquitoes repelled with time, with an optimal % repellency activity attained at approximately 45.6-49.9 min.
However, the sensitivities of the A. gambiae to the test compounds vary slightly, with A. gambiae being most and least sensitive to sulcatone at 15 min and β-selinene at about 26.7 min. Further, the time required to repel 50% of the A. gambiae varies from one pure compound to another. Sulcatone and β-caryophellene showed 50% repellency activity in the shortest time of 23.9 min, contrary to the observed repellency of α-ylangene at 25.3 min. Sulcatone has been reported to be potent against Drosophila melanogaster and Toxorhynchites amboinensis [57]. This inference is in consonance with the studies of Angraeni et al. [58] and Nararak et al. [59,60]. Most insect repellents, including DEET, have a species-specific method of action that is sometimes unknown. Activation of specific odorant receptors (ORs), inhibition of specific ORs, and modulation of multiple ORs, according to Tsitoura et al. [34] and Degennaro et al. [61], may be suggestive of functional blocking of the olfactory receptor co-receptor, resulting in A. gambiae susceptibility to the molecules in the repellent.

Conclusions
The essential oil of V. negundo showed significant mosquito repellent efficacy (p < 0.05) when compared to DEET, with optimal repellency achieved in 29.8 min versus 35.7 min for DEET. According to the docking studies, four phytocompounds (β-selinene, βcaryophellene, sulcatone, and α-ylangene) showed decreased activity when tested individually, with an average optimal repellency time of 48.25 min. This decrease in activity is due to the loss of synergism. The ligand efficiency values were all within the expected ranges, implying that the ligands were quantitatively hit and thus qualify as a potential odorant-binding protein repellent lead. Among the lead phytochemicals, β-selinene displayed a multiplicity of highest binding tendencies to three of the OBPs (OBP 1, OBP 4, and OBP). β-selinene, in part or in association with the four lead phytochemicals, may be involved in the functional blocking of the olfactory receptors and co-receptors, which un- In 60 min, sulcatone, α-ylangene, β-caryophellene, and β-selinene, repelled 75%, 77%, 85%, and 83%, while the oil and DEET repelled 100% of the A. gambiae. There is a significant difference in the responses of the A. gambiae to the oil and DEET compared to the pure compounds (Figure 11c-f). The optimal repellencies of the oil and DEET were recorded at 29.8 min and 35.9 min. Pure compounds showed an increase in the amount of mosquitoes repelled with time, with an optimal % repellency activity attained at approximately 45.6-49.9 min.
However, the sensitivities of the A. gambiae to the test compounds vary slightly, with A. gambiae being most and least sensitive to sulcatone at 15 min and β-selinene at about 26.7 min. Further, the time required to repel 50% of the A. gambiae varies from one pure compound to another. Sulcatone and β-caryophellene showed 50% repellency activity in the shortest time of 23.9 min, contrary to the observed repellency of α-ylangene at 25.3 min. Sulcatone has been reported to be potent against Drosophila melanogaster and Toxorhynchites amboinensis [57]. This inference is in consonance with the studies of Angraeni et al. [58] and Nararak et al. [59,60].
Most insect repellents, including DEET, have a species-specific method of action that is sometimes unknown. Activation of specific odorant receptors (ORs), inhibition of specific ORs, and modulation of multiple ORs, according to Tsitoura et al. [34] and Degennaro et al. [61], may be suggestive of functional blocking of the olfactory receptor co-receptor, resulting in A. gambiae susceptibility to the molecules in the repellent.

Conclusions
The essential oil of V. negundo showed significant mosquito repellent efficacy (p < 0.05) when compared to DEET, with optimal repellency achieved in 29.8 min versus 35.7 min for DEET. According to the docking studies, four phytocompounds (β-selinene, β-caryophellene, sulcatone, and α-ylangene) showed decreased activity when tested individually, with an average optimal repellency time of 48.25 min. This decrease in activity is due to the loss of synergism. The ligand efficiency values were all within the expected ranges, implying that the ligands were quantitatively hit and thus qualify as a potential odorant-binding protein repellent lead. Among the lead phytochemicals, β-selinene displayed a multiplicity of highest binding tendencies to three of the OBPs (OBP 1, OBP 4, and OBP). β-selinene, in part or in association with the four lead phytochemicals, may be involved in the functional blocking of the olfactory receptors and co-receptors, which underlines the repellant activities of the essential oil from A. gambiae. The complexes formed between the lead phytocompounds and the OBP 7 were stable in the dynamic environment, hence can be adapted for further wet experimental studies. Based on these findings, V. negundo oil can be used in a variety of ways to enhance or replace repellent agents in traditional repellents. It can also be used as a natural anti-mosquito agent in the preparation of aerosol, repellent lotion, and repellent fabrics.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/app12157500/s1, Plate S1. Lane M is the 100bp marker, Lane 1-20 are randomly selected Anopheles samples. Lane 21 = negative Sample. Distinguishing band size: Anopheles gambiae s.l at 390bp; Anopheles arabiensis 315 bp; Plate S2. Agarose gel 1.5% for distinguishing An. gambiae s.s and coluzzi after PCR with primers (R3, R5, B/Sint & MoPint). Lane M = 100bp molecular weight marker; visible at 500 bp and 1000 bp. Samples 1 to 7 are Anopheles gambiae s.l randomly picked from a group of 20 in Kaduna and shows DNA band sizes of 475bp to authenticate species to be Anopheles gambiae s.s, Figure S1: Chromatogram of V. negundo essential oil, Figure S2: Ligand properties of β-selinene, Figure S3: Ligand properties of β-caryophellene, Figure  S4: Ligand properties of sulcatone, and Figure Table S1: Molecular docking results for the interaction between the ligands and the target proteins, and Table S2: Global reactivity descriptors of the FLP.