Abstract

Schinus molle is a medicinal plant used as an anti-inflammatory and for rheumatic pain in the traditional medicine of Peru. On the other hand, Aedes aegypti is the main vector of several tropical diseases and the transmitter of yellow fever, chikungunya, malaria, dengue, and Zika virus. In this study, the aim was to investigate the antioxidant activity in vitro and the insecticidal activity in silico, in the presence of the mosquito juvenile hormone-binding protein (mJHBP) from Aedes aegypti, of the essential oil from S. molle leaves. The volatile phytochemicals were analyzed by gas chromatography-mass spectrometry (GC-MS), and the profile antioxidants were examined by DPPH, ABTS, and FRAP assays. The evaluation in silico was carried out on mJHBP (PDB: 5V13) with an insecticidal approach. The results revealed that EO presented as the main volatile components to alpha-phellandrene (32.68%), D-limonene (12.59%), and beta-phellandrene (12.24%). The antioxidant activity showed values for , , and . Regarding the insecticidal approach in silico, alpha-muurolene and gamma-cadinene had the best biding energy on mJHBP (), followed by beta-cadinene (). Additionally, the volatile components did not reveal antioxidant activity, and its potential insecticidal effect would be acting on mJHBP from A. aegypti.

1. Introduction

Essential oils (EOs) are a class of natural products characterized by their volatile components, especially terpenes and sesquiterpenes [1]. EOs are obtained by several methods such as hydrodistillation, supercritical fluids, microwaves, and ultrasound [2]. Although the yield rate is very low, a majority are used in perfumery [3], aromatherapy [4], and the pharmaceutical [5] and food industries [6]. Their biological activities have been reported as insecticidal, antiviral, antibacterial, antifungal, antimalarial, antidepressive, anticancer, antimutagenic, hepatoprotective, anti-inflammatory, antioxidant, analgesic, and antipyretic [5, 7]. Schinus molle L., internationally commonly known as peruvian pepper, false pepper, American pepper, or pink pepper, is a perennial tree belonging to the Anacardiaceae family, which is native to subtropical regions of South and Central America [8]. In Peru, it is popularly known as molle (in Spanish) and its leaves and fruits exhibit anti-inflammatory, analgesic, antirheumatic, antibacterial, antiseptic, and repellent activity [9]. The essential oil of S. molle (leaves and fruits) exhibits biological activities such as antibacterial (gram-positive and gram-negative bacteria), antifungal, insecticidal, repellent, and cytotoxic activities [10]. Locally, the plant is used as a rheumatic reliever, analgesic, antiseptic, purgative, stomach cramp reliever, and diuretic. Moreover, leaves, fruits, and latex are used for many conditions, including menstrual, respiratory, and urinary disorders, and as an antidepressant, digestive stimulant, and astringent [11].

Regarding insecticidal activity, several synthetized chemicals are used as repellents, vector control, or insecticides, but during the last year, some of them have been retired due to their carcinogenic, genotoxic, or harmful-to-the-environment effects by bioaccumulation [12]. The application of natural products as insecticides might be highly effective, less expensive, biodegradable, and safer than that of synthetic insecticides [13]. The harmful effect of essential oils or their isolated compounds against insects can be manifested in various ways, including mortality, toxicity, inhibiting growth, the suppression of reproductive behavior, and reducing fertility and fecundity [14]. Hence, the use of bioinsecticides based on natural products might be a powerful tool to combat some insects who are responsible for several tropical diseases (Zika, malaria, dengue, and yellow fever) transmitted mainly by the mosquito Aedes aegypti [15].

Recently, a group of salivary D7 proteins known as mosquito juvenile hormone-binding protein (mJHBP) has been identified in A. aegypti, which consist of two modified odorant-binding protein domains with its respective ligand in the N-terminal domain named JH-III [16]. This hormone is found in pupae and adults; furthermore, it regulates the larval development in the mosquito. On the other hand, the juvenile hormone (JH) plays an essential role in adult female promoting the ovary maturation before blood feeding; additionally, it is involved with the nutritional state and its relationship with blood meal-dependent reproductive development [17]. In recent years, a class of insecticides named JH analogues have been designed to disrupt this endocrine process and affect the normal development in mosquitoes. A synthetic JH analogue is pyriproxyfen, which is a phenyl carbonyl derivative, and its main effect is to produce an imbalance in the mosquito hormonal system, inhibiting the embryogenesis, adult metamorphosis, and development of the adult mosquito [18].

