Vitex negundo L. Essential Oil: Odorant Binding Protein Efficiency Using Molecular Docking Approach and Studies of the Mosquito Repellent

Simple Summary Malaria fever kills millions of people annually in the tropical and subtropical countries of Africa and Asia. Because there is no effective vaccine, malaria prevention is exclusively dependent on avoiding human-vector interaction. The interaction of Vitex negundo essential oil constituents with Anopheles gambiae Odorant Binding Proteins (OBP), as well as its compositional variation, repellent efficacy, and toxicity profile, are investigated in this work. The oils were subjected to GC-MS analysis, a mosquito behavioral test, OBP-ligand interactions, Anopheles species authentication, and toxicity profile. Docking protocol validation was achieved by redocking the co-crystallized ligands and root mean square deviation (RMSD) calculation. The oil yields and compositions are climate–soil dependent with ≈71.39% monoterpenes and ≈16.32% sesquiterpene. Optimal repellency is achieved at 15 min at ED50 0.08–0.48% v/v while the RMSD was estimated to be within 0.24–1.35 Å. Strong affinities, −6.4 to −5.4 kcal/mol, were demonstrated by α-pinene, citronellal, linalool, and myrcene for OBP1, OBP7, OBP4, and OBP. respectively. The hydrophobic interactions involve Leu17, Cys35, ALA52, Leu73, Leu76, Ala88, Met91, Lys93, Trp114, Phe123, and Leu124 receptors on α-helixes 1–7 within the binding cavities, and may block the olfactory receptors resulting in disorientation. α-pinene, linalool, and myrcene are safe and suitable for use in the development of green and innovative repellents because their ligand efficiency metrics, ADME/tox, and repellency screening are all within the threshold values. Abstract (1) Background: Malaria fever affects millions of people yearly in Africa and Asia’s tropical and subtropical areas. Because there is no effective vaccine, malaria prevention is solely dependent on avoiding human-vector interaction. (2) Aim: This study examines the interaction between the constituents of Vitex negundo essential oil and Anopheles gambiae Odorant Binding Proteins (OBP) as well as the compositional variation, repellent efficacy, and toxicity profile. (3) Methods: The oils were subjected to GC-MS and mosquito behavioral analysis. OBP–ligand interactions, Anopheles species authentication, and the toxicity profile were determined by molecular docking, PCR assay and in silico ADME/tox tool. Docking protocol validation was achieved by redocking the co-crystallized ligands into the protein binding pocket and root mean square deviation (RMSD) calculation. (4) Results: The oil yields and compositions are climate–soil dependent with ≈71.39% monoterpenes and ≈16.32% sesquiterpene. Optimal repellency is achieved at 15 min at ED50 0.08–0.48% v/v while the RMSD was estimated to be within 0.24–1.35 Å. Strong affinities were demonstrated by α-pinene (−6.4 kcal/mol), citronellal (−5.5 kcal/mol), linalool (−5.4 kcal/mol), and myrcene (−5.8 kcal/mol) for OBP1, OBP7, OBP4, and OBP; respectively. The hydrophobic interactions involve Leu17 (α-helix 1), Cys35 (α-helix 2), ALA52 (α-helix 3), Leu73, Leu76 (α-helix 4), Ala88, Met91, Lys93, Trp114 (α-helix 5), Phe123 (α-helix 6), and Leu124 (α-helix 7) receptors within the binding cavities, and may cause blocking of the olfactory receptors resulting in disorientation. (5) Conclusion: The ligand efficiency metrics, ADME/tox and repellency screening are within the threshold values; hence, α-pinene, linalool, and myrcene are safe and fit-to-use in the development of a green and novel repellent.


