Laurus nobilis L. Essential Oil-Loaded PLGA as a Nanoformulation Candidate for Cancer Treatment

The aim of this study was to obtain essential oil (LNEO) from the Laurus nobilis L. plant, and to prepare LNEO-loaded poly lactic-co-glycolic acid (PLGA) nanoparticles (NPs) as an approach in cancer treatment. The components of the obtained LNEO were analyzed using GC-MS. The LNEO-NPs were synthesized by the single-emulsion method. The LNEO-NPs were characterized using UV-Vis spectrometry, Dynamic Light Scattering (DLS), Scanning Electron Microscopy (SEM), and a DNA binding assay, which was performed via the UV-Vis titration method. According to the results, the LNEO-NPs had a 211.4 ± 4.031 nm average particle size, 0.068 ± 0.016 PdI, and −7.87 ± 1.15 mV zeta potential. The encapsulation efficiency and loading capacity were calculated as 59.25% and 25.65%, respectively, and the in vitro drug release study showed an LNEO release of 93.97 ± 3.78% over the 72 h period. Moreover, the LNEO was intercalatively bound to CT-DNA. In addition, the mechanism of action of LNEO on a dual PI3K/mTOR inhibitor was predicted, and its antiproliferative activity and mechanism were determined using molecular docking analysis. It was concluded that LNEO-loaded PLGA NPs may be used for cancer treatment as a novel phytotherapeutic agent-based controlled-release system.


Introduction
According to the World Health Organization (WHO), cancer is the second leading cause of death worldwide, accounting for around 10 million deaths annually. Globally, one in six deaths are associated with cancer. It is estimated that cancer-related deaths will increase by 70% in 2027 [1]. Cancer-causing infections, such as hepatitis and human papilloma virus (HPV), are responsible for approximately 30% of cancer cases in low-and middle-income countries [2]. More than 90% of high-income countries have comprehensive treatment, while less than 15% of low-income countries have this opportunity [3].
Traditional cancer treatment includes different interventions such as chemotherapy, surgery and radiation therapy; however, these treatment modalities also damage the adjacent healthy cells [4]. New anti-proliferative drugs are used for cancer treatment, but they come with severe adverse effects, toxicity and have high costs. Therefore, costeffective and biocompatible therapeutic applications are needed [5]. New natural products will play an important role in this period. Many of the antitumor agents currently used in clinical practice are of natural origin. Among the great diversity of phytochemicals,

LNEO Composition
Dry Laurus nobilis leaf yielded 1.98 ± 0.16 mL of essential oil/100 g on average after three rounds of hydrodistillation. A total of 11 compounds in laurel leaf essential oils were detected and quantified. Table 1 presents the constituents in terms of the retention indices' order. The components of Laurus nobilis leaf essential oil were found to be 64.556% 1.8 cineole as the major compound, 10.099% sabinene, and 6.180% α-pinene, respectively. Oxygenated Monoterpenes (7,(9)(10)(11)(12)

Average Particle Size, PdI and Zeta Potential Analysis Results
The average particle size, PdI, and zeta potential of nanoparticles play an important role in the interaction of bioactive substances with the target tissue in drug delivery systems to produce a therapeutic effect [17][18][19]. DLS is one of the widely used methods to determine the average particle size, PdI, and zeta potential of synthesized nanoparticles [19][20][21]. In this study, the average particle size, PdI and zeta potential of the nanoparticles were determined using DLS. The DLS data showed that the mean particle size of the blank NPs was 182.4 ± 0.499 nm, and the PdI value was 0.071 ± 0.006 (Figure 1a). The zeta potential of the empty NPs was −10.5 ± 0.4 mV (Figure 1b). The mean particle size of the LNEO-NPs was 211.4 ± 4.031 nm, and their PdI was 0.068 ± 0.016 (Figure 2a). The zeta potential of the LNEO-NPs was −7.87 ± 1.15 mV (Figure 2b).  Furthermore, the band edge positions of these MoTe 2 /PtS 2 vdWHs are obtained b HSE06 calculations in Figure 3, and compared with the HER and OER potential energ for the water-splitting at pH 0. One can see that the P-4 and P-6 MoTe 2 /PtS 2 vdWH possess a decent band alignment, which can promote the HER and OER for water-splitting at pH 0. Additionally, together with a type-II band structure, the P-6 MoTe 2 /PtS 2 vdWH can be considered as an advantageous photocatalyst to decompose water. Although th P-4 MoTe 2 /PtS 2 vdWH has a type-II band alignment to separate the photogenerated electrons and holes, the suitable band edge positions still have the ability to induce redo for water-splitting. Furthermore, the type-I heterostructure is also reported to show a nove photocatalytic performance [41,42].