Although the insecticidal activity of S. molle EO was evaluated in several species such as Trogoderma granarium, Tribolium castaneum [19], Sitophilus oryzae [20], Haematobia irritans [21], Ctenocephalides felis felis [22], and Gonipterus platensis [23], currently, there is not scientific literature on A. aegypti. Thus, we investigated the potential insecticidal effect in silico of the volatile components of S. molle EO to find any responsible molecule as a JH analogue, which might combat the presence of mosquitoes in tropical regions and reduce the prevalence of diseases transmitted by them. The aim in this study was (1) to determine the total volatile component of the EO from S. molle leaves by gas chromatography-mass spectrometry (GC-MS), (2) to evaluate the antioxidant activity using the DPPH, ABTS, and FRAP methods, and (3) to determine the insecticidal activity using a virtual screening of the EO from S. molle leaves on the mosquito juvenile hormone-binding protein from Aedes aegypti.

2. Materials and Methods

2.1. Chemicals

All the solvents (dichloromethane, chloroform, methanol, and hydrogen peroxide), of analytical grade (99.5%), were purchased from Merck (Darmstadt, Germany). 2,2-Diphenyl-1-picrylhydrazyl (DPPH), 2,2-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) (ABTS), 2,4,6-tri(2-pyridyl)-s-triazine (TPTZ), and Trolox were purchased from Sigma-Aldrich (St. Louis, MO, USA).

2.2. Plant Material

A quantity of 4700 g of Schinus molle (leaves) cultivated in Tinguiña, Ica Region, Peru (406 m.a.s.l.) in December 2020 was received. Leaves were cleaned and peeled to be incorporated in a Clevenger equipment and to obtain essential oil by hydrodistillation for 2 h [24]. The essential oil was separated by decantation; then, anhydrous Na2SO4 was added to eliminate any water drops. Finally, the EO was stored in a sealed amber vial until further use.

2.3. Identification of Volatile Compounds by Gas Chromatography-Mass Spectrometry (GC-MS)

Volatile chemicals were determined with a GC-MS system (Agilent Technologies 7890 Gas Detector and Agilent Technologies 5975C Mass Spectrometer Detector, Santa Clara, CA, USA). Then, 20 μL of EO was mixed with 1.0 mL of dichloromethane. Next, 1.0 μL of the working solution was injected into the equipment in splitless mode (split: 20 : 1). The EO was run on a J&W 122-1545.67659 DB-5 ms column, () (Agilent Technologies, Santa Clara, CA, USA). The working conditions were as follows: the temperature program was 40°C, starting with increments of 5°C/min up to 180°C, followed by increases of 2.5°C/min up to 200°C for 5 min, and finally 10°C/min up to 300°C, remaining for 3 min. The helium flow rate was at 1 mL/min. Volatile components’ identification was based on a comparison of relative retention indices (RIs) and mass spectra data with the NIST20 library data and the published literature [25]. Each RI was calculated compared with a homologous series of n-alkanes C9–C25 (C9, BHD purity 99%; C10–C25, Fluka purity 99%). The relative amount (expressed as a percentage) of each compound identified in the EO was calculated by comparing the area of the corresponding peak in the chromatogram with the total area of identified peaks.

2.4. Determination of the Antioxidant Capacity by the Free Radical 2,2-Diphenyl-1-picrylhydrazyl (DPPH)

A 150 μL of EO (10 mg/mL) was mixed with 2850 μL of a methanolic solution of DPPH (20 mg/L) with an absorbance adjusted to . After mixing, it was incubated in the dark for 30 minutes, and the absorbance reading was carried out at 515 nm. The standard curve was elaborated with Trolox at concentrations of 0 to 800 μmol/mL. The Trolox equivalent antioxidant capacity (TEAC) was expressed as μmol TE/g of essential oil [26]. To calculate the half inhibitory concentration (IC50), a linear regression method was used based on the different concentrations of the EO.

