In vitro and in silico studies of terpenes, terpenoids and related compounds with larvicidal and pupaecidal activity against Culex quinquefasciatus Say (Diptera: Culicidae)
BackgroundIn order to develop new larvicidal agents derived from phytochemicals, the larvicidal activity of fifty molecules that are constituent of essential oils was evaluated against Culex quinquefasciatus Say. Terpenes, terpenoids and phenylpropanoids molecules were included in the in vitro evaluation, and QSAR models using genetic algorithms were built to identify molecular and structural properties of biological interest. Further, to obtain structural details on the possible mechanism of action, selected compounds were submitted to docking studies on sterol carrier protein-2 (SCP-2) as possible target.ResultsResults showed high larvicidal activity of carvacrol and thymol on the third and fourth larval stage with a median lethal concentration (LC50) of 5.5 and 11.1 µg/mL respectively. Myrcene and carvacrol were highly toxic for pupae, with LC50 values of 31.8 and 53.2 µg/mL. Structure–activity models showed that the structural property π-bonds is the largest contributor of larvicidal activity while ketone groups should be avoided. Similarly, property–activity models attributed to the molecular descriptor LogP the most contribution to larvicidal activity, followed by the absolute total charge (Qtot) and molar refractivity (AMR). The models were statistically significant; thus the information contributes to the design of new larvicidal agents. Docking studies show that all molecules tested have the ability to interact with the SCP-2 protein, wherein α-humulene and β-caryophyllene were the compounds with higher binding energy.ConclusionsThe description of the molecular properties and the structural characteristics responsible for larvicidal activity of the tested compounds were used for the development of mathematical models of structure–activity relationship. The identification of molecular and structural descriptors, as well as studies of molecular docking on the SCP-2 protein, provide insight on the mechanism of action of the active molecules, and the information can be used for the design of new structures for synthesis as potential new larvicidal agents.Electronic supplementary materialThe online version of this article (10.1186/s13065-018-0425-2) contains supplementary material, which is available to authorized users.
Introduction<!><!>Insect cultures and rearing conditions<!>Larvicidal activity bioassays and statistical analysis<!>DFT study and descriptors calculations<!>Structure–property–larvicidal activity models<!>Molecular docking studies on protein SCP-2<!><!>Larvicidal activity and quantitative structure–larvicidal activity relationship<!><!>Larvicidal activity and quantitative structure–larvicidal activity relationship<!><!>Quantitative property–larvicidal activity relationship and DFT study<!><!>Quantitative property–larvicidal activity relationship and DFT study<!><!>Docking studies on sterol carrier protein-2 (SCP-2)<!><!>Docking studies on sterol carrier protein-2 (SCP-2)<!><!>Docking studies on sterol carrier protein-2 (SCP-2)<!>Conclusions<!>
<p>More than half of the global human population is exposed to the risk of infection spread by mosquitoes; including Culex spp., Anopheles spp. and Aedes spp. that are considered a public health problem, sin are vectors of pathogenic parasites. Lymphatic filariasis uses Culex quinquefasciatus Say (Diptera: Culicidae) as vector; it is one of the leading causes of global morbidity, with close to 150 million infected, especially in tropical climates [1]. Culex quinquefasciatus is present in most tropical regions of the world; it is commonly found in many urban areas and has been reported as resistant to registered insecticides [2].</p><p>The control of mosquito larvae and pupae currently relies on the use of synthetic chemical insecticides [3]. However, prolonged use of these synthetic pesticides has caused numerous problems, such as the development of resistance [4], undesirable effects on non-target organisms, effects on wildlife, damage to human health and other negative impacts on the environment [5–7]. Several studies have searched for natural products derived from plants as possible mosquito control environmentally-friendly strategy; reports include the larvicidal action of essential oils (EOs) and their constituents [8, 9]. EOs can be alternative pest control agents, because some of their compounds have proven to be highly selective, easily removable, biodegradable, with low or no toxicity against mammals and are effective against a full spectrum of mosquito pests [10, 11]. Also EOs are characterized by reduced effects on non target organisms and minimal environmental persistence [12]. With few exceptions, some of the purified terpenoid constituents of EOs are moderately toxic to mammals, but the oils themselves or their compounds are mostly non toxic to mammals, birds, and fish [12].</p><p>EOs are heterogeneous mixtures of organic chemical compounds [13] mainly terpenoids and phenylpropanes, but low molecular weight aliphatic compounds, acyclic esters or lactones may also be present [14]. The EOs chemical composition is affected by diverse factors, including plant species and subspecies, geographical location, harvested time, the part of the plant used and the extraction methods employed to obtain the EO [15]. In spite of several studies on the larvicidal activity of EOs and their constituents, little is known on the mechanism of action exerted by terpenoids and phenylpropanoids on mosquito larvae. This has motivated the study of the molecular properties, reactivity or structural modulation of essential oil chemical components in order to minimize synthetic and biological evaluation effort for the development of new compounds with potential larvicidal activity.</p><p>Computer assisted prediction of the biological activity of specific chemical compounds considering their chemical structure is now a common technique used in drug discovery [16, 17]. Quantitative structure–activity relationship (QSAR) and quantitative property–activity relationship (QPAR) studies can provide information to understand the relationship between molecule's chemical structure and biological activity [18]. Also, molecular docking is an in silico technique used to estimate the strength of the protein–ligand interaction, to determine biding poses and free energy values [19]. Docking describe ligand binding to a receptor through noncovalent interactions which is commonly used to explore the ligand recognition on targets for new drug development [20].</p><p>This article describes the larvicidal activity of fifty compounds against larvae and pupae of Culex quinquefasciatus (Diptera: Culicidae). Terpenes, terpenoids and others related compounds constituents of different EOs were evaluated in this work. Likewise, the present work reports the theoretical characterization of the molecular and electronic properties of experimentally tested molecules. QSAR/QPAR models and docking studies are also included to emphasize the molecular and structural properties that are essential in the larvicidal activity.</p><!><p>(1) p-Anisaldehyde, (2) Canphor, (3) (3) Carene, (4) Carvacrol, (5) Carveol, (6) Carvomenthol, (7) Carvone, (8) Carvotanacetol, (9) β-Caryophyllene, (10) Citronellal, (11) β-Citronellol, (12) m-Cresol, (13) o-Cresol, (14) Cuminaldehyde, (15) p-Cimene, (16) t-Dihydrocarvone, (17) 3,4-Dimethylcumene, (18) Eucalyptol, (19) Geranial, (20) Geraniol, (21) Germacrene-D, (22) α-Humulene, (23) Hydrocarvone, (24) Hydrodihydrocarvone, (25) 3-Isopropylphenol, (26) Isoborneol, (27) Isopulegol, (28) t-Isopulegone, (29) Lavandullol, (30) Limonene, (31) Linalool, (32) Menthol, (33) Menthone, (34) Myrcene, (35) Neoisopulegol, (36) Perillaldehyde, (37) β-Phellandrene, (38) α-Pinene, (39) β-Pinene, (40) Pulegone, (41) Rotundifolone, (42) Sabinene (43) α-Terpinene, (44) γ-Terpinene, (45) 4-Terpineol, (46) α-Terpineol, 47) β-Terpineol, 48) γ-Terpineol, (49) Terpinolene, (50) Thymol</p><!><p>Larvae of Cx. quinquefasciatus were collected from water tanks in the Sanctorum Cemetery in Mexico City, Mexico (19°27′17″N, 99°12′47″W) and identified using Harwood and James descriptions [21]. Groups of 50 individuals of first and second instar larvae were placed in glass bottles with purified water, maintained at 26 ± 2° C with a natural photoperiod and supplied with 3:1 powdered mixture of dog food and baking powder. The third instar emerging larvae were then separated by groups of 10 individuals in 100 mL tubes with distilled water [22].</p><!><p>Bioassays were done according to the World Health Organization (WHO) protocol with few modifications [23]. Third and fourth instar larvae as well as pupae, were used for testing. Five groups of 20 larvae were isolated in beakers of 250 mL, exposed to different concentrations of the tested compounds and maintained in starvation throughout the experimental period; the surviving larvae were counted in order to record larval mortality. The compounds were diluted in dimethyl sulfoxide (DMSO) (Sigma, 472301) before being added to the aqueous medium which contained the larvae. Temephos H at 0.1 ppm (commercial concentration) was used as a standard for comparison. Larvae were considered dead if they were immobile and unable to reach the water surface [24]. Lethal concentrations (LC50) was calculated using Probit analysis. Data were processed using MS Excel 2010 and SAS v. 9 (Proc Probit) computer programs.</p><!><p>Computational studies were carried out using the Spartan 03 [25] and Gaussian 09 quantum chemistry computer programs [26]. The molecular structures were analyzed by a conformational analysis of each molecule in gas phase using the mechanics force field SYBYL [27]. The minimum energy conformation was selected in order to obtain the geometry optimization using the density functional theory (DFT). The equilibrium geometries of the molecules in the electronic ground state were determined with the Becke three-parameter hybrid functional combined with Lee–Yang–Parr correlation functional (B3LYP) [28, 29]. The basis set 6-311G(d,p) was used for the geometry optimization and vibrational frequency calculations and the 6-311+G(d,p) was applied for vertical excitation energy calculations [30–32]. Analytical frequency calculations were carried out, where the absence of imaginary frequencies confirmed that the stationary points correspond to the global minima of the potential energy hypersurfaces.</p><p>The Koopmans theorem [33] was applied for calculations of the chemical reactivity descriptors such as: the ionization potential (I), electron affinity (A), electronegativity (χ), chemical potential (μ), hardness (ɳ), softness (σ), global electrophilicity (ω), as well as the electronic parameters of, EHOMO (energy of highest occupied molecular orbital), ELUMO (energy of the lowest unoccupied molecular orbital) and band gap (GAPE) were calculated. All molecules were analyzed in the gas and aqueous phase. The polarizable continuum model (PCM) was used to model the solvent effects [34].</p><p>Structure, constitutional, physicochemical and topological descriptors were generated using Dragon 5.0 software [35] using the optimized structure in the aqueous phase.</p><!><p>QSAR/QPAR studies was carried out using all biological activities obtained in vitro and the calculated theoretical descriptors; the analysis was carried out using genetic algorithms with the Mobydigs Software [36]. The quality of the model was considered statistically satisfactory based on the determination coefficient (R2), leave-one-out cross-validated explained variance (Q2), standard deviation (s) and the ANOVA (F) of the model.</p><!><p>The sequence of sterol carrier protein (SCP-2) of Cx. quinquefasciatus (GenBank: AAO43438.1) was obtained from the database of the National Center for Biotechnology Information (NCBI). The protein was modeled through Swiss-Model server [37, 38], using as template the sterol carrier protein of Aedes aegypti (PDB: 1PZ4) [39] reported in the RCSB Protein Data Bank. The final model was subjected to Ramachandran analysis using the Rampage server [40]. Docking analysis was done using the AutoDock4 software [41]. For the docking the active site was defined considering the residues within a grid of 60 A° × 60 A° × 60 A° centered in the active site, with an initial population of 100 randomly placed individuals and a maximum number of 1.0 × 107 energy evaluations. Active site was determined under the description made by Dyer et al. [39]. Compounds for docking were drawn in Gauss view before docking, the compounds were subjected to energy minimization using the hybrid functional B3LYP with a 6, 311G(d,p) basis set. The Kd and ΔG (Kcal/mol) values were obtained from the conformation with the lowest minimum free energy of the ligand coupled on the protein targets. The figures were prepared with ChemBioOffice [42] for the structures and Chimera [43] for the proteins and ligands.</p><!><p>Larvicidal activity of the terpenes, terpenoids and related compounds against Cx. quinquefasciatus</p><p>In parenthesis, 95% confidence intervals, compounds activity is considered significantly different when the 95% CI fail to overlap</p><!><p>EOs are aromatic extracts obtained from plant material that are complex mixtures of volatile secondary metabolites [44]. Some of the compounds present in EOs are terpenes (molecules formed of isoprene units) [45], terpenoids (terpenes with oxygen on its structure) [45] and phenylpropanoids [47]. In the present report, carvacrol and thymol (terpenoids found mainly in the EO of oregano) were the most active molecules with a LC50 of 7.7 and 8.4 μg/mL respectively, against larvae at fourth stage. Myrcene presented a relevant activity against pupae with a LC50 of 31.8 μg/mL. Cheng et al. reported the results of screening EOs and suggested that oils with LC50 values >  100 ppm should not be considered active, whereas those with LC50 values < 50 ppm could be regarded as highly active [48]. Our results agree with reports of the larvicidal activity of constituents of oregano EO; the reports demonstrate that these compounds have fumigant and repellent activity [49–53].</p><p>In relation to chemical structure and larvicidal activity, results have been grouped considering the main chemical moiety of the tested compounds in monocyclic-terpenes, monocyclic-terpenoids, bicyclic-terpenes and bicyclic-terpenoids, and phenylpropanes. β-Caryophyllene, a bicyclic sesquiterpene, showed the lower larvicidal activity with a LC50 of 57.7 μg/mL against fourth instar and 222.3 μg/mL against pupae, Doria et al. also report low larvicidal activity of β-caryophyllene against Aedes aegypti [54]. Sabineno, a bicyclic monoterpene, also had a low activity, with LC50 values of 59.0 μg/mL for fourth instar and 258 μg/mL against pupae. β-Pinene and 3-carene presented a LC50 of 19.6 and 24.7 μg/mL respectively against the fourth stair being the most active of the bicyclic terpenes. Eucalyptol was the bicyclic terpenoid most active against pupae, the only activity lower than 100 μg/mL of all bicyclic compounds evaluated.</p><!><p>Summary of the statistics quantitative structure–larvicidal activity relationship models for activity against fourth instar of Cx. quinquefasciatus</p><p>n, number of systems evaluated; Q2, the square of the coefficient of cross-validation; R2, the square of the correlation coefficient; s, standard deviation; F, Fisher statistic; WC, without contribution; nCt, number of total tertiary C (sp3); nCconj, number of non-aromatic conjugated C (sp2); nR = Cp, number of terminal primary C (sp2); nRCO, number of ketones (aliphatic); nROR, number of ethers; nArOH, number of phenolic groups; nOH, number of a hydroxyls</p><!><p>Models 3 and 4 (Table 2) were constructed based on the larvicidal activity of 47 evaluated molecules, excluding the sesquiterpenes β-caryophyllene (9), germacrene (21) and α-humulene (22) from the analysis. The models showed the same relationship with the nCconj, nRCO, nROR descriptors and the number of phenolic groups (nArOH) and the number of hydroxyl groups (nOH) as descriptors directly related to the biological activity. This is consistent with the most biologically-active molecules: carvacrol and thymol. In monocyclic terpenoids and monocyclic terpenes, increasing the number of double bonds also increased the larvicidal activity. Menthol has a LC50 of 38.1 μg/mL against fourth instar larvae, while thymol had an activity of 12.2 μg/mL. The structural difference between these two compounds is the phenolic group in thymol as compared to menthol that only has the hydroxyl group; p-Cymene has the benzene group without hydroxyl group with an activity of 24.0 μg/mL; this demonstrate the importance of the phenolic group in the larvicidal activity. Carvacrol, an isomer of thymol, has a LC50 of 7.7 μg/mL; therefore, the position of the hydroxyl group plays an important role in the larvicidal activity.</p><p>For acyclic terpenes and terpenoids, higher larvicidal activity was observed in compounds with a higher number of double bonds and increased lipophilicity. Ketone acyclic terpenes were the compounds with lowest larvicidal activity; substitution of the ketone group by the hydroxyl group increased the biological activity considerably. Citronellol was the alcohol terpene with lower activity against the fourth instar and pupae. Geraniol has one double bond more than citronellol, and this structural difference increase the larvicidal activity. Linalool also presents one double bond more than citronellol, and this differential structure is reflected in an increase in larvicidal activity, however the position of the hydroxyl group changes from a primary to a secondary alcohol; this difference could be responsible for the lower biological activity shown. On the other hand, myrcene, an acyclic molecule with no oxygen in its structure, has the highest larvicidal activity and is the only compund with significant activity against pupae, with a LD50 of 19.1 μg/mL against the fourth instar larvae and 31.8 μg/mL against pupae. Myrcene has three double bonds in its structure, and since the lipophilicity is increased in the absence of oxygen, these is an important trait for their potential activity. Accordingly, Lucia et al. consider that the octanol water partition coefficient (LogP) is an important molecular property in the larvicidal activity of monoterpenes [55].</p><!><p>Structural descriptors calculated</p><p>nCs, Number of total secondary C (sp3); nCt, number of total tertiary C (sp3); nCconj, Number of non-aromatic conjugated C (sp2); nR = Cp, number of terminal primary C (sp2); nR = Cs, number of aliphatic secondary C(sp2); nR = Ct, number of aliphatic tertiary C(sp2); nRCO, number of ketones (aliphatic); nArOH, number of aromatic hydroxyls; nOH, number of a hydroxyls; nHDon, number of donor atoms for H-bonds; nHAcc, number of acceptor atoms for H-bonds</p><p>Predicted versus experimental larvicidal activity from structural–activity relationship models. a Model 1, b model 3, c model 5</p><!><p>The models that describe the relationship between the molecular properties and biological activity demonstrate that the octanol–water partition coefficient (MlogP) descriptor is the largest contributor to the larvicidal activity. Lucia et al. developed a QSAR model based on six monoterpenes and they found that when vapor pressure and lipophilicity values decreased, the larvicidal activity against A. aegypti also diminished. The strong effect of the octanol–water partition coefficient can be explained considering that the main conduit for component entrance to the organism is tactile (external cuticle) [55]. Therefore, the partition occurs between the hydrophilic environment (water) and the lipophilic environment (larval epicuticle); therefore, molecule hydrophobicity plays an important role in the intoxication of the larva [56].</p><!><p>Summary of the statistics quantitative property–larvicidal activity relationship models for activity against fourth instar of Cx. quinquefasciatus</p><p>n, number of systems evaluated; Q2, the square of the coefficient of cross-validation; R2, the square of the correlation coefficient; s, standard deviation; F, Fisher statistic; WC, without contribution; J, Balaban-like index; MlogP, Moriguchi octanol–water partition coeff. (logP); TIE, E-state topological parameter; AMR, Ghose–Crippen molar refractivity; Qtot, total absolute charge; BAC, Balaban centric index; Hy, hydrophilic factor; η, chemical hardness; EHOMO, energy of the HOMO orbital; m, dipole moment</p><p>Predicted versus experimental larvicidal activity from property–activity relationship models. a Model 1, b model 3, c model 5</p><!><p>The lipophilic character of terpenes and their derivatives have been widely discussed as a key factor in the antimicrobial and larvicidal properties of these compounds [14–16, 44, 45]; however, it does not finish describing their larvicidal behavior. Sesquiterpenes, for example, have high MlogP values and are not the most active compounds.</p><!><p>Molecular and physicochemical descriptors calculated</p><p>Qpos, total positive charge; Qneg, total negative charge; Qtot, total absolute charge (electronic charge index-ECI); Ui, unsaturation index; Hy, hydrophilic factor; AMR, Ghose–Crippen molar refractivity; TPSA, topological polar surface area; MlogP, Moriguchi octanol–water partition coeff.; AlogP, Ghose–Crippen octanol–water partition coeff</p><p>Chemical reactivity descriptors calculated</p><p>EHOMO, energy of the HOMO orbital; ELUMO, energy of the LUMO orbital; GapE, ELUMO–EHOMO; I, ionization potential; A, electron affinity; μ, chemical potential; χ, electronegativity; ƞ, Chemical hardness σ, chemical softness; m, dipole moment</p><p>The contour plots of LUMO orbitals of the most active molecules. (1) p-Anisaldehyde, (2) Canphor, (3) 3-Carene, (4) Carvacrol, (9) β-Caryophyllene, (10) Citronellal, (11) β-Citronellol, (12) m-Cresol, (13) o-Cresol, (14) Cuminaldehyde, (15) p-Cimene, (17) 3,4-Dimethylcumene, (18) Eucalyptol, (19) Geranial, (20) Geraniol</p><p>The contour plots of LUMO orbitals of the most active molecules (conti…). (21) Germacrene-D, (22) α-Humulene, (25) 3-Isopropylphenol, 26) Isoborneol, (29) Lavandullol, (30) Limonene, (31) Linalool, (34) Myrcene, (38) α-Pinene, (41) Rotundifolone, (42) Sabinene, (43) α-Terpinene, (46) α-Terpineol, (48) γ-Terpineol, (49) Terpinolene, (50) Thymol</p><p>Topological descriptors calculated</p><p>J, Balaban-like index; TIE, E-state topological parameter; UNIP, unipolarity; CENT, centralization; VAR, variation; BAC, Balaban centric index; LOP, lopping centric index; ICR, radial centric; CSI, eccentric connectivity index; ECC, eccentricity; PHI, Kier flexibility index</p><!><p>The mechanism of action of the larvicidal and repellent activity exerted by EOs and their constituents is not fully described. Inhibition of the acetylcholinesterase (AChE) enzyme has been frequently proposed, a similar neurotoxic effect produced by organophosphorus and carbamate incesticides [69, 70]. Similar results have been reported when flies and cockroaches are exposed to eugenol and α-terpineol [71]. However, some authors agree that in most cases there is no relationship between inhibition of AChE and larvicidal effects of terpenes and derivatives [72, 73].</p><p>Priestley et al. proposes that EOs and their constituents act on GABA receptors, as indicated by their results when exposing Drosophila melanogaster to thymol [74]. In addition, Kumar et al. have reported that terpenes present in Calotropis gigantea have larvicidal activity due to the ability to block the sterol carrier protein (AeSCP-2) [75], which is partially responsible for intracellular cholesterol transport in insects [76]. The larvaes during the feeding step contain high concentrations of SCP-2 because they depend on exogenous sources of cholesterol for biosynthesis of steroid derivatives [77]. Therefore, compounds that can inhibit this protein have a high potential as vector control agents.</p><!><p>Results of the construction homology of the sterol carrier protein (SCP-2). a Model of sterol carrier protein (SCP-2). b Ramachandran plot corresponding for the model of SCP-2</p><p>Docking results by SCP-2 from Culex quinquefasciatus</p><p>* Hydrogen bonds interaction</p><!><p>Results showed that monoterpenes and monoterpenoids with the highest larvicidal activity were also the compounds with better binding energy values, being carvacrol the most active followed by α-terpinene and terpinolene. Another important observation is that monoterpenes and monoterpenoides with the highest larvicidal activity are capable of interact with the Phe105 residue.</p><!><p>Interaction of cyclic terpenes and terpenoids with SCP-2. a Interaction of carvacrol (green), α-Terpinene (orange) and terpinolene (yellow). b Interaction of limonene (brown) and carvona (green). c Interaction of p-cimene (yellow), menthol (pink) and thymol (purple)</p><!><p>The relevance of the phenol group is observed when the binding energies of cymene, menthol, thymol and carvacrol are compared. Cymene binding energy is − 5.28 kcal, while menthol is − 5.69 kcal, this energy difference can be attributed to the hydroxyl group; on the other hand, thymol has a binding energy of − 6.66 kcal, which shows that the phenolic group is also important. This characteristics are also observed when comparing the bonding energy of carvomenthol and carvacrol. These results are consistent with the QSAR models also included in this work. The structural difference between the aromatic ring present in thymol and menthol without π bonds, generates a change in the arrangement of the later in the SCP-2 protein active site. It can be observed that the larger aliphatic chain in para position of cymene and thymol is in the direction of Phe105 residue, but does not interact with it, while the menthol is in the opposite position; however the hydroxyl group is kept in the same coordinates as for thymol (Fig. 7c). This is because the hydroxyl group of thymol and menthol are capable of forming hydrogen bonds with the amino group of Arg24 residue.</p><p>The position of the hydroxyl group in the phenolic group is also relevant. The hydroxyl group in the meta position of carvacrol leaves more exposed to larger aliphatic chain, which interacts with the Phe105 residue; the results is an increased biological activity as well as a more favorable binding energy as compared to thymol. The hydroxyl group of carvacrol can form hydrogen bonds with the amino group in Arg24 and with the amino group of the peptide bond between Gln25 and Val26 residues. The isopropyl group, on the other hand, also plays a fundamental role in the recognition of monoterpenes; for example, m-cresol and o-cresol, does not have the isopropyl residue and have no affinity on the SCP-2. This observation also agrees with the QSAR models, which propose that nCt are important in biological activity.</p><p>The results on acyclic terpenes denote the importance of π bonds despite not being aromatic moieties. Citronellol, the molecule with lower number of π bonds, is an acyclic terpene less able to interact with SCP-2 and is also the molecule with lower larvicidal activity. On the other hand, myrcene has the highest number of π bonds, presented the highest larvicidal activity and is also the best to interact with SCP-2. Geraniol and myrcene are the acyclic terpenes with the higher larvicidal activity and both interact with the Phe105 residue (Additional file 1: Figure S2). All acyclic terpenes, except those with ketone groups, are capable of interact with residues Ile19, Asn23, Arg24 and Gln25. Geraniol has the ability to form a hydrogen bond with the amino group of the backbone between the Ile19 and Asn20.</p><p>Anisaldehyde presented a binding energy of − 5.72 kcal/mol and was able to interact with the Phe105 residue and form a hydrogen bond with the amino group of Arg24 (Additional file 1: Figure S3a). The cuminaldehyde does not interact with the Phe105 residue and was not able to form hydrogen bonds. Sesquiterpenes presented the highest affinity on the SCP-2 active site, presented interactions with the Phe105 residue and with the hydrophobic pocket (Additional file 1: Figure S3b).</p><!><p>The larvicidal activity of terpenes and terpenoids was analyzed by LC50 determination for different stairs of Culex quinquefasciatus Say. The description of the molecular properties and the structural characteristics responsible for larvicidal activity of the tested compounds, were used for the development of mathematical models of structure–property–activity relationship. The docking studies were able to show that molecular and structural descriptors provide evidence of SCP-2 as a possible biological target, an important protein in cholesterol and fatty acid catabolism, which cleaves the 3-oxoacyl-CoAs of methyl-branched fatty acid and bile acid intermediates. However experimental studies should be conducted to elucidate this effect.</p><!><p>Additional file 1. Additional tables and figures.</p><p>Electronic supplementary material</p><p>The online version of this article (10.1186/s13065-018-0425-2) contains supplementary material, which is available to authorized users.</p>