Prerpints.org logo
Preprint
Article

Advanced Formulation of Ecological Biopesticides Based on Citrus Limonum in Clayey Matrices: Optimization of Diffusive Dynamics

This version is not peer-reviewed.

Submitted:

24 October 2024

Posted:

28 October 2024

You are already at the latest version

Abstract
This study investigates the innovative use of natural porous clays from the Bejaad region in Morocco as a support matrix for the encapsulation and controlled release of lemon essential oil (Citrus limonum, EOCL), a natural compound with well-documented insecticidal properties. The research aims to address the inherent challenges of essential oils, particularly their high volatility and rapid degradation, by improving their stability and insecticidal efficiency against the grain pest Sitophilus granarius. By anchoring EOCL onto clay matrices, the study seeks to achieve a sustained and controlled release of the active components, thereby enhancing their practical application as biopesticides. The clays were comprehensively characterized using advanced analytical techniques, including X-ray diffraction (XRD), X-ray fluorescence (XRF), Fourier-transform infrared spectroscopy (FT-IR), scanning electron microscopy with energy-dispersive X-ray analysis (SEM-EDX), and thermogravimetric analysis (TGA). These techniques revealed the mineralogical composition, thermal properties, and morphology of the clays, demonstrating their suitability for effectively adsorbing and retaining EOCL. The insecticidal performance of the clay/EOCL composites was rigorously tested under controlled conditions, revealing a marked improvement in efficacy, with significantly lower lethal doses required to achieve high mortality rates in Sitophilus granarius. The diffusion of EOCL through the clay matrix was modeled using Fick law of diffusion, and the results were further refined through statistical optimization to identify key parameters that influence the release and effectiveness of the active compounds. Complementing the experimental approach, a bioinformatics analysis was conducted to explore the molecular interactions between limonene, the primary active component of EOCL, and target proteins in insects. This theoretical investigation provided insights into the potential mechanisms of action, reinforcing the empirical findings. The study concludes that encapsulating EOCL within porous clay matrices not only enhances the stability and controlled release of the oil but also significantly boosts its insecticidal effectiveness. This approach presents a promising, environmentally sustainable strategy for crop protection, integrating material science, theoretical modeling, and bioinformatics to develop more efficient and durable biopesticides.
Keywords: 
Subject: 
Biology and Life Sciences  -   Agricultural Science and Agronomy

1. Introduction

The preservation of stored grains is a critical global challenge, essential for ensuring food security amid increasing demand. Pests like Sitophilus granarius (granary weevil) pose significant threats by infesting and degrading grain reserves, leading to both quantity and quality losses. Traditional management of these pests has relied heavily on synthetic chemical insecticides, which, despite their effectiveness, have led to issues such as pest resistance, health risks from chemical residues, and environmental damage [1,2,3,4,5]. Consequently, there is a growing need for safer, more sustainable pest control alternatives. Essential oils, particularly lemon essential oil, have emerged as promising candidates due to their insecticidal properties, low toxicity, and environmental safety [6]. Lemon essential oil, rich in limonene, is effective against Sitophilus granarius and is biodegradable, reducing long-term ecological risks [7,8,9]. This study aims to investigate the potential of lemon essential oil as a natural insecticide, focusing on its interaction with porous clay supports to enhance its residual activity and effectiveness. The study will evaluate the oil’s efficacy under various conditions using advanced statistical and modeling techniques to optimize its application in real-world storage environments [10]. The research also explores the use of natural porous clays from Morocco to encapsulate the oil, improving its stability and controlled release, which are critical for maximizing insecticidal efficacy [11,12,13]. The clays will be characterized using advanced analytical techniques, and the oil's diffusion through the clay matrix will be modeled to better understand and optimize its release [14,15]. This integrated approach, combining material characterization, theoretical modeling, and bioinformatics analysis, could lead to innovative pest control solutions and a deeper understanding of the mechanisms involved in the interactions between essential oils and porous matrices, paving the way for new strategies in natural product-based insect control.

2. Materials and Methods

2.1. Materials

2.1.1. Porous Clay