2.5. Determination of the Antioxidant Capacity by the Method of the Radical 2,2-Azinobis-(3-ethylbenzothiazoline)-6-sulfonic Acid (ABTS.+)

To carry out the antioxidant activity using the ABTS radical, a solution was prepared using the mixing of 10 mL of ABTS (4.06 mg/mL) and 10 mL of potassium persulfate (0.7 mg/mL), both reacted for 16 hours. Then, 150 μL of the EO (5 mg/mL, diluted with methanol) was mixed with 2850 μL of ABTS radical and incubated for 7 min, the same procedure was used with Trolox standard ranging between 0 and 400 μmol/mL. The absorbance was read to at a wavelength of 734 nm. The Trolox equivalent antioxidant capacity (TEAC) was expressed as μmol TE/g of essential oil [27]. To calculate the half inhibitory concentration (IC50), a linear regression method was used based on the different concentrations of the EO.

2.6. Determination of Antioxidant Capacity by the Ferric Reducing/Antioxidant Power (FRAP)

In the determination of FRAP of EO, a reagent battery was made using 25 mL of acetate buffer pH 3.6; 2.5 mL of 20 mM TPTZ was dissolved in 40 mM HCl and 2.5 mL of 20 mM ferric chloride hexahydrate. The mixing of this reagents constituted the FRAP solution which reacted with the EO at different concentrations as well as Trolox standard. 150 μL of EO (1 mg/mL) was mixed with 2850 μL of FRAP reagent and reacted for 4 min at room temperature. The absorbance reading was carried out at 593 nm. A standard curve was prepared with Trolox (50-800 μM). The results were expressed as μmol equivalents of Trolox per gram of essential oil (μmol TE/g EO) [28]. To calculate the half inhibitory concentration (IC50), the linear regression method was used based on the different concentrations of the EO.

2.7. Molecular Docking Studies

The molecules isolated from phytochemicals of essential oil were docked against the mosquito juvenile hormone-binding protein (PDB id: 5V13). In order to validate the docking, the crystal structures were docked with the native ligand pyriproxifen and JH3 bound to X-ray structures of mosquito juvenile hormone-binding protein. In addition to the former for comparison, better validation of docking at the binding cavity of mosquito juvenile hormone-binding protein was determined surrounding the 3 Å distance from the bound ligand JH3 using the DoGSiteScorer server. The binding cavity residues comprised mainly TYR33, LEU37, TRP50, AL51, TRP53, TYR64, SER69, TYR129, TYR133, ILE140, PHE269, TRP278, and ALA281. Protein and ligand preparations were performed using AutoDock Tools (v. 1.5.6) (Forli et al., 2016). Gasteiger charges were added to the ligand molecules prior to converting to PDBQT format. The online server DoGSiteScorer and the information about the binding site residues of the native ligand were used to construct the grid box. The grid box of dimensions for the mosquito juvenile hormone-binding protein with 0.375 Å grid spacing was constructed using AutoGrid 4.2. Semiflexible docking was performed keeping the receptor molecule rigid and ligands flexible. AutoDock 4.2 using the Lamarckian genetic algorithm (LGA) scoring function with number of , , and maximum number of was used to develop the molecular docking of all volatile components of EO. Then, the RMSD clustering maps were obtained by a reclustering command with a clustering tolerance of 0.25 Å, 0.5 Å, and 1 Å, respectively, in order to obtain the best cluster having the lowest energy score with a high number of populations [18].

2.8. Molecular Dynamics Simulation (MD)