Introduction
Malaria has been and continues to be one of the world's leading causes of mortality. Despite the fact that the disease has been eradicated in the United States, it remains a global health threat. Malaria affects an estimated 300 to 500 million people worldwide each year, resulting in 1.5 to 2.7 million fatalities. Sub-Saharan Africa, on the other hand, accounts for about 90% of all malaria cases worldwide [1], followed by the World Health Organization (WHO) South-East Asia Region (3.4%) and Eastern Mediterranean Region (2.1%). Nearly 85% of the worldwide malaria burden was carried by 19 nations in Sub-Saharan Africa and India. Nigeria (25%), the Democratic Republic of the Congo (12%), Uganda (5%), Côte d'Ivoire (4%), Mozambique (4%), and Niger (4%) accounted for more than half of all malaria cases globally [2]. In addition, one to two million people suffering from malaria die from the disease, the majority of whom are children under the age of five and pregnant women [3]. Unfortunately, 90% of malaria deaths occur in Sub-Saharan Africa, where approximately 3000 people die each year. The Nigerian situation is consistent with the global trend, as malaria is the leading cause of death among pregnant women in the country, accounting for one out of every ten deaths [4,5].
Malaria is caused by a species of Plasmodium which is vectored by the adult female Anopheles mosquito that thrives in hot and humid climates around the world [6]. The Anopheles gambiae complex consists of at least seven morphologically similar but six genetically and behaviorally distinct mosquito species from the genus Anopheles. The complex contains the most important malaria vectors in Sub-Saharan Africa, especially those that transmit Plasmodium vivax, Plasmodium ovale, Plasmodium malariae, Plasmodium falciparum, and Plasmodium knowlesi [7]. A. gambiae is more than just a pest and a nuisance; it is also responsible for the spread of malaria and other deadly diseases across Africa.
Quinoline-based drugs were the mainstay of malaria treatment and prevention for decades. Unfortunately, the emergence of drug-resistant Plasmodium species as a result of mutation has rendered conventional therapeutic treatment ineffective [8,9]. WHO, on the other hand, has approved insecticides, fumigation, air shields, ultrasonic rays, pesticide spraying, insecticidal treated nets, and insecticidal-treated clothing as malaria prevention measures [10]. For a variety of reasons, insect repellents such as lotions, coils, mosquito-treated nets, and liquidators have limited efficiency [11,12]. However, consumers have recently become more interested in cheap plant-based commercial repellents, which are frequently considered as "safe" when compared to synthetic repellents; this is not always the case depending on the secondary metabolites present [13]. Plant-based metabolites have been shown to be an effective mosquito control alternative to synthetic insecticides or when used in conjunction with other pesticides as part of an integrated vector control strategy [14].
Several plant extracts have been reported to have mosquitocidal growth regulator activity, as well as insecticidal properties in the elimination of larval or adult mosquitos or as mosquito repellents to protect against mosquito bites [13,15,16]. Vitex negundo L. (Verbenaceae) is a well-known medicinal herb [17][18][19][20], which has been linked to a variety of pharmacological properties, including enzyme inhibition [21,22], antifeeding [23], larvicidal [24], and mosquito repellent activities [15]. Research on synthetic insect repellents has developed due to breakthroughs in our understanding of mosquito behavior and olfactory receptors [25], but there have been insufficient studies on the safety of plant-based mosquito repellents and their interaction with olfactory receptors. The olfactory system of insects is made up of a number of transmembrane odorant receptor proteins that are expressed in olfactory membrane neurons all throughout their bodies [26]. Odorant binding proteins (OBPs) are potent mosquito biosensors, according to Di Pietrantonio et al. [27], Sankaran et al. [28], and Possas-Abreu et al. [29]. Odorant-binding protein 1 (OBP1) (PDB ID 3N7H) and (OBP:PDB ID 2ERB), odorantbinding protein 7 (OBP7: PDB ID 3R1O), and odorant-binding protein 4 (OBP4: PDB ID 3Q8I) are the key olfactory proteins involved in signals for host recognition process, repellent alarm pheromone, and diphenyl ester-specific binding protein, respectively.
More systematic research is needed to better evaluate plant-based repellents and produce novel solutions that are safe for consumers. This study investigates the major constituents of V. negundo essential oil in order to establish repellent efficacy, predict their in-silico toxicity profile, and determine the interactions with Anopheles odorant binding proteins using a molecular docking-based technique.

Collection Sites and Identification of V. negundo Leaves
The leaves of V. negundo were harvested in September 2020 from six states of the North-Central Geopolitical Zone of Nigeria with the climatic condition and major soil type presented in Table 1

Leaf Processing and Extraction of Essential Oils
Collected fresh V. negundo leaves were washed with tap water and extracted within 12 h of collection using a 25 kg capacity fabricated Essential oil Distillation System (EDS) based on the steam distillation principle (Figure 1). The EDS steam generator was filled with 50 L of distilled water while the sample container was loaded to capacity and distilled over a period of 45 min. The distillate was recovered and separated in batches using a 2 L separatory funnel into essential oil and aqueous distillate (hydrosol), after which the essential oils were dried over anhydrous Na 2 SO 4 and stored for further analysis. Finally, the oil yield was calculated relative to the fresh matter and the result presented as the mean ± standard deviation of triplicate extractions.  Figure 1. Schematic of the Essential oil Distillation System (EDS).