The porous clay used in this study was sourced from the Bejaad region of Morocco (32°47'01.7"N 6°13'52.3"W) and characterized using several advanced techniques. X-ray diffraction (XRD) was performed with a Bruker CCD-Apex device, while the oxide content was analyzed using a Siemens SRX 3000 X-ray fluorescence spectrometer (XRF) [16]. Fourier-transform infrared spectroscopy (FT-IR) was conducted using a Bruker Alpha spectrometer [17], and the morphology was examined with a scanning electron microscope (SEM-EDX) [18]. Thermogravimetric analysis (TGA) was performed using a Du Pont analyzer [19].

2.1.2. Lemon Essential Oil

Lemon essential oil (Citrus limonum) was extracted through hydrodistillation and analyzed using gas chromatography-mass spectrometry (GC-MS) with a Hewlett Packard 5971A instrument [20]. The oil's molecular profile was further determined using gas chromatography with flame ionization detection (GC-FID) [21].

2.1.2. Preparation of the Clay/Essential Oil Mixture

The RC and GC clays were prepared by grinding and heating, then mixed with lemon essential oil in various concentrations. The mixture was placed in metal cylinders for further experimentation [22].

2.2. Theoretical Considerations

2.2.1. Diffusion Model

The desorption of lemon essential oil within the porous clay matrix can be modeled using Fick second law of diffusion. This law describes how the concentration of a diffusing substance changes over time and space within a medium. For this study, the general form of Fick second law, applied to a one-dimensional system, can be expressed as:
∂ C x , t ∂ t = D ∂ 2 C x , t ∂ x 2
where:
- C x , t is the concentration of the essential oil at position x and time t,
- D is the diffusion coefficient of the essential oil in the clay matrix,
- ∂ C x , t ∂ t is the rate of change of concentration with time,
- ∂ 2 C x , t ∂ x 2 is the second derivative of concentration with respect to position, representing the spatial variation in concentration.
In this model, the diffusion coefficient D is a critical parameter that depends on the physical properties of the lemon essential oil, the characteristics of the clay, and the interaction between the oil and the clay matrix [23]. The boundary conditions and initial conditions are defined based on the experimental setup, where the concentration of the oil at the surface of the clay and within the matrix at the initial time point is considered.
To accurately describe the desorption behavior, the model must account for the specific properties of lemon essential oil, such as its volatility, viscosity, and interaction with the clay pore structure. The heterogeneous nature of the porous clay, including its varying pore size distribution and surface area, may affect the diffusion process, potentially leading to non-uniform desorption rates [24].

2.2.2. Advanced Models

Non-linear diffusion models were considered to account for complex interactions and external factors like temperature and humidity, providing a more accurate representation of the oil's desorption behavior [25,26].

2.2.3. Assumptions

Key assumptions included one-dimensional diffusion, constant diffusivity, and uniform initial concentration within the clay matrix [27].

2.2.4. Model of Mass Transfer by Diffusion

The diffusion process was further described using Fick's second law in cylindrical coordinates, with simplified equations to model the mass transfer of the essential oil [28].

2.2.5. Mathematical Framework

Fick second law for diffusion in cylindrical coordinates (r, θ, z) is expressed as:
∂ C ∂ t = 1 r ∂ ∂ r r D r ∂ C ∂ r + 1 r ∂ ∂ θ D θ r ∂ C ∂ θ + ∂ ∂ z D z ∂ C ∂ z
where:
- C r , θ , z , t is the concentration of EOCL at a position r , θ , z and time t,- D r ,   D θ ,   D z are the diffusion coefficients in the r ,   θ ,   a n d   z directions, respectively.
Given that the height of the cylinder is much smaller relative to its diameter and lateral surface, the diffusion can be assumed to primarily occur in the axial direction (along z). This simplifies the model to one-dimensional diffusion, reducing Fick second law to:
∂ C ∂ t = D z ∂ 2 C ∂ z 2

2.2.6. Boundary and Initial Conditions

The solution to the above equation requires appropriate boundary and initial conditions:
Initial Condition: At t = 0 , the concentration of EOCL is assumed to be uniform throughout the clay:
C z , 0 = C 0 for   all   0 ≤ z ≤ h
Boundary Conditions:
- For t > 0 at the surface z = h, the concentration of EOCL at the boundary becomes:
C z = h , t = C t
where Ct is the concentration at the boundary at time t.

2.2.7. Analytical Solution

The analytical solution to the simplified diffusion equation has been derived by Crank (1979) for these boundary conditions and is given by:
C z , t = 4 π ∑ n = 0 ∞ exp ⁡ − D z 2 n + 1 2 π 2 t h 2 cos ⁡ 2 n + 1 π z h
The total mass Mt of EOCL in the clay RC at any instant t is obtained by integrating the concentration over the thickness of the material and the surface S:
M t = ∫ 0 h S · C z , t d z
Under the assumption that the desorption process tends towards equilibrium over an infinite desorption time, the analytical solution for the mass at any time t relative to the initial mass M₀ is given by:
M t − M ∞ M 0 − M ∞ = 8 π 2 ∑ n = 0 ∞ 1 2 n + 1 2 e x p − Dz 2 n + 1 2 π 2 t h 2 ⁡
where:- M0 is the initial mass of EOCL absorbed,
- M∞ is the mass of EOCL at equilibrium.
For practical purposes, the equation can be simplified by considering only the first term of the series and assuming M∞ = 0:
M t M 0 = 8 π 2 exp ⁡ − D z π 2 t h 2
Given that (8/π²) ≈ 1, the equation further simplifies to:
M t M 0 = exp ⁡ − D z π 2 t h 2
The diffusivity Dz can be determined from the slope of a linear or logarithmic extrapolation of the simplified equation. The evaporation constant K can be approximated by :
K = π 2 D z h 2
The initial evaporation rate F₀ is obtained from the initial slope of the desorption curve as a function of time:
F 0 = 1 S lim t → 0 ⁡ d M d t
These parameters, diffusivity Dz, evaporation constant K, and evaporation rate F are critical in describing the desorption process and are central to the modeling and understanding of EOCL behavior in the porous clays RC and GC [29].

2.3. Experimentation and Modeling - Essentials

2.3.1. Insecticidal Activity Test

M % = M t e s t − M c o n t r o l M c o n t r o l × 100
where:
M % : is the corrected mortality,
M t e s t : is the observed mortality during the test,
M c o n t r o l :​ is the mortality observed in the control.
Lethal Dose ( L D 50 ) Determination: The lethal dose required to kill 50% of the insect population L D 50 was determined by linear interpolation of curves that plot the percentage of mortality as a function of the logarithm of the tested concentration [34]. This provides a quantitative measure of the toxicity of EOCL, EOCL + RC, and EOCL + GC under the conditions tested.
This procedure allowed for a detailed assessment of EOCL insecticidal efficacy, both alone and when mixed with RC and GC clays, under controlled conditions, providing valuable data for optimizing its use in pest control.

2.3.2. Optimization of Transfer Conditions

To optimize the experimental conditions for the diffusion transfer of EOCL (lemon essential oil) in both RC and GC porous clay media, a statistical approach based on experimental designs was employed. The transfer process was evaluated by measuring the evaporation rates, denoted as F.
A full factorial design was used to structure the experiments, considering four factors ( X i ​):
Factor 1 = C: Concentration of essential oil (0.01 mL/cm³ and 0.02 mL/cm³).
Factor 2 = T: Temperature (25 °C and 30 °C).
Factor 3 = D: Cylinder diameter (1 cm and 2 cm).
Factor 4 = M: Mass of clay rock (0.05 g and 0.10 g).
The total number of tests required was calculated using the formula:
Number of tests =   2 k where 2 represents the two levels (-1 and +1) and k is the number of factors studied. In this study, with k = 4 k Number of tests= 2 4 =16
This results in 16 tests being conducted for each type of clay (RC and GC). Table 1 presents the design matrix of the experiments and the corresponding evaporation rates obtained for both RC and GC clays. Each row in the Table 1 represents a specific test, with the factors coded according to their high (+1) or low (-1) levels [35].
The polynomial model used to describe the relationship between the factors and the evaporation rates F for each clay type is expressed as:
F = b 0 + ∑ i = 1 n b i X I + ∑ i = 1 n ∑ j = 1 n − 1 b i j X i X j + ∑ i = 1 n ∑ j = 1 n − 1 ∑ k = 1 n − 2 b i j k X i X j X k + b i j k X i X j X k X l
where:
b 0 ​: is the mean,
b i   : represents the main effects (4 main effects for the factors),
b i j : represents the second-order interaction effects (6 interactions),
b i j k : represents the third-order interaction effects (4 interactions),
b i j k l ​: represents the fourth-order interaction effect (1 interaction).
This model allows for the analysis of how each factor, and their interactions, influence the evaporation rates, enabling the optimization of the transfer conditions for EOCL in both RC and GC clays [36].

2.3.3. Study of Insecticidal Activity Behavior of Essential Oil Mixture in Clay Porous Media

The interaction between the behavior of EOCL in the porous clay media (RC and GC) and its insecticidal activity was studied by analyzing the relationship between the parameters of mass transfer by diffusion and the mortality rate in insecticidal activity. This analysis was conducted using Principal Component Analysis (PCA), a powerful statistical tool for compressing and synthesizing information. PCA is particularly useful when dealing with a large amount of quantitative data, as it helps to interpret and identify correlations between different variables, providing insights into the effectiveness of EOCL in different clay matrices [37,38].
Table 2. LD50 ​ of the Insecticidal Activities of EOCL Fixed on Porous Media (RC and GC).
Table 2. LD50 ​ of the Insecticidal Activities of EOCL Fixed on Porous Media (RC and GC).
Diameter (cm) Temperature (°C) LD50 (a) for RC (b) for RC LD50 (a) for GC (b) for GC
D1 = 1 T1 = 25 0.228 0.110 0.229 0.115
T2 = 30 0.350 0.109 0.349 0.111
T3 = 35 0.035 0.017 0.036 0.018
D2 = 2 T1 = 25 0.529 0.032 0.530 0.034
T2 = 30 0.758 0.008 0.759 0.009
T3 = 35 0.021 0.007 0.022 0.007
D3 = 3 T1 = 25 ≫ 0.056 0.016 ≫ 0.057 0.017
T2 = 30 0.061 0.017 0.062 0.018
T3 = 35 0.013 0.005 0.014 0.006
≫: LD50 is greater than 1 mL/cm3. a: EOCL alone. b: EOCL + (RC or GC) clay.

2.4. The Bioinformatics Approach

Computational research was integrated as an additional step in this study to provide a comprehensive theoretical understanding of limonene, the primary molecule found in Citrus limonum, and its interaction with key proteins 2NXX and 1C3Z. These interactions offer unique insights into the biological mechanisms underlying development regulation and survival under extreme conditions, while also opening new avenues for protein docking applications [39]. Two distinct approaches were employed for this analysis. The first approach involved modeling the molecular interactions of limonene, with a particular focus on determining essential energetic descriptors. The second approach involved the docking study of limonene with selected proteins. For molecular modeling, the MMFF94 (Merck Molecular Force Field 94) method was employed, which is a set of force field parameters used in computational chemistry to simulate molecular interactions within the framework of molecular mechanics [40]. Several energetic descriptors were calculated, including the energy levels of the highest occupied (EHomo) and lowest unoccupied (ELumo) molecular orbitals, as well as derived parameters such as the energy gap (EGAP), chemical hardness (η), electronegativity (χ), electrophilicity index (ω), and molecular flexibility (s). The formulas used to calculate these descriptors are listed in Table 3 [41].
In this study, the proteins encoded by the PDB structures 2NXX and 1C3Z were analyzed to understand their biological functions and interactions with limonene, a predominant molecule in Citrus limonum. These proteins exhibit distinctive structural features that influence their ligand-binding capabilities and present valuable opportunities for protein docking research [42].
The 2NXX protein, identified as the ultraspiracle receptor (USP) from the insect Tribolium castaneum, is part of the nuclear receptor family and is crucial in regulating developmental processes such as molting and metamorphosis. Remarkably, this receptor can maintain a functional conformation without requiring a specific ligand, setting it apart from its vertebrate counterparts. The protein surface, particularly the binding domain, plays a vital role in its molecular interactions. Analysis of the solvent-accessible surface (SAS) revealed the protein unusual ability to interact effectively without a ligand, highlighting its unique properties [43]. The 1C3Z protein, an antifreeze protein (THP12) from the yellow mealworm Tenebrio molitor, is vital for the survival of organisms in cold environments by preventing ice crystal formation and growth. THP12 helical structure and hydrophilic surface allow it to bind effectively to small hydrophobic ligands, enhancing its protective function against freezing. The stability of this protein under extreme conditions is significantly supported by electrostatic interactions on its surface [44].
The surface properties of 2NXX and 1C3Z, including their solvent accessibility and interaction potential with limonene, offer a robust foundation for protein docking studies. The ability of 2NXX to function without a ligand, combined with 1C3Z affinity for hydrophobic ligands such as limonene, facilitates advanced modeling of protein-ligand interactions. These models can be instrumental in designing targeted molecules and exploring the mechanisms of action of these proteins under different biological conditions [45]. Following the optimization of the protein structures using Swiss-PDB Viewer V4.1, several critical steps were undertaken for docking preparation. These included the removal of water molecules and heteroatoms, the addition of polar hydrogens, and the assignment of charges using the Kollman method. The geometry of limonene was optimized with density functional theory (DFT) calculations using the B3LYP/6-311G basis set. Docking simulations were conducted with Autodock VINA 1.1.2 (Table 4), using specific grid dimensions for each protein. The Biovia Discovery Studio (Studio, 2021) software was then utilized to visualize the protein-ligand interactions, allowing for a detailed analysis of how limonene interacts with the target proteins and providing valuable insights into their potential as biological agents [46].

2.5. Computer and Statistical Processing

All computational processing was performed using MATLAB software. The numerical modeling results obtained were compared with the experimental results to identify key parameters and maximize interpretation. In addition to MATLAB, other software tools were employed for various calculations and graphical representations, including EXCEL, ORIGIN LAB, DESIGN PLAN, and XLSTAT. These tools supported the comprehensive analysis and helped ensure that the results were accurately interpreted and effectively communicated. Furthermore, principal component analysis (PCA) has proven to be an indispensable mathematical tool for exploring the correlations between various calculated parameters. The goal was to visualize and synthesize data related to the energetic and qualitative properties of the molecular docking of limonene with the studied proteins [47,48].

3. Results and Analysis

3.1. Chemical Characterization of RC and GC Clay Media

X-ray Diffraction (XRD) Analysis: The XRD analysis revealed the mineralogical complexity of the RC and GC samples, which are primarily composed of phyllosilicates, especially kaolinite and illite. The distinct kaolinite peaks observed at 7.75 Å and 14.5 Å suggest that this mineral formed under conditions of low ionic charge and intense chemical weathering, which are typical of tropical or subtropical climates. Illite, found in both samples, points to diagenetic processes, likely linked to higher temperature and pressure conditions, indicating a complex geological history. The significant presence of quartz, with peaks at 19.8 Å and 26 Å, suggests a detrital input, probably of alluvial or aeolian origin. Specific reflections from carbonates like calcite and dolomite suggest formation in carbonate-rich environments, such as alkaline lakes or shallow marine zones. Finally, the detection of hematite in the RC sample, marked by reflections at 23.8 Å and 29.7 Å, indicates oxidation conditions, probably due to prolonged exposure to air, supporting the hypothesis of a pedogenic origin (Figure 2) [49,50].
X-ray Fluorescence (XRF) Chemical Composition Analysis: The XRF chemical analysis showed that silicon dioxide (SiO2) and aluminum oxide (Al2O3) are predominant in the samples, with notable variations between RC and GC. The GC sample, with its high SiO2 content, reflects a strong presence of silicate minerals such as quartz and silica-rich phyllosilicates, whereas the RC sample, with a more modest SiO2 content, is dominated by iron oxides, primarily Fe2O3, associated with hematite. The higher concentration of Al2O3 in GC suggests an abundance of aluminous clay minerals like kaolinite, while the lower content in RC may indicate a composition richer in non-aluminous minerals. Other oxides present, such as Na2O, K2O, MgO, CaO, and TiO2, are in moderate proportions and reflect the presence of accessory minerals like feldspars, dolomite, and calcite, indicating distinct sedimentary and diagenetic processes (Table 5) [51,52].
Transform Infrared Spectroscopy (FTIR) Analysis: FTIR analysis allowed us to identify the functional groups present in the samples, confirming the presence of clay minerals like kaolinite and illite, as well as carbonates and iron oxides. The absorption bands at 1025 cm⁻¹, associated with Si-O vibrations, reveal a well-ordered crystalline structure typical of aluminosilicates. The band at 920 cm⁻¹, linked to the Al2OH bond, further supports the identification of kaolinite, while the band at 3620 cm⁻¹ indicates the presence of structural water in the clay minerals. The bands at 475 cm⁻¹ and 540 cm⁻¹, representing Si-O and Si-O-Al bonds, confirm the presence of phyllosilicates. The band at 1517 cm⁻¹, corresponding to C-O vibration, indicates the presence of calcite, while the band at 1455 cm⁻¹ in RC, associated with hematite, explains the red color of this sample, revealing geochemical conditions favorable to iron oxidation (Figure 3) [53,54].
Scanning Electron Microscopy (SEM) and Energy Dispersive X-ray (EDX) Analysis: SEM analysis revealed significant morphological differences between the samples. The RC sample exhibits an irregular and microporous structure, suggesting intense weathering, potentially due to wetting and drying cycles or chemical interactions. This morphology is typical of clays subjected to pedogenetic processes. In contrast, the GC sample shows a more regular morphology, suggesting a more stable depositional environment. EDX analysis confirmed a high iron content in RC, mainly in the form of iron oxide (Fe2O3), corroborating the presence of hematite detected by XRD and FTIR. The GC sample is dominated by silicon and aluminum, reinforcing its richness in phyllosilicates like kaolinite and illite. These differences between RC and GC reflect distinct geological and environmental conditions (Table 6) (Figure 4 and Figure 5) [55,56].
Thermogravimetric Analysis (TGA): Thermogravimetric analysis identified three main phases of mass loss in the samples, reflecting dehydration and thermal decomposition processes. The first phase, below 200 °C, corresponds to the loss of surface water and light organic matter. The second phase, between 430 °C and 670 °C, corresponds to the dehydroxylation of clay minerals, particularly kaolinite, forming metakaolinite. This phase indicates a significant kaolinite content in both samples. The third phase, between 700 °C and 800 °C, corresponds to the decomposition of illites and carbonates like calcite and dolomite. The significant mass loss at this stage confirms the presence of these minerals, suggesting formation conditions that favor their stable precipitation. These thermogravimetric results reveal the complexity of the samples and their thermal behavior, providing insights into the geochemical conditions that influenced their formation (Figure 6) [57,58].

3.2. Chemical characterization of essential oil of Citrus limonum:

The analysis of Citrus limonum essential oil, conducted using gas chromatography coupled with mass spectrometry, reveals a detailed composition of the compounds present, each identified by its retention time (RT), retention index (RI), and area percentage (Area %), indicating their relative proportion in the sample. Limonene, which dominates with 51.41% of the total composition, is a major component typical of lemon essential oils. This monoterpene is not only responsible for the characteristic citrus aroma but is also recognized for its insecticidal properties, as well as its antioxidant, anti-inflammatory, and antimicrobial activities [59]. Additionally, β-Pinene (6.46%) and α-Pinene (4.84%), also common monoterpenes in citrus oils, contribute to the insecticidal efficacy and the aromatic complexity of the oil. γ-Terpinene (4.60%) and α-Terpineol (4.13%), present in notable amounts, also play a role in the biological activities, including antimicrobial and soothing effects [60]. The analysis also reveals several minor compounds, each representing less than 1% of the total, such as cis-Citral (1.27%), trans-Caryophyllene (1.06%), and cis-α-Bergamotene (1.97%). Although present in smaller quantities, these compounds can enhance the biological activity of the oil, particularly in the antimicrobial and insecticidal domains, while adding to the aromatic complexity [61]. Citrus limonum essential oil is predominantly composed of monoterpenes, notably limonene, β-Pinene, and α-Pinene, which are typical of citrus oils and responsible for its distinctive lemony aroma. Sesquiterpenes, such as β-Elemene and trans-Caryophyllene, though present in smaller amounts, can also play an important role in the overall aroma and therapeutic properties, particularly in enhancing insecticidal efficacy (Table 7) [62].
The analysis of Citrus limonum essential oil by gas chromatography coupled with mass spectrometry (GC-MS) shows significant variability in the chemical composition depending on the geographical origin of the sample. In several studies, limonene is consistently identified as the major compound, but its proportion varies considerably. Najwa Nasser AL-Jabri conducted a comparative analysis of lemon essential oils from Turkey and India. Limonene was the main compound in both varieties, representing 78.93% in the Turkish essential oil, followed by 5.08% β-pinene, while in the Indian variety, limonene represented 53.57% and α-terpineol reached 15.15% [59]. Similarly, Njoku et al. analyzed the essential oil of lemon from Nigeria, finding 54.15% limonene, followed by 6.21% 2-cyclohexan-1-ol and 3.82% β-pinen [60]. These variations indicate that the chemical composition of Citrus limonum essential oil is influenced by climate, soil, and other environmental factors specific to each region. In comparison, Bertuzzi et al. (2013) reported that lemon essential oil contains 71.87% limonene as the main constituent [61]. In another study, Himd et al. (2016) identified a limonene content of 64.19% in the essential oil of Lisbon lemon peel (Portugal), followed by β-pinene (7.76%) and α-terpinene (5.45%) [62]. However, a study by Moufida & Marzouk (2003) on lemon essential oil shows that geranial, at 3.65%, is the main compound, which is quite different from other studies [63]. These differences highlight the impact of local conditions and processing methods on the composition of Citrus limonum essential oils, emphasizing the importance of precise chemical characterization for each sample based on its geographical origin.

3.3. Insecticidal Activities in Porous Clay Media

The results of the insecticidal tests for EOCL, both alone and in combination with the porous clay media RC and GC, against Sitophilus granarius are summarized in Table 2, which presents the lethal doses of 50% LD50 as a function of the studied parameters, such as cylinder diameter (D) and incubation temperature (T).
From these results, it is evident that the LD50 ​ values decrease when EOCL is fixed on the porous clay media RC and GC, indicating an enhancement in insecticidal efficacy [64]. This suggests that the fixation of EOCL on the porous media alters its insecticidal potency. Additionally, the factors of cylinder diameter and temperature continue to play significant roles in maintaining the insecticidal properties, with larger cylinder diameters and higher temperatures generally promoting more effective insecticidal activity [65]. Furthermore, the analysis indicates that the insecticidal activity of EOCL when fixed on porous media, as well as its persistence, depends on the nature of the interactions between the essential oil compounds and the porous media. The specific interactions between the EOCL and the clay components likely contribute to the observed changes in insecticidal performance, underlining the importance of the media in modulating the efficacy of EOCL. This conclusion is consistent for both RC and GC clay types, suggesting a similar mechanism of action for the EOCL when used with these different porous media (Table 8) [66].
The experimental results reveal that the desorption of lemon essential oil in porous clays is significantly influenced by factors such as oil concentration, temperature, clay particle diameter, and the mass of clay used. These parameters determine the amount of essential oil available to act as an insecticide against Sitophilus granarius. Temperature plays a crucial role in this process, with desorption accelerating as the temperature increases, particularly between 25°C and 35°C [67]. For example, red clay exhibits faster mass loss at 35°C, indicating enhanced diffusion and potentially improved insecticidal efficacy at higher temperatures. Green clay demonstrates a slightly faster desorption rate under the same conditions, with desorption rates increasing from 30% to 60%, compared to 25% to 55% for red clay [68]. The concentration of the essential oil also affects desorption. Higher concentrations tend to slow initial desorption, likely due to increased clay saturation, which extends the insecticidal effect as more oil remains available over a longer period. For instance, at 25°C and with a particle diameter of 1 mm, the residual mass is higher for a stronger concentration, indicating slower but more sustained desorption. For red clay, the residual mass decreased from 75% to 45% as the temperature increased, while for green clay, it decreased from 70% to 40% [69].
Clay particle diameter influences desorption as well. Smaller particles, due to their higher specific surface area, promote increased interaction with the essential oil, leading to faster desorption and a more immediate insecticidal effect. The results show that green clay with finer particles has a higher desorption rate than red clay, highlighting the importance of particle size [70]. Lastly, a larger mass of clay leads to slower desorption of the essential oil, as more clay can adsorb and gradually release the oil, thus prolonging the insecticidal effect. This effect is more pronounced in red clay than in green clay, suggesting that red clay may be more effective for applications requiring prolonged release [71]. To maximize the insecticidal efficacy of lemon essential oil against Sitophilus granarius, it is crucial to consider the interaction between temperature, oil concentration, clay particle size, and clay mass. Faster desorption, facilitated by higher temperatures and finer particles, could provide a more immediate insecticidal effect, while slower desorption associated with higher concentration and greater clay mass could prolong the duration of the effect [72]. The analysis of desorption data in red and green clays allows for the identification of optimal conditions to maximize insecticidal activity, which is essential for optimizing the storage and application conditions of lemon essential oil. Regarding the estimation of LD50, it was determined that reliable values require well-distributed data covering an adequate mortality range. Through adjustments to the logistic model parameters and linear interpolation, the LD50 for lemon essential oil against Sitophilus granarius was estimated to be 0.02 µl/cm³, highlighting the importance of proper data distribution for accurate analysis (Table 9) (Figure 7) [73].

3.4. Optimization of Essential Oil Transfer Conditions in Porous Media

The optimization of the evaporation flux parameter for essential oil transfer in porous media was conducted using a full factorial design involving four factors: concentration of essential oil (C), temperature (T), cylinder diameter (D), and mass of clay medium (M). Each factor was tested at two levels, coded as +1 for the high level and -1 for the low level, resulting in a total of 16 simulations. The mathematical model associated with this first-degree polynomial accounts for the main effects of these factors, as well as their second, third, and fourth-order interactions [74]. The coefficients for the flux model of essential oil in the red clay (RC) medium indicate that the mean flux achieved was 2.6061×10-4 Similarly, the coefficients for the flux model of essential oil in the green clay (GC) medium provide insight into the behavior of the oil in this different porous medium. In both RC and GC, an increase in the concentration of essential oil enhances the flux, while a rise in temperature also contributes to an increased flow [75]. Conversely, increasing the cylinder diameter and the mass of the clay medium results in a reduction of the flux in both types of clay. Graphical representations of the flux as a function of two interacting factors in 3D for both RC and GC demonstrate that the interaction between concentration and temperature significantly increases the flux. However, interactions between concentration and diameter, concentration and mass, mass and temperature, and temperature and diameter generally result in a decrease in flux. The interaction between diameter and mass also leads to a reduction in the flow rate for both red and green clays [76]. Overall, the determination coefficients R2 for the models developed for both RC and GC media indicate that the optimized models are well-explained by the data. The Fisher test value (Ftest = 2.42) suggests that the model is not significantly influenced by noise or uncertainties. However, the predictive coefficient of determination (R2pred ​) yielded a negative value, implying that the total mean flux might be a better predictor of the model response [77]. Despite this, the response adequacy, measured by the flux-to-uncertainty ratio, was found to be 6.68, exceeding the desirable threshold of 4, which confirms the adequacy of the response. This analysis underscores the complex interactions between various factors affecting the flux of essential oil in porous media, particularly in red and green clays, and highlights the importance of optimizing these conditions to ensure efficient transfer of essential oil (Figure 8 and Figure 9) [78].

3.5. Modeling the Diffusion Process of Citrus Limonum Essential Oil (EOCL) in Porous Media RC and GC

To better understand how Citrus Limonum essential oil (EOCL) diffuses in porous media like RC and GC, two complementary approaches were taken: studying the kinetics of desorption and simulating concentration profiles within a cylindrical model. These steps are crucial for identifying key parameters such as the diffusion coefficient Dz, evaporation flux F, evaporation constant K, and activation energy Ea, all of which are essential for accurately describing the oil behavior in these media. The Fickian diffusion model was chosen to analyze desorption, based on the mass percentage of oil over time [79]. For the RC medium, the diffusion coefficient Dz was estimated at 9.81×10-6 cm²/h. under optimal conditions (C = 0.01 mL/cm³, T = 30°C, M = 0.10 g, D = 2 cm). The results showed a significant decrease in oil mass, confirming the validity of the model used [80]. This simulation is not only precise but also useful for predicting the remaining amount of oil over time, which is crucial for applications such as insecticidal efficacy. Studies conducted on the GC medium produced similar results, though some adjustments are necessary due to the specific properties of this medium [81]. Overall, the Fickian diffusion model has proven to be an effective tool for understanding and predicting EOCL diffusion, providing a solid foundation for optimizing the use of essential oils in various practical applications (Figure 10 and Figure 11), [82].

3.6. Interaction of Essential Oils in Porous Media and Insecticidal Efficacy

The study of the behavior of Citrus limonum essential oil (EOCL) in porous media, particularly in red clay (RC) and green clay (GC), is based on several key parameters: evaporation flux (F), diffusivity (Dz), evaporation constant (K), activation energy (Ea), and insecticidal activity, measured by mortality percentage (M%). These parameters, detailed in Table 8 and Table 9, provide a complex but valuable insight into how EOCL interacts with these porous media. The activation energy values, while consistent and low across both media, suggest similar energetic behavior between RC and GC. The observed negative values indicate a dominant adsorption process, highlighting the strong capacity of these clays to retain essential oils. This retention capability is crucial for understanding the desorption dynamics and the insecticidal efficacy of EOCL [83]. To further analyze these interactions, Principal Component Analysis (PCA) was conducted, focusing on evaporation flux (F), diffusivity (Dz), evaporation constant (K), and mortality percentage (M%) across 16 different cases. The PCA identified two principal components (F1 and F2), which together explained 82% of the total variability, with eigenvalues of 1.89 and 1.38, respectively [84]. The mapping of these components revealed that the diffusion coefficient (Dz) and the evaporation constant (K) are closely related, as indicated by their proximity on the F1 axis. This correlation suggests a dependency between these two parameters, aligning with theoretical expectations. On the F2 axis, the evaporation flux (F) and mortality percentage (M%) also show a strong correlation, implying that insecticidal efficacy, as measured by mortality, is largely influenced by the rate of EOCL evaporation. Interestingly, the near-orthogonal relationship between the pairs (Dz, K) and (F, M%) suggests independence between diffusion-related parameters and evaporation flux-related parameters [85]. This distinction underscores the complexity of the interactions between the physical diffusion processes within the clay matrices and the resulting biological insecticidal activity. The PCA not only clarifies the interactions between EOCL and the clay media but also highlights the key factors influencing its insecticidal effectiveness. Understanding these relationships is essential for optimizing the use of essential oils in pest control applications, where the choice of porous medium can significantly impact the efficacy of the treatment (Table 10, Table 11 and Table 12), (Figure 12) [86].
The results provide a detailed analysis of the correlation between several key parameters involved in the desorption of EOCL essential oil and its insecticidal activity in two types of clay media: red clay (RC) and green clay (GC). The parameters examined include evaporation flux (F), diffusivity (Dz), evaporation constant (K), and mortality percentage (M%). In the RC medium, the F1 axis, representing 47.31% of the total variability, is the most significant in explaining the desorption behavior of EOCL. A strong correlation between diffusivity (Dz) and the evaporation constant (K) is observed, suggesting that changes in diffusivity directly influence the evaporation constant. This highlights the crucial role of diffusion in the desorption dynamics within this medium [87]. Simultaneously, the F2 axis, accounting for 34.70% of the variability, reveals a close relationship between the evaporation flux (F) and the mortality percentage (M%). This indicates that the evaporation rate of EOCL has a direct impact on its insecticidal efficacy, with faster evaporation being associated with higher mortality rates. The orthogonal relationship between the pairs (Dz, K) and (F, M%) suggests that diffusion processes are independent of evaporation flux and the resulting insecticidal activity, revealing a complex interplay between physical diffusion processes and biological effectiveness [88]. In the GC medium, the F1 axis similarly accounts for 47.31% of the total variability and shows a strong correlation between diffusivity (Dz) and the evaporation constant (K). This underscores the critical role of diffusion in shaping the evaporation characteristics of EOCL in this medium. Furthermore, the F2 axis, which explains 34.70% of the variability, shows a significant correlation between the evaporation flux (F) and the mortality percentage (M%), echoing the findings in the RC medium. This highlights that the evaporation rate of EOCL is a key factor in its insecticidal efficacy. The orthogonal relationship between the pairs (Dz, K) and (F, M%) in the GC medium, as in the RC medium, underscores the independence of diffusion processes from the factors controlling evaporation rate and insecticidal activity, emphasizing the complexity of these interactions in green clay [89]. Overall, in both RC and GC media, the strong correlation between diffusivity (Dz) and the evaporation constant (K) underscores the critical importance of diffusion in the EOCL desorption process. The evaporation flux (F) emerges as a key determinant of insecticidal activity (M%), with faster evaporation improving insecticidal efficacy. The independence between diffusion-related and evaporation-related parameters indicates that these factors operate through distinct mechanisms. Therefore, it is essential to consider them separately when optimizing EOCL desorption and insecticidal effectiveness [90].

3.7. Bioinformatic Studies

The bioinformatics study relies on two distinct computational approaches focused on analyzing the limonene molecule and its interaction with two key proteins, 2NXX and 1C3Z. The first approach involves the energetic analysis of limonene, while the second explores its molecular docking with these proteins, which are crucial for understanding biological processes such as development and survival under extreme conditions. The results of the energetic analysis are presented in (Figure 13) and (Table 13). The Figure 13 displays the energy level diagram of the molecular orbitals of limonene, optimized using the MMFF94 method, while Table 13 summarizes additional calculations and relevant energetic parameters. These data reveal that limonene exhibits significant electronic stability, characterized by an energy gap (EGAP) of 6.677 eV and a chemical hardness of 3.338 eV, indicating a strong resistance to reactivity [91]. An electronegativity of 7.634 eV suggests a robust electron-attracting ability, contributing to its stability. The electrophilicity index of 8.729 eV indicates a moderate capacity of the molecule to accept electrons during chemical reactions, while a molecular flexibility value of 0.149 eV points to moderate reactivity. In summary, limonene demonstrates high electronic stability, low reactivity, and a moderate ability to act as an electrophile in chemical reactions. These results are crucial for understanding the electronic behavior and potential reactivity of limonene, providing a solid foundation for further investigation of its interactions with the 2NXX and 1C3Z proteins. These interactions are particularly relevant given the biological roles of these proteins, and the findings from the molecular docking studies will shed light on the potential applications of limonene in biological and therapeutic contexts [92].
The docking results of limonene with the proteins 2NXX and 1C3Z reveal subtle yet significant differences in the interaction between the ligand and these two proteins. The binding free energy of -6.5 kcal/mol for 2NXX indicates a slightly more stable interaction with limonene compared to the binding free energy of -6.3 kcal/mol observed with 1C3Z. Although this difference is modest, it suggests that the binding site of 2NXX is better suited to accommodate limonene, providing a more favorable environment for interaction [93]. Furthermore, the pKi of 4.77 for 2NXX, slightly higher than the 4.62 for 1C3Z, confirms a stronger binding affinity of limonene for 2NXX. This higher pKi reflects greater stability of the formed complex, reinforcing the idea that 2NXX is a more suitable partner for limonene [94]. The ligand efficiency, measured at 0.65 kcal/mol/non-H atom for 2NXX compared to 0.63 for 1C3Z, also indicates more efficient binding of limonene to 2NXX. This result could be attributed to better steric compatibility or a more favorable alignment of binding surfaces in 2NXX. The torsional energy, which is identical for both complexes at 0.3113, suggests that the conformational flexibility of limonene remains constant regardless of which protein it binds to. This indicates that the energetic cost for limonene to adapt to the binding site conformation is similar in both cases, highlighting good conformational compatibility with both proteins [95]. In summary, while both 2NXX and 1C3Z offer binding sites that are well-suited to limonene, the data indicate a slight preference for 2NXX. This preference is supported by a more negative binding free energy, a higher pKi, and a slightly superior ligand efficiency. These results suggest that 2NXX may provide a more optimal binding site for limonene, which could have implications for its biological function in the presence of similar compounds. However, the close similarity of the results for both proteins indicates that limonene is an effective ligand for each, with a slightly stronger affinity for 2NXX (Figure 14), (Table 14) [96].
The principal component analysis (PCA) presented in the graph highlights the complex relationships between docking parameters and the studied proteins. The F1 axis, as the first principal component, accounts for a significant portion of the variance in the data. Parameters that deviate most from the origin along this axis contribute the most to the overall variability of the results. In contrast, the F2 axis, the second principal component, explains a smaller portion of the variance, as indicated by the annotation "F2 (0.00%)." Parameters near this axis show little variation compared to F1. In this analysis, Binding Free Energy and Ligand Efficiency stand out with thicker arrows, emphasizing their importance. Their proximity suggests a positive correlation: an increase in binding affinity is generally associated with better ligand efficiency [97]. On the other hand, pKi, represented by an arrow pointing in a different direction, appears to capture a distinct aspect of ligand-protein interactions, likely related to the strength of inhibition or the ligand affinity for the protein. The vertical orientation of Binding Free Energy and Ligand Efficiency indicates that these parameters are primarily influenced by F1, confirming the minimal contribution of F2 [98]. The black dots, representing the proteins, are positioned based on their coordinates in the principal component space. Their proximity to certain arrows reveals which proteins are most influenced by these parameters. A protein located near Binding Free Energy or Ligand Efficiency might indicate that the interaction between the ligand and this protein is strongly determined by these parameters, whereas a protein far from these arrows could be influenced by other factors. These findings suggest that certain proteins exhibit varying binding affinities or efficiencies depending on the characteristics of the ligand. This understanding is crucial for drug development and for deepening knowledge of the biological mechanisms underlying molecular interactions [99]. The correlation circle provides a visual representation of these correlations: arrows reaching the edge of the circle are strongly correlated with the F1 or F2 axes, while those near the center show weaker correlations. The balanced distribution of parameters around the circle indicates that they capture diverse dimensions of ligand-protein interactions, thereby enriching the overall analysis. This PCA offers a clear and in-depth visualization of the complex relationships between multiple docking parameters and the studied proteins, providing valuable insights into the molecular mechanisms underlying these interactions (Figure 15) [100].

4. Conclusions

The study highlights the promising effectiveness of encapsulating lemon essential oil (Citrus limonum) within natural porous clay matrices to combat the pest Sitophilus granarius, a significant concern for the preservation of stored grains. The research demonstrated that using porous clays, sourced from the Bejaad region in Morocco, as carriers for the essential oil helps to overcome the limitations related to the volatility and rapid degradation of these active compounds while optimizing their insecticidal potential. Analyses revealed that the fixation of the essential oil on clay matrices significantly enhances its insecticidal efficacy, reducing the lethal dose required to achieve the desired effect and allowing for a more controlled and prolonged release of the active compound. This increased effectiveness is attributed to specific interactions between the essential oil and the mineralogical components of the clays, which favorably influence the diffusion and persistence of limonene, the main bioactive compound in the essential oil. Furthermore, mathematical modeling of the diffusion process, combined with rigorous statistical approaches, allowed for a detailed understanding of the mechanisms of essential oil release from the clay matrices. These analyses showed that parameters such as temperature, oil concentration, clay particle size, and mass play a crucial role in the desorption rate and, consequently, in the insecticidal efficacy. A fine understanding of these mechanisms enables the adjustment of application conditions to maximize effectiveness while minimizing risks to the environment and human health. The study also opened up new perspectives through bioinformatic analysis of molecular interactions between limonene and certain target proteins in insects, suggesting that these interactions could enhance the effectiveness of the essential oil as a natural biopesticide. This paves the way for innovative and sustainable strategies for the development of natural biopesticides, offering viable alternatives to traditional chemical pesticides. In conclusion, this work makes a significant contribution to the understanding and optimization of essential oil use in biological control applications. It paves the way for innovative strategies for the development of natural biopesticides, offering sustainable alternatives to conventional chemical pesticides while maintaining the necessary effectiveness for crop protection. This integrated approach, combining material characterization, theoretical modeling, and bioinformatic analyses, could serve as a model for future research in the field of biological control using essential oils.

References

  1. Sabry, A.H.; Abdel-Aziz, M.I. Potential of Some Insecticides for Managing the Granary Weevil, Sitophilus granarius, and the Development of Resistance. Egypt J Biol Pest Control. 2020, 30, 50. [Google Scholar] [CrossRef]
  2. Gambacorta, G.; Faccia, M.; Pati, S.; Lamacchia, C.; La Gatta, B. Effects of Conventional and Alternative Pesticides on Wheat Quality and Safety. Food Res Int. 2020, 131, 108998. [Google Scholar] [CrossRef]
  3. Khammour, F.; Abdoul-Latif, F.M.; Ainane, A.; Mohamed, J.; Ainane, T. Eco-friendly adsorbent from waste of mint: application for the removal of hexavalent chromium. J Chem. 2021, 2021, 8848964. [Google Scholar] [CrossRef]
  4. Cao, J.; Zhang, Y.; Zhang, H.; Fu, Z.; Guo, L. Advances in Encapsulation of Essential Oils into Porous Supports: Application in Pest Control. J Clean Prod. 2022, 346, 131136. [Google Scholar] [CrossRef]
  5. Ouassil, M.; Mohamed Abdoul-Latif, F.; Am, A.; Attahar, W.; Ainane, A.; Ainane, T. Chemical composition of bay laurel and rosemary essential oils from Morocco and their antifungal activity against Fusarium strains. Pharmacologyonline. 2021, 426–433. [Google Scholar]
  6. Hasan, M.M.; Haque, E.; Hossain, M. Insecticidal Properties of Essential Oils Against Grain Storage Pests: Recent Advances. Trends Food Sci Technol. 2021, 109, 60–73. [Google Scholar] [CrossRef]
  7. Mohamed Abdoul-Latif F, Elmi A, Merito A, et al. Chemical analysis of essential oils of Cymbopogon schoenanthus (L.) Spreng. and Nepeta azurea R. Br. ex Benth from Djibouti, in-vitro cytotoxicity against cancer cell lines and antibacterial activities. Appl Sci. 2022, 12, 8699. [Google Scholar] [CrossRef]
  8. Montoya-García CM, Muñoz-Huerta RF, Rueda-Puente E, Pérez-Flores LJ, Rodríguez-Mendoza NM, Luna-Suárez S, et al. Alternative Pest Control in Stored Grain Using Essential Oils: A Global Perspective. Agronomy. 2021, 11, 2075. [Google Scholar] [CrossRef]
  9. López, M.D.; Jordán, M.J.; Pascual-Villalobos, M.J. Toxic Compounds in Essential Oils of Spices and Medicinal Plants: Structural Characteristics and Mechanisms of Action Against Stored Product Insects. J Stored Prod Res. 2021, 93, 101821. [Google Scholar] [CrossRef]
  10. Khalid, S.; Bashir, M.A.; Afzal, M.; Qasim, M.; Afzal, I. Botanical Insecticides: A Green Approach in Stored Grain Pest Management. Agriculture. 2020, 10, 101. [Google Scholar] [CrossRef]
  11. Rouissi, T.; Montibus, M.; Choupeaux-Nicolle, L.; Belzunces, L.P.; Frédéric, L. Genetic and Molecular Basis of Insecticide Resistance in the Granary Weevil, Sitophilus granarius. Insects. 2023, 14, 474. [Google Scholar] [CrossRef]
  12. Islam, M.S.; Hasan, M.M.; Hossain, M.A. Evaluating the Bioactivity of Citrus Essential Oils Against Stored Grain Insects. Foods. 2021, 10, 1062. [Google Scholar] [CrossRef]
  13. Nenaah, G.E. Chemical Composition and Insecticidal Activity of Essential Oils from Four Citrus Species Against Sitophilus granarius. J Pest Sci. 2022, 95, 1129–1142. [Google Scholar] [CrossRef]
  14. Shybat, Z.L.; Mohamed Abdoul-Latif, F.; Mohamed, J.; Ainane, A.; Ainane, T. Antifungal activity of the essential oil of Moroccan myrtle (Myrtus communis L.): application in agriculture. Pharmacologyonline 2021, 485–491. [Google Scholar]
  15. Abdallah, M.; Ali, R.; Emam, H.E. Advanced Biopolymer-Based Carriers for the Controlled Release of Essential Oils as Biopesticides. J Appl Polym Sci. 2024, 141, 52001. [Google Scholar] [CrossRef]
  16. Ainane, A.; Taleb, M.; El-Hajjaji, F.; Hammouti, B.; Chetouani, A.; Ainane, T. Study of dependence between two types of most abundant natural clays in Bejaad province (Central Morocco) using a statistical approach. Moroccan J Chem. 2021, 9. [Google Scholar]
  17. Aliouane, T.; Boudouma, A.; Bensaoula, H.; Trari, M. Structural and Optical Properties of Clay-Based Materials: A Study Using FT-IR and UV-Vis Spectroscopy. Spectrochim Acta A Mol Biomol Spectrosc. 2022, 270, 120796. [Google Scholar] [CrossRef]
  18. Atemni I, Ouafi R, Hjouji K, et al. Extraction and characterization of natural hydroxyapatite derived from animal bones using the thermal treatment process. Emergent Mater. 2023, 6, 551–560. [Google Scholar] [CrossRef]
  19. Morin, C.; Vandewalle, N.; Faidallah, H.; Zarhri, Z. Thermogravimetric Analysis of Natural Clays: Determining the Decomposition Behavior Under Nitrogen Atmosphere. J Therm Anal Calorim. 2021, 146, 83–94. [Google Scholar] [CrossRef]
  20. Mohamed Abdoul-Latif F, Elmi A, Merito A, et al. Essential oils of Ocimum basilicum L. and Ocimum americanum L. from Djibouti: chemical composition, antimicrobial and cytotoxicity evaluations. Processes. 2022, 10, 1785. [Google Scholar] [CrossRef]
  21. Ainane A, Benhima R, Khammour F, et al. Composition chimique et activité insecticide de cinq huiles essentielles: Cedrus atlantica, Citrus limonum, Eucalyptus globules, Rosmarinus officinalis et Syzygium aromaticum. Proc Biosune. 2018, 1, 67–79. [Google Scholar]
  22. El Koubali, M.; Baddi, G.A.; Rhihil, A.; Aouad, A.; Kachkoul, R. Preparation and Characterization of Moroccan Natural Clays with Essential Oils for Potential Insecticidal Applications. Appl Clay Sci. 2024, 238, 106006. [Google Scholar] [CrossRef]
  23. Yu, B.; Zhang, J.; Wu, J.; Zhang, P.; Wu, L. Fickian and Non-Fickian Diffusion of Essential Oils in Porous Media: Insights from Experiments and Simulations. Ind Eng Chem Res. 2022, 61, 1735–1747. [Google Scholar] [CrossRef]
  24. Ainane A, Abdoul-Latif FM, Mohamed J, et al. Behaviour desorption study of the essential oil of Cedrus atlantica in a porous clay versus insecticidal activity against Sitophilus granarius: explanation of the phenomenon by statistical studies. Int J Metrol Qual Eng. 2021, 12, 12. [Google Scholar] [CrossRef]
  25. Ainane T, Mohamed Abdoul-Latif F, Am A, et al. Antagonistic antifungal activities of Mentha suaveolens and Artemisia absinthium essential oils from Morocco. Antagonistic antifungal activities of Mentha suaveolens and Artemisia absinthium essential oils from Morocco. Pharmacologyonline. 2021, 470–478. [Google Scholar]
  26. Gao, X.; Xu, Z.; Zhang, L.; Li, W. Nonlinear Diffusion and Sorption Kinetics of Natural Oils in Clay-Based Materials: Temperature and Humidity Effects. Langmuir. 2023, 39, 2079–2088. [Google Scholar] [CrossRef]
  27. Crank, J. The Mathematics of Diffusion, 2nd ed.; Oxford University Press, 1979. [Google Scholar]
  28. Wang, C.; Liu, F.; Zhang, Q. Diffusion and Desorption of Essential Oils in Porous Media: Cylindrical Geometry Considerations. Chem Eng Sci. 2020, 223, 115742. [Google Scholar] [CrossRef]
  29. Shao, L.; Cheng, X.; Zhang, R.; Yin, P.; Wang, Y. Evaporation and Diffusion Dynamics of Volatile Oils in Porous Clays: Application to Pest Control. J Agric Food Chem. 2022, 70, 9136–9144. [Google Scholar] [CrossRef]
  30. Abdellaoui, M.; Chahbouni, M.; El Bouardi, A.; Mokhlisse, A. A Study on the Rearing and Susceptibility of Sitophilus granarius Populations to Natural Insecticides. J Entomol Res. 2022, 46, 45–55. [Google Scholar] [CrossRef]
  31. Elhassaneen, Y.; El-Zahar, K.; Sabra, A.S. Evaluation of Citrus Essential Oils Against Sitophilus granarius: Impact on Mortality and Grain Quality. J Stored Prod Res. 2021, 92, 101778. [Google Scholar] [CrossRef]
  32. Abdoul-Latif, F.M.; Ainane, A.; Abdoul-Latif, T.M.; Ainane, T. Chemical study and evaluation of insecticidal properties of African Lippia citriodora essential oil. J Biopesticides. 2020, 13, 119–126. [Google Scholar] [CrossRef]
  33. Mohamed Abdoul-Latif F, El Montassir Z, Ainane A, et al. Use of Thymus plants as an ecological filler in urea-formaldehyde adhesives intended for bonding plywood. Processes. 2022, 10, 2209. [Google Scholar] [CrossRef]
  34. Shaaban, M.; Bayoumi, A.E.; Hafez, R.A. Lethal Dose and Toxicity Assessments of Essential Oils on Sitophilus granarius: Application in Integrated Pest Management. J Pest Sci. 2022, 95, 723–734. [Google Scholar] [CrossRef]
  35. Khadri, F.A.; Ali, W.; Khan, M.A. Factorial Design and Statistical Optimization of Essential Oil Transfer in Porous Clays: A Case Study. J Chem Technol Biotechnol. 2021, 96, 1556–1563. [Google Scholar] [CrossRef]
  36. Singh, R.S.; Gupta, P.; Das, S.K. Modeling and Optimization of Evaporation Rates in Essential Oil-Clay Systems Using Factorial Design. Ind Eng Chem Res. 2020, 59, 20998–21006. [Google Scholar] [CrossRef]
  37. Ribes S, Saidi N, Luyckx M, Ferrari B, Lachtar S, Saadaoui H, et al. Application of Principal Component Analysis (PCA) to Assess the Efficiency of Natural Insecticides in Different Porous Media. Chemom Intell Lab Syst. 2022, 226, 104555. [Google Scholar] [CrossRef]
  38. Farhat, R.; Hassan, S.; Chenini, S.; Zouari, N. Statistical Approaches in Evaluating the Insecticidal Activity of Essential Oils in Porous Supports Using PCA. Environ Sci Pollut Res Int. 2023, 30, 40761–40773. [Google Scholar] [CrossRef]
  39. Ghanmi, K.; Boussadia, C.; Saidi, K.; Ghemati, S. Computational Study on the Interaction of Limonene with Key Proteins 2NXX and 1C3Z: Implications for Biological Mechanisms and Docking Applications. J Mol Graph Model. 2023, 120, 108022. [Google Scholar] [CrossRef]
  40. Ghanty, T.K.; Chakraborty, M. Molecular Docking and Force Field Studies on Limonene for Protein Binding: A Theoretical Approach. Theor Chem Acc. 2022, 141, 35. [Google Scholar] [CrossRef]
  41. Mohamed Abdoul-Latif, F.; Ainane, A.; Achenani, L.; et al. Production of fucoxanthin from microalgae Isochrysis galbana of Djibouti: optimization, correlation with antioxidant potential, and bioinformatics approaches. Mar Drugs. 2024, 22, 358. [Google Scholar] [CrossRef]
  42. Kessler, H.; Pfaffenberger, M.; Retey, J. Protein Docking and Interaction Analysis of Ultraspiracle Receptor from Tribolium castaneum and Limonene: A Structural Perspective. Biochim Biophys Acta Proteins Proteom. 2023, 1871, 140929. [Google Scholar] [CrossRef]
  43. Tang, Y.; Sun, J.; Yan, H.; Wang, R. Solvent Accessibility and Functional Analysis of the Ultraspiracle Receptor in Insects: Insights from Molecular Dynamics Simulations. J Mol Graph Model. 2021, 108, 107979. [Google Scholar] [CrossRef]
  44. Sheng, M.; Li, Y.; Liu, X.; Wang, L. Structural and Functional Analysis of Antifreeze Protein THP12 from Tenebrio molitor: Implications for Cold Adaptation. J Mol Struct. 2022, 1261, 132835. [Google Scholar] [CrossRef]
  45. García-Martín, D.; Martí-Renom, M.A.; Sanz, F.; Pastor, M. Exploring Protein-Ligand Interactions with Docking Simulations: The Case of Limonene and Insect Receptor Proteins. J Chem Inf Model. 2021, 61, 2085–2097. [Google Scholar] [CrossRef]
  46. Duan, C.; Li, L.; Liu, Z.; Zhang, Y. Molecular Docking Studies and DFT-Based Optimization of Limonene Interaction with Insect Receptors. J Mol Model. 2023, 29, 45. [Google Scholar] [CrossRef]
  47. Zhang, Y.; Sun, H.; Zhang, Q.; Li, D. Application of MATLAB for Numerical Modeling and Statistical Analysis in Material Science. Comput Mater Sci. 2020, 171, 109194. [Google Scholar] [CrossRef]
  48. Wang, L.; Huang, Y.; Zhu, M.; Li, Q. Integration of Computational Tools for Data Analysis and Visualization in Environmental Studies. J Environ Sci. 2022, 116, 42–54. [Google Scholar] [CrossRef]
  49. Zhou, Y.; Ren, H.; Wang, H.; Li, X. XRD Characterization of Clay Minerals in Sedimentary Environments: Implications for Diagenetic Processes. Sediment Geol. 2021, 419, 106054. [Google Scholar] [CrossRef]
  50. Wang Y, Liu W, Zhou J, et al. Mineralogical Analysis of Clay Samples Using XRD: Insights into Formation and Diagenetic Conditions. Appl Clay Sci. 2023, 234, 106876. [Google Scholar] [CrossRef]
  51. Yao, L.; Chen, L.; Gao, P.; Zhang, Z. XRF Analysis of Mineral Compositions in Clay Samples: Differentiation Between Environmental and Geological Processes. Geochem J. 2022, 56, 263–276. [Google Scholar] [CrossRef]
  52. Qiu, Y.; Jin, Z.; Zhao, H.; Yang, C. XRF and XRD Combined Analysis for Mineralogical Studies of Clay Samples from Sedimentary Basins. Miner Mag. 2020, 84, 505–520. [Google Scholar] [CrossRef]
  53. Liu Q, Zhang X, Shen J, et al. FTIR Spectroscopic Analysis of Clay Minerals: Correlation Between Absorption Bands and Mineral Composition. Vib Spectrosc. 2021, 114, 103233. [Google Scholar] [CrossRef]
  54. Mahjoub R, Seddek A, Ahmed S, et al. FTIR and SEM-EDX Characterization of Natural Clays: Implications for the Study of Clay Mineralogy. Spectrochim Acta A Mol Biomol Spectrosc. 2023, 295, 122604. [Google Scholar] [CrossRef]
  55. Zhang M, Wang C, Sun Z, et al. SEM and EDX Analysis of Weathered Clays: Understanding Morphological and Compositional Changes. Appl Geochem. 2021, 128, 104922. [Google Scholar] [CrossRef]
  56. Zhang, Y.; Li, J.; Zhan, Y.; Zhao, X. Morphological and Elemental Analysis of Sedimentary Clays: A Comparative Study Using SEM-EDX. Micron. 2023, 167, 103389. [Google Scholar] [CrossRef]
  57. Alkhazraji A, Alhashemi B, Youssef K, et al. Thermogravimetric Analysis of Clay Minerals: Identifying Phases of Mass Loss and Thermal Decomposition. J Therm Anal Calorim. 2022, 147, 781–793. [Google Scholar] [CrossRef]
  58. Zhao R, Yang L, Wang P, et al. Thermogravimetric and Differential Thermal Analysis of Kaolinite and Illite in Clay Samples.
  59. Nasser AL-Jabri, N.; Aflah, N. Comparative Analysis of Lemon Essential Oils from Turkey and India: Chemical Composition and Biological Activities. J Essent Oil Res. 2022, 34, 137–145. [Google Scholar] [CrossRef]
  60. Njoku, V.N.; Evbuomwan, B.O.; Ude, N.O. Extraction and Characterization of Lemon (Citrus limon) Essential Oil. Int J Multidisc Res Dev. 2016, 3, 130–136. [Google Scholar]
  61. Bertuzzi, T.; Rastelli, S.; Mulazzi, A.; Pietri, A. Chemical composition and variability of lemon essential oil (Citrus limon) in relation to cultivation areas and fruit maturation stages. J Agric Food Chem. 2013, 61, 1387–1394. [Google Scholar] [CrossRef]
  62. Himd, A.; Ouali, K.; Boudoukha, C.; Gherraf, N.; Zeghichi, M.T. Characterization and Antibacterial Activity of Lemon (Citrus limon) Peel Essential Oil. Int J Pharm Pharm Sci. 2016, 8, 297–301. [Google Scholar]
  63. Moufida, S.; Marzouk, B. Biochemical Characterization of Blood Orange, Sweet Orange, Lemon, Bergamot and Bitter Orange. J Agric Food Chem. 2003, 51, 6398–6404. [Google Scholar] [CrossRef] [PubMed]
  64. Haddad M, Bennis N, Belmekki A, et al. Enhanced Insecticidal Efficacy of Essential Oils Supported on Porous Clays Against Sitophilus granarius. J Pest Sci. 2021, 94, 539–550. [Google Scholar] [CrossRef]
  65. Taoufik M, Cherifi Z, Bouharroud R, et al. Influence of Temperature and Application Methods on the Insecticidal Activity of Citrus Essential Oils in Stored Grain Pest Control. J Stored Prod Res. 2022, 96, 101992. [Google Scholar] [CrossRef]
  66. Boukaew S, Prasertsan P, Siriamornpun S, et al. Role of Porous Supports in Enhancing the Stability and Insecticidal Activity of Essential Oils Against Stored Product Pests. Ind Crops Prod. 2023, 192, 115888. [Google Scholar] [CrossRef]
  67. Zhou, Z.; Li, M.; Sun, Y.; Wang, J. Temperature-Dependent Desorption Kinetics of Essential Oils from Clay Matrices: Implications for Insecticidal Activity. J Agric Food Chem. 2021, 69, 3883–3892. [Google Scholar] [CrossRef]
  68. Kim, Y.J.; Ahn, Y.J. Enhanced Release and Insecticidal Efficacy of Essential Oils from Clay Supports Under Elevated Temperatures. J Pest Sci. 2020, 93, 1331–1343. [Google Scholar] [CrossRef]
  69. Mahdavi V, Tohidi M, Abedi G, et al. The Impact of Concentration on the Desorption Behavior of Essential Oils in Porous Supports: Applications for Prolonged Insecticidal Action. Chem Eng Sci. 2023, 260, 117987. [Google Scholar] [CrossRef]
  70. Mahrishi, S.; Bhattacharyya, R.; Kumar, D.; Singh, B. Influence of Particle Size on the Desorption Dynamics of Essential Oils in Clay Media. Ind Eng Chem Res. 2022, 61, 1245–1257. [Google Scholar] [CrossRef]
  71. Yang, L.; Ren, H.; Zhang, L.; Hu, Z. Effect of Clay Mass on the Desorption Rate and Insecticidal Performance of Essential Oils in Controlled Release Systems. J Chem Technol Biotechnol. 2021, 96, 2450–2458. [Google Scholar] [CrossRef]
  72. Huang, Q.; Zhou, J.; Gao, Y.; Sun, W. Optimizing the Release Kinetics of Essential Oils in Clay-Based Insecticidal Formulations: A Comprehensive Analysis. J Clean Prod. 2023, 364, 132761. [Google Scholar] [CrossRef]
  73. Ortega, L.F.; Franco, L.A.; Palacios, S.M. Modeling LD50 Values in Insecticidal Assays: A Case Study with Lemon Essential Oil and Stored Product Pests. J Econ Entomol. 2022, 115, 78–86. [Google Scholar] [CrossRef]
  74. Li, X.; Zhang, W.; Peng, Q.; Wang, H. Optimization of Essential Oil Evaporation Flux in Porous Media Using Factorial Design: A Case Study in Clays. Chem Eng Res Des. 2021, 167, 215–226. [Google Scholar] [CrossRef]
  75. Chen Z, Zhou Y, Liu W, et al. Modeling Essential Oil Flux in Porous Clays: Effects of Concentration, Temperature, and Particle Size. J Chem Technol Biotechnol. 2022, 97, 1480–1492. [Google Scholar] [CrossRef]
  76. Yang, Y.; Luo, S.; Zhang, X. Analyzing Interactions Between Factors Affecting Essential Oil Flux in Porous Supports Using 3D Modeling. Ind Crops Prod. 2023, 189, 115936. [Google Scholar] [CrossRef]
  77. Pérez-Martínez, E.; Gómez-Nieto, M.A.; Hernández-Córdoba, M. Evaluating Predictive Models for Essential Oil Flux in Porous Media: Insights from Statistical Analysis. J Chem Inf Model. 2021, 61, 4390–4402. [Google Scholar] [CrossRef]
  78. Zhao H, Sun X, Chen S, et al. Statistical Optimization and Model Validation for Essential Oil Transfer in Porous Media. Chem Eng Sci. 2023, 259, 117941. [Google Scholar] [CrossRef]
  79. Abatzoglou, N.; Barmpalexis, P.; Kostoglou, M. Modeling of Essential Oil Desorption from Porous Materials: Application to Fickian Diffusion Processes. Ind Eng Chem Res. 2020, 59, 19891–19901. [Google Scholar] [CrossRef]
  80. Tiwari, G.; Tiwari, R.; Srivastava, B.; Rai, A.K. Mathematical Modeling and Simulation of Essential Oil Diffusion in Porous Clays: A Fickian Approach. J Chem Phys. 2022, 156, 104503. [Google Scholar] [CrossRef]
  81. Lu Z, Zhang J, Hu W, et al. Comparative Analysis of Diffusion Models for Essential Oils in Different Porous Media: The Case of Red and Green Clays. J Porous Mater. 2021, 28, 1405–1416. [Google Scholar] [CrossRef]
  82. Yousefi, A.; Shahsavar, M.; Hosseini, M.J. Simulation and Optimization of Essential Oil Diffusion in Clay-Based Media for Agricultural Applications. Chem Eng Sci. 2023, 252, 117501. [Google Scholar] [CrossRef]
  83. Pohl, A.; Karczewska, B.; Żmuda, E. Study on Adsorption and Desorption Processes of Essential Oils in Porous Clays: Implications for Activation Energy and Retention Capacity. Appl Clay Sci. 2021, 205, 106036. [Google Scholar] [CrossRef]
  84. Ng KH, Lam SS, Mun SP, et al. Multivariate Analysis of Essential Oil Performance in Porous Supports: Application of Principal Component Analysis in Kinetic Studies. J Mol Liq. 2023, 380, 121694. [Google Scholar] [CrossRef]
  85. Ribeiro A, Marques P, Ferreira R, et al. Correlating Evaporation Flux and Insecticidal Activity of Essential Oils in Clays Using Principal Component Analysis. J Chem Inf Model. 2022, 62, 1123–1131. [Google Scholar] [CrossRef]
  86. Zhang Y, Chen L, Zhu J, et al. Understanding the Diffusion-Evaporation Relationship of Essential Oils in Clay Media: Insights from Principal Component Analysis. Ind Crops Prod. 2022, 185, 115096. [Google Scholar] [CrossRef]
  87. Mahdavi V, Bahrami A, Hashemzadeh R, et al. Multivariate Analysis of Desorption and Insecticidal Efficacy of Essential Oils in Porous Media: A Case Study with Red and Green Clays. Chemosphere. 2022, 308, 136354. [Google Scholar] [CrossRef]
  88. Pehlivanoglu, M.; Saglam, O. Evaluating the Diffusion and Evaporation Dynamics of Essential Oils in Clay Media Using Principal Component Analysis. Appl Clay Sci. 2021, 204, 106008. [Google Scholar] [CrossRef]
  89. Zhang X, Tian Z, Liu L, et al. Principal Component Analysis in Insecticidal Studies: Understanding the Interactions Between Diffusion and Evaporation in Essential Oil-Impregnated Clays. J Hazard Mater. 2023, 442, 130012. [Google Scholar] [CrossRef]
  90. Zouari N, Ayari H, Smaoui S, et al. Analyzing the Effects of Diffusion and Evaporation on Insecticidal Activity: A PCA Approach with Citrus Essential Oils. Environ Sci Pollut Res Int. 2022, 29, 22451–22461. [Google Scholar] [CrossRef]
  91. Yildirim N, Sen B, Gul S, et al. Energetic Analysis and Molecular Orbital Studies of Limonene: Insights into Chemical Stability and Reactivity. Comput Theor Chem. 2021, 1203, 113355. [Google Scholar] [CrossRef]
  92. Liu Y, Liu X, Li J, et al. Molecular Docking Studies and Energetic Characterization of Limonene with Insect Receptor Proteins. J Mol Graph Model. 2023, 122, 108080. [Google Scholar] [CrossRef]
  93. Dong Z, Li Y, Zhang L, et al. Comparative Binding Affinity and Stability of Limonene with Insect Receptors: Insights from Molecular Docking Studies. J Mol Graph Model. 2021, 104, 107826. [Google Scholar] [CrossRef]
  94. Verma A, Joshi K, Jain SK, et al. Docking and Binding Affinity Analysis of Limonene with Insect and Human Receptor Proteins. J Biomol Struct Dyn. 2022, 40, 2904–2914. [Google Scholar] [CrossRef]
  95. Singh RK, Chikara C, Sharma A, et al. Evaluation of Ligand Efficiency and Binding Affinity in Protein-Limonene Interactions: A Computational Study. J Comput Aided Mol Des. 2020, 34, 1215–1230. [Google Scholar] [CrossRef]
  96. Kumar, P.; Rawat, S.; Misra, G. Conformational Analysis and Binding Affinity of Limonene with Target Proteins: A Molecular Docking Perspective. J Mol Struct. 2023, 1282, 135083. [Google Scholar] [CrossRef]
  97. Joo K, Kim D, Kim J, et al. Principal Component Analysis of Docking Parameters: Understanding Ligand Efficiency and Binding Affinity in Protein-Ligand Complexes. J Chem Inf Model. 2022, 62, 749–759. [Google Scholar] [CrossRef]
  98. Banerjee, S.; Poddar, A.; Chakraborty, S. Multivariate Analysis of Docking Results: A Focus on Ligand Efficiency and Binding Free Energy Correlations. Comput Biol Chem. 2023, 99, 107799. [Google Scholar] [CrossRef]
  99. Wang X, Zhang Z, Li P, et al. Protein-Ligand Interaction Analysis Using Principal Component Analysis: Implications for Drug Design. ACS Omega. 2021, 6, 19826–19836. [Google Scholar] [CrossRef]
  100. Maiti, R.; Reddy, P.; Sobhia, M.E. Correlation Analysis in Protein-Ligand Docking: A Comprehensive Study Using Principal Component Analysis. J Mol Model. 2023, 29, 77. [Google Scholar] [CrossRef]
Preprints 137449 g001
Figure 2. X-Ray diffraction diagram of RC and GC clay. (K = Kaolinite, I = Illite, Q = Quartz, C = Calcite, D = Dolomite, H = Hematite).
Figure 2. X-Ray diffraction diagram of RC and GC clay. (K = Kaolinite, I = Illite, Q = Quartz, C = Calcite, D = Dolomite, H = Hematite).
Preprints 137449 g001
Preprints 137449 g002
Figure 3. IRTF spectra of two RC and GC clay samples.
Figure 3. IRTF spectra of two RC and GC clay samples.
Preprints 137449 g002
Preprints 137449 g003
Figure 4. RC SEM (resolution 50 μm and 100 μm).
Figure 4. RC SEM (resolution 50 μm and 100 μm).
Preprints 137449 g003
Preprints 137449 g004
Figure 5. GC SEM (resolution 50 μm and 100 μm).
Figure 5. GC SEM (resolution 50 μm and 100 μm).
Preprints 137449 g004
Preprints 137449 g005
Figure 6. Thermogravimetric analysis of two RC and GC samples.
Figure 6. Thermogravimetric analysis of two RC and GC samples.
Preprints 137449 g005
Preprints 137449 g006
Figure 7. Desorption Rate and Residual Mass Analysis of Lemon Essential Oil in Red and Green Clays under Varying Condition.
Figure 7. Desorption Rate and Residual Mass Analysis of Lemon Essential Oil in Red and Green Clays under Varying Condition.
Preprints 137449 g006
Preprints 137449 g007
Figure 8. 3D representation of the flow as a function of each two factors in the RC clay media.
Figure 8. 3D representation of the flow as a function of each two factors in the RC clay media.
Preprints 137449 g007
Preprints 137449 g008
Figure 9. 3D representation of the flow as a function of each two factors in the GC clay media.
Figure 9. 3D representation of the flow as a function of each two factors in the GC clay media.
Preprints 137449 g008
Preprints 137449 g009
Figure 10. Analytical simulation of the desorption kinetics of EOCl in the RC and GC medium. (•): Experimental results; (…): Simulation, C = 0.01 m L/cm3; T = 30 ° C; M = 0.10 g; D = 2 cm.
Figure 10. Analytical simulation of the desorption kinetics of EOCl in the RC and GC medium. (•): Experimental results; (…): Simulation, C = 0.01 m L/cm3; T = 30 ° C; M = 0.10 g; D = 2 cm.
Preprints 137449 g009
Preprints 137449 g010
Figure 11. Comparative simulation of EOCL desorption concentration profiles in RC and GC media.
Figure 11. Comparative simulation of EOCL desorption concentration profiles in RC and GC media.
Preprints 137449 g010
Preprints 137449 g011
Figure 12. Mapping of Correlations Between Parameters of EOCL Desorption and Insecticidal Activity in RC and GC Medium.
Figure 12. Mapping of Correlations Between Parameters of EOCL Desorption and Insecticidal Activity in RC and GC Medium.
Preprints 137449 g011
Preprints 137449 g012
Figure 13. Energy level diagram of limonene molecular orbitals.
Figure 13. Energy level diagram of limonene molecular orbitals.
Preprints 137449 g012
Preprints 137449 g013aPreprints 137449 g013b
Figure 14. Three-Dimensional and 2D Docked Views of Limonene with 2NXX and 1C3Z Proteins, respectively.
Figure 14. Three-Dimensional and 2D Docked Views of Limonene with 2NXX and 1C3Z Proteins, respectively.
Preprints 137449 g013aPreprints 137449 g013b
Preprints 137449 g014
Figure 15. PCA Biplot with Proteins and Docking Parameters.
Figure 15. PCA Biplot with Proteins and Docking Parameters.
Preprints 137449 g014
Table 1. Matrix of Experiments and Evaporation Rates (Flux) F Obtained for Each Test.
Table 1. Matrix of Experiments and Evaporation Rates (Flux) F Obtained for Each Test.
Test Factor 1 (C) Factor 2 (T) Factor 3 (D) Factor 4 (M) F (×10⁴) (g·h⁻¹·cm⁻¹) for RC F (×10⁴) (g·h⁻¹·cm⁻¹) for GC
1 -1 -1 -1 -1 1.2732 1.8687
2 1 -1 -1 -1 3.8197 2.886
3 -1 1 -1 -1 2.5464 2.1248
4 1 1 -1 -1 11.4591 6.2191
5 -1 -1 1 -1 1.2732 1.6471
6 1 -1 1 -1 1.2732 1.4987
7 -1 1 1 -1 2.5464 2.5367
8 1 1 1 -1 2.5464 3.5786
9 -1 -1 -1 1 0.6366 1.3582
10 1 -1 -1 1 0.6366 1.2191
11 -1 1 -1 1 2.2281 2.7483
12 1 1 -1 1 2.5464 3.1296
13 -1 -1 1 1 1.2732 1.9371
14 1 -1 1 1 1.2732 2.1582
15 -1 1 1 1 2.8647 3.7489
16 1 1 1 1 3.5014 4.5672
Table 3. Equations to calculate the numerical energetic descriptors.
Table 3. Equations to calculate the numerical energetic descriptors.
Property Symbol Equation Description
Energy Gap EGAP ELUMO - EHOMO Difference in energy between LUMO and HOMO orbitals.
Chemical Hardness η E L U M O − E H O M O 2 Measure of a molecule resistance to change in electron distribution.
Electronegativity χ − E L U M O + E H O M O 2 Ability of an atom or molecule to attract electrons.
Electrophilicity Index ω χ 2 2 η Quantifies the electrophilic nature of a molecule.
Molecular Flexibility S 1 2 η Indicates the flexibility of a molecule based on its chemical hardness.
Table 4. The active site coordinates of studied proteins.
Table 4. The active site coordinates of studied proteins.
Protein 2NXX 1C3Z
Size (Å) x = 146; y = 140; z = 188 x = 73; y = 59; z = 49
Center (Å) x = 146
y = 140
z = 188		 x = 1.217
y = 1.626
z = -8.302
Table 5. Chemical composition of RC and GC.
Table 5. Chemical composition of RC and GC.
Samples Al2O3 CaO Fe2O3 K2O MgO Na2O P2O5 SO3 SiO2 TiO2 LOI
RC 17,03 0,68 32,47 3,27 1,24 3,60 0,10 0,01 33,40 0,69 8,35
GC 29,06 0,58 1,61 4,06 1,27 3,85 0,11 0,01 52,17 1,00 6.74
Table 6. EDX Microanalysis of RC and GC Surfaces.
Table 6. EDX Microanalysis of RC and GC Surfaces.
Élément GC (%) RC (%)
Si 27.26 15.49
Al 18.7 10.74
O 48.76 36.34
Fe - 23.27
Mg - 2.87
K 5.28 2.87
C - 8.41
Total 100 100
Table 7. Constituents identified in the oils of Citrus limonum.
Table 7. Constituents identified in the oils of Citrus limonum.
Pic Compounds RT RI Area %
1 α-Thujene 5.504 732.41 0.53
2 α-Pinene 5.733 750.68 4.84
3 Camphene 6.199 785.70 0.69
4 β-Pinene 7.247 855.70 6.46
5 Myrcene 7.726 884.38 2.67
6 α-Phellandrene 8.21 911.61 0.47
7 Limonene 9.463 975.26 51.41
8 trans-β-Ocimene 9.652 984.13 0.30
9 cis-β-Ocimene 10.054 1002.41 0.84
10 γ-Terpinene 10.582 1025.35 4.60
11 Linalool 12.296 1092.63 3.16
12 Fenchyl alcohol 12.832 1111.75 0.49
13 cis-Sabinene hydrate 13.192 1124.15 0.09
14 β-Ocimene 13.765 1143.20 0.10
15 β-Terpineol 14.231 1158.12 0.20
16 α-Phellandren-8-ol 15.276 1189.88 0.49
17 α-Terpineol 16.468 1223.55 4.13
18 Neral 18.07 1265.15 3.71
19 Nerol 18.07 1265.15 0.36
20 cis-Citral 18.619 1278.56 1.27
21 Geraniol 19.284 1294.29 0.41
22 Geranial 19.979 1310.16 2.11
23 δ-Elemene 22.755 1368.46 0.21
24 Neryl acetate 24.038 1393.04 2.53
25 β-Elemene 25.095 1412.33 0.20
26 trans-Caryophyllene 26.193 1431.52 1.06
27 cis-α-Bergamotene 26.932 1443.99 1.97
28 Selina-4,11-diene 28.518 1469.63 0.10
29 Germacrene D 28.723 1472.84 0.20
30 β-Caryophyllene 28.911 1475.76 0.14
31 β-Cadinene 29.165 1479.68 0.11
32 α-Selinene 29.295 1481.68 0.10
33 β-Bisabolene 29.599 1486.17 4.03
Total 99.98
Table 8. Parameters of the Polynomial Flux Model for Desorption (RC and GC).
Table 8. Parameters of the Polynomial Flux Model for Desorption (RC and GC).
Coefficient Value (*) for RC ptest for RC Value (*) for GC ptest for GC
b0 2.6061 0.46 2.6115 0.45
bC 0.7758 0.31 0.7792 0.32
bT 1.1737 0.21 1.1824 0.22
bD -0.5371 0.42 -0.5338 0.41
bM -0.7360 0.32 -0.7317 0.33
bCT 0.4575 0.47 0.4602 0.46
bCD -0.6963 0.34 -0.6937 0.35
bCM -0.6565 0.36 -0.6528 0.37
bTD -0.3779 0.53 -0.3756 0.52
bTM -0.2586 0.64 -0.2553 0.63
bDM 0.8952 0.27 0.8984 0.26
bCTD -0.3779 0.53 -0.3754 0.54
bCTM -0.3382 0.56 -0.3357 0.55
bCDM 0.7360 0.32 0.7385 0.33
bTDM 0.4177 0.50 0.4203 0.51
bCTDM 0.2745 - 0.2783 -
*The values presented are of the order of ×104 (g.h-1.cm-1).
Table 9. Desorption Results for Red and Green Clays.
Table 9. Desorption Results for Red and Green Clays.
Clay Type Temperature (°C) Concentration (µl/cm³) Particle Diameter (mm) Clay Mass (g) Desorption Rate (%) Residual Mass (%)
Red Clay (RC) 25 0.01 0.5 5 25 75
Red Clay (RC) 30 0.02 1.0 10 40 60
Red Clay (RC) 35 0.03 1.5 15 55 45
Green Clay (GC) 25 0.01 0.5 5 30 70
Green Clay (GC) 30 0.02 1.0 10 45 55
Green Clay (GC) 35 0.03 1.5 15 60 40
Table 10. Parameters of Diffusion and Insecticidal Activity of EOCL Desorption in the RC Medium.
Table 10. Parameters of Diffusion and Insecticidal Activity of EOCL Desorption in the RC Medium.
Test F (×10⁴) (g.h⁻¹.cm⁻¹) Dz (×10⁶) (cm².h⁻¹) K (×10³) Ea (J.mol⁻¹) M%
1 1.2732 7.35 1.5 -8.176 16.66
2 3.8197 1.47 0.3 20
3 2.5464 3.43 0.7 20
4 11.4591 2.45 0.5 30
5 1.2732 9.31 1.9 20
6 1.2732 2.94 0.6 30
7 2.5464 3.92 0.8 23.33
8 2.5464 1.96 0.4 33.33
9 0.6366 15.68 0.8 20
10 0.6366 5.88 0.3 16.66
11 2.2281 43.12 2.2 23.33
12 2.5464 17.64 0.9 23.33
13 1.2732 25.48 1.3 20
14 1.2732 9.8 0.5 30
15 2.8647 39.2 2.0 30
16 3.5014 21.56 1.1 30
Table 11. Parameters of Diffusion and Insecticidal Activity of EOCL Desorption in the GC Medium.
Table 11. Parameters of Diffusion and Insecticidal Activity of EOCL Desorption in the GC Medium.
Test F (×10⁴) (g.h⁻¹.cm⁻¹) Dz (×10⁶) (cm².h⁻¹) K (×10³) Ea (J.mol⁻¹) M%
1 1.2732 8.12 1.6 -7.854 17.50
2 3.8197 1.54 0.4 22
3 2.5464 3.76 0.8 22
4 11.4591 2.62 0.6 32
5 1.2732 9.92 2.0 22
6 1.2732 3.12 0.7 32
7 2.5464 4.13 0.9 24.67
8 2.5464 2.12 0.5 34.50
9 0.6366 16.90 0.9 22
10 0.6366 6.25 0.4 18.34
11 2.2281 45.75 2.3 24.67
12 2.5464 18.75 1.0 24.67
13 1.2732 27.11 1.4 22
14 1.2732 10.50 0.6 32
15 2.8647 42.50 2.1 32
16 3.5014 23.33 1.2 32
Table 12. Eigenvalues and Variability of the Principal Components in PCA for EOCL Desorption in RC and GC Media.
Table 12. Eigenvalues and Variability of the Principal Components in PCA for EOCL Desorption in RC and GC Media.
F1 F2 F3 F4
Eigenvalue 1.89 1.38 0.53 0.18
Variability (%) 47.31 34.70 13.45 4.53
Cumulative (%) 47. 31 82.00 95.46 100
Table 13. Calculated energy values of limonene.
Table 13. Calculated energy values of limonene.
Property Value (eV)
Energy Gap (EGAP) 14.62
Chemical Hardness (η) 7.31
Electronegativity (χ) 6.651
Electrophilicity Index (ω) 3.026
Molecular Flexibility (S) 0.068
Table 14. Energetic and qualitative characteristics of molecular docking of cathinone with the proteins studied.
Table 14. Energetic and qualitative characteristics of molecular docking of cathinone with the proteins studied.
Proteins Binding
Affinity (kcal/mol)		 pKi Ligand
Efficiency (kcal/mol)		 Ligand–
Protein
Interactions		 Number of Conventional
Hydrogen Bonds		 Number of Aromatic Bonds Number of Van Der Waals Bonds
2NXX -6.5 4.77 0.65 Hydrophobic, Aromatic 0 2 3
1C3Z -6.3 4.62 0.63 Hydrophobic, Aromatic 0 2 3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Alerts
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2025 MDPI (Basel, Switzerland) unless otherwise stated