This evaluation was carried on the dock complexes for α-phellandrene (the most abundant in GC-MS analysis) and α-muurolene (the most active in molecular docking analysis) with the mosquito juvenile hormone-binding protein (mJHBP) using the Desmond 2020.1 from Schrödinger, LLC. The reproducibility was completed using three replicates for each MD run. The OPLS-2005 force field and explicit solvent model with the SPC water molecules were used in this system. To neutralize and simulate the physiological conditions, the charge Na+ ions and 0.15 M NaCl were added, respectively. The system was equilibrated using constant-temperature and constant-volume ensemble (NVT) for 100 ps. This was followed by a short run equilibration and minimization using constant-temperature and constant-pressure ensemble (NPT) for 12 ps and was set up using the Nose-Hoover chain coupling scheme with temperature 27°C, the relaxation time of 1.0 ps, pressure 1 bar, and a time of 2 fs. The Martyna-Tuckerman-Klein chain coupling scheme barostat method was used for pressure control using a relaxation time of 2 ps. The particle mesh Ewald method was used for calculating long-range electrostatic interactions with a radius of 9 Å for the coulomb interactions. Bonded forces were calculated using a RESPA integrator with a time step of 2 fs for each trajectory. The root mean square deviation (RMSD), radius of gyration (Rg), root mean square fluctuation (RMSF), and solvent accessible surface area (SAS area) were calculated to monitor the stability of the MD simulation.

2.9. Statistical Analysis

In the analysis of the antioxidant profile of essential oil from S. molle, the IC50 values were estimated by linear regression statistics. A Spearman’s Rho coefficient was calculated to establish a correlation between the evaluated concentrations and the antioxidant response. values less than 0.01 were considered statistically significant. GraphPad Prism program version 6.0 (La Jolla, CA, USA) was used to carry out the statistical analysis.

3. Results and Discussion

3.1. Chemical Profile of the Essential Oil of S. molle

The obtained EO of fresh leaves showed a light-yellow color, a density of at 20°C, and an extraction yield of 0.73% during 2 hours of distillation. The profiles of volatile constituents of S. molle essential oil were analyzed by GC-MS and are presented in Figure 1 and Table 1. Regarding the chemical profile, the EO showed 34 compounds (Figure 1), two of which are of unknown structures (compounds 27 and 29), which accounted for 99.7% of the total composition. Those unknown chemical structures were classified as oxygenated sesquiterpenes with the following formula C15H26O and a molecular weight of 222.37 g/mol. The analysis identified alpha-phellandrene as the main volatile chemical with 32.68%, followed by D-limonene (12.59%) and beta-phellandrene (12.24%). According to Figure 1, the retention time at 17.80 min corresponded to the major component in the EO (compound 7).

In the phytochemical analysis, our results differed from those of other investigations and could be explained by different types of factor such as type of extraction, temperature conditions, storage conditions, edaphic conditions, plant part extracted, ecosystem, environmental, harvesting time, and the presence of environmental pollution [29]. In this study, the most representative molecule was alpha-phellandrene with 32.68%. However, in recent studies, as reported by Morales-Rabanales et al., β-phellandrene (15.7%) and α-phellandrene (12.1%) were the main components of the EO from S. molle leaves cultivated in Tepetitla de Lardizábal, Mexico [30]. In an analysis of S. molle EO obtained of dried leaves from Seropédica, Rio de Janeiro, Brazil, epi-α-cadinol (22.85%) was the major component [31]. Additionally, in S. molle from Caxias do Sul, Brazil, alpha-pinene () was the most abundant constituent [32]. In Naviraí-MS, Brazil, the major component was epi-α-cadinol with , followed by myrcene () and sabinene () [33]. S. molle from North Cyprus showed values of alpha-phellandrene, limonene, and beta-phellandrene equivalent to 31.5%, 10.1%, and 9.9%, respectively [34]. On the other hand, S. molle from Sonora, Mexico, revealed 60.8% of D-limonene (17.8%), o-cymene (16.4%), β-phellandrene (12.6%), δ-cadinene (7.5%), and caryophyllene (6.5%) [35]. However, in Jordan and Turkey, alpha-phellandrene was the most abundant constituent with 48.2% and 29.0%, respectively.

3.2. Antioxidant Profile of S. molle Essential Oil

