key: cord-0069841-l4yulfp0
authors: Pereira, Renato B.; Pinto, Nuno F. S.; Fernandes, Maria José G.; Vieira, Tatiana F.; Rodrigues, Ana Rita O.; Pereira, David M.; Sousa, Sérgio F.; Castanheira, Elisabete M. S.; Fortes, A. Gil; Gonçalves, M. Sameiro T.
title: Amino Alcohols from Eugenol as Potential Semisynthetic Insecticides: Chemical, Biological, and Computational Insights
date: 2021-10-31
journal: Molecules
DOI: 10.3390/molecules26216616
sha: 68c783b930bd088cbe5aae100fa12bc8620d4b2c
doc_id: 69841
cord_uid: l4yulfp0

A series of β-amino alcohols were prepared by the reaction of eugenol epoxide with aliphatic and aromatic amine nucleophiles. The synthesized compounds were fully characterized and evaluated as potential insecticides through the assessment of their biological activity against Sf9 insect cells, compared with a commercial synthetic pesticide (chlorpyrifos, CHPY). Three derivatives bearing a terminal benzene ring, either substituted or unsubstituted, were identified as the most potent molecules, two of them displaying higher toxicity to insect cells than CHPY. In addition, the most promising molecules were able to increase the activity of serine proteases (caspases) pivotal to apoptosis and were more toxic to insect cells than human cells. Structure-based inverted virtual screening and molecular dynamics simulations demonstrate that these molecules likely target acetylcholinesterase and/or the insect odorant-binding proteins and are able to form stable complexes with these proteins. Encapsulation assays in liposomes of DMPG and DPPC/DMPG (1:1) were performed for the most active compound, and high encapsulation efficiencies were obtained. A thermosensitive formulation was achieved with the compound release being more efficient at higher temperatures.

The use of synthetic pesticides for decades to manage pest control in crops has resulted in an accumulation of various residues with adverse effects on many organisms and potential negative impact in human health. At the same time, crop destruction by pests, mainly by insects, is one of the main problems responsible for losses in agricultural production. Pesticides from natural sources are an effective alternative to synthetic pesticides and are becoming more important for pest management in agriculture and also public health. In this respect, plants offer a wide variety of secondary metabolites with efficacy against insects [1, 2] . In recent years, essential oils (EOs) became an important natural source of pesticides. Many EOs present insecticidal, repellent, fumigant, and antifeedant activities against a wide variety of insects [3, 4] . Essential oil components and their derivatives are Eugenol 1 is easily obtained by hydrodistillation from clove, and is known for its various biological activities, as mentioned above, namely insecticidal. Following our recent interests in finding new biopesticides [22] , the present work describes a strategy consisting of structural changes of eugenol in an attempt to obtain semisynthetic alternatives with improved insecticidal activity. Eugenol epoxide 2 was prepared from eugenol through reaction with m-CPBA in DCM using a known procedure [11, 24] in 48% yield. Then, the epoxide was further reacted at 50 • C with a series of aliphatic and aromatic amine nucleophiles in ethanol/water as solvent [25] , which is followed Molecules 2021, 26 , 6616 3 of 20 by column chromatography purification on silica gel using dichloromethane/methanol, mixtures of increasing polarity (3b, 3d-f) as the eluent, or by evaporation of solvents under reduced pressure (3a and 3c) to afford the corresponding β-amino alcohol derivatives 3a-f as oil materials. Thus, the reaction of 2-methoxy-4-(oxiran-2-ylmethyl)phenol 2 with 2-methylpropan-2-amine, octan-2-amine, piperidine, aniline, 3-methoxyaniline, and 4-cyanoaniline gave 4-(3-(tert-butylamino)-2-hydroxypropyl)-2-methoxyphenol 3a, 4-(2-hydroxy-3-(octan-2-ylamino)propyl)-2-methoxyphenol 3b, (4-(2-hydroxy-3-(piperidin-1-yl)-propyl)-2-methoxyphenol 3c, 4-(2-hydroxy-3-(methyl(phenyl)amino)propyl)-2methoxyphenol 3d, 4-(2-hydroxy-3-((3-methoxyphenyl)amino)propyl)-2-methoxy-phenol 3e, and 4-(2-hydroxy-3-(4-hydroxy-3-methoxyphenyl)propyl)amino)benzonitrile 3f, respectively, in 9-97% yield, and they were fully characterized by 1 H and 13 C NMR spectroscopy and HRMS (Scheme 1).

amine nucleophiles in ethanol/water as solvent [25] , which is followed by column chromatography purification on silica gel using dichloromethane/methanol, mixtures of increasing polarity (3b, 3d-f) as the eluent, or by evaporation of solvents under reduced pressure (3a and 3c) to afford the corresponding β-amino alcohol derivatives 3a-f as oil materials. Thus, the reaction of 2-methoxy-4-(oxiran-2-ylmethyl)phenol 2 with 2methylpropan-2-amine, octan-2-amine, piperidine, aniline, 3-methoxyaniline, and 4cyanoaniline gave 4-(3-(tert-butylamino)-2-hydroxypropyl)-2-methoxyphenol 3a, 4-(2hydroxy-3-(octan-2-ylamino)propyl)-2-methoxyphenol 3b, (4-(2-hydroxy-3-(piperi-din-1-yl)-propyl)-2-methoxyphenol 3c, 4-(2-hydroxy-3-(methyl(phenyl)amino)propyl)-2methoxyphenol 3d, 4-(2-hydroxy-3-((3-methoxyphenyl)amino)propyl)-2-methoxyphenol 3e, and 4-(2-hydroxy-3-(4-hydroxy-3-methoxyphenyl)propyl)amino)benzonitrile 3f, respectively, in 9-97% yield, and they were fully characterized by 1 H and 13 C NMR spectroscopy and HRMS (Scheme 1).

The main 1 H-NMR features of compounds 3a-f are the signals for protons of the CHOH, CH2N, and OCH3 groups. The CH2N protons show up as two distinct signals as doublets of doublets or multiplets (δ 3.45-2.53 ppm); the CHOH proton displays as a multiplet in all compounds (δ 4.91-3.89 ppm); and the OCH3 group shows up as a singlet (δ 3.89-3.81 ppm). The tert-butyl group in 3a shows up as a singlet (δ 1.16 ppm), while in 3b and 3c, the N-alkyl chain corresponds to a series of multiplets (δ 2.88-0.88 ppm). The 13 The main 1 H-NMR features of compounds 3a-f are the signals for protons of the CHOH, CH 2 N, and OCH 3 groups. The CH 2 N protons show up as two distinct signals as doublets of doublets or multiplets (δ 3.45-2.53 ppm); the CHOH proton displays as a multiplet in all compounds (δ 4.91-3.89 ppm); and the OCH 3 group shows up as a singlet (δ 3.89-3.81 ppm). The tert-butyl group in 3a shows up as a singlet (δ 1.16 ppm), while in 3b and 3c, the N-alkyl chain corresponds to a series of multiplets (δ 2.88-0.88 ppm). The 13 C main features are the CH 2 N (δ 41.24-41.10 ppm), CHOH carbon (δ 71.1-67.01 ppm), OCH 3 carbon (δ 55.93-55.72 ppm), and the additional OCH 3 in 3e (δ 54.96 ppm). In addition, the tert-butyl carbons in 3a shows up (δ 26.79 and 24.58 ppm), while in 3b, the methyl terminal is highlighted (δ 13.98 ppm).