SEM Micrograph of the LNEO-NPs
The morphology of the LNEO-NPs was determined using SEM. The SEM image ( Figure 3) confirmed the synthesis of spherical LNEO-NPs [22], and showed that the LNEO-NPs were in a homogeneous size distribution.

Determination of the Encapsulation Efficiency and Loading Capacity
The standard curve of the LNEO was prepared in order to determine the encapsulation efficiency and loading capacity of the LNEO-NPs (Figure 4). The encapsulation efficiency and loading capacity of the LNEO-loaded PLGA nanoparticles were calculated as 59.25% and 25.65%, respectively.

In Vitro Release Profile of the LNEO
The release profile of the LNEO-NP was evaluated as a function of time. Figure 5 shows the 72 h release results of the LNEO. The results showed that 47.76% ± 2.47% of the LNEO was released in the first six hours, and at the end of 72 h, 93.89% ± 4.93% of the LNEO had been released from the PLGA nanoparticles to the medium.

DNA Binding
UV-Vis absorption spectrophotometry is a useful technique to investigate the interaction between anticancer drug molecules and DNA [23][24][25]. Therefore, we investigated the potential of LNEO as an anticancer drug molecule using UV-Vis absorption spectrophotometry. In the study, CT-DNA was used in the interaction of the LNEO with DNA. The absorption spectra of the LNEO in the presence and absence of CT-DNA are given in Figure 6. LNEO was found to exhibit a 93.80% hypochromic effect at the 201 nm wavelength, and 6 nm of bathochromic (red) shift.