GC-MS Profiling of the Essential Oils
The GC-MS analyses of the essential oils were performed with a Varian CP-3800 gaschromatograph equipped with a HP-5 capillary column (30 mm × 0.25 mm; coating thickness 0.25 μm), carrier gas nitrogen at 1.2 mL/min, and a Varian Saturn 2000 ion trap mass detector. The oven temperature was programmed from 50 to 280 at 3 °C/min. Analytical conditions: injector and transfer line temperatures were 220 and 240 °C, respectively. Volume injected: 0.2 μL of 10% hexane solution, split ratio 1:30. Co-injection of the essential oil with a solution containing a similar series of C8-C22 n-alkanes yielded linear retention indices for all molecules. Retention indices were used to identify the individual components, which were then compared to compounds previously reported in the literature [30,31]. Further, the identification of the compounds was made using data of a computer library (Wiley 275L) connected to the GC-MS, Adams library (https://bok.cc/book/3506611/3b1f4f (accessed on 15 October 2021)), the NIST website (https://webbook.nist.gov/chemistry/(accessed on 15 October 2021)) using RI values from comparable polarity columns, and/or the Mondello library (https://www.sisweb.com/software/wileyffnsc.htm (accessed on 15 October 2021)).

Sampling, Rearing, and Identification of Mosquitoes
In the months of March and April 2021, mosquito larvae were collected from selected localities in Kaduna metropolitan located between 10°33′ N of the equator and 07°27′ E of the Greenwich Meridian (Ungwan Gwari; 10°35.71′; 07°27.17′ Ungwan Romi 10°25.19′; 07°25.20′ Kamanzou 10°46.24′; 07°49.41′). The larvae were collected using a 7 cm diameter, 5 cm deep, and 30 cm long handle plastic standard dipper from a 0.12 m × 2.5 m deep temporary pool with grass vegetation. At the breeding sites, larvae identification and morphological classification were carried out. The absence of a siphon, the parallel swimming pattern on the water surface, and the morphology of the combs were used to sort larvae into the Anophelinae and Culicinae subfamilies under a compound microscope, compared to the Culex larva with a long siphon, lighter color, and "hairy" body, as well as the identification key of Gillies and Coetzee [32]. The immature larval stages were carefully transported in vials to the Insectary at the Biological Sciences laboratory, Kaduna State University, Nigeria. The larvae were then placed in an open plastic container 29 cm × 21 cm × 30 cm containing 1 L of ground water and allowed to acclimate for 2 h before being fed finely ground low-fat flour-baked food product [33].
The larvae were batched in separate breeding containers and were reared to adults in separate 30 cm × 30 cm wooden made net chambers for three weeks under controlled

GC-MS Profiling of the Essential Oils
The GC-MS analyses of the essential oils were performed with a Varian CP-3800 gaschromatograph equipped with a HP-5 capillary column (30 mm × 0.25 mm; coating thickness 0.25 µm), carrier gas nitrogen at 1.2 mL/min, and a Varian Saturn 2000 ion trap mass detector. The oven temperature was programmed from 50 to 280 at 3 • C/min. Analytical conditions: injector and transfer line temperatures were 220 and 240 • C, respectively. Volume injected: 0.2 µL of 10% hexane solution, split ratio 1:30. Co-injection of the essential oil with a solution containing a similar series of C8-C22 n-alkanes yielded linear retention indices for all molecules. Retention indices were used to identify the individual components, which were then compared to compounds previously reported in the literature [30,31]. Further, the identification of the compounds was made using data of a computer library (Wiley 275L) connected to the GC-MS, Adams library (https://b-ok.cc/book/3506611/3b1f4f (accessed on 15 October 2021)), the NIST website (https://webbook.nist.gov/chemistry/ (accessed on 15 October 2021)) using RI values from comparable polarity columns, and/or the Mondello library (https://www.sisweb.com/software/wiley-ffnsc.htm (accessed on 15 October 2021)).

Sampling, Rearing, and Identification of Mosquitoes
In the months of March and April 2021, mosquito larvae were collected from selected localities in Kaduna metropolitan located between 10 • 33 N of the equator and 07 • 49.41 ). The larvae were collected using a 7 cm diameter, 5 cm deep, and 30 cm long handle plastic standard dipper from a 0.12 m × 2.5 m deep temporary pool with grass vegetation. At the breeding sites, larvae identification and morphological classification were carried out. The absence of a siphon, the parallel swimming pattern on the water surface, and the morphology of the combs were used to sort larvae into the Anophelinae and Culicinae subfamilies under a compound microscope, compared to the Culex larva with a long siphon, lighter color, and "hairy" body, as well as the identification key of Gillies and Coetzee [32]. The immature larval stages were carefully transported in vials to the Insectary at the Biological Sciences laboratory, Kaduna State University, Nigeria. The larvae were then placed in an open plastic container 29 cm × 21 cm × 30 cm containing 1 L of ground water and allowed to acclimate for 2 h before being fed finely ground low-fat flour-baked food product [33].
The larvae were batched in separate breeding containers and were reared to adults in separate 30 cm × 30 cm wooden made net chambers for three weeks under controlled optimum conditions of 25 • C, 65% relative humidity, and regulated light/dark (14/10 h) cycle. The emerged adults were identified morphologically using taxonomic characters such as the palps, proboscis, wing venation, and markings or tuffs on legs or abdomen as provided by the dichotomous keys employed by Coetzee and Gillies [34]. This was performed using the simple Olympus light microscope to genera and species level. The adults in the cages were fed a 10% sucrose solution after eclosion from their pupal cases and allowed to rest and mature for 2 to 3 days. Only newly emerged adult females A. gambiae were manually aspirated into a 200 mL perforated plastic container and allowed to rest for 1 hr before exposure to the essential oils ( Figure 2). optimum conditions of 25 °C, 65% relative humidity, and regulated light/dark (14/10 h) cycle. The emerged adults were identified morphologically using taxonomic characters such as the palps, proboscis, wing venation, and markings or tuffs on legs or abdomen as provided by the dichotomous keys employed by Coetzee and Gillies [34]. This was performed using the simple Olympus light microscope to genera and species level. The adults in the cages were fed a 10% sucrose solution after eclosion from their pupal cases and allowed to rest and mature for 2 to 3 days. Only newly emerged adult females A. gambiae were manually aspirated into a 200 mL perforated plastic container and allowed to rest for 1 hr before exposure to the essential oils ( Figure 2).

Anopheles Species Authentication: Genomic DNA Extraction and PCR Amplification
The emerged adult mosquitoes belonging to the A. gambiae (s.l) complex were subjected to PCR and genomic DNA assays designed for species, molecular form identification, and molecular analysis to fully authenticate the adult Anopheles species. The molecular analyses were conducted at the Centre for Biotechnology Research and Training, Ahmadu Bello University, Zaria. Genomic DNA (gDNA) was extracted from 20 Anopheles

Anopheles Species Authentication: Genomic DNA Extraction and PCR Amplification
The emerged adult mosquitoes belonging to the A. gambiae (s.l) complex were subjected to PCR and genomic DNA assays designed for species, molecular form identification, and molecular analysis to fully authenticate the adult Anopheles species. The molecular analyses were conducted at the Centre for Biotechnology Research and Training, Ahmadu Bello University, Zaria. Genomic DNA (gDNA) was extracted from 20 Anopheles adult mosquitoes using the Quick-DNA™ Miniprep Plus Kit (D4069) product by ZYMO research company according to the protocol of the manufacturer.
Genomic mosquito DNA (2 µL) was extracted and placed in a 0.2 mL thin-walled Eppendorf tube and 23 µL of the PCR reaction mixture containing species specific primers for A. gambiae (GA = 5 -CTGGTTTGGTCGGCACGTTT-3 ), A. Arabiensis (AR = 5 -AAGTGTCC TTCTCCATCCTA-3 ), and a unique primer for all species (UN = 5 -GTGTGCCCC TTC-CTCGAT GT-3 ) according to the protocols of Scott et al. [35] and Favia et al. [36]. Then, deoxynucleotide triphosphates, magnesium chloride, PCR buffer, distilled water, and recombinant Taq DNA polymerase were all added. Amplification was performed in an initial denaturation step at 94 • C for two minutes, then a 30-cycle denaturation at 94 • C for 30 s, annealing at 49 • C for 30 s and elongation at 68 • C for 5 min in a thermal cycler machine. The PCR amplicons ran on a 1.5% agarose gel electrophoresis gel tank and were viewed under the trans-illuminator UV light for the characteristic A. gambiae positive band sizes at 390 bp.
Again, from the samples that showed base pairs of 390 band sizes (A. gambiae s.l), seven out of a group of 20 were randomly selected and 1 µL picked and mixed in 24 µL of PCR reaction mixture containing species specific primers for A. gambiae s.s and coluzzi. Species identification of A. gambiae (s.l.) was performed by PCR according to Favia et al. [36]. The PCR conditions were 10 min at 94 • C as the initial step, followed by 30 cycles (94 • C for 30 s, 53 • C for 30 s, and 68 • C for 30 s). After the last cycle, the products were finally extended for 5 min at 68 • C. Primers used in the PCR were: R5 (5 -GCC AAT CCG AGC TGA TAG CGC-3 ), R3 (5 -CGA ATT CTA GGG AGC TCC AG-3 ), Mop int (5 -GCC CCT TCC TCG ATG GCA T-3 ), and B/S int (5 -ACC AAG ATG GTT CGT TGC-3 ). Amplified fragments were analyzed on a 1.5% agarose gel for the characteristic band size of 475 for A. gambiae s.s.

Mosquito Behavioral Study Repellency Y-Tube Olfactometer Test
This experiment used an olfactometer with a glass Y-tube of 1.6 cm diameter, 12-cm base, and two 40.62 cm arms at a 45 • angle to one another [37]. V. negundo essential oil diluted with paraffin oil (analytical purity, Sigma-Aldrich) was added in a gradient sequence of 0.1, 0.25, 0.5, 0.75, and 1.0% v/v to a 5 g cotton ball, which was then mounted on the cotton ball holder of the test arm.
The cotton ball holder on the control arm of the glass Y-tube supported the negative control cotton ball (paraffin oil without the essential oil). At a rate of 180 mL/min, fresh air was provided by an electric air pump and filtered on a carbon and silicone base with air flowing into the respective arm of the Y-tube.
The essential oil was combined with the purified air and went through one arm of the Y-tube whereas purified air with no essential oil went through the opposite arm. Precisely, 100 female A. gambiae mosquitoes were discharged into the glass Y-tube, where behavioral responses were monitored and recorded for 1800 s [38]. In addition, to establish the mosquito behavioral activity and the degree of synergism between the pure constituents of the essential oil, selected commercially available pure constituents based on the absorption, distribution, metabolism, excretion, and toxicity (ADME/tox) and docking studies were evaluated. After each investigation, the Y-tubes were air cleaned with a stream of hot air (>60 • C), the cotton ball was removed, and the holder was cleaned. The olfactory test was repeated thrice. The repellent rate was calculated according to the following Equation (1) [39,40]: %Mosquitoe Repelled = 100 − mean number of mosquitoes selecting essential oil 100 − mean number of mosquitoes not selecting essential oil 100 (1) The % mosquito repelled of 50% individuals (ED 50 ) was estimated using the Probit analysis model in the IBM SPSS v.25 statistical software.

Target Protein Selection and Preparation
The odorant binding proteins (OBPs) were selected as a target based on their application as bio-recognition elements and biosensors for small ligands. The three-dimensional (3D) structures of four A. gambiae OBPs; OBP 1 (PDB ID 3N7H), OBP 7 (PDB ID 3R1O), Insects 2021, 12, 1061 7 of 26 OBP 4 (PDB ID 3Q8I), and OBP (PDB ID 2ERB) were retrieved from the Protein Data-bank (http://www.rcsb.org (accessed on 12 February 2021)) ( Figure 3). The crystal structures of the of OBPs were processed by removing existing ligands and water molecules while missing hydrogen atoms were added according to the amino acid protonation state at pH 7.0 utilizing Autodock 4.2 (Molecular Graphics Laboratory, Scripps Research Institute, La Jolla, CA, USA). Thereafter, non-polar hydrogens were merged while polar hydrogens were added to each protein. The process was repeated for each protein and subsequently saved into a dockable format.

Target Protein Selection and Preparation
The odorant binding proteins (OBPs) were selected as a target based on their application as bio-recognition elements and biosensors for small ligands. The three-dimensional (3D) structures of four A. gambiae OBPs; OBP 1 (PDB ID 3N7H), OBP 7 (PDB ID 3R1O), OBP 4 (PDB ID 3Q8I), and OBP (PDB ID 2ERB) were retrieved from the Protein Data-bank (http://www.rcsb.org (accessed on 12 February 2021)) ( Figure 3). The crystal structures of the of OBPs were processed by removing existing ligands and water molecules while missing hydrogen atoms were added according to the amino acid protonation state at pH 7.0 utilizing Autodock 4.2 (Molecular Graphics Laboratory, Scripps Research Institute, La Jolla, CA, USA). Thereafter, non-polar hydrogens were merged while polar hydrogens were added to each protein. The process was repeated for each protein and subsequently saved into a dockable format.

Molecular Docking Studies
Molecular docking was carried out using PyRx-Python Prescription 0.8 software (Hangzhou, Zhejiang, China). The input file was in the form of PDB code of the receptor or PDB file format and the molecules were in PDB file format. The output file was a docking report. The docked image was viewed by "BIOVIA Discovery Studio Visualizer" software (Waltham, MA, USA) to review the interactions between ligands and proteins, and the length of the interaction along with amino acids. The ligands were imported into PyRx 0.8 through the OpenBabel plug-in tool for each docking phase, with the Universal Force Field (UFF) as the energy minimization parameter, and conjugate gradient descent as the optimization algorithm. Table 1 shows the coordinate of the active sites of the four A.

Molecular Docking Studies
Molecular docking was carried out using PyRx-Python Prescription 0.8 software (Hangzhou, Zhejiang, China). The input file was in the form of PDB code of the receptor or PDB file format and the molecules were in PDB file format. The output file was a docking report. The docked image was viewed by "BIOVIA Discovery Studio Visualizer" software (Waltham, MA, USA) to review the interactions between ligands and proteins, and the length of the interaction along with amino acids. The ligands were imported into PyRx 0.8 through the OpenBabel plug-in tool for each docking phase, with the Universal Force Field (UFF) as the energy minimization parameter, and conjugate gradient descent as the optimization algorithm. Table 1 shows the coordinate of the active sites of the four A. gambiae odorant binding proteins as determined by the grid boxes utilized in the docking studies.
The inhibition constant and ligand efficiency metrics were quantitatively estimated using Equations (2)-(7) [49][50][51][52][53]. Ligand Ligand Efficiency Scaled (LE scale ) = 0.873e −0.026(HA) − 0.064 (4) Ligand Lipophilic Efficiency (LLE) = −(logk i + logP) Fit Quality (FQ) = LE LE scale (6) Ligand Efficiency Lipophilic Price(LELP) = Log P LE (7) 2.10. Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADME/tox) Investigation The ADME/tox filtering analysis of all selected ligands for docking including the physiochemical were predicted using the ADMETlab 2.0 webserver (https://admetmesh.scbdd. com/service/evaluation/index (accessed on 15 October 2021)). The canonical SMILES of the ligands downloaded from the PubChem Database were used for the calculation of the ADME/tox parameters in default mode. Table 2 indicates the general soil types, climatic conditions, temperature, and annual rainfall range of the collection sites, as reported in previous literature [54][55][56][57][58][59], and essential oil yield of the six states. Ferralsols and hydromorphic tropical soils were the predominant soil type of Kwara, while the soils of the Plateau are poorly drained sandy clay loam surfaces with clay subsurface. The Niger, Nasarawa, and Benue are dominated by sandy loamy soil types with a clay subsoil as well, except for Kogi soil which is richly sandy and loamy. The oil yields were generally high for all sampled states; however, V. negundo from Niger and Kogi produced the highest and lowest essential oil yields, respectively. Besides, there was no significant difference (p > 0.05) in the yields from Nasarawa, Plateau, and Kwara. In accordance with the findings of Tirillini et al., [60], the essential oil yield and composition were influenced by Ca 2+ and K + concentrations, percentage of organic matter, and temperature. The high yields in Niger, Nasarawa, Benue, Kwara, and Plateau may be attributable to their proximity to a region with a moderately low to low annual rainfall, low temperature, and predominant loamy soil with clay substructure with significant concentrations of K + , Na + , and Ca 2+ [61]. However, the significantly high yield obtained in Niger is justifiably due to the low annual rainfall (59.9-274.2 mm) compared to other regions with the same soil substructure (Table 2). On the other hand, the low yield recorded in the Kogi sample may be directly linked to the lack of clay minerals in the soil. Table 3 shows the results of the V. negundo essential oils GC-MS analysis collected from six states and the chromatograms are presented as Supplementary Figures S1-S6. The essential oil from the various collection sites showed compositional variation. Niger, Kwara, Benue, Plateau, Kogi, and Nasarawa oil, respectively, contained 16,18,30,24,15, and 28 identified compounds. Monoterpenes made up the majority of the constituents of the essential oils across all study states with about 74.65-96.23%, followed by about 0.75-16.32% sesquiterpene content. The rest were other compounds of about 3.55-10.88%. Essential oils from Niger, Kwara, and Kogi had more than 90.04% monoterpene while oils from Plateau, Nasarawa, and Benue demonstrated the highest sesquiterpene content (7.63-16.32%). Sesquiterpene content was found to be quite low in Niger and Kwara (0.75-1.23%), with sesquiterpene completely absent in the Kogi sample. The observed compounds are in consonant with the reports of Hebbalkar et al. [17], Huang et al. [62], and Kumar et al. [63].  Major components (≥10%) such as α-phellandrene were found in Niger (34.65%) and Kwara (20.27%), sabinene in Plateau (12.21%) and Nasarawa (11.31%), β-pinene (42.04%) and p-cymene (16.47%) in Kogi, and myrcene (16.78%) in Benue. As a result, the predominant monoterpenes (≥10%) present in V. negundo essential oil throughout the six states were α-phellandrene (20.27-34.65%), sabinene (11.31-12.21%), β-pinene (42.04%), p-cymene (16.47%), and myrcene (16.78%). In the chemotaxonomical classification of V. negundo, unique compounds such as α-pinene, linalool, cis-sabinene hydrate, citronellal, verbenone, and bornyl acetate found in all samples irrespective of the collection site could be used in the fingerprint of the essential oil. Several factors, including growth conditions, temperature, altitude, soil type, agricultural methods and practices, developmental stage, plant part extracted, and harvesting period, are strong factors that determine the presence and absence of certain terpenes, according to Moghaddam and Mehdizadeh [64]. Few monoterpenes and sesquiterpenes reported by Issa et al. [65] and Khokra et al. [31] were noticeably absent across all six states; however, compounds such as p-cymene [66], α-terpinolene [67], and citronellal [68] with reported potent insecticidal properties were present in the oils.

Mosquito Behavioural Study
The exposure of adult female A. gambiae to the essential oil from all six states and N-diethyl-3-methylbenzamide (DEET) for a period of 30 min at doses ranging from 1% v/v was investigated and reported in Figures 7 and 8.

Mosquito Behavioural Study
The exposure of adult female A. gambiae to the essential oil from all six states and N,N-diethyl-3-methylbenzamide (DEET) for a period of 30 min at doses ranging from 0.1-1% v/v was investigated and reported in Figures 7 and 8.
There was an increase in the number of mosquitos repelled with time with an optimal % repellency activity attained at approximately 15 min. All the essential oil samples showed a significant increase in the percentage of mosquitoes repelled within the period of investigation with no significant difference in the percentage of mosquitoes repelled (p > 0.05) between the essential oils and the N,N-diethyl-3-methylbenzamide. As the doses increased, the repellency activity increased to a concentration where there are no observable changes in activity. However, the optimal concentration varies from state to state as a result of its compositional variation (Table 3). Essential oils from Niger, Kwara, Plateau and Nasarawa showed optimal repellency at a concentration of 0.5% v/v while Niger, Benue, and Kogi oil samples showed an optimal effect at 0.75% v/v. On the contrary, DEET showed no significant difference (p < 0.05) in repellency activity as the concentration changed. This inference is in consonance with the studies of Cárdenas-Ortega et al. [69] and Senthil-Nathan [70], which emphasize the slight variation in the repellency activities of samples due to the presence and percentage composition of unique compounds. Insects 2021, 12, x FOR PEER REVIEW 13 of 27   There was an increase in the number of mosquitos repelled with time with an optimal % repellency activity attained at approximately 15 min. All the essential oil samples showed a significant increase in the percentage of mosquitoes repelled within the period of investigation with no significant difference in the percentage of mosquitoes repelled (p > 0.05) between the essential oils and the N, N-diethyl-3-methylbenzamide. As the doses increased, the repellency activity increased to a concentration where there are no observable changes in activity. However, the optimal concentration varies from state to state as a result of its compositional variation (Table 3). Essential oils from Niger, Kwara, Plateau and Nasarawa showed optimal repellency at a concentration of 0.5% v/v while Niger, Benue, and Kogi oil samples showed an optimal effect at 0.75% v/v. On the contrary, DEET showed no significant difference (p < 0.05) in repellency activity as the concentration changed. This inference is in consonance with the studies of Cárdenas-Ortega et al. [69] and Senthil-Nathan [70], which emphasize the slight variation in the repellency activities of samples due to the presence and percentage composition of unique compounds.

Effective Dose of the Essential Oils from the North-Central Geopolitical Zone
Using the Probit analysis model, the effective dose (ED50) that would repel 50% of the mosquito population is presented in Table 4. The ED50 of the oils and positive control are in the order of DEET > Kwara > Niger > Plateau and Nasarawa > Benue > Kogi. There is a significant different (p < 0.05) in the repellency of the oils from Kwara, Niger, Plateau, and Nasarawa compared to Benue and Kogi. The repellency property of the oils showed a composition-concentration dependent activity, which is not in variance with the result obtained in the mosquito behavioral investigation (Figures 7 and 8). However, DEET showed very potent repellency at ED50 of 0.01%v/v compared to the oils from all states. This observation is due to the variation in the composition.

Effective Dose of the Essential Oils from the North-Central Geopolitical Zone
Using the Probit analysis model, the effective dose (ED 50 ) that would repel 50% of the mosquito population is presented in Table 4. The ED 50 of the oils and positive control are in the order of DEET > Kwara > Niger > Plateau and Nasarawa > Benue > Kogi. There is a significant different (p < 0.05) in the repellency of the oils from Kwara, Niger, Plateau, and Nasarawa compared to Benue and Kogi. The repellency property of the oils showed a composition-concentration dependent activity, which is not in variance with the result obtained in the mosquito behavioral investigation (Figures 7 and 8). However, DEET showed very potent repellency at ED 50 of 0.01%v/v compared to the oils from all states. This observation is due to the variation in the composition.

Validation of Molecular Docking Protocol
According to the literature, a validated protocol must have a RMSD value < 2.0 in the binding mode prediction, when superimposed on the crystallographic pose of the ligand [71,72]. To establish that the conformation of the interaction between co-crystallized ligands and OBPs can be replicated in silico to validate our docking method, the cocrystallized ligands were redocked in the protein binding pocket and the root mean square deviation (RMSD) data were used to evaluate the fitness of each redocked pose. Figure 9 illustrates the poses estimated in relation to the deposited PDB complexes, with the RMSD of 0.67 Å, 0.24 Å, 0.71 Å, and 1.35 Å for OBP1, OBP 7, OBP 4, and OBP; respectively. binding mode prediction, when superimposed on the crystallographic pose of the ligand [71,72]. To establish that the conformation of the interaction between co-crystallized ligands and OBPs can be replicated in silico to validate our docking method, the co-crystallized ligands were redocked in the protein binding pocket and the root mean square deviation (RMSD) data were used to evaluate the fitness of each redocked pose. Figure 9 illustrates the poses estimated in relation to the deposited PDB complexes, with the RMSD of 0.67 Å, 0.24 Å, 0.71 Å, and 1.35 Å for OBP1, OBP 7, OBP 4, and OBP; respectively.

Molecular Docking
The binding energies and inhibition constants of the proteins with the selected ligands are reported in Table 5.

Molecular Docking
The binding energies and inhibition constants of the proteins with the selected ligands are reported in Table 5. All the selected ligands demonstrated multiplicity of binding properties and varying degrees of interaction within the active pockets of the proteins with the exception of OBP 4 which only has affinity for α-pinene, linalool, verbenone, and β-pinene. In Table 5, α-pinene and myrcene (−6.4 kcal/mol), citronellal (−5.5 kcal/mol), linalool (−5.4 kcal/mol), and myrcene (−5.8 kcal/mol) demonstrated the strongest affinities for OBP1, OBP7, OBP4, and OBP, respectively. As a result, myrcene has been identified as an inhibitor of the main olfactory proteins involved in the signals for host recognition processes (OBP1 and OBP), which aid in the collection and transport of hydrophobic odorants into and through the fluid. Repellent alarm pheromone and diphenyl ester-specific binding proteins were inhibited by citronellal (0.1528 mM) and linalool (0.1089 mM), respectively. According to Bohacek et al. [73] and Hughes et al. [74], a molecule with a low Ki value in the micromolar range but less than 10 −4 mM qualifies as a lead. The estimated Ki values of α-pinene and myrcene, citronellal, linalool, and myrcene against the four targeted OBPs are within the predicted range which qualifies them as a lead mosquito repellent agent. By implication, the strong binding energies and inhibition constants observed for α-pinene and myrcene, citronellal, linalool, and myrcene when compared to other investigated ligands may be suggestive of functional blocking of the olfactory receptor coreceptor, activation of specific odorant receptors (ORs), inhibition of specific ORs responding to attractants, and/or modulation of multiple ORs causing olfactory confusion, according to Tsitoura et al. [46] and Degennaro et al. [75]. Figure 10 illustrates the 3D images of the active pockets of the four selected OBPs while Table 6 reports the list of the active pockets of the four selected OBPs from A. gambiae. The structures of OBP, OBP1, and OBP7 are made up of two monomers, with OBP and OBP1 each having six α-helices and OBP7 having seven α-helices, with the odorant binding pocket placed in the center of a hydrophobic tunnel that runs through the dimeric interface. However, OBP4 is made up of a single monomer. In this investigation, the active pockets of odorant binding proteins were identified by removing the ligands that were previously linked to the receptors before targeting these cavities. The experimental ligands all docked at the same pockets as the native ligands, validating the docking protocol adopted in this study. Figures 11-14 depict the 3D interaction between OBPs and ligands, whereas Table 7  The structures of OBP, OBP1, and OBP7 are made up of two monomers, with OBP and OBP1 each having six α-helices and OBP7 having seven α-helices, with the odorant binding pocket placed in the center of a hydrophobic tunnel that runs through the dimeric interface. However, OBP4 is made up of a single monomer. In this investigation, the active pockets of odorant binding proteins were identified by removing the ligands that were previously linked to the receptors before targeting these cavities. The experimental ligands all docked at the same pockets as the native ligands, validating the docking protocol adopted in this study. Table 6. Active pockets on the four selected odorant binding proteins of A. gambiae.
The mechanisms of interaction between the various ligands differ and will most likely result in a variety of activities ranging from functional blocking of the olfactory receptor coreceptor due to repression of Leu73 in OBP1, inhibition of specific ORs responding to attractants, and/or modulation of multiple Ors causing disorientation, as reported by Murphy et al. [76]. A strong affinity of OBP7 for citronellal and myrcene, according to Sun et al. [77], could create disturbance in the insect's chemical information decoding potential. These rare interactions of α-pinene, linalool, verbenone, and β-pinene with OBP4 are strongly associated with their spatial orientation of the dialkyl and π-alkyl groups; with the likelihood of blocking the olfactory receptor coreceptor and modulate OBP4, resulting in the susceptibility of A. gambiae to these molecules in the repellent.

Efficiency Metrics of Selected Ligands
Supplementary Tables S1-S4 show the ligand efficiency metrics of the selected ligands, which were calculated using Equations (3)- (7). 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 hit quantitatively [78]. During lead discovery, the Ligand Efficiency Lipophilic Price (LELP) is estimated to be between −10 and 10 [79]. The ligand efficiency metrics against the four OBPs are within the criteria, qualifying them as a possible odorant binding protein repellent lead. The ADMET properties of all selected ligands were carried out to determine the molecules as safe potential OBP inhibitors and the results are presented in Tables S5-S16. Even though cis-sabinene hydrate, citronellal, sabinene, and verbenone failed the human oral bioavailability test (values were 0.7-1.0), this pharmacokinetic parameter is less of a concern regarding skin sensitization and eye irritation for dermally applied products such as repellent lotion or aerosols. In this investigation the empirical decision for skin sensitization and eye irritation tests for linalool, cis-sabinene hydrate, citronellal, sabinene, verbenone, α-terpinene, bornyl acetate, β-pinene, and α-phellandrene are >0.8, which is within the rejection zone because such molecules could induce allergic contact dermatitis, cornea, and conjunctiva tissue damage [80]. Furthermore, citronellal and α-phellandrene have been discovered to be respiratory and human hepatotoxicants, respectively, with high morbidity and mortality potential [81]. The plasma protein binding and blood-brain barrier penetration of cis-sabinene hydrate, α-phellandrene, and α-terpinene were found to be greater than 90%, indicating that these compounds have a low therapeutic index [82]. The metabolic profile of the ligands indicated that they are all either substrates or inhibitors of human cytochrome P450 based on chemical biotransformation reactions [83].

In-silico Environmental Toxicity
To estimate the environmental impact of the essential oil, and the bioconcentration factor (BCF), the concentration of the selected ligands in water in mg/L that causes 50% growth inhibition of Tetrahymena pyriformis after 48 h (IGC 50 ), 50% of fathead minnow to die after 96 h (LC 50 ), and 50% of Daphnia magna to die after 48 h (LC 50 DM) were evaluated. The result of the analysis is presented in Table 8. Table 8. In-silico environmental toxicity profile of the selected ligands.  Bioconcentration factors range from 0.553-3.003 L/kg, reflecting the very low potential for the ligands to bioaccumulate in the environment. Citronellal, verbenone (with high human oral bioavailability) (Tables S10 and S12), linalool, and myrcene demonstrated the lowest BCF compared to β-pinene, p-cymene, sabinene, α-pinene, and cis-sabinene hydrate ( Table 8)

Repellence Study of the Pure Selected Ligands
Pure α-pinene, linalool, and myrcene could be employed as a safe active ingredient in the development of a new mosquito repellent, according to ADME/tox and docking studies. Exactly 1.0% v/v of pure α-pinene, linalool, myrcene, equivalent mixture of α-pinene, linalool, myrcene (PLM), and N-diethyl-3-methylbenzamide (DEET) were evaluated for repellent activities and reported after exposure for 30 min (Figure 15). Bioconcentration factors range from 0.553-3.003 L/kg, reflecting the very low potential for the ligands to bioaccumulate in the environment. Citronellal, verbenone (with high human oral bioavailability) (Tables S10 and S12), linalool, and myrcene demonstrated the lowest BCF compared to β-pinene, p-cymene, sabinene, α-pinene, and cis-sabinene hydrate ( Table 8)

Repellence Study of the Pure Selected Ligands
Pure α-pinene, linalool, and myrcene could be employed as a safe active ingredient in the development of a new mosquito repellent, according to ADME/tox and docking studies. Exactly 1.0% v/v of pure α-pinene, linalool, myrcene, equivalent mixture of αpinene, linalool, myrcene (PLM), and N-diethyl-3-methylbenzamide (DEET) were evaluated for repellent activities and reported after exposure for 30 min (Figure 15). The pure linalool and myrcene demonstrated above 50% mosquito repellency activity except for α-pinene. This repellence activity is comparatively lower than the observed The pure linalool and myrcene demonstrated above 50% mosquito repellency activity except for α-pinene. This repellence activity is comparatively lower than the observed activity of the essential oil after 30 min exposure (Figure 7). However, there is a significant difference in the mosquito repellence activity of DEET (100%) compared to the pure ligands (42-60%). Despite the strong interactions of the ligands with the OBPs the observed repellency activity of the essential oil is significantly synergistic in nature. The synergistic activity was validated by comparing the repellence activity of PLM (73%) with DEET (100%) after 30 min of exposure.

Conclusions
The essential oils of V. negundo from six states in the North Central Geopolitical Zone were found to have composition-dependent mosquito repellent efficacy against A. gambiae. In comparison to the other states, the essential oil derived from Niger, Kwara Plateau, and Nasarawa demonstrated substantial repellency with an ED 50 of 0.14-0.08% v/v. αphellandrene, sabinene, β-pinene, p-cymene, and myrcene were the most common terpenes found in the essential oil throughout the six states. Regardless of the collection site, all essential oils contained α-pinene, linalool, cis-sabinene hydrate, citronellal, verbenone, and bornyl acetate. Linalool, α-pinene, verbenone, β-pinene, myrcene, and citronellal had the strongest affinity for OBPs, while -pinene, citronellal, linalool, and myrcene inhibited strongly by creating hydrophobic interactions at the binding pocket. The LE, LLE, FQ, and LELP values were all within the predicted ranges, indicating that the ligands are quantitatively hit and so qualify as a potential odorant binding protein repellent lead. Linalool, cis-sabinene hydrate, citronellal, sabinene, verbenone, α-terpinene, bornyl acetate, α-pinene, and α-phellandrene all had a low ecotoxicological profile, while linalool, cissabinene hydrate, citronellal, sabinene, verbenone, and α-terpinene did not. According to ADME/tox and docking results, α-pinene, linalool, and myrcene could be used as safe active components in the development of an environmentally friendly new mosquito repellent. Commercial standards of α-pinene, linalool, and myrcene were themselves active in mosquito repellent assays, and a mixture containing these compounds in equivalent proportions was found to be as significantly active as DEET, suggestive of a synergistic activity itself. Docking showed that these ligands bind to OBPs and may play an important role in blocking the olfactory receptor (ORs) coreceptor and inhibition of specific ORs causing disorientation and confusion in A. gambiae.