All molecules obtained were evaluated for their impact in the viability of the Sf9 cells at 100 µg/mL (i.e., 1-6.09 × 10 −4 M; 2-5. 55 CHPY-2.85 × 10 −4 M) by the means of a resazurin-based method. For benchmarking purposes, the insecticide chlorpyrifos was used at the same concentration. As shown in Figure 1 , the only molecule devoid of toxicity was 3c, which incidentally was also the only one bearing a piperidine ring. A second group of molecules, which elicited residual toxicity (under 25% of viability loss), was 1, 2, 3a, and 3b. Eugenol 1 was the starting material, and the results show that the replacement of the terminal methylene group by the epoxide had no effect upon the biological activity of the molecule. Finally, the most potent molecules were 3d, 3e, and 3f, which caused losses of 40%, 30%, and 50% viability in insect cells, respectively. These three molecules were also the only ones bearing a benzene ring next to the nitrogen atom. Considering the unsubstituted ring, 3d, its methoxylation resulted in decreased potency, while the presence of the cyanide group increased it. In light of these results, we decided to advance our studies solely with the two most potent molecules, 3d and 3f, as they were more potent than the benchmark used, chlorpyrifos. 

All molecules obtained were evaluated for their impact in the viability of the Sf9 cells at 100 µg/mL (i.e., 1-6.09×10 −4 M; 2-5.55×10 −4 M; 3a-3.95×10 −4 M; 3b-3.23×10 −4 M; 3c-3.77×10 −4 M; 3d-3.66×10 −4 M; 3e-3.30×10 −4 M; 3f-3.35×10 −4 M; CHPY-2.85×10 −4 M) by the means of a resazurin-based method. For benchmarking purposes, the insecticide chlorpyrifos was used at the same concentration. As shown in Figure 1 , the only molecule devoid of toxicity was 3c, which incidentally was also the only one bearing a piperidine ring. A second group of molecules, which elicited residual toxicity (under 25% of viability loss), was 1, 2, 3a, and 3b. Eugenol 1 was the starting material, and the results show that the replacement of the terminal methylene group by the epoxide had no effect upon the biological activity of the molecule. Finally, the most potent molecules were 3d, 3e, and 3f, which caused losses of 40%, 30%, and 50% viability in insect cells, respectively. These three molecules were also the only ones bearing a benzene ring next to the nitrogen atom. Considering the unsubstituted ring, 3d, its methoxylation resulted in decreased potency, while the presence of the cyanide group increased it. In light of these results, we decided to advance our studies solely with the two most potent molecules, 3d and 3f, as they were more potent than the benchmark used, chlorpyrifos. Figure 1 . Viability of the Sf9 cells after incubation with the presented molecules (100 µg/mL), medium (control), or the reference insecticide chlorpyrifos (CHPY, 100 µg/mL). Cells were incubated for 24 h, after which viability was evaluated. ** p < 0.01, *** p < 0.001.

After establishing the toxicity of the selected molecules toward insect cells, we investigated the mechanism of action behind this effect. In fact, the loss of viability could be a consequence of an array of different biological processes, from necrosis to cell cycle arrest and apoptosis, among others. Necrosis is a process of uncontrolled cell death that encompasses the destruction of cell membranes, with consequent leakage of cytoplasmic content to the surrounding tissues; for this reason, it is usually avoided in biological contexts [26] . To assess the potential unfolding of this event, we assessed the levels of leaked lactate dehydrogenase (LDH) in cells incubated with the selected molecules. Being Viability of the Sf9 cells after incubation with the presented molecules (100 µg/mL), medium (control), or the reference insecticide chlorpyrifos (CHPY, 100 µg/mL). Cells were incubated for 24 h, after which viability was evaluated. ** p < 0.01, *** p < 0.001.

After establishing the toxicity of the selected molecules toward insect cells, we investigated the mechanism of action behind this effect. In fact, the loss of viability could be a consequence of an array of different biological processes, from necrosis to cell cycle arrest and apoptosis, among others. Necrosis is a process of uncontrolled cell death that encompasses the destruction of cell membranes, with consequent leakage of cytoplasmic content to the surrounding tissues; for this reason, it is usually avoided in biological contexts [26] . To assess the potential unfolding of this event, we assessed the levels of leaked lactate dehydrogenase (LDH) in cells incubated with the selected molecules. Being a cytoplasmic enzyme, the finding of extracellular LDH is widely used as a marker of necrosis. As shown in Figure 2A , the incubation of cells with a lysis solution (LS) resulted in a three to four-fold increase in extracellular LDH. Conversely, the incubation of cells with 3d and 3f had no detectable impact in LDH levels in culture media. In light of this, we concluded that the impact of these molecules in the viability of the Sf9 cells was not a consequence of an ongoing necrotic process. Next, we assessed if a process of organized cell death, such as apoptosis, could be taking place. Given the pivotal role of cysteine-aspartic proteases in most forms of apoptosis, we investigated the impact of 3d and 3f in the insect equivalent of mammal caspases, in this case DRACE, using a proluminescent substrate of this target.

As shown in Figure 2B , both 3d and 3f significantly increased the caspase-like activity in treated cells, the latter having a more pronounced effect. This result suggests that both 3d and 3f elicit their cytotoxic effect toward the Sf9 cells by triggering an organized process of cell death with the involvement of cysteine-aspartic proteases. a cytoplasmic enzyme, the finding of extracellular LDH is widely used as a marker of necrosis. As shown in Figure 2A , the incubation of cells with a lysis solution (LS) resulted in a three to four-fold increase in extracellular LDH. Conversely, the incubation of cells with 3d and 3f had no detectable impact in LDH levels in culture media. In light of this, we concluded that the impact of these molecules in the viability of the Sf9 cells was not a consequence of an ongoing necrotic process. Next, we assessed if a process of organized cell death, such as apoptosis, could be taking place. Given the pivotal role of cysteineaspartic proteases in most forms of apoptosis, we investigated the impact of 3d and 3f in the insect equivalent of mammal caspases, in this case DRACE, using a proluminescent substrate of this target. As shown in Figure 2B , both 3d and 3f significantly increased the caspase-like activity in treated cells, the latter having a more pronounced effect. This result suggests that both 3d and 3f elicit their cytotoxic effect toward the Sf9 cells by triggering an organized process of cell death with the involvement of cysteine-aspartic proteases. 