Molecular Docking and ADME Analysis Results
The inhibitory effects of the PI3K/mTOR target receptor by 1,8-cineole, α-terpinyl acetate, methyleugenol and sabinene molecules-which are the most dominant components in laurel essential oil (LNEO)-were analyzed, and their antiproliferative activity was evaluated by molecular docking analysis ( Table 2). Table 2 also shows the ligand-receptor interactions within the target's binding site, as well as the binding energy with each ligand. The docked poses of four major ingredients of LNEO are also shown in Figure 7    The active binding residues in the PI3K/mTOR target receptor (PDB 4FA6) were visualized for four major components in Figure 8 and Table 2.
The amino acid Val 882, located in the binding pocket of 4FA6, is crucial for the inhibition of PI3Kα-mTOR; as can be seen from the interactions, all four major compounds are bound to the protein from this binding pocket, even by two separate hydrogen bonds (2.30Å and 2.43 Å) of methyleugenol with Val 882 (Figure 9c), which is the most stable binding component in this pocket, with a docking score energy of −5.97 kcal/mol. In addition, these four major compounds also showed hydrophobic interactions with Met 953, Ile 831, Ile 879 and Ile 963 by the co-crystallized ligand [26,27] interacting with the target receptor. α-Terpinyl acetate also has one hydrogen bonding interaction with Asp 950 in the PI3K/mTOR target receptor ( Figure 9b).
The hydrogen bonding interactions of methyleugenol with Val 882 ( Figure 10a) and α-terpinyl acetate with Asp 950 in the PI3K/mTOR target receptor are illustrated in Figure 10. The electrostatic potentials of the binding pocket of the PI3K/mTOR target receptor and 1,8-cineole, α-terpinyl acetate, methyleugenol and sabinene molecules are also given in Figure 11, respectively.    Furthermore, the band edge positions of these MoTe 2 /PtS 2 vdWHs are obtained by HSE06 calculations in Figure 3, and compared with the HER and OER potential energy for the water-splitting at pH 0. One can see that the P-4 and P-6 MoTe 2 /PtS 2 vdWHs possess a decent band alignment, which can promote the HER and OER for water-splitting at pH 0. Additionally, together with a type-II band structure, the P-6 MoTe 2 /PtS 2 vdWH can be considered as an advantageous photocatalyst to decompose water. Although the P-4 MoTe 2 /PtS 2 vdWH has a type-II band alignment to separate the photogenerated electrons and holes, the suitable band edge positions still have the ability to induce redox for water-splitting. Furthermore, the type-I heterostructure is also reported to show a novel photocatalytic performance [41,42].  The interfacial characteristics of the MoTe 2 /PtS 2 vdWHs are investigated using the charge density difference and potential drop across the interface. Using Equation (2), the charge density difference in these MoTe 2 /PtS 2 vdWHs is calculated, showing that the PtS 2 layer obtained the electrons from the MoTe 2 layer ( Figure 4). Additionally, the charge transfer between the MoTe 2 and PtS 2 was investigated by Bader charge analysis [43], as shown in Table 1, which suggests a maximum charge transfer in P6-MoTe 2 /PtS 2 vdWH o about 0.047 electrons. Furthermore, the potential drop (shown by Figure 5) in the interface for MoTe 2 /PtS 2 vdWHs is addressed in Table 1. One can see that MoTe 2 /PtS 2 vdWHs possesses a pronounced potential drop across the interface, ranging from 4.41 eV to 4.67 eV which is larger than that in the WSSe/BSe vdW heterostructure [44]. Importantly, this Among the four major components, the 1,8-cineole molecule contains only negative electrostatic potential regions, which are shown in red; α-terpinyl acetate and methyleugenol both have negative electrostatic regions, shown in red, and positive potential regions shown in blue, while sabinene has no negative or positive regions (Table 3). For sabinene, the number of donor and acceptor hydrogen bonds was zero. As can be seen in Table 3, 1,8-cineole and sabinene, which have no rotatable single bonds, unfortunately failed to bind to the active binding site of the target protein with hydrogen bonds. The log P for octanol/water was given as 2.417, 2.919, 1.892 and 5.106 values for 1,8-cineole, α-terpinyl acetate, methyleugenol and sabinene, respectively, and was found to be in the range of 95% of drugs (−2.0/6.5). We also calculated the log BB for the brain/blood values of four major constituents, which were determined to be in the range from −3.0/to 1.2. All four compounds indicated very high Caco-2 and MDCK permeability properties ( Table 3). Each of these four compounds, which are in compliance with the Lipinski Rule of Five, have a high human oral absorption percentage, have low molecular weights, and are promising as drug candidates; therefore, the possess robust ADME properties, and laurel essential oil may also have such an activity.