S. molle EO exhibited a low antioxidant activity, as is shown in Table 2. On the other hand, there was a significant difference between the different methods used in the antioxidant activity (). Other reports have shown different values; according to Eryigit et al., EO showed a Trolox of S. molle grown in Turkey [36]. EO from Portugal reported an inhibition value of 4.8% at 16 mg/mL and was obtained by hydrodistillation. Although alpha-phellandrene, myrcene, limonene, and beta-phellandrene were the main components in this EO, its low antioxidant activity could be associated with the low ability of monoterpene hydrocarbons for DPPH scavenging activity. However, in other methods, such as beta-carotene bleaching, the antioxidant activity was high and might be justified by the presence of the monoterpenes α and β-phellandrene, α-pinene, sabinene, limonene, β-myrcene, and others [10]. In this study, the IC50 against DPPH was 41.84 μg/mL, which was higher compared to the findings of the EO from southeast Portugal [10]. In another study, S. molle from Egypt showed an IC50 at DPPH of 172.45 μg/mL and a high percentage of oxygenated sesquiterpene (17.73%) [37]. According to Table 2, the best effect was found in the FRAP assay, followed by ABTS and DPPH. It is known that DPPH and ABTS methods are based on electron and H atom transfer reaction, while the FRAP method is based on electron transfer reaction. Thus, the volatile phytochemicals from S. molle EO could be acted by electron transfer to reduce the free radicals.

3.3. Molecular Docking of the Essential Oil from S. molle in the Presence of the Mosquito Juvenile Hormone-Binding Protein (mJHBP) from Aedes aegypti

The virtual screening of the 32 chemical components, additionally JH3 and the pyriproxyfen (synthetic insecticide), was carried out in order to understand the interaction profile of various volatile compounds present in S. molle EO leaves with mJHBP from A. aegypti. Out of 32 specific compounds found abundantly in chromatography, gamma-cadinene and alpha-muurolene displayed a lower binding energy () of -9.7 kcal/mol and predicted inhibitory concentration (Ki) of 0.134 μM (Table 3). Although the major component (alpha-phellandrene) had a high binding free energy (), compared to alpha-muurolene, the total components might be synergizing the insecticidal effect.

The principal residues of mJHBP, TYR133, TRP53, TYR33, and PHE144 were involved in Pi-Sigma bond formation with alpha-muurolene (Figure 2). We observed that TYR133 was the main residue of mJHBP in which the molecules of the EO established a Pi-sigma, Van der Walls, and Pi-alkyl bonds (Supplementary Materials S1–S32). Interestingly, elixene, germacrene D, beta-cadinene, gamma-muurolene, germacrene-D, humulene, and gamma-eludesmol exhibited activity and binding energies close to those of the known insecticide pyriproxyfen and the cocrystallized ligand JH3. Alpha-phellandrene, beta-caryophyllene, beta-elemene, alpha-pinene, beta-phellandrene, beta-pinene, elemol, etc. have exhibited significant binding and inhibition of the mosquito juvenile hormone-binding protein. Therefore, from the docking study, it can be predicted that the molecules from this essential oil have great potential as inhibitors of the mosquito juvenile hormone-binding protein. The interactions of JH3 in this study are according to the results of Ramos et al., which the epoxy group forms a conventional hydrogen bond with the phenolic hydroxyl of Tyr-129, and other interactions showed on the isoprenoid chain were Val65,Val68, Pro55, Phe144, Tyr64, Tyr33, Ala281, Trp53, and Phe269 [18]. Furthermore, the low binding free energy () of some volatile components was evidenced in those interactions with the Tyr133, Trp53, Tyr33, and Phe144 residues as shown in Supplementary Materials Figures S1–S32. These could infer that the insecticidal activity is related to these residues mentioned above.

3.4. Molecular Dynamics of the Phytoconstituents of the Essential Oil from S. molle