Up to this point, we had already identified two molecules that presented higher potency than the commercial insecticide chlorpyrifos and that were shown the be nonnecrotic and pro-apoptotic. In addition to these traits reported herein, it is also important that prospective new insecticides present some degree of selectivity, specifically low toxicity to human cells. To this end, we assessed the impact of the 3d and 3f in 2D models of human cells. We chose human keratinocytes (HaCaT cell line), as one of the most relevant routes of human contact with pesticides is usually via the skin, where keratinocytes are the first population of living cells in the skin. As shown in Figure 3 , both molecules elicited a weak loss of viability, around 20%. Importantly, both molecules were less toxic than the benchmark chlorpyrifos and, relevantly, they were less toxic to human cells than insect cells.

These results are promising and pave the way for further developments in the field, as the chemical diversity obtained allows drawing some structure-activity relationships, as addressed above. 

Up to this point, we had already identified two molecules that presented higher potency than the commercial insecticide chlorpyrifos and that were shown the be nonnecrotic and pro-apoptotic. In addition to these traits reported herein, it is also important that prospective new insecticides present some degree of selectivity, specifically low toxicity to human cells. To this end, we assessed the impact of the 3d and 3f in 2D models of human cells. We chose human keratinocytes (HaCaT cell line), as one of the most relevant routes of human contact with pesticides is usually via the skin, where keratinocytes are the first population of living cells in the skin. As shown in Figure 3 , both molecules elicited a weak loss of viability, around 20%. Importantly, both molecules were less toxic than the benchmark chlorpyrifos and, relevantly, they were less toxic to human cells than insect cells.

Molecules 2021, 26, 6616 6 of 20 Figure 3 . Viability of HaCaT cells exposed to compounds 3d and 3f (100 µg/mL), medium (control), or the reference insecticide chlorpyrifos (CHPY, 100 µg/mL). Cells were incubated for 24 h, after which viability was evaluated. *** p < 0.001.

After identifying the most promising molecules and the biological processes involved in their cytotoxic effect, we were interested in shedding light on the possible molecular targets. To this end, an array of computational methods was used. Table 1 presents the average scores obtained for compounds 3d and 3f for each potential target with each scoring function. Regarding the different scoring functions, it is Figure 3 . Viability of HaCaT cells exposed to compounds 3d and 3f (100 µg/mL), medium (control), or the reference insecticide chlorpyrifos (CHPY, 100 µg/mL). Cells were incubated for 24 h, after which viability was evaluated. *** p < 0.001. These results are promising and pave the way for further developments in the field, as the chemical diversity obtained allows drawing some structure-activity relationships, as addressed above.

After identifying the most promising molecules and the biological processes involved in their cytotoxic effect, we were interested in shedding light on the possible molecular targets. To this end, an array of computational methods was used. Table 1 presents the average scores obtained for compounds 3d and 3f for each potential target with each scoring function. Regarding the different scoring functions, it is important to mention that they are based on different metrics and scales. The score for all the GOLD scoring functions is dimensionless, and the higher the score, the better the binding affinity. The Vina scoring function, on the other hand, uses a metric that approximates that of binding free energies, so a more negative value means better affinity. The PDB structure with the best score was selected for each potential target, and they were ranked from the best target to worst, according to the predictions of the different docking programs/scoring functions. Globally, considering the results obtained with the different scoring functions, the odorant binding proteins class (OBP) and acetylcholinesterase (AChE) are the protein targets with the highest affinity toward compounds 3d and 3f. This tendency was quite clear with all the different scoring functions evaluated, which further strengthens our results.

To validate the inverted screening results, we evaluate the protein flexibility and characterize the molecular interactions formed, and molecular dynamics simulations were performed for the complexes formed with compound 3d and compound 3f and the two groups of targets predicted at the inverted VS stage: OBP and AChE. Structures with the best score from each group were selected (3K1E for OBP and 1QON for AChE). The stability of AChE: compound 3d, AChE: compound 3f, OBP: compound 3d, and OBP: compound 3f complexes was evaluated using RMSD calculations for the Cα atoms of each complex and for the ligands, Solvent-Accessible Surface Area (SASA) analysis, and hydrogen bonding analysis.

All systems and ligands presented relatively low RMSD values, as seen in Table 2 (and Figure S1) , showing that the target-ligand complexes are well equilibrated and that the eugenol derivatives evaluated maintain their binding conformation predicted from the docking. Table 2 . Average protein and ligand RMSD values (Å), ligand RMSD (Å), average ligand SASA (Å), percentage of potential ligand SASA buried, and average number of ligand-target hydrogen bonds obtained from the MD simulations. ∆G binding energy determined using MM/GBSA and per-residue decomposition, which were calculated for the last 90 ns of the simulation. 

When analyzing the percentage of potential SASA area buried for compound 3d and compound 3f when complexed with AChE and OBP, it can be seen that the two molecules remain tightly bound to the two targets evaluated and well protected from the solvent with average buried areas over 90% (Table 2 ). A small decrease was noticed for compound 3d bound to AChE, in relation to the initial configuration predicted from docking, with an average buried area oscillating between 80 and 90%. These results demonstrate that compounds 3d and 3f remain well bound to the two targets evaluated, even after 100 ns. In particular, the eugenol derivatives evaluated in complex with OBP remain very well protected from the solvent throughout time.

Hydrogen bonding analysis is important to understand the stability of the interactions between the targets and ligands throughout time. The results presented in Table 2 show that both ligands maintain a stable hydrogen bonding profile with the targets evaluated, maintaining between one and three hydrogen bonds with AChE and one and four hydrogen bonds with OBP. Globally, the profile observed shows that compounds 3d and 3f establish more hydrogen bonds with OBP and with AChE. Table 2 summarizes the results discussed so far and presents the values for the Gibbs binding free energy calculated using MM/GBSA. The analysis of the residue contribution to the eugenol derivatives' binding free energy to the two protein targets evaluated highlights the interaction profile of compounds 3d and 3f against AChE and OBP, showing the most important amino acid residues involved in ligand stabilization.

AChE is a serine hydrolase, and it is a very common target for pesticides as it is an enzyme vital for the regulation of acetylcholine in several organisms, from insects to mammals. Since this is an enzyme transversal to many species, the use of anticholinesterase insecticides can cause serious health and environmental problems. In addition, there are reports of insect resistance due to mutation of the AChE gene [27] .