Discussion
Traditional cancer treatment methods such as chemotherapy and radiotherapy have a wide range of side effects, such as cardiotoxicity, myelosuppression, hepatotoxicity and neurotoxicity [28][29][30]. The efforts of modern medicine are directed towards the development of drug delivery systems (DDSs) to increase the effectiveness of these treatments, and to minimize their systemic side effects. The milestone for DDSs is the use of Doxil ® , which was the first commercially available DDS. Doxorubicin (DOX) should be masked due to its toxic properties to healthy cells. In addition, the therapeutic effect of DOX can be enhanced by encapsulation. Doxil ® , a DOX-containing liposomal shell, minimizes side effects and increases the efficacy of treatment [31]. In addition to Doxil ® , Abraxane ® [32] and Onivyde ® [33] are commercially used anticancer drugs. Nanoparticles of various structures, especially polymeric ones, are the most promising platform for the creation of such DDSs [34,35]. PLGA polymer is one of the most widely used polymers in the development of nanoparticles as anticancer agents [36][37][38][39][40]. Laurel essential oil has been used as an antimicrobial, antioxidant [41,42] diuretic [43], and for the relief of rheumatic pains [44]. Additionally, the essential oil of bay leaves is widely used in the perfume and soap industries [45,46]. However, there are no studies on their anticancer activity in the literature. In this study, we developed PLGA nanoparticles loaded with laurel essential oil for cancer therapy.
The DLS analysis results showed that the mean droplet size of the developed nanoparticles was 211.4 ± 4.031 nm, their PdI value was 0.068 ± 0.016 (Figure 2a), and their zeta potential was −7.87 ± 1.15 mV (Figure 2b).
The PLGA nanoparticle study conducted by Maksimenko et al. reported that the particle size values ranged from 102 to 253 nm, and the zeta potential ranged from −10.1 mV to −11.1 mV [47]. Pereira et al. synthesized PLGA nanoparticles, and they found that the average particle size was 202.5 ± 50.8nm and the PdI was 0.37 ± 0.04 [13]. In a Cymbopogon citratus essential oil-loaded PLGA study, the mean particle size, PdI and zeta potential were found to be 277 nm, 0.18, and -16mV, respectively [19]. Average particle size changes between 204 ± 41.3 and 356 ± 54.9 nm, and PdI value changes between 0.25 ± 0.02 and 0.32 ± 0.04 were reported in a clove bud essential oil-loaded nanoparticle assay [48].
In another study on PLGA, it was reported that the average particle size was 226.9 nm, the PdI value was 0.004, and the zeta value was −7.41 mV [49]. The average particle size, PdI and zeta potential values obtained in our study are compatible with the previous studies.
Next, the encapsulation efficiency and loading capacity of the developed nanoparticles were calculated. The encapsulation efficiency and loading capacity of synthesized nanoparticles are strategic parameters for the development of effective nanoformulation [50]. The encapsulation efficiency was calculated as 59.25% using Equation (1). In a study conducted with anethole-and carvone-loaded PLGA nanoparticles, it was found that the encapsulation efficiency ranged from 2.12% ± 0.098 to 87.31% ± 5.84% [51]. In another study on PLGA, the encapsulation efficiency was calculated to be 50% [52]. In a study conducted with bergamot essential oil, the encapsulation efficiency was found to be between 28% and 84% [53]. The percentage value we obtained in our study is consistent with the values obtained in the literature [54]. The loading efficiency of the LNEO-loaded PLGA nanoparticles was calculated as 25.65% using Equation (2). In a study conducted with PLGA nanoparticles loaded with carvacrol, the encapsulation efficiency and loading capacity were found to be 26% and 21%, respectively. These results are lower than our results [55]. Compared to the study by Fonte et al., it was found that the encapsulation efficiency is lower and the loading capacity is higher [56].
The in vitro release profile of the LNEO from PLGA nanoparticles showed a biphasic release pattern. The initial rapid release can be explained by the rapid release of the active compound near the surface of the nanoparticles [57,58]. However, the 72-h sustained release may be related to the nature of the LNEO trapped in the core of the PLGA nanoparticles [59]. In general, PLGA degradation is slow. Therefore, the active compound release is also related to the diffusion constant, PLGA swelling and surface PLGA erosion [58].
DNA is the primary pharmacological target of anticancer drugs. Therefore, it is very important to determine the interactions between the new anticancer drug molecule and DNA [25]. Therefore, in this study, we demonstrated the potential of LNEO as a potential anticancer molecule by evaluating its interaction with DNA. The UV-Vis spectra results showed that the interaction of the LNEO with CT-DNA resulted in a 6 nm bathochromic (red) shift and a 93.80% hypochromic effect.
In spectral effects, the vacant π* orbital of the molecule pairs with the π* orbital of the DNA base pairs, resulting in an energy reduction and a reduction of the π-π* transition energy. This is detected by the redshift of the absorption in molecular DNA interaction. At the same time, the empty π* orbital is partially filled with electrons in order to reduce the transition probability, resulting in hypochromism [60][61][62][63]. The hypochromic effect and redshift are typical of the binding of the anticancer molecule to DNA in the intercalation mode [63]. Molecules that intercalatively bind to DNA are used in cancer treatment because they inhibit DNA replication in cancer cells [25,64]. Our results showed that LNEO could be an effective anticancer agent for cancer treatment.
In addition to DNA, PI3K/mTOR signaling pathways, which are accepted as the main regulators for cancer, are also important in the pharmacological targeting of anticancer drugs. PI3Ks are members of the intracellular lipid kinases that regulate the cell metabolism, survival, growth and profile, and mTOR is also a class IV PI3K kinase; both are involved in the regulation of cellular growth and proliferation [26]. Compounds that inhibit both PI3Ka and mTOR, which may offer great potential in cancer therapy, are designed by many companies, including Pfizer, and their PI3K/mTOR dual inhibitory activities are currently being investigated for clinical trials. The antiproliferative activity of laurel essential oil was correlated with interactions on this dual inhibitor of PI3K/mTOR as a receptor [26,27] for molecular docking calculations. The molecular docking results indicated that the amino acid Val 882 in the binding pocket of the PI3K/mTOR receptor, which was chosen in the demonstration of the antiproliferative activity of LNEO, has an important place; the compounds in LNEO interacted with this amino acid residue, and even methyleugenol formed two hydrogen bonds with this residue and achieved the most stable binding to the binding pocket. In addition, the four most abundant compounds also showed hydrophobic interactions with the residues Met 953, Ile 831, Ile 879 and Ile 963, with which the molecule interacts with the target receptor and has antiproliferative activity [64,65]. The ADME results showed that the most abundant compounds in LNEO have low molecular weights (complying with Lipinski's Rule of Five), have a high percentage of human oral absorption, and ultimately each have promising pharmacokinetic properties to be a drug candidate. Therefore, laurel essential oil also gives such activity.
The interactions of LNEO with DNA and the PI3K/mTOR dual inhibitor indicate that LNEO-loaded PLGA NPs have the potential to be used in cancer therapy as a novel phytotherapeutic agent-based controlled-release system.  Air-dried Laurus nobilis L. leaves (1000 g) were subjected to hydrodistillation using a Clevenger apparatus, as described in the European Pharmacopoeia, for 3 h in triplicate. The essential oils were recovered and dried over anhydrous sodium sulphate. The samples were kept in brown vials at 4 • C until further analysis.

Analysis of the Essential Oil
Solutions of 10% (v/v) essential oil in n-hexane were subjected to GC-FID/MS analysis. An Agilent 7890B GC-FID (Santa Clara, CA, USA) coupled with an Agilent 5977E electron impact mass spectrometer (Santa Clara, CA, USA) via a two-way capillary splitter was used to identify and quantify the essential oil components. An Agilent G4513A (Santa Clara, CA, USA) auto injector was employed for the injections of 1 µL of the sample solutions. A DB-WAX column (60 m, 0.25 mm, 0.25 µm) was operated with the following temperature program: 70 • C for 15 min, raised to 180 • C at a rate of 2 • C/min. The column temperature was kept isothermal at 180 • C for 5 min, and then increased to 230 • C at a rate of 5 • C/min. Finally, the column temperature was set isothermally at 230 • C for 15 min. The total analysis time was 100 min. Helium was used as a carrier gas with a constant flow of 1.5 mL/min. The split ratio was set to 1:50. The temperatures of the injector port, ion source, quadrupole, MSD transfer line and FID were 250 • C, 230 • C, 150 • C, 250 • C and 220 • C, respectively. The FID air flow was 400 mL/min, and the H 2 flow was adjusted to 30 mL/min. The mass detector scan range was set to 45-450 m/z. The compounds were identified by comparing their spectral data obtained from the Wiley Registry of Mass Spectral Data 9th edition (April 2011) with the NIST 11 Mass Spectral Library (NIST11/2011/EPA/NIH), and by using authentic reference samples. The retention indices were calculated from a co-injected alkane series (C7-C40) compared with previous studies and the NIST online webbook. A quantification was performed by the external standard method, using calibration curves generated by running a GC-FID analysis of representative compounds. All of the analyses were performed in triplicate.

Preparation of the LNEO-NPs
The LNEO-NPs were synthesized using the single-emulsion method [65][66][67][68]. Briefly, 50 mg PLGA was dissolved in 2 mL dichloromethane (DCM). In total, 10 mg LNEO was added to the PLGA solution and mixed on a magnetic stirrer. Then, the LNEO PLGA solution was added dropwise to the 3% PVA solution with an injector. The solution obtained was homogenized by sonication (Bandelin, Sonopuls, Berlin, Germany) for 3 min at 70 W energy in an ice bath. It was stirred on a magnetic stirrer for 16 h in order to remove the solvent from the emulsion. The LNEO-NPs were then purified by applying a three-cycle centrifugation step (at 8000 rpm 40 min). Finally, the LNEO-NPs were filtered through a 0.45 µm cellulose membrane.

Spectrophotometric Analysis of the LNEO
The LNEO analysis was performed using a UV-Vis Spectrophotometer. The maximum absorbance values at 224 nm were obtained for seven different concentrations of the LNEO in ethanol (100, 50, 25, 12.5, 6.25, 3.125, and 1.5625 µg/mL), and the calibration curve of the LNEO was plotted [69]. The calibration curve was used to determine the LNEO's encapsulation efficiency, loading capacity, and in vitro release.

Average Particle Size, Polydispersity Index (PdI) and Zeta Potential Analyses of the LNEO-NPs
The DLS analyses, including the average particle size, PdI and zeta potential of the LNEO-NPs, were performed using a Zeta Sizer Nano ZS device operating at 25 • C, equipped with a 4.0 mV He-Ne laser (633 nm). The 60 nm and 200 nm polystyrene latex particles were used as references to determine the particle size of the LNEO-NPs. Every sample was freshly prepared. The suspension of the LNEO-NPs was diluted with distilled water, and then measured in triplicate.

Scanning Electron Microscopy (SEM)
The morphology of the LNEO-NPs was determined using SEM. The LNEO-NPs were dissolved in ultra water by mixing in an ultrasonic bath for 15 min. They were then prepared by dropping 10 µL of the sample onto glass, and were dried at room temperature for 4 h. The SEM images were obtained at 20.00 KX magnification, 10.00 kV electron high tension, and a 15.0-mm working distance with an in-lens detector.

Determination of the Encapsulation Efficiency and Loading Capacity
The encapsulation yield was achieved by the solvent extraction method [19]. Briefly, LNEO-NPs were dissolved in DMSO, vortexed and centrifuged for 15 min (8000). The supernatants were analyzed by UV-Vis spectrophotometry at 224 nm. Empty NPs were used as reference. The total LNEO content was calculated using the curve equation obtained from the standard curve of the LNEO using the absorbance value obtained from the UV-Vis spectrometer. The encapsulation efficiency (EE) was calculated using Equation (1), and the loading capacity (LC) was calculated using Equation (2 In order to determine the in vitro release profile of the LNEO, 5 mg of the LNEO-NPs were dispersed in distilled water and placed in a dialysis capsule. In total, 75 mL phosphate-buffered saline (60%; pH 7.4) + 40% ethanol were used for the release medium. The samples were incubated at 120 rpm in a horizontal shaking water bath at ambient temperature [69]. The samples were taken at specified time intervals and replaced with an equivalent volume of fresh medium. The samples from the release medium were analyzed using a UV-Vis spectrophotometer. The release amount (%) of the LNEO from the polymeric nanoparticles was calculated using Equation (3).

DNA Binding
Calf Thymus DNA (CT-DNA) was used for the determination of the interaction of the LNEO with DNA [25]. The CT-DNA solution in Tris-HCl/NaCl (pH 7.2) buffer gave a UV absorbance ratio of 1.9 at A 260 /A 280 at 260 and 280 nm wavelengths. This ratio indicates that the CT-DNA does not contain protein [70,71]. The LNEO stock solution was dissolved in ethanol and the required concentration for the experiment, and was prepared by dilution with a buffer. The experiment was performed by keeping the LNEO concentration (40 µg) constant in Tris-HCl/NaCl buffer (pH 7.2) and adding increasing concentrations (0-120 µM) of CT-DNA. After each addition of DNA to the Tris-HCl/NaCl buffer containing LNEO, it was incubated for 5 min, and the absorbance values were recorded.

Molecular Docking and ADME Analysis
In this study, the major components of LNEO from Laurus nobilis L. were chosen as potential ligands. As stated in the introduction, laurel leaf and flower essential oils contain 1,8-cineol, α-terpinyl acetate and methyleugenol as their main components [8]. In addition, the leaf essential oil contains at least 48% 1,8 cineole and 70% α-terpinyl acetate [9]. According to the results obtained by GC-MS, the major components of the LNEO from Laurus nobilis L. were noted as 62.37% 1,8-cineole, 4.48% α-terpinyl acetate, 0.201% methyleugenol, and 9.76% sabinene. The 1,8-Cineole (PubChem CID2758), α-terpinyl acetate (PubChem CID111037), methyleugenol (PubChem CID7127) and sabinene (Pub-Chem CID18818) molecules were preferred as possible ligands for the molecular docking calculations, and their three-dimensional molecular structures were downloaded from the PubChem site [72], in the light of the data obtained from the analysis results of the essential oil components by GC-MS. Molecular structures downloaded as 3D .sdf files were introduced one by one to the Gaussian package program [73], and were optimized using the DFT/B3LYP-631G(d,p) basis set; the optimized geometries were used for the molecular docking analysis. Using the Lig Prep tool in Maestro version 11.4 in the Schrodinger Software program [74][75][76], the ionization and tautomeric states of the optimized molecular structures were created using the Epic process, and the optimization process was performed using the OPLS force field [77]; possible stereoisomers were produced for each ligand. In the molecular docking analysis, the antiproliferative activity and mechanism of laurel essential oil were supported by their mechanism of action on a dual PI3K/mTOR inhibitor. The phosphatidylinositol 4,5-bisphosphate 3-kinase catalytic subunit gamma isoform with 2.7 Å resolution (PDB 4FA6) was chosen as a potential target for molecular docking from the Protein Data Bank (PDB) [26]. All of the water and ions were deleted, polar hydrogens were added, bond orders were assigned, the charges were defined using PROPKA [78] at pH 7.0, and optimization and minimization processes were applied after preprocess analysis via the Protein Preparation Wizard tool. Using the receptor grid tool, the receptor active site to which the ligands bind was obtained by constructing a grid with a cubic box formed by positioning ligands of certain sizes in the center. After the identification of the atomic groups that could rotate in the binding site of the receptor, ligand-receptor docking was performed. In addition, by using the Qik-Prop module [79][80][81] of the same program, the pharmacokinetic and physicochemical properties of the major chemical compounds in laurel essential oil, such as the molecular weights (M W ), percent human oral absorp-tion, estimated octanol/water partition coefficient (QPlogPo/w) and Lipinski Rule of Five compliance of each ligand were evaluated.

Conclusions
In this study, LNEO-NPs were synthesized and characterized. Characterization studies showed that LNEO-NPs have a spherical shape and good distribution in suspension. Moreover, the LNEO-NPs provided the controlled release of LNEO. The usage potential of the LNEO in cancer treatment was investigated using both in vitro and in silico methods.
The antiproliferative activity of the LNEO was achieved by revealing interactions on the dual inhibitor of PI3K/mTOR as a receptor for molecular docking calculations. Based on the strong ADME properties of drug candidate compounds, it may be possible to say that LNEO can be a drug candidate, taking into account the remarkable ADME properties of each of the compounds in the content of laurel essential oil.

Recommendations and Future Works
LNEO is an essential oil that has biological activities: namely antioxidant, antimicrobial, and analgesic activities. Although some components of the LNEOs-such as monoterpenes, sesquiterpenes, oxygenated monoterpenes, and phenolic sesquiterpenesplay an important role in cancer treatment, the usage of LNEO is restricted because it is degraded easily by external factors. This restriction can be overcome using controlled release systems. The approach in our present study might be applied to different phytotherapeutic agents. After in vivo applications are conducted, they might become anticancer drug candidates in the future.  Acknowledgments: The authors thank Anupriya Kumar for the permission to use the Schrödinger Maestro tool for the docking studies. Additionally, in this study, the infrastructure of the Applied Nanotechnology and Antibody Production Laboratory established with TUBITAK support (project numbers: 115S132 and 117S097) was used. The authors would also like to thank TUBITAK for their support.

Conflicts of Interest:
The authors declare no conflict of interest.