Molecular dynamics and simulation (MD) studies were carried out in order to determine the stability and convergence of mJHBP+α-phellandrene and mJHBP+α-muurolene complexes. Each simulation of 100 ns displayed stable conformation comparing the root mean square deviation (RMSD) values. The Cα backbone of mJHBP bound to α-phellandrene exhibited a deviation of 0.5 Å (Figure 3(a)), while mJHBP bound to α-muurolene exhibited a deviation of 1.0 Å (Figure 3(f)). RMSD plots are within the acceptable range signifying the stability of mJHBP in the ligand-bound state before and after simulation, and it can also be suggested that ligand α-phellandrene and α-muurolene-bound mJHBP are quite stable in complex due to the higher affinity of the ligand. The plots for root mean square fluctuations (mMSF) displayed a significant spike in fluctuation (3.5 Å) at amino acid residue 220 in mJHBP, while the rest of the residues fluctuated less during the entire 100 ns simulation (Figures 3(b) and 3(g)) indicating the stable amino acid conformations during the simulation time. Therefore, from the RMSF plots, it can be suggested that the structures of mJHBP were stable during simulation in α-phellandrene- and α-muurolene-bound conformations. The radius of gyration is the measure of compactness of the protein. In this study, the mJHBP Cα backbone displayed a lowering of the radius of gyration (Rg) from 19.6 Å to 19.4 Å, and the lowering of the Rg indicates the compactness of the complex (Figure 3(c)). On the other hand, the lowering of Rg was observed for the mJHBP+α-muurolene complex till 35 ns, and a later increment of the peak was observed. This indicates the less stable conformation of mJHBP+α-muurolene as compared to mJHBP+α-phellandrene. However, the stable Rg peak confirmed the significant compactness of the protein in the α-muurolene-bound state (Figure 3(h)). The overall quality analysis from RMSD and Rg suggests that α-phellandrene and α-muurolene bound to the protein targets posthumously in the binding cavities and played a significant role in stability of the proteins. Solvent accessible surface area provides the information about the compactness of protein complex with the ligand. The lowering of SASA in the case of α-phellandrene and α-muurolene bound to mJHBP as compared to the unbound state signifies the achievement of stable converged structures due to the high compactness of both the systems (Figures 3(d) and 3(i)). The interaction plots of both α-phellandrene and α-muurolene bound to mJHBP displayed no involvement of conventional hydrogen bonds, while hydrophobic interactions played a major role in ligand stabilization at the binding cavity of mJHBP (Figures 3(e) and 3(j)).

4. Conclusions

The essential oil of S. molle leaves presented 34 volatile constituents determined by gas chromatography-mass spectrometry (GC-MS). Alpha-phellandrene was the major component, which represented 32.68% of the total composition; furthermore, beta-phellandrene and D-limonene were the following more abundant metabolites with 12.24% and 12.59%, respectively. Additionally, the antioxidant activity was evaluated in order to determine the antioxidant profile of EO using the three methods of DPPH, ABTS, and FRAP. According to the results, the EO showed better affinity and a good effect on FRAP assay and IC50 equivalent to . However, this study revealed that S. molle EO is not a good antioxidant. Regarding the insecticidal activity in silico based on a virtual screening on mosquito juvenile hormone-binding protein, several volatile compounds were active against this target such as alpha-muurolene and gamma-cadinene, being similar to pyriproxyfen, which is a synthetic insecticide analogue of JH. The molecular dynamics carried out for alpha muurolene (the best result in molecular docking) and alpha-phellandrene (the most abundant molecule) were very stable during 100 ns of evaluation. In the future, S. molle EO might be used as bioinsecticide on A aegypti, but an in vitro and in vivo assay has to be evaluated to validate our findings.

Data Availability

The data used to support the findings of this study are included within the supplementary information file.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Conceptualization was supervised by O.H.-C., H.C., G.E.-S.B., and N.H.A. Methodology was supervised by S.A. and H.C. Formal analysis was supervised by E.C.E.-R. Investigation was supervised by O.H.-C., P.W.C.-V., and R.D.H.-Q. Writing the original draft preparation was supervised by O.H.-C., J.B.P.-O., and E.L.G. Writing the review and editing was supervised by O.H.-C. and L.F.-S. Funding acquisition was supervised by N.H.A. and S.A.

Acknowledgments

This work was funded by Princess Nourah Bint Abdulrahman University Researchers Supporting Project number (PNURSP2022R62), Princess Nourah Bint Abdulrahman University, Riyadh, Saudi Arabia. Researchers Supporting Project number (RSP-2021/26), King Saud University, Riyadh, Saudi Arabia.

Supplementary Materials

Figures S1–S32: molecular interaction studies of the phytochemical components from S. molle and mosquito juvenile hormone-binding protein (mJHBP) from Aedes aegypti. (Supplementary Materials)