For AChE, compound 3d binding is stabilized mostly by residues Trp83 (−2.4 ± 0.8), Tyr370 (−1.3 ± 0.4), and His480 (−1.3 ± 0.6) through non-polar interactions. For compound 3f, the residues contributing more toward AChE binding are Tyr370 (−2.4 ± 0.1), Tyr372 (−2.5 ± 0.8), and Trp83 (−1.9 ± 0.4), with non-polar interactions playing an important role and π-π stacking with Trp83. Figure 4 represents the average structure of the dominant cluster of the AChE-eugenol derivatives complexes obtained from the analysis of the MD trajectory, illustrating the binding pocket and main interactions formed.

Molecules 2021, 26, 6616 9 respectively. This is consistent with the results presented so far and indicates that eug derivatives have indeed a high affinity toward OBP. In fact, there is a structure depo in the PDB of a bee OBP14 from bound to eugenol [28] . OBP are a large and diverse family of insect proteins. They are involved in transport of hydrophobic odorant and pheromone molecules toward the olfac receptors. They are abundant in the insect family and different in structure but carry similar roles. In the Drosophila melanogaster, there are 52 different types of OBPs a Even though diverse in number and sequence, they present some common features. are small, have six conserved cysteine residues joined by three disulfide bridges, and six alpha-helical domains [29] [30] [31] [32] .

Compounds 3d and 3f have a higher molecular weight than eugenol (273.33 g/ 298.34 g/mol, respectively, versus 164.20 g/mol), but they are also lipophilic and if vol they can in fact be capable of binding OBP. The precise mechanism of action still nee be further validated.

For OBP1, compound 3d binding is stabilized mostly by residues Leu67 (−2.5 ± Trp105 (−2.1 ± 0.4), and Ala79 (−1.7 ± 0.5), through van der Waals interactions compound 3f, the residues contributing more toward OBP binding are Met75 (−2.9 ± and Phe114 (−1.8 ± 0.8) through van der Waals interactions and Trp105 (−2.5 ± 0.4) thro a hydrogen bond with the backbone. Figure 5 illustrates the average structure o dominant cluster of the OBP1 binding pocket and main interactions formed betw compound 3d-OBP1 and compound 3f-OBP1, respectively. For the sake of warranting the potential off-target effect that could pose a toxicity risk, human AChE was also analyzed, the docking scores (Table S1 ) being inferior to the ones obtained for insect AChE, hence suggesting that the eugenol derivatives evaluated favor binding to the insect AChE considered over that of human AChE. When comparing the sequence of amino acids between insect and human AChE, there is only 53-54% sequence identity, even though their 3D structures are very similar. The active-site gorge in the insect enzyme is narrower, and the amino acid residues are different. Moreover, the residues in the opening of the gorge are also different [27] , which might explain the difference in affinity of the eugenol derivatives.

The results show that the most stable complexes are OBP-compound 3d and OBPcompound 3f, with binding free energy values of −31.7 ± 0.2 and −41.6 ± 0.2 kcal/mol, respectively. This is consistent with the results presented so far and indicates that eugenol derivatives have indeed a high affinity toward OBP. In fact, there is a structure deposited in the PDB of a bee OBP14 from bound to eugenol [28] .

OBP are a large and diverse family of insect proteins. They are involved in the transport of hydrophobic odorant and pheromone molecules toward the olfactory receptors. They are abundant in the insect family and different in structure but carry out similar roles. In the Drosophila melanogaster, there are 52 different types of OBPs alone. Even though diverse in number and sequence, they present some common features. They are small, have six conserved cysteine residues joined by three disulfide bridges, and have six alpha-helical domains [29] [30] [31] [32] .

Compounds 3d and 3f have a higher molecular weight than eugenol (273.33 g/mol, 298.34 g/mol, respectively, versus 164.20 g/mol), but they are also lipophilic and if volatile, they can in fact be capable of binding OBP. The precise mechanism of action still needs to be further validated.

For OBP1, compound 3d binding is stabilized mostly by residues Leu67 (−2.5 ± 0.5), Trp105 (−2.1 ± 0.4), and Ala79 (−1.7 ± 0.5), through van der Waals interactions. For compound 3f, the residues contributing more toward OBP binding are Met75 (−2.9 ± 0.4) and Phe114 (−1.8 ± 0.8) through van der Waals interactions and Trp105 (−2.5 ± 0.4) through a hydrogen bond with the backbone. Figure 5 illustrates the average structure of the dominant cluster of the OBP1 binding pocket and main interactions formed between compound 3d-OBP1 and compound 3f-OBP1, respectively. 

The most active compound against the Sf9 cells, compound 3f, was encapsula liposomal systems of the phospholipids DMPG (100%) and DMPG/DPPC (1:1). Th (hydrodynamic diameter), polydispersity index, and zeta potential of the comp loaded liposomes were determined by Dynamic and Electrophoretic Light Scat (Table 3) . These properties can affect the bulk properties, performance, processabilit stability of a nanoformulation. Particularly, the surface charge highly influence stability of the liposomes, and zeta potential values more negative than -30 mV or positive than +30 mV are considered optimal values for good stabilization nanodispersion [33] . In view of this fact, a negatively charged lipid, phosphatidylgly was chosen for the liposomal formulation. Specifically, the phospholipid DMPG phase transition temperature (23 °C) near room temperature [23] , allowing obtaini enhanced release at summer temperatures (around or above 30 °C), where the lipi the fluid (liquid-crystalline) phase. However, the relatively short hydrocarbon ch DMPG and the tendency to form leaky vesicles [34] can hamper a high encapsu efficiency of compound 3f. Therefore, the DPPC/DPPG (1:1) formulation was also t Table 3 . Hydrodynamic diameter (Dh), polydispersity index (PDI), and zeta (ζ) potential of D 

The most active compound against the Sf9 cells, compound 3f, was encapsulated in liposomal systems of the phospholipids DMPG (100%) and DMPG/DPPC (1:1). The size (hydrodynamic diameter), polydispersity index, and zeta potential of the compoundloaded liposomes were determined by Dynamic and Electrophoretic Light Scattering (Table 3) . These properties can affect the bulk properties, performance, processability, and stability of a nanoformulation. Particularly, the surface charge highly influences the stability of the liposomes, and zeta potential values more negative than −30 mV or more positive than +30 mV are considered optimal values for good stabilization of a nanodispersion [33] . In view of this fact, a negatively charged lipid, phosphatidylglycerol, was chosen for the liposomal formulation. Specifically, the phospholipid DMPG has a phase transition temperature (23 • C) near room temperature [23] , allowing obtaining an enhanced release at summer temperatures (around or above 30 • C), where the lipid is in the fluid (liquid-crystalline) phase. However, the relatively short hydrocarbon chain of DMPG and the tendency to form leaky vesicles [34] can hamper a high encapsulation efficiency of compound 3f. Therefore, the DPPC/DPPG (1:1) formulation was also tested. Table 3 . Hydrodynamic diameter (D h ), polydispersity index (PDI), and zeta (ζ) potential of DMPG (100%) and DMPG:DPPC (1:1) liposomes (SD from three independent measurements).

PDI ± SD ζ-Potential ± SD (mV) The values in Table 3 show that both liposome formulations are small in size, with hydrodynamic diameters around 200 nm, presenting also a low polydispersity. A PDI value below 0.3 is considered to be acceptable, indicating a homogenous population of phospholipid vesicles [35] . The zeta potential values indicate a highly negative surface charge, anticipating a low aggregation (due to the electrostatic repulsion) and high colloidal stability.

The compound 3f is a fluorescent molecule ( Figure 6 ) in several solvents and in liposomes. This is a great advantage for the determination of the encapsulation efficiency and drug release, due to the high sensitivity (and selectivity) of fluorescence spectroscopy. The compound 3f is a fluorescent molecule ( Figure 6 ) in several solvents and in liposomes. This is a great advantage for the determination of the encapsulation efficiency and drug release, due to the high sensitivity (and selectivity) of fluorescence spectroscopy.

The encapsulation efficiency (EE%) of compound 3f in liposomes was determined by fluorescence measurements (Table 4 ). In these assays, it was kept in mind that the compound is active against the Sf9 cells at a concentration of 100 µg/mL (3.35 × 10 −4 M) and, therefore, at least this concentration must be encapsulated. The encapsulation efficiencies are high, the system DPPC/DMPG being the most advantageous for 3f encapsulation. Nevertheless, the EE% values show that both formulations are able to encapsulate 3f at concentrations that may guarantee an insecticidal activity (if compound release is effective). Compound release from both liposomal formulations was studied at 20 ºC and 35 ºC, to investigate the temperature dependence of the release profile (Figure 7) . The experimental data were analyzed with the Weibull model (Table S2 and Figure S4 in Supporting Information) and the cumulative concentration released was compared in terms of liposome formulation and temperature. The encapsulation efficiency (EE%) of compound 3f in liposomes was determined by fluorescence measurements (Table 4 ). In these assays, it was kept in mind that the compound is active against the Sf9 cells at a concentration of 100 µg/mL (3.35 × 10 −4 M) and, therefore, at least this concentration must be encapsulated. The encapsulation efficiencies are high, the system DPPC/DMPG being the most advantageous for 3f encapsulation. Nevertheless, the EE% values show that both formulations are able to encapsulate 3f at concentrations that may guarantee an insecticidal activity (if compound release is effective). Compound release from both liposomal formulations was studied at 20 • C and 35 • C, to investigate the temperature dependence of the release profile (Figure 7) . The experimental data were analyzed with the Weibull model (Table S2 and Figure S4 in Supporting Information) and the cumulative concentration released was compared in terms of liposome formulation and temperature. An enhanced release of 3f was observed from liposomes of DMPG, reaching a cumulative release of 62% in 24 h at 35 °C, while, at 20 °C, a 36% release was attained. This is a result of the higher membrane fluidity at 35 °C (above transition temperature) for DMPG liposomes, which provide a thermosensitive formulation. Moreover, at both temperatures, the compound released is higher than 100 µg/mL in 24 h. The rigidity of DPPC at temperatures below its gel to liquid-crystalline phase transition (41 °C) justifies the much lower compound release from DPPC/DMPG liposomes, with 14% and 16% of released compound at 20 °C and 35 °C, respectively, not displaying a significant sensitivity to temperature.

The parameter b of the Weibull model can be related to the release mechanism [36] . If b ≤ 0.75, the release is due to a Fickian diffusion, which is the case for DMPG liposomes at 20 °C (Table S2 ). If b > 1, a complex release mechanism takes place, with multiple mechanisms contributions, which is verified in the other cases.

Dichloromethane, ethanol, methanol, ethyl acetate, light petroleum, and mchloroperbenzoic acid were purchased from Fisher Scientific (Geel, Belgium). The anhydrous magnesium sulfate was PanReac Applichem (Barcelona, Spain) products. Chloroform-d was produced by Eurisotop (Cambridge, England). Thin-layer chromatography (TLC) analyses were carried out on 0.25 mm thick, precoated silica plates (Merck Fertigplatten Kieselgel 60F254, Germany), and spots were visualized under UVlight. Chromatography on silica gel was carried out on Merck Kieselgel (230-240 mesh).

NMR spectra were obtained on a Bruker Avance III (Bruker Corporation, Billerica, An enhanced release of 3f was observed from liposomes of DMPG, reaching a cumulative release of 62% in 24 h at 35 • C, while, at 20 • C, a 36% release was attained. This is a result of the higher membrane fluidity at 35 • C (above transition temperature) for DMPG liposomes, which provide a thermosensitive formulation. Moreover, at both temperatures, the compound released is higher than 100 µg/mL in 24 h. The rigidity of DPPC at temperatures below its gel to liquid-crystalline phase transition (41 • C) justifies the much lower compound release from DPPC/DMPG liposomes, with 14% and 16% of released compound at 20 • C and 35 • C, respectively, not displaying a significant sensitivity to temperature.

The parameter b of the Weibull model can be related to the release mechanism [36] . If b ≤ 0.75, the release is due to a Fickian diffusion, which is the case for DMPG liposomes at 20 • C (Table S2 ). If b > 1, a complex release mechanism takes place, with multiple mechanisms contributions, which is verified in the other cases.

Dichloromethane, ethanol, methanol, ethyl acetate, light petroleum, and m-chloroperbenzoic acid were purchased from Fisher Scientific (Geel, Belgium). The anhydrous magnesium sulfate was PanReac Applichem (Barcelona, Spain) products. Chloroform-d was produced by Eurisotop (Cambridge, England). Thin-layer chromatography (TLC) analyses were carried out on 0.25 mm thick, precoated silica plates (Merck Fertigplatten Kieselgel 60F254, Germany), and spots were visualized under UVlight. Chromatography on silica gel was carried out on Merck Kieselgel (230-240 mesh).

NMR spectra were obtained on a Bruker Avance III (Bruker Corporation, Billerica, MA, USA) at an operating frequency of 400 MHz for 1 H NMR and 100.6 MHz for 13 C NMR using the solvent peak as internal reference at 25 • C. All chemical shifts are given in ppm using δ Me 4 Si = 0 ppm as reference, and J values are given in hertz. Assignments were made by comparison of chemical shifts, peak multiplicities, and J values and were supported by spin decoupling-double resonance and bidimensional heteronuclear correlation techniques. High-resolution mass spectrometry analyses were performed at the "CACTI-Unidade de Masas e Proteómica, at University of Santiago de Compostela", Spain.

A solution of eugenol 1 (0.500 g 3.0 mmol; 1 equiv) dissolved in dichloromethane (18 mL) was added dropwise to a suspension of 70% m-chloroperbenzoic acid (0.750 g; 4.3 mmol; 1 equiv) in dichloromethane (10 mL) at 0 • C. After stirring for 1 h, m-chloroperbenzoic acid was again added (1 equiv), and the reaction mixture was stirred for another 24 h at room temperature. A 10% aqueous solution of sodium sulfate (2 × 20 mL) was added, and the resulting solution was washed with 5% aqueous solution of sodium hydrogen carbonate (2 × 20 mL). The organic phase was dried with anhydrous magnesium sulfate, and the solvent was evaporated to afford the compound 2 as a dark yellow oil 

To a suspension of 2-methoxy-4-(oxiran-2-ilmethyl)phenol 2 (0.163 g; 0.90 mmol; 1 equiv) in H 2 O/EtOH 2:1 (2 mL) was added octan-2-amine (0.502 g; 3.89 mmol) and the resulting mixture was heated at 50 • C for 4 h. Then, water (2 mL) was added, and the resulting mixture extracted with EtOAc (2 mL), The organic phase was collected, dried with anhydrous MgSO 4 , and the solvent evaporated to afford an oil (0.202 g), which was subjected to column chromatography using DCM/MeOH as eluent of increasing polarity to give the compound 3b as a brown oil (0.165 g; 0.53 mmol; 59%). Rf 

To a suspension of 2-methoxy-4-(oxiran-2-ilmethyl)phenol 2 (0.1 g; 0.56 mmol; 1 equiv) in H 2 O/EtOH 2:1 (2 mL) was added piperidine (0.047 g; 0.56 mmol), and the resulting mixture was heated at 50 • C for 5 h. The solvent was evaporated under reduced pressure to afford compound 3c as a brown oil (0.142 g; 0.54 mmol; 97%). Rf = 0.35 (MeOH/DCM 10:90). 1 

To a suspension of 2-methoxy-4-(oxiran-2-ilmethyl)phenol 2 (1 equiv) in H 2 O/EtOH 2:1 (2 mL) was added aniline (0.4 mL; 3.9 equiv), and the resulting mixture was heated at 50 • C for 5.5 h. Then, water (2 mL) was added, and the resulting mixture was extracted with EtOAc (2 mL). The organic phase was collected, dried with anhydrous MgSO 4 , and the solvent was evaporated to afford an oil (0.284 g), which was subjected to column chromatography using DCM/MeOH as eluent of increasing polarity to give compound 3d as a yellow oil (0.095 g; 0.35 mmol; 36%). Rf = 0.7 (MeOH/DCM 5:95). 1 

To a suspension of 2-methoxy-4-(oxiran-2-ilmethyl)phenol 2 (0.162 g; 1 equiv) in H 2 O/EtOH 2:1 (2 mL) was added 3-methoxyaniline (0.544 mg; 4.42 mmol), and the resulting mixture was heated at 50 • C for 4 h. Then, water (2 mL) was added, and the resulting mixture was extracted with DCM (2 mL). The organic phase was collected, dried with anhydrous MgSO 4 , and the solvent was evaporated to afford an oil (0.579 g), which was subjected to column chromatography using light petroleum/EtOAc as an eluent of increasing polarity to give compound 3e as a brown oil (0.098 g; 0.32 mmol; 36%). Rf = 0.35 (ethyl acetate/light petroleum 1:1). 1 

To a suspension of 2-methoxy-4-(oxiran-2-ilmethyl)phenol 2 (0.162 g; 0.90 mmol; 1 equiv) in H 2 O/EtOH 2:1 (2 mL) was added 4-aminobenzonitrile (0.528 mg; 4.47 mmol), and the resulting mixture was heated at 50 • C for 37 h. Then, water (2 mL) was added, and the resulting mixture was extracted with EtOAc (2 mL), The organic phase was collected, dried with anhydrous MgSO4, and the solvent was evaporated to afford an oil (0.265 g), which was subjected to column chromatography using DCM/MeOH as an eluent of increasing polarity to give compound 3f as a dark yellow oil (0.025 g; 0.08 mmol; 9%). Rf 

Insect cells (Sf9, Spodoptera frugiperda) cells were maintained as a suspension culture and cultivated in Grace's medium with 10% FBS and 1% penicillin/streptomycin, at 28 • C with agitation. Cells were used in experiments while in the exponential phase of growth. On the other hand, HaCaT (human keratinocytes) cells were culture in DMEM supplemented with 10% FBS and 1% penicillin/streptomycin at 37 • C in a humidified atmosphere of 5% CO 2 .

For the assessment of viability, a resazurin-based method was used. The Sf9 and HaCaT cells were plated at a density of 3.0 × 10 4 and 1.5 × 10 4 cells/well, respectively, incubated for 24 h, and then exposed to the molecules under study (at 100 µg/mL in Grace's medium) for 24 h. After this period, a commercial solution of resazurin was added (1:10), and the kinetic reaction of fluorescence increase was monitored at 560/590 nm. For HaCaT and the Sf9 cells, 30 and 60 min of incubation were used, respectively.

The Sf9 cells were cultured at the same density described above for the viability assessment. To assess the release of the stable cytosolic enzyme lactate dehydrogenase (LDH) into the media, 24 h after the incubation of cells with the molecules under study (at 100 µg/mL in Grace's medium), 50 µL of culture media were removed to a 96-well plate. The LDH released was determined using a CytoTox 96 ® assay kit (Promega; Madison, WI, USA) according to the manufacturer's instructions. A lysis solution (LS) was used as a positive control to generate a maximum LDH release (45 min). Absorbances were measured at 490 nm in a Multiskan GO plate reader (Thermo Fisher Scientific; Waltham, MA, USA), and results correspond to the fold increase of absorbance in treated vs. untreated cells of four independent experiments performed in duplicate.

The Sf9 cells were plated at the same density described for viability studies and exposed to the molecules under study for the designated time. Generally, the same method described before by some of us was used [37] ; however, it has been adapted toward insect cells, as previously reported [21, 22] . After the incubation period, caspase-3/7 substrate was added to wells, and cells were incubated for 20 min at 22 • C. The luminescent signal was measured in a microplate reader (Cytation™ 3, BioTek, Winooski, VT, USA), and three independent experiments were performed in duplicate. Then, to normalize the results, DNA quantification was performed in a triplicate pool using a Qubit™ 1X dsDNA HS Assay Kit according to a previously described procedure [38] and manufacturer's instructions.

For biological assays, the Shapiro-Wilks normality test was performed in the data to ensure that it followed a normal distribution. Comparison between the means of controls and each experimental condition was performed using one-way ANOVA. Outliers were identified by the Grubbs' test. Data were expressed as the mean ± standard deviation (SD) of at least three independent experiments. GraphPad Prism 7.0 software was used, and values were considered statistically significant with a p < 0.05.

To identify possible molecular targets of the amino alcohols derived from eugenol, an inverted virtual screening protocol was applied. A search on Scopus was performed for papers describing virtual screening (VS) studies involving targets and molecules with insecticidal activity using the keywords: "virtual screening" and "pesticides". The selection criteria were relevance of the target and year of publication. In the 18 studies found, 23 PDB structures were identified and downloaded, enabling the creation of a structural database of putative insecticide targets. These are listed in Table S3 .

The 23 PBD structures were prepared for docking starting using the Autodock Vina plugin for Pymol [39] . Crystallographic waters were removed. Then, the crystallographic ligands were saved in separate files and used as reference for active site coordinates as well as for validation in the re-docking steps. In the absence of ligands, the active site coordinates were based on the most important residues described in the literature. Redocking was used to evaluate the ability of the docking software to reproduce the geometry and orientation of the crystallographic pose, as well as the quality of the docking protocol, and to optimize the docking protocol.

The docking programs/scoring functions used were GOLD [40] (PLP, ASP, Chem-Score, and GoldScore scoring functions), and AutoDock Vina [41] . With each docking program/scoring function, the protocol was optimized for each protein target, to minimize the RMSD in the docking predictions of the reference ligand in re-docking, by comparison with the crystallographic structure of the corresponding complex.

The optimized parameters for each program/scoring function were as follows: Vinadocking box position, docking box dimension, exhaustiveness; GOLD (PLP, ASP, Chem-Score, GoldScore)-binding pocket center, docking region radius, search efficiency, number of runs. The final optimized conditions were used for the subsequent stages. Structures for the two eugenol amino alcohol derivatives with the highest insecticide activity were prepared for docking using Datawarrior [42] and OpenBabel [43] and were docked into each structure with the optimized protocol with all the five scoring functions. A ranked list of most likely targets was prepared based on the average scores obtained for each target with the different scoring functions.

Molecular dynamics simulations were performed using the Amber18 software (University of California, San Francisco, USA) for the two compounds identified from the experimental studies to have the highest insecticide activity (compounds 3d and 3f) , which are bound to the two most promising targets identified from the inverted virtual screening study (odorant binding protein 1-3KIE and acetylcholinesterase-1QON). Since 1QON presented a gap in the structure, a homology model was generated using SWISS-MODEL [44] . A total of 1466 templates were found to match the original sequence, but only the top 50 were used to build the model ( Figure S5 in Supplementary Information) .

Models for the MD simulations were prepared starting from the pose predicted for these complexes in the docking experiments during the inverted virtual screening protocol with GOLD/PLP and treated with the Leap module of AMBER [45] . The protein targets were described with the ff14SB force field [46] , while the eugenol derivatives were parameterized using ANTECHAMBER, with RESP HF/6-31G (d) charges calculated with Gaussian16 [47] and the General Amber Force Field (GAFF) [48] . The overall charge on the system was neutralized through the addition of counter-ions (Na + or Cl − ), and the systems were placed in TIP3P water boxes with a minimum distance of 12 Å between the protein surface and the side of the box.

In order to remove the clashes, the systems were submitted to four consecutive minimizations stages, which were followed by an equilibration and production. In the first four minimization stages, the procedure was applied to (1) water molecules (2500 steps);

(2) hydrogens atoms (2500 steps); (3) side chains of all the amino acid residues (2500 steps); and (4) the full system (10,000 steps). After the complete minimization, the systems were equilibrated by a procedure, which was divided into two stages: in the first stage, NVT ensemble, the systems were gradually heated to 298 K using a Langevin thermostat at constant volume (50 ps); in the second stage, the density of the systems was further equilibrated at 298 K (subsequent 50 ps). Finally, the productions runs were performed during 100 ns. Production was executed with an NPT ensemble at constant temperature (298 K, Langevin thermostat) and pressure (1 bar, Berendsen barostat), with periodic boundary conditions. An integration time of 2.0 fs using the SHAKE algorithm was used to constrain all covalent bonds involving hydrogen atoms. The nonbonded interactions were cut off at 10 Å throughout the entire molecular simulation procedure. The final trajectories were analyzed in terms of RMSD to obtain confirmation that both systems were well equilibrated after the initial 10 ns. The last 90 ns of the simulation were considered for hydrogen bonding analysis, and cluster analysis of the conformations was generated. This overall procedure has been previously used with success in the treatment of several biomolecular systems [49] [50] [51] [52] [53] [54] [55] [56] [57] .

The Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) method [27] was applied to estimate the binding free energies of compounds 3d and 3f to the odorant binding protein 1 and to acetylcholinesterase, considering a salt concentration of 0.100 mol.dm −3 . In addition, the energy decomposition method was employed to estimate the contribution of all the amino acid residues for each of these binding free energies. From each MD trajectory, a total of 1400 conformations taken from the last 70 ns of simulation were considered for the MM-GBSA calculations.

The most active compound against the Sf9 cells, compound 3f, was encapsulated in liposomes composed of the phospholipids 1,2-dimyristoyl-sn-glycero-3-phospho-(1 -racglycerol) (sodium salt) (DMPG) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), which are either composed of DMPG (100%) or DPPC/DMPG (1:1). The liposomes (2 mM total lipid concentration) were prepared by the ethanolic injection method [58] above the transition temperature of each lipid, as previously described [21] . The compound (at an initial concentration of 6.7 × 10 −4 M) was encapsulated by co-injection with the ethanolic lipid solution. The size, polydispersity, and zeta potential of compound-loaded liposomes were measured in a Litesizer 500 Dynamic Light Scattering apparatus from Anton Paar (Anton Paar GmbH, Graz, Austria).

The encapsulation efficiency (EE%) was determined as previously reported [21] and calculated through the Equation (1):

using a calibration curve of fluorescence intensity vs. concentration (C), taking advantage of the fluorescence emission of the compound. The compound release was followed during 24 h at 20 • C and 35 • C, below and above the phase transition temperature of DMPG. The Weibull model (a distribution function) was used to study the transport mechanism involved in the compound release [36] , being expressed in terms of the compound fraction accumulated (m) in solution at time t (Equation (2)):

where a is a scale parameter that defines the timescale of the process, T i represents the latency time of the release process (often being zero), and b is a formal parameter that characterizes the type of curve (b = 1 is exponential; b > 1 is sigmoid, with ascendant curvature delimited by an inflection point; and b < 1 is parabolic, displaying high initial slope and a consistent exponential character).

A series of β-amino alcohols were prepared by reaction of eugenol epoxide with various aliphatic and aromatic amines. The obtained eugenol derivatives were subjected to biological activity evaluation in the Sf9 cell line, in comparison with the corresponding precursors, in order to evaluate their application as potential natural based insecticides.

We identified that the three derivatives bearing a terminal benzene ring, either substituted or unsubstituted, were those showing higher potency, in some cases higher than the benchmark used. We further clarified that the molecules were eliciting their effect by triggering organized cell death, and they were selective for insect cells.

Inverted virtual screening studies with five independent methods suggest that these molecules display their insecticide activity most likely by targeting the insect acetylcholinesterase and/or the insect odorant binding proteins. Molecular dynamics simulations and free energy calculations confirm that these two molecules bind strongly to both targets forming very stable complexes with well-defined molecular interactions that are maintained through time.

Nanoencapsulation studies allow obtaining very reasonable encapsulation efficiencies and a controlled release. Liposomes of DMPG provide a temperature-sensitive compound release, which is more effective than the DPPC/DMPG (1:1) formulation.

The following are available online, Figure S1 : Protein and ligand RMSD (Å) of the AChE and OBP-ligand complexes, Figure S2 : Percentage of the potential solvent accessible surface area of the ligands that is buried by the protein targets evaluated, Figure S3 : Number of ligand-target hydrogen bonds formed during the simulations for compound 3d and 3f when complexed with AchE and OBP, Figure S4 : Fitting of the release profiles to the Weibull model, Figure  S5 : Homology model built for 1QON, Table S1 : Docking scores for compounds 3d and 3f in complex with human and insect AChE, Table S2 : Parameters of the Weibull model for the release of compound 3f from liposomes and corresponding coefficients of determination (R 2 ), Table S3 : List of targets selected for the inverted virtual screening study.

Plant extracts for developing mosquito larvicides: From laboratory to the field, with insights on the modes of action

Repellent and acaricidal activity of essential oils and their components against Rhipicephalus ticks in cattle

Green pesticides: Essential oils as biopesticides in insect-pest management

A review on prospects of essential oils as biopesticide in insect-pest management

Syzygium aromaticum): A precious spice

Antimicrobial activity of eugenol derivatives

Eugenol confers resistance to tomato yellow leaf curl virus (TYLCV) by regulating the expression of SlPer1 in tomato plants

Ortho-eugenol exhibits anti-nociceptive and anti-inflammatory activities

Eugenol derivatives as future potential drugs

Molecular and pharmacological analysis of an octopamine receptor from American cockroach and fruit fly in response to plant essential oils

Eugenol derivatives: Synthesis, characterization, and evaluation of antibacterial and antioxidant activities

re-and si-Face-selective nitroaldol reactions catalyzed by a rigid chiral quaternary ammonium salt: A highly stereoselective synthesis of the HIV protease inhibitor amprenavir (Vertex 478)

1,2-Amino alcohols and their heterocyclic derivatives as chiral auxiliaries in asymmetric synthesis

Sharpless asymmetric aminohydroxylation: Scope, limitations, and use in synthesis

Catalytic asymmetric aminohydroxylation (AA) of olefins

Synthesis, antibacterial and antitubercular evaluation of cardanol and glycerol based β-amino alcohol derivatives

Synthesis of β-amino alcohols using the tandem reduction and ring-opening reaction of nitroarenes and epoxides

Discovery of novel bis-oxazolidinone compounds as potential potent and selective antitubercular agents

One-pot synthesis of N-substituted β-amino alcohols from aldehydes and isocyanides

Cytotoxic plant extracts towards insect cells: Bioactivity and nanoencapsulation studies for application as biopesticides

New eugenol derivatives with enhanced insecticidal activity

Thermotropic phase transitions of pure lipids in model membranes and their modifications by membrane proteins

Antibacterial and antioxidant photoinitiated epoxy co-networks of resorcinol and eugenol derivatives

Highly chemoselective addition of amines to epoxides in water

Apoptosis, pyroptosis, and necrosis: Mechanistic description of dead and dying eukaryotic cells

Can acetylcholinesterase serve as a target for developing more selective insecticides?

Crystal structure of Apis mellifera OBP14, a C-minus odorant-binding protein, and its complexes with odorant molecules

Virtual screening of plant volatile compounds reveals a high affinity of Hylamorpha elegans (Coleoptera: Scarabaeidae) odorant-binding proteins for sesquiterpenes from Its native host

The diverse small proteins called odorant-binding proteins

Beyond chemoreception: Diverse tasks of soluble olfactory proteins in insects

Structure and biotechnological applications of odorantbinding proteins

Importance of physicochemical characterization of nanoparticles in pharmaceutical product development. In Advances in Pharmaceutical Product Development and Research: Basic Fundamentals of Drug Delivery

Aqueous dispersions of DMPG in low salt contain leaky vesicles

Impact of particle size and polydispersity index on the clinical applications of lipidic nanocarrier systems

On the use of the Weibull function for the discernment of drug release mechanisms

New chalcone-type compounds and 2-pyrazoline derivatives: Synthesis and caspase-dependent anticancer activity

Bolaamphiphilic bis-dehydropeptide hydrogels as potential drug release systems

Ligand docking and binding site analysis with PyMOL and Autodock/Vina

Development and validation of a genetic algorithm for flexible docking

AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading

DataWarrior: An open-source program for chemistry aware data visualization and analysis

Open Babel: An open chemical toolbox

SWISS-MODEL: Homology modelling of protein structures and complexes

Development and testing of a general amber force field

Improving the accuracy of protein side chain and backbone parameters from ff99SB

Automatic atom type and bond type perception in molecular mechanical calculations

In silico, in vitro and in cellulo models for monitoring SARS-CoV-2 Spike/Human ACE2 complex, viral entry and cell fusion

Identification of new potential inhibitors of quorum sensing through a specialized multi-level computational approach

Molecular dynamic simulations and structure-based pharmacophore development for farnesyltransferase inhibitors discovery

Enzymatic "tricks": Carboxylate shift and sulfur shift

Molecular dynamics analysis of a series of 22 potential farnesyltransferase substrates containing a CaaX-motif

Chemically modified tetracyclines as inhibitors of MMP-2 matrix metalloproteinase: A molecular and structural study

Molecular dynamics simulations on the critical states of the farnesyltransferase enzyme

Enzyme flexibility and the catalytic mechanism of farnesyltransferase: Targeting the relation

End-point binding free energy calculation with MM/PBSA and MM/GBSA: Strategies and applications in drug design

Ethanol injection method for hydrophilic and lipophilic drug-loaded liposome preparation

The authors declare no conflict of interest.Sample Availability: Not